<?xml version="1.0" encoding="UTF-8"?>
<article xmlns="http://specifications.silverchair.com/xsd/1/18/SCJATS-journalpublishing.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://specifications.silverchair.com/xsd/1/18/SCJATS-journalpublishing.xsd 1/18/SCJATS-journalpublishing.xsd" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">ptep</journal-id>
<journal-title-group>
<journal-title>Progress of Theoretical and Experimental Physics</journal-title>
</journal-title-group>
<issn pub-type="epub">2050-3911</issn>
<publisher>
<publisher-name>Oxford University Press</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.1093/ptep/ptz084</article-id>
<article-id pub-id-type="publisher-id">ptz084</article-id>
<article-id pub-id-type="arxiv">arXiv:1805.00262</article-id>
<article-categories>
<subj-group subj-group-type="category-toc-heading">
<subject>Papers</subject>
<subj-group subj-group-type="category-toc-heading">
<subject>Theoretical Particle Physics</subject>
</subj-group>
</subj-group>
<subj-group subj-group-type="category-journal-collection">
<subject>PTEP/B32</subject>
<subject>PTEP/E03</subject>
<subject>PTEP/E05</subject>
</subj-group>
</article-categories>
<title-group>
<article-title><inline-formula><tex-math notation="LaTeX" id="ImEquation1"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability in higher-derivative theories</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Abe</surname> <given-names>Yugo</given-names></name>
<xref ref-type="aff" rid="AFF1"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Inami</surname> <given-names>Takeo</given-names></name>
<xref ref-type="aff" rid="AFF2"/>
<xref ref-type="aff" rid="AFF3"/>
<xref ref-type="aff" rid="AFF4"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Izumi</surname> <given-names>Keisuke</given-names></name>
<xref ref-type="aff" rid="AFF5"/>
<xref ref-type="aff" rid="AFF6"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Kitamura</surname> <given-names>Tomotaka</given-names></name>
<xref ref-type="aff" rid="AFF7"/>
<xref ref-type="aff" rid="AFF8"/>
<xref ref-type="corresp" rid="COR1"/>
<email xlink:type="simple">t.k.dirac@gmail.com</email></contrib>
<contrib contrib-type="author">
<name><surname>Noumi</surname> <given-names>Toshifumi</given-names></name>
<xref ref-type="aff" rid="AFF9"/>
</contrib>
</contrib-group>
<aff id="AFF1"><italic>National Institute of Technology, Miyakonojo College, Miyakonojo 885-8567, Japan</italic></aff>
<aff id="AFF2"><italic>Theoretical Research Division, Nishina Center, RIKEN, Wako 351-0198, Japan</italic></aff>
<aff id="AFF3"><italic>Institute of Physics, VAST, Bo Ho, Hanoi, Vietnam</italic></aff>
<aff id="AFF4"><italic>Department of Physics, Sungkyunkwan University, Suwon 16419, Republic of Korea</italic></aff>
<aff id="AFF5"><italic>Kobayashi-Maskawa Institute, Nagoya University, Nagoya 464-8602, Japan</italic></aff>
<aff id="AFF6"><italic>Department of Mathematics, Nagoya University, Nagoya 464-8602, Japan</italic></aff>
<aff id="AFF7"><italic>Department of Physics, Rikkyo University, Toshima-ku, Tokyo 171-8501, Japan</italic></aff>
<aff id="AFF8"><italic>Department of Physics, Waseda University, Shinjyuku-ku, Tokyo 169-8555, Japan</italic></aff>
<aff id="AFF9"><italic>Department of Physics, Kobe University, Kobe 657-8501, Japan</italic></aff>
<author-notes>
<corresp id="COR1">E-mail: <email>t.k.dirac@gmail.com</email></corresp>
</author-notes>
<pub-date pub-type="cover">
<month>08</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection">
<day>01</day>
<month>08</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="epub" iso-8601-date="2019-08-23">
<day>23</day>
<month>08</month>
<year>2019</year>
</pub-date>
<volume>2019</volume>
<issue>8</issue>
<elocation-id>083B06</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>02</month>
<year>2019</year>
</date>
<date date-type="rev-recd">
<day>20</day>
<month>05</month>
<year>2019</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>06</month>
<year>2019</year>
</date>
</history>
<permissions>
<copyright-statement>&#x00A9; The Author(s) 2019. Published by Oxford University Press on behalf of the Physical Society of Japan.</copyright-statement>
<copyright-year>2019</copyright-year>
<license license-type="cc-by" xlink:href="http://creativecommons.org/licenses/by/4.0/">
<license-p>This is an OpenAccess article distributed under the terms of the Creative CommonsAttribution License (<ext-link ext-link-type="uri" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.</license-p>
<license-p>Funded by SCOAP<sup>3</sup></license-p>
</license>
</permissions>
<self-uri xlink:href="ptz084.pdf"/>
<abstract abstract-type="abstract"><title>Abstract</title>
<p>We investigate the relation between <inline-formula><tex-math notation="LaTeX" id="ImEquation2"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity (<inline-formula><tex-math notation="LaTeX" id="ImEquation3"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula>) and renormalizability in theories with negative-norm states. The relation has been confirmed in many field theories, including gauge theories and Einstein gravity, by analyzing the unitarity bound, which follows from the <inline-formula><tex-math notation="LaTeX" id="ImEquation4"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and the norm positivity. On the other hand, renormalizable theories with a higher-derivative kinetic term do not necessarily satisfy the unitarity bound because of the negative-norm states. In these theories it is not known whether the <inline-formula><tex-math notation="LaTeX" id="ImEquation5"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity provides a nontrivial constraint related to the renormalizability. In this paper, by relaxing the assumption of norm positivity we derive a bound on scattering amplitudes weaker than the unitarity bound, which may be used as a consistency requirement for <inline-formula><tex-math notation="LaTeX" id="ImEquation6"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity. We demonstrate in scalar field models with a higher-derivative kinetic term that the weaker bound and the renormalizability imply identical constraints.</p>
</abstract>
<kwd-group kwd-group-type="jel">
<kwd>B32</kwd>
<kwd>E03</kwd>
<kwd>E05</kwd>
</kwd-group>
<funding-group>
<award-group award-type="grant">
<funding-source><named-content content-type="funder-name">Japan Society for the Promotion of Science</named-content>
<named-content content-type="funder-identifier">10.13039/501100001691</named-content>
</funding-source>
</award-group>
<award-group award-type="grant">
<funding-source><named-content content-type="funder-name">Grants-in-Aid for Young Scientists</named-content>
</funding-source>
</award-group>
<award-group award-type="grant">
<funding-source><named-content content-type="funder-name">Japan&#x2013;Korea Bilateral Joint Research</named-content>
</funding-source>
</award-group>
</funding-group>
<counts>
<page-count count="18"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="SEC1"><title>1. Introduction</title>
<p>Unitarity and renormalizability are fundamental principles in quantum field theories (QFTs). Renormalizability means that a theory can be described by a finite number of parameters. Fixing all of these parameters by experiments, we can make nontrivial predictions. On the other hand, unitarity consists of two independent conditions: the <inline-formula><tex-math notation="LaTeX" id="ImEquation7"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity,<xref ref-type="fn" rid="FN1"><sup>1</sup></xref> i.e. <inline-formula><tex-math notation="LaTeX" id="ImEquation8"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula>, and the positivity of physical state norms (norm positivity). <inline-formula><tex-math notation="LaTeX" id="ImEquation9"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity guarantees a unique time evolution, whereas the positivity is to make the probability positive. We emphasize that <inline-formula><tex-math notation="LaTeX" id="ImEquation10"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and unitarity are carefully distinguished in this paper. Both unitarity and renormalizability are necessary to make the theory predictive, hence it is interesting to explore their possible connections.</p>
<p>Indeed, there are a variety of known examples where the unitarity bound, derived from the <inline-formula><tex-math notation="LaTeX" id="ImEquation11"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and the norm positivity, and renormalizability imply the same constraints on interactions of the theory. Historical arguments include gauge theories and Einstein gravity [<xref ref-type="bibr" rid="B3">3</xref>&#x2013;<xref ref-type="bibr" rid="B6">6</xref>].<xref ref-type="fn" rid="FN2"><sup>2</sup></xref> More recently, it was shown that this equivalence holds true in more generic QFTs such as Lifshitz-type nonrelativistic scalar field theories [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>]. If this equivalence is universal, one may expect, e.g., that quantum gravity theories consistent with the unitarity bound are automatically renormalizable. Since the unitarity bound significantly constrains the particle spectrum and interactions of quantum gravity theories (see, e.g., Ref. [<xref ref-type="bibr" rid="B9">9</xref>] for recent discussions), it would be helpful if we could use the unitarity bound as a probe of renormalizability in this way. We would like to explore this possibility by further investigating the possible equivalence.</p>
<p>Let us recall here that the models [<xref ref-type="bibr" rid="B3">3</xref>&#x2013;<xref ref-type="bibr" rid="B8">8</xref>] mentioned above do not involve negative-norm states and thus the unitarity bound directly follows from the <inline-formula><tex-math notation="LaTeX" id="ImEquation12"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity. On the other hand, one may easily show by explicit calculations that renormalizable theories with higher-derivative terms, such as quadratic gravity, do not necessarily satisfy the unitarity bound. This is because the unitarity bound does not hold in these theories due to the negative-norm states stemming from the higher derivatives. So does this mean that the unitarity bound is not helpful anymore in investigating renormalizability? To answer this question, we raise the following question: Does <inline-formula><tex-math notation="LaTeX" id="ImEquation13"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity, without the positive norm assumption, provide nontrivial constraints? In this paper we argue that <inline-formula><tex-math notation="LaTeX" id="ImEquation14"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity indeed gives a nontrivial constraint even in theories with negative-norm states. In particular, we demonstrate in a concrete model that the obtained constraint is identical to the one implied by renormalizability.</p>
<p>Among various models with negative-norm states we would like to focus on quadratic gravity, whose action contains quadratic terms in the curvature tensor on top of the Einstein&#x2013;Hilbert term. In quadratic gravity the graviton propagator is modified as <inline-formula><tex-math notation="LaTeX" id="ImEquation15"><![CDATA[${\sim} p^{-4}$]]></tex-math></inline-formula> in the UV region and thus renormalizability is achieved at the cost of norm positivity [<xref ref-type="bibr" rid="B10">10</xref>,<xref ref-type="bibr" rid="B11">11</xref>]. Even though it is not yet clear how to deal with the negative-norm states, the theory is predictive as long as the <inline-formula><tex-math notation="LaTeX" id="ImEquation16"><![CDATA[$S$]]></tex-math></inline-formula>-matrix is unitary, so that quadratic gravity could offer an interesting approach to UV completion of gravity. In the present paper we consider a scalar field model with a modified propagator as a toy model for quadratic gravity and demonstrate that <inline-formula><tex-math notation="LaTeX" id="ImEquation17"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability lead to identical constraints on this model. Our result provides new evidence for the connection between <inline-formula><tex-math notation="LaTeX" id="ImEquation18"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability. It will also be useful for a better understanding of quadratic gravity. In particular, our analysis implies that the unitarity bound is a useful probe for renormalizability in theories without negative-norm states, whereas <inline-formula><tex-math notation="LaTeX" id="ImEquation19"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity has to be used instead if the theory contains negative-norm states.</p>
<p>The rest of this paper is organized as follows. In <xref ref-type="sec" rid="SEC2">Sect. 2</xref> we introduce a scalar field model with a higher-derivative kinetic term. The scalar propagator in this model is modified as <inline-formula><tex-math notation="LaTeX" id="ImEquation20"><![CDATA[${\sim} p^{-4}$]]></tex-math></inline-formula> in the UV, so that we can think of it as a toy model for quadratic gravity. We also elaborate on the renormalizability conditions in our model based on Refs. [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>]. In <xref ref-type="sec" rid="SEC3.1">Sect. 3.1</xref>, we review the optical theorem and the unitarity bound. In particular, we argue that the lowest-order approximation of <inline-formula><tex-math notation="LaTeX" id="ImEquation21"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula> may be used to constrain the theory, just as the tree-level unitarity of high-energy scattering constrains interactions in unitary theories. The lowest-order approximation of <inline-formula><tex-math notation="LaTeX" id="ImEquation22"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula> in the scalar field model is then investigated in <xref ref-type="sec" rid="SEC3.2">Sect. 3.2</xref>. There we find that <inline-formula><tex-math notation="LaTeX" id="ImEquation23"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability imply identical constraints. A summary and discussion are given in <xref ref-type="sec" rid="SEC4">Sect. 4</xref>. Technical details are collected in the appendices.</p>
</sec>
<sec id="SEC2"><title>2. Higher-derivative scalar field theory</title>
<p>In this section we introduce a scalar field model with a higher-derivative kinetic term as the simplest example of theories with negative-norm states. Due to the fourth-order derivative term in the kinetic term of the scalar field <inline-formula><tex-math notation="LaTeX" id="ImEquation24"><![CDATA[$\phi$]]></tex-math></inline-formula>, its propagator is modified as <inline-formula><tex-math notation="LaTeX" id="ImEquation25"><![CDATA[${\sim} p^{-4}$]]></tex-math></inline-formula> at UV and thus negative-norm (ghost) states appear. This model can then be thought of as a toy model for quadratic gravity. We also elaborate on the renormalizability conditions in this higher-derivative theory based on the criterion introduced in Refs. [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>].<xref ref-type="fn" rid="FN3"><sup>3</sup></xref></p>
<sec id="SEC2.1"><title>2.1. Quadratic action</title>
<p>The quadratic action for the scalar field <inline-formula><tex-math notation="LaTeX" id="ImEquation26"><![CDATA[$\phi$]]></tex-math></inline-formula> with the fourth-order derivative in four-dimensional spacetime is given by
<disp-formula id="ptz084-M1"><label>(1)</label><tex-math notation="LaTeX" id="Equation1"><![CDATA[
\begin{eqnarray}
\mathcal{S}_2=\int d^4x \left[-\frac{1}{2}\phi(\Box-m_{1}^2)(\Box-m_{2}^2)\phi \right]\!,\label{HDSFT1}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation27"><![CDATA[$\Box=\partial^{\mu}\partial_{\mu}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation28"><![CDATA[$m_{2}>m_{1}\ (>0)$]]></tex-math></inline-formula>. The propagator of <inline-formula><tex-math notation="LaTeX" id="ImEquation29"><![CDATA[$\phi$]]></tex-math></inline-formula> is given by <inline-formula><tex-math notation="LaTeX" id="ImEquation30"><![CDATA[$\displaystyle\frac{1}{(p^2+m_1^2)(p^2+m_2^2)}$]]></tex-math></inline-formula>, which has a more convergent UV behavior <inline-formula><tex-math notation="LaTeX" id="ImEquation31"><![CDATA[${\sim} p^{-4}$]]></tex-math></inline-formula>. Also, it has poles at <inline-formula><tex-math notation="LaTeX" id="ImEquation32"><![CDATA[$p^2=-m_1^2$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation33"><![CDATA[$p^2=-m_2^2$]]></tex-math></inline-formula> representing two species of particles. In order to construct the two sets of creation and annihilation operators, it is better to rewrite the action into that of two scalars with second-order derivatives. This can be done by the following redefinition of the fields (a careful analysis using the Lagrange multiplier is given in <xref ref-type="sec" rid="SECA">Appendix A</xref>):
<disp-formula id="ptz084-M2"><label>(2)</label><tex-math notation="LaTeX" id="Equation2"><![CDATA[
\begin{eqnarray}
\psi_{1}:=\frac{(\Box-m_{2}^2)\phi}{M}\hspace{1cm}\mbox{and}\hspace{1cm}
\psi_{2}:=\frac{(\Box-m_{1}^2)\phi}{M}\hspace{1cm}
\left( \phi=\frac{1}{M}(\psi_2-\psi_1) \right)\!, \label{phipsi}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation34"><![CDATA[$M:=\sqrt{m_2^2-m_1^2}$]]></tex-math></inline-formula>. Then, the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-M1">1</xref>) becomes
<disp-formula id="ptz084-M3"><label>(3)</label><tex-math notation="LaTeX" id="Equation3"><![CDATA[
\begin{eqnarray}
\mathcal{S}_2= \int d^4x \left[\frac{1}{2}\psi_{1}(\Box-m_{1}^2)\psi_{1}-\frac{1}{2}\psi_{2}(\Box-m_{2}^2)\psi_{2}\right]\!.
\label{HDSFT2}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The kinetic term of <inline-formula><tex-math notation="LaTeX" id="ImEquation35"><![CDATA[$\psi_1$]]></tex-math></inline-formula> has a positive sign, whereas that of <inline-formula><tex-math notation="LaTeX" id="ImEquation36"><![CDATA[$\psi_2$]]></tex-math></inline-formula> has a negative sign. Thus, the one-particle state of <inline-formula><tex-math notation="LaTeX" id="ImEquation37"><![CDATA[$\psi_{2}$]]></tex-math></inline-formula> becomes a negative-norm state (see <xref ref-type="sec" rid="SECC">Appendix C</xref>). The appearance of negative-norm states is inevitable if the original action has higher derivatives. Note that <inline-formula><tex-math notation="LaTeX" id="ImEquation38"><![CDATA[$\phi$]]></tex-math></inline-formula> is dimensionless, <inline-formula><tex-math notation="LaTeX" id="ImEquation39"><![CDATA[$[\phi]=0$]]></tex-math></inline-formula>, and <inline-formula><tex-math notation="LaTeX" id="ImEquation40"><![CDATA[$\psi_{i}$]]></tex-math></inline-formula> has the canonical dimension, <inline-formula><tex-math notation="LaTeX" id="ImEquation41"><![CDATA[$[\psi_i]=1$]]></tex-math></inline-formula>, where the bracket means the mass dimension of a parameter or field. We also introduce a canonical scalar field <inline-formula><tex-math notation="LaTeX" id="ImEquation42"><![CDATA[$\sigma$]]></tex-math></inline-formula> by
<disp-formula id="ptz084-M4"><label>(4)</label><tex-math notation="LaTeX" id="Equation4"><![CDATA[
\begin{eqnarray}
\mathcal{S}_{2\sigma}=\int d^4x \left[\frac{1}{2}\sigma(\Box-m_{\sigma}^2)\sigma \right] \label{HDSFTsigma}
\end{eqnarray}
]]></tex-math></disp-formula>
to construct a nonrenormalizable coupling term in the next subsection. This field <inline-formula><tex-math notation="LaTeX" id="ImEquation43"><![CDATA[$\sigma$]]></tex-math></inline-formula> has dimension <inline-formula><tex-math notation="LaTeX" id="ImEquation44"><![CDATA[$[\sigma]=1$]]></tex-math></inline-formula>.</p>
</sec>
<sec id="SEC2.2"><title>2.2. Interaction terms</title>
<p>We then introduce interaction terms. To discuss the relation between the tree-level <inline-formula><tex-math notation="LaTeX" id="ImEquation45"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability, we consider both renormalizable and nonrenormalizable interactions, and clarify whether they are consistent with <inline-formula><tex-math notation="LaTeX" id="ImEquation46"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula> or not in the perturbative expansion. Among various interactions we focus on the marginal ones, namely <inline-formula><tex-math notation="LaTeX" id="ImEquation47"><![CDATA[$[\lambda]=0$]]></tex-math></inline-formula>. It is first because other cases can easily be inferred from the results in the marginal cases. Moreover, if the theory contains a field with a nonpositive dimension (our <inline-formula><tex-math notation="LaTeX" id="ImEquation48"><![CDATA[$\phi$]]></tex-math></inline-formula> field has <inline-formula><tex-math notation="LaTeX" id="ImEquation49"><![CDATA[$[\phi]=0$]]></tex-math></inline-formula>), the renormalizability condition is not as simple as in ordinary QFTs [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>]. As we will explain shortly, the renormalizability condition of marginal interactions does not simply follow from ordinary power-counting, but rather requires an additional condition in our setup. Therefore, we expect to see similar nontrivial conditions for the <inline-formula><tex-math notation="LaTeX" id="ImEquation50"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity, the main objective of this paper. We also concentrate on tree-level four-point scattering because it is simple, but it already provides interesting information on the unitarity and renormalizability issues. Thus, we will consider marginal four-point vertex operators.</p>
<p>Before introducing interaction terms, we briefly summarize the results in Refs. [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>]. Let us start from a general discussion of the following <inline-formula><tex-math notation="LaTeX" id="ImEquation51"><![CDATA[$n$]]></tex-math></inline-formula>th-order interaction in four-dimensional spacetime,
<disp-formula id="ptz084-M5"><label>(5)</label><tex-math notation="LaTeX" id="Equation5"><![CDATA[
\begin{eqnarray}
S_{\rm int} = \lambda \int d^{4} x \, (\partial_x^{a_1} \phi_1) \, (\partial_x^{a_2} \phi_2) \cdots (\partial_x^{a_n} \phi_n),
\label{eq:vertex}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation52"><![CDATA[$a_1,\ldots, a_n$]]></tex-math></inline-formula> are nonnegative integers. Here, <inline-formula><tex-math notation="LaTeX" id="ImEquation53"><![CDATA[$\phi_i$]]></tex-math></inline-formula> is either <inline-formula><tex-math notation="LaTeX" id="ImEquation54"><![CDATA[$\phi$]]></tex-math></inline-formula> or <inline-formula><tex-math notation="LaTeX" id="ImEquation55"><![CDATA[$\sigma$]]></tex-math></inline-formula> as introduced in the previous subsection. In particular, their kinetic terms are nondegenerate. Each <inline-formula><tex-math notation="LaTeX" id="ImEquation56"><![CDATA[$\partial_x$]]></tex-math></inline-formula> denotes spacetime derivatives <inline-formula><tex-math notation="LaTeX" id="ImEquation57"><![CDATA[$\partial_\mu$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation58"><![CDATA[$(\mu=0, \ldots, 3)$]]></tex-math></inline-formula> whose indices are contracted in a Lorentz-invariant way. In the interaction action of Eq. (<xref ref-type="disp-formula" rid="ptz084-M5">5</xref>), the dimension of the coupling constant <inline-formula><tex-math notation="LaTeX" id="ImEquation59"><![CDATA[$\lambda$]]></tex-math></inline-formula> reads
<disp-formula id="ptz084-M6"><label>(6)</label><tex-math notation="LaTeX" id="Equation6"><![CDATA[
\begin{eqnarray}
[\lambda] = 4 - \sum_{l=1}^n ( a_l + [\phi_l] ).
\label{eqlam}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The usual power-counting renormalizability (PCR) condition requires that the mass dimension of the coupling constant is nonnegative, <inline-formula><tex-math notation="LaTeX" id="ImEquation60"><![CDATA[$[\lambda]\geq0$]]></tex-math></inline-formula>. However, an infinite number of interaction terms satisfy this condition if there exists a field with a nonpositive dimension <inline-formula><tex-math notation="LaTeX" id="ImEquation61"><![CDATA[$[\phi_l]\le0$]]></tex-math></inline-formula>. In other words, an infinite number of parameters are required, and thus the theory fails to be renormalizable. For instance, let us look at
<disp-formula id="ptz084-M7"><label>(7)</label><tex-math notation="LaTeX" id="Equation7"><![CDATA[
\begin{eqnarray}
\lambda \int d^4 x \phi^4 (\Box \phi)^2,
\label{ex:ePCR}
\end{eqnarray}
]]></tex-math></disp-formula>
which satisfies the conventional PCR condition <inline-formula><tex-math notation="LaTeX" id="ImEquation62"><![CDATA[$[\lambda]\geq0$]]></tex-math></inline-formula>. This term generates <inline-formula><tex-math notation="LaTeX" id="ImEquation63"><![CDATA[$\int d^4 x \phi^n (\Box \phi)^2$]]></tex-math></inline-formula> for any <inline-formula><tex-math notation="LaTeX" id="ImEquation64"><![CDATA[$n$]]></tex-math></inline-formula> with a divergent coefficient in the one-loop effective action (for a detailed discussion, see Refs. [<xref ref-type="bibr" rid="B7">7</xref>,<xref ref-type="bibr" rid="B8">8</xref>]). This means that an infinite number of counterterms are required, and thus (7) is not renormalizable. To make the number of parameters finite, we introduce extended PCR conditions given by
<disp-formula id="ptz084-M8"><label>(8)</label><tex-math notation="LaTeX" id="Equation8"><![CDATA[
\begin{eqnarray}
\label{eq:ePCR-1}
\sum_{~l=1~}^n (a_l + [\phi_l]) &\le& 4, \\
\end{eqnarray}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M9"><label>(9)</label><tex-math notation="LaTeX" id="Equation9"><![CDATA[
\begin{eqnarray}
\sum (a_l + [\phi_l]) &<& 4
\quad
\mbox{for an arbitrary partial sum}.
\label{eq:ePCR}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The former is nothing but the usual PCR condition, while the latter additional conditions mean that the dimension of <inline-formula><tex-math notation="LaTeX" id="ImEquation65"><![CDATA[$\partial_x^{a_{i_1}} \phi_{a_{i_1}} \cdots \partial_x^{a_{i_k}} \phi_{a_{i_k}}$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation66"><![CDATA[$(k\le n-1)$]]></tex-math></inline-formula> should be less than the spacetime dimension, <inline-formula><tex-math notation="LaTeX" id="ImEquation67"><![CDATA[$4$]]></tex-math></inline-formula>. In particular, we can see that the second conditions provide nontrivial constraints only on marginal interactions in our setup, which contains a field <inline-formula><tex-math notation="LaTeX" id="ImEquation68"><![CDATA[$\phi$]]></tex-math></inline-formula> with vanishing dimension, but does not contain fields with negative dimensions. We can see that the example in (7) does not satisfy the extended PCR condition of Eq. (<xref ref-type="disp-formula" rid="ptz084-M9">9</xref>). Based on the extended PCR conditions, we introduce two marginal interaction terms.</p>
<sec id="SEC2.2.1"><title>2.2.1. Renormalizable interaction terms</title>
<p>We consider the simplest renormalizable marginal interaction term:
<disp-formula id="ptz084-M10"><label>(10)</label><tex-math notation="LaTeX" id="Equation10"><![CDATA[
\begin{eqnarray}
S_{\rm ren} = \lambda \int d^4x 
\left\{( \partial_\mu \phi)^2\right\}^2 . 
\label{4th}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>This interaction has <inline-formula><tex-math notation="LaTeX" id="ImEquation69"><![CDATA[$a_1=a_2=a_3=a_4=1$]]></tex-math></inline-formula>, and satisfies the extended PCR conditions in Eq. (<xref ref-type="disp-formula" rid="ptz084-M8">8</xref>). For the later calculations of the scattering amplitudes, it is better to rewrite it in terms of <inline-formula><tex-math notation="LaTeX" id="ImEquation70"><![CDATA[$\psi_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation71"><![CDATA[$\psi_2$]]></tex-math></inline-formula> as defined in Eqs. (<xref ref-type="disp-formula" rid="ptz084-M2">2</xref>) (the reformulation using the Lagrange multiplier in <xref ref-type="sec" rid="SECA">Appendix A</xref> gives the same result):
<disp-formula id="ptz084-M11"><label>(11)</label><tex-math notation="LaTeX" id="Equation11"><![CDATA[
\begin{eqnarray}
S_{\rm ren} & = & \alpha \int d^4x \, \left\{\big( \partial_\mu (\psi_2-\psi_1)\big)^2\right\}^2,
\label{4thpsi} \\
\end{eqnarray}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M12"><label>(12)</label><tex-math notation="LaTeX" id="Equation12"><![CDATA[
\begin{eqnarray}
\alpha & := & \frac{\lambda}{M^4}.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>This gives the four-point interaction terms for <inline-formula><tex-math notation="LaTeX" id="ImEquation72"><![CDATA[$\psi_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation73"><![CDATA[$\psi_2$]]></tex-math></inline-formula> with derivatives, and they have the same coupling constant <inline-formula><tex-math notation="LaTeX" id="ImEquation74"><![CDATA[$\alpha$]]></tex-math></inline-formula>, with <inline-formula><tex-math notation="LaTeX" id="ImEquation75"><![CDATA[$[\alpha]=-4$]]></tex-math></inline-formula>. Although <inline-formula><tex-math notation="LaTeX" id="ImEquation76"><![CDATA[$[\alpha]<0$]]></tex-math></inline-formula> does not satisfy the PCR condition, the theory turns out to be renormalizable because the common coupling <inline-formula><tex-math notation="LaTeX" id="ImEquation77"><![CDATA[$\alpha$]]></tex-math></inline-formula> brings about a cancellation among the <inline-formula><tex-math notation="LaTeX" id="ImEquation78"><![CDATA[$\psi_{1}$]]></tex-math></inline-formula> loop and the <inline-formula><tex-math notation="LaTeX" id="ImEquation79"><![CDATA[$\psi_{2}$]]></tex-math></inline-formula> loop. In the next section we use the vertex functions (a) <inline-formula><tex-math notation="LaTeX" id="ImEquation80"><![CDATA[$\psi_{1}\psi_{1}\psi_{1}\psi_{1}$]]></tex-math></inline-formula>, (b) <inline-formula><tex-math notation="LaTeX" id="ImEquation81"><![CDATA[$\psi_{1}\psi_{1}\psi_{1}\psi_{2}$]]></tex-math></inline-formula>, and (c) <inline-formula><tex-math notation="LaTeX" id="ImEquation82"><![CDATA[$\psi_{1}\psi_{1}\psi_{2}\psi_{2}$]]></tex-math></inline-formula> given by (see <xref ref-type="fig" rid="F1">Fig. 1</xref>)</p>
<fig id="F1" orientation="portrait" position="float"><label>Fig. 1.</label><caption><p>The scalar four-point vertex functions (a) <inline-formula><tex-math notation="LaTeX" id="ImEquation83"><![CDATA[$\psi_{1}\psi_{1}\psi_{1}\psi_{1}$]]></tex-math></inline-formula>, (b) <inline-formula><tex-math notation="LaTeX" id="ImEquation84"><![CDATA[$\psi_{1}\psi_{1}\psi_{1}\psi_{2}$]]></tex-math></inline-formula>, and (c) <inline-formula><tex-math notation="LaTeX" id="ImEquation85"><![CDATA[$\psi_{1}\psi_{1}\psi_{2}\psi_{2}$]]></tex-math></inline-formula>. The broken and dotted lines stand for the scalar fields <inline-formula><tex-math notation="LaTeX" id="ImEquation86"><![CDATA[$\psi_{1}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation87"><![CDATA[$\psi_{2}$]]></tex-math></inline-formula>, respectively.</p></caption>
<graphic xmlns:xlink="http://www.w3.org/1999/xlink" orientation="portrait" position="float" xlink:href="ptz084f1.tif"/></fig>
<p><disp-formula id="ptz084-M13"><label>(13)</label><tex-math notation="LaTeX" id="Equation13"><![CDATA[
\begin{equation}
\begin{split}
\mathrm{(a)}&=8i\alpha\left[\left(p_{1}\cdot p_{2}\right)\left(p_{3}\cdot p_{4}\right)+\left(p_{1}\cdot p_{3}\right)\left(p_{2}\cdot p_{4}\right)+\left(p_{1}\cdot p_{4}\right)\left(p_{2}\cdot p_{3}\right)\right]\!,\\
\mathrm{(b)}&=-8i\alpha\left[\left(p_{1}\cdot p_{2}\right)\left(p_{3}\cdot p_{4}\right)+\left(p_{1}\cdot p_{3}\right)\left(p_{2}\cdot p_{4}\right)+\left(p_{1}\cdot p_{4}\right)\left(p_{2}\cdot p_{3}\right)\right]\!,\\
\mathrm{(c)}&=8i\alpha\left[\left(p_{1}\cdot p_{2}\right)\left(p_{3}\cdot p_{4}\right)+\left(p_{1}\cdot p_{3}\right)\left(p_{2}\cdot p_{4}\right)+\left(p_{1}\cdot p_{4}\right)\left(p_{2}\cdot p_{3}\right)\right]\!,
\label{4pointvertex11XXa}
\end{split}
\end{equation}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation88"><![CDATA[$p_{i}\ (i=1, 2, 3, 4)$]]></tex-math></inline-formula> are the four-momenta of scalar fields.</p>
</sec>
<sec id="SEC2.2.2"><title>2.2.2. Non-renormalizable interaction terms</title>
<p>We wish to consider the simplest examples of nonrenormalizable marginal interaction terms.<xref ref-type="fn" rid="FN4"><sup>4</sup></xref> One is renormalizable, while the other is not:<xref ref-type="fn" rid="FN5"><sup>5</sup></xref>
<disp-formula id="ptz084-M14"><label>(14)</label><tex-math notation="LaTeX" id="Equation14"><![CDATA[
\begin{eqnarray}
\int d^4 x \phi(\partial_{\mu}\phi)^2(\Box\phi)
\hspace{1cm}\mbox{or}\hspace{1cm}
\int d^4 x \phi^2(\Box\phi)^2
\label{4DNSI}.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>However, it turns out that we have to analyze six-point scattering amplitudes to check their consistency with the <inline-formula><tex-math notation="LaTeX" id="ImEquation89"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity <inline-formula><tex-math notation="LaTeX" id="ImEquation90"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula>, essentially because these interactions contain the <inline-formula><tex-math notation="LaTeX" id="ImEquation91"><![CDATA[$\Box$]]></tex-math></inline-formula> operator (as shown in <xref ref-type="sec" rid="SECB">Appendix B</xref>). They are not suitable for our analysis based on four-point scattering amplitudes. Notice here that any marginal fourth-order interaction with four derivatives can be reformulated into the three interaction terms in Eqs. (<xref ref-type="disp-formula" rid="ptz084-M10">10</xref>) and (<xref ref-type="disp-formula" rid="ptz084-M14">14</xref>) by partial integrals. We then introduce another scalar field <inline-formula><tex-math notation="LaTeX" id="ImEquation92"><![CDATA[$\sigma$]]></tex-math></inline-formula> defined in Eq. (<xref ref-type="disp-formula" rid="ptz084-M4">4</xref>) and consider the following interaction as a simple example of nonrenormalizable marginal interactions:
<disp-formula id="ptz084-M15"><label>(15)</label><tex-math notation="LaTeX" id="Equation15"><![CDATA[
\begin{eqnarray}
S_{\rm non}=\lambda' \int d^4 x \phi^2(\partial_{\mu}\sigma)^2\label{4th-non},
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation93"><![CDATA[$[(\partial_{\mu}\sigma)^2]=4$]]></tex-math></inline-formula>, and thus this coupling term does not satisfy the extended PCR conditions in Eq. (<xref ref-type="disp-formula" rid="ptz084-M8">8</xref>). This interaction term can be expressed with <inline-formula><tex-math notation="LaTeX" id="ImEquation94"><![CDATA[$\psi_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation95"><![CDATA[$\psi_2$]]></tex-math></inline-formula> as
<disp-formula id="ptz084-M16"><label>(16)</label><tex-math notation="LaTeX" id="Equation16"><![CDATA[
\begin{eqnarray}
S_{\rm non} & = & \alpha' \int d^4 x (\psi_2-\psi_1)^2(\partial_{\mu}\sigma)^2\label{4thpsi-non}, \\
\end{eqnarray}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M17"><label>(17)</label><tex-math notation="LaTeX" id="Equation17"><![CDATA[
\begin{eqnarray}
\alpha' & := & \frac{\lambda'}{M^2}.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>In the next section we use the vertex functions <inline-formula><tex-math notation="LaTeX" id="ImEquation96"><![CDATA[$\sigma\psi_{1}\sigma\psi_{1}$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation97"><![CDATA[$\sigma\psi_{1}\sigma\psi_{2}$]]></tex-math></inline-formula> (see <xref ref-type="fig" rid="F2">Fig. 2</xref>) given by</p>
<fig id="F2" orientation="portrait" position="float"><label>Fig. 2.</label><caption><p>The scalar four-point vertex functions (d) <inline-formula><tex-math notation="LaTeX" id="ImEquation98"><![CDATA[$\sigma\psi_{1}\sigma\psi_{1}$]]></tex-math></inline-formula> and (e) <inline-formula><tex-math notation="LaTeX" id="ImEquation99"><![CDATA[$\sigma\psi_{1}\sigma\psi_{2}$]]></tex-math></inline-formula>. The solid, broken, and dotted lines stand for the scalar fields <inline-formula><tex-math notation="LaTeX" id="ImEquation100"><![CDATA[$\sigma$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation101"><![CDATA[$\psi_{1}$]]></tex-math></inline-formula>, and <inline-formula><tex-math notation="LaTeX" id="ImEquation102"><![CDATA[$\psi_{2}$]]></tex-math></inline-formula>, respectively.</p></caption>
<graphic xmlns:xlink="http://www.w3.org/1999/xlink" orientation="portrait" position="float" xlink:href="ptz084f2.tif"/></fig>
<p><disp-formula id="ptz084-M18"><label>(18)</label><tex-math notation="LaTeX" id="Equation18"><![CDATA[
\begin{equation}
\begin{split}
\mathrm{(d)}&=-4i\alpha'\left(p_{1}\cdot p_{3}\right)\!,\\
\mathrm{(e)}&=4i\alpha'\left(p_{1}\cdot p_{3}\right)\!,
\label{4pointvertexS1SX}
\end{split}
\end{equation}]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation103"><![CDATA[$p_{i}\ (i=1, 3)$]]></tex-math></inline-formula> are the four-momenta of scalar fields.</p>
</sec>
</sec>
</sec>
<sec id="SEC3"><title>3. <inline-formula><tex-math notation="LaTeX" id="ImEquation104"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability</title>
<p>We now investigate the relation between renormalizability and <inline-formula><tex-math notation="LaTeX" id="ImEquation105"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity for high-energy scattering. In theories without negative-norm states, <inline-formula><tex-math notation="LaTeX" id="ImEquation106"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity implies the unitarity bound. We know that the unitarity bound gives identical constraints to renormalizability in various theories. On the other hand, if there exist negative-norm states, the unitarity bound no longer follows from <inline-formula><tex-math notation="LaTeX" id="ImEquation107"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity. However, in this section we argue that the optical theorem (which is based on <inline-formula><tex-math notation="LaTeX" id="ImEquation108"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity, but does not rely on norm positivity) provides nontrivial consistency conditions even in theories with negative-norm states. Interestingly, we find that the constraint obtained is identical to that required by renormalizability.</p>
<p>In <xref ref-type="sec" rid="SEC3.1">Sect. 3.1</xref> we firstly see the relation between <inline-formula><tex-math notation="LaTeX" id="ImEquation109"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity, the optical theorem, and the unitarity bound by carefully reviewing the derivation of the unitarity bound and locating the step in which the positivity of all norms is required. There, we will see that the optical theorem may be used to constrain theories with negative norms.<xref ref-type="fn" rid="FN6"><sup>6</sup></xref> In <xref ref-type="sec" rid="SEC3.2">Sect. 3.2</xref> we then concretely investigate the relation between renormalizability and the optical theorem in the simple scalar field model introduced in <xref ref-type="sec" rid="SEC2">Sect. 2</xref>. In particular, we will see that the (non)renormalizable interaction term is (in)consistent with the optical theorem in the UV limit.</p>
<sec id="SEC3.1"><title>3.1. Optical theorem and unitarity bound</title>
<p>We begin with the <inline-formula><tex-math notation="LaTeX" id="ImEquation110"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity,
<disp-formula id="ptz084-M19"><label>(19)</label><tex-math notation="LaTeX" id="Equation19"><![CDATA[
\begin{eqnarray}
SS^\dagger = 1.
\label{US}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Decomposing the <inline-formula><tex-math notation="LaTeX" id="ImEquation111"><![CDATA[$S$]]></tex-math></inline-formula>-matrix as <inline-formula><tex-math notation="LaTeX" id="ImEquation112"><![CDATA[$S = 1 + i T$]]></tex-math></inline-formula>, we can rewrite the above equation as
<disp-formula id="ptz084-M20"><label>(20)</label><tex-math notation="LaTeX" id="Equation20"><![CDATA[
\begin{eqnarray}
-i(T - T^{\dagger})=TT^{\dagger}.
\label{T1}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Let us also suppose that we can take a complete orthonormal basis of the Hilbert space <inline-formula><tex-math notation="LaTeX" id="ImEquation113"><![CDATA[$\{ |{X} \rangle \}$]]></tex-math></inline-formula> as
<disp-formula id="ptz084-M21"><label>(21)</label><tex-math notation="LaTeX" id="Equation21"><![CDATA[
\begin{eqnarray}
\sum_{X}|X\rangle C_X \langle X|=1 \hspace{1cm} \mbox{with} \hspace{1cm}
\langle X|Y\rangle=C_X^{-1}\delta_{XY},
\label{defCX}
\end{eqnarray}
]]></tex-math></disp-formula>
where we have assumed that <inline-formula><tex-math notation="LaTeX" id="ImEquation114"><![CDATA[$\langle X|X\rangle\neq 0$]]></tex-math></inline-formula>, but it can be either positive or negative. Note that if the theory contains a negative-norm state, its normalization factor <inline-formula><tex-math notation="LaTeX" id="ImEquation115"><![CDATA[$C_X$]]></tex-math></inline-formula> is inevitably negative, because the sign of <inline-formula><tex-math notation="LaTeX" id="ImEquation116"><![CDATA[$C_X$]]></tex-math></inline-formula> cannot be flipped by normalization of the state <inline-formula><tex-math notation="LaTeX" id="ImEquation117"><![CDATA[$|X\rangle$]]></tex-math></inline-formula>. We would then like to reformulate the relation in Eq. (<xref ref-type="disp-formula" rid="ptz084-M20">20</xref>) in terms of the covariant scattering amplitude <inline-formula><tex-math notation="LaTeX" id="ImEquation118"><![CDATA[${\cal{M}}{\left(i \to f\right)}$]]></tex-math></inline-formula> defined as
<disp-formula id="ptz084-M22"><label>(22)</label><tex-math notation="LaTeX" id="Equation22"><![CDATA[
\begin{eqnarray}
\langle f |T| i \rangle = \delta^{4} \left( {p}_i - {p}_f \right) {\cal{M}} \left( i \to f \right) . 
\label{defM}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>By substituting Eq. (<xref ref-type="disp-formula" rid="ptz084-M22">22</xref>) into Eq. (<xref ref-type="disp-formula" rid="ptz084-M20">20</xref>), we arrive at the cutting rule in terms of covariant amplitudes,
<disp-formula id="ptz084-M23"><label>(23)</label><tex-math notation="LaTeX" id="Equation23"><![CDATA[
\begin{eqnarray}
-i[{\cal{M}}\left(i\to f\right)-{\cal{M}}\left(f\to i \right)^{\ast}] = \sum_{X}C_X\delta^{4} \left( {p} - {p}_X \right){\cal{M}}\left(i\to X\right){\cal{M}}\left(f\to X\right)^{\ast},
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation119"><![CDATA[${p}$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation120"><![CDATA[$(={p}_{i}={p}_{f})$]]></tex-math></inline-formula> is the total energy-momentum. By considering forward scattering (<inline-formula><tex-math notation="LaTeX" id="ImEquation121"><![CDATA[$| i \rangle=| f \rangle$]]></tex-math></inline-formula>) the relation is reduced to the optical theorem,
<disp-formula id="ptz084-M24"><label>(24)</label><tex-math notation="LaTeX" id="Equation24"><![CDATA[
\begin{eqnarray}
2{\rm{Im}}\,{\cal{M}}\left(i\to i\right)=\sum_{X} C_X \delta^{4} \left( {p} - {p}_X \right)|{\cal{M}}\left(i\to X\right)|^2,
\label{op}
\end{eqnarray}
]]></tex-math></disp-formula>
where the summation is over all possible intermediate on-shell states. In particular, the optical theorem implies that
<disp-formula id="ptz084-M25"><label>(25)</label><tex-math notation="LaTeX" id="Equation25"><![CDATA[
\begin{eqnarray}
|{\cal{M}}\left(i\to i\right)|\ge|{\rm{Im}}\,{\cal{M}}\left(i\to i\right)|
=\frac{1}{2}\bigg|\sum_{X} C_X \delta^{4} \left( {p} - {p}_X \right)|{\cal{M}}\left(i\to X\right)|^2\bigg|. \label{ImM}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The discussion so far is applicable even if the theory contains negative-norm states.</p>
<p>If the theory does not contain negative-norm states, the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) implies the unitarity bound. To explain this, let us consider theories without negative-norm states and set <inline-formula><tex-math notation="LaTeX" id="ImEquation122"><![CDATA[$C_X=1$]]></tex-math></inline-formula> by taking the normal basis. In this case, the right-hand side (RHS) of Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) is bounded as
<disp-formula id="ptz084-M26"><label>(26)</label><tex-math notation="LaTeX" id="Equation26"><![CDATA[
\begin{align}
\text{RHS of (25)}&=\frac{1}{2}\sum_{X}\delta^{4} \left( {p} - {p}_X \right)|{\cal{M}}\left(i\to X\right)|^2\geq\frac{1}{2}|{\cal{M}}\left(i\to j\right)|^2
\end{align}
]]></tex-math></disp-formula>
for an arbitrary on-shell intermediate state <inline-formula><tex-math notation="LaTeX" id="ImEquation123"><![CDATA[$j$]]></tex-math></inline-formula>. By setting <inline-formula><tex-math notation="LaTeX" id="ImEquation124"><![CDATA[$j=i$]]></tex-math></inline-formula>, we obtain the inequality 
<disp-formula id="ptz084-M27"><label>(27)</label><tex-math notation="LaTeX" id="Equation27"><![CDATA[
\begin{eqnarray}
|{\cal{M}}\left(i\to i\right)|\le \mathrm{const.}, \label{UB1}
\end{eqnarray}
]]></tex-math></disp-formula>
which further implies that the left-hand side of Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) is bounded from above by a constant. By combining Eqs. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>)&#x2013;(<xref ref-type="disp-formula" rid="ptz084-M27">27</xref>), we arrive at the following bound on scattering amplitudes with an arbitrary on-shell state <inline-formula><tex-math notation="LaTeX" id="ImEquation125"><![CDATA[$j$]]></tex-math></inline-formula>:
<disp-formula id="ptz084-M28"><label>(28)</label><tex-math notation="LaTeX" id="Equation28"><![CDATA[
\begin{eqnarray}
|{\cal{M}}\left(i\to j\right)|\le \mathrm{const}. \label{UB}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>This is the unitarity bound.<xref ref-type="fn" rid="FN7"><sup>7</sup></xref> In particular, this bound for high-energy scattering at the tree level (that is, the lowest order in coupling constants) is called the tree-level unitarity condition.</p>
<p>It is known that the unitarity bound in Eq. (<xref ref-type="disp-formula" rid="ptz084-M28">28</xref>), especially the tree-level unitarity, gives an identical condition to renormalizability in various theories without negative-norm states. A natural question will be if the weaker bound, Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>), which holds true even in theories with negative-norm states, may provide nontrivial constraints useful in analyzing the UV properties such as renormalizability. In the next subsection we use the scalar field model introduced in <xref ref-type="sec" rid="SEC2">Sect. 2</xref> to demonstrate that it is indeed the case. In particular, we show that tree-level high-energy scattering induced by the (non)renormalizable interaction is (in)consistent with the bound in Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) implied by the optical theorem.</p>
</sec>
<sec id="SEC3.2"><title>3.2. Weaker bound for amplitudes in higher-derivative scalar models</title>
<p>We finally come to our main question of whether the weaker unitarity bound of Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) is satisfied for amplitudes in higher-derivative scalar models of <xref ref-type="sec" rid="SEC2">Sect. 2</xref>. From the two ways of quantization for <inline-formula><tex-math notation="LaTeX" id="ImEquation126"><![CDATA[$\psi_2$]]></tex-math></inline-formula> (see <xref ref-type="sec" rid="SECC">Appendix C</xref>), we choose the one for which negative-norm states exist. <xref ref-type="fn" rid="FN8"><sup>8</sup></xref> For the reason mentioned in Sec.t <xref ref-type="sec" rid="SEC3.1">3.1</xref>, we consider two&#x2013;two scattering:
<disp-formula id="ptz084-M29"><label>(29)</label><tex-math notation="LaTeX" id="Equation29"><![CDATA[
\begin{eqnarray}
{p}_1 + {p}_2 \to {p}_3 + {p}_4 \label{1} . 
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>It is customary to work in the center-of-mass frame, in which the three momenta satisfy
<disp-formula id="ptz084-M30"><label>(30)</label><tex-math notation="LaTeX" id="Equation30"><![CDATA[
\begin{eqnarray}
{\bf p}_2 = - {\bf p}_1, 
\qquad {\bf p}_4 = - {\bf p}_3 .\label{2}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>We denote the scattering angle by <inline-formula><tex-math notation="LaTeX" id="ImEquation127"><![CDATA[$\theta$]]></tex-math></inline-formula>. We consider two&#x2013;two scattering of <inline-formula><tex-math notation="LaTeX" id="ImEquation128"><![CDATA[$\psi_{1}, \psi_{2}, \sigma$]]></tex-math></inline-formula>:
<disp-formula id="ptz084-M31a"><label>(31a)</label><tex-math notation="LaTeX" id="Equation31"><![CDATA[
\begin{align}
\psi_1({\bf p}_1)\psi_1({\bf p}_2) &\to\psi_1({\bf p}_3)\psi_1({\bf p}_4)\label{3a} , 
\\[1mm]
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M31b"><label>(31b)</label><tex-math notation="LaTeX" id="Equation32"><![CDATA[
\begin{align}
\psi_1({\bf p}_1)\psi_1({\bf p}_2) &\to\psi_1({\bf p}_3)\psi_2({\bf p}_4)\label{3b} , 
\\[1mm]
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M31c"><label>(31c)</label><tex-math notation="LaTeX" id="Equation33"><![CDATA[
\begin{align}
\psi_1({\bf p}_1)\psi_1({\bf p}_2) &\to \psi_2({\bf p}_3)\psi_2({\bf p}_4) \label{3c} , 
\\[1mm]
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M31d"><label>(31d)</label><tex-math notation="LaTeX" id="Equation34"><![CDATA[
\begin{align}
\sigma({\bf p}_1)\psi_1({\bf p}_2) &\to \sigma({\bf p}_3)\psi_1({\bf p}_4)\label{3d} , 
\\[1mm]
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M31e"><label>(31e)</label><tex-math notation="LaTeX" id="Equation35"><![CDATA[
\begin{align}
\sigma({\bf p}_1)\psi_1({\bf p}_2) &\to \sigma({\bf p}_3)\psi_2({\bf p}_4) \label{3e} . 
\end{align}
]]></tex-math></disp-formula></p>
<p>Note that (31b) and (31e) contain a single negative-norm particle <inline-formula><tex-math notation="LaTeX" id="ImEquation129"><![CDATA[$\psi_{2}$]]></tex-math></inline-formula>. The amplitudes are denoted by
<disp-formula id="ptz084-M31"><label>(32)</label><tex-math notation="LaTeX" id="Equation36"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1({\bf p}_2)\to\psi_1({\bf p}_3)\psi_1({\bf p}_4)\big) \label{4a}
\end{eqnarray}
]]></tex-math></disp-formula>
for the process in (31a), and similarly for the other four. The calculation of the amplitudes is somewhat technical and is given in <xref ref-type="sec" rid="SECE">Appendix E</xref>.</p>
<p>We first give the results for the elastic amplitudes for positive-norm particles given in Eqs. (<xref ref-type="disp-formula" rid="ptz084-ME-7">E.7</xref>) and (<xref ref-type="disp-formula" rid="ptz084-ME-15">E.15</xref>) as
<disp-formula id="ptz084-M32"><label>(33)</label><tex-math notation="LaTeX" id="Equation37"><![CDATA[
\begin{align}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1({\bf p}_2)\to\psi_1({\bf p}_3)\psi_1({\bf p}_4)\big) & = 8\alpha\big( (6+ 2\cos^2\theta)|{\bf p}_1|^4+8|{\bf p}_1|^2m_1^2+3m_1^4\big), \label{SA1} \\
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-M33"><label>(34)</label><tex-math notation="LaTeX" id="Equation38"><![CDATA[
\begin{align}
{\cal M}\big(\sigma({\bf p}_1)\psi_1({\bf p}_2)\to\sigma({\bf p}_3)\psi_1({\bf p}_4)\big) & = -4\alpha'\big( (1-\cos\theta)|{\bf p}_1|^2+m_\sigma^2\big),\label{SA2}
\end{align}
]]></tex-math></disp-formula>
where the former is the renormalizable interaction and the latter the nonrenormalizable one. We immediately note that both amplitudes diverge in the limit of <inline-formula><tex-math notation="LaTeX" id="ImEquation130"><![CDATA[$|{\bf p}_1| \to \infty$]]></tex-math></inline-formula>. Hence the unitarity bound of Eq. (<xref ref-type="disp-formula" rid="ptz084-M28">28</xref>) is apparently violated for both renormalizable and nonrenormalizable marginal interactions, Eqs. (<xref ref-type="disp-formula" rid="ptz084-M10">10</xref>) and (<xref ref-type="disp-formula" rid="ptz084-M14">14</xref>), respectively, of <xref ref-type="sec" rid="SEC2">Sect. 2</xref>. This result is not unexpected because our model contains negative-norm states.</p>
<p>We then discuss whether these two amplitudes are consistent with the weaker bound of Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>) implied by the optical theorem. Similarly to the tree-level unitarity argument, we work in the lowest-order approximation in the coupling constants. In this approximation, only the intermediate two-particle states contribute to the right-hand side of the relation in Eq. (<xref ref-type="disp-formula" rid="ptz084-M25">25</xref>):
<disp-formula id="ptz084-M34"><label>(35)</label><tex-math notation="LaTeX" id="Equation39"><![CDATA[
\begin{align}
\Big|&{\cal M}\left(\Psi_1({\bf{p}}_{1}),\Psi_2({\bf{p}}_{2}) \to \Psi_1({\bf{p}}_{1}),\Psi_2({\bf{p}}_{2})\right)\Big| \nonumber \\
& \ge \Big|\rm{Im}\,{\cal M}\left(\Psi_1({\bf{p}}_{1}),\Psi_2({\bf{p}}_{2}) \to \Psi_1({\bf{p}}_{1}),\Psi_2({\bf{p}}_{2})\right) \Big|\nonumber \\
& = \frac{1}{2}\Bigg| \sum_{{i}{j}}\int \frac{d^3{\bf p}_i}{ 2E_i} \frac{d^3{\bf p}_j}{ 2E_j} \delta^4 (p_1+p_2-p_i-p_j) (-1)^{n_{ij}} \Big| {\cal M} \left( {\Psi_1(\bf{p}}_{1}),\Psi_2({\bf{p}}_{2})\to\Psi_i({\bf{p}}_{i}),\Psi_j({\bf{p}}_{j})\right) \Big|^2\Bigg| \label{op2},
\end{align}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation131"><![CDATA[$\Psi_i({\bf p}_i)$]]></tex-math></inline-formula> means a one-particle state of <inline-formula><tex-math notation="LaTeX" id="ImEquation132"><![CDATA[$\psi_1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation133"><![CDATA[$\psi_2$]]></tex-math></inline-formula>, or <inline-formula><tex-math notation="LaTeX" id="ImEquation134"><![CDATA[$\sigma$]]></tex-math></inline-formula> with the three-dimensional momentum <inline-formula><tex-math notation="LaTeX" id="ImEquation135"><![CDATA[${\bf p}_i$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation136"><![CDATA[$n_{ij}$]]></tex-math></inline-formula> is the number of negative-norm particles <inline-formula><tex-math notation="LaTeX" id="ImEquation137"><![CDATA[$\psi_2$]]></tex-math></inline-formula> in the state <inline-formula><tex-math notation="LaTeX" id="ImEquation138"><![CDATA[$|\Psi_i({\bf{p}}_{i}),\Psi_j({\bf{p}}_{j})\rangle$]]></tex-math></inline-formula>, and we take all possible two-particle states in the sum with respect to <inline-formula><tex-math notation="LaTeX" id="ImEquation139"><![CDATA[$i$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation140"><![CDATA[$j$]]></tex-math></inline-formula>. <inline-formula><tex-math notation="LaTeX" id="ImEquation141"><![CDATA[$E_i$]]></tex-math></inline-formula> is the energy of the intermediate on-shell particles <inline-formula><tex-math notation="LaTeX" id="ImEquation142"><![CDATA[$i$]]></tex-math></inline-formula>, and similarly for other on-shell momenta. The measure of the integral is fixed by the inner product in Eq. (<xref ref-type="disp-formula" rid="ptz084-MC-8">C.8</xref>) in <xref ref-type="sec" rid="SECC">Appendix C</xref>. In the center-of-mass frame (<inline-formula><tex-math notation="LaTeX" id="ImEquation143"><![CDATA[${\bf p}_1=-{\bf p}_2$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation144"><![CDATA[${\bf p}_i=-{\bf p}_j$]]></tex-math></inline-formula>), the last line of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) can be written as (see <xref ref-type="sec" rid="SECD">Appendix D</xref>)
<disp-formula id="ptz084-M35"><label>(36)</label><tex-math notation="LaTeX" id="Equation40"><![CDATA[
\begin{eqnarray}
\propto\Bigg|\sum_{ij} (-1)^{n_{ij}} \frac{ |\hat{{\bf p}}_i|}{(E_1+E_2)} \int_0^\pi \sin \theta \Big| {\cal M} \left( \Psi_1({\bf{p}}_{1}),\Psi_2(-{\bf{p}}_{1})\to\Psi_i(\hat {\bf{p}}_{i}),\Psi_j( -\hat{\bf{p}}_{i})\right) \Big|^2 d\theta \Bigg| \label{op3},
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation145"><![CDATA[$\theta$]]></tex-math></inline-formula> is the angle between <inline-formula><tex-math notation="LaTeX" id="ImEquation146"><![CDATA[${\bf p}_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation147"><![CDATA[${\bf p}_i$]]></tex-math></inline-formula>. Also, <inline-formula><tex-math notation="LaTeX" id="ImEquation148"><![CDATA[$\hat{\bf p}_i$]]></tex-math></inline-formula> is the on-shell momentum satisfying energy-momentum conservation, so that its absolute value <inline-formula><tex-math notation="LaTeX" id="ImEquation149"><![CDATA[$|\hat{\bf p}_i|$]]></tex-math></inline-formula> is a function of <inline-formula><tex-math notation="LaTeX" id="ImEquation150"><![CDATA[$|{\bf p}_1|$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation151"><![CDATA[$\theta$]]></tex-math></inline-formula>. The proportionality constant is an unimportant numerical factor.</p>
<p>We can now compare the energy dependence of both sides in the inequality of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) in the high-energy limit. Since we saw in Eq. (<xref ref-type="disp-formula" rid="ptz084-M33">33</xref>) that <inline-formula><tex-math notation="LaTeX" id="ImEquation152"><![CDATA[${\cal M}\big(\psi_1({\bf p}_1)\psi_1({\bf p}_2)\to\psi_1({\bf p}_3)\psi_1({\bf p}_4)\big)$]]></tex-math></inline-formula> behaves as <inline-formula><tex-math notation="LaTeX" id="ImEquation153"><![CDATA[$|{\bf p}_1|^4$]]></tex-math></inline-formula> in the UV limit, naive comparison of the energy dependence in each term seems to imply violation of the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) even for the renormalizable interaction. However, we will find that negative contributions in the sum of the optical theorem, Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>), bring about a cancellation that restores the optical theorem for the renormalizable interaction. On the other hand, we will find that the cancellation for the nonrenormalizable case is not enough to restore the inequality of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>). The argument below is based on the calculation of each scattering amplitude <inline-formula><tex-math notation="LaTeX" id="ImEquation154"><![CDATA[${\cal M} \left( \Psi_1({\bf{p}}_{1}),\Psi_2(-{\bf{p}}_{1})\to\Psi_i(\hat {\bf{p}}_{i}),\Psi_j( -\hat{\bf{p}}_{i})\right)$]]></tex-math></inline-formula> in <xref ref-type="sec" rid="SECE">Appendix E</xref>.</p>
<sec id="SEC3.2.1"><title>3.2.1. Renormalizable interaction term</title>
<p>We first consider the renormalizable interaction in Eq. (<xref ref-type="disp-formula" rid="ptz084-M11">11</xref>) and investigate the inequality of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>). The calculation is simple if we take the center-of-mass initial state with two <inline-formula><tex-math notation="LaTeX" id="ImEquation155"><![CDATA[$\psi_1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation156"><![CDATA[$|i\rangle=|\psi_1,\psi_1 \rangle$]]></tex-math></inline-formula>. The left-hand side of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is obtained by setting <inline-formula><tex-math notation="LaTeX" id="ImEquation157"><![CDATA[$\cos \theta =1$]]></tex-math></inline-formula> in the amplitude in Eq. (<xref ref-type="disp-formula" rid="ptz084-M33">33</xref>),
<disp-formula id="ptz084-UM1"><tex-math notation="LaTeX" id="Equation41"><![CDATA[
\begin{equation*}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1(-{\bf p}_1)\to\psi_1({\bf p}_1)\psi_1(-{\bf p}_1)\big) = 8\alpha\big( 8|{\bf p}_1|^4+8|{\bf p}_1|^2m_1^2+3m_1^4\big)
= \alpha {\cal O}\left(|{\bf p}_1|^4\right)\!.
\end{equation*}
]]></tex-math></disp-formula></p>
<p>The right-hand side of the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is the sum of Eqs. (<xref ref-type="disp-formula" rid="ptz084-ME-8">E.8</xref>), (<xref ref-type="disp-formula" rid="ptz084-ME-11">E.11</xref>), and (<xref ref-type="disp-formula" rid="ptz084-ME-14">E.14</xref>). Here, we have to multiply Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-11">E.11</xref>) by <inline-formula><tex-math notation="LaTeX" id="ImEquation158"><![CDATA[$2$]]></tex-math></inline-formula> because there are two cases, <inline-formula><tex-math notation="LaTeX" id="ImEquation159"><![CDATA[$(\Psi_i,\Psi_j)=(\psi_1,\psi_2),(\psi_2,\psi_1)$]]></tex-math></inline-formula>. Then, it can easily be seen that the cancellation occurs in the leading order and the next-to-leading order. As a result, the right-hand side of the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is
<disp-formula id="ptz084-M36"><label>(37)</label><tex-math notation="LaTeX" id="Equation42"><![CDATA[
\begin{eqnarray}
&&64 \alpha^2\frac{m_1^4|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi {\cal O}\left( \left( \frac{m_1}{|{\bf p}_1|}\right)^0 \right) d \theta= \alpha^2 {\cal O} \left(|{\bf p}_1|^4\right)\!.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Hence, the energy dependence of both sides of the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) at high energy is the same. Therefore, the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is satisfied at high energy, provided we take the coupling constant <inline-formula><tex-math notation="LaTeX" id="ImEquation160"><![CDATA[$\alpha$]]></tex-math></inline-formula> sufficiently small.</p>
</sec>
<sec id="SEC3.2.2"><title>3.2.2. Nonrenormalizable interaction term</title>
<p>Next, we will see that the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) does not hold true for the nonrenormalizable interaction term of Eq. (<xref ref-type="disp-formula" rid="ptz084-M15">15</xref>). The left-hand side of Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is the amplitude for the forward scattering and it can be obtained by setting <inline-formula><tex-math notation="LaTeX" id="ImEquation161"><![CDATA[$\cos \theta=1$]]></tex-math></inline-formula> in Eq. (<xref ref-type="disp-formula" rid="ptz084-M34">34</xref>),
<disp-formula id="ptz084-M37"><label>(38)</label><tex-math notation="LaTeX" id="Equation43"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\sigma({\bf p}_1)\psi_1(-{\bf p}_1)\to\sigma({\bf p}_1)\psi_1(-{\bf p}_1)\big) =4\alpha' m_\sigma^2.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Note that the leading-order term in the high-energy limit, that is proportional to <inline-formula><tex-math notation="LaTeX" id="ImEquation162"><![CDATA[$\left|{\bf p}_1\right|^2$]]></tex-math></inline-formula>, disappears, and that this approaches a constant in the high-energy limit. The right-hand side is the sum of Eqs. (<xref ref-type="disp-formula" rid="ptz084-ME-16">E.16</xref>) and (<xref ref-type="disp-formula" rid="ptz084-ME-19">E.19</xref>). The leading order of the sum is cancelled, but the next-to-leading order is not. The exact form of the sum is
<disp-formula id="ptz084-UM2"><tex-math notation="LaTeX" id="Equation44"><![CDATA[
\begin{equation*}
16 {\alpha'}^2\frac{|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi \sin\theta \left[ (1-\cos\theta)^2 \frac{3M^2}{4|{\bf p}_1|^2}+ {\cal O}\left(\left(\frac{m_1}{|{\bf p}_1|}\right)^4\right) \right] d \theta = \alpha'^2 {\cal O}(|{\bf p}_1|^2).
\end{equation*}
]]></tex-math></disp-formula></p>
<p>This means that the right-hand side diverges in the high-energy limit. Therefore, the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) is not satisfied in the high-energy limit, i.e. the optical theorem is violated at some energy.</p>
</sec>
</sec>
</sec>
<sec id="SEC4"><title>4. Summary</title>
<p>We have studied the relation between <inline-formula><tex-math notation="LaTeX" id="ImEquation163"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity, <inline-formula><tex-math notation="LaTeX" id="ImEquation164"><![CDATA[$SS^\dagger =1$]]></tex-math></inline-formula>, and renormalizability in scalar field models with a higher-derivative kinetic term, which gives negative-norm states. Two cases, a theory with a renormalizable interaction and another with a nonrenormalizable interaction, have been investigated. Tree-level <inline-formula><tex-math notation="LaTeX" id="ImEquation165"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity is satisfied in the former, while it is violated in the latter. This provides evidence that <inline-formula><tex-math notation="LaTeX" id="ImEquation166"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability are related to each other even in theories with negative-norm states.</p>
<p>Combining our result with past analyses of the relation between the unitarity bound and renormalizability without negative-norm states [<xref ref-type="bibr" rid="B3">3</xref>&#x2013;<xref ref-type="bibr" rid="B8">8</xref>], the relation between <inline-formula><tex-math notation="LaTeX" id="ImEquation167"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity and renormalizability is expected to hold true in generic QFTs, regardless of whether there are negative-norm states or not. Our results therefore support an expectation that the unitarity bound is useful in investigating the renormalizability of theories with norm positivity. We believe that it is particularly powerful when applied to quantum gravity theories.</p>
<p>The relation to the renormalizability of <inline-formula><tex-math notation="LaTeX" id="ImEquation168"><![CDATA[$R_{\mu\nu}^2$]]></tex-math></inline-formula> gravity is a particularly interesting objective. <inline-formula><tex-math notation="LaTeX" id="ImEquation169"><![CDATA[$R_{\mu\nu}^2$]]></tex-math></inline-formula> gravity is known to be renormalizable. The matter scattering of canonical fields (i.e. non-ghost fields) with a graviton propagator has been confirmed to satisfy the unitarity bound [<xref ref-type="bibr" rid="B14">14</xref>]. It will be interesting to see the optical theorem for graviton scattering. The counter terms in the renormalization can be combined due to the symmetry (i.e. the general covariance), and thus the number of required parameters of the counter terms is finite. Our result suggests that <inline-formula><tex-math notation="LaTeX" id="ImEquation170"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity also requires the symmetry. It should be possible to be seen in graviton scattering. From this point of view, it will also be very interesting to analyze the graviton loops to the scalar potential in <inline-formula><tex-math notation="LaTeX" id="ImEquation171"><![CDATA[$R^{2}_{\mu\nu}$]]></tex-math></inline-formula> gravity theory. Actually, in Einstein gravity we know that gravitational corrections give rise to nonrenormalizable divergent terms [<xref ref-type="bibr" rid="B15">15</xref>&#x2013;<xref ref-type="bibr" rid="B17">17</xref>]. It is suggested that <inline-formula><tex-math notation="LaTeX" id="ImEquation172"><![CDATA[$R^{2}_{\mu\nu}$]]></tex-math></inline-formula> gravity is a useful attempt at this problem [<xref ref-type="bibr" rid="B18">18</xref>,<xref ref-type="bibr" rid="B19">19</xref>].</p>
</sec>
</body>
<back>
<ack id="ack1">
<title>Acknowledgements</title>
<p>The authors would like to thank Y.-T. Huang for valuable discussions. T.I. thanks Yoonbai Kim and S. Kawai for their kind hospitality during his one and a half year stay at Sungkunkwan University, where he benefitted much from discussions with the theory people and the students. T.K. thanks K. Maeda, S. Yamada, and S. Miyashita for valuable discussions. K.I. is supported by Japan Society for the Promotion of Science (JSPS) Grants-in-Aid for Young Scientists (B) (No. 17K14281) and for Scientific Research (A) (No. 17H01091). T.K. is supported by a Waseda University Grant-in-Aid for Special Research Projects (2017K-233). T.N. is in part supported by JSPS KAKENHI Grant Numbers JP17H02894 and JP18K13539, and MEXT KAKENHI Grant Number JP18H04352. This work is partially supported by Japan&#x2013;Korea Bilateral Joint Research Projects (JSPS-NRF collaboration) String Axion Cosmology. The authors wish to thank T. Hatsuda of RIKEN and Y. Tsuboi of Chuo University for their kind hospitality.</p>
</ack>
<sec><title>Funding</title>
<p>Open Access funding: SCOAP<inline-formula><tex-math notation="LaTeX" id="ImEquation173"><![CDATA[$^3$]]></tex-math></inline-formula>.</p>
</sec>
<app-group>
<app id="APP1"><title/>
<sec id="SECA"><title>Appendix A. Decomposition to two scalar modes</title>
<p>In <xref ref-type="sec" rid="SEC2">Sect. 2</xref> we rewrite the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-M1">1</xref>) into the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-M3">3</xref>) by the redefinition of fields in Eq. (<xref ref-type="disp-formula" rid="ptz084-M2">2</xref>). One may wonder whether the transformation involving derivatives is allowed. We show in this appendix that this redefinition is consistent with the reformulation by introducing a Lagrange multiplier and an auxiliary field.</p>
<p>We introduce the auxiliary field <inline-formula><tex-math notation="LaTeX" id="ImEquation174"><![CDATA[$\chi$]]></tex-math></inline-formula> as
<disp-formula id="ptz084-MA-1"><label>(A.1)</label><tex-math notation="LaTeX" id="Equation45"><![CDATA[
\begin{eqnarray}
\chi = \Box \phi.\label{chi}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>This can be done by introducing the Lagrange multiplier <inline-formula><tex-math notation="LaTeX" id="ImEquation175"><![CDATA[$\lambda$]]></tex-math></inline-formula>, and the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-M1">1</xref>) can be expressed as
<disp-formula id="ptz084-MA-2"><label>(A.2)</label><tex-math notation="LaTeX" id="Equation46"><![CDATA[
\begin{eqnarray}
\mathcal{S}_2=\int d^4x \left[-\frac{1}{2}(\chi-m_{1}^2\phi)(\chi-m_{2}^2\phi) -\lambda (\chi-\Box \phi)\right]\!.\label{HDSFT3}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>We can easily check that the variation of the above action with respect to <inline-formula><tex-math notation="LaTeX" id="ImEquation176"><![CDATA[$\lambda$]]></tex-math></inline-formula> gives Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-1">A.1</xref>) and that, substituting Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-1">A.1</xref>), the above action is reduced to the original action of Eq. (<xref ref-type="disp-formula" rid="ptz084-M1">1</xref>). In the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-2">A.2</xref>) <inline-formula><tex-math notation="LaTeX" id="ImEquation177"><![CDATA[$\chi$]]></tex-math></inline-formula> does not have the kinetic term, and thus the variation with respect to <inline-formula><tex-math notation="LaTeX" id="ImEquation178"><![CDATA[$\chi$]]></tex-math></inline-formula> gives a constraint equation,
<disp-formula id="ptz084-MA-3"><label>(A.3)</label><tex-math notation="LaTeX" id="Equation47"><![CDATA[
\begin{eqnarray}
\chi=\frac{1}{2}\left(m_1^2+m_2^2\right)\phi - \lambda.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Substituting this into the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-2">A.2</xref>), we have
<disp-formula id="ptz084-MA-4"><label>(A.4)</label><tex-math notation="LaTeX" id="Equation48"><![CDATA[
\begin{eqnarray}
\mathcal{S}_2=\int d^4x \left[\lambda \Box \phi + \frac{1}{2}\lambda^2 -\frac{1}{2}\left(m_1^2+m_2^2\right)\lambda \phi+\frac{1}{8}\left(m_1^2-m_2^2\right)^2 \phi^2\right]\!.\label{HDSFT4}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>This action can be diagonalized by the following field redefinition:
<disp-formula id="ptz084-MA-5"><label>(A.5)</label><tex-math notation="LaTeX" id="Equation49"><![CDATA[
\begin{eqnarray}
\lambda= -\frac{M}{2} (\psi_1+\psi_2) \hspace{1cm} \mbox{and}\hspace{1cm}
\phi=\frac{(\psi_2-\psi_1)}{M},\label{lambdaphi2}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation179"><![CDATA[$M=\sqrt{m_2^2-m_1^2}$]]></tex-math></inline-formula>. We can see that the form of <inline-formula><tex-math notation="LaTeX" id="ImEquation180"><![CDATA[$\phi$]]></tex-math></inline-formula> is the same as that in Eq. (<xref ref-type="disp-formula" rid="ptz084-M2">2</xref>). Then, the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-4">A.4</xref>) becomes
<disp-formula id="ptz084-MA-6"><label>(A.6)</label><tex-math notation="LaTeX" id="Equation50"><![CDATA[
\begin{eqnarray}
\mathcal{S}_2= \int d^4x \left[\frac{1}{2}\psi_{1}(\Box-m_{1}^2)\psi_{1}-\frac{1}{2}\psi_{2}(\Box-m_{2}^2)\psi_{2}\right]\!, \label{HDSFT5}
\end{eqnarray}
]]></tex-math></disp-formula>
which is the same as the action in Eq. (<xref ref-type="disp-formula" rid="ptz084-M3">3</xref>).</p>
<p>For the interaction term <inline-formula><tex-math notation="LaTeX" id="ImEquation181"><![CDATA[${\mathcal L}_{\rm int} (\partial, \phi)$]]></tex-math></inline-formula>, we can substitute Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-5">A.5</xref>). The reason is as follows. When we introduce the auxiliary field <inline-formula><tex-math notation="LaTeX" id="ImEquation182"><![CDATA[$\chi$]]></tex-math></inline-formula>, we can leave <inline-formula><tex-math notation="LaTeX" id="ImEquation183"><![CDATA[$\phi$]]></tex-math></inline-formula> in the interaction term untouched. Then, the elimination of <inline-formula><tex-math notation="LaTeX" id="ImEquation184"><![CDATA[$\chi$]]></tex-math></inline-formula> can be done in the same way and we do the field redefinition of Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-5">A.5</xref>). This is nothing but substituting Eq. (<xref ref-type="disp-formula" rid="ptz084-MA-5">A.5</xref>).</p>
</sec>
<sec id="SECB"><title>Appendix B. Interaction with <inline-formula><tex-math notation="LaTeX" id="ImEquation185"><![CDATA[$\Box$]]></tex-math></inline-formula> operator</title>
<p>In this paper we discuss the relation between renormalizability and <inline-formula><tex-math notation="LaTeX" id="ImEquation186"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity. However, in the latter analysis, if the interaction term has <inline-formula><tex-math notation="LaTeX" id="ImEquation187"><![CDATA[$\Box$]]></tex-math></inline-formula> operators attached to external lines then the amplitude is suppressed. Therefore, in this order, we cannot see the relation. If we go to the higher-point amplitude we would see the correspondence, but calculation of the higher-point amplitude is involved. We briefly show why <inline-formula><tex-math notation="LaTeX" id="ImEquation188"><![CDATA[$\Box$]]></tex-math></inline-formula> operators make the discussion different.</p>
<p>We naively expect that the <inline-formula><tex-math notation="LaTeX" id="ImEquation189"><![CDATA[$\Box$]]></tex-math></inline-formula> operator gives the <inline-formula><tex-math notation="LaTeX" id="ImEquation190"><![CDATA[${\cal O}(E^2)$]]></tex-math></inline-formula> contribution. However, if the <inline-formula><tex-math notation="LaTeX" id="ImEquation191"><![CDATA[$\Box$]]></tex-math></inline-formula> operator acts on an external line, it becomes <inline-formula><tex-math notation="LaTeX" id="ImEquation192"><![CDATA[$m^2$]]></tex-math></inline-formula> <inline-formula><tex-math notation="LaTeX" id="ImEquation193"><![CDATA[$(={\cal O}(E^0))$]]></tex-math></inline-formula> by the on-shell condition; that is, the UV behavior becomes mild. This is also true even in the canonical theory (without negative-norm states). For instance, we consider the following action for a canonical scalar field <inline-formula><tex-math notation="LaTeX" id="ImEquation194"><![CDATA[$\zeta$]]></tex-math></inline-formula>:
<disp-formula id="ptz084-MB-1"><label>(B.1)</label><tex-math notation="LaTeX" id="Equation51"><![CDATA[
\begin{eqnarray}
{\mathcal S} = \int d^4x \left[ \frac{1}{2} \zeta \left( \Box-m^2\right) \zeta + \alpha_\zeta \zeta^3 \Box \zeta \right]\!.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Based on the power-counting argument, the interaction term is not renormalizable. Meanwhile, the tree-level four-point scattering amplitude is constant, and thus it satisfies the unitarity bound. To check the violation of the unitarity bound, we have to consider a tree diagram that has internal lines where the <inline-formula><tex-math notation="LaTeX" id="ImEquation195"><![CDATA[$\Box$]]></tex-math></inline-formula> operator is assigned. If we consider the six-point amplitude an internal line appears, and thus it is expected to violate the unitarity bound.</p>
<p>The situation is the same in the case with higher derivatives that we analyze in this paper. We have seven kinds of four-point marginal vertices. Nevertheless, they are related to the integration by parts, and actually all of them can be expressed by the following three vertex functions:
<disp-formula id="ptz084-MB-2"><label>(B.2)</label><tex-math notation="LaTeX" id="Equation52"><![CDATA[
\begin{eqnarray}
(\Box^2 \phi)\phi^3, \hspace{1cm} (\Box \phi)^2 \phi^2, \hspace{1cm} \left((\partial_\mu \phi)^2\right)^2 \label{4pt}.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The last term is considered in this paper as a renormalizable vertex function. The other two are nonrenormalizable vertex terms, but because of the <inline-formula><tex-math notation="LaTeX" id="ImEquation196"><![CDATA[$\Box$]]></tex-math></inline-formula> operator we have to go the six-point scattering amplitude to see violation of the <inline-formula><tex-math notation="LaTeX" id="ImEquation197"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity.</p>
<p>Another possible way to see the correspondence for a nonrenormalizable operator would be considering a three-point marginal vertex term and constructing the four-point scattering amplitude with two vertices. However, the situation is the same as the case with the four-point vertex. All three-point marginal operators can be expressed by
<disp-formula id="ptz084-MB-3"><label>(B.3)</label><tex-math notation="LaTeX" id="Equation53"><![CDATA[
\begin{eqnarray}
(\Box \phi)^2 \phi, \hspace{1cm} (\Box \phi)(\partial_\mu \phi)^2. \label{3pt}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The nonrenormalizable term is the former, but the <inline-formula><tex-math notation="LaTeX" id="ImEquation198"><![CDATA[$\Box$]]></tex-math></inline-formula> operator operates on a different <inline-formula><tex-math notation="LaTeX" id="ImEquation199"><![CDATA[$\phi$]]></tex-math></inline-formula>. In each vertex of tree-level four-point amplitudes only one line is internal, and thus we still have one <inline-formula><tex-math notation="LaTeX" id="ImEquation200"><![CDATA[$\Box$]]></tex-math></inline-formula> operator in the external line. Although it would be interesting to check the six-point amplitude, the calculation is complicated. Meanwhile, we can see the relation by the four-point vertex in the simple model that we introduce in <xref ref-type="sec" rid="SEC2">Sect. 2</xref>. Therefore, we do not analyze the nonrenormalizable interaction terms appearing in (B.3) and (B.2).</p>
</sec>
<sec id="SECC"><title>Appendix C. Quantization</title>
<p>We show the quantization for the fields <inline-formula><tex-math notation="LaTeX" id="ImEquation201"><![CDATA[$\psi_1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation202"><![CDATA[$\psi_2$]]></tex-math></inline-formula>, and <inline-formula><tex-math notation="LaTeX" id="ImEquation203"><![CDATA[$\sigma$]]></tex-math></inline-formula>. The quantization for the canonical fields <inline-formula><tex-math notation="LaTeX" id="ImEquation204"><![CDATA[$\psi_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation205"><![CDATA[$\sigma$]]></tex-math></inline-formula> is done in the usual manner:
<disp-formula id="ptz084-MC-1"><label>(C.1)</label><tex-math notation="LaTeX" id="Equation54"><![CDATA[
\begin{align}
\psi_1 & = \int \frac{d^3p}{\sqrt{(2\pi)^3 E^{(1)}}} \left\{ a_1({\bf p}) e^{ipx}+a_1^\dagger({\bf p}) e^{-ipx} \right\}\!,
\quad \left( E^{(1)}:=\sqrt{{\bf p}^2 +m_1^2}(=-p^{(1)}_0) \right)\!, \\
\end{align}
]]></tex-math></disp-formula>
<disp-formula id="ptz084-MC-2"><label>(C.2)</label><tex-math notation="LaTeX" id="Equation55"><![CDATA[
\begin{align}
\sigma & = \int \frac{d^3p}{\sqrt{(2\pi)^3 E^{(\sigma)}}} \left\{ a_\sigma({\bf p}) e^{ipx}+a_\sigma^\dagger({\bf p}) e^{-ipx} \right\}\!, 
\quad \left( E^{(\sigma)}:=\sqrt{{\bf p}^2 +m_\sigma^2}(=-p^{(\sigma)}_0) \right)\!.
\end{align}
]]></tex-math></disp-formula></p>
<p>The commutation relations of the creation and annihilation operators for the positive-norm fields <inline-formula><tex-math notation="LaTeX" id="ImEquation206"><![CDATA[$\psi_1$]]></tex-math></inline-formula> and <inline-formula><tex-math notation="LaTeX" id="ImEquation207"><![CDATA[$\sigma$]]></tex-math></inline-formula> are
<disp-formula id="ptz084-MC-3"><label>(C.3)</label><tex-math notation="LaTeX" id="Equation56"><![CDATA[
\begin{eqnarray}
\left[ a_1({\bf p}), a_1^\dagger({\bf k}) \right]= \delta^3({\bf p}-{\bf k}) ,
\hspace{1cm}
\left[ a_\sigma({\bf p}), a_\sigma^\dagger({\bf k}) \right]= \delta^3({\bf p}-{\bf k}) .
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The quantization of <inline-formula><tex-math notation="LaTeX" id="ImEquation208"><![CDATA[$\psi_2$]]></tex-math></inline-formula> is unsettled due to the negative sign of the kinetic term. We have two natural choices to define creation and annihilation operators for <inline-formula><tex-math notation="LaTeX" id="ImEquation209"><![CDATA[$\psi_2$]]></tex-math></inline-formula>:
<disp-formula id="ptz084-MC-4"><label>(C.4)</label><tex-math notation="LaTeX" id="Equation57"><![CDATA[
\begin{eqnarray}
\psi_2 = \int \frac{d^3p}{\sqrt{(2\pi)^3 E^{(2)}}} \left\{ a_2({\bf p}) e^{ipx}+a_2^\dagger({\bf p}) e^{-ipx} \right\}\!, 
\quad \left( E^{(2)}:=\sqrt{{\bf p}^2 +m_2^2}(=-p^{(2)}_0) \right) 
\end{eqnarray}
]]></tex-math></disp-formula>
or
<disp-formula id="ptz084-MC-5"><label>(C.5)</label><tex-math notation="LaTeX" id="Equation58"><![CDATA[
\begin{eqnarray}
\psi_2 = \int \frac{d^3p}{\sqrt{(2\pi)^3 E^{(2)}}} \left\{ \tilde a_2^\dagger ({\bf p}) e^{ipx}+ \tilde a_2({\bf p}) e^{-ipx} \right\}\!, 
\quad \left( E^{(2)}:=\sqrt{{\bf p}^2 +m_2^2}(=-p^{(2)}_0) \right)\!.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>With the former definition, the creation operator creates a positive-energy particle, but the sign of the commutation relation is flipped, <inline-formula><tex-math notation="LaTeX" id="ImEquation210"><![CDATA[$\big[ a_2({\bf p}), a_2^\dagger({\bf k}) \big] = -\delta^3({\bf p}-{\bf k})$]]></tex-math></inline-formula>. On the other hand, the latter definition gives the usual commutation relation, <inline-formula><tex-math notation="LaTeX" id="ImEquation211"><![CDATA[$\big[ \tilde a_2({\bf p}), \tilde a_2^\dagger({\bf k}) \big]= \delta^3({\bf p}-{\bf k})$]]></tex-math></inline-formula>, but particle states have negative energy. We often choose the latter definition for physical fields (especially in cosmology), because the flipped sign of the kinetic term would mean negativity of the kinetic energy and any state satisfies positivity of the norm. However, the negative energy leads to instability of the vacuum state (the ghost instability), and its decay rate is indeed infinity! Thus, the perturbative approach breaks down. <xref ref-type="fn" rid="FN9"><sup>9</sup></xref> The renormalizability of the theory with higher-order derivative is discussed with the former definition. Although it is unclear how to interpret the negative-norm state, without the negative-energy state the perturbative approach is applicable because of the absence of the ghost instability. Therefore, we choose the former definition here.</p>
<p>The particle states are usually expressed by the covariant creation operators,
<disp-formula id="ptz084-MC-6"><label>(C.6)</label><tex-math notation="LaTeX" id="Equation59"><![CDATA[
\begin{eqnarray}
A^\dagger_i ({\bf p}) = a^\dagger_i ({\bf p}) (2\pi)^{3/2} \sqrt{2 E^{(i)} ({\bf k})} \hspace{1cm} \left( i=1,2,\sigma \right)\!,
\end{eqnarray}
]]></tex-math></disp-formula>
as
<disp-formula id="ptz084-UM3"><tex-math notation="LaTeX" id="Equation60"><![CDATA[
\[\begin{align}
& |{{\mathbf{p}}_{{{1}_{1}}}},\ldots ,{{\mathbf{p}}_{{{1}_{l}}}};{{\mathbf{p}}_{{{2}_{1}}}},\ldots ,{{\mathbf{p}}_{{{2}_{m}}}};{{\mathbf{p}}_{{{\sigma }_{1}}}},\ldots ,{{\mathbf{p}}_{{{\sigma }_{n}}}}\rangle := \\ 
 & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,A_{1}^{\dagger }({{\mathbf{p}}_{{{1}_{1}}}})\cdots A_{1}^{\dagger }({{\mathbf{p}}_{{{1}_{l}}}})A_{2}^{\dagger }({{\mathbf{p}}_{{{2}_{1}}}})\cdots A_{2}^{\dagger }({{\mathbf{p}}_{{{2}_{m}}}})A_{\sigma }^{\dagger }({{\mathbf{p}}_{{{\sigma }_{1}}}})\cdots A_{\sigma }^{\dagger }({{\mathbf{p}}_{{{\sigma }_{m}}}})|0\rangle . \\ 
\end{align}\]
]]></tex-math></disp-formula></p>
<p>Based on the covariant creation operators <inline-formula><tex-math notation="LaTeX" id="ImEquation212"><![CDATA[$A^\dagger_i ({\bf p})$]]></tex-math></inline-formula>, the usual Feynman rule gives the scattering amplitude. Then, the norm of the one-particle state becomes
<disp-formula id="ptz084-MC-7"><label>(C.7)</label><tex-math notation="LaTeX" id="Equation61"><![CDATA[
\begin{eqnarray}
\langle 0| A_i ({\bf k}) A_j^\dagger ({\bf p})|0 \rangle = \eta_i (2\pi)^3 2 E^{(i)} \delta^3({\bf k}-{\bf p}) \delta_{ij}, \label{innerP}
\end{eqnarray}
]]></tex-math></disp-formula>
with
<disp-formula id="ptz084-MC-8"><label>(C.8)</label><tex-math notation="LaTeX" id="Equation62"><![CDATA[
\begin{eqnarray}
\eta_1=1=\eta_\sigma, \hspace{1cm} \eta_2=-1.
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The normalization factor <inline-formula><tex-math notation="LaTeX" id="ImEquation213"><![CDATA[$C_X$]]></tex-math></inline-formula> for particle states defined in Eq. (<xref ref-type="disp-formula" rid="ptz084-M21">21</xref>) is fixed by Eq. (<xref ref-type="disp-formula" rid="ptz084-MC-8">C.8</xref>).</p>
</sec>
<sec id="SECD"><title>Appendix D. Derivation of Eq. (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>)</title>
<p>Here, we show the derivation of Eq. (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) from the right-hand side of the inequality in Eq. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>). Three of the four delta functions represent the momentum conservation law, and there each component of <inline-formula><tex-math notation="LaTeX" id="ImEquation214"><![CDATA[${\bf p}_j$]]></tex-math></inline-formula> appears linearly. Therefore, the integration with respect to <inline-formula><tex-math notation="LaTeX" id="ImEquation215"><![CDATA[${\bf p}_j$]]></tex-math></inline-formula> can be done easily and in the center-of-mass frame it gives the replacement of <inline-formula><tex-math notation="LaTeX" id="ImEquation216"><![CDATA[${\bf p}_j$]]></tex-math></inline-formula> with <inline-formula><tex-math notation="LaTeX" id="ImEquation217"><![CDATA[$-{\bf p}_i$]]></tex-math></inline-formula> in the integrand. The remaining delta function, namely the energy conservation, can be deformed as follows: 
<disp-formula id="ptz084-MD-1"><label>(D.1)</label><tex-math notation="LaTeX" id="Equation63"><![CDATA[
\begin{eqnarray}
\delta(E_1+E_2-E_i-E_j)& =&\delta\left(\sqrt{|{\bf p}_1|^2 + \mu_1^2}+\sqrt{|{\bf p}_1|^2 + \mu_2^2}-\sqrt{|{\bf p}_i|^2 + \mu_i^2}-\sqrt{|{\bf p}_i|^2 + \mu_j^2}\right) \nonumber\\
&=&
\frac{E_iE_j}{|{\bf p}_1|(E_i+E_j)} \delta \big(|{\bf p}_i|-|\hat{{\bf p}}_i| (|{\bf p}_1|)\big)\nonumber\\
&=&\frac{E_iE_j}{|{\bf p}_i|(E_1+E_2)} \delta \big(|{\bf p}_i|-|\hat{{\bf p}}_i|(|{\bf p}_1|)\big),\label{EC}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation218"><![CDATA[$\mu_k$]]></tex-math></inline-formula> (<inline-formula><tex-math notation="LaTeX" id="ImEquation219"><![CDATA[$k=1,2,i,j$]]></tex-math></inline-formula>) shows the mass corresponding to the particle with momentum <inline-formula><tex-math notation="LaTeX" id="ImEquation220"><![CDATA[${\bf p}_k$]]></tex-math></inline-formula>. Taking spherical coordinates for the momentum space of <inline-formula><tex-math notation="LaTeX" id="ImEquation221"><![CDATA[${\bf p}_i$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation222"><![CDATA[$d^3{\bf p}_i$]]></tex-math></inline-formula> is written as <inline-formula><tex-math notation="LaTeX" id="ImEquation223"><![CDATA[$|{\bf p}_i|^2 \sin \theta_p d|{\bf p}_i| d\theta_p d \phi_p$]]></tex-math></inline-formula>. Substituting this and Eq. (<xref ref-type="disp-formula" rid="ptz084-MD-1">D.1</xref>), integrating with respect to <inline-formula><tex-math notation="LaTeX" id="ImEquation224"><![CDATA[$\phi_p$]]></tex-math></inline-formula>, and ignoring the unimportant numerical factor, we can obtain the form in Eq. (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>).</p>
</sec>
<sec id="SECE"><title>Appendix E. Calculation of the scattering amplitudes</title>
<p>In <xref ref-type="sec" rid="SEC3.2">Sect. 3.2</xref> we considered five amplitudes of two&#x2013;two scattering of <inline-formula><tex-math notation="LaTeX" id="ImEquation225"><![CDATA[$\psi_1$]]></tex-math></inline-formula>, <inline-formula><tex-math notation="LaTeX" id="ImEquation226"><![CDATA[$\psi_2$]]></tex-math></inline-formula>, and <inline-formula><tex-math notation="LaTeX" id="ImEquation227"><![CDATA[$\sigma$]]></tex-math></inline-formula>: (31a)&#x2013;(31e). The momentum assignment of scattering is as
<disp-formula id="ptz084-ME-1"><label>(E.1)</label><tex-math notation="LaTeX" id="Equation64"><![CDATA[
\begin{eqnarray}
{p}_1 + {p}_2 \to {p}_3 + {p}_4 \label{E1}
\end{eqnarray}
]]></tex-math></disp-formula>
for all five processes. Here, <inline-formula><tex-math notation="LaTeX" id="ImEquation228"><![CDATA[$p_i$]]></tex-math></inline-formula> are four-momenta. The amplitudes are denoted by
<disp-formula id="ptz084-ME-2"><label>(E.2)</label><tex-math notation="LaTeX" id="Equation65"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1({\bf p}_2)\to\psi_1({\bf p}_3)\psi_1({\bf p}_4)\big) \label{E2}
\end{eqnarray}
]]></tex-math></disp-formula>
for the process in (31a), and similarly for the other four. After taking account of the momentum conservation we write
<disp-formula id="ptz084-ME-3"><label>(E.3)</label><tex-math notation="LaTeX" id="Equation66"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1({\bf p}_2)\to\psi_1({\hat{\bf p}}_3)\psi_1({\hat{\bf p}}_4)\big) \label{E3} ,
\end{eqnarray}
]]></tex-math></disp-formula>
and similarly for the other processes. Here, <inline-formula><tex-math notation="LaTeX" id="ImEquation229"><![CDATA[$\hat{\bf p}_i$]]></tex-math></inline-formula> is the momentum satisfying both the on-shell condition and energy-momentum conservation. The momentum variables <inline-formula><tex-math notation="LaTeX" id="ImEquation230"><![CDATA[${\bf p}_1, \ldots , {\bf p}_4$]]></tex-math></inline-formula> are often omitted when it is understood. Throughout this paper we take the center-of-mass frame, hence
<disp-formula id="ptz084-ME-4"><label>(E.4)</label><tex-math notation="LaTeX" id="Equation67"><![CDATA[
\begin{eqnarray}
{\bf p}_2 = - {\bf p}_1, 
\qquad {\bf p}_4 = - {\bf p}_3 .\label{E4}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>The scattering angle is denoted by <inline-formula><tex-math notation="LaTeX" id="ImEquation231"><![CDATA[$\theta$]]></tex-math></inline-formula>. <inline-formula><tex-math notation="LaTeX" id="ImEquation232"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> is expressed in terms of <inline-formula><tex-math notation="LaTeX" id="ImEquation233"><![CDATA[$|{\bf p}_1|$]]></tex-math></inline-formula>, but its expression is different depending on the processes in (31a)&#x2013;(31e). The first three processes arise from the interaction term in Eq. (<xref ref-type="disp-formula" rid="ptz084-M14">14</xref>) and hence their amplitudes are the same except for the all-over sign <inline-formula><tex-math notation="LaTeX" id="ImEquation234"><![CDATA[$\epsilon$]]></tex-math></inline-formula>,
<disp-formula id="ptz084-ME-5"><label>(E.5)</label><tex-math notation="LaTeX" id="Equation68"><![CDATA[
\begin{eqnarray}
{\cal M}=8 \epsilon \alpha\left[\left(p_{1}\cdot p_{2}\right)\left(p_{3}\cdot p_{4}\right)+\left(p_{1}\cdot p_{3}\right)\left(p_{2}\cdot p_{4}\right)+\left(p_{1}\cdot p_{4}\right)\left(p_{2}\cdot p_{3}\right)\right]\!,\label{E5}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation235"><![CDATA[$\epsilon =1$]]></tex-math></inline-formula> if the final state gives a positive norm, otherwise <inline-formula><tex-math notation="LaTeX" id="ImEquation236"><![CDATA[$\epsilon = -1$]]></tex-math></inline-formula>. The last two processes arise from the interaction term in Eq. (<xref ref-type="disp-formula" rid="ptz084-M16">16</xref>), and their amplitudes are given by
<disp-formula id="ptz084-ME-6"><label>(E.6)</label><tex-math notation="LaTeX" id="Equation69"><![CDATA[
\begin{eqnarray}
{\cal M}=4 \epsilon \alpha'\left(p_{1}\cdot p_{3}\right)\!,\label{E6}
\end{eqnarray}
]]></tex-math></disp-formula>
where <inline-formula><tex-math notation="LaTeX" id="ImEquation237"><![CDATA[$\epsilon =1$]]></tex-math></inline-formula> if the final state gives a positive norm, otherwise <inline-formula><tex-math notation="LaTeX" id="ImEquation238"><![CDATA[$\epsilon = -1$]]></tex-math></inline-formula>. However, the sign <inline-formula><tex-math notation="LaTeX" id="ImEquation239"><![CDATA[$\epsilon$]]></tex-math></inline-formula> is immaterial in computing Eqs. (<xref ref-type="disp-formula" rid="ptz084-M35">35</xref>) and (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>), since the absolute values <inline-formula><tex-math notation="LaTeX" id="ImEquation240"><![CDATA[$|{\cal M}|^2$]]></tex-math></inline-formula> are used there.</p>
<sec id="SECE.1"><title>E.1. <inline-formula><tex-math notation="LaTeX" id="ImEquation241"><![CDATA[$\psi_1+\psi_1\to\psi_1+\psi_1$]]></tex-math></inline-formula></title>
<p>The two particles of the initial state and the two in the final state all have the same mass <inline-formula><tex-math notation="LaTeX" id="ImEquation242"><![CDATA[$m_1$]]></tex-math></inline-formula> and hence <inline-formula><tex-math notation="LaTeX" id="ImEquation243"><![CDATA[$|\hat{{\bf p}}_i|=|{\bf p}_1|$]]></tex-math></inline-formula>. Using this in the <inline-formula><tex-math notation="LaTeX" id="ImEquation244"><![CDATA[${\cal{M}}$]]></tex-math></inline-formula> of Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-5">E.5</xref>), we have
<disp-formula id="ptz084-ME-7"><label>(E.7)</label><tex-math notation="LaTeX" id="Equation70"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1(-{\bf p}_2)\to\psi_1(\hat {\bf p}_i)\psi_1(-\hat {\bf p}_i)\big)
=8\alpha\big( (6+ 2\cos^2\theta)|{\bf p}_1|^4+8m_1^2|{\bf p}_1|^2+3m_1^4\big) . 
\label{E7}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) of <xref ref-type="sec" rid="SEC3.2">Sect. 3.2</xref> now takes the form
<disp-formula id="ptz084-ME-8"><label>(E.8)</label><tex-math notation="LaTeX" id="Equation71"><![CDATA[
\begin{equation}
64 \alpha^2\frac{m_1^4|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi \sin\theta (6+2\cos^2 \theta)^2 \left[ \left(\frac{|{\bf p}_1|}{m_1}\right)^4 + \frac{16}{6+2\cos^2\theta}\left( \frac{|{\bf p}_1|}{m_1}\right)^2 + {\cal O}\left( \left( \frac{m_1}{|{\bf p}_1|}\right)^0 \right) \right] d \theta.\label{1111}
\end{equation}]]></tex-math></disp-formula></p>
</sec>
<sec id="SECE.2"><title>E.2. <inline-formula><tex-math notation="LaTeX" id="ImEquation245"><![CDATA[$\psi_1+\psi_1\to\psi_1+\psi_2$]]></tex-math></inline-formula></title>
<p>The expression for <inline-formula><tex-math notation="LaTeX" id="ImEquation246"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> is slightly complicated, containing <inline-formula><tex-math notation="LaTeX" id="ImEquation247"><![CDATA[$M_2 = {m_2}^2 - {m_1}^2$]]></tex-math></inline-formula>:
<disp-formula id="ptz084-ME-9"><label>(E.9)</label><tex-math notation="LaTeX" id="Equation72"><![CDATA[
\begin{eqnarray}
|\hat{{\bf p}}_i| &=& |{\bf p}_1| \sqrt{ 1-\frac{M^2}{2|{\bf p}_1|^2}+\frac{M^4}{16|{\bf p}_1|^2(|{\bf p}_1|^2+m_1^2)}} \nonumber \\
&=& |{\bf p}_1| \left[ 1- \frac{M^2}{4|{\bf p}_1|^2} + {\cal O} \left(\left(\frac{M}{|{\bf p}_1|}\right)^4\right) \right] . 
\label{E9}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Using this <inline-formula><tex-math notation="LaTeX" id="ImEquation248"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> in the <inline-formula><tex-math notation="LaTeX" id="ImEquation249"><![CDATA[${\cal M}$]]></tex-math></inline-formula> of Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-5">E.5</xref>), we have
<disp-formula id="ptz084-ME-10"><label>(E.10)</label><tex-math notation="LaTeX" id="Equation73"><![CDATA[
\begin{multline}
{\cal M}\big(\psi_1({\bf p}_1)\psi_1(-{\bf p}_2)\to\psi_1(\hat {\bf p}_i)\psi_2(-\hat {\bf p}_i)\big) 
\\
=8\alpha\left( (6+ 2\cos^2\theta)|{\bf p}_1|^4+\left[8m_1^2-M^2(1+\cos^2\theta) \right]|{\bf p}_1|^2+ m_1^2\bigg(3m_1^2-\frac{M^2}{ 2}\bigg) \right. \\
\left. -\frac{M^4}{8}+ \frac{M^4|{\bf p}_1|^2}{8(|{\bf p}_1|^2+m_1^2)} \cos^2\theta \vphantom{\frac{M^2}{ 2}} \right) . 
\label{E10}
\end{multline}
]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) now takes the form
<disp-formula id="ptz084-ME-11"><label>(E.11)</label><tex-math notation="LaTeX" id="Equation74"><![CDATA[
\begin{multline}
-64 \alpha^2\frac{m_1^4|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi \sin\theta (6+2\cos^2 \theta)^2 
\\
\times
\left[ \left(\frac{|{\bf p}_1|}{m_1}\right)^4 + \frac{1}{6+2\cos^2\theta}\left(16-\frac{M^2}{2m_1^2}(7+5\cos^2 \theta) \right) \left( \frac{|{\bf p}_1|}{m_1}\right)^2 + {\cal O}\left( \left( \frac{m_1}{|{\bf p}_1|}\right)^0 \right) \right] d \theta. \label{1112}
\end{multline}
]]></tex-math></disp-formula></p>
</sec>
<sec id="SECE.3"><title>E.3. <inline-formula><tex-math notation="LaTeX" id="ImEquation250"><![CDATA[$\psi_1+\psi_1\to\psi_2+\psi_2$]]></tex-math></inline-formula></title>
<p>The two particles in the final state have the same mass <inline-formula><tex-math notation="LaTeX" id="ImEquation251"><![CDATA[$m_2$]]></tex-math></inline-formula> and hence <inline-formula><tex-math notation="LaTeX" id="ImEquation252"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> is easily found:
<disp-formula id="ptz084-ME-12"><label>(E.12)</label><tex-math notation="LaTeX" id="Equation75"><![CDATA[
\begin{eqnarray}
|\hat{{\bf p}}_i| &=& |{\bf p}_1| \sqrt{ 1 - \frac{M^2}{|{\bf p}_1|^2} } \nonumber \\
&=& |{\bf p}_1|\left[ \left( 1 - \frac{M^2}{2|{\bf p}_1|^2}\right) + {\cal O} \left(\left( \frac{M}{|{\bf p}_1|}\right)^4 \right)\right]\!. 
\label{E12}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Using this <inline-formula><tex-math notation="LaTeX" id="ImEquation253"><![CDATA[$|\hat{{\bf p}}_1|$]]></tex-math></inline-formula> in the <inline-formula><tex-math notation="LaTeX" id="ImEquation254"><![CDATA[${\cal M}$]]></tex-math></inline-formula> of Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-5">E.5</xref>), we have
<disp-formula id="ptz084-ME-13"><label>(E.13)</label><tex-math notation="LaTeX" id="Equation76"><![CDATA[
\[\begin{align}
& M({{\psi }_{1}}({{\mathbf{p}}_{1}}){{\psi }_{1}}(-{{\mathbf{p}}_{2}})\to {{\psi }_{2}}({{\widehat{\mathbf{p}}}_{i}}){{\psi }_{2}}(-{{\widehat{\mathbf{p}}}_{i}})) \\ 
 & \,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,=8\alpha \left( (6+2{{\cos }^{2}}\theta )|{{\mathbf{p}}_{1}}{{|}^{4}}+\left[ 8m_{1}^{2}-2{{M}^{2}}(1+{{\cos }^{2}}\theta ) \right]|{{\mathbf{p}}_{1}}{{|}^{2}}+m_{1}^{2}(3m_{1}^{2}-{{M}^{2}}) \right). \\ 
\end{align}\label{E13}\]
]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) now takes the form
<disp-formula id="ptz084-ME-14"><label>(E.14)</label><tex-math notation="LaTeX" id="Equation77"><![CDATA[
\[\begin{align}
& 64{{\alpha }^{2}}\frac{m_{1}^{4}|{{\mathbf{p}}_{1}}{{|}^{5}}}{{{E}_{1}}+{{E}_{2}}}\int_{0}^{\pi }{\sin }\theta {{(6+2{{\cos }^{2}}\theta )}^{2}} \\ 
 & \,\,\,\,\,\times \left[ {{\left( \frac{|{{\mathbf{p}}_{1}}|}{{{m}_{1}}} \right)}^{4}}+\frac{1}{6+2{{\cos }^{2}}\theta }\left( 16-\frac{{{M}^{2}}}{m_{1}^{2}}(7+5{{\cos }^{2}}\theta ) \right){{\left( \frac{|{{\mathbf{p}}_{1}}|}{{{m}_{1}}} \right)}^{2}}+O\left( {{\left( \frac{{{m}_{1}}}{|{{\mathbf{p}}_{1}}|} \right)}^{0}} \right) \right]d\theta . \\ 
\end{align}\label{1122}\]
]]></tex-math></disp-formula></p>
<p>We note that if we set <inline-formula><tex-math notation="LaTeX" id="ImEquation255"><![CDATA[${M}^2 =0$]]></tex-math></inline-formula>, all three amplitudes in Eqs. (<xref ref-type="disp-formula" rid="ptz084-ME-7">E.7</xref>), (<xref ref-type="disp-formula" rid="ptz084-ME-10">E.10</xref>), and (<xref ref-type="disp-formula" rid="ptz084-ME-13">E.13</xref>) coincide.</p>
</sec>
<sec id="SECE.4"><title>E.4. <inline-formula><tex-math notation="LaTeX" id="ImEquation256"><![CDATA[$\sigma+\psi_1\to\sigma+\psi_1$]]></tex-math></inline-formula></title>
<p>Since the two final-state particles are the same as the two initial ones, <inline-formula><tex-math notation="LaTeX" id="ImEquation257"><![CDATA[$|\hat{{\bf p}}_i| = |{\bf p}_1|$]]></tex-math></inline-formula>. Using this <inline-formula><tex-math notation="LaTeX" id="ImEquation258"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> in the <inline-formula><tex-math notation="LaTeX" id="ImEquation259"><![CDATA[${\cal M}$]]></tex-math></inline-formula> of Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-6">E.6</xref>), we have
<disp-formula id="ptz084-ME-15"><label>(E.15)</label><tex-math notation="LaTeX" id="Equation78"><![CDATA[
\begin{eqnarray}
{\cal M}\big(\sigma({\bf p}_1)\psi_1({\bf p}_2)\to\sigma(\hat {\bf p}_i)\psi_1(\hat {\bf p}_j)\big) =-4\alpha'\big( (1-\cos\theta)|{\bf p}_1|^2+m_\sigma^2\big).\label{E15}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) now takes the form
<disp-formula id="ptz084-ME-16"><label>(E.16)</label><tex-math notation="LaTeX" id="Equation79"><![CDATA[
\begin{equation}
16 {\alpha'}^2\frac{|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi \sin\theta \left[ (1-\cos\theta)^2 + (1-\cos\theta) \frac{2m_\sigma^2}{|{\bf p}_1|^2}+ {\cal O}\left(\left(\frac{m_1}{|{\bf p}_1|}\right)^4\right) \right] d \theta.\label{s1s1}
\end{equation}]]></tex-math></disp-formula></p>
</sec>
<sec id="SECE.5"><title>E.5. <inline-formula><tex-math notation="LaTeX" id="ImEquation260"><![CDATA[$\sigma+\psi_1\to\sigma+\psi_2$]]></tex-math></inline-formula></title>
<p>This process involves three unequal masses, and <inline-formula><tex-math notation="LaTeX" id="ImEquation261"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> in terms of <inline-formula><tex-math notation="LaTeX" id="ImEquation262"><![CDATA[$|{\bf p}_1|$]]></tex-math></inline-formula> is more involved; we only show the high-energy approximation,
<disp-formula id="ptz084-ME-17"><label>(E.17)</label><tex-math notation="LaTeX" id="Equation80"><![CDATA[
\begin{eqnarray}
|\hat{{\bf p}}_1|= |{\bf p}_1|\left[ 1 - \frac{M^2}{4|{\bf p}_1|^2} + {\cal O} \left(\left(\frac{m_1}{|{\bf p}_1|}\right)^4\right)\right]\!. \label{E17}
\end{eqnarray}
]]></tex-math></disp-formula></p>
<p>Using this <inline-formula><tex-math notation="LaTeX" id="ImEquation263"><![CDATA[$|\hat{{\bf p}}_i|$]]></tex-math></inline-formula> in the <inline-formula><tex-math notation="LaTeX" id="ImEquation264"><![CDATA[${\cal M}$]]></tex-math></inline-formula> of Eq. (<xref ref-type="disp-formula" rid="ptz084-ME-6">E.6</xref>), we have
<disp-formula id="ptz084-ME-18"><label>(E.18)</label><tex-math notation="LaTeX" id="Equation81"><![CDATA[
\begin{align}
{\cal M}\big(\sigma({\bf p}_1)\psi_1({\bf p}_2)\to & \, \sigma(\hat {\bf p}_i)\psi_2(\hat {\bf p}_j)\big) \nonumber \\
& = 4\alpha'\big( \sqrt{|{\bf p}_1|^2+m_\sigma^2}\sqrt{|\hat{{\bf p}}_i|^2+m_\sigma^2}-|{\bf p}_1||\hat{{\bf p}}_i| \cos \theta \big) 
\nonumber \\
& = 4\alpha' |{\bf p}_1|^2 \left[ (1-\cos\theta)\left( 1-\frac{M^2}{4|{\bf p}_1|^2}\right) \frac{m_\sigma^2}{|{\bf p}_1|^2}+ {\cal O}\left(\left(\frac{m_1}{|{\bf p}_1|}\right)^4\right) \right]\!.\label{E18}
\end{align}
]]></tex-math></disp-formula></p>
<p>Equation (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>) now takes the form
<disp-formula id="ptz084-ME-19"><label>(E.19)</label><tex-math notation="LaTeX" id="Equation82"><![CDATA[
\begin{equation}
-16 {\alpha'}^2\frac{|{\bf p}_1|^5}{E_1+E_2}\int_0^\pi \sin\theta \left[ (1-\cos\theta)^2 \left( 1-\frac{3M^2}{4|{\bf p}_1|^2}\right) + \frac{2m_\sigma^2}{|{\bf p}_1|^2}+ {\cal O}\left(\left(\frac{m_1}{|{\bf p}_1|}\right)^4\right) \right] d \theta.\label{s1s2}
\end{equation}]]></tex-math></disp-formula></p>
</sec>
</sec>
</app>
</app-group>
<fn-group>
<title>Footnotes</title>
<fn id="FN1"><p><sup>1</sup> The condition <inline-formula><tex-math notation="LaTeX" id="ImEquation265"><![CDATA[$SS^\dagger=1$]]></tex-math></inline-formula> in field theories containing negative-norm states is often called &#x201C;pseudo-unitarity&#x201D; [<xref ref-type="bibr" rid="B1">1</xref>,<xref ref-type="bibr" rid="B2">2</xref>]. While the word &#x201C;pseudo-unitarity&#x201D; might give the impression that we consider theories with negative norms, the purpose of this paper is to show that the relation to renormalizability is independent of the norm positivity. To avoid giving such a misleading impression, we use &#x201C;<inline-formula><tex-math notation="LaTeX" id="ImEquation266"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity&#x201D; instead of &#x201C;pseudo-unitarity.&#x201D;</p></fn>
<fn id="FN2"><p><sup>2</sup> Here, we mention Einstein gravity as one of the examples that possesses neither unitarity nor renormalizability.</p></fn>
<fn id="FN3"><p><sup>3</sup> The renormalization group flow is discussed in Ref. [<xref ref-type="bibr" rid="B12">12</xref>].</p></fn>
<fn id="FN4"><p><sup>4</sup> Recall that the ordinary PCR condition is not necessarily equivalent to renormalizablitiy in the presence of extra symmetries and/or constraints. The same caveat applies to the extended PCR condition in general. However, the interactions considered in this paper have no such extra symmetries or constraints, so that the extended PCR condition can be used to dictate renormalizability.</p></fn>
<fn id="FN5"><p><sup>5</sup> The difference between the renormalizable interaction in Eq. (<xref ref-type="disp-formula" rid="ptz084-M10">10</xref>) and the non-renormalizable interaction in (14) can be also explained by shift symmetry. However, some of the nonrenormalizable interaction terms with shift symmetry could not be excluded by only the nonnegativity of their coefficients. The dimensions of each factor in the interaction terms are indeed crucial to distinguish the renormalizability. Detailed discussion of this is presented in Ref. [<xref ref-type="bibr" rid="B8">8</xref>].</p></fn>
<fn id="FN6"><p><sup>6</sup> An <inline-formula><tex-math notation="LaTeX" id="ImEquation267"><![CDATA[$S$]]></tex-math></inline-formula>-matrix containing a negative-norm state is also discussed in Ref. [<xref ref-type="bibr" rid="B13">13</xref>]</p></fn>
<fn id="FN7"><p><sup>7</sup> To determine the precise value of the upper bound (which was simply denoted by &#x201C;const.&#x201D;), we need to carefully evaluate the normalization factor originating from the delta function in Eq. (<xref ref-type="disp-formula" rid="ptz084-M22">22</xref>). See, e.g., Ref. [<xref ref-type="bibr" rid="B8">8</xref>] for details.</p></fn>
<fn id="FN8"><p><sup>8</sup> Choosing the other quantization creates negative-energy particles, and no negative-norm state appears. The cancellation in Eq. (<xref ref-type="disp-formula" rid="ptz084-M36">36</xref>), as we will explain, never occurs, and thus the theory does not possess the perturbative <inline-formula><tex-math notation="LaTeX" id="ImEquation268"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitary. The relation between renormalizability and <inline-formula><tex-math notation="LaTeX" id="ImEquation269"><![CDATA[$S$]]></tex-math></inline-formula>-matrix unitarity suggests that the perturbative renormalizability property is never better in such quantization with higher-order derivatives.</p></fn>
<fn id="FN9"><p><sup>9</sup> If Lorentz symmetry is broken, the decay rate of vacuum can be mild and we would be able to live with the ghost modes [<xref ref-type="bibr" rid="B20">20</xref>,<xref ref-type="bibr" rid="B21">21</xref>].</p></fn>
</fn-group>
<ref-list id="ref1">
<title>References</title>
<ref id="B1"><label>[1]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Anselmi</surname> <given-names>D.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>94</volume>, <fpage>025028</fpage> (<year>2016</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1606.06348">arXiv:1606.06348</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1606.06348">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.94.025028">https://doi.org/10.1103/PhysRevD.94.025028</ext-link></comment>)</mixed-citation></ref>
<ref id="B2"><label>[2]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Karowski</surname> <given-names>M.</given-names></string-name></person-group>, <source>Nuovo Cim. A</source> <volume>23</volume>, <fpage>126</fpage> (<year>1974</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1007/BF02748296">https://doi.org/10.1007/BF02748296</ext-link></comment>)</mixed-citation></ref>
<ref id="B3"><label>[3]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Llewellyn Smith</surname> <given-names>C. H.</given-names></string-name></person-group>, <source>Phys. Lett. B</source> <volume>46</volume>, <fpage>233</fpage> (<year>1973</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/0370-2693(73)90692-8">https://doi.org/10.1016/0370-2693(73)90692-8</ext-link></comment>)</mixed-citation></ref>
<ref id="B4"><label>[4]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Cornwall</surname> <given-names>J. M.</given-names></string-name>, <string-name name-style="western"><surname>Levin</surname> <given-names>D. N.</given-names></string-name>, and <string-name name-style="western"><surname>Tiktopoulos</surname> <given-names>G.</given-names></string-name></person-group>, <source>Phys. Rev. Lett.</source> <volume>30</volume>, <fpage>1268</fpage> (<year>1973</year>); <volume>31</volume>, <fpage>572</fpage> (<year>1973</year>) [erratum]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevLett.30.1268">https://doi.org/10.1103/PhysRevLett.30.1268</ext-link>; <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevLett.31.572">https://doi.org/10.1103/PhysRevLett.31.572</ext-link></comment>)</mixed-citation></ref>
<ref id="B5"><label>[5]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Cornwall</surname> <given-names>J. M.</given-names></string-name>, <string-name name-style="western"><surname>Levin</surname> <given-names>D. N.</given-names></string-name>, and <string-name name-style="western"><surname>Tiktopoulos</surname> <given-names>G.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>10</volume>, <fpage>1145</fpage> (<year>1974</year>); <volume>11</volume>, <fpage>972</fpage> (<year>1975</year>) [erratum]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.10.1145">https://doi.org/10.1103/PhysRevD.10.1145</ext-link>; <ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.11.972">https://doi.org/10.1103/PhysRevD.11.972</ext-link></comment>)</mixed-citation></ref>
<ref id="B6"><label>[6]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Berends</surname> <given-names>F. A.</given-names></string-name> and <string-name name-style="western"><surname>Gastmans</surname> <given-names>R.</given-names></string-name></person-group>, <source>Nucl. Phys. B</source> <volume>88</volume>, <fpage>99</fpage> (<year>1975</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/0550-3213(75)90528-3">https://doi.org/10.1016/0550-3213(75)90528-3</ext-link></comment>)</mixed-citation></ref>
<ref id="B7"><label>[7]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Fujimori</surname> <given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Inami</surname> <given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Izumi</surname> <given-names>K.</given-names></string-name>, and <string-name name-style="western"><surname>Kitamura</surname> <given-names>T.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>91</volume>, <fpage>125007</fpage> (<year>2015</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1502.01820">arXiv:1502.01820</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1502.01820">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.91.125007">https://doi.org/10.1103/PhysRevD.91.125007</ext-link></comment>)</mixed-citation></ref>
<ref id="B8"><label>[8]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Fujimori</surname> <given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Inami</surname> <given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Izumi</surname> <given-names>K.</given-names></string-name>, and <string-name name-style="western"><surname>Kitamura</surname> <given-names>T.</given-names></string-name></person-group>, <source>Prog. Theor. Exp. Phys.</source> <volume>2016</volume>, <fpage>013B08</fpage> (<year>2016</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1510.07237">arXiv:1510.07237</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1510.07237">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1093/ptep/ptv185">https://doi.org/10.1093/ptep/ptv185</ext-link></comment>)</mixed-citation></ref>
<ref id="B9"><label>[9]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Arkani-Hamed</surname> <given-names>N.</given-names></string-name></person-group>, <source>Towards deriving string theory as the weakly coupled UV completion of gravity</source>. <comment>Talk at Strings 2016</comment>.</mixed-citation></ref>
<ref id="B10"><label>[10]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Stelle</surname> <given-names>K. S.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>16</volume>, <fpage>953</fpage> (<year>1977</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.16.953">https://doi.org/10.1103/PhysRevD.16.953</ext-link></comment>)</mixed-citation></ref>
<ref id="B11"><label>[11]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Fradkin</surname> <given-names>E. S.</given-names></string-name> and <string-name name-style="western"><surname>Tseytlin</surname> <given-names>A. A.</given-names></string-name></person-group>, <source>Nucl. Phys. B</source> <volume>201</volume>, <fpage>469</fpage> (<year>1982</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/0550-3213(82)90444-8">https://doi.org/10.1016/0550-3213(82)90444-8</ext-link></comment>)</mixed-citation></ref>
<ref id="B12"><label>[12]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Elizalde</surname> <given-names>E.</given-names></string-name>, <string-name name-style="western"><surname>Jacksenaev</surname> <given-names>A. G.</given-names></string-name>, <string-name name-style="western"><surname>Odintsov</surname> <given-names>S. D.</given-names></string-name>, and <string-name name-style="western"><surname>Shapiro</surname> <given-names>I. L.</given-names></string-name></person-group>, <source>Phys. Lett. B</source> <volume>328</volume>, <fpage>297</fpage> (<year>1994</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/0370-2693(94)91483-4">https://doi.org/10.1016/0370-2693(94)91483-4</ext-link></comment>)</mixed-citation></ref>
<ref id="B13"><label>[13]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Salvio</surname> <given-names>A.</given-names></string-name></person-group>, <source>Phys. Rev. D</source> <volume>94</volume>, <fpage>096007</fpage> (<year>2016</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1608.01194">arXiv:1608.01194</ext-link> [hep-ph]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1608.01194">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1103/PhysRevD.94.096007">https://doi.org/10.1103/PhysRevD.94.096007</ext-link></comment>)</mixed-citation></ref>
<ref id="B14"><label>[14]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Abe</surname> <given-names>Y.</given-names></string-name>, <string-name name-style="western"><surname>Inami</surname> <given-names>T.</given-names></string-name>, <string-name name-style="western"><surname>Izumi</surname> <given-names>K.</given-names></string-name>, and <string-name name-style="western"><surname>Kitamura</surname> <given-names>T.</given-names></string-name></person-group>, <source>Prog. Theor. Exp. Phys.</source> <volume>2018</volume>, <fpage>031E01</fpage> (<year>2018</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1712.06305">arXiv:1712.06305</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1712.06305">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1093/ptep/pty010">https://doi.org/10.1093/ptep/pty010</ext-link></comment>)</mixed-citation></ref>
<ref id="B15"><label>[15]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Smolin</surname> <given-names>L.</given-names></string-name></person-group>, <source>Phys. Lett. B</source> <volume>93</volume>, <fpage>95</fpage> (<year>1980</year>). (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/0370-2693(80)90103-3">https://doi.org/10.1016/0370-2693(80)90103-3</ext-link></comment>)</mixed-citation></ref>
<ref id="B16"><label>[16]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Bhattacharjee</surname> <given-names>S.</given-names></string-name> and <string-name name-style="western"><surname>Majumdar</surname> <given-names>P.</given-names></string-name></person-group>, <source>Nucl. Phys. B</source> <volume>885</volume>, <fpage>481</fpage> (<year>2014</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1210.0497">arXiv:1210.0497</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1210.0497">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1016/j.nuclphysb.2014.05.031">https://doi.org/10.1016/j.nuclphysb.2014.05.031</ext-link></comment>)</mixed-citation></ref>
<ref id="B17"><label>[17]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Abe</surname> <given-names>Y.</given-names></string-name>, <string-name name-style="western"><surname>Horikoshi</surname> <given-names>M.</given-names></string-name>, and <string-name name-style="western"><surname>Inami</surname> <given-names>T.</given-names></string-name></person-group>, <source>Commun. Phys.</source> <volume>26</volume>, <fpage>229</fpage> (<year>2016</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1602.03792">arXiv:1602.03792</ext-link> [hep-ph]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1602.03792">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.15625/0868-3166/26/3/9033">https://doi.org/10.15625/0868-3166/26/3/9033</ext-link></comment>)</mixed-citation></ref>
<ref id="B18"><label>[18]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Salvio</surname> <given-names>A.</given-names></string-name> and <string-name name-style="western"><surname>Strumia</surname> <given-names>A.</given-names></string-name></person-group>, <source>J. High Energy Phys.</source> <volume>1406</volume>, <fpage>080</fpage> (<year>2014</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1403.4226">arXiv:1403.4226</ext-link> [hep-ph]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1403.4226">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1007/JHEP06(2014)080">https://doi.org/10.1007/JHEP06(2014)080</ext-link></comment>)</mixed-citation></ref>
<ref id="B19"><label>[19]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Abe</surname> <given-names>Y.</given-names></string-name>, <string-name name-style="western"><surname>Horikoshi</surname> <given-names>M.</given-names></string-name>, and <string-name name-style="western"><surname>Inami</surname> <given-names>T.</given-names></string-name></person-group>, <source>Proc. 2nd LeCosPA Symp.: Everything about Gravity, Celebrating the Centenary of Einstein&#x2019;s General Relativity</source> <ext-link ext-link-type="doi" xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="https://doi.org/10.1142/9789813203952_0052">https://doi.org/10.1142/9789813203952_0052</ext-link>(LeCosPA2015), p. <fpage>392</fpage> [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/1603.00192">arXiv:1603.00192</ext-link> [hep-ph]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+1603.00192">Search <sc>in</sc>SPIRE</ext-link>].</mixed-citation></ref>
<ref id="B20"><label>[20]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Izumi</surname> <given-names>K.</given-names></string-name> and <string-name name-style="western"><surname>Tanaka</surname> <given-names>T.</given-names></string-name></person-group>, <source>Prog. Theor. Phys.</source> <volume>121</volume>, <fpage>427</fpage> (<year>2009</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/0810.4811">arXiv:0810.4811</ext-link> [hep-th]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+0810.4811">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1143/PTP.121.427">https://doi.org/10.1143/PTP.121.427</ext-link></comment>)</mixed-citation></ref>
<ref id="B21"><label>[21]</label><mixed-citation publication-type="journal"><person-group person-group-type="author"><string-name name-style="western"><surname>Izumi</surname> <given-names>K.</given-names></string-name> and <string-name name-style="western"><surname>Tanaka</surname> <given-names>T.</given-names></string-name></person-group>, <source>Prog. Theor. Phys.</source> <volume>121</volume>, <fpage>419</fpage> (<year>2009</year>) [<ext-link ext-link-type="uri" xlink:href="http://www.arxiv.org/abs/0709.0199">arXiv:0709.0199</ext-link> [gr-qc]] [<ext-link ext-link-type="uri" xlink:href="http://www.inspirehep.net/search?p=find+EPRINT+0709.0199">Search <sc>in</sc>SPIRE</ext-link>]. (<comment><ext-link ext-link-type="doi" xlink:href="https://doi.org/10.1143/PTP.121.419">https://doi.org/10.1143/PTP.121.419</ext-link></comment>)</mixed-citation></ref>
</ref-list>
</back>
</article>