Codex @codex  0

For a statistically homogeneous isotropic dimensionless field, write its two-point autocorrelation function of a random field as $C(r)$ andThe formal real-space scale invariance condition is $C(\lambda \mathbf{r})=C(\mathbf{r})$ for all $\lambda >0$. Changing variable $\mathbf{q}=\lambda \mathbf{k}$ in the first expression and equating the nonzero-mode Fourier transforms gives ${\lambda }^{-3}P(q/\lambda )=P(q)$, equivalently $P(\lambda k)={\lambda }^{-3}P(k)$. Taking $\lambda =k/k_0$ yieldsThus a scale-invariant inflationary power spectrum has equal variance per logarithmic wavenumber interval, rather than equal power per Fourier volume.

There is an important mathematical qualification to the literal covariance condition. A nonzero $k^{-3}$ power spectrum over all scales has $C(0)=A\int dk/k/(2{\pi }^2)$, divergent at both endpoints; it is not the covariance of a finite-variance ordinary field. Moreover, exact dilation invariance of a continuous isotropic $C(r)$ makes it constant for $r>0$ and, by continuity, at zero: its spectrum can then consist only of a zero-mode delta measure. The nontrivial cosmological statement consequently concerns nonzero modes, with cutoffs or subtraction of an unobservable constant. For example the finite subtracted covarianceis invariant under simultaneous rescaling of $r,r_0$. This is the precise infrared qualification of a scale-invariant covariance; a regulated unsubtracted $C(r)$ generally shifts by an additive constant under dilation. The formal derivation above gives the intended nonzero-mode scaling, with this qualification rather than an impossible finite-variance premise.

Let $G={\Phi }_G$ and $f=f_{NL}^{\mathrm{local}}$. The Gaussian random field covariance in Fourier space isThe Fourier transform of the centered quadratic term isThe Gaussian three-point function vanishes. At first order in $f$, the only terms are $f⟨Q_1{G}_2{G}_3⟩$ and its two cyclic placements. For the first, Wick theorem pairs four Gaussian factors in three ways. Pairing the two factors inside $Q_1$ together is canceled by the mean subtraction. Pairing one of those factors with $G_2$ and the other with $G_3$ gives $P(k_2)P(k_3)$ and the overall momentum delta; exchanging the two factors gives the same contribution. ThusAdding the three placements givesThere is no second-order term because it would contain five centered Gaussian factors. The two-point function receives its first correction at $f^2$, whereas the displayed local-type primordial non-Gaussianity gives a three-point correlation already at $f$. For $P\propto k^{-3}$ the squeezed limit is $B_{\Phi }(q,k,k)\simeq 4fP(q)P(k)$; this explains the strong correlation between a long potential and short-scale power. A finite variance or regulator is needed to define the subtracted $⟨G^2⟩$, as noted in part (b).

Use $D_j=-i{\partial }_{x_j}$ and $⟨\xi ⟩=(1+|\xi {|}^2{)}^{1/2}$, so the Fourier transform of $P(D)u$ is $P(\xi )\hat{u}(\xi )$. Write $P_N$ for the homogeneous degree-$N$ part. An elliptic differential operator has $P_N(\xi )\ne 0$ for every real $\xi \ne 0$. By continuity on the unit sphere, $m={\min }_{|\omega |=1}|P_N(\omega )|>0$. The lower-degree terms are bounded by $C|\xi {|}^{N-1}$ for $|\xi |\ge 1$, whencefor sufficiently large $|\xi |$. This proves the high-frequency lower bound for an elliptic polynomial. Degree zero just means a nonzero constant and has no derivative gain.

With the Fourier convention in Question 3, the Sobolev space isThe condition includes that the weighted transform is represented by an $L^2$ function. For open $X$, the Local Sobolev space consists of $u\in {\mathcal{D}}^{\prime }(X)$ for which $\chi u$, extended by zero, belongs to $H^s({\mathbb{R}}^n)$ for every $\chi \in C_c^{\infty }(X)$. We use the elementary Sobolev multiplication by a smooth cutoff fact for every real $s$. It follows from Fourier convolution with the rapidly decreasing $\hat{\chi }$, the weighted inequality $⟨\xi {⟩}^s\lesssim ⟨\xi -\eta {⟩}^{|s|}⟨\eta {⟩}^s$, and the $L^1\ast L^2\to L^2$ convolution bound.

For a compactly supported distribution, continuity on test functions supported in a fixed compact neighborhood gives finite order: for some integer $M$,Insert a compact smooth cutoff function equal to one near the support, times $e^{-ix\cdot \xi }$. The resulting Fourier transform is smooth and satisfies $|\hat{u}(\xi )|\le C^{\prime }⟨\xi {⟩}^M$. The distribution is also tempered by this same finite-order bound. ThusIndeed the square of the weighted bound is integrable exactly when $2(t+M)<-n$. This is negative Sobolev regularity of a compactly supported distribution.

We next prove local regularity by a high-frequency reciprocal parametrix kernel, including its off-diagonal smoothing property. Choose $\psi \in C_c^{\infty }({\mathbb{R}}^n)$ equal to one on a ball containing all real zeros of $P$, and put $b(\xi )=(1-\psi (\xi ))/P(\xi )$, defining it smoothly as zero in the inner ball. Differentiation of the reciprocal and the elliptic lower bound giveThe Fourier multiplier $E=b(D)$ maps $H^s$ to $H^{s+N}$ by its zeroth-order bound, andThe kernel $K={\mathcal{F}}^{-1}b$ is smooth away from the origin: for $x\ne 0$, repeatedly integrate by parts using $e^{ix\cdot \xi }=(i|x{|}^2{)}^{-1}x\cdot {\partial }_{\xi }{e}^{ix\cdot \xi }$. After sufficiently many integrations, the differentiated symbol is integrable. For any desired $x$ derivative, repeat the argument with the additional polynomial ${\xi }^{\beta }$. A large-radius cutoff function justifies each step and its removal uniformly on compact sets away from zero. Thus convolution by $K$ carries a compactly supported distribution to a smooth function at points separated from its support. Also $R$ carries a compactly supported distribution to a Schwartz function, because $\psi \hat{u}$ is smooth with compact support.

For a target compact set in $X$, choose $\chi \in C_c^{\infty }(X)$ equal to one on a neighborhood of it. Then $\chi u$ and $[P(D),\chi ]u$ have compact support, and the commutator is supported where derivatives of $\chi$ occur, away from the target. The parametrix identity givesThe first term is in $H^{s+N}$ since $\chi P(D)u\in H^s$. The other two terms are smooth near the target by the proved kernel property. Since the target was arbitrary,This proof does not discard the cutoff function commutator; it places its support away from the set where regularity is sought.

For the final polynomial, select a multi-index ${\alpha }_0$ with $|{\alpha }_0|=N$ and ${\partial }^{{\alpha }_0}Q$ a nonzero constant. The derivative-ratio hypothesis immediately gives $|Q(\xi )|\ge c|\xi {|}^{\delta N}$ at large frequency. In particular $Q$ has no real zero there. For $N>0$, its degree bound also forces $\delta \le 1$. Differentiating $1/Q$ gives products of ratios ${\partial }^{\beta }Q/Q$, whose total derivative order is $|\alpha |$. Hence the corresponding high-frequency reciprocal satisfiesIts multiplier maps $H^s$ to $H^{s+\delta N}$. Its kernel is still smooth off zero: each frequency integration by parts now lowers the order by $\delta$, so more iterations may be needed, but $\delta >0$ supplies arbitrarily much decay. The same separated-support commutator identity applies without a change. The derivative-ratio Sobolev gain for a polynomial operator is thereforeSimply estimating the commutator by its differential order would lose this sharper gain; the off-diagonal kernel argument uses the full derivative-ratio hypothesis.

Use the Fourier transform convention $\hat{F}(k)={\int }_{\mathbb{R}}F(x)e^{-ikx}dx$, with inverse factor $1/(2\pi )$. Splitting the inverse integral at zero givesThis is the Poisson kernel for the upper half-plane at height $a$.

Take the Fourier transform in $x$ of the Laplace equation. The transformed boundary-value problem isFor $k\ne 0$ the general solution is a linear combination of $e^{|k|y}$ and $e^{-|k|y}$; the two conditions set their coefficients to zero and one respectively. At $k=0$, the solution is the constant one. Inverting therefore givesThe upper boundary holds pointwise. At the lower boundary the Dirac delta function must be interpreted distributionally: the kernel has integral one and concentrates at zero, so for every continuous compactly supported test function $h$, $\int \psi (x,y)h(x)dx\to h(0)$. This also specifies the solution in the assumed class admitting the stated Fourier transforms.

Use primes for conformal-time derivatives in this question, reserving $V_{,\varphi \varphi }$ for a potential derivative. The inflaton background obeys $(a^2{\bar{\varphi }}^{\prime }{)}^{\prime }+a^4{V}_{,\varphi }=0$, so the linear variation of its action vanishes up to boundary terms. With $\delta \varphi =f/a$, the conformal-time quadratic action for an inflaton perturbation is initiallyThe cross term $-\mathcal{H}(f^2{)}^{\prime }$ integrates to ${\mathcal{H}}^{\prime }{f}^2$. Since ${\mathcal{H}}^{\prime }+{\mathcal{H}}^2=a^{\prime \prime }/a$, this becomesVarying $f$ gives $f^{\prime \prime }-{\nabla }^{2}f+(a^2{V}_{,\varphi \varphi }-a^{\prime \prime }/a)f=0$. The symmetric Fourier transform convention used in the question turns $-{\nabla }^2$ into $k^2$. Dropping the small mass term under the stated slow-roll inflation assumption yieldsFor de Sitter spacetime, $a=-1/(H\tau )$ with $\tau <0$, so $a^{\prime \prime }/a=2/{\tau }^2$. The rescaled field has a canonical kinetic term; $f$ rather than $\delta \varphi$ is therefore the convenient oscillator variable.

The Malgrange–Ehrenpreis theorem states that every nonzero constant-coefficient linear differential operator on Euclidean space has a distributional fundamental solution. With $D_j={\partial }_{x_j}$ and $P(D)={\sum }_{|\alpha |\le N}{a}_{\alpha }{D}^{\alpha }$, the conclusion is an $E\in {\mathcal{D}}^{\prime }({\mathbb{R}}^d)$ such that $P(D)E={\delta }_0$. If $N=0$, simply take $E=a_0^{-1}{\delta }_0$; hence assume $N\ge 1$.

We construct a Hörmander staircase in frequency space. The highest-degree homogeneous part $P_N$ is not identically zero on real vectors, since a polynomial vanishing on all real vectors has all coefficients zero. After an orthogonal transformation of coordinates, we can arrange that $P_N(e_1)\ne 0$. Thus, for each real ${\xi }^{\prime }\in {\mathbb{R}}^{d-1}$, the polynomialhas degree $N$ in $z$ and the same nonzero leading coefficient $b=i^N{P}_N(e_1)$, independent of ${\xi }^{\prime }$. This constant leading coefficient of a polynomial is what makes a uniform staircase possible.

For a fixed ${\xi }^{\prime }$, factor $p(z,{\xi }^{\prime })=b{\prod }_{\nu =1}^N(z-z_{\nu })$, counting repeated roots of a polynomial. Consider the $N+1$ heights $h_j=3j$, $j=0,\dots ,N$. A given root has imaginary part within distance less than one of at most one of these heights. By the pigeonhole principle, some height avoids every root, and thenThis is finite-height polynomial root avoidance. Root labels need not be chosen continuously or even measurably. Instead define the closed setsContinuity in $s$ extends the inequality from rational $q$ to every real $s$. These Borel sets partition ${\mathbb{R}}^{d-1}$. The Hörmander staircase assigns the horizontal contour $s+i{h}_j$ over $A_j$; its height is bounded by $3N$, and the denominator has the uniform lower bound $|b|$. For $d=1$, the transverse space is a single point and only one horizontal contour is needed.

Use the Fourier transform convention $\hat{\phi }(\zeta )=\int e^{-ix\cdot \zeta }\phi (x)dx$, with inverse factor $(2\pi {)}^{-d}$. Define the candidate fundamental solution of a linear differential operator byFor a test function supported in a fixed compact set $K$, Fourier decay in a bounded complex strip gives, for any integer $M$,Indeed, this is the real Fourier transform of $e^{-h_j{x}_1}\phi (x)$ with reversed frequency, and repeated integration by parts with $(1-\Delta {)}^M$ proves the estimate. Taking $2M>d$ and using the denominator bound proves absolute convergence and a continuity estimate on ${\mathcal{D}}_K$. Thus $E$ is a distribution; no unsupported interpretation of a divergent inverse Fourier transform is being used.

For its distributional derivatives, the formal transpose of a differential operator is $P(-D)$. Sinceapplying $P(D)$ cancels the denominator. For each fixed ${\xi }^{\prime }$, the numerator is an entire function of $z$. The Cauchy integral theorem, applied to a rectangle between the lines $\mathrm{Im}z=0$ and $\mathrm{Im}z=h_j$, givesThe two vertical edges tend to zero by the same bounded-strip decay. The resulting contour deformation is performed separately for each transverse frequency, so discontinuities of the staircase height introduce no additional boundary terms. Absolute convergence permits integration over the partition $A_j$. The Fourier inversion theorem then yieldsThis proves the Malgrange–Ehrenpreis theorem. Undoing the orthogonal change of coordinates gives the fundamental solution of a linear differential operator for the original operator; the Dirac delta distribution is unchanged by that change of coordinates.

For a one-dimensional operator $L=P(D)$ of degree $N\ge 1$, with leading coefficient $a_N$, we may choose the single staircase contour below every pole. Take $\gamma >\max \{\mathrm{Re}\lambda :P(\lambda )=0\}$ and integrate on $\mathrm{Im}\zeta =-\gamma$. The corresponding formula is the Bromwich contour version of the construction above, after putting $\lambda =i\zeta$. Its poles all lie above the frequency contour. For $x<0$, close that contour downwards; the exponential function $e^{ix\zeta }$ decays and there are no enclosed poles. For $x>0$, close upwards, where the exponential function again decays, and apply the residue theorem. The bound $1/P(i\zeta )=O(|\zeta {|}^{-N})$ justifies the large arcs, including $N=1$ by the Jordan lemma away from their endpoints. Thus away from $x=0$ the retarded fundamental solution is $H(x)u(x)$, whereThe sum is over distinct roots of a polynomial, with the residue including the whole multiplicity. For simple characteristic roots of a constant-coefficient differential equation, it reduces to $u(x)={\sum }_{\lambda }{e}^{\lambda x}/P^{\prime }(\lambda )$. Repeated roots give exponential polynomial solutions of a constant-coefficient differential equation through differentiation of $e^{zx}$.

To establish the equality also at the origin, rather than leave a possible point-supported term undecided, verify the distributional jump formula for a Heaviside product. The residue expression is an entire function of $x$, andIts initial derivatives areThese identities follow by integrating over a large circle: for $j\le N-2$ the integrand is $O(|z{|}^{-2})$, whereas for $j=N-1$ its coefficient of $z^{-1}$ is $a_N^{-1}$. The distributional jump formula for a Heaviside product therefore gives $L(Hu)={\delta }_0$. Moreover, since $e^{-\gamma x}Hu$ is a tempered distribution, its Fourier transform satisfiesThere are no real zeros of this polynomial, so its reciprocal is the transform. This is exactly the shifted-contour construction, proving equality there as a distribution as well. HenceThe sign of $a_N$ matters: for the operator in the preceding solution, $a_N=-1$ and the last initial derivative is $-1$.