{'fraction_non_alphanumeric': 0.042173840020507926, 'fraction_numerical': 0.017263783181733803, 'mean_word_length': 4.004777915413201, 'pattern_counts': {'":': 0, '<': 13, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 44, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
String Field Theory as World-sheet UV Regulator

Even at tree level, the first quantized string theory suffers from apparent short distance singularities associated with collision of vertex operators that prevent us from straightforward numerical computation of various quantities. Examples include string theory S-matrix for generic external momenta and computation of the spectrum of string theory under a marginal deformation of the world-sheet theory. The former requires us to define the S-matrix via analytic continuation or as limits of contour integrals in complexified moduli space, while the latter requires us to use an ultraviolet cut-off at intermediate steps. In contrast, string field theory does not suffer from such divergences. In this paper we show how string field theory can be used to generate an explicit algorithm for computing tree level amplitudes in any string theory that does not suffer from any short distance divergence from integration over the world-sheet variables. We also use string field theory to compute second order mass shift of string states under a marginal deformation without having to use any cut-off at intermediate steps. We carry out the analysis in a broad class of string field theories, thereby making it manifest that the final results are independent of the extra data that go into the formulation of string field theory. We also comment on the generalization of this analysis to higher genus amplitudes.


Introduction and summary
String field theory was originally formulated with the hope of getting non-perturbative results in string theory. With the exception of some non-perturbative classical solutions in open string field theory [1], this hope has not been realized so far. However string field theory can play another useful role -in getting a better description of perturbative string theory. Indeed, superstring field theory has been useful for giving a definition of perturbative string amplitudes free from all divergences other than the ones expected on physical grounds, e.g. the usual infrared divergences in four space-time dimensions. It has also been useful for proving general properties of the perturbative S-matrix e.g. unitarity, analyticity etc. A review of these developments can be found in [2].
The goal of this paper will be to use string field theory to remove some of the inconvenient features of tree level string theory -divergences associated with short distance singularities on the world-sheet. We analyze two classes of problems -computation of on-shell amplitudes and studying the effect of marginal deformation on the world-sheet (super-) conformal field theory (CFT) on the spectrum of string theory. Since at tree level, theories with tachyons e.g. bosonic string theory and superstring theory formulated around certain non-supersymmetric backgrounds, also give sensible results, our analysis will also be valid for the corresponding string field theories. For this reason we shall use the phrase 'string field theory' instead of 'superstring field theory', with the understanding that superstring field theory should be considered as a special case of string field theory. We shall see that all the apparent short distance divergences on the world-sheet arise from wrong treatment of (L 0 +L 0 ) −1 where L n ,L n are the usual Virasoro generators of matter ghost CFT and, for superstring theory, also from wrong treatment of picture changing operators (PCOs) [3].
There are many versions of string field theory -for our analysis we shall work with a broad class of string field theories whose interaction vertices are constructed in such a way that the equality of the S-matrix computed from string field theory and first quantized string theory is manifest. This includes for example the original open bosonic string field theory of Witten [4], its generalization to closed bosonic string field theory in [5][6][7][8][9] and closed and open superstring field theories described e.g. in [10][11][12][13][14][15], but does not include some versions of string field theory where this equality is not manifest e.g. those described in [16][17][18][19][20]. Although our analysis can be applied to both open and closed string theories, for definiteness we shall focus on closed string theories. For open string theories, the L 0 +L 0 factors in the subsequent discussion will be replaced by L 0 and there will be no analog of the projection to L 0 =L 0 states. Let us first consider tree level scattering amplitudes in string theory. The usual definition of these amplitudes is in terms of integrals of certain correlation functions of vertex operators in a CFT, but this often suffers from divergences from regions of integration where the locations of the vertex operators come close to each other. A simple example of this is the Koba-Nielsen formula for multi-tachyon amplitudes in bosonic string theory. Usually we avoid this problem by defining these integrals via analytic continuation, but this prevents us from directly evaluating the integrals numerically when an analytic formula is absent. Refs. [21][22][23] suggest alternative approaches by regarding the integrals as (limits of) integrals in the complexified moduli space -some aspects of this will be discussed later.
In string field theory the tree level S-matrix is evaluated by summing over Feynman diagrams. Part of the contribution to the S-matrix with N external legs comes from the diagram involving an N-point interaction vertex to which the external legs connect, without any internal propagator. 1 This diagram gives the same world-sheet integral as the first quantized string theory, except that the integration is over a restricted region, avoiding regions where two or more vertex operators come close to each other. These missing contributions are given by Feynman diagrams with one or more internal propagators and interaction vertices with less number of external legs. Formally the contribution from the latter diagrams also can be represented as the same world-sheet integral as the first quantized string theory, but the integration runs over regions where one or more vertex operators come close, complementing the region associated with the the N-string interaction vertex. However this correspondence is only formal. To understand the difference between the formal expression and the actual contribution from the Feynman diagram of string field theory, we first note that in Siegel gauge [24] the propagator of the string field takes the form with some additional numerator factors for the Ramond sector states in superstring theories.
Here b 0 ,b 0 denote the b,b ghost zero modes and L n ,L n denote the total Virasoro generators.
In order to express the contribution from the Feynman diagrams in the usual form of the string amplitude in the first quantized formalism, we need to express (1.1) as We use the convention that ds ∧ dθ describes positive integration measure. The collection of the variables (s, θ) for different propagators, together with the integration variables that enter in the definition of the interaction vertices, become the coordinates of the vertex operators on the world-sheet after some change of variables.
Note however that the equality of (1.1) and (1.2) holds only when both sides act on states with positive L 0 +L 0 eigenvalues. Even tachyon free superstring field theories have off-shell internal states with negative L 0 +L 0 eigenvalues since they carry Lorentzian space-time momentum. Acting on such states (1.1) is perfectly well-defined, but (1.2) is divergent. It is these divergences that show up as apparent divergences in the integration over the locations of the vertex operators (punctures) in the first quantized formalism.
String field theory suggests different (but equivalent) ways of resolving this problem. The first approach will be to directly use the expressions for amplitudes given as sum over Feynman diagrams without using (1.2). It is possible to formulate string field theory in a way that the sum over intermediate states converges rapidly even though the sum runs over infinite number of states [2]. The second approach will be to rewrite (1.1) by: where Λ is a large number and in the last term the integration over q runs in the anti-clockwise direction. In the second expression in (1.3) we have used the fact that in the second term the integration over q projects to L 0 =L 0 states, allowing us to replace the L 0 +L 0 factor in the denominator by 2L 0 . The total expression is independent of Λ and is equal to (1.1) both for positive and negative values of L 0 +L 0 . However the advantage of this representation is that for large Λ the contribution from the first term can be represented as the usual string amplitude with integration range over the world-sheet variables having certain excluded domains corresponding to the |q| < e −Λ regions, while the second term inside the square bracket receives appreciable contribution only from finite number of states with L 0 +L 0 ≤ 0. 2 A third representation of the propagator (1.1) is: where ǫ is a small positive parameter which we eventually take to 0. This representation of the propagator was used in [22] to get finite results for string amplitudes. This is equivalent to the 2 As discussed in the last paragraph of §3.1, for on-shell amplitudes the organization of various terms based on the representation (1.3) can be reinterpreted as the one based on the representation (1.1) by redefining the interactions vertices of string field theory. This can be achieved by absorbing the first term inside the square bracket in (1.3) and the e −Λ(L0+L0) factor from the second term into the definition of the interaction vertices, leaving only the propagator factor proportional to (L 0 +L 0 ) −1 . In the string field theory literature, this operation is known as 'adding stubs' [25].
Feynman iǫ prescription, but applied to tree level diagrams. In order to use (1.4) for numerical evaluation of the amplitude, we must work with a finite ǫ since for ǫ = 0 the integrand is oscillatory and not damped. The actual answer is then obtained by evaluating this expression for various values of ǫ and then extrapolating the result to ǫ → 0.
Of course in all these approaches, we have a genuine divergence when an intermediate state has exactly vanishing L 0 +L 0 . This represents the physical poles of the tree level S-matrix, and happens when we choose the external momenta such that an intermediate state goes on-shell.
In the computation of on-shell amplitudes in superstring theory, the apparent divergence arising from intermediate states with negative (L 0 +L 0 ) eigenvalue can some time be enhanced due to the use of vertex operators of wrong picture number. String field theory gives a definite prescription for the choice of locations of PCOs when two or more vertex operators come close.
This in particular guarantees that only −1 and −3/2 picture states propagate as intermediate states. If we use the wrong picture number, e.g. by converting some of the NS sector vertex operators into zero picture, states in 0 picture propagate in the intermediate state when a pair of zero picture vertex operators approach each other. This sector has states with negative conformal weights leading to apparent divergences which are not present in the correct treatment of string field theory.
In this paper we show that using the representation (1.3) of the propagator (equivalently representation (1.1) following footnote 2) and the fact that the change in the locations of PCOs results in addition of total derivative terms to the integrand, one can arrive at the following prescription for computing (super-)string tree amplitudes that avoids all divergences: 1. Let us suppose that for an (n + 3)-point amplitude of an arbitrary set of vertex operators in an arbitrary (compactified) string theory, I (0) (σ 1 , · · · , σ n ) is the integrand obtained by computing the relevant correlation functions of (n + 3) vertex operators on the sphere, with three of the vertex operators at fixed positions and the n vertex operators inserted at σ 1 , · · · , σ n . Naively the amplitude is given by the integral of the 2n-form I (0) over the whole moduli space parametrized by σ i 's, but this integral diverges from the regions where the vertex operators come close, or equivalently, the sphere degenerates. labelling the moduli of the pair of punctured spheres into which the original sphere degenerates.
3. Multiple degenerations correspond to the region of the moduli space where the original sphere degenerates into several spheres. A degeneration into (k + 1) spheres is characterized uniquely by a set of k unordered labels s 1 , · · · , s k where each s i corresponds to one of the labels characterizing single degenerations. Near such degenerations we can construct from the {σ i }'s a set of complex coordinates u s 1 , · · · , u s k that vanish at the degeneration and another set of n − k complex coordinates m (s 1 ,···,s k ) , labelling the moduli of the punctured spheres to which the original sphere degenerates.
4. We now denote by D s an open tubular neighborhood of the subspace u s = 0 and define C (0) to be the region of the moduli space that excludes D s for all s. 3 Therefore by construction C (0) excludes all degenerations. We also denote by C s 1 ···s k is fixed by the equation where the sum over s runs over all labels for which C (k+1) s 1 ···s k s exists.
5. Then the full amplitude (for generic external momenta) is given by 3 For example in a local patch we could take |u s | < ǫ for some small number ǫ as the definition of D s . The choice of u s 's is of course not unique, but the final result is independent of this choice. Furthermore the coordinates {u s1 , · · · , u s k } may need to be defined differently in different local patches of the space spanned by m (s1,···,s k ) . In this case on the overlap one needs to choose appropriate interpolation in the definition of D s between the two patches. An explicit example elaborating the choice of the coordinates {u s }, {m (s1,···,s k ) } and the definition of D s for a five punctured sphere has been described in appendix A.
where the (2n − k) forms I (k) s 1 ···s k , defined in a neighborhood of C (k) s 1 ···s k , are obtained by solving the equations Beginning with the known expression for I (0) , we could solve for I (k) s 1 ···s k iteratively in k by expanding both sides of (1.7) in a power series in u s 1 , · · · , u s k and their complex conjugates and matching the coefficients. 6. The solution to (1.7) is not unique but the final result is not affected by this ambiguity, as long as we work with generic external momenta so that the exponents in the power series expansion in u s ,ū s are not integers. Similarly the choice of C Note that for k = 0 term in (1.6), the integrand is the same as the original integrand but the integration runs over the region C (0) that avoids all degenerations. The missing regions are compensated for by the remaining integrals involving k ≥ 1 terms. Also note that the final formula for the amplitude encoded in (1.6) and (1.7) does not require any detailed knowledge of string field theory even though we use string field theory to arrive at these formulae.
We can regard C (k) s 1 ···s k I (k) s 1 ···s k as the result of using the second term in (1.3) -or equivalently (1.4) -for the propagators labelled by s 1 , · · · , s k if we identify the variable u s with the variable q for the s-th propagator and the tubular neighborhood D s as the region |u s | < e −Λ . As we shall see, this choice of {u s } and {D s } is not strictly necessary. Nevertheless this makes the connection to the iǫ prescription of [22] clear: C (k) s 1 ···s k I (k) s 1 ···s k represents part of the contribution where the integrals over ln |u s | runs parallel to the imaginary axis for s = s 1 , · · · , s k . We note however that the prescription given in (1.6) works equally well for complex external momenta, while the iǫ prescription will typically require different choices of the contour of integration in the complexified moduli space for different complex external momenta to ensure convergence. The procedure described above works even if the original integrand I (0) is computed with wrong choice of PCO locations, e.g. by taking the PCO locations to coincide with some of the vertex operators. This adds a total derivative to I (0) and could introduce additional divergences when the vertex operators collide. However the addition of such terms to I (0) also requires us to add corrections to I (1) s when we solve (1.7), and these cancel in (1.6). Therefore while implementing this procedure, one does not need to be careful about the choice of PCO locations.
As already mentioned, during the analysis described above we do not assume any particular choice of string background or any specific form of the external states. Our analysis holds for generic external states in a generic (compactified) string theory.
The second application of classical string field theory that we discuss is in the study of marginal deformations of the world-sheet CFT describing the background space-time. String field theory is formulated around some specific classical solution, corresponding to a CFT of matter and ghost fields on the world-sheet. If this CFT has marginal deformations, then one can in principle compute the spectrum and correlations functions of all the operators in the deformed theory using conformal perturbation theory. This in turn will determine the spectrum of physical states and the S-matrix of string theory formulated in the deformed background, in terms of the spectrum and correlation functions of the original CFT. In practice this requires carefully regulating the ultraviolet divergences on the world-sheet since the marginal operator needs to be integrated on the world-sheet together with other vertex operators and we have to regulate the divergences when the locations of the operators come close to each other.
In particular at higher order in perturbation theory when there are several insertions of the marginal operator we need to carefully remove all the divergences, while making sure not to remove the finite pieces. Now given the deformed CFT, one can formulate a string field theory around this deformed CFT. It is known that this new string field theory is related to the string field theory formulated around the original CFT by a field redefinition that includes a shift [26][27][28]. This means that the original string field theory has a classical solution that describes string field theory around the deformed background. 4 Furthermore the spectrum and S-matrix of the string field theory around the deformed CFT can be computed by expanding the original string field theory action around the classical solution describing the deformed background. Since sting field theory has no divergences, this procedure also does not suffer from any divergences. In other worlds string field theory automatically provides an ultraviolet regulator for the worldsheet theory. A fully systematic procedure for constructing the classical solution in string field theory to any given order in the expansion in terms of the deformation parameters, and computing the spectrum and S-matrix of string theory in the deformed background, can be found in [2,34]. Its application to a particular class of examples has also been described in [35] (see also [36,37] for construction of the solution). The apparent ultraviolet divergences in the conformal perturbation theory can be traced either to wrong use of Schwinger parametrization as in (1.2), or applying (L 0 +L 0 ) −1 on states that have zero (L 0 +L 0 ) eigenvalue. The first case is dealt with by using (1.1) or (1.3), while the second case requires a more elaborate treatment as explained in sections 3.3 and 4.1 of [34].
We apply this method to the special case of string compactification on a circle, and consider the marginal deformation to be the one that changes the radius of the circle. The mass of a state carrying momentum n/R along the circle gets shifted under this deformation. We find that both for the bosonic and the heterotic string theory the mass 2 shift to second order in the deformation parameter µ is given by Since the leading order expression for m 2 has an additive term given by n 2 R −2 , (1.8) corresponds to a scaling of R −2 by (1 − µ + µ 2 /2). This is consistent since the marginal deformation is expected to change the radius. If ∆m 2 had not been proportional to n 2 then this would not be possible. Indeed at the intermediate stages of the calculation there are other terms proportional to n 4 R −4 and n 0 R 0 , but they cancel at the end. Nowhere at any stage of the calculation we need an ultraviolet regulator on the world-sheet. Furthermore the method we use is completely general and can be used to compute the shift in the spectrum under a general marginal deformation where the answer may not be a priori known, e.g. shift in the masses of heavy string states under a blowup of the orbifold singularity.
Another possible application of string field theory is in the study of Ramond-Ramond background [38]. Once a string field theory is formulated around a given world-sheet CFT, one may have a family of consistent classical solutions that involve switching on Ramond-Ramond background. In this case one can use the same procedure used for marginal deformations to systematically construct the classical solution in string field theory to any order in deformation parameter, and compute the spectrum and S-matrix of the theory around the deformed background [38]. As in the case of marginal deformations, this procedure never requires an ultraviolet regulator on the world-sheet if we apply the formalism of string field theory systematically.
As is well known, formulation of string field theory requires us to specify certain datathe choice of local coordinates at the punctures of the Riemann surface and locations of the PCOs. The results we have quoted are independent of the data. Indeed, one of the goals in our analysis will be to manipulate the various expressions in such a way that even though at the intermediate stages of the calculation the result depends on these additional data (which we leave unspecified), this dependence cancels at the end. This is manifest in our expressions (1.6) and (1.7). One finds for example that different choice of local coordinate system at the punctures lead to different choices of the subspaces C s 's. 5 Similarly different choice of PCO locations can be shown to induce an additive total derivative term in the expression for I (0) , but (1.6) can be shown to remain unchanged under addition of such a term. In essence, such an additive term in I (0) forces us to also add certain terms in I (1) s in order to satisfy (1.7). These extra terms cancel the effect of additive term in I (0) when we evaluate (1.6).
Similarly we see that (1.8) does not depend on the choice of local coordinate systems or PCO locations. At intermediate stages of the analysis various expressions we get do depend on the additional data, but this dependence cancels at the end. We must note however that the total independence of the result (1.8) of the additional data is accidental since different choice of local coordinates and PCO locations lead to string field theories that are related by field redefinition and we expect that under such a field redefinition the deformation parameter µ -which is the component of a field -will also get redefined. Therefore at higher order in the expansion in powers of µ, we do expect the result to depend on the additional data involved in the construction of string field theory, but this dependence should be removable by a redefinition of µ. The analysis leading to (1.8) is close to the spirit of the analysis in [38] in that the latter paper computed the shift in mass under deformation involving RR background to second order in the deformation parameter µ. There is however one important technical difference between the two analysis. In the case of RR background, the leading contribution to the mass shift appears at second order in µ. Therefore for computing the mass shift, the vertex operators of the states whose mass shift is being calculated could be taken to be dimension zero primaries. For the deformation we consider in this paper, the first correction to the mass already appears at order µ, and therefore to compute the mass at second order in µ we need to use vertex operators whose dimensions differ from zero by order µ. This introduces non-trivial dependence on the choice of local coordinate system at intermediate steps of the calculation even though the dependence cancels at the end. This difference is similar to the difference between the calculation of one loop mass renormalization and two loop mass renormalization. The former requires computing the torus two point function of a pair of vertex operators that satisfy tree level on-shell condition. However the latter requires use of vertex operators that take into account one loop mass renormalization, and therefore do not satisfy the tree level on-shell condition.
The rest of the paper is organized as follows. In §2 we review some conventions we use to compute correlation functions in the world-sheet CFT and briefly review some aspects of string field theory that we need for our analysis. In §3 we discuss how choice of local coordinate system and PCO locations affect the results for Feynman diagrams in string field theory. In §4 we make use of the representation (1.3) of the propagator to give a manifestly finite expression for the four tachyon amplitude in bosonic string theory. In §5 we use (1.3) to give a general algorithm for getting manifestly finite expressions for general amplitudes in bosonic string theory involving arbitrary number of external states carrying arbitrary quantum numbers. Our analysis in this section leads to eqs.(1.6), (1.7). In §6 we use superstring field theory to arrive at the same formulae (1.6), (1.7), giving manifestly finite expressions for the superstring tree amplitudes. In particular we show that (1.6), (1.7) give the correct amplitude even if I (0) is computed using wrong choice of PCO locations. In §7 we apply string field theory to study the effect of marginal deformations in bosonic and heterotic string theory and arrive at (1.8) in both theories. We conclude in §8 with some comments on higher genus amplitudes. Appendix A contains examples of the choice of the coordinate {u s }, {m (s 1 ···s k ) } for a five punctured sphere. Appendix B contains some technical results needed to complete the analysis in §7.2.

Conventions
In this section we shall briefly review some aspects of the world-sheet theory and the string field theory that we shall need for our analysis . More details can be found in [2]. We begin by describing the normalization conventions for the vacuum of the world-sheet (super-)conformal field theory of the matter ghost system. For bosonic string theory the ghost system has b, c,b andc ghosts with the usual mode expansion. We normalize the SL(2,C) invariant vacuum of the bosonic string as up to a factor given by the overall volume of space-time which eventually generates the momentum conserving delta function in a correlator. We shall not write this factor explicitly. For heterotic string theory, besides the b, c,b,c ghost fields we also have the β, γ ghosts, related to the ξ, η, φ system via the relations We choose the normalization of the vacuum of the heterotic string such that For type II string theories we also have anti-holomorphicβ,γ system and the normalization of the vacuum will be chosen as We shall denote by X µ the world-sheet fields corresponding to non-compact space time coordinates, and, for the analysis in §7, by Y the world-sheet field corresponding to a compact space direction of radius R. We also denote by ψ µ and χ their holomorphic superpartners on the world-sheet -in type II theories we also have anti-holomorphic fieldsψ µ andχ. Their operator product expansions have the form: where · · · denote less singular terms whose knowledge will not be needed for our analysis. There may also be additional component of the matter CFT describing other compact directions that will not be relevant for our analysis. The operator product expansion of the ghost fields take the form e q 1φ (z) e q 2φ (w) = (z −w) −q 1 q 2 e (q 1 +q 2 )φ(w) + · · · ,∂φ(z)∂φ(w) = − 1 (z −w) 2 + · · · ,(2.6) where · · · denote less singular terms.
In the heterotic string theory we have holomorphic PCO given by where Q B is the BRST charge: T m ,T m denote components of the matter stress tensor and T β,γ is the stress tensor of the β, γ system. z includes the 1/2πi factor for holomorphic integral and −1/2πi factor for the anti-holomorphic integral. T F is the super-stress tensor of the matter SCFT, given by Here (T F ) int denotes the contribution from the additional compact target space directions other than the Y -χ system. In the bosonic string theory the contribution from the β, γ system will be absent, while in type II theory there will be additional contribution involvingβ,γ. In type II theory we also have the anti-holomorphic PCO, obtained by replacing the holomorphic fields by anti-holomorphic fields in (2.7).
In the bosonic string theory the physical unintegrated vertex operators take the form ccV where V is a dimension (1,1) primary in the matter CFT. From this one can construct the integrated vertex operator: where z denotes a contour around z.
In the heterotic string theory the unintegrated −1 picture NS sector vertex operator takes the form cce −φ V , (2.11) where V is a dimension (1, 1/2) superconformal primary in the matter SCFT. The unintegrated zero picture vertex operator takes the form where W is a dimension (1,1) matter sector vertex operator defined via The integrated −1 picture vertex operator takes the form  14) and the integrated 0 picture vertex operator is given by  In type II string theory we can similarly define integrated and unintegrated NSNS sector vertex operators carrying picture numbers (−1, −1), (0, −1), (−1, 0) and (0, 0). We shall not write down the explicit form of the Ramond sector vertex operators since they will not be needed for our analysis. Tree level (n + 3)-point amplitude of vertex operators ccV i for 1 ≤ i ≤ (n + 3) in bosonic string theory is given by converting n of them to integrated vertex operators and integrating the resulting correlation function over the locations of the integrated vertex operators: (2.16) Here y n+1 , y n+2 and y n+3 represent arbitrary points in the complex plane. The (−1/2πi) n factor is a normalization factor that appears in the definition of the interaction vertices of string field theory and is related to the −1/2πi factor in (1.2). A given in (2.16) can be regarded as a contribution to the T -matrix, related to the S-matrix by S = 1 + i T .
The results for the tree level (n+3) point function of NS sector states in the heterotic string theory and NSNS sector states in type II string theories are similar, except that the vertex operators are taken in the −1 picture and the amplitude has insertion of (n + 1) of PCOs in the heterotic theory and (n + 1) holomorphic PCO's and (n + 1) pair of anti-holomorphic PCO's in type II theory. The PCO locations are arbitrary when the vertex operators are well separated from each other but need to satisfy certain relations when two or more vertex operators approach each other. These rules are induced from superstring field theory and essentially tell us that when n vertex operators approach each other, we must also have (n − 1) PCOs approaching them so that the picture numbers of all the operators add up to −1. If there are R sector vertex operators present, then they are taken in the −1/2 picture. The number of PCO's need to be adjusted so that the total picture number of all vertex operators and PCOs add up to −2 in the heterotic string theory and (−2, −2) in the type II string theory.
We shall now briefly review some aspects of (super-)string field theory. In bosonic string theory we denote by H the Hilbert space of matter-ghost CFT satisfying the conditions: For NS sector of heterotic string theory and NSNS sector of type II string theory we have similar constraints except that we also require the states in H to carry picture number −1 in heterotic string theory and picture number (−1, −1) in type II string theory. In all the theories the string field |Ψ is an arbitrary element of H. We shall denote by Ψ the corresponding vertex operator in the CFT.
The classical action of bosonic string field theory and the NS sector fields in superstring field theory takes the form: 6 where for |A i ∈ H, {A 1 · · · A n } is a multilinear function of the A i 's obtained by first constructing the sphere correlation functions of certain ghost operators and PCOs and the vertex operators A i in specified coordinate system, and then integrating the result over an appropriate subspace of the moduli space that excludes all the singular regions where two or more vertex operators come close. The precise choice of these subspaces, or the coordinate system in which the vertex operators are inserted, or the PCO locations are not fixed completely but are subject to stringent constraints, and different choices lead to different string field theories which are related by field redefinition. Note that in this definition we do not require the A i 's to be BRST invariant.
We also define [A 2 · · · A n ] ∈ H such that In terms of this, the equations of motion of the string field |Ψ , derived from the action (2.18), takes the form: Inclusion of the Ramond (R) sector states requires additional structure. In heterotic string theory we introduce a pair of string fields Ψ R and Ψ R belonging to −1/2 and −3/2 picture number sectors satisfying (2.17). In superstring theory Ψ R will include fields in the (−1, −1/2), (−1/2, −1) and (−1/2, −1/2) sectors while Ψ R will include fields in the (−1, −3/2), (−3/2, −1) and (−3/2, −3/2) sectors. The kinetic term takes the form In heterotic string theory G is the zero mode of the PCO. In type II string theory G is the zero mode of the holomorphic (anti-holomorphic) PCOs in the NSR (RNS) sectors and product of zero modes of holomorphic and anti-holomorphic PCOs in the RR sector. The interaction term involving the NS and R sector fields have the form In particular the field Ψ R does not enter the interaction term. Therefore the Feynman rules only require the Ψ − Ψ propagator. In Siegel gauge this has the form similar to (1.1) except for an extra factor of G: The relevant equations of motion in the Ramond sector, obtained by taking a linear combination of the equations of motion of Ψ R and Ψ R is given by The other linear combination satisfies free field equations of motion and decouples from the theory. Of course we must also replace Ψ by Ψ + Ψ R in the second term in (2.20).

Local coordinates and picture changing operators
As mentioned earlier, for the construction of a string field theory action, we need to define the interaction vertices involving off-shell string fields. For n-point interaction vertex this requires identifying certain subspaces of the moduli space of a sphere with n punctures, and a choice of local coordinates at the punctures, satisfying certain consistency conditions. For superstring theory we also need to specify choice of the locations of the picture changing operators. In this section we shall describe some general criteria for choosing the local coordinate system at the punctures of a three punctured sphere involved in the definition of 3-string interaction vertex and the corresponding choice of a subspace of the moduli space of a sphere with 4 punctures, involved in the definition of the 4-point interaction vertex. We shall also discuss the general criteria for the choice of locations of the PCOs. However we shall refrain from committing ourselves to any particular choice of local coordinates or PCO locations since one of our goals will be to demonstrate that the final result is independent of these choices.

Local coordinate system
On a sphere with three punctures, we can use SL(2,C) transformation to choose the locations of the punctures to be at We shall choose the local coordinates w i at these punctures to be related to the global coordinate z by up to arbitrary phases. h i (z) satisfies h i (z i ) = 0 for 1 ≤ i ≤ 3, so that w i vanishes at z = z i . λ is an arbitrary positive real number that could have been included in the definitions of h i (z), but we have displayed it explicitly since at the end we shall try to simplify our analysis by taking the large λ limit. h i (z) is an analytic map between an open neighborhood around z = z i and the unit disk |w i | < 1, but may have singularities outside this domain. We require the images of the |w i | ≤ 1 regions in the z-plane to be non-overlapping -this can be achieved by taking λ sufficiently large. Furthermore, up to overall phases, this choice of local coordinates should be invariant under permutation of the punctures. For example, since the transformation z → 1/z exchanges the punctures 1 and 2 and leaves z 3 invariant, we must have up to phases. Similarly we should have and up to phases. 7 Let us now consider another sphere with three punctures, carrying global coordinate z ′ and the punctures situated at carrying local coordinates around the three punctures. We can construct a two parameter family of spheres with four punctures by sewing the two spheres at their third punctures via the relation This family of four punctured spheres is what we shall obtain from the s-channel Feynman diagram of string field theory. Using (3.2) and (3.7) we get the relation between the global coordinates z and z ′ on the two spheres: inside open neighborhoods of z = z 3 and z ′ = z ′ 3 . Since the glued Riemann surface is a sphere with four punctures, we can introduce a global coordinate y on the four punctured sphere. If the h i 's had been SL(2,C) transformations as in footnote 7 then y could be taken to be either z or z ′ or related to these by an SL(2,C) transformation, but in general the relation between y and the original coordinates z, z ′ is more complicated. By an SL(2,C) transformation we can ensure that in the y plane the original punctures at z = 0 and ∞ are located at y = 0 and ∞, and the puncture at z ′ = ∞ is located at y = 1. We shall denote by σ the location of the puncture at z ′ = 0 in the y plane. It takes the form of a function of q/λ 2 due to (3.9). Therefore in the y plane the four punctures are located at: for some function g. y 4 = σ is a holomorphic function of q for |q| < 1, but we shall not make use of this information in an essential way. We also define v 1 , v 2 , v 3 and v 4 to be the original local coordinates w 1 , w 2 , w ′ 1 and w ′ 2 , but now expressed as function of y. These take the form for some functionsĥ i . The overall multiplicative factors of λ comes from the multiplicative factors of λ in the definitions of w 1 , w 2 , w ′ 1 and w ′ 2 while the dependence on q/λ 2 enters through (3.9).
It follows from (3.9) that as q → 0, the four punctured sphere degenerates into two spheres, with the original punctures at z = 0 and z = ∞ on one sphere and the original punctures at z ′ = 0 and z ′ = ∞ on the other sphere. In the y plane this will correspond to the punctures at σ and 1 coming together. To find the behavior of g(q/λ 2 ) for small q, we note that for small q/λ 2 we can take z, z ′ close to 1 and express the relation (3.9) as Since in this case the local coordinates are related to global coordinates z and z ′ via S(2,C) transformation, we can identify z (and z ′ ) with y up to SL(2,C) transformation. Now under the identification (3.12) the puncture at z ′ = ∞ is mapped to z = 1 and the puncture z ′ = 0 is mapped to z = 1 − q λ −2 (h ′ 3 (1)) −2 . Therefore we can make the identification y = z and Comparing this with (3.10) we get (3.14) The region |q| ≤ 1 corresponds to a neighborhood of the point 1 in the σ plane. In string field theory, if we use the local coordinates given in (3.2) to define the 3-point interaction vertex, then the family of four punctured spheres (3.10) corresponding to |q| ≤ 1, and the local coordinates given in (3.11), describe the contribution to the four point Green's function due to 's-channel diagrams'. In this the external states represented by the vertex operators inserted at the punctures at y 1 and y 2 merge to form an intermediate state which then splits into the states represented by the vertex operators inserted at the punctures y 3 and y 4 .
Contribution from the u-channel diagram is obtained by exchanging 1 and 3. Denoting the global coordinate on the plane byỹ, we get, by exchanging 1 and 3 in (3.10), 15) and the local coordinates at the punctures are We now make a change of variable so that in the y coordinate system the locations of the punctures are given by 18) and the local coordinates around the punctures are given by, (3.19) Contribution from the t-channel diagram is given by exchanging 2 and 3 in the s-channel contribution. Denoting the global coordinate on the plane byŷ, we get We make a change of variables In this coordinate system the punctures are located at 23) and the local coordinates are given by Using (3.13) we see that in the y 4 = σ plane, the region |q| ≤ 1 in (3.10), (3.18) and (3.23) map respectively to some regions around 1, 0 and ∞ which we shall denote by R s , R u and R t . Integrals over these regions describe the contributions from s, u and t-channel Feynman diagrams. The rest of the region in the σ plane, which we shall call R, must come from the elementary four point interaction vertex.
In order to define the four point interaction vertex for off-shell string states, we need to look for local coordinates of the form satisfying the following conditions: 2.h i (y; σ) must match the results (3.11), (3.19) and (3.24) of s, u and t channel diagrams when σ takes value at the common boundary of R and R s , R and R u and R and R t respectively, corresponding to setting q = e −iθ in (3.15), (3.18) and (3.23).
3. v i (σ) must obey permutation symmetry. This means the following. Let us suppose that we exchange the punctures k and ℓ for some fixed k and ℓ. This exchanges y k with y ℓ and also v k with v ℓ . Now by an SL(2,C) transformation we can bring the locations of the punctures 1, 2 and 3 to their original values 0, ∞ and 1, but in that process y 4 = σ will typically map to a different point σ ′ and the local coordinates v i , expressed in terms of the new complex coordinate, which we shall still denote by y, will change to v ′ i . The requirement that we would like to impose is that v ′ i should agree with v i at σ ′ . This will guarantee that interaction vertex is invariant under permutation of the external states.
While it will be nice to have an interaction vertex of this type, the third condition is not strictly necessary. Given a choice of local coordinates that satisfy conditions 1 and 2 but not 3, we can define the interaction vertex by taking average of the result computed from this choice of local coordinates and the images of these local coordinates under permutation. Similar procedure can be followed for the construction of the local coordinates needed for higher point interaction vertices of string field theory.
We shall end this section with a few comments on the parameter λ.
1. While we expect the on-shell amplitudes to be independent of the parameter λ (and the functions h i ), the off-shell amplitudes depend on λ. Since dw i /dz computed from Therefore for large λ the contribution from the modes with large L 0 +L 0 eigenvalue is highly suppressed. This suppression factor in turn makes the sum over massive states in the intermediate state finite even though the number of such modes has a Hagedorn growth.
2. As can be seen from (3.10), (3.18) and (3.23), for large λ the s, u and t channel diagrams cover small regions of the σ ≡ y 4 plane around 1, 0 and ∞ respectively, given by the |q| ≤ 1 regions in (3.10), (3.18) and (3.23). Therefore most of the integration over σ = y 4 is generated by the four point interaction vertex. This is related to the fact that the 2) cover small regions in the z plane.
3. Note that q always appears in the combination q/λ 2 . This means in particular that if we change the cut-off Λ in (1.3) and simultaneously change λ so that λ 2 e Λ remains fixed, then the cut-off |q| = e −s ≥ e −Λ remains unchanged in the σ = y 4 plane. Using this we can set Λ = 0 by scaling λ by e Λ/2 . In this case the first term in (1.3) vanishes and the second term becomes (1.1). Therefore use of (1.3) to deal with the negative (L 0 +L 0 ) states is equivalent to the use of (1.1), with a rescaled value of λ.

Locations of picture changing operators
For constructing the interaction vertices of superstring theories, we also need to make a choice of PCO locations. We shall now briefly describe the procedure. For simplicity we focus on the heterotic string theory. For type II string theories we have to repeat the procedure for the left and the right-moving sectors of the world-sheet theory separately.
We begin with the cubic interaction vertex. There are two types of vertices: R-R-NS and NS-NS-NS. For the R-R-NS vertex the total picture number of the vertex operators is −2 and we do not need to add any PCO. Therefore the vertex is the same as that in bosonic string theory. For the NS-NS-NS vertex, the picture numbers of the three NS-sector vertex operators add up to −3 and we need to insert one PCO. The location of the PCO needs to be invariant under permutation symmetry, i.e. under the SL(2,C) transformations that exchange the positions 0, 1 and ∞ of the vertex operators. One can either try to choose the location to be invariant under this symmetry group, or pick an arbitrary location and average over all its images. It is easy to verify that there is no single location that is invariant under the subgroup of SL(2,C) generating arbitrary permutation of 0, 1 and ∞. The brute force way will be to choose an arbitrary point p 1 and its images p 2 , · · · p 6 under the symmetry generated by z → 1 − z and z → z/(z − 1) and take the PCO insertion to be (3.27) A more economical procedure will be to pick a pair of points, each of which is invariant under the subgroup that generates cyclic permutation of 0, 1 and ∞, and the pair gets exchanged under a Z 2 transformations that exchange 0 and 1, leaving ∞ fixed. These pair of points are: We can now insert (X (p 1 ) + X (p 2 ))/2 to get a symmetric 3-point interaction vertex. In the following we shall denote the PCO location for a 3-point interaction vertex by p, keeping in mind that p may stand for averages over several locations.
Let us now consider the contribution to the 4-point amplitudes associated with s, t and u channel diagrams. First consider the case of four R-sector external states. In this case the intermediate state in each channel is an NS sector state. No PCO is needed at any stage of this calculation, reducing the analysis to that of bosonic string field theory.
Next consider the case of four NS sector external states. In the s-channel diagrams the PCOs are located at z = p, z ′ = p , (3.29) where p corresponds to average of several insertions as in (3.27). Their images W 1 and W 2 in the y plane may be found using the known map between z, z ′ and y. The information that will be useful to us later is that in the small q or large λ limit, one of the PCO's (say W 1 ) is finite distance away from the punctures 1 and σ while the other one (say W 2 ) is close to the punctures 1 and σ. For t and u channel diagrams the locations of the PCOs are related to W 1 , W 2 by the transformations y → y/(y − 1) and y → (1 − y) respectively. While constructing the 4-point interaction vertex of four NS sector states, one can choose the local coordinates at the punctures for σ ∈ R as in the case of bosonic string theory, but we also need to fix the location of two PCOs as a function of the modulus σ. The PCO locations (W 1 , W 2 ) on ∂R s , ∂R u and ∂R t induced from the s, u and t-channel diagrams provide appropriate boundary condition that the PCO locations must satisfy. They also must respect the permutation symmetry. This can be achieved by brute force by choosing one configuration satisfying the boundary condition and then averaging over all its images under permutation. However, one may be able to reduce the number of configurations to average over by making judicious choice of the initial configuration so that it is invariant under a subgroup of the permutation group. Other than these constraints the PCO locations on the four punctured sphere can be chosen arbitrarily for σ ∈ R. If necessary (e.g. to avoid collision with each other) we can even allow the PCO locations to jump discontinuously across codimension one subspaces inside R and add correction terms [44,45], so there is no obstruction to choosing the PCOs satisfying the desired conditions.
The case of two NS and two R vertex operators can be analyzed similarly. In this case in the channel where an R sector state propagates as intermediate state, the PCO insertion comes from the factor of G in the R sector propagator (2.23). Therefore the PCO insertion involves an average not over discrete number of possibilities but a continuous set of possibilities. As before, for σ ∈ R the PCO locations can be chosen arbitrarily subject to the boundary conditions on ∂R s , ∂R t and ∂R u and the symmetry requirement (which in this case corresponds to the exchange of the two NS punctures and (independently) the two R punctures.

Warm up with four point amplitude
We shall now show how bosonic string field theory can be used to compute the four point amplitude over the full range of external momenta without any need for analytic continuation. In §4.1 we shall outline the general procedure for determining the boundary term that arises when we use (1.3) to deal with the contribution from states with negative L 0 +L 0 eigenvalue. In §4.2 we shall use this to express the four tachyon amplitude as an integral over the world-sheet that does not have any divergence.

A concise description of the boundary terms
Let us consider either the s, t or u-channel Feynman diagram of string field theory four point amplitude. Each of these diagrams has a single propagator. Using the representation (1.3) of the propagator we can express the contribution to the diagram from the first term on the right hand side of (1.3) as a two dimensional integral over the complex variable q, with the integrand given by where γ i and δ i are the L 0 andL 0 eigenvalues of the i-th internal state and the constants A i are given by the products of the vertices of the Feynman diagram. We shall work with generic external momenta so that there are no terms with integer γ i , δ i . On the other hand the contribution to the Feynman diagram from the second term on the right hand side of (1.3) may be represented as a boundary integral of the form: 2) The choice (4.3) of I (1) is not unique since we could have taken this to be proportional to dq or a linear combination of dq and dq which all satisfy (4.4). This ambiguity corresponds to the freedom of adding exact forms to I (1) .
By the standard procedure in string field theory illustrated in §3.1 one can convert the integration over the parameters q,q associated with the propagators to the integration over the standard world-sheet variables σ,σ denoting the world-sheet coordinate of one of the vertex operators, keeping fixed the positions of the other three vertex operators. In the convention of §3.1 the region |q| ≤ 1 gets mapped to the regions R s , R t and R u in the σ-plane for the s, t and u-channel diagrams. Let us suppose that the region |q| ≤ e −Λ gets mapped to the regions R s , R t and R u for these diagrams. Then the boundary integrals (4.2) will run over ∂ R s , ∂ R t and ∂ R u . On the other hand the bulk integral, after combining the contributions from the s, t and u-channel Feynman diagrams and the diagram involving the four point interaction vertex, will run over the full complex σ-plane except the excluded regions R s , R t and R u .
(4.4) gives a simple way of determining I (1) from the original integrand I (0) without knowing the relation between q and σ. Let us suppose that σ approaches σ 0 as q approaches 0. Here σ 0 is the location of one of the other vertex operators. Then near q = 0 the relation between σ and q takes the form In terms of the variables σ andσ we can express (4.1) as where B j are new coefficients and α j , β j take values that differ from γ j , δ j at most by integers.
We can now solve (4.4) to get up to closed differential forms. Now one can show that a closed differential form with an expansion in non-integer powers of (σ − σ 0 ) and (σ −σ 0 ) is also exact, with an expansion involving non-integer powers of (σ − σ 0 ) and (σ −σ 0 ). Therefore I (1) given in (4.7) differs from the one in (4.3) at most by an exact form and its integral over the cycle |q| = e −Λ gives the same result as that of (4.3). When σ 0 = ∞, one has to use a slight variant of this procedure, with In this procedure the only part that requires the knowledge of the relation between q and σ is the curve C = ∂ R u , ∂ R s or ∂ R u -the image of the curve |q| = e −Λ in the σ plane. However we shall now show that the sum of the bulk and the boundary terms is independent of the form of C. For this let us recall that the bulk integration runs over the region outside the curve C in the σ plane. Therefore if we change C to C ′ , the bulk integration region will change by a region D with ∂D = C − C ′ . The net change in the bulk and the boundary integral is given by This shows that even the knowledge of C is not needed and (4.7) gives us complete information about the boundary term in terms of the bulk integrand (4.6). This in turn shows that in order to evaluate the sum of the boundary terms and the bulk integral that defines the amplitude, we do not need any knowledge of the local coordinates used in defining the interaction vertices of string field theory. While for computing amplitudes with generic external momenta, we can assume that γ j , δ j and hence α j , β j are not integers, in §7, where we analyze the effect of marginal deformations, we shall encounter situations where α j , β j are integers. In this case there may be an additional contribution to the boundary term involving the integral of a one form that is closed but not exact: Here K is a constant. As we shall see in §7, possible additional contributions from such terms can be calculated by working directly with (4.1).

Four tachyon amplitude
We shall now apply the procedure described in §4.1 to give an expression for the four tachyon amplitude as a world-sheet integral that does not have any divergence. First let us describe the conventional world-sheet computation in the α ′ = 1 unit. Let us denote by k 1 , k 2 , k 3 and k 4 the momenta of the tachyons, all counted as positive if ingoing, and satisfying the on-shell condition where m 2 is the mass 2 of the tachyon. Then according to (2.16), the amplitude takes the form where y 1 , y 2 , y 3 are arbitrary fixed points, which we can take to be 0, ∞ and 1 respectively, After evaluating the correlation function using (2.1) the result takes the form: up to overall momentum conserving delta function. Defining (4.14) we can express (4.13) as This integral has divergences from the region σ → 0, σ → 1 and σ → ∞ if, respectively, Our goal will be to modify (4.15) suitably by drawing insights from string field theory so that we can evaluate the integral without having to invoke analytic continuation. For generic s, t and u we can invoke the results of §4.1 to arrive at such an expression. We cut out small regions R s , R u and R t around 1, 0 and ∞ respectively, and denote by R the left over region in the complex plane. Then we can express the amplitude as where the last three terms are boundary terms which can be evaluated from the bulk integrand using (4.7) as follows. Suppose that near σ = 0 where α i , β i are given by −u/4 − 1 plus non-negative integers and C i are constants. Then (4.7) gives Similarly if we denote the power series expansion of the integrand in the region around σ = 1 as then the required boundary terms on ∂ R s can be expressed as Finally if we expand the integrand in a power series expansion in |σ| −1 around σ = ∞: then the required boundary terms on ∂ R t are given by The expressions for B u , B s and B t given above involve infinite sum over states. However the sum converges and there is no difficulty in numerical evaluation of this expression. We have tested (4.18) by numerically evaluating this and comparing this with the exact result (4.17) in the regime in (s, t, u) space where the original integral (4.15) diverges. By taking the contours around 0 and 1 to be sufficiently small and the contour around ∞ to be sufficiently large, we can ensure that only finite number of terms with negative α i + β i , α ′ i + β ′ i and α ′′ i + β ′′ i give appreciable contribution to (4.20), (4.22) and (4.24).
As a simple example we can consider the domain s < −4, t < −4 and −4 < u < −2. In this case the divergences in the integral (4.19) will come only from the region near σ = 0. Therefore we can take the limit in which R s shrinks to σ = 1 and R t recedes to infinity, making B s and B t vanish. To determine B u we expand the bulk integrand in a power series in σ to get This gives, from (4.20), Since (4.18) has been shown to be independent of the choice of contours ∂ R u , ∂ R s and ∂ R t , we can take R u to be a rectangular region around the origin. As the size of R u shrinks, contribution to B u from the first term inside the square bracket increases and the other terms decrease for −4 < u < −2. By taking the contour to be sufficiently small, we can ignore the contribution from all terms other than the first term in the square bracket. Adding this to the bulk integral given by the first term in (4.18), with R given by the complement of R u in the whole complex plane, we can evaluate the full amplitude. With relatively little effort we can reduce the error in computation -measured by comparing the result to the exact expression given in (4.17) -to less than .1%. The error comes from the higher order terms inside the square bracket in (4.26). We can further reduce the error either by including these terms or by shrinking the size of R u , at the same time making sure that the error in the numerical evaluation of the integrals are not significant. 8 We shall now give an alternate expression that does not require infinite sum. This requires taking the limit in which the regions R s , R t and R u shrink to points 1, ∞ and 0 respectively.
In this limit only the terms with negative α i + β i , α ′ i + β ′ i and α ′′ i + β ′′ i survive in the boundary terms. We now use where R c u is the complement of the region R u . In the limit when R u shrinks to a point R c u will expand into the whole complex plane, but we do not want to take this limit yet since the integral (4.27) will diverge in this limit. Using (4.27) and analogous results for the s and t-channel boundary terms, we can express (4.18) as (4.28) The integral now runs over the full complex σ-plane, giving a finite expression for the amplitude since all the possible divergences near 0, 1 and ∞ have been subtracted. In a related case, this form of the amplitude was used in [38].

Comments
We end this section with a few comments on our result.
1. One could ask whether the amplitudes defined via (4.18) or (4.28) agree with the ones defined by analytic continuation. Again string field theory can be invoked to prove the equality of these different procedures. Each tree level Feynman diagram constructed from string field theory is manifestly an analytic function of external momenta except for poles when an on-shell particle propagates in the intermediate state. Therefore if we begin in a region of the external momenta where there are no divergences in the world-sheet description and the string field theory and world-sheet descriptions coincide, and then analytically continue the result to the region where the direct world-sheet description gives divergent results, the result of analytic continuation will coincide with the one calculated from string field theory ı.e. eqs.(4.18) and (4.28).
2. Witten gave an operational procedure for performing the world-sheet integrals that coincides with the analytic continuation of the world-sheet description [22].

Higher point amplitudes
We shall now generalize the procedure described above to higher point amplitudes of generic vertex operators for generic compactification. As in §4 we shall decompose each propagator into two parts according to (1.3) and include the first part inside the square bracket as part of bulk integration while the second term inside the square bracket will provide a boundary term at |q| = e −Λ . The main strategy will be to begin with the boundary terms that are implied by string field theory and then convert them to a general coordinate system in the moduli space. Our goal will be to arrive at an algorithm to compute the boundary terms in a general coordinate system without explicitly knowing the relation between the variables q and the general coordinates used to parametrize the moduli space.
We shall begin by introducing a convenient way of labelling the propagators of a Feynman diagram. For a tree diagram, cutting a propagator divides the external states into two disconnected sets, with each set containing at least two external states. We associate to each propagator an index s that specifies this division. Therefore s can take N possible values, where N is the total number of ways in which the set of (n + 3) external states can be divided into two sets, with each set containing at least two external states. In a given Feynman diagram different propagators carry different labels, but two propagators from two different Feynman diagrams will carry the same label if cutting them leads to the same division of external par- The − sign in front of C (1) corresponds to the fact that the normal to C (1) is taken to be along the direction of increasing |q s |, pointing into C (0) . The region inside C (0) corresponds to |q s | > ǫ.
Even though for numerical evaluation it may be better to work with small ǫ, all our formulae will be valid for finite ǫ ≤ 1.
We now denote by C   We can regard C (2) s 1 s 2 as an oriented subspace that is anti-symmetric under the exchange of s 1 and s 2 since this exchanges their transverse coordinates |q s 1 | and |q s 2 |. Generalizing this we define C (k) s 1 ···s k to be the codimension k subspace obtained by setting |q s | = ǫ for s = s 1 , · · · , s k , |q s | ≥ ǫ for s = s 1 , s 2 , · · · , s k , (5.3) in all Feynman diagrams that contain propagators with labels s 1 , · · · s k . It is easy to verify the relations and where the sum over s runs over all labels for which C where m (s 1 ,···,s k ) denote the coordinates of the 2(n − k) dimensional moduli space of the punctured spheres into which the original sphere degenerates when we set q s 1 = · · · = q s k = 0. In a Feynman diagram of string field theory, these moduli come from the moduli integration appearing in the definition of the interaction vertices, as well as the parameters q associated with the propagators other than those carrying labels s 1 , · · · , s k . These, together with the phases of q s 1 , · · · , q s k form the coordinates of C We can now express the total contribution to the amplitude as: The sum runs over all inequivalent combinations {s 1 , · · · , s k } for which C (k) s 1 ···s k exists. Using (5.6) with k replaced by k + j, one also finds that near C Note that the C i 's in (5.7) and (5.9) are not the same, but the latter are obtained by expanding the former near C (k+j) s 1 ···s k r 1 ···r j . It is easy to verify from (5.7), and (5.9) with j = 1, k replaced by (k − 1), that in the neighborhood of C (k) Note that (5.8) and (5.10) are written in the coordinate free notation. Therefore we can use these to give an expression for the amplitude in an arbitrary coordinate system in the moduli space without referring to the coordinates q s ,q s . However we need to first check if these equations fix the amplitude uniquely. First we see from (5.10) that for non-integer γ there is an ambiguity in determining I (k) for given I (k−1) of the form: where J (k) s 1 ···s k . Due to (5.10) with k replaced by k + 1, this requires us to change I (k+1) s 1 ···s k s k+1 as Therefore the net ambiguity in the whole set {I (k) s 1 ···s k } for all k and {s 1 , · · · , s k } satisfying (5.10) takes the form We shall now compute its effect on the total integral I defined in (5.8). We get (5.14) Using (5.5) we get On the other hand, replacing the summation variable k by k + 1 we can express the second term in (5.14) as where in the last step we have rearranged the sum by relabelling {s 1 , · · · , s ℓ−1 , s ℓ+1 , · · · , s k+1 } as {s 1 , · · · , s k } and s ℓ as s, and used the antisymmetry of C (k+1) to write (−1) k+1−ℓ C (k+1) s 1 ···s k s k+1 as C (k+1) s 1 ···s k s . Substituting (5.15) and (5.16) into (5.14) we see that the two terms cancel and we get ∆I = 0 . (5.17) Therefore the ambiguity (5.13) in the determination of I (k) by solving (5.10) does not affect the integral I given in (5.8). This means that we can solve (5.10) to determine the boundary corrections directly by examining the form of the bulk integrand in any given coordinate system, without having to know the relation between the chosen coordinates and the q s 's near the boundaries. As already emphasized before, this requires us to work with generic external momenta so that the exponents γ are not integers. Otherwise the determination of I (1) s could have an additive ambiguity proportional to q −1 s dq s d 2(n−1) m (s) C(m (s) ) that is not an exact differential dJ (1) s . Similar ambiguities could be present in the other I (k) s 1 ···s k 's as well. Now the expression for the amplitude given in (5.8) depends not only on the differential forms I (k) s 1 ···s k but also on the boundaries C (k) s 1 ···s k on which these forms are to be integrated. Since these are given by |q s 1 | = · · · = |q s k | = ǫ, the shapes of these boundaries in a generic coordinate system depend on the relation between the special coordinates {q s ℓ } induced from string field theory and the coordinates we are using (which could be taken to be the positions σ 1 , · · · , σ n of n of the vertex operators keeping three of them at fixed positions). We shall now demonstrate that (5.8) is actually invariant under arbitrary deformations of C  By making a k → (k+1) shift in the first term and relabelling s ℓ as r, {s 1 , · · · s ℓ−1 , s ℓ+1 , · · · , s k+1 } as {s 1 , · · · , s k }, we see that the two terms in (5.22) cancel and we have δ s I = 0 .
s . Using these results we arrive at the following algorithm for evaluating an (n + 3)-point amplitude: 1. We can take three of the vertex operators at fixed locations y n+1 , y n+2 and y n+3 and take the locations σ 1 , · · · σ n of the other vertex operators as coordinates of the moduli space.
2. We now list all possible ways of dividing the set 1, · · · , (n + 3) into two subsets, with each subset containing at least two particles. Associated with each such decomposition, labelled by s, there is a possible degeneration in which the original sphere degenerates into two spheres, each carrying external punctures belonging to one of the two subsets.
For every such degeneration, we can construct, in terms original moduli {σ i }, a new complex parameter u s that vanishes at the degeneration and a set of (n − 1) complex moduli m (s) labelling the moduli of the spheres left after degeneration. The change of variables from {σ i } to u s , m (s) may not be globally defined over the whole moduli space spanned by m (s) and we typically will have to use different u s in different patches.
3. Multiple degenerations will correspond to several of the u s 's vanishing simultaneously.
Near a degeneration where u s 1 , · · · , u s k vanish, we can construct, in terms of the {σ i }'s, a set of (n − k) complex coordinates m (s 1 ,···,s k ) which remain finite at the degeneration, labelling the moduli of the spheres into which the original sphere degenerates. We can use u s 1 , · · · , u s k and m (s 1 ,···,s k ) to parametrize the moduli space near such degenerations.
5. Once we know the relation between the coordinates σ 1 , · · · , σ n and {u s } and {m (s 1 ,···,s k ) }, we can determine I for appropriate constants α s 1 s 2 as a power series in u s 1 , u s 2 and their complex conjugates, with the coefficients of expansion given as functions of m (s 1 ,s 2 ) . This can be done by using the known relations between the coordinate systems {u s 1 ,ū s 1 , m (s 1 ) }, {u s 1 ,ū s 1 , u s 2 ,ū s 2 , m (s 1 ,s 2 ) } and {u s 2 ,ū s 2 , m (s 2 ) }, each of which can be related to the original coordinates {σ 1 , · · · , σ n }. Given I

Amplitudes in superstring theory
In superstring theory the computation of amplitudes suffers from divergences similar to the ones described above for bosonic string theory -appearing from propagation of internal states carrying negative L 0 +L 0 eigenvalue. However there is another source of divergences in superstring theory -from the wrong choice of PCOs. Consider for example a pair of integrated −1 picture vertex operators e −φ ψ µ∂ X ν (z,z) and e −φ ψ ρ∂ X σ (w,w). The rule that follows from superstring field theory is that when we bring them close, we must also bring a PCO close to them so that their product generates a −1 picture state. This can be achieved e.g. by converting the second vertex operator to a zero picture vertex operator ∂X ρ∂ X σ (w,w). The leading singularity in the collision of the two vertex operators now is proportional to (z −w) −2 e −φ ψ µ ∂X ρ η νσ . If instead we had taken the product of the two −1 picture vertex operators, then we would have gotten a stronger leading singularity proportional to (z − w) −2 (z −w) −2 e −2φ η µρ η νσ . A similar situation would arise if we had converted both vertex operators to zero picture before bringing them close to each other. Now often it is convenient to take the PCO's to coincide with some vertex operators before computing the correlation function and use the same arrangement everywhere in the moduli space. According to the discussion in the preceding paragraph this would give the wrong integrand I (0) near many degenerations, e.g. when two or more zero picture vertex operators come together their product would carry a net picture number 0 instead of −1. Therefore this would be the wrong starting point for implementing the procedure described in §5. However one can show that the error that one makes by taking the wrong choice of PCO locations adds to I (0) an exact differential in C (0) . We shall now argue that the addition of such exact differentials does not affect the amplitude I defined in (5.8). For this let us suppose that we change I (0) to I (0) + dL (1) where L (1) is a globally defined 2n − 1 form on C (0) . Then the equation dI We can now study the effect of the change in I (0) and I (1) on I defined in (5.8). We have Therefore the net change in the sum of the left hand sides of (6.2) and (6.

Mass shift under marginal deformation
We shall now consider a different situation in tree level string theory where short distance divergences appear in the world-sheet of the string, and show how string field theory removes the divergence. The case that we shall describe involves marginal deformation of the matter CFT that is used in describing the target space geometry. Since the deformed background is also a conformal field theory, one can formulate string theory around the deformed background as well. The spectrum and S-matrix of the new theory is in principle computable from the data in the original CFT before the deformation using conformal perturbation theory. This however requires ultraviolet regularization on the world-sheet, since the naive computation requires integrating correlation functions of marginal operators and other operators over the locations of the marginal operators and they diverge when the marginal operators collide with each other or with other operators. Our goal will be to show how in string field theory we can carry out the computation without encountering any divergence.
Although the procedure we shall describe is valid for any marginal deformation, in order to get concrete results we shall focus on a particular class of examples where the target space includes a compact circle and the marginal deformation corresponds to deforming the radius of the compact circle. We shall denote the world-sheet scalar field corresponding to the compact circle by Y . In string field theory the marginal deformation that changes the radius of the circle can be represented as the effect of switching on a background string field solving the classical equations of motion. We shall use string field theory to compute the shift in the masses of various states under this deformation to second order in the deformation parameter and compare the result to known results. In doing this we shall make use of the general strategy described in [34] for solving classical equations of motion of string field theory. During this analysis we shall see that while in the intermediate stages of calculation the results depend on the choice of local coordinates used to define the theory, the final result is independent of this choice.

Bosonic string theory
The equations of motion of bosonic string field theory is given by (2.20): We shall solve this equation in a power series in the parameter µ: by starting with a solution It is clear that this solves (7.1) to order µ since |Ψ cl 1 is annihilated by Q B . The second order correction µ 2 |Ψ cl 2 to the solution can be expressed as [34] |Ψ cl 2 = − where P is the projection operator into the L + 0 = 0 states and |ψ cl 2 ∈ P H satisfies: It is easy to see however that for the choice given in ( For the first two states in (7.6), {φΨ cl 1 Ψ cl 1 } vanishes since in the three point function the total ghost number has to add to 3 separately in the left and the right sector. For the last state in (7.6), {φΨ cl 1 Ψ cl 1 } vanishes since the three point function of three ∂Y∂Y 's vanishes due to separate Y → −Y symmetry in the left and the right sector. 9 Therefore we can take Since the states listed in (7.6) are all BRST invariant, we could include an arbitrary linear combination of these states in the definition of |ψ cl 2 . This will correspond to a redefinition of µ or the string coupling constant, or a gauge transformation.
Even though we shall use the solution to order µ 2 , it is instructive to examine how the solution can be extended to order µ 3 . We can solve (7.1) by taking [34] |Ψ cl Possible obstruction to extending the solution to order µ 3 can arise from the failure to solve (7.9). Using (7.4) and (7.7), we can express this as A consistency condition for the existence of a solution to this equation is the vanishing of the inner product of both sides with φ|c − 0 for the three states |φ listed in (7.6). The inner product of the left hand side with all the states vanish since all of these states are BRST invariant. The inner product of the right hand side with the first two states in (7.6) vanish due to the vanishing of the correlation functions of odd number of ∂Y 's and / or odd number of∂Y 's. Since the last state in (7.6) is given by |Ψ cl 1 , we need to check the vanishing of Formally this can be regarded as the four point function of four Ψ cl 1 states, with the first term representing the contribution from the elementary vertex and the second term representing the contribution from the s, t and u channel diagrams all of which are equal. However unlike the four point function which has divergences from collision of vertex operators, (7.11) is manifestly finite. To evaluate it we shall choose a local coordinate system of the type described in §2 with large λ. We emphasize however that the choice of large λ is a matter of convenience, but is in no way necessary.
Since the external momenta all vanish and are not generic, we can have ambiguities in applying the procedure of §4,5. Therefore we shall evaluate (7.11) directly. First we have where we have used ∂R = −∂R s − ∂R u − ∂R t . Since all the external states are identical, the contribution from ∂R u , ∂R s and ∂R t are identical. So we shall calculate the contribution from one of the boundaries ∂R u around σ = 0 and multiply the result by 3. Now we see from (3.13) and (3.18) that in the large λ limit, ∂R u given by |q| = 1 corresponds approximately to a circle around σ = 0 of radius of order λ −2 . If we drop terms proportional to negative powers of λ, which will drop out when we take the λ → ∞ limit, then the only contribution to (7.12) comes from the σ −2σ−1 term in the integrand. This leads to Since the ∂R u is not strictly a circle in the σ-plane, it will be more convenient to express the boundary term in the q plane via the relation (3.18). This gives , (7.14) where ξ = q/λ 2 . Since the integral has to be evaluated at |ξ| = 1/λ 2 , we can expand g(ξ) in a power series in ξ for large λ. Using (3.13) we have Substituting this into (7.14) we get where ≃ means that we have dropped terms involving inverse powers of λ.
The second term on the right hand side of (7.11) is given by where |ξ p denote a complete basis of ghost number 2 states in H and |ξ c p is a conjugate basis of ghost number 4 states in c − 0 H satisfying In arriving at the right hand side of (7.17), we have replaced [Ψ cl using the fact that [Ψ cl 1 Ψ cl 1 ] ∈ H, and then inserted a sum over complete set of states between b − 0 and c − 0 . Now from (3.2) we see that if ξ p has L + 0 eigenvalue ∆ p , then for dimension 0 primaries A, B, the {ABξ p } carries a factor of |dw 3 /dz| −∆p Since Ψ cl 1 is a dimension 0 primary, it follows that for large λ we can restrict the sum over p, q in (7.17) to states with ∆ p ≤ 0. Now the contributions from the ∆ p = 0 states are already removed by the projection operator (1 − P ), therefore we have to focus on states with ∆ p < 0. The only such state is the ground state with ∆ p = −2, giving Therefore we have This gives (7.21) (7.21) exactly cancels (7.16). Therefore (7.11) gives This in turn shows that there is no obstruction to extending the solution to order µ 3 . This is equivalent to proving the vanishing of the β-function for the marginal deformation V to order µ 3 without having to introduce ultraviolet cut-off at any step. Earlier explicit analysis of marginal deformations in closed string field theory [37] examined the solution to second order and therefore did not encounter dependence on the local coordinate λ h 3 (z) in the intermediate steps of the calculation. As stated below (7.11), and noted in [37], at order λ 3 the existence of the solution formally requires the vanishing of the four point function of marginal operators. However the corresponding integrand diverges when vertex operators collide and one has to regularize it. In contrast, string field theory gives finite result by treating (L + 0 ) −1 correctly. We now turn to the computation of the spectrum of the deformed theory to order µ 2 . For this we need to study the fluctuation of the string field around the new background. Let us define: Then the equation of motion (7.1) to linear order in Φ is given by: Let us for definiteness consider the tachyon state carrying momentum n/R along the compact direction and k along the non-compact directions. Starting with the leading order solution to (7.24), we shall compute systematic corrections to the solution following the procedure described in [34]. We shall denote by |Φ ℓ the full solution to order µ ℓ . To order µ 0 the solution to (7.24) takes the form: For subsequent analysis we need to introduce a projection operator P that projects to states which carry momentum k = k (0) + O(µ) along non-compact directions, n/R along the compact direction and have L + 0 = O(µ). It is easy to verify that in this case only vertex operators invariant under P are those of the form cce ik.X e inY /R . We normalize |Φ ℓ such that |φ ℓ ≡ P|Φ ℓ = c 1c1 e ik (ℓ) .X (0)e inY /R (0)|0 , k (ℓ) = k (0) + O(µ) . (7.26) As discussed in [34], since k (ℓ) is expected to be different from k (0) due to a change in the mass, it is not convenient to use the solution |Φ 0 and correct it to obtain |Φ ℓ . Instead we begin with |φ ℓ given in (7.26) as the seed solution and correct it to order µ ℓ to determine |Φ ℓ and k (ℓ) . We begin with the ansatz for the solution |Φ 1 : (1) .X e inY /R , k 2 (1) = k 2 (0) + a 1 µ , (7.27) where the constant a 1 will be determined shortly. Taking (7.27) as the leading order solution, we can substitute this into the right hand side of (7.24) to get a solution to order µ: provided |φ 1 satisfies: If we defineφ (1) .X e −inY /R , (7.30) satisfying the normalization φ 1 |c 0c0 |φ 1 = −1 , (7.31) and take the inner product of φ 1 |c − 0 with (7.29), we get where in the last step we have used the fact that Ψ cl 1 is a dimension 0 primary and φ 1 andφ 1 are primaries of dimension of order µ. Only the (c 0 L 0 +c 0L0 ) term in Q B contributes to the left hand side. This gives We now turn to order µ 2 computation. We begin with the ansatz with |φ 2 as defined in (7.26), and solve the equations of motion (7.24) iteratively to get successive order solutions: and provided |φ 2 satisfies, to order µ 2 , In the second line we have used the expression for |Ψ cl 2 given in (7.4), (7.7). Defining (2) .X e −inY /R , (7.38) and taking the inner product of φ 2 |c − 0 with (7.37), using the analog of (7.31) with φ 1 ,φ 1 replaced by φ 2 ,φ 2 , we get We shall now evaluate the different terms appearing on the right hand side of (7.39). First of all we have: where the (λ|h ′ 3 (1)|) −k 2 (2) +4−n 2 R −2 factor arises from the fact that the states φ 2 ,φ 2 have L + 0 eigenvalues ∆ = (k 2 (2) − 4 + n 2 R −2 )/2 and therefore {φ 2 Ψ cl 1 φ 2 } will carry factors of |dw 3 /dz| −∆ z=z i for i = 1, 3. Using (3.1), (3.2) and (3.4) we get (dw 1 /dz) z=z 1 = (dw 3 /dz)| z=z 3 = λh ′ 3 (1). To evaluate the right hand side of (7.40) to order µ 2 , we can replace k (2) by k (1) at the cost of making an error of order µ 2 . Using (7.33) we get The evaluation of the last three terms on the right hand side of (7.39) can be simplified by noting that since all these terms already have explicit factors of µ 2 , we can replace the momenta k (2) by k (0) given in (7.25). This makes all the external states on-shell.
Let us begin with the evaluation of the last term in (7.39). This has the same expression as the on-shell four point function of the states φ,φ, Ψ cl 1 and Ψ cl 1 , except that the integration over the moduli runs over the region R introduced in §3.1 instead of the whole complex plane. The integrand for the on-shell amplitude is given by where and While the analysis can be carried out for any choice of local coordinate system encoded in the functions λ h i (z) introduced in §3.1, we shall simplify our analysis by taking λ to be large.
In this case it follows from (3.10), (3.18) and (3.23) that the excluded regions |q| ≤ 1, denoted by R s , R u and R t in §3.1, correspond to small regions around 1, 0 and a region outside a large radius respectively. Therefore for large λ we can ignore the excluded regions R s , R t and/or R u as long as the integrands do not encounter any divergence. 10 For example in I 1 we can ignore the excluded regions R s and R t since the integrand does not have any divergence from 1 and ∞, but cannot ignore the excluded region R u around the origin. On the other hand in I 3 we can ignore the excluded region R u around the origin, but cannot ignore R s and R t .
Let us begin with the evaluation of I 1 . Since the only excluded region in this expression is a small region R u around 0, we can express this as where we have used ∂R = −∂R u . This has the same structure as (7.13) and can be evaluated using identical procedure, leading to the analog of (7.16): where ≃ denotes equality up to terms containing inverse powers of λ.
I 2 can be evaluated by expressing this as We can evaluate this using integration by parts, picking up boundary contributions from ∂R.
The only boundary that contributes is the boundary −∂R s around σ = 1. This gives I 3 can be evaluated by noting that the integral has logarithmic divergence both near 1 and ∞, and therefore we need to use eqs. (3.10) and (3.23) to determine the cut-off on σ integral by identifying the |q| = 1 curves. Using (3.13) we see that in the σ = y 4 plane these curves are and respectively. Since the integrals are at most logarithmically divergent we do not need to know the corrections to these curves. This gives Using the values of I 1 , I 2 , I 3 determined above, we finally get from (7.43) Let us now turn to the contribution from the third term on the right hand side of (7.39).
Evaluation of this proceeds exactly as that of (7.17) with the only difference that the first two factors of Ψ cl 1 are replaced byφ 2 φ 2 . The first equation in (7.20) is replaced by This gives the analog of (7.21): The contribution from the second term on the right hand side of (7.39) vanishes in the large λ limit due to the following reason. As already argued before, for evaluation of this term we can set k (2) = k (0) . In this case momentum conservation forces the intermediate state to carry a factor of e ik (0) .X+inY /R which has L + 0 = 2. Therefore the lowest L + 0 eigenvalue state propagating in this channel is cce ik (0) .X+inY /R with L + 0 = 0. The contribution from this state is projected out by the (1 − P) operator insertion. Therefore the contribution comes only from states with positive L + 0 eigenvalue, leading to terms with negative powers of λ. Such terms will vanish in the large λ limit. Therefore We now have from (7.39), (7.41), (7.54), (7.56) and (7.57), This gives Therefore the mass 2 of the state to order µ 2 is given by This is consistent with the expectation that the marginal deformation induces a deformation of the radius R of the compact direction. Note that we have arrived at this result without encountering any divergence from collision of the pair of marginal operators or of the marginal operator with the vertex operator of the tachyon. The µ dependent terms can be regarded as an expansion of e −µ but we do not have a compelling reason to believe that this pattern will continue to hold at higher order.

Heterotic string theory
We shall now repeat the analysis of the previous section for heterotic string theory. The marginal deformation corresponds to switching on NS sector string field. For definiteness we shall take the states, whose mass shift we compute, also to be in the NS sector, but the generalization to Ramond sector states is straightforward. Since in the NS sector the structure of heterotic string field theory is identical to that of bosonic string field theory, most of the analysis takes identical form. For this reason, we shall mention only the differences.
The first difference is in the form of the leading order classical solution describing the shifted vacuum. We take this to be of the form: where we have used the notation of §2 for various fields. The form of (7.4) and (7.5) remain unchanged. The possible list of L + 0 = 0 states in H which have the correct total ghost and picture numbers and (Y → −Y, χ → −χ) symmetry for having non-zero inner product with However all of these actually have vanishing inner product with c − 0 [Ψ cl 1 Ψ cl 1 ] due to separate ghost charge conservation in holomorphic and anti-holomorphic sectors and / or the χ → −χ symmetry. Therefore we can take |ψ cl 2 = 0 as in bosonic string theory. Showing that the solution extends to order µ 3 is also straightforward, but we shall not describe it here since we only use the solution to order µ 2 .
Analysis of the fluctuations around the deformed vacuum also proceeds as in the case of bosonic string theory. We replace the tachyon carrying spatial momentum k and internal momentum n/R by a massless field with the same momenta. The ansatz for φ ℓ andφ ℓ , replacing (7.26), (7.30) and (7.38) takes the form: To first order in µ we get the analog of (7.32): where p is the PCO location on the NS-NS-NS interaction vertex chosen according to the rules discussed in §3.2. Explicit evaluation gives a result independent of p: At order µ 2 we have the analog of (7.39): The first term on the right hand side can be evaluated exactly as in bosonic string theory, leading to the analog of (7.41): In this case the contribution from the second and the third terms on the right hand side of (7.67) carry only negative powers of λ since there are no states with L + 0 < 0 and the contributions from the L + 0 = 0 states are removed by the projection operators (1 −P ) and (1 −P). Therefore we are left to evaluate the last term.
As in the bosonic string theory, while evaluating φ 2 φ 2 Ψ cl 1 Ψ cl 1 we can take the external states to be on-shell. The integrand depends on the choice of PCO locations but under a change in PCO locations the integrand changes by a total derivative. Therefore we can adjust them at will in the interior of R, as long as on ∂R they coincide with the PCO locations in ∂R s , ∂R t and ∂R u , fixed by the arrangements described in §3.2. Our strategy will be to take the PCO locations to coincide with the locations of the vertex operators Ψ cl 1 in the interior of R, and at the boundary ∂R make them jump to the values they take inside R s , R u and R t . The effect of this jump can be computed via vertical integration [44,45]. If we had chosen a different PCO assignment inside R, the bulk integrand will change by a total derivative. However the result of vertical integration will also change, precisely cancelling this effect.
The effect of taking the PCO locations to the locations of Ψ cl 1 is to convert the Ψ cl 1 inserted at 0 to an unintegrated 0 picture vertex operator given in (2.12) Using these we get where and B v is the result of vertical integration that moves the PCOs from their locations at 0 and σ to the required values on ∂R so that they coincide with the positions of the PCOs on ∂R s , ∂R t and ∂R u given in §3.2. If we denote by W 1 and W 2 the final locations of the PCOs, and follow the convention that we first move the PCO at 0 to W 1 and then move the PCO at σ to W 2 , then this has the effect of replacing inside the correlation function [44,45] the factor Only the terms involving σb (w)dw survive after imposing ghost charge conservation in the anti-holomorphic sector. This gives First let us compute the contribution from the first term on the right hand side of (7.71). The φ-charge conservation (or equivalently ξ-η charge conservation) tells us that the term proportional to η in the second line of (7.72) does not contribute. Evaluation of the rest of the correlator gives This is identical to the integrand in (7.42). Therefore we get, using (7.54) Next we turn to the analysis of B v given in (7.75). Using the relation we can express B v as We shall furthermore express B u as and We have shown in appendix B that in the large λ limit, B s , B t and B ′′ u vanish so that Intuitively the vanishing of B s , B t and B ′′ u may be understood as follows. Let us start with B s . For large λ the integration region is near a degeneration where the points 0 and ∞ are on one sphere and the points σ and 1 are on another sphere, with the two spheres connected by a narrow neck. In this case the dominant contribution comes from nearly on-shell states propagating along the neck. In the initial configuration the PCO at 0 lies on the first sphere and the PCO at σ lies on the second sphere. In the final configuration one of the PCOs (W 1 ) lies on the first sphere and the other PCO (W 2 ) lies on the second sphere. B s describes the effect of moving the first PCO from 0 to W 1 and moving the second PCO from σ to W 2 . Each of the PCOs remains on its own sphere. Since for on-shell three point function, moving the PCO on the sphere does not have any effect, we expect B s to vanish in the large λ limit. Similar argument can be given for B t . This does not apply to B u since the initial position of the two PCOs, at 0 and σ, lie on the same sphere while the final arrangements W 1 and W 2 must lie on different spheres. We analyze this as a combination of two moves: first move one of the PCO's from 0 to 1, and call this contribution B ′ u , and then move the pair of PCO's at 1 and σ, which are now on different spheres, to W 1 and W 2 . The latter contribution, called B ′′ u vanishes in the large λ limit due to the same arguments as for B s and B t .
We shall now analyze B ′ u . Due to ξ-η charge conservation, only the term proportional to η inside the curly bracket contributes. This gives: For large λ, ∂R u is a small contour around the origin. In this case the contribution from the term proportional to (σ − 1) −1 is suppressed by inverse power of λ and can be neglected. The remaining term has the same structure as the right hand side of (7.47) with the roles of σ and σ interchanged. This generates an extra minus sign, leading to (7.48) multiplied by −1:  (7.87), and the fact that the contribution from the second and third terms on the right hand side of (7.67) vanishes, we get This gives Although the µ dependent corrections have the same form as (7.59), there is a subtle difference in the analysis. In the case of bosonic string theory the potential short distance divergence, reflected in the term I 1 in (7.48) that grows as λ 4 in the large λ limit, can be identified to the contribution from the intermediate tachyon state. The naive divergence in the world-sheet integral over σ can be traced to the wrong treatment of the tachyon propagator as in (1.2). In contrast the potential short distance divergence in the heterotic string theory, reflected in the term proportional to λ 4 in (7.77), can be traced to the wrong assignment of PCO locations, and is cancelled by the boundary term B ′ u that corrects the PCO location via vertical integration.

Comments of higher genus amplitudes
As we have emphasized earlier, there is a subtle difference between the analysis in §4,5,6 and that in §7. The analysis in sections 4, 5 and 6 was carried out in a manner that is manifestly independent of the string field theory data -choice of local coordinate system and PCO locations. On the other hand the analysis in §7 required, in the intermediate steps, use of the local coordinate system, e.g. the function h 3 (z) in (7.41) and g 4 (z) in (7.15), although at the end the dependence cancelled. There is an intrinsic reason for this difference. The on-shell amplitudes discussed in §4,5,6 are expected to be genuinely independent of the choice of string field theory data, and so it is not surprising that the analysis can be made independent of these choices. The result of §7, describing the effect of marginal deformation, is not expected to be manifestly independent of the string field theory data although the dependence on these data is expected to be removable by a redefinition of the deformation parameter µ. This is due to the fact that a change in the string field theory data causes a redefinition of the string fields including the one corresponding to the marginal operator, and such field redefinitions will induce a redefinition of the deformation parameter. Therefore the total independence of the final results (7.59), (7.89) on the string field theory data is accidental, and we expect that in general there will be such dependence of the result on these data.
It is natural to ask if the analysis can be generalized to give a systematic procedure for computing higher genus amplitudes that gives manifestly finite results, and yet minimizes the dependence on the choice of string field theory data. To this end note that for higher genus amplitudes, the dependence on the string field theory data of the kind presented in §7 will always be present even in the absence of marginal deformations. This is due to the fact that at higher genus, under quantum corrections there will be mass and wave-function renormalizations of all the external states. Now since different choices of string field theory data lead to string field theories that are related by field redefinition, there will be two effects [34,46,47]. First the amplitudes computed in different string field theories will differ due to different wave-function renormalizations of external states. Second the definition of the moduli fields, e.g. the string coupling constant encoded in the dilaton, will differ in different string field theories, causing a change in the amplitude. These will lead to ambiguities in the final result that depend on the string field theory data. However such ambiguities can be absorbed into a finite renormalization of external states and the values of the moduli. Therefore, as in this paper, one could proceed with the computation assuming the existence of a consistent set of string field theory data without making a particular choice and then verify at the end that the final result depends on these data only through the normalization of the external states and definition of the moduli fields.
We shall end this section by describing the different kinds of degenerations that we need to deal with for higher genus amplitudes: 1. The first type of degeneration is separating type degeneration with generic momentum flowing across the degenerating punctures. These correspond to the original Riemann surface degenerating into a pair of Riemann surfaces, each of which carries two or more external punctures. These can be treated in the same way as in §5,6 and do not introduce any ambiguity in the final result.
2. The second type of degenerations, analyzed extensively in [47], involve separating type degenerations with special momentum flowing across the degenerating punctures. Examples of these involve degenerations where the original Riemann surface degenerates into a pair of Riemann surfaces, one of which carries either no puncture or one puncture. The first one represents tadpole type diagrams with zero momentum flowing across the degenerating puncture while the second one represents mass and wave-function renormalization diagrams with on-shell momentum flowing across the degenerating puncture.
Near these degenerations the expansion of the integrand in powers of the variables u,ū that vanish at the degeneration involves integer exponents and we cannot apply the general trick of §5,6 to remove these divergences. There may be genuine divergences of the form d 2 u/|u| 2 signaling the presence of massless tadpoles and renormalization of physical masses. There may also be ambiguities in determining the I (k) 's of the form du/u encoding redefinitions of massless moduli fields and/or external states. We need to use the analog of the procedure described in §7 to address these cases. The general procedure based on string field theory can be found in [34]. In simple cases, one can follow the procedure described in §7 to minimize the dependence on the explicit knowledge of string field theory, reproducing the results in [48][49][50].
3. The final category of degenerations involves non-separating type degenerations -degenerations where one pinches the handle of a Riemann surface but the Riemann surface still remains connected after degeneration. This case lies in between the two cases described earlier, in that the momentum flowing across the degenerating puncture is a loop momentum that needs to be integrated. For most of the range of integration over momenta the momentum is generic, but on codimension one subspaces of the loop momentum space the L 0 +L 0 eigenvalue corresponding to that momentum may vanish, causing the integrand to diverge. In four or less dimensions the momentum integrals themselves are divergent reflecting the presence of infrared divergences in the theory and the procedure for getting finite result is complicated -requiring the same methods that are normally used in quantum field theories with massless fields. However in higher than four dimensions there is no genuine divergence and one should be able to extract finite results. Nevertheless the answer is not free from ambiguity without additional input, since the iǫ prescription in the integration over the loop momenta is hidden in the prescription of how we deal with the pole in the 1/(L 0 +L 0 ) factor. This can be done following the procedure described in [22], where we replace the (L 0 +L 0 ) −1 factor by the representation (1.4). This means that if the original integral has the form where we have suppressed the integration over the other moduli, we replace it by 11 11 Even though the Cutkosky rules are not manifest in the procedure described in [22], this has been shown −2i The relation between the variables q,q and some predetermined moduli parameters on the punctured Riemann surface can be determined if we know the string field theory data. However this is not necessary. Let us suppose that u,ū correspond to some other set of variables with the property that u vanishes linearly with q near the degeneration.
Parametrizing u as e −t−iφ , we can conclude that if the original integrand has the form The equivalence between these two procedures -one in q,q variable and the other in the u,ū variable -follows from the fact that near the degeneration q and u are related to each other in a one to one fashion. Therefore if we map the integration contour used in (8.4) to the (s, θ) plane, the resulting contour can be smoothly deformed to the contour in (8.2) without passing through any singularity. In dimensions larger than four, and for generic momenta of external states, the integrand F (u,ū) has power law suppression in inverse powers of − ln |u| that makes the integral convergent from the large Im (t) region even without the e iǫt factor. Therefore we can compute the integral directly, without having to take limits of integrals. This remains true even for multiple integrals of this kind where more than one handle degenerates simultaneously. An explicit example of this for one loop amplitude can be found in [53] (see e.g. eq.(3.16)). A Choice of local coordinates on a five punctured sphere In this appendix we shall describe, for the five punctured sphere, the relation between the global coordinates (σ 1 , σ 2 ) and the coordinates u s 1 , · · · , u s k , m (s 1 ,···,s k ) for some choices of {s 1 , · · · s k }.
We begin by making a specific choice of global coordinates. This will be done by fixing the puncture locations y 1 , · · · y 5 in the complex plane to be at y 1 = σ 1 , y 2 = σ 2 , y 3 = 1, y 4 = 2, y 5 = 0 . (A.1) The boundaries we shall consider are shown in Fig. 2. These are similar to the Feynman diagrams shown in Fig. 1, but Fig. 2 should be regarded as depicting regions near those represented by the Feynman diagrams in Fig. 1 for |q| ∼ e −Λ -not necessarily the precise regions that follow from a string field theory. For example Fig. 2(i) represents region near the boundary C (1) a where the original sphere is near a degeneration into a four punctured sphere carrying the original punctures 1, 2 and 5 and a three punctured sphere carrying the original punctures 3 and 4. Our goal will be to introduce the coordinates u a , m (a) near this boundary in terms of the global coordinates σ 1 , σ 2 . This is done as follows. Let us take the three punctured sphere on the left, carrying global coordinate z, to have puncture 3 at z = 1, puncture 4 at z = 2 and the sewing puncture at z = 0. We also take the four punctured sphere on the right, carrying global coordinate z ′ , to have puncture 1 at z ′ = m (a) , puncture 2 at z ′ = 2, puncture 5 at z ′ = 1 and the sewing puncture at z ′ = 0. m (a) should keep a finite distance away from 0, 1 and 2 so that the four punctured sphere is not close to degeneration -as will be discussed later, for m (a) close to 0, 1 or 2, we need to choose the coordinate systems differently. We now to be equivalent to the procedure of [51] and therefore satisfies the Cutkosky rules [52]. We have usedû a instead of u a to take into account the fact that the correct candidate for the coordinate u a may have different form in different domains in the moduli space. We shall see that while theû a appearing in the above equation is the correct choice of u a as long as m (a) is kept away from 0, 1 and 2, we need modifications when m (a) approaches any of these points.
In order to bring the 5th puncture at 0 leaving the third and fourth punctures at 1 and 2 respectively, we make a change of coordinates: The location of the punctures 1 and 2 in the y plane are now given by, respectively, σ 1 = 2(1 − m (a) )û â u a (û a − 3m (a) ) + 2 m (a) , σ 2 = − 2û â u a (û a − 6) + 4 . (A.5) Since we shall use this formula only for smallû a , we shall replace (A.5) by a simpler equation: a . These can be taken to represent the coordinates (u a , m (a) ) in the notation of §5 when m (a) is not close to 0, 1 or 2.
Let us now turn to the region near the boundary C (1) b represented by Fig. 2(ii). We take the four punctured sphere on the left, carrying global coordinate z, to have puncture 3 at z = 1, puncture 4 at z = 2, puncture 1 at z = m (b) and the sewing puncture at z = 0. m (b) needs to keep finite distance away from 0, 1 and 2 so that this sphere is not close to degeneration. We also take the three punctured sphere on the right, carrying global coordinate z ′ , to have puncture 2 at z ′ = 2, puncture 5 at z ′ = 1 and the sewing puncture at z ′ = 0. We now sew the two Riemann surfaces via the relation (A.7) In the z coordinate the punctures are located at: 3 : z = 1, 4 : z = 2, 5 : z =û b , 1 : z = m (b) , 2 : z =û b /2 . (A. 8) In order to bring the 5th puncture to 0 leaving the third and fourth punctures at 1 and 2 respectively, we make a change of coordinates: The location of the punctures 1 and 2 in the y plane are now given by, respectively, Again, since we shall use this coordinate system for smallû b , we replace this by a simpler set of equations: (A.11) (A.11) defines the coordinatesû b , m (b) appropriate near C (b) in terms of the global coordinates σ 1 , σ 2 of the moduli space. These can be taken to represent the coordinates (u b , m (b) ) in the notation of §5 when m (b) is not close to 0, 1 or 2.
Next we turn to the region near C (2) ab represented by the diagram 2(iii). Keeping in mind that the variablesû a andû b introduced earlier may not exactly match with the variables suitable for parametrizing the region near C (2) ab , we shall denote the new parameters byũ a andũ b . We take the left sphere carrying global coordinate z to have puncture 3 at z = 1, puncture 4 at z = 2 and the sewing puncture at z = 0, the middle sphere carrying global coordinate z ′ to have puncture 1 at z ′ = 1, the left sewing puncture at z ′ = 0 and the right sewing puncture at z ′ = ∞ and the right sphere carrying global coordinate z ′′ to have the sewing puncture at z ′′ = 0, puncture 2 at z ′′ = 2 and puncture 5 at z ′′ = 1. We now sew the three spheres via the relation zz ′ =ũ a , z ′′ /z ′ =ũ b .
(A. 12) In the z coordinate the punctures are located at: 3 : z = 1, 4 : z = 2, 5 : z =ũ aũb , 1 : z =ũ a , 2 : z =ũ aũb /2 . We now recall that for Fig. 2(i), identification of the coordinate system (û a , m (a) ) with the coordinates (u a , m (a) ) introduced in §5 breaks down for m (a) close to 0 since the right sphere degenerates in this limit. By examining the choice of coordinates of the punctures on the original sphere one can see that this degeneration is precisely the one depicted in Fig.(2)(iii). Therefore in this region we can identify (ũ a , m (a) ), instead of (û a , m (a) ), with the coordinates (u a , m (a) ) introduced in §5. Similar modifications must also be made when m (a) approaches 1 and 2 by analyzing good coordinate systems near other degenerations. The choice of the coordinate system (u b , m (b) ) needs to be similarly modified when m (b) approaches 0, 1 and 2.
For example when m (b) approaches 0 we can use (ũ b , m (b) ) to label coordinates near C (b) . Using this coordinate system we can also define the boundaries C to be given by |û b | = ǫ when m (b) is away from 0, 1 and 2 but |ũ b | = ǫ when m (b) is close to zero. These can be formally stated as follows. Let H(x) be a smooth function of a complex variable x that approaches 1 for large |x| and zero for small |x|, e.g.
H(x) ≡ |x| 2 |x| 2 + η 2 , (A. 19) where η is a fixed number. Then we define: C (1) a : |û a | = ǫH(m (a) ) + ǫ|m (a) |(1 − H(m (a) )), C When m (a) (m (b) ) is close to 1 or 2, the definitions of C (1) a (C (2) b ) need to be further modified along the same line. C (1) ab is simply the intersection of these two subspaces, given approximately by |ũ a | ≃ |ũ b | ≃ ǫ for small ǫ. Note however that we do not need to take the ǫ → 0 limit since (5.8) gives the correct result even when the C (A.21) We now obtain I where in the second line we have displayed its behavior in the coordinate system appropriate near C (2) ab using the coordinate transformations (A.18). Similarly we have We can now expand C i and C i in the second lines of (A.22) and (A.23) in power series expansion inũ b andũ a respectively to find the expressions for I In this appendix we shall show that B s , B t and B ′′ u , defined in (7.80), (7.83), vanish in the large λ limit. We begin with B s . Keeping only the part of the correlator in (7.80) that does not vanish by ghost charge conservation, we have B s = −2 µ 2 2πi ∂Rs dσ cce −φ ψ 1∂ X 2 e −ik (0) .X e −inY /R (1)cce −φ ψ 1∂ X 2 e ik (0) .X e inY /R (∞) For large λ, ∂R s represents a contour around 1 of linear size of order λ −2 . Also W 1 is a point at finite distance away from 1 and W 2 is a point within distance λ −2 of 1. We shall now estimate the integrand for σ ∈ ∂R s and show that the integral is suppressed in the large λ limit. The integration measure dσ contributes a factor of λ −2 since the contour has linear size λ −2 . Therefore the integrand must grow as λ 2 for getting a finite contribution. For the term in the second line of (B.1) the b, c,b,c correlators give finite contribution, the φ correlator gives a contribution of order (σ − 1) ∼ λ −2 , the ξ, η correlator gives finite contribution and the matter correlator gives a contribution of order (σ − 1) −1 ∼ λ 2 from the Y contribution. Therefore the correlator grows as order λ 0 and its contribution to the integral vanishes for large λ. For the term in the third line of (B.1) the b, c,b,c correlators give a contribution of order (σ − 1) ∼ λ −2 from the c-c operator product, the φ correlator gives a contribution of order (σ − 1) −1 ∼ λ 2 , the matter correlator gives a contribution of order (σ − 1) −1 ∼ λ 2 from the Y contribution and the ξ, η correlator gives a contribution of order (σ − W 2 ) ∼ λ −2 due to the fact that the points σ and W 2 where ξ is inserted are within a distance λ −2 of each other. This again makes the integrand is of order λ 0 in the large λ limit, making the integral vanish in this limit. This shows that B s vanishes.
A similar analysis can be carried out for B t which has the same expression as (B.1) except that the integration contour ∂R t lies in the large σ region (σ ∼ λ 2 ), W 1 lies at a finite point and W 2 also becomes large of order λ 2 . However this case is related to that for B s via a z → z/(z − 1) transformation accompanied by a reversal in sign of k (0) and n. Therefore vanishing of B s in the large λ limit also implies vanishing of B t .
Let us now turn to B ′′ u which, according to (7.83), is given by (0) .X e −inY /R (1)cce −φ ψ 1∂ X 2 e ik (0) .X e inY /R (∞) In this case the integration contour ∂R u encloses 0 and, according to (3.13), (3.18), represents approximately a circle of radius ∝ λ −2 around the origin, with corrections of order λ −4 . The point W 1 is at finite distance away from 0 and the point W 2 is within a distance of order λ −2 of 0. Now it follows from the analysis in §3.1, 3.2 that W 1 and W 2 are holomorphic functions of q and hence of σ. Evaluating the correlator in (B.2) we find that B ′′ u has the form: for some holomorphic function f (σ) with a regular Taylor series expansion around σ = 0. It follows from this that for large λ, B ′′ u is suppressed by inverse powers of λ.