The positive geometry for $\phi^{p}$ interactions

Starting with the seminal work of Arkani-Hamed et al arXiv:1711.09102, in arXiv:1811.05904, the"Amplituhedron program"was extended to analyzing (planar) amplitudes in massless $\phi^{4}$ theory. In this paper we show that the program can be further extended to include $\phi^{p}$ ($p>4$) interactions. We show that tree-level planar amplitudes in these theories can be obtained from geometry of polytopes called accordiohedron which naturally sits inside kinematic space. As in the case of quartic interactions the accordiohedron of a given dimension is not unique, and we show that a weighted sum of residues of the canonical form on these polytopes can be used to compute scattering amplitudes. We finally provide a prescription to compute the weights and demonstrate how it works in various examples.


Introduction
In [1,3] a space-time independent formulation for scattering amplitudes in N = 4 SYM was proposed.
A key feature of this formulation was that scattering amplitudes were to be thought of as differential forms rather than functions and it was shown that there is a precise connection between an object called the amplituhedron living in an auxiliary space and planar scattering amplitudes in N = 4 SYM. It was further shown that unitarity and locality emerge quite naturally from geometric properties of the amplituhedron. Understanding scattering amplitudes in terms on canonical form on the amplituhedron makes the dual conformal symmetry manifest and has already produced striking results regarding positivity properties of the amplitudes as well as computational simplification for higher loop amplitudes [4,5].
Recently Arkani-Hamed et al [1] extended this formulation to non-supersymmetric domain by analysing tree-level amplitudes of bi-adjoint scalar φ 3 theory. A precise connection was established between a polytope called the associahedron and scattering amplitude of the theory. It was shown that associahedron (which is a combinatorial polytope) can be embedded as a convex polytope in the kinematic space and the canonical form associated to it is the scattering form for φ 3 theory. It was also shown that various properties like soft limits and recursion relations of scattering amplitudes could be obtained from geometric properties of the associahedron. This formulation was further extended to include 1-loop amplitudes in bi-adjoint φ 3 theory [6,7].
A natural question to ask is for what class of theories such a formulation exists. In particular for theories containing independent quartic and higher order vertices, does this picture relating scattering amplitudes to differential forms and convex polytopes continues to hold. First step in this direction was taken in [2] where it was shown that for tree-level planar amplitudes in φ 4 theory such a formulation does indeed exist. It was established in [2] that there is a relationship between the scattering amplitude, scattering forms and objects called Stokes polytope. There were however important structural differences between the amplituhedron/associahedron picture and the picture that emerged in the case of quartic interactions. Namely, for an n-particle scattering, the corresponding geometry was a union of many polytopes of a given dimension and one had to account for all of them to capture complete tree-level planar scattering amplitude.
In this paper, we extend these ideas to all φ p theories by establishing a connection between unique planar scattering forms and convex polytopes called accordiohedron. We also show that (just as in the case of cubic and quartic interactions), geometric realisation of accordiohedron inside the kinematic space leads one to identify scattering amplitudes in such theories with canonical top forms of the accordiohedron.
One of the key motivations for our work (apart from it's intrinsic interest in developing the Amplituhedron program further) was already outlined in [2] and arises from trying to better understand the CHY formulation for φ p theories [8,9]. Unlike bi-adjoint scalar field theory, or Yang-Mills theory the CHY integrands associated to φ p theories are not d-log forms and hence do not seem to have a clear geometric meaning in terms of forms on the moduli space. Our hope is that by first understanding the structure of such forms in kinematic space, we may be able to better understand CHY formula for such theories.
The rest of the paper is organised as follows. In section (2) we define the notion Q-Compatibility and use it to define the accordiohedron for φ p interactions. In section (3) we explain how to embed the accordiohedron in kinematic space and how to obtain the scattering amplitude from the accordiohedron by intorducing primitives and weights to simplify our computations. In section(4) we derive a formula for the number of primitves of a given dimension n and provide a classification all primitves up to n ≤ 3, we also provide general prescription for determining the weights. In section (5) we prove factorization for accordiohedra and also discuss the role of factorisation in determining the weights. We finally end with some conclusions and future directions.

Amplitudes and Accordiohedron
We shall begin by defining an object called the accordiohedron [10,11] associated with general dissections of polygons. We shall then show how this reduces to the associahedron and Stokes polytopes [1,2] for cubic and quartic interactions respectively. We shall then argue that the accordiohedron is thus the natural candidate for the positive geometry of all φ p interactions with p ≥ 3.

Accordion lattices and Accordiohedra
Let A be a convex polygon. Let us consider the division of A into identical p-gons which we call p-angulation of A. We can represent A as a set of points on the unit circle oriented clockwise where the arcs represent edges of A and chords represent diagonals of A. The simplest example is the case where we divide (2p−2)-gon A into two p-gons (see figure (1)). There are p−1 possible p-angulations which correspond to having the diagonals {(1, p), (2, p + 1), · · · , (p − 1, 2p − 2)}. We define a notion We can use this rule to define accordion lattices AL P p,n of dimension n associated with a reference p -angulation P 2 as follows: • We can start with any p-angulation P of a convex polygon with n diagonals.
• In the first step for each of the n diagonals, we go to the unique (2p − 2)-gon which contains it and replace it with its Q-compatible diagonal.
• In the second step for each of the n p-angulations at the end of step one we choose one of the original n − 1 diagonals and replace it with its Q-compatible diagonal as in step one.
• We repeat this till none of the original n diagonals remain in step n.
This generates a graph 3 which is the 1-skeleton of a convex polytope called the Accordiohedron [10,11], which we shall also call AC P p,n . 4 1 In [10,11] there is different definition of compatability, but these two definitions can be shown to be equivalent to each other and we shall use the definition (1) as its most suited for our purposes. We thank Alok Laddha for explaining this fact to us. 2 We consider only the case where we divide the polygon into p-gons in this paper, but accordion lattices are defined for arbitrary dissections. 3 More precisely one can define a partial order on the the set of all p-angulations that can be obtained from a particular p-angulation P, which has a maximal and a minimal element (P itself) thus making it lattice and this procedure generates the Hasse diagram of the lattice 4 Interestingly enough, there is yet another, "complimentary" family of polytopes known as graph associahedra [12] which also had deep connections with geometry of scattering amplitudes. In contrast to accordiohedron, graph associahedra can not always be obtained by considering dissections of polygons. graph associahedra is a set of polytopes which includes, associahedron, permutahedron, halohedron etc. Many of these members, e.g. permutahedron and halohedron are associated to amplitudes in bi-adjoint scalar theory with non-planar [13] and 1-loop amplitudes [7] respectively. It is intriguing that one class of polytopes helps one to move beyond tree-level and planarity in bi-adjoint φ 3 theories and the other class helps one move beyond cubic vertices.
Thus the Accordiohedra AC P p,n are a general class of polytopes which contain both associahedron and Stokes polytopes as special cases when the p-angulations corresponds to triangulations and quadrangulations respectively. The accordiohedra AC P p,n with p > 4 also retains many of the features of the Stokes polytopes we had discussed earlier in [2] including the fact that the accordiohedron AC P p,n of a given dimension n is not unique and depends on the reference p-angulation P, which is due to the fact that (1) is not an equivalence relation as The case p = 3 is special in this sense as for AC P 3,n is independent of P , as every diagonal is Q-compatible with every other diagonal and thus we could start with any triangulation P and we would generate all possible triangulations.
The accordiohedron obtained by starting with a particular p-angulation is also completely determined by the relative configuration of diagonals. 5 The n = 1 accordiohedron AC The n = 3 case the accordiohedra continue to be one of the four n = 3 Stokes polytopes with different multiplicities, i.e., AC 4,3 for all p provided P and Q have the same configuration of diagonals. We elaborate on this in section (4.4). We expect that at higher n new polytopes which are not one of the Stokes polytopes will be eventually generated.

Positive geometry for φ p interactions
We would like to show that the accordiohedron AC (P ) p,n is the positive geometry associated to φ p interactions. We shall do this by first embedding the accordiohedron into kinematic space and then showing that the canonical form of the accordiohedron when pulled back gives the right planar scattering amplitude for φ p interactions. We start by noting the following facts: • The only tree level amplitudes consistent with φ p interactions have p + (p − 2)n external legs for n + 1 vertices.
• Analogously to the cubic and quartic cases there is a 1-1 correspondence between planar tree level Feynman graphs and dissections of p + (p − 2)n-gon into p-gons.
• We also require the accordiohedron AC (P ) p,n to have dimension n, which is the number of propagators. 6 We shall first introduce some notations and variables that we shall use to describe the kinematic space.

Kinematic space
We shall primarily follow the conventions of [1] which we briefly summarise here.
The kinematic space K n of n massless momenta in d dimensions (with n < d − 1) is spanned by the Mandelstram variables: For any set of particle labels I ⊂ {1, 2, ..., n} we can define generalised Mandelstram variables as: A more convenient basis for the kinematic K n for our purposes are the planar kinematic variables, It is clear that X i i+1 = 0 and X 1n = 0 due to masslessness and momentum conservation respectively.
The mandelstram variables s ij can be written interns of planar variables X ij as The planar variables X ij can be interpreted as the the length of a diagonal between vertices i and j of an n-gon which has the momenta p i as sides.

Planar scattering form for φ p interactions
We would like to define a planar scattering form for φ p interactions. We can associate to each planar graph g with propagators X i 1 j 1 , X i 2 j 2 · · · X injn a scattering form: Thus, when we sum over all planar graphs g we have several possible scattering forms. Since we do not have a notion of projectivity except in the case of p = 3 which helps us fix a unique scattering form [1]. We can choose a particular reference graph g (equivalently a p-angulation P ) and look at only those subset of graphs which are related to this graph by a sequence of Q-flips namely all the vertices of the accordiohedron. If a graph g is related to g by an odd (even) number of Q-flips we can associate −(+) sign to it. Thus, we can define a p-angulation P dependent planar scattering form Ω P n : Since, the Q-compatible p-angulations corresponding to any reference p-angulation P does not exhaust all the p-angulations, we need to define such a planar scattering form for each P .
We now turn to embedding the accordiohedron in kinematic space and showing that when the planar scattering form is pulled back onto the accordiohedron it gives the canonical form of the accordiohedron.

Locating the accordiohedron inside kinematic space
We now define the kinematic accordiohedron AC (P ) p,n . We locate the accordiohedron inside the positive region of kinematic space X ij > 0 for all 1 ≤ i < j ≤ p + (p − 2)n by imposing the following constraints: where c ij , d r i s i are positive constants. 8 Physically we choose the above set of constraints as they do not appear as propagators of any φ p graph. The first constraint above is the famous associahedron embedding [1]. We have thus embedded the accordiohedron inside the associahedron. The positivity of X ij 's, the above constraints along with the equation (2) are a set of inequalities satisfied by the X ij which makes the convexity of the accodiohedron manifest.
We first consider the n = 1 case with the reference p-angulation to be P = {(1, p)} for p = 5, 6. The above constraints then translate to: p=5: 16 which a line with boundaries at X 16 , X 510 = 0 provided the following are satisfied In the case when p = 4 that is when the Accordiohedron is a Stokes polytope, there is a canonical choice for the additional constraints if the Stokes polytope is itself not an Associahedron [2]. However for p > 4 we do not have any canonical choice of these constraints. As we show in this section, there is at least one choice which consistently embeds the Accordiohedron in the Kinematic space. 9 Since we are slicing the associahedron using some hyperplanes X risi = d risi to get the accordiohedron these constraints tell us how the slicing should be made. For higher n we shall not state these constraints for brevity but we shall assume that they are satisfied.
The above equations define lines with Q-compatible vertices {15, 48} and {16, 510} for p = 5 and p = 6 respectively. We can trivially repeat this exercise for any other reference p-angulation P = {i, i + p}, the results of which can be obtained by taking k → k + i − 1 in the above equations. We can now pull back the scattering form onto the accordiohedron AC P p,n as: As before to get the full amplitude M n we consider a weighted sumM n of m P p,n over all P .
It is clear that M n =M n if and only if α i = 1 2 for all i = 0, · · · , p − 1. Thus, we can simplify our computation by considering a subset of p-angulations {P 1 , . . . , P I } called primitive p-angulations for which: The primitives are the class of rotationally inequivalent diagrams. Since, a rotation does not change the relative configuration of diagonals it is clear that accodiohedra remain the same for all the diagrams that belong to a primitive class and that the weights depend only on primitives. We shall say more about primitives in section (4).
The positivity of X 610 , X 411 carves out a square region in the X 15 , X 711 space. The accordiohedron in this case is also a square as we had emphasised in section (2.1).

Analysing the combinatorics of Accordiohedra
Both in [2] and in the present case, a complete computation of the amplitude from the geometry of the polytope requires determination of all the primitives of a given dimension n and computation of the corresponding weights. We shall address the problem in this section. We emphasise that this a purely combinatorial problem and hence does not depend on the construction of kinematic space accordiohedron. In sections (4.1) ,(4.2) we shall first derive formulae to count the number of primitive accordiohedra of a given dimension n. Then in sections (4.3) ,(4.4) we provide a complete classification of primitive accordiohedra for n ≤ 3 and compute the corresponding weights for any φ p interactions. Let us first consider the quartic case.

Counting primitives for the quartic case
In this section we shall address the quartic case first and provide a formula for the number of primitive Stokes polytopes p n of a given dimension n. The main result of this section is: We shall now prove this result.
We shall consider a (4 + 2n)-gon as equally spaced points a i on the circle i.e. a i = exp 2πi 4+2n with i = 0, ..., 2n + 3. The edges of the polygon correspond to arcs a i a i+1 on the circle and the diagonals of a quadrangulation correspond to chords.
There is a natural action of the dihedral group D 2n+4 on any given quadrangulation which is generated by rotation and reflections about a given diagonal. We are interested in counting primitive quadrangulations, no two of which are related to each other by a cyclic permutation, which corresponds to a rotation on the circle. Thus, it is sufficient to consider only the cyclic group Z 2n+4 for our purposes. The problem of counting primitive quadrangulations is thus equivalent to finding the number orbits of the set of all quadrangulations of a (4 + 2n)-gon under the action of the cyclic group Z 2n+4 . We shall do this by using the celebrated Burnside's lemma, which is the standard way to count the number of orbits G/G X for the action of any finite group G on a set X. It states that the number of orbits is equal to the average number of points that remain invariant when acted on by elements of G.
Thus, to count the number of primitive quadrangulations we just need to find the subset of quadrangulations that are invariant under some rotation. This problem has been addressed by [15] using the method of generating functions, but we shall take a simpler approach here following [16].
We can consider the division of (2n + 4)-gon into n + 1 quadrilaterals. We first note that the centre of the circle is left invariant by the action of the cyclic group Z 2n+4 . The centre of the circle can lie on: (1) A diameter. This can only happen when n is odd since, the relative angle between the end points a i , a j of this diameter has to π. In case (2) the diagonals can either be rotated to themselves or into each other. This can only be accomplished by rotations of π, ± π 2 and the corresponding quadrangulations are shown in the figure (3). The number of quadrantulations of type (2) depends on whether n is divisible by 2 or 4 and is given by (n+2) 2F n 2 and (n + 2)F n 4 respectively. In the case where n = 2k we can divide k into k 1 , k 2 which we call A and B in the third figure of (3) and the number of such quadrangulations would correspond to F k 1 and F k−k 1 respectively. The total number of such quadrangulations would then be k k 1 =0 F k 1 F k−k 1 and since there are n+2 2 ways we can relabel the invariant square. Using the following combinatorial identity (see appendix B).
We have (n+2) 2F n 2 invariant quadrangulations under a rotation by π. When n = 4k we have F n 4 subquadrangulations A as shown in figure (3). There are also (n + 2) ways to relabel the invariant cell and thus there are a total of (n + 2)F n 4 quadrangulations that are invariant under a rotation by ± π 2 . Thus, after also including the identity rotation which leaves all the elements invariant we get the total number of primitive quadrangulations p n is given by: We can easily check the above formula for n = 1, 2, 3 cases by using

counting primitives for φ p case
We shall now extend our analysis for the quartic case to any general p and provide a formula for the number of primitive accordiohedra of dimension n. The number of primitives p-angulations of an (p − 2)n + p-gon is the same as the number of orbits of the cyclic group Z (p−2)n+p when it acts on the set of all p-angulations. There number of such orbits can be straightforwardly computed from Burnside's lemma just as we had done in (4.1). We proceed analogously to the quartic case (4.1) by noting that the centre of the circle is invariant under any rotation and can lie: (1) On a diameter, this happens only when n is odd and leaves the p-angulation invariant under a rotation by π (see figure (5)).
(2) Inside an invariant cell, in this case we have p-angulations for every d | Gcd(p, n) which is invariant under rotation by 2π d (see figure (5) ). The total number of p-angulations of an (p − 2)n + p-gon into (n + 1) p-gons is given by the Fuss Thus, the total number of such p-angulations once we also include identity rotation is: p n = withF n,p,1 = F n,p which we shall prove in appendix B.
To determine the weights we would need to actually identify the primitive accordiohedra, then we need to find all the vertices of the accordiohedra starting with these primitives as the reference p-angulations. There is no general classification for primitives of an arbitrary dimension n to our knowledge since they grow as p n . We provide a complete classification for n ≤ 3 and give the compute the corresponding weights.

primitives and weights for n = 2 case
We would like to provide the details of the primitives and weights for φ p interactions for the n = 2 ( 2d case). In this case there are 3 vertices with 3p − 4 legs. We could try to recursively construct these graphs from the n = 1 graphs. So without loss of generality we consider only Feynman graphs in which two of the vertices lie in a line as shown in the figure (6). The 3rd vertex can then be made to lie on the central line connecting the first two vertices and it can be either above or below this line. These graphs can be denoted as (k 1 , k 2 ) such that k 1 + k 2 = p − 2, where the 3rd vertex has k 1 legs above and k 2 legs below the central line.
Since, we are only interested in primitives the graphs for which 3rd vertices are (k 1 , k 2 ) and ). We shall call them [1], [2],...,[k], where k = p 2 . The primitives are all shown in the above figure (6). We shall now show that these are the only primitives. It is clear that all the graphs above are inequivalent under cyclic permutations. As explained earlier the total number of such graphs is given by the Fuss-Catalan number We shall show that if we perform the channel sum starting with these primitives we generate all the graphs. A cyclic permutation corresponds to a clockwise rotation of the labels and every graph returns to itself after a rotation of period 3p − 4 but a graph for which the 3rd vertex is symmetric about the central line returns to itself after only half a rotation i.e has a period 3p−4 2 . The only such graph is the last graph in the case when p = 2k. Thus, the total number of graphs generated by performing sum over all the channels is : we get the total number of graphs to be (p−1)(3p−4) 2 which agrees with our results from (7). As explained in the previous subsection the accordiohedra generated by starting with a particular graph (or p-angulation ) depends only on the relative configuration of the diagonals and in this case since there are 2 diagonals the only possibilities are: 1. The diagonals meet as in [1], in this case the accordiohedron is an associahedron A n .(see (2)) 2. The diagonals do not meet as in [2],... , [k], in all these other cases the accordiohedron turns out to be a square.
We can provide a mapping between the vertices of the Stokes polytope AC P 4,2 and AC P p,2 for the n = 2 case as follows: When the two diagonals do not meet then we notice that there are several possible choices of the diagonals for p > 5, but we notice that the maximum number of i's which can appear for any choice of P is 8, since each diagonal appears twice thus there are only 4 possible ij's. We can thus identify these i's with i = 1, ..., 8 and define a mapping.
For example we provide such a mapping between the Stokes polytope corresponding to Q = {(14, 58)} and the accordiohedra corresponding to p = 5, 6 (see (2) )below: It is straightforward to find such a mapping for general p > 4 and a general choice of P .

primitives for n = 3 case
We would now try to find all primitive graphs for the n = 3 case. In this case there are 4 vertices with 4p − 6 legs. As before we can try to recursively construct primitive graphs from n = 2 case.
,if p is odd ,if p is even There are two possible ways we add the 4th vertex which could be any of the k = p 2 types in the n = 2 case : 1. We can add another central vertex either above or below the central line to a n = 2 graph.
2. We can add the vertex to any one of the external legs of the 3rd vertex. We denote these    We could also consider graphs of the type (k 1 , k 2 , k 3 ) such that k 1 + k 2 + k 3 = p − 3. Since there are n−1 r−1 non zero solutions of x 1 + ... + x r = n, in this case we have the following possibilities: • We have a graph of the type (0,0,p-3) with period 4p − 6.
• We can have a two graphs (k 1 , 0, k 3 ) and (0, k 1 , k 3 ) for each k 1 + k 3 = p − 3 with k 1 = k 3 (which are inequivalent since they are reflections of each other) and when p is odd we also have one graph (k 1 , 0, k 3 ) with k 1 = k 3 with period 4p − 6. Thus, there are • We have one graph for each (k 1 , k 2 , k 3 ) with k 1 = k 2 = k 3 with period 4p − 6. In this case we have such diagrams.
• When p = 3k we have exactly one graph with k 1 = k 2 = k 3 which has a period 4p−6 3 .
• We have two graphs (k 1 , k 2 , k 3 ) and (k 2 , k 1 , k 3 ) for each k 1 = k 2 = k 3 . In this case we have These possibilities are summarised in the table 2 below.
We can now find the total number of p-angulations by summing over all channels by multiplying columns two and three of the tables above and adding them up.
The result of this exercise turns out to be (p−1)(2p−3)(4p−5) 3 which matches with the expected Fuss-Catalan number which agrees with equation (7) 1 4(p − 2) + 1 Type of primitive number of primitives period of the primitive Since there are 3 diagonals now the relative configuration of diagonals can be of one of the following types: 1. None of the diagonals meet -in this case the corresponding accordiohedron is a cube. There are 1 + 2(p−3) 2 Figure 9: The primitive graphs of the type (k 1 , k 2 , k 3 ).
Thus the total number of primitives is : which agrees with we our general formula (7).
The accordiohedra for n = 3 we get continue to be one of the four kinds of Stokes polytopes. We could define a function from vertices of the Stokes polytopes to that of the accordiohedra as we had done in the n = 2 case to establish that this is indeed the case. We can thus continue to use the same names Lucas, Mixed etc for the n = 3 Stokes polytopes for accordiohedra as well. We expect that at sufficiently higher n, accordiohedra will be generated which do not correspond to any Stokes polytope.

Determination of the weights
In this section we shall provide a simple method to determine the weights for the general case and demonstrate the method in a few examples. We recall that we had the reduced amplitudeM n which is a weighted sum of canonical forms of all the primitive accordiohedra of a given dimension n. We would like to determine the weights such that this gives the full amplitude i.e.M n = M n .
The full amplitude M n is given by : where, the sum is over all (i 1 j 1 , ..., i n j n ) that form a complete p-angulation.
Thus, to get the full amplitude from the partial amplitude we need to impose the constraint that each n k=1 1 X i k j k appears exactly once.
But as we had emphasised before the accordiohedron depends only on the relative configuration of diagonals of the reference p-angulation which does not change under rotations and thus it is sufficient to impose these constraints for the primitive p-angulations.
where, n i p is number of times primitive i appears in the vertices of all accodiohedra, α i p are the corresponding weights.
Since, we have managed to classify all the primitives unto n = 3 we should be able to implement this straightforward procedure to get all the weights and we shall now discuss our results.
We shall first see what these conditions are for n = 2 in the p = 5, 6 cases. p = 5 : In this case there are two primitives as we had explained in the section (4.3) and we get: which can be solved to give: In this case there are 3 primitives and we get: 2α 1 6 + α 2 6 + 2α 3 6 = 1 α 1 6 + 4α 2 6 = 1 α 1 6 + 2α 3 6 = 1 which can be solved to give:

Factorization
One of the remarkable consequences of relating tree level scattering amplitudes to positive geometries like associahedron, Stokes polytope is the fact that geometric factorization implied physical factorization of scattering amplitude. This in turn implied that tree-level unitarity and locality are emergent properties of the positive geometry [1,2]. In this section we will try to argue that this is indeed the case even for planar amplitudes in massless φ p theory. We shall first argue that the geometric factorisation of accordiohedron holds and then show how this leads to the factorisation of the amplitude.
In [2] it was shown how the factorization of Stokes polytope leads to a recursions relation on α's. We shall see that even for more general Accordiohedra α's are required to satisfy analogous recursion relations. Our first assertion is the following. Given any diagonal (ij), consider all P which contains ij and the consider all the corresponding kinematic accordiohedron AC P p,n . We contend that for each accordiohedron, the corresponding facet X ij = 0 is a product of lower dimensional accordiohedra.
where P 1 and P 2 are such that P 1 ∪ P 2 ∪ (ij) = P . P 1 is the p-angulation of the polygon {i, i + 1, . . . , j} and P 2 is the p-angulation of {j, j + 1, . . . , n, . . . , i}. Now we know that, on AC P p,n any planar scattering variable X kl is a linear combination of X ij and remaining X's which constitute P . Hence in order to prove this assertion we need to show that any X kl with i ≤ k < l ≤ j can be written as a linear combination of X ij and elements of P 1 and similarly any variable in the complimentary set can be written in terms of X ij and elements of P 2 .
However this is immediate since we know from the factorization property of associahedron proven in [1] that any X kl = X ij + i<m<n<j X mn , some of these X mn ∈ P 1 and the others are constrained via X mn = d mn . This proves our assertion. Thus X ij = 0 facet factorizes into two lower dimensional accordiohedra.
Our second assertion is that the geometric factorization implies amplitude factorization of φ p theory. This assertion is based on the following two facts (For details, we refer the reader to appendix A of [1] and [4]).
(1) As the accordiohedra is a positive geometry, we know that it's canonical form satisfies the following properties satisfed by canonical form on any positive geometry A where we think of ω A as defined on the embedding space and H is any subspace in the embedding space which contains the face B.
(2) It is also known that if B = B 1 × B 2 then Thus we immediately see that where m = j − i + 1.
We thus see that residue over each accordiohedron which contains a boundary X ij → 0 factorizes into residues over lower dimensional accordiohedra. This factorization property naturally implies factorization of amplitudes as follows. Consider the n-gon with a diagonal (ij) (with i, j such that this diagonal can be part of a p-angulation). This diagonal subdivides the n-gon into a two polygons with vertices {i, . . . , j} and {j, . . . , n, 1, . . . i} respectively. By considering all the kinematic accordiohedra associated to these polygons, we can evaluate M |j−i+1| , M n+2−(|j−i+1|) which correspond to left and right sub-amplitudes respectively. This immediately implies that This proves physical factorization. We also note that, eqns. (9) and (13) imply following constraints on α's.
The left hand side of the above equation involves sum over all accordiohedra AC P p,n for which (ij) ∈ P and the right hand side involves sum over P L and P R which range over all the p-angulations of the two polygons to the left and right of the diagonal (ij) respectively.
It can be verified that in all the examples up to p = 12 and n = 3 the α P 's do indeed satisfy these constraints.
For n = 1 there is only one diagonal which we can to be (1, p). The accordiohedra is always a line as we had emphasised in section (3.3) and (1, p) appears in the vertex of exactly two of these lattices namely {1 p, p − 1 2p − 2} and {p + 1 2p + 1, 1 p}. There is only way to divide P into P L and P R and both these are trivial have α = 1. Thus, the above equation (14) gives : We could expect that the set of equations (14) would help in determining all the weights [2].
But, as we shall now show this is not the case as (14) provides too few equations. In other words the set of equations (14) provide a set of necessary but not sufficient conditions.
Let, us consider the case n = 2 for p = 5 and the diagonal (15), in this case the 11-gon gets divided into a 5-gon and 8-gon and we have weights α P L = 1 2 and α P R = 1. There are exactly 4 squares and 5 pentagons which contain the diagonal (15) in their vertices. Thus, we have We can check that for any other choice of diagonal (ij) we get the same equation. It is clear that this is not sufficient to solve for α 1 5 , α 2 5 . The solution we had obtained using our prescription in section (4.5) namely α 1 5 = 3 11 , α 2 5 = 2 11 does indeed satisfy (15).

Conclusions and Future work
Re-formulating scattering amplitudes as differential forms on positive geometries (succinctly called the Amplituhedron program) has had profound impact on how we understand Quantum field theories and how properties like unitarity and locality are a natural consequence of the positive geometries. In theories like N = 4 Super Yang Mills theory, the Amplituhedron program offers conceptual as well as striking technical advancements in the understanding of planar S-matrix. In the non supersymmetric world, these ideas were extended to bi-adjoint scalar theory in [1] where it was shown that the corresponding amplituhedron is an associahedron in Kinematic space and the canonical form on this associahedron was proportional to the scattering amplitude.
These ideas were extended from cubic to quartic interactions in [2] where the underlying positive geometry was Stokes polytope. However unlike Associahedron, which is unique (in a given dimension), there are several Stokes polytopes in any given dimension and it was shown that one had to sum over canonical forms of all such polytopes to obtain scattering amplitude of φ 4 theory.
Not all Stokes polytopes contributed equally but one had to assign different weights to each Stokes polytope. In [2] it was argued that these weights were not assigned to a given Stokes polytope but to an equivalence class of such polytopes which were related to each other by cyclic permutations and that in each such class, one could choose a representative that we called primitive. Whence the computation of scattering amplitude reduced to the problem of finding all the primitives and assigning weights to them. In this paper, continuing along the lines of [1,2] we extend the Amplituhedron program to (treelevel) planar amplitudes for massless scalar field theories with φ p interactions. We have shown that the positive geometry underlying scattering amplitudes in this theory is a class of polytopes called accordiohedron. Accordiahedron is a family of polytopes whose members include associahedron and Stokes Polytope. 11 Just as in the case of quartic interactions there exists no single accordiohedron of a given dimension n and a weighted sum of canonical forms of all the accordiohedra of a given dimension n does indeed produce the full planar amplitude. This re-affirms and generalises the result we had obtained in in the case of quartic interactions.
In [2] the problem of counting all primitives (corresponding to Stokes polytopes) of a given dimension had remained open. In this paper we fill this gap and in fact extend this notion to primitive accordioheron. We give an enumeration of the number of primitives at arbitrary dimension n and a complete classification of primitive diagrams for n ≤ 3.
We then gave a prescription to compute the weights and provided the results for the weights obtained by using our prescription for all p > 4 in n = 1, 2 dimensions and p ≤ 12 for n = 3.
Accordioheron is a very general polytope and one may wonder if they can be used to extend the Amplituhedron program to Scalar theories with mixed-vertices (e.g. theories with cubic as well as quartic interactions). It turns out that this is indeed the case. [17] There are several outstanding questions that arise out of our analysis. In [6] it was shown that the 1-loop integrand of φ 3 theory also corresponds to canonical form on a polytope which is well known in mathematical literature called Halohedron. Whether this idea can be extended to 1-loop integrand of φ p theories remains to be seen.
One of the most striking results obtained in [1] was the derivation of CHY formula for biadjoint scalar φ 3 interactions from the canonical form on kinematic space associahedron. More in detail, it was shown that the (tree-level) moduli space of CHY scattering equations admits a compactification which is nothing but an Associahedron (called worldsheet Associahedron) and the Scattering equations can be understood as diffeomorphisms between the world-sheet Associahedron and the Associahedron in Kinematic space. It was also shown that the canonical form on Kinematic space Associahedron is nothing but a push-forward of the so-called Park Taylor form which is the canonical form on world-sheet associahedron. Although CHY integrands exist for φ p interactions for p > 3 they do not admit any such geometric interactions. Our hope is that understanding of kinematic space Accordiohedron is the first step in "geometrizing" the CHY formula for φ p theories.

Acknowledgements
I would like to thank Alok Laddha for many illuminating discussions, constant guidance and various valuable comments on improving the manuscript. I would also like to thank Nemani Suyanaryana, Sujay Ashok, Nima-Arkani Hamed, Song He and Guilio Salvatori and Pinaki Banerjee for insightful discussions and constant encouragement. I would also like to thank the participants of Soft Holography conference in pune and the string group at university of Torino where a preliminary version of this work was presented.

A Counting planar diagrams
In this appendix we want to derive an explicit formula for the number a n of planar tree level Feynman diagrams with n vertices and arbitrary interactions. We begin by considering the "toy" equations of motion which are equations of motion with a source where we set the d'Alembertian and couplings to one [18]. In other words this is equivalent to setting all propagators and vertices to one. Since, the perturbative solutions to equations of motion with a source are in 1-1 correspondence with connected tree amplitudes, the solutions of these equations gives us the generating functional for the number of Feynman diagrams.
Planar diagrams in φ p theory: The toy equations for this case are We want to φ(J) from above which would be the generating function of Planar diagrams. Note that we have φ p−1 instead of φ p−1 (p−1)! since we only want to count planar diagrams which correspond to diagrams that have a fixed ordering. We can invert (16) by treating it as a formal power series and using the Lagrange series inversion formula.
If z = f (w) which is analytic about w = 0 and f (0) = 0 then w = g(z) in the neighbourhood of z = 0 with g(z) given by the power series: where, Using this for (16) with f (z) = z − z p−1 we get The number of planar trees a n with n vertices is the coefficient of J n−1 in the series above which is the Fuss-Catalan number F n,p .