the rapid neutron capture process ( r - process ) is one of main processes to synthesize elements heavier than the iron group .
the dominant r - process element source in the universe is still matter of debate .
two possible sites have been proposed ; core - collapse supernovae ( ccsne , e.g. , * ? ? ?
* ; * ? ? ?
* ) , and coalescence of binaries of double neutron star ( ns - ns ) and black hole - neutron star ( bh - ns ) systems ( e.g. , * ? ? ?
* ; * ? ? ?
we refer to scenarios based on chemical evolution models of r - process elements exclusively supplied from binary neutron star mergers ( nsms ) as the nsm scenarios in the following .
one promising way to understand the enrichment history of heavy elements in the universe is observations of extremely metal - poor ( emp ) stars .
these stars are the living fossils born in the early universe and can be probes for the galactic chemical evolution
. europium ( eu ) is commonly used as an element representative of r - process elements but we also discuss with barium ( ba ) abundance as another representative element since it is more easily observed in emp stars than eu .
though ba in the solar system is mainly produced by s - process in intermediate mass stars on the asymptotic giant branch ( agb ) , the contribution from s - process is thought to be negligible for emp stars since their formation epoch is in prior to the period when the first intermediate mass stars in the universe evolve to agb stars .
we note that the measured abundance ratios between eu and ba in emp stars are similar to the ratio of the r - process components of these elements in the solar system , except for carbon - enhanced emp stars , of which the surface abundances have probably been changed by binary mass transfer .
it is known that the r - process abundances of emp stars with @xmath0}\lesssim -2.5 $ ] show a very large star - to - star scatter of @xmath1 ( e.g. , * ? ? ?
* ; * ? ? ?
the average abundance of ba shows a decreasing trend as metallicity decreases in the metallicity range of @xmath2 } < -2.6 $ ] , while it seems to reach a plateau at @xmath0 } < -3.3 $ ] @xcite . based on this trend ,
@xcite and @xcite pointed out that we can detect ba on almost all the emp stars .
it has been pointed out from the viewpoint of chemical evolution that r - process elements originating from neutron star mergers ( nsms ) have difficulties in realizing the observed abundance features in emp stars .
@xcite argued that the model with nsm as the dominant r - process source predicts too large abundance scatter since nsms have a very low event rate ( @xmath3 in the milky way ( mw ) , * ? ? ? * ) and very large r - process yield ( @xmath4 , e.g. , * ? ? ?
* ) per event .
in addition , it has been argued that the long delay - time of nsm ( @xmath5yr ) is incompatible with the presence of r - process enhanced stars at @xmath0 } < -2.5 $ ] .
on the other hand , from other points of view , there is growing evidence that the nsms are major r - process sources .
theoretical studies of nucleosynthesis in ccsne show difficulties in synthesizing elements of the second and third r - process peaks , in the framework of neutrino driven explosion ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , while these elements are successfully synthesized in the nsm ejecta ( e.g. , * ? ?
? * ; * ? ? ?
observationally , a short gamma - ray burst grb  130603b showed enhancement of infrared radiation in its afterglow @xcite .
this event agrees with the predicted properties of a  kilonova " , radiation from matter heated by radioactive decay of r - process ejecta from a nsm @xcite .
@xcite pointed out that the discrepancy in abundances of the radioactive r - process element @xmath6pu between the early solar system and the current solar vicinity may indicate a very low event rate for the r - process source .
the discovery of an r - process rich ultra - faint dwarf galaxy @xcite also favors the r - process source with a large yield and rare event rate .
recently , @xcite pointed out that the ejecta from a nsm can spread into very large volume since the nsm ejecta expand at very large velocities of @xmath7 of the speed of light ( e.g. , * ? ? ?
* ; * ? ? ?
the stopping length , @xmath8 , of @xmath9eu with speed of 0.2@xmath10 is @xmath11kpc , which can be larger than the radius of a dwarf spheroidal galaxy .
( see section 2.3.1 for the derivation of @xmath8 .
note that @xcite evaluate this value as @xmath12 , which is in error by a factor of 153 . )
we investigate effects of this large - scale spreading of the nsm ejecta on the galactic chemical evolution .
recent chemical evolution studies confirmed that the emergence of eu around @xmath0}\sim -3 $ ] can not be explained by the nsm scenario , unless binary nss coalesce on very short timescales of 10 myr or less @xcite , while both observations of binary pulsars @xcite and models of binary population synthesis @xcite indicate average coalescence timescales beyond 100 myr .
in addition , even with the short delay time , the secondary r - process source such as ccsne is invoked to reproduce the abundance distribution at @xmath0 } < -3 $ ] @xcite .
@xcite show that the appearance of r - process elements at @xmath0}\sim -3 $ ] can be explained by the combination with the coalescence timescale of 100 myr for the majority ( 95% ) of nsms and 1 myr for the rest ( 5% ) , if the star formation efficiencies are lower for less - massive sub - halos .
but they use a one - zone model and do not account for the scatter of @xmath13 .
@xcite follow the enrichment history of r - process elements by chemo - dynamical simulations but only for a dwarf galaxy . @xcite and
@xcite compute the formation of the mw by chemo - dynamical simulations but their models lack the resolution to correctly follow the formation of pop iii and early pop ii stars .
we note that the results of these simulation studies are dependent on assumptions about the sub - grid mixing process and numerical resolution .
all these previous studies do not take into account the difference of the propagation process between nsm ejecta and sn ejecta . in our previous study , we investigated the enrichment of r - process elements using a chemical evolution model with merging history of proto - galaxies , but without considering a difference in propagation @xcite .
we pointed out that the observed very low frequency of ba lacking stars is in contradiction with the nsm scenario even when we assume a significant surface pollution . in this paper
, we revisit this problem by considering the possibility that nsm ejecta can escape from its host proto - galaxy and pollute neighboring proto - galaxies and intergalactic medium ( igm ) .
the structure of this paper is as follows .
next section describes our model and its ingredients .
section  [ results ] presents results and compares them with observations .
section  [ summarys ] summarizes the paper .
we use the hierarchical chemical evolution model presented in @xcite . in the following ,
we give a brief summary of the model and changes from the previous works .
our model treats individual metal poor stars in the framework of the hierarchical galaxy formation .
this model also considers gas infall and outflow , metal pre - enrichment in igm , and surface pollution of stars by accretion of interstellar matter ( ism ) .
we build merger trees of dark - matter halos based on the extended press - schechter theory @xcite .
the first stars are formed in halos with @xmath14 , where @xmath15 is the virial temperature , and @xmath16 corresponds to the halo mass of @xmath17 at a redshift of @xmath18 .
we refer to the baryonic component in the branches of the merger tree as proto - galaxies , which are low - mass galaxies as building blocks of the milky way .
each halo hosts one proto - galaxy at the beginning .
we assume that when two dark halos merge , the proto - galaxy in the more massive halo is placed at the center of the newly formed halo and the proto - galaxy in the smaller halo becomes a satellite in a sub - halo and that the mass of the smaller halo remains unchanged .
we neglect tidal stripping and ram pressure stripping .
the orbital energy of the satellite decays through dynamical friction @xcite .
the satellite merges with the central galaxy on the dynamical friction timescale , @xmath19 , given by the standard chandrasekhar formula , and stars in the satellite are added to the stellar halo of the central galaxy ( e.g. , * ? ? ?
the star formation rate , @xmath20 , is expressed as the product of the gas mass , @xmath21 , and the star formation efficiency ( sfe ) , @xmath22 .
the masses of individual stars are set randomly following the initial mass function ( imf ) described in section  [ s : imf ] .
metal yields of supernovae ( sne ) are adopted from theoretical results of @xcite for type ii sne and @xcite for type ia sne .
we assume homogeneous mixing of metal in each proto - galaxy .
we consider gas outflow triggered by individual sne .
the first sn in a proto - galaxy form a galactic wind , and subsequent sne add gas and energy to the wind .
for the outflow energy , @xmath23 , we use the following formula as a function of the sn kinetic energy , @xmath24 , and the gas binding energy , @xmath25 , of a proto - galaxy , @xmath26 which supposes that the outflow carries away all the sn kinetic energy if @xmath27 but a fraction , @xmath28 , of the kinetic energy turns to the outflow if @xmath29 .
we compute inhomogeneous metal enrichment in the igm by following the evolution of the winds from proto - galaxies assuming the momentum conserving snow - plow shell model .    in the fiducial model
, we adopt the star formation efficiency depending on the mass of a halo , @xmath30 , motivated by the observational mass - metallicity relation , while a constant efficiency was adopted in our previous studies . the observed correlation between stellar masses , @xmath31 , and average metallicities , @xmath32}\right>$ ] , of galaxies
is expressed as @xmath33}\right > } \propto   m_*^{0.3}$ ] @xcite , and may indicate a lower star formation efficiency in a lower mass galaxy @xcite .
we use @xmath34 instead of @xmath35 to avoid @xmath36 at the beginning .
measurements of rotation curves of low mass star forming galaxies showed that the halo mass is almost proportional to the stellar mass @xcite .
we also present results of the model with a constant sfe but with a higher outflow rate .
we adopt @xmath37 in the constant sfe model rather than @xmath38 as the fiducial value in our previous study .
this model yields a slope of mass - metallicity relation very similar to the fiducial model with @xmath39 and @xmath38 .
we evaluate the proportional constant of the sfe to yield the solar metallicity in the current mw , @xmath40 in the fiducial model and @xmath41 in the constant sfe model .
for the initial mass function ( imf ) , @xmath42 , we use the log - normal function with a power - law tail at higher mass part as @xcite adopted , @xmath43 @xcite yielded @xmath44 , @xmath45 , and @xmath46 for the galactic spheroid . for emp stars , however , we have shown that the low - mass peaked imf obtains a much smaller frequency of carbon enhanced stars and more numerous emp stars than observed by the he s survey @xcite .
a high or intermediate - mass peaked imf shows a better consistency with observations . here , we compute models with the imfs peaked at @xmath47 , and @xmath48 ( @xmath49 is also shifted in proportion to @xmath50 ) .
the binary fraction is set to be 0.5 , and we adopt the flat mass ratio distribution .
we note that the predicted r - process abundance distribution is not sensitive to the choice of the peak mass .
( in previous studies , we may have overestimated the volume observed by the he s survey due to the assumption of @xmath51 for giant stars and applying an incorrect limiting magnitude . the correction slightly lower the estimated value of @xmath50 but it still remains in the intermediate mass range . )
we consider that the surface of a star is polluted by accreting ism .
we trace changes of surface abundance of stars with @xmath0 } < -2.5 $ ] using the bondi accretion formula
. the surface pollution can be significant on the surface abundances of emp stars , as shown in @xcite .
} $ ] of proto - galaxies against redshift , with the gas mass of proto - galaxies color coded .
bottom panel : metallicity distribution functions ( mdfs ) of metal - poor stars in the mw halo .
the green , red , and blue lines show the predicted number of giant stars in the survey volume of the he s survey for models with @xmath52 , respectively .
the dashed red line shows the mdf after surface pollution .
the shaded histogram is the raw ( dark - gray ) and bias - corrected ( light - gray ) mdfs by the he s survey @xcite .
the solid black line is the high - resolution sample biased toward lower metallicity by the saga database @xcite and the dashed black line is the sdss / segue sample by @xcite , scaled to fit the he s data .
the color scale shows the predicted number density and the black crosses are the observational sample from the saga database .
, title="fig : " ] + } $ ] of proto - galaxies against redshift , with the gas mass of proto - galaxies color coded .
bottom panel : metallicity distribution functions ( mdfs ) of metal - poor stars in the mw halo .
the green , red , and blue lines show the predicted number of giant stars in the survey volume of the he s survey for models with @xmath52 , respectively .
the dashed red line shows the mdf after surface pollution .
the shaded histogram is the raw ( dark - gray ) and bias - corrected ( light - gray ) mdfs by the he s survey @xcite .
the solid black line is the high - resolution sample biased toward lower metallicity by the saga database @xcite and the dashed black line is the sdss / segue sample by @xcite , scaled to fit the he s data .
the color scale shows the predicted number density and the black crosses are the observational sample from the saga database .
, title="fig : " ] +    the top panel of figure  [ mdf ] shows the metal enrichment history in our fiducial model .
we plot @xmath0}$ ] of sampled proto - galaxies against redshift .
color denotes the gas mass of proto - galaxies . in our model ,
@xmath0}$ ] can decrease when a proto - galaxy merges with a more metal deficient proto - galaxy .
most of stars with @xmath0}<-3 $ ] are formed at @xmath53 , and stars with @xmath54 } < -2 $ ] are formed at @xmath55 .
we show the resultant metallicity distribution functions ( mdfs ) in the bottom panel .
the green , red , and blue lines show the predicted mdfs in the cases of @xmath52 , respectively . the solid and dashed
red lines are the intrinsic metallicity and the surface metallicity after pollution , respectively .
shaded histograms show the raw ( dark - shaded ) and bias - corrected ( light - shaded ) mdfs of the he s sample @xcite .
we scale the predicted mdfs considering the field of view ( @xmath56 ) and the magnitude limit ( @xmath57 ) of the he s survey @xcite .
since most of the sample stars in @xcite are around the base of the red giant branch , we count the predicted number of giant stars in the survey volume assuming that the typical luminosity is @xmath58 .
we employ the de vaucouleur density profile for the mw halo .
we also show the mdf of the sdss / segue sample ( * ? ? ?
* dashed black line ) scaled to match the he s data and the mdf by the saga database ( solid black line ) , which collects stellar abundances by high - resolution spectroscopy of metal - poor stars from literatures @xcite .    as seen in the bottom panel , though all the three imfs lead to a similar form of the mdf
, the imf with @xmath59 solely reproduces the total number .
the imf with a higher @xmath50 produces a fewer number of low - mass stars which can survive to date , and at the same time form more numerous massive stars producing iron to shorten the timescale for low - metallicity star formation .
we show results of the models with @xmath59 in section  [ results ] .
the main results about r - process elements remain true for the other two @xmath50 .
the mdf with @xmath59 slightly underestimates the number of stars at @xmath0}\gtrsim -2.5 $ ] .
it may be due to the change of the imf to the low - mass peaked one @xcite .
our model obtains the abundance distribution of alpha elements consistent with observations .
we present the correlation of the abundance ratios of magnesium to iron with metallicity in figure  [ mgfe ] , and oxygen and silicon shows the similar abundance distributions .
the mass range for progenitor stars of nss is @xmath60 .
we assume that @xmath61 of binary systems in this mass range form ns binaries which coalesce within the age of the universe .
this yields one coalescing ns binary out of 1000 massive stars .
we adopt the delay time distribution of the nsm , @xmath62 for @xmath63 following results of binary population synthesis by @xcite . in this case , the averaged delay time is @xmath64 .
the nsm event rate is proportional to the following integration , @xmath65    the ejecta mass of a nsm is assumed to be @xmath66 , based on the estimation of ejecta mass of the observed kilonova candidate grb 130603b @xcite .
we assume the scaled solar r - process abundance pattern @xcite for nsm ejecta .
the elemental yields of ba and eu are @xmath67 and @xmath68 , respectively .
we do not distinguish ns - ns mergers and bh - ns mergers in our computation , while recent numerical studies show that a bh - ns merger can eject lager mass than a ns - ns merger @xcite .
in this paper , we do not consider other heavy element sources such as s - process in intermediate massive stars or r - process in sne  ii .
a significant fraction of the nsm ejecta can escape from proto - galaxies .
nearby proto - galaxies can capture the escaped r - process elements . the rest of r - process ejecta travel away into the intergalactic space with high velocities .
the r - process elements cool and are coupled with igm , and are incorporated into galaxies along gas infall as described in the following .
r - process elements from a nsm propagate until they lose kinetic energies by interaction with ism .
the dominant energy deposition process is the ionization of neutral hydrogen in hi gas , and the coulomb scattering in the ionized gas .
the energy loss rate is @xmath69   \nonumber \\ & \times n_{\rm hi } ( 1 + 0.0185 \ln \beta ) \beta^{-1 } \hbox { ev s}^{-1 }      \end{aligned}\ ] ] in the neutral gas and @xmath70 in the ionized plasma @xcite , where @xmath71 and @xmath72 is the number density of neutral hydrogen and free electron , respectively , @xmath73 is the velocity in units of the speed of light and set at @xmath74 in this study , and @xmath75 is the atomic number .
the stopping length , @xmath76 for @xmath77 in the neutral gas , where @xmath78 is the proton mass , @xmath79 is the mass number , and @xmath80 is the lorentz factor .
we assume that r - process elements ejected from an nsm travel the stopping length of @xmath8 on a straight line at a constant velocity .
the ejected particles stop rather suddenly once they travels the stopping length since the energy loss rate increases as the energy decreases .
the fraction , @xmath81 , of nsm ejecta escaping from the proto - galaxy is estimated as @xmath82 where @xmath83 is the radius of the host proto - galaxy . for simplicity , we assume that a proto - galaxy is a homogeneous sphere of neutral gas , and an nsm takes place at the center of the proto - galaxy .
we give the gas density and radius of proto - galaxies assuming isobaric cooling of the virialized gas @xcite , i.e. , @xmath84 , and @xmath85 , where @xmath86 is the averaged density of the virialized dark - matter halos as a function of redshift , @xmath87 is the mean molecular weight , @xmath88 is the virial radius of the host halo .
gas can be cooled to @xmath89 k in the mw but @xmath90 k in the primordial mini - halo , and @xmath91 k in relic hii regions without metal .
the gas temperature , @xmath92 , is set at the maximum of the temperature of the cosmic microwave background or @xmath93 .
these assumptions yield @xmath94 and @xmath95 pc for the primordial proto - galaxy with @xmath96 at @xmath97 , and @xmath98 and @xmath99 kpc for the current mw ( @xmath100 ) .      in this study
, we assume isotropic mass ejection from a nsm @xcite .
the fraction of escaped nsm ejecta encountering with a proto - galaxy with radius @xmath101 is @xmath102 , where @xmath103 is the distance to the proto - galaxy .    for a particle which runs into a spherical proto - galaxy with an impact parameter of @xmath104 ,
the length of the pass crossing the proto - galaxy is @xmath105 . under the same assumptions as in eq .
( [ eq : fesc ] ) , the fraction of particles that are not captured by the proto - galaxy is @xmath106 .
the fraction of captured high - energy r - process elements out of all particles encountering the proto - galaxy is @xmath107 we mention that , in the limit of @xmath108 , the capture rate approaches to @xmath109 and the total amount of captured particles is proportional to @xmath110 since @xmath111 ( eq .  [ eq : ls ] ) , i.e. , the capture rate is simply proportional to the gas mass but independent of radius and density .
we estimate the distance between the @xmath112-th and the @xmath113-th proto - galaxies using the following equation ,    d_n , m(t ) = r(t|m_t , n , m , t_c , n , m ) [ eq : dhalo ] + ( t < t_c ,
n , m ) + d_sat , m(t ) = r_vir , n ( 1- ) [ eq : dsat ] + ( t_c , n ,
m t < t_c , n
, m + _ df ) + [ eq : d2sat ] +    in the case of two proto - galaxies in different dark halos ( eq .  [ eq : dhalo ] ) , we estimate the distance between them under the spherical collapse approximation of halos following @xcite . here
two halos are supposed to merge at time @xmath114 to become a single halo with the mass of @xmath115 , and @xmath104 is the distance between the two halos ( see section 4.1 in @xcite for details ) .
after the two dark halos merge ( eq .  [ eq : dsat ] ) , the proto - galaxy in the more massive halo ( labeled @xmath112 ) is placed at the center of the merged halo and the other one ( labeled @xmath113 ) becomes a satellite .
the satellite galaxy is assumed to approach the center as the distance , @xmath116 , from the center linearly decreases as a function of time .
@xmath117 is the merger timescale by the dynamical friction and we use the formula given in @xcite . at @xmath118 ,
the satellite merges to the central galaxy . in the case of two satellites in the same dark halo ( eq .  [ eq : d2sat ] ) , we estimate the distance assuming the random distribution for the orbital angles , @xmath119 , between the satellite galaxies .
the energy loss timescale , @xmath120 , of @xmath9eu propagating in the igm is given as @xmath121 where we use the energy loss rate for neutral gas ( eq .
[ eq : omeganeutral ] ) before reionization and one for ionized gas ( eq .  [ eq : omegaplasma ] ) after that .
we assume the reionization redshift @xmath122 . since the stopping length in the igm , @xmath123 , is larger than the computation box , we assume homogeneous mixing of the escaped r - process elements in the intergalactic space , while we follow the inhomogeneous metal pollution of the igm for elements provided by sne as mentioned above .
r - process elements having escaped from the host galaxy lose their kinetic energies on the timescale of @xmath124 unless they are captured by another proto - galaxy .
such r - process elements are eventually settled in the igm and some of them are attracted by the gravity of proto - galaxies and accreted onto them along with gas infall .
we can describe the evolution of mass @xmath125 of the high - energy particle in the intergalactic space and the abundance , @xmath126 , of element @xmath127 in the igm as follows , @xmath128 @xmath129 where @xmath130 is the elemental yields from one nsm event , @xmath131 is the time when the @xmath132-th nsm event takes place in the @xmath112-th proto - galaxy . the second term in the right hand side of eq .
( [ eq : he ] ) describes capture of particles incoming from outside of the computation box , and @xmath133 is the volume of the box , which is a comoving sphere enclosing the total mass of @xmath134 .
@xmath135 is the gas accretion rate onto the @xmath112-th proto - galaxy .    here
, we assume that the escaping of the high energy r - process elements from the simulation box balances the incoming from galaxies outside the box .
we discuss the validity of this assumption in appendix .      in summary , evolution of r - process abundances of a proto - galaxy
is described as follows , @xmath136 where the first term on the right - hand side is the contribution from nsms in the proto - galaxy , the second term is the increment by capturing ejecta from nsm events in other proto - galaxies in the simulation box , the third term is high - energy r - process background from galaxies outside the computation box , and the fourth term is the infall of igm .
of proto - galaxies , with the number of nsm events in them color coded .
we randomly sample @xmath137 of proto - galaxies and plot their abundances .
the dash - dotted blue line denote the abundance of high - energy eu in intergalactic space , @xmath138 , and the dashed black line is the low - energy eu abundance in igm . ]
figure  [ ce ] shows the enrichment history of eu abundances for sampled proto - galaxies and the igm ( dashed black line ) .
the number of nsm events in each proto - galaxy is color coded .
the dash - dotted blue line denotes the abundance of high energy r - process elements in intergalactic space , @xmath139 .
the first nsm in this computation run takes place at @xmath140 .
most of the r - process ejecta escape from its host proto - galaxy to enhance @xmath138 to @xmath141 , \end{aligned}\ ] ] where @xmath142 is the solar eu abundance @xcite , and @xmath143 and @xmath144 are densities of matter and baryon relative to the critical density , respectively . in the limit of @xmath108 ,
the capture rate of high energy r - process particles is proportional to the total gas mass and independent of the density as mentioned , i.e. , the capture rate per unit mass is almost the same for all proto - galaxies with @xmath108 and for the igm . of the igm ( black line ) and proto - galaxies ( cyan lines ) increase in a similar way by capturing high - energy r - process elements at @xmath145 . at such a high redshift ,
majority of baryon is still in the igm , and most of the escaped nsm ejecta are mixed with the igm on the timescale of @xmath146 .    at @xmath122
, @xmath138 suddenly decreases because the energy loss rate in the igm increases due to cosmic reionization .
the capture rate per unit mass by the ionized igm becomes several tens of times higher than the one by neutral proto - galaxies .
therefore , the abundance of proto - galaxies without nsm events ( cyan lines ) is @xmath147 dex or more smaller than the igm abundance at @xmath148 . after the reionization , pre - enrichment of igm is mainly responsible for the eu abundance of proto - galaxies without nsm events . at @xmath149
, @xmath138 overwhelms @xmath150 because the energy loss timescale in the igm increases as the igm density decreases .
a single nsm event enriches its host proto - galaxy with r - process elements to @xmath151      \end{aligned}\ ] ] in the limit of @xmath152 ( @xmath108 ) .
when we adopt assumptions for @xmath112 and @xmath101 in section  [ s : fesc ] , we can describe this equation as a function of @xmath153 and the cosmological parameters as follows , @xmath154 ,     \ ] ] where @xmath155 is the overdensity of the virialized halo relative to the critical density , and @xmath156 is the dimensionless hubble parameter .
these proto - galaxies become much more r - process rich than the intergalactic space , and the contribution from other proto - galaxies is negligible for them . at @xmath157 ,
the most massive photo - galaxy shows the highest r - process abundance because of a large event rate and a small escape probability .
figures  [ rfe ] and [ outflow ] show the r - process abundance of stars against metallicity for the fiducial model and the constant sfe model , respectively .
the left and right panels denote results without escape of nsm ejecta ( @xmath158 ) and with escape , respectively .    the constant sfe significantly delays the enrichment of r - process elements compared to iron and the predicted abundances of emp stars do not match observational results , as mentioned in previous chemical evolution studies .
the predicted abundances of r - process elements are significantly lower than observations at @xmath0}\lesssim -2 $ ] , and r - ii stars with @xmath159 around @xmath0}\sim -3 $ ] is not reproduced . on the other hand , if the sfe is lower at lower mass proto - galaxies , r - process elements emerge at @xmath0}\sim -3 $ ] , even with the long average delay time ( @xmath160 yr ) .
this is because that enrichment timescale of iron is longer for low - mass proto - galaxies in which emp stars were formed and some nsms can take place within the formation epoch of emp stars .
these are consistent with the result of the one zone model by @xcite .
when we do not consider the escape of nsm ejecta from the host proto - galaxies , however , the majority of stars with @xmath0 } < -3 $ ] show very low r - process abundances with [ ba / fe ] @xmath161 , which is apparently at odds with observations ( see the top left panel of fig .
3 ) . on the other hand ,
our escape model predicts @xmath162 or @xmath163 for the majority of stars at @xmath0}\lesssim -3 $ ] , and well reproduces observations of halo stars , as shown in the right top panel .
proto - galaxies lacking nsm events are polluted by r - process elements ejected from nsms in other proto - galaxies .
as mentioned , igm pre - enrichment is mostly responsible for the r - process elements in these stars , and capture of high - energy r - process elements also contribute at @xmath164 .
stars in proto - galaxies without a nsm are distributed at @xmath165 .
nsm events of short delay time , @xmath166 , are the dominant sources of r - process elements for these stars . on the other hand , proto - galaxies in which nsm events take place form stars with @xmath167 .
surface pollution can play an important role for stars in proto - galaxies without a nsm . without surface pollution , there emerges a depression in the number density of stars at @xmath0}\sim -3 $ ] in the range of @xmath168 where the first nsm occurs .
the surface pollution alleviates the gap . in the model with @xmath158 ,
the surfaces of most emp stars have a tiny amount of ba such as @xmath169 even after surface pollution .    at @xmath0 } > -2 $ ] , the r - process abundance is averaged .
the predicted ba abundance is lower than observations since we neglect the s - process contribution .    at @xmath0}\gtrsim -1 $ ] , while we can see a clear decreasing trend of @xmath170 as metallicity increases , we do not see a similar  knee " in the predicted r - process abundance distributions .
this is because the increase of the nsm rate due to the long delay time compensates the increasing production rate of iron enhanced by sne ia .
one possible explanation for the observed depression of @xmath13 is a lower frequency of nsms with high metallicity progenitor stars , although we assume constant frequency in the computation for simplicity .
the event rate of nsms in this model is @xmath171 in the current mw , as shown in figure  [ nsmrate ] .
it is also slightly higher than the estimation based on the observations of binary pulsars in the mw ( @xmath172 , * ? ? ?
both the binary pulsars and the abundances of the metal rich stars indicate a factor of @xmath173 smaller nsm fraction than used in this paper at @xmath174 . at the right axis of figure  [ nsmrate ] ,
the event rate in the mw is converted into the expected detection rate of the gravitational waves ( gws ) by the advanced ligo - virgo network .
ongoing observations of gws will provide better constraints on the nsm rate .    for the metallicity dependence of the nsm rate ,
studies of binary population synthesis predict a lower ns - ns merger rate for lower metallicity stars @xcite .
for the bh - ns merger rate , on the other hand , @xcite predict that a rate for stars with the metallicity of @xmath175 is higher than the population i ( pop i ) stars while @xcite predict that a rate for pop iii stars is lower than that for pop i stars .
further studies are required for the event rate , delay time distribution , and r - process yield of nsms , and their metallicity dependences .    in this paper , we assume that high - energy r - process elements go straight until they lose energies by interaction with gas particles . if proto - galaxies have magnetic field , however , particles can travel along a winding path , which decreases the escape probability and increases the capture rate . on the other hand , the escape probability can be higher in a disk - like gas density profile than in a gas sphere assumed in this work .
in addition , nsm events at the outer region of a proto - galaxy will effectively transfer r - process elements into intergalactic space , although we only consider nsm events at the center of proto - galaxies . furthermore , a ns binary itself can go away from a proto - galaxy if it can get enough velocity from sn explosions .
the velocity of a ns binary is depending on its binary orbit and sn kick .
although we still have no reliable estimation for the velocity distribution of ns binaries , two observed binary pulsars show a proper motion of more than @xmath176 km / s @xcite , which exceeds the typical escape velocity of a proto - galaxy hosting emp stars .
the igm pre - enrichment and transfer of r - process elements between proto - galaxies are depending on these details , but a model taking them into account is beyond the scope of this paper .
coalescence of ns binaries is one of the most plausible sites for r - process nucleosynthesis but results of previous chemical evolution models with nsms as major r - process sources is inconsistent with the observed r - process abundances of metal - poor stars .
we revisit the galactic chemical evolution of r - process elements following the nsm scenario .
the novelty of this paper is in considering the effect of propagation of nsm ejecta across proto - galaxies .
as pointed out by @xcite , r - process elements from nsms have very large stopping lengths due to their high velocities .
we compute the escape of r - process elements from proto - galaxies , pre - enrichment of igm , and capture of the escaped elements by other proto - galaxies in the chemical evolution model based on the framework of the hierarchical galaxy formation .
we have confirmed that the nsm scenario can reproduce the emergence of r - process elements at @xmath0}\sim -3 $ ] in the case of lower sfe for less massive proto - galaxies , while the constant sfe delays the enrichment of r - process elements , which results in abundance correlations incompatible with observations .    while previous chemical evolution models with the nsm scenario predict many stars lacking r - process elements at very low metallicities , our model with the long - distance propagation of nsm ejecta predicts @xmath165 for stars formed in proto - galaxies without a nsm .
our model succeeds in reproducing the r - process abundance distribution of emp stars .
this result indicates that the nsm scenario is not rejected by the galactic chemical evolution , when we consider the propagation of nsm ejecta beyond the proto - galaxies .
our current model predicts a slightly higher event rate of nsms and results in @xmath177 higher than observations in the current mw while the good consistency with observations at very low metallicity .
it may indicate the metallicity dependence of the nsm rate .
future observations of gws and their electromagnetic counterparts will give us information about the nsm rate and its contribution to the galactic chemical evolution .
we adopt a periodic boundary condition which assumes that the high energy r - process elements incoming from outside the simulation box balance the elements escaping from the box .
this seems to be a reasonable assumption at very high redshift since large scale density fluctuation tends to grow at low redshift .
for example , at @xmath97 , only @xmath178 fluctuation of @xmath179 collapses , and the box size is much larger . on the other hand ,
it is not a good approximation at low redshift since the environment around the current mw is significantly different from the average of the universe .    in the current universe , there are @xmath180 mw equivalent galaxies ( mwegs ) on average @xcite , but the number density in our computation box corresponds to @xmath181 ( @xmath182 ) .
it may indicate that the production rate of high energy r - process elements ( the first term of eq .
[ eq : he ] ) may be overestimated by factor @xmath183 at @xmath174 . the capture rate by proto - galaxies ( the second term of eq .
[ eq : he ] ) may also be overestimated since the capture rate in this computation box is dominated by the mw at low redshift . in order to mimic this difference , we compute the chemical evolution model with replacing the eq .
( [ eq : he ] ) to the following equation , @xmath184 where @xmath185 is the relative number density of mwegs .
but for the different boundary condition . incoming rate of high energy r - process particles from boundaries are smaller than outgoing rate by proportional to the relative density of milky way equivalent galaxies .
see text for details . ]
figure  [ boundary ] show the result .
the predicted abundance distribution is indistinguishable from the fiducial model at high metallicity or high r - process abundance ( @xmath167 ) .
this is because r - process elements of proto - galaxies in which these stars are formed are dominated by their own nsms . on the other hand , the abundance of proto - galaxies without nsm events decreases since the r - process element abundances of the igm decrease .
the abundance of emp stars at @xmath186 formed before the first nsm in their host proto - galaxies decreases by @xmath187 dex .
we note that the difference can be overestimated at low metallicity since the region around the mw should be more similar to the average of the universe at high redshift .
abadie , j. , abbott , b. p. , abbott , r. , et al . , 2010 ,
cqgra , 27 , 3001 andrievsky , s. m. , spite , m. , korotin , s. a. , 2009 , , 494 , 1083 arlandini , c. , kappeler , f. , wisshak , k. , gallino , r. , lugaro , m. , busso , m. , & straniero , o. 1999 , , 525 , 886 argast , d. , samland , m , thielemann , f .- k . , & qian , y .- z .
2004 , , 416 , 997 bauswein , a. , goriely , s. , & janka , h .-
, 2013 , , 773 , 78 barnes , j. , & kasen , d. , 2013 , , 775 , 18 beniamini , p. , & piran , t. 2016 , , 456 , 4089 belczynski , k. , kalogera , v. , & bulik , t. 2002 , , 572 , 407 berger , e. , fong , w. , chornock , r. 2013 , , 774 , l23 burris , d.  l. , pilachowski , c.  a. , armandroff , t.  e. , sneden , c. , cowan , j.  j. , & roe , h. 2000 , , 544 , 302 carollo , d. , beers , t. c. , chiba , m. , et al .
2010 , , 712 , 692 cescutti , g. , romano , d. , matteucci , f. , chiappini , c. , & hirschi , r. 2015 , , 577 , 139 chabrier , g.  2003 , , 115 , 763 christlieb , n. , schrck , t. , frebel , a. , et al .
2008 , , 484 , 721 cole , s. , lacey , c. g. , baugh , c. m. , & frenk , c. s. , 319 , 1 dominik , m. , belczynski , k. , fryer , c. , et al .
2012 , , 759 , 52 goriely , s. , bauswein , a. , & janka , h .- t .
2011 , , 738 , l32 grevesse , n. , noels , a. , & sauval , a. j. 1996 , aspc , 99 , 117 hirai , y. , ishimaru , y. , saitoh , t. r. , fujii , m. s. , hidaka , j. , & kajino , t. 2015 , , 814 , 41 honda , s. , aoki , w. , kajino , t. , ando , h. , beers , t.  c. , izumiura , h. , sadakane , k. , takada - hida , m. 2004 , , 607 , 474 hotokezaka , k. , kiuchi , k. , kyutoku , k. , et al .
2013 , phrvd , 87 , 4001 hotokezaka , k. , kyutoku , k. , tanaka , m. , et al .
2013 , , 778 , l16 hotokezaka , k. , piran , t. , & paul , m. natph , 11 , 1042 lacey , c. , & cole , s.  1993 , , 262 , 627 ishimaru , y. , wanajo , s. , & prantzos , n. , 2015 , , 804 , l35 ji , a. p. , frebel , a. , chiti , a. , & simon , j. d. , nature , 531 , 610 kim , c. , perera , b.  b.  p. , mclaughlin , m.  a. , 2015 , , 448 , 928 kinugawa , t. , inayoshi , k. , hotokezaka , k. , et al .
2014 , , 442 , 2963 kirby , e. n. , cohen , j. g. , guhathakurta , p. , cheng , l. , bullock , j. s. , & gallazzi , a. 2013 , , 779 , 102 kobayashi , c. , umeda , h. , nomoto , k. , tominaga , n. , & ohkubo , t. 2006 , , 653 , 1145 komiya , y. , habe , a. , suda , t. , & fujimoto , y. m. 2009 , , 696l , 79 komiya , y. , suda , t. , & fujimoto , y. m. 2009 , , 694 , 1577 komiya , y. , yamada , s. , suda , t. , & fujimoto , m.  y. 2014 , , 783 , 132 komiya , y. , suda , t. , & fujimoto , y. m. 2015 , , 808 , l47 komiya , y. , suda , t. , & fujimoto , y. m. 2015 , , accepted .
lattimer , j. m. , & schramm , d. n. , 1974 , , 192 , 145 lorimer , d. r. 2008 , llr , 11 , 8 matteucci , f. , romano , d. , arcones , a. , korobkin , o. , & rosswog , s. 2014 , , 438 , 2177 miller , s. h. , ellis , r. s. , newman , a. b. , & benson , a. 2014 , , 782 , 115 nomoto , k. , thielemann , f .- k . , yokoi , k. 1984 , , 286 , 644 roederer , i. u. 2013 , , 145 , 26 roederer , i. u. , karakas , a. i. , pignatari , & herwig , f. , 2016 , , 821 , 37 rosswog , s. , liebendrfer , m. , thielemann , f .- k . , davies , m. b. , benz , w. , & piran , t. 1999 , , 341 , 499 sato , k. , 1974 , pthph , 51 , 726 schlickeiser , r. 2002 , cosmic ray astrophysics ( berlin : springer ) schrck , t. , christlieb , n. , et al .
2009 , , 507 , 817 shen , s. , cooke , r. j. , ramirez - ruiz , e. , madau , p. , mayer , l. , guedes , j. 2015 , , 807 , 115 somerville , r.  s. , & kolatt , t.  s.  1999 , , 305 , 1 suda , t. , et al .
2008 , pasj , 60 , 1159 suda , t. , et al .
2013 , , 432 , l46 tanaka , m. , & hotokezaka , k. , , 775 , 113 tanvir , n. r. , levan , a. j. , fruchter , a. s. , et al .
2013 , nature , 500 , 547 tsujimoto , t. , & shigeyama , t. 1999 , , 565 , l5 van de voort , f . ,
quataert , e. , hopkins , p. f. , keres , d. , & faucher - gigure , c - a . , 2015 , , 447 , 140 wanajo , s. , 2013 , , 770 , l22 wanajo , s. , sekiguchi , y. , nishimura , n. , et al . 2014 , , 789 , l39 wanajo , s. , janka , h - t . , & mller , b. 2011 , , 726 , l15 wehmeyer , b. , pignatari , m. , & thielemann , f .- k . , 2015 , , 452 , 1970 woosley , s.  e. & hoffman , r.  d. 1992 , , 395 , 202 yamada , s. , suda , t. , komiya , y. , aoki , w. , & fujimoto , m.  y. , , 436 , 1362


light - cone gauge string field theory provides a useful way to define string theories @xcite@xcite@xcite .
given the action , it is possible to define the amplitudes and calculate them perturbatively , although we should check if they are well - defined .
since it is a gauge fixed theory , there is no problem in considering a light - cone gauge string field theory in noncritical space - time dimensions .
being noncritical , it does not possess the space - time lorentz invariance .
it should correspond to a string theory in a lorentz noninvariant background .
in other words , we should be able to find a brst invariant worldsheet theory , with a nonstandard @xmath0 part .
what we would like to do in this paper is to study this @xmath0 theory .
we give the energy - momentum tensor and the action of the theory and calculate the correlation functions .
we show that the energy - momentum tensor possesses desired properties , and construct a brst invariant formulation of the worldsheet theory .
the reason why we are interested in this theory is that it can be used to regularize the string field theory .
the dimensional regularization is one of the most powerful regularizations for ordinary quantum field theory and it may be useful also in string theory . in particular , in the light - cone gauge superstring field theory , unwanted divergences occur even at the tree level , because of the transverse supercurrent insertions at the interaction points of the vertices . in ref .
@xcite , we have proposed a dimensional regularization scheme to deal with these divergences .
we have found that the divergences of the tree amplitudes can be regularized by shifting the number of the space - time dimensions .
we have checked that the results of the first quantized formalism are reproduced without any counterterms for the four - point case . in order to proceed
further we need to show that the dimensional regularization preserves important symmetries of the theory .
if the light - cone gauge string field theory corresponds to a brst invariant formulation even in noncritical dimensions , it means that the dimensional regularization preserves the brst symmetry .    in this paper
, we deal with the case of closed bosonic string field theory .
we restrict ourselves to the tree amplitude and consider the worldsheet theory on the complex plane .
the organization of this paper is as follows . in section [ sec : lccov ]
, we present a way to rewrite light - cone gauge string amplitudes in a brst invariant form .
we follow the procedure in noncritical dimensions , and show what kind of worldsheet theory of @xmath0 should appear in the end . in section [ sec : xpmcft ] , we study the theory for @xmath3 and show that it is a cft with the right virasoro central charge .
namely , combining the cft for @xmath3 constructed here with the worldsheet theory for the light - cone gauge strings in noncritical dimensions , we obtain a cft with the central charge @xmath4 by which we can define a brst invariant worldsheet theory . in section [ sec : brst ] , we show that the tree amplitudes in the noncritical case can be written in a brst invariant form .
section [ sec : discussions ] is devoted to discussions . in appendix
[ sec : action ] , we present the action of the light - cone gauge bosonic string field theory in @xmath5 dimensions . in appendix
[ sec : mandelstam ] , the mandelstam mapping is given . in appendix
[ sec : gamma ] , we present a way to calculate @xmath6 $ ] , which is used in sections [ sec : lccov ] and [ sec : xpmcft ] . in appendix
[ sec : relations ] , the details of the calculations in section [ sec : xpmcft ] are presented .
in order to study the noncritical case , it is useful to consider the relation between light - cone gauge amplitudes and covariant ones for critical strings .
one can calculate the amplitudes starting from the light - cone gauge string field theory action given in appendix [ sec : action ] .
the tree amplitudes can be expressed by path integral on the light - cone diagrams , on which the complex @xmath7 coordinate is introduced as usual . via the mandelstam mapping @xmath8 , which is given in appendix [ sec : mandelstam ]
, one can express them using correlation functions of vertex operators on the @xmath9-plane endowed with the metric @xmath10 the correlation functions we should consider are @xmath11_{\phi }   e^{-s_{x^{i } } } \prod_{r=1}^{n}v_{r}^{\mathrm{lc}}\;. \label{eq : correlation}\ ] ] here @xmath12 denotes the free action for the @xmath13 transverse coordinates @xmath14 and @xmath15 denotes the vertex operator .
we assume that the @xmath16-th external state is of the form @xmath17 and the vertex operator should be @xmath18 where @xmath19 is defined in eq.([eq : rho - zir ] ) , @xmath20 is the coordinate of the unit disk of the @xmath16-th external string given in eq.([eq : rho - wr ] ) and @xmath21 .
the on - shell and the level - matching conditions require that @xmath22    the path integral measure @xmath23_{\phi}$ ] should be defined using the metric ( [ eq : metric ] ) .
it is related to the measure @xmath23 $ ] which is defined with the flat metric @xmath24 as @xmath25_{\phi }   & \propto &   \left[dx^{i}\right]e^{-\gamma\left[\phi\right]}\;,\end{aligned}\ ] ] roughly speaking , @xmath26 $ ] is given by the liouville action @xmath27
\sim   -\frac{24}{24\pi}\int d^{2}z \ ,   \partial \phi \bar{\partial}\phi ~ , \label{eq : liouville}\ ] ] where @xmath28 .
since @xmath29 diverges at the poles and zeros of @xmath30 by the definition ( [ eq : metric ] ) , the liouville action is not well - defined at these points .
$ ] can be evaluated by regularizing the divergences and carefully taking various effects into account @xcite@xcite .
we present an alternative derivation of @xmath31 $ ] in appendix [ sec : gamma ] and give an explicit form in eq.([eq : gamma - c ] ) .
, this form of @xmath31 $ ] is only up to a constant which can be fixed by factorization . with the string field action in eq.([eq : sftaction ] ) and the @xmath31 $ ] in eq.([eq : gamma - c ] ) , one can see that the factor is fixed as @xcite @xmath25_{\phi }   & \sim &   \left[dx^{i}\right ] \mathrm{sgn}\left(\prod_{r=1}^n\alpha_r\right ) e^{-\gamma\left[\phi\right]}\;,\end{aligned}\ ] ] up to a numerical factor .
we ignore the phase factor @xmath32 in the following , because it does not play any important roles . ]
thus eq.([eq : correlation ] ) can be rewritten as @xmath33 e^{-s_{x^{i}}-\gamma[\phi ] } \prod_{r}v_{r}^{\mathrm{lc}}\;. \label{eq : correlation2}\ ] ]      in order to covariantize the amplitudes , we need to introduce the longitudinal coordinates .
the light - cone gauge condition implies that @xmath34 equals the lorentzian time on the light - cone diagram , which means on the @xmath9-plane , @xmath35 therefore we introduce the variable @xmath34 with the delta function @xmath36 , which can be expressed as @xmath37 e^{-\frac{1}{\pi}\int d^{2}z   x^{\prime -}\partial\bar{\partial}x^+ } \prod_{r=1}^{n } e^{-ip_{r}^{+}x^{\prime-}}\left(z_{r},\bar{z}_{r}\right ) \nonumber \\ & & \hspace{5 mm } \times \int \left [ db^{\prime}dc^{\prime}d\tilde{b}^{\prime}d\tilde{c}^{\prime } \right ] e^{-\frac{1}{\pi}\int d^{2}z     \left ( b^{\prime}\bar{\partial}c^{\prime }            + \tilde{b}^{\prime } \partial \tilde{c}^{\prime}\right ) } c'(\infty ) \tilde{c}'(\infty ) ~.    \label{eq : x+delta}\end{aligned}\ ] ] eq.([eq : x+delta ] ) should be considered as the formal euclideanized version of a lorentzian path integral .
the grassmann odd fields @xmath38 of conformal weights @xmath39 , @xmath40 , @xmath41 , @xmath42 are introduced to cancel the determinant factor @xmath43 .
with these variables , we can rewrite the right hand side of eq.([eq : correlation2 ] ) as @xmath44 e^{-s_{x^{i}}-s_\pm^\prime -s_{b^{\prime}c^{\prime } } } \nonumber\\ & & \hspace{10em } \times c^\prime ( \infty ) \tilde{c}^\prime ( \infty ) \prod_{r=1}^{n }   \left(v_{r}^{\mathrm{lc}}e^{-ip_{r}^{+}x^{\prime- } } \left(z_{r},\bar{z}_{r}\right ) \right ) \ ; , \label{eq : correlation3}\end{aligned}\ ] ] where @xmath45 \;. \label{eq : xprimeaction}\end{aligned}\ ] ] we consider @xmath46 as the worldsheet action for these variables .
the action @xmath47 is with the interaction term which depends only on @xmath34 .
this interaction term does not affect the correlation functions with less than two insertions of @xmath48 , and they coincide with those of the free theory .
thus one can show the ope @xmath49    from the worldsheet action , the energy - momentum tensor can be obtained as @xmath50 when @xmath51 . here
@xmath52 denotes the schwarzian derivative .
this term can be derived from the variation of @xmath53 which coincides with the liouville action except for the singular points .
the covariant expression for the amplitudes can be obtained by introducing @xcite @xmath54 and their anti - holomorphic counterparts @xmath55 and @xmath56 .
the fields @xmath57 and @xmath58 have now weights @xmath59 and @xmath60 respectively .
we should also introduce @xmath61 so that the ope s between @xmath62 and the ghosts are regular and the energy - momentum tensor ( [ eq : emtensor ] ) takes the form @xmath63 which coincides with the energy - momentum tensor in the conformal gauge . in the following , we would like to rewrite the correlation function ( [ eq : correlation3 ] ) using these new variables .
let us first consider the vertex operator . from eq.([eq :
emtensor2 ] ) , one can see that the action for the new variables should be the free action .
@xmath64 appears essentially in the same way as @xmath65 does in eq.([eq : x+delta ] ) , and we obtain the delta function . therefore @xmath8 and @xmath66 which appear in the integrand can be replaced by @xmath67 and @xmath68 , where @xmath69 and @xmath70 are the holomorphic and the anti - holomorphic part of @xmath34 respectively .
we will denote the equality which holds under this identification by @xmath71 .
thus the factor @xmath72 in the definition  ( [ eq : vrlc ] ) of @xmath73 can be rewritten as @xmath74 where @xmath75 is the ddf operator given by @xmath76 one can also use @xmath77 for @xmath78 , where @xmath79 is defined in eq.([eq : nbarrr00 ] ) .
using these , one can show that @xmath80 in eq.([eq : correlation3 ] ) subject to the on - shell condition ( [ eq : on - shell - cov ] ) can be rewritten as @xmath81 where @xmath82 is the vertex operator corresponding to the ddf state .
one subtle point to notice in eq.([eq : defbc ] ) is that the variables @xmath83 should have zeros and poles at the zeros and poles of @xmath84 . since the integration over @xmath62 leads to the identification @xmath85 , the zeros and poles of @xmath84 become those of @xmath30 , namely @xmath86 and @xmath87 .
accordingly @xmath57 should be inserted at @xmath86 and @xmath58 should be inserted at @xmath87 , if we wish to rewrite eq.([eq : correlation3 ] ) @xcite@xcite@xcite@xcite .
these insertions should come with appropriate factors made from @xmath88 or @xmath34 in order that eq.([eq : correlation3 ] ) is reproduced .
using eqs.([eq : lcddf ] ) and ( [ eq : exp - gamma2 ] ) , one can show that eq.([eq : correlation3 ] ) can be rewritten as @xmath89    e^{-s_{x}-s_{bc } } \ ,    \left|\sum_{r}\alpha_{r}z_{r}\right|^{2 }   \left (   \lim_{z\to\infty}\frac{1 }   { \left|z\right|^{4 } }   c\left(z\right)\tilde{c}\left(\bar{z}\right)\right ) \nonumber \\   &   & \hspace{1 cm }   \times\prod_{i }   \left (   \frac{b}{\partial^{2}\rho}(z_{i } )   \frac{\tilde{b}}{\bar{\partial}^{2}\bar{\rho}}(\bar{z}_{i})\right )   \prod_{r=1}^{n } \left(c\tilde{c } v_{r}^{\mathrm{ddf}}\right )            \left(z_{r},\bar{z}_{r}\right ) \ ; ,   \label{eq : correlation4}\end{aligned}\ ] ] where @xmath90 and @xmath91 are the free actions for @xmath92 and for @xmath83 respectively .
one can show that eq.([eq : correlation4 ] ) yields the expression for the amplitudes in the conformal gauge . using @xmath93 and deforming the contour
, one can recast eq.([eq : correlation4 ] ) into the form @xmath89
e^{-s_{x}-s_{bc}}\nonumber \\   &   & \hspace{1 cm }   \times\prod_{i }   \left (     \oint_{c_{i}}\frac{dz}{2\pi i}\frac{b}{\partial\rho}(z )     \oint_{c_{i } }     \frac{d\bar{z}}{2\pi i }       \frac{\tilde{b}}{\bar{\partial}\bar{\rho } } ( \bar{z } )   \right )   \prod_{r=1}^{n }     \left(c\tilde{c}v_{r}^{\mathrm{ddf}}\right )            \left(z_{r},\bar{z}_{r}\right)\ ; ,   \label{eq : correlation5}\end{aligned}\ ] ] where the integration contour @xmath94 is depicted in fig .
[ fig : npt ] in appendix [ sec : gamma ] on the @xmath7-plane .
the amplitudes can be obtained by integrating the correlation function @xmath95 over the @xmath96 moduli parameters @xmath97 defined in eq.([eq : mathcalt ] ) as @xmath98 one can see that the antighost insertion @xmath99 corresponds to the quasiconformal vector field associated with the moduli parameter @xmath100 .
the form of the amplitude is brst invariant , because the vertex operator @xmath101 is invariant and @xmath102 yields a total derivative with respect to @xmath100 .
we can change the integration variables to @xmath103 and @xmath104 with the expression ( [ eq : correlation5 ] ) for @xmath95 , the determinant factor can be combined with the antighost factors as @xmath105 performing the contour integrals around @xmath106 in the last line , we eventually obtain the familiar expression @xmath107    e^{-s_{x}-s_{bc}}\nonumber \\   &   & \hspace{1 cm }   \times   \prod_{s=1,2,n }
\left(c\tilde{c}v_{s}^{\mathrm{ddf}}\right )      \left(z_{s},\bar{z}_{s}\right )   \prod_{r=3}^{n-1 }   \int d^2z_r   v_{r}^{\mathrm{ddf}}\left(z_{r},\bar{z}_{r}\right)\;.\end{aligned}\ ] ]    hence the variables @xmath108 can be identified with those in the covariant formulation . from eq.([eq
: newvariables ] ) , one can find the ope s of the variables @xmath109 to be @xmath110 and regular otherwise .
these can also be deduced from the action ( [ eq : xprimeaction ] ) .
using these ope s , one can easily show that the energy - momentum tensor ( [ eq : emtensor ] ) satisfies the virasoro algebra of central charge @xmath111 .      the light - cone gauge string field theory in @xmath5
@xmath112 space - time dimensions can be defined with the action given in appendix [ sec : action ] .
this time , the correlation functions we should consider is_{\phi }   & \sim &   \left[dx^{i}\right ] \mathrm{sgn}\left(\prod_{r=1}^n\alpha_r\right ) e^{-\frac{d-2}{24}\gamma\left[\phi\right]}\;,\end{aligned}\ ] ] and we ignore the phase @xmath32 in the following . ]
@xmath33 e^{-s_{x^{i}}-\frac{d-2}{24}\gamma[\phi ] } \prod_{r}v_{r}^{\mathrm{lc}}\;. \label{eq : noncriticalf}\ ] ] we can follow the above procedure and introduce the variables @xmath113 without any problem . in this case , eq.([eq : emtensor ] ) should be @xmath114 by using @xmath108 defined in eq.([eq : newvariables ] ) , @xmath115 can be rewritten as @xmath116 and the action @xmath117 should be @xmath118\;. \label{eq : sxpm}\ ] ] therefore in noncritical dimensions the worldsheet theory for @xmath3 is different from the usual free theory , and obviously the lorentz symmetry is broken .
what we should study is this theory for @xmath3 .
it is a conformal field theory similar to the one for @xmath119 and @xmath65 in the previous subsection .
we will study its properties in the next section .    before closing this subsection ,
a comment is in order . for @xmath120
, one can define @xmath121 with @xmath122 using these variables , the energy - momentum tensor ( [ eq : nenergymomentum ] ) can be written as @xmath123 therefore the worldsheet theory is a free theory with ghosts of noninteger spins
. we may be able to study the theory using these variables , although we need to figure out the way to deal with the ghost zero - modes .
the theory we would like to consider is with the action ( [ eq : sxpm ] ) and the energy - momentum tensor @xmath124 this theory will be well - defined if @xmath84 has a nonvanishing expectation value . as in the case of the theory for @xmath119 and @xmath65
, we always consider this theory in the presence of the insertions of the vertex operators @xmath125 ( @xmath126 ) with @xmath127 , so that the classical equation of motion implies @xmath128 thus the quantities we would like to calculate are the expectation values @xmath129 \prod_{r=1}^{n }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle } \nonumber\\ & & \equiv     \int\left[dx^{+}dx^{-}\right ]      e^{-s_{x^\pm}}f\left[x^{+},x^{-}\right ]      \prod_{r=1}^{n }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right),\end{aligned}\ ] ] for functionals @xmath130 $ ] which satisfy @xmath131 = f\left[x^{+},x^{-}\right]~,\ ] ] for arbitrary constants @xmath132 .    for functionals @xmath133 $ ] which do not depend on @xmath64 , it is formally possible to perform the path integral as in eq.([eq : x+delta ] ) and obtain @xmath134 \prod_{r=1}^{n }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle
\sim f\left[-\frac{i}{2}\left(\rho+\bar{\rho}\right)\right ] \exp \left(-\frac{d-26}{24}\gamma     \left[\ln\left(\partial \rho \bar{\partial}\bar{\rho}\right)\right ] \right ) ~ , \label{eq : fx+}\ ] ] up to the factor coming from the integration over the zero - modes of @xmath0 .
it is convenient to define @xmath135\right\rangle_{\rho } \equiv \frac { \left\langle   f\left[x^{+},x^{-}\right ] \prod_{r=1}^{n }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle } { \left\langle   \prod_{r=1}^{n }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle } ~,\ ] ] and we obtain @xmath136\right\rangle _ { \rho }    =   f\left[-\frac{i}{2}\left(\rho+\bar{\rho}\right)\right]\;. \label{eq : expx+}\ ] ]      the manipulation to obtain eq.([eq : fx+ ] ) is rather formal because @xmath138 possesses zeros and poles and we need to specify the regularization procedure to define @xmath139 $ ] . as we mentioned in the last section , this was done in ref .
@xcite and @xmath140 $ ] is given for arbitrary mandelstam mapping @xmath7 .
assuming that eq.([eq : fx+ ] ) is true with @xmath53 derived in ref .
@xcite ( and also in appendix [ sec : gamma ] ) , one can calculate correlation functions which involve @xmath64 . roughly speaking , differentiating the left hand side of eq.([eq : fx+ ] ) with respect to @xmath141 and setting @xmath142
, we obtain @xmath134x^-(z_n,\bar{z}_n ) \prod_{r=1}^{n-1 }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle ~,\ ] ] although the momentum conservation condition complicates the procedure a bit . in the same way , we can in principle
obtain arbitrary correlation functions , treating the right hand side of eq.([eq : fx+ ] ) as a kind of generating functional .
since all the correlation functions on the complex plane are given this way , we define the theory based on eq.([eq : fx+ ] ) , rather than starting from the action ( [ eq : sxpm ] ) .
if one define the theory in this way , it is not a priori clear if @xmath143 in eq.([eq : emtensorxpm ] ) can be considered as the energy - momentum tensor . in the rest of this section
, we would like to calculate the correlation functions of @xmath143 and examine if it can be regarded as the energy - momentum tensor of the theory .
one comment is in order .
the operator @xmath144 can be considered to create a hole of length @xmath145 in the light - cone diagram .
therefore it is similar to the macroscopic observables in the old matrix models @xcite .
these operators are nonlocal objects and give rise to singularities at @xmath146 where no operators are inserted .
taking @xmath147 limit to obtain local operators is exactly what was done in the old matrix models .      in order to calculate the correlation functions of @xmath143 , we need correlation functions of @xmath64 .
let us follow the above - mentioned procedure and calculate them .
first we consider the simplest example , @xmath148 . in order to calculate this quantity
we start from @xmath149 \right ) ~,\ ] ] where @xmath150 and @xmath151 then the one point function @xmath148 can be given as @xmath152                    \right )                \right|_{\alpha_{0}=0 }         \;. \label{eq : expdx-}\ ] ]    it is straightforward to evaluate the right hand side of eq.([eq : expdx- ] ) by using the expression  ( [ eq : gamma - c ] ) for @xmath6 $ ] . in the calculations , one should notice the following points . while there are @xmath153 interaction points @xmath154 for @xmath155
, there are only @xmath156 interaction points @xmath86 for @xmath157 . since @xmath158 what happens in the limit @xmath159 is as follows : one of @xmath154 s , which we denote by @xmath160 , goes to @xmath161 ; one of @xmath154 s , which we denote by @xmath162 , goes to @xmath163 ; the other interaction points @xmath164 go to the interaction points @xmath86 for @xmath7 , which are denoted with the same subscripts . as we mentioned earlier , @xmath165 with finite @xmath166 should be considered as a nonlocal operator , which induces singularities at the interaction points where there are no operator insertions .
@xmath167 in the limit @xmath168 is consistent with the fact that @xmath165 tends to a local operator .
one subtle point to notice is that if @xmath169 is close to one of @xmath170 and the interaction point nearest to @xmath171 coincides with @xmath160 , the neumann coefficient @xmath172 does not go to @xmath173 in the limit @xmath168 .
since @xmath174 $ ] depends on @xmath172 , it implies that @xmath165 has a nonlocal effect even if we take @xmath168 in such a configuration .
then we can not expect to get correlation functions of local operators from @xmath174 $ ] . in the following
, we will assume that @xmath169 and @xmath163 are not close to any of @xmath170 , to avoid such a situation .
namely we assume that @xmath175 and @xmath176 are not in the regions of external propagators in the @xmath7 plane .
we calculate @xmath177 for such @xmath169 , and analytically continue the result to the whole complex plane .
one can generalize eq.([eq : expdx- ] ) and calculate the correlation functions with insertions of @xmath34 as @xmath180
\right\rangle _ { \rho } \nonumber \\ & & \hspace{5 mm } =           \left .
\frac { 2i \partial_{z_0 } \partial_{\alpha_0 } \left\langle   f\left[x^{+}\right ] \prod_{r=0}^{n+1 }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle } { \left\langle   \prod_{r=0}^{n+1 }          e^{-ip_{r}^{+}x^{-}}\left(z_{r } , \bar{z}_{r}\right ) \right\rangle }                           \right|_{\alpha_{0}=0 }     \nonumber \\   & &   \hspace{5 mm }   =      \left\langle           \partial x^{-}\left(z_{0}\right )      \right\rangle _ { \rho }      f
\left [ -\frac{i}{2 } \left ( \rho+\bar{\rho } \right ) \right ]    + \int d^{2}z \
, \frac{1}{z_{0}-z }      \left .
\frac{\delta f\left[x^{+}\right ] }                  { \delta x^{+ } ( z,\bar{z } ) }      \right|_{x^{+}=-\frac{i}{2}\left(\rho+\bar{\rho}\right)}\;. \label{eq : expxpm}\end{aligned}\ ] ] since we are dealing with the correlation functions with source terms for @xmath64 , we can read off the operator relations from these correlation functions . from eqs.([eq :
expx+ ] ) and ( [ eq : expxpm ] ) , we obtain the ope s @xmath181 which are valid if @xmath9 and @xmath182 are away from the singularities @xmath183 .
these are consistent with eq.([eq : ope ] ) .
the two point function for @xmath184 can be calculated by using eq.([eq : dx- ] ) as @xmath185 here we are interested in the singularity at @xmath186 .
it is straightforward to calculate the right hand side of eq.([eq : twopoint ] ) and obtain @xmath187   + \mbox{regular terms}\;. \label{eq : dx-2}\ ] ] from this , we deduce the ope @xmath188   \nonumber\\   & \sim &   { } -\frac{d-26}{12 }     \left[- \frac{1}{(z - z')^{4 } }            \frac{6}{\left(\partial x^{+}(z')\right)^{2 } }          - \frac{1}{(z - z')^{3 } }            3 \partial \left (                   \frac{1}{\left ( \partial x^{+ } ( z ' ) \right)^{2 } }                 \right )     \right .
\nonumber\\   & & \hspace{5em } \left .          { } -\frac{1}{(z - z')^{2 } } \frac{1}{2 }            \partial^{2 } \left (                   \frac{1}{\left ( \partial x^{+ } ( z ' ) \right)^{2 } }                   \right )     \right]~ , \label{eq : ope3}\end{aligned}\ ] ] which is valid if @xmath9 and @xmath182 are away from the singularities @xmath183 .
this is consistent with eq.([eq : ope ] ) .
we can also obtain an expression for the correlation functions with @xmath34 insertions as we did in eq.([eq : expxpm ] ) : @xmath189    \right\rangle _ { \rho } } \nonumber\\   & = & \left\langle \partial x^{-}\left(z\right )              \partial x^{-}\left(w\right)\right\rangle _ { \rho }        f\left[-\frac{i}{2}\left(\rho+\bar{\rho}\right)\right ] \nonumber \\   &   & \quad { } + \left\langle \partial x^{-}\left(z\right )                 \right\rangle _ { \rho }         \int d^{2}w ' \ , \frac{1}{w - w ' }             \left.\frac{\delta f\left[x^{+}\right ] }                        { \delta x^{+ } ( w',\bar{w}^\prime ) }             \right|_{x^{+}=-\frac{i}{2}\left(\rho+\bar{\rho}\right ) }
\nonumber \\   &   & \quad { }     + \left\langle \partial x^{-}\left(w\right )                 \right\rangle _ { \rho }         \int d^{2}z ' \ , \frac{1}{z - z ' }              \left.\frac{\delta f\left[x^{+}\right ] }                         { \delta x^{+ } ( z',\bar{z}^\prime ) }              \right|_{x^{+}=-\frac{i}{2}\left(\rho+\bar{\rho}\right ) } \nonumber \\   &   & \quad { } + \int d^{2}z ' \int d^{2}w ' \ ,       \frac{1}{z - z ' }   \frac{1}{w - w ' }       \left.\frac{\delta^{2}f\left[x^{+}\right ] }                  { \delta x^{+}(z',\bar{z } ' )                     \delta x^{+ } ( w',\bar{w } ' ) }       \right|_{x^{+}=-\frac{i}{2}\left(\rho+\bar{\rho}\right)}\;. \label{eq : expdx - dx - f}\end{aligned}\ ] ]        one can evaluate the right hand side of eq.([eq : vevtpm ] ) and examine how it behaves around the possible singularities : @xmath194 \right ) }              { \partial z_{r } }       & z\sim z_{r }      \\[1ex ]       \displaystyle         \mathrm{regular } & z\sim z_{i }        \\[0.5ex ]         \mathcal{o } \left(\frac{1}{z^{4 } } \right )         & z\sim \infty    \end{array } \right . .
\label{eq : vevtpm2}\ ] ] eqs.([eq : vevtpm2 ] ) and ( [ eq : exptf ] ) imply that @xmath195 is regular at @xmath196 and @xmath197 , if no operators are inserted there . therefore , although @xmath198 is singular at @xmath199 without any operator insertions , the energy momentum tensor is regular and conserved .
this property is essential for constructing the brst charge . from eq.([eq : vevtpm2 ] ) for @xmath200 we can read off the ope @xmath201 although @xmath202 is a nonlocal operator , it behaves as a primary field of weight @xmath203 .
using the ope s ( [ eq : ope - xpm ] ) and ( [ eq : ope3 ] ) , we can show that @xmath143 satisfies @xmath204 therefore the central charge of the virasoro algebra in the @xmath3 cft is @xmath205 .
we can also find that @xmath206 and thus @xmath207 are primary fields of weight @xmath208 .
since the worldsheet theory for @xmath0 is a cft with virasoro central charge @xmath205 , with the transverse coordinates @xmath14 added the total central charge of the system for @xmath0 and @xmath14 is @xmath4 .
therefore with ghosts @xmath57 and @xmath58 , we can construct a nilpotent brst charge @xmath209 .    as we have shown in section [ sec : lccov ] , the amplitude for the light - cone gauge string field theory can be rewritten by using these variables .
we start from the correlation function given in eq.([eq : noncriticalf ] ) , where the vertex operator @xmath73 is of the form ( [ eq : vrlc ] ) but with the on - shell condition @xmath210 we can proceed in the same way as in section [ sec : lccov ] and eventually obtain @xmath211   e^{-s_{x^{i}}-s_{x^{\pm}}-s_{bc } }   \nonumber\\   &   & \hspace{1 cm }   \times   \prod_{i }   \left (   \oint_{c_{i}}\frac{dz}{2\pi i}\frac{b}{\partial\rho}(z )   \oint_{c_{i } }   \frac{d\bar{z}}{2\pi i}\frac{\tilde{b}}{\bar{\partial}\bar{\rho } }    ( \bar{z } )   \right )   \nonumber   \\   & &   \hspace{1 cm }   \times   \prod_{r=1}^{n } \left(c\tilde{c}v_{r}^{\mathrm{ddf } }   \exp \left(-i\frac{d-26}{24}\frac{x^+}{p_r^+}\right )   \right )   \left(z_{r},\bar{z}_{r}\right )   \nonumber   \\   & &   \hspace{1 cm }   \times   \prod_{r=1}^{n }   \exp \left(i\frac{d-26}{24}\frac{x^+}{p_r^+}\right )   \left ( z_i^{(r)},\bar{z}_i^{(r ) } \right)~ , \label{eq : brstinv}\end{aligned}\ ] ] where @xmath212 is defined in the same way as in eq.([eq : vddf ] ) .
it is easy to show that @xmath213 is a primary field of weight @xmath214 and @xmath143 is regular even with the insertion @xmath215 .
therefore eq.([eq : brstinv ] ) gives a brst invariant expression for the amplitude .
@xmath213 may look as a vertex operator with momentum @xmath216 instead of @xmath217 , but the momentum which is conserved is @xmath218 .
the conserved momentum can be identified with the operator @xmath219 on the worldsheet , which is conserved at the interaction points with insertions @xmath220 .
in this paper , we have constructed a brst invariant worldsheet theory which corresponds to the light - cone gauge string field theory in @xmath221 space - time dimensions . the worldsheet theory for the longitudinal coordinate variables
@xmath0 is different from the usual free theory , but it is a cft with @xmath222 .
our results provide yet another way to construct string theories in noncritical dimensions
. the brst invariant formulation will be useful to study d - branes for such string theories .
now that the cft is given , we can at least formally construct the interaction vertices of the string field theory based on this cft through the prescription of ref .
@xcite . since we have constructed the cft on the worldsheet of the light - cone string diagram , the gauge unfixed version of the string field theory
is supposed to possess the joining - splitting type of interactions .
such a theory looks similar to the @xmath223 hikko theory given in ref .
@xcite .
the results in this paper should be generalized to be used in regularizing string field theory .
one should consider the @xmath0 on the riemann surfaces with higher genera , in order to check if it works for regularizing the uv and ir divergences .
one should also construct a supersymmetric version of the cft . the light - cone gauge superstring field theory in noncritical dimensions
can be used to dimensionally regularize the tree amplitudes of the critical theory , as was discussed @xcite .
it is possible to generalize the calculations performed in this paper into the superstring case , although they are much more complicated .
we will present these results elsewhere .
this work was supported in part by grant - in - aid for scientific research  ( c ) ( 20540247 ) and grant - in - aid for young scientists  ( b ) ( 19740164 ) from the ministry of education , culture , sports , science and technology ( mext ) .
in order to fix the notation , we present the action for the light - cone gauge bosonic string field theory in @xmath5 space - time dimensions .
the worldsheet variables are @xmath224 .
the action of the string field theory takes the form @xmath225~. \label{eq : sftaction}\end{aligned}\ ] ] here @xmath226 is the zero - mode of the transverse virasoro generators for the @xmath16-th string , @xmath227 is the coupling constant , and @xmath228 denotes the integration measure for the momentum zero - modes of the @xmath16-th string defined as @xmath229 where @xmath230 is the string - length parameter of the @xmath16-th string .
@xmath231 is the reflector given by @xmath232 @xmath233 denotes the three - string interaction vertex defined as @xmath234 here @xmath235 is the three - string lpp vertex @xcite for the transverse coordinate variables @xmath14 and @xmath236 .
the string field @xmath237 is taken to obey the reality and the level - matching conditions .
let @xmath7 be the standard complex coordinate on the @xmath153-string tree diagram with the joining - splitting type interaction .
the portion on the @xmath7-plane corresponding to the @xmath16-th external string @xmath238 is mapped to the unit disk , @xmath239 , of the @xmath16-th string as @xmath240 where @xmath241 is the coordinate on the @xmath7-plane at which the @xmath16-th string interacts .
the @xmath153-string tree diagram is mapped to the complex @xmath9-plane with @xmath153 punctures by the mandelstam mapping @xcite @xmath242 where the puncture @xmath243 correspond to the origin of the unit disk @xmath244 .
the @xmath156 interaction points @xmath245 are determined by @xmath246 .
these are related to the interaction points on the @xmath7-plane by @xmath247 where @xmath248 denotes one of @xmath86 at which the @xmath16-th external string interacts .
the fact that @xmath86 s are the zeros of @xmath249 yields @xmath250
@xmath6 $ ] can be obtained by evaluating the liouville action for the metric ( [ eq : metric ] ) on the @xmath9-plane  @xcite@xcite . here
we will present an alternative derivation .
@xmath251}$ ] can be regarded as the partition function of the light - cone gauge string theory in @xmath4 space - time dimensions from the view point of the worldsheet cft on the light - cone diagram .
therefore the variation @xmath252 $ ] under @xmath253 can be given by using the expectation value of the transverse energy - momentum tensor @xmath254 .
we can obtain @xmath255 $ ] by integrating @xmath256 $ ] .
the result should be up to a factor which does not change under the variation eq.([eq : deltazr ] ) .
such a factor can be fixed by imposing the factorization condition in each case @xcite .    before we begin the calculation , two comments are in order .
because of the constraint @xmath257 , all the @xmath258 s can not be treated as independent variations . here
we think of @xmath258 with @xmath259 as independent of each other and @xmath260 as being determined by the relation @xmath261 we also note that @xmath255 $ ] involves divergence originating from the infinite length of the external lines .
we regularize it by cutting off the @xmath16-th external line so that its length becomes @xmath262 , where @xmath263 is a large positive constant . by doing so
, the contribution from the @xmath16-th external line to the partition function becomes @xmath264 , and the divergences of @xmath255 $ ] become independent of the parameters @xmath265 . on the @xmath9-plane
, this corresponds to cutting a hole of radius @xmath266 around @xmath171 , where @xmath266 and @xmath263 are related as @xmath267      the variations ( [ eq : deltazr ] ) correspond to the variations of the following parameters of the light - cone diagram : @xmath269 the moduli parameters @xmath270 @xmath271 the heights @xmath272 of the external cylinders ; @xmath273 the circumferences @xmath274 of the external cylinders .
thus the variation @xmath275)$ ] is expressed as @xmath276 \right ) = ( \mathrm{i } ) + ( \mathrm{ii } ) + ( \mathrm{iii})~ , \label{eq : variation0}\ ] ] where @xmath277     + \mathrm{c.c.}~ ,    \nonumber\\   ( \mathrm{iii } )    & = & \sum_{r=1}^{n-1 }       i 2\pi \delta \alpha_{r } \int_{l_{rn } } \frac{d\rho}{2\pi i }      \left\langle t_{x^{i } } ( \rho ) \right\rangle
+ \mathrm{c.c.}~. \label{eq : variation1}\end{aligned}\ ] ] here @xmath278 stands for the complex conjugate .
the integration contour @xmath94 of the term @xmath269 lies between the consecutive interaction points @xmath279 and @xmath280 as depicted in fig .  [ fig : npt ] .
the integration path @xmath281 of term @xmath273 in eq.([eq : variation1 ] ) is a line stretching from the asymptotic region of the @xmath16-th external string to that of the @xmath153-th string on the @xmath7-plane . as an example , the path @xmath282 is depicted in fig .  [ fig : npt ] . on the @xmath9-plane ,
the path @xmath281 becomes a segment connecting the two punctures @xmath87 and @xmath283 with the orientation from @xmath87 to @xmath283 .
-string tree diagram with @xmath284 .
the contours @xmath94 and the path @xmath282 on the @xmath7-plane of the integrals in eq.([eq : variation1 ] ) are depicted for this case . ]
the expectation value of @xmath285 can be evaluated by going to the @xmath9-plane . since @xmath286
we obtain @xmath287 using eq.([eq : vevt ] ) , one can rearrange the term @xmath269 in eq.([eq : variation1 ] ) as @xmath288 where @xmath289 denotes the interaction point closest to @xmath290 .    from eq.([eq :
delta - alphan ] ) , one can easily find that term @xmath271 becomes @xmath291    we can recast the @xmath16-th term of @xmath273 in eq.([eq : variation1 ] ) into @xmath292 here we have taken the cut of the function @xmath293 in the integrand on the right hand side to be the segment @xmath281 and the integration path @xmath294 to be the sum of the two oriented segments connecting two punctures @xmath87 and @xmath283 as depicted in fig .  [
fig : cut ] . deforming the contours to rearrange eq.([eq
: log - int ] ) further and summing over @xmath16 , we obtain @xmath295 where @xmath296 and @xmath297 are depicted in fig .
[ fig : cut ] .     consists of the two oriented segments described by the arrowed bold lines .
the dashed line denotes the cut of the integrand on the right hand side of eq.([eq : log - int ] ) . ]
eventually we obtain the expression @xmath298 \right ) & = & \left [ { } -\sum_{r=1}^{n } \int_{c_{r } } \frac{dz}{2\pi i }       \frac{\delta \rho ( z ) - \delta \rho ( z^{(r)}_{i } )              + \delta \alpha_{r } \ell_{r } }            { \partial \rho ( z ) }       \left(- 2   \ { \rho , z\ } \right )      \right .
\nonumber\\ & & \quad { }   - \sum_{i } \oint_{z_{i } } \frac{dz}{2\pi i }       \frac{\delta \rho ( z ) - \delta \rho ( z_{i } ) }            { \partial \rho ( z ) }       \left(- 2   \ { \rho ,
z\ } \right ) \nonumber\\ & & \left .
\quad { } - \oint_{\infty } \frac{dz}{2\pi i }     \frac{\delta \rho ( z ) - \delta \rho ( z^{\infty}_{i } ) }          { \partial \rho ( z ) }     \left ( - 2 \ { \rho , z\ } \right )     \right ] + \mathrm{c.c.}~. \label{eq : variation2}\end{aligned}\ ] ]      let us evaluate the contour integrals on the right hand side of eq.([eq : variation2 ] ) . for @xmath299
, the schwarzian derivative @xmath300 behaves as @xmath301 \displaystyle   \frac{3}{(z - z_{i})^{2 } }       + \frac{1}{z - z_{i } }          \frac{\partial \left(-w\right)}{\partial z_{i } }       + \mathcal{o } \left ( ( z - z_{i})^{0 } \right ) & z \sim z_{i } \\[2ex ] \mathcal{o } \left(\frac{1}{z^{4 } } \right ) & z \sim \infty \end{array } \right.~ , \label{eq : schwarzian - around}\ ] ] where @xmath179 is a function of @xmath302 defined as @xmath303    for the contour integral along @xmath296 , one should notice that @xmath304 involves a term @xmath305 .
this part can be integrated as @xmath306 for @xmath307 .
the other terms can be evaluated by taking the residues of the simple poles and one can show @xmath308 where @xmath79 is a neumann coefficient given by @xmath309 and @xmath310 is defined in eq.([eq : rho - zir ] ) .    for @xmath311
, one can show @xmath312 using @xmath313 which is derived by varying the equation @xmath314 under eq.([eq : deltazr ] ) .
the contour integral around @xmath146 can be obtained as @xmath315    it is easy to see that the integral around @xmath197 vanishes .
putting all the pieces together , we obtain @xmath316 = - w      -2 \sum_{r=1}^{n }         \mathop{\mathrm{re } } \bar{n}^{rr}_{00 }        -3\sum_{i}\ln         \left| \partial^{2 } \rho ( z_{i } ) \right| ~. \label{eq : gamma - c}\ ] ] this form of @xmath6 $ ] is useful in the calculations in section [ sec : xpmcft ] .
we note that @xmath317 can be described as @xmath318 which follows from eqs.([eq : deldelrho ] ) and ( [ eq : zr - zs ] ) .
this yields @xmath319 }   =    \prod_{r=1}^{n } \left| \alpha_{r }   \right|^{-2 }       \left| \sum_{s=1}^{n } \alpha_{s } z_{s }       \right|^{4 }       e^{-2 \sum_{r=1}^{n }            \mathop{\mathrm{re } } \bar{n}^{rr}_{00 } }      \prod_{i } \left| \partial^{2 } \rho ( z_{i } )                \right|^{-1}~. \label{eq : exp - gamma}\ ] ] this expression is the one obtained by the method in ref .
we also note that @xmath319 }   =    \prod_{r=1}^{n } \left| \alpha_{r }   \right|^{-1 }    e^{-\frac{1}{2}w }       \left| \sum_{s=1}^{n } \alpha_{s } z_{s }       \right|^{2 }       e^{-2 \sum_{r=1}^{n }            \mathop{\mathrm{re } } \bar{n}^{rr}_{00 } }      \prod_{i } \left| \partial^{2 } \rho ( z_{i } )                \right|^{-2}~ , \label{eq : exp - gamma2}\ ] ] which is used in sections [ sec : lccov ] and [ sec : brst ] .
in this appendix , we present the details of the calculations to obtain eqs.([eq : dx- ] ) and ( [ eq : dx-2 ] ) .
let us consider the interaction points @xmath321 and @xmath162 for @xmath322 defined in eq.([eq : drhohat ] ) , which tend to @xmath323 and @xmath163 as @xmath159 respectively . from eq.([eq
: drhohat ] ) , we can obtain the expansion of @xmath324 and @xmath162 in terms of @xmath325 , @xmath326~ ,   \nonumber \\ & & z'_{i^{(n+1)}}-z_{n+1 } \nonumber\\ & & \qquad \quad   =   \frac{\alpha_{0}}{\partial\rho\left(z_{n+1}\right ) }   \left[1 -\left ( \frac{\partial^{2}\rho\left(z_{n+1}\right ) }                { \left(\partial\rho\left(z_{n+1}\right)\right)^{2 } }          + \frac{1}{\partial\rho\left(z_{n+1}\right ) }            \frac{1}{z_{n+1}-z_{0 } }    \right )    \alpha_{0 }    + \mathcal{o } \left ( \alpha_{0}^{2 } \right)\right].~~~~~~~ \label{eq : zi0-z0}\end{aligned}\ ] ] the neumann coefficients @xmath327 for the mandelstam mapping @xmath155 behave as @xmath328 using eqs.([eq : delrho2 ] ) and ( [ eq : zi0-z0 ] ) , we obtain @xmath329\alpha_0 + \mathcal{o}\left ( \alpha_{0}^2 \right),\end{aligned}\ ] ] and @xmath330 \ ; , \nonumber\\     \partial^{2 } \rho '             \left(z'_{i^{\left(n+1\right)}}\right )   & = & \frac{\left(\partial\rho\left(z_{n+1}\right)\right)^{2 } }      { \alpha_{0 } } \nonumber\\ & & \ \times \left [ 1      + \left ( \frac{3\partial^{2}\rho\left(z_{n+1}\right ) }           { \left(\partial\rho\left(z_{n+1}\right)\right)^{2 } }     + \frac{2}{\partial\rho\left(z_{n+1}\right)\left(z_{n+1}-z_{0}\right ) } \right )   \alpha_0     + \mathcal{o } \left(\alpha_{0}^2 \right)\right]~.~~~~~ \ ] ] gathering all the relations obtained above , we have @xmath331   = 6 \ln |\alpha_0 | - 4   + \gamma [ \ln\left(\partial\rho\bar{\partial}\bar{\rho}\right ) ] } \nonumber\\ & & \hspace{4em}{}+ 2 \mathop{\mathrm{re } }    \left
[ \sum_{r=1}^n \frac{1}{\alpha_r }      \ln \left|\frac{\left(z_i^{(r ) } - z_0\right )                       \left(z_r- z_{n+1}\right ) }                 { \left(z_r - z_0\right )                        \left(z_i^{(r ) } - z_{n+1}\right ) } \right|\right . \nonumber\\ & & \hspace{8em } { }         -\sum_{i } \frac{1 } { \partial^{2 } \rho ( z_{i } ) }                 \left ( \frac{1}{z_{i } -z_{0 } }                        - \frac{1}{z_{i } - z_{n+1 } }                 \right )                 \frac{\partial w } { \partial z_{i } } \nonumber\\ & & \hspace{8em } { }   -\frac{3}{2}\sum_{i }             \frac{1}{\partial^2 \rho(z_i ) }   \left (             \frac{1}{\left(z_{i}-z_{0}\right)^{2 } }             -\frac{1}{\left(z_{i}-z_{n+1}\right)^{2 } }
\right ) \nonumber\\ & & \hspace{8em }   \left . { } -\frac{3}{2}\sum_{i }             \frac{\partial^3 \rho\left(z_i\right ) }                 { \left ( \partial^{2}\rho\left(z_{i}\right ) \right)^2 }            \left (               \frac{1}{z_{i}-z_{0}}-\frac{1}{z_{i}-z_{n+1 } }            \right )     \right]\alpha_0 + \mathcal{o}\left ( \alpha_{0}^2 \right).~~~~\end{aligned}\ ] ] we can see that @xmath332 $ ] is divergent and does not coincide with @xmath333 $ ] . this singularity can be avoided by modifying @xmath334 , which corresponds to a renormalization of the operator @xmath335 .
such a renormalization is irrelevant to the calculation of the right hand side of eq.([eq : expdx- ] ) and we obtain @xmath336    \right|_{\alpha_0=0 }   = \partial_{z_0 } \left [ \sum_{r=1}^n \frac{1}{\alpha_r }   \ln \left(\frac{z_i^{(r ) } - z_0 }                     { z_r - z_0 } \right)\right .
-\sum_{i } \frac{1 } { \partial^{2 } \rho ( z_{i } ) }                  \frac{1}{z_{i } -z_{0 } }                 \frac{\partial w } { \partial z_{i } } } \nonumber\\ & & { } \hspace{12em }   \left.-\frac{3}{2}\sum_{i }             \frac{1}{\partial^2 \rho(z_i ) }              \left (              \frac{1}{\left(z_{i}-z_{0}\right)^{2 } }             + \frac{\partial^3
\rho\left(z_i\right ) }                 { \partial^{2}\rho\left(z_{i}\right ) }               \frac{1}{z_{i}-z_{0 } }             \right )     \right ] .
\label{eq : dgamma}\end{aligned}\ ] ] from this equation , we can compute the right hand side of eq.([eq : expdx- ] ) and obtain eq.([eq : dx- ] ) .
let us evaluate the right hand side of eq.([eq : twopoint ] ) .
this can be evaluated by using eq.([eq : dgamma ] ) with @xmath161 , @xmath337 and @xmath338 replaced by @xmath9 , @xmath339 and @xmath340 respectively and with the range of the index @xmath16 taken to be from @xmath203 to @xmath341 .
the terms in which we are interested are the @xmath342 contribution of the first term and the @xmath343 case of the second and the third terms in the square brackets on the right hand side of eq.([eq : dgamma ] ) with the replacements mentioned above . in order to evaluate the first term
, we use @xmath344\alpha_{0 } \nonumber \\   &   & \qquad      { } + \mathcal{o } \left ( \alpha_{0}^{2 } \right)\ ; , \label{eq : for - rprime0}\end{aligned}\ ] ] which follows from eq.([eq : zi0-z0 ] ) . for the computation of the second and the third terms , we use @xmath345       \alpha_{0}^{2 }    + \mathcal{o } \left ( \alpha_{0}^{3 } \right)\ ; , \nonumber \\ \frac{1}{z - z'_{i^{(0 ) } } }   & = &     \frac{1}{z - z_{0 } }     - \frac{1}{\partial\rho ( z_{0 } ) }      \frac{1 } { ( z - z_{0})^{2 } } \alpha_{0 }     + \mathcal{o } \left ( \alpha_{0}^{2 } \right)\ ; , \nonumber \\ \frac{\partial w^{\prime } }       { \partial z'_{i^{(0 ) } } }   & = &   - \frac{\partial^3\rho'(z'_{i^{(0)}})}{\partial^2\rho'(z'_{i^{(0 ) } } ) }   = - 2     \partial \rho ( z_{0 } )     \frac{1}{\alpha_{0 } }     + \frac{2}{z_{0}-z_{n+1 } }    + \mathcal{o } ( \alpha_{0})~. \label{eq : for - i0}\end{aligned}\ ] ] combining these relations , we obtain eq.([eq : dx-2 ] ) .
m.  kaku and k.  kikkawa , `` the field theory of relativistic strings , pt .
1 . trees , '' phys .
d * 10 * , 1110 ( 1974 ) ; `` the field theory of relativistic strings . 2 .
loops and pomerons , '' phys .
d * 10 * , 1823 ( 1974 ) .
y.  baba , n.  ishibashi and k.  murakami , `` light - cone gauge superstring field theory and dimensional regularization , '' arxiv:0906.3577 [ hep - th ] .
s.  mandelstam , `` the interacting string picture and functional integration , '' in _ workshop on unified string theories , _ eds .
green and d.j .  gross ( world scientific , singapore , 1986 ) , p46 .
m.  b.  green , j.  h.  schwarz and e.  witten , `` superstring theory .
2 : loop amplitudes , anomalies and phenomenology , '' _ cambridge , uk : univ . pr . ( 1987 ) 596 p. ( cambridge monographs on mathematical physics ) _
e.  dhoker and s.  b.  giddings , `` unitary of the closed bosonic polyakov string , '' nucl .
b * 291 * , 90 ( 1987 ) .
h.  kunitomo and k.  suehiro , `` a new method of constructing gauge string field theory , '' nucl .
b * 289 * , 157 ( 1987 ) .
h.  sonoda , `` functional determinants on punctured riemann surfaces and their application to string theory , '' nucl .
b * 294 * , 157 ( 1987 ) .
e.  dhoker and d.  h.  phong , `` functional determinants on mandelstam diagrams , '' commun .
* 124 * , 629 ( 1989 ) .
t.  banks , m.  r.  douglas , n.  seiberg and s.  h.  shenker , `` microscopic and macroscopic loops in nonperturbative two - dimensional gravity , '' phys .  lett .
b * 238 * , 279 ( 1990 ) .