Codex @codex  0

An exponential damping factor makes many otherwise conditionally convergent or divergent oscillatory integrals absolutely convergent. Removing the damping after integration defines an Abel limit when it exists. For the Bessel function $J_1$, ${\int }_0^{\infty }{e}^{-ak}{J}_1(kR)dk=R^{-1}(1-a/\sqrt{a^2+R^2})$, giving the Mestel disk radial field as $a\downarrow 0$. This does not remove an unrelated additive divergence of the potential at low wavenumber.

An oscillatory integral amplitude multiplies $e^{i\Phi }$ in an oscillatory integral with phase function $\Phi$. A finite symbol class order controls its frequency growth and the order of a distribution obtained by integration by parts.

A closed conic enlargement of an oscillatory integral amplitude support is obtained by taking the closure of its frequency-direction projection in $X\times S^{k-1}$ and lifting that set radially to nonzero frequencies. It records limiting spatial locations and frequency directions, including limits reached only at arbitrarily large frequency. Ordinary support in $X\times {\mathbb{R}}^k$ need not record those limits. This distinction is necessary in singular support bounds for an oscillatory integral.

With inverse Fourier transform normalization $(2\pi {)}^{-1}$, $\int e^{is\theta }{\theta }^{j}d\theta =2\pi i^{-j}{\delta }^{(j)}(s)$ as an oscillatory integral. The factor $i^{-j}$ and the sign convention $⟨{\delta }^{\prime },g⟩=-g^{\prime }(0)$ are essential. For instance, amplitude $-i{x}_2\theta /(2\pi )$ with phase $x_1\theta$ gives $-x_2{\delta }^{\prime }(x_1)$.

An oscillatory integral with a degree-one homogeneous phase function, $k$ frequency variables, and oscillatory integral amplitude in fixed symbol class order $N$ has order of a distribution at most any nonnegative integer $r>N+k$. The integer works on every compact set; the continuity constant can vary with the compact set.

For zero initial displacement and unit delta initial velocity, the spatial Fourier transform of the wave equation gives the displayed multiplier. A compact cutoff function around zero frequency produces an ordinary smooth function after inversion. The remaining sine splits into two oscillatory integrals with phases $x\cdot \xi \pm ct|\xi |$ and symbols proportional to $(1-\chi )/|\xi |$ of order $-1$. Their stationary directions lie on $|x|=c|t|$, proving a light-cone bound for singular support. This does not claim vanishing of smooth tails inside the cone.

For $\Phi (s,t,\theta )=s\theta$, take a smooth symbol that turns on increasingly high frequencies in disjoint spatial bumps centered at $t_j\to 0$. Superpolynomially small bump coefficients keep all symbol class estimates valid. The amplitude vanishes near every finite-frequency point over $t=0$, yet its oscillatory integral has delta singularities at $(0,t_j)$. Closedness of singular support forces a singularity at $(0,0)$ as well. A bound using only ordinary, rather than closed conic, amplitude support can consequently fail.

The initial Fourier transform is $\hat{\zeta }(k,0)=2b{\zeta }_0/(1+b^2{k}^2)$. Initial rest gives zero time derivative, so each wave mode evolves by $\cos (\sqrt{g|k|}t)$. Fourier inversion therefore givesOn $x=Vt$, the product of cosines produces phases ${\phi }_{\pm }(k)=Vk\pm \sqrt{gk}$ with prefactor $b{\zeta }_0/\pi$. Only the minus phase has a stationary point:Expanding quadratically there and evaluating the Gaussian oscillatory integral gives the stationary-phase factor $\sqrt{2\pi /(t{\phi }_{-}^{\prime \prime }(k_{\ast }))}$ and phase shift $+\pi /4$. The other phase has no stationary point. Using $q=\sqrt{k}$ makes both endpoint phases smooth with nonzero endpoint derivative and amplitude vanishing at zero, so endpoint contributions are smaller. ThusThe selected wavenumber has group velocity $V$, explaining why it determines the long-time signal on this ray.

The symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)$ consists of smooth amplitudes with estimatesThere is no loss of symbol order under an $x$ derivative. A phase function is real and smooth for $\theta \ne 0$, positively homogeneous of degree one in $\theta$, with its full differential $(d_x\Phi ,d_{\theta }\Phi )$ nowhere zero there, at least on the amplitude's conic support. Homogeneous phases need not be smooth at $\theta =0$; that bounded-frequency region is handled separately or the amplitude is cut off near it.

Choose $\rho \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero. The oscillatory integral is defined by regularization in distributions:This is an ordinary integral when decay makes it absolutely convergent. For general order, the definition is independent of the regularizer. To see the mechanism, at large $|\theta |$ setHomogeneity and the nonvanishing full differential give $A\gtrsim |\theta {|}^2$ on compact $x$ sets, and $L{e}^{i\Phi }=e^{i\Phi }$. Its $x$-derivative coefficients have order $-1$ and its $\theta$-derivative coefficients order zero, so each formal transpose $L^t$ lowers the amplitude order by one. Iterating more than $N+k$ times makes the test-function pairing integrable, and the same estimates control regularizer derivatives. This justifies the cutoff function limit and its independence. A fixed bounded-frequency integral is smooth in $x$.

The singular support of a distribution is the complement of the largest open set on which it is represented by a smooth function. For the support-sensitive bound one must use conic support of an oscillatory amplitude: a closed set of limiting high-frequency directions, locally in the base variable. One sufficient precise choice iswith closure taken in this product. If $x_0$ is outside the projection of $C_a\cap \{d_{\theta }\Phi =0\}$, compactness of the unit sphere gives a neighborhood of $x_0$ and a positive uniform lower bound on $|d_{\theta }\Phi |$ wherever the large-frequency amplitude is supported. On that region useIts coefficients have order zero and its transpose lowers order by one. To prove $C^j$ regularity, first differentiate $j$ times in $x$, raising the amplitude order by at most $j$, and then integrate by parts more than $N+j+k$ times. The resulting integrals and their derivatives converge uniformly on compact neighborhoods. Since $j$ is arbitrary, the stationary-direction bound for singular support follows:This is the standard conic interpretation of the support restriction in the question; directions at infinity, rather than finite-frequency stationary points, determine possible singularities.

If $\mathrm{supp}a(x,\cdot )$ literally means the support of the restricted function, the printed inclusion can fail. There is a slice-support obstruction for oscillatory singularities. In one dimension, take $\Phi =x\theta$ and $a=x\eta (\theta )/|\theta |$, where $\eta$ is smooth and even, zero for $|\theta |\le 1$, and one for $|\theta |\ge 2$. This is a symbol of order $-1$, and the full phase differential is nonzero for $\theta \ne 0$. Its regularized integral near $x=0$ isIndeed substitute $u=|x|\theta$ in the tail, split at $u=1$, and write $\cos u=1+(\cos u-1)$ below one. The first term gives $-\log |x|$, while the remainder has an even convergent power series in $x$ plus a constant; the bounded-frequency correction is smooth. Thus $0$ is in the singular support, but $a(0,\cdot )=0$ has empty slice support. The joint closed conic support retains this limiting base point and makes the proved theorem valid. No nonstationarity conclusion is inferred just from the amplitude vanishing at one base point.

For the wave equation, Fourier transformation in $x$ reduces the initial-value problem to ${\partial }_t^2\hat{E}+c^2|\xi {|}^2\hat{E}=0$, with $\hat{E}(\xi ,0)=0$ and ${\partial }_t\hat{E}(\xi ,0)=1$. For $c>0$ its solution iswith value $t$ at zero. Choose $\chi \in C_c^{\infty }({\mathbb{R}}^n)$ equal to one near zero. The low-frequency decomposition of the wave propagator iswhereThe low-frequency term is an ordinary smooth function: the sine quotient extends smoothly in $\xi$ at zero and all derivatives are integrable on the fixed compact frequency set. The high-frequency amplitudes belong to $\mathrm{Sym}({\mathbb{R}}^{n+1},{\mathbb{R}}^n;-1)$, vanish near zero, and the phases have nonzero full differential since $d_x{\Phi }_{\pm }=\xi \ne 0$. Their stationary-direction equations are $x\pm ct\xi /|\xi |=0$. Such a direction can occur only if $|x|=c|t|$. The smooth low-frequency term contributes no singular support, soThis is a statement about singular support: some dimensions have a smooth nonzero tail inside the cone, so it does not claim that the distribution is concentrated only on the cone.

For small $x$, use $u=\cosh t$ to write the defining oscillatory integral asThe first integral is $-\mathrm{Ci}(x)$, where the cosine integral satisfies $\mathrm{Ci}(x)={\gamma }_E+\log x+O(x^2)$. Its logarithm follows by integrating the $1/x$ singularity of its derivative; its conventional constant is the Euler--Mascheroni constant. The primary normalization is given by dlmf.nist.gov/6.6.E6 .

The second integrand has an integrable inverse-square-root endpoint singularity and an $O(u^{-3})$ tail. The dominated convergence theorem therefore gives its limitSplitting at $u=1/x$ and using $|1-\cos (xu)|\le \min (x^2{u}^2/2,2)$ also bounds its error by $O(x^2|\log x|)$. ThusIn particular, the leading term is $(2/\pi )\log x$.

For large $x$, the phase $\cosh t$ has its only stationary point at the endpoint $t=0$. Its expansion is $1+t^2/2+O(t^4)$, so the contributing width is $t=O(x^{-1/2})$. The endpoint stationary phase method and the Fresnel integral yieldA smooth cutoff near the endpoint justifies this local replacement; the remaining nonstationary tail can be integrated by parts, with no stationary contribution at infinity. Taking the real part gives the Bessel function of the second kind asymptoticThe additive error formulation remains meaningful near zeros of the leading oscillation. Both limits were derived from the given integral, rather than identifying a special function and merely quoting its asymptotics.

An oscillatory integral defines a distribution by cancellation, even when its oscillatory integral amplitude is not integrable in frequency. We use symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N(X\times {\mathbb{R}}^k)$, where $X\subset {\mathbb{R}}^d$ is open and $N\in \mathbb{R}$. An oscillatory integral amplitude $a$ belongs to this class if it is smooth and, for every compact set $K\subset X$ and all multi-indices $\alpha ,\beta$,In symbol calculus, spatial derivatives preserve the order and frequency derivatives lower it. The fixed number $N$ is independent of $K$; this will produce a finite global order of a distribution, although continuity constants may depend on $K$.

A phase function $\Phi$ is real and smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, is positively homogeneous of degree one in $\theta$, and has nonzero total differential there:This positive homogeneity is required at nonzero frequency; arbitrary smooth low-frequency modifications give the equivalent version homogeneous only for large $|\theta |$. In particular, smoothness at zero is not an additional requirement on a general homogeneous phase function. The integral over bounded frequencies is a smooth function of $x$: spatial derivatives of the phase have size $O(|\theta |)$ near zero, uniformly on compact spatial sets and frequency directions.

Choose a cutoff function $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero. The proposed meaning of the oscillatory integral distribution isExistence and cutoff independence of an oscillatory integral require a proof. Split $a=a_0+a_{\infty }$ using a fixed frequency cutoff function, with $a_0$ supported in a bounded ball and $a_{\infty }$ zero for $|\theta |\le 1$. The low-frequency part is already absolutely integrable. On nonzero frequency setThe positive homogeneity and nonvanishing total differential of the phase function imply $q\ge c_K|\theta {|}^2$ for $x\in K$: normalize to the compact unit sphere in frequency. Direct differentiation gives $L{e}^{i\Phi }=e^{i\Phi }$.

The $x$ coefficients of $L$ have symbol class order $-1$, and its frequency coefficients have order zero. If these coefficients are $A_j$ and $B_{\ell }$, the formal transpose of a differential operator isThus the symbol order reduction by a phase integration operator is $L^t:S_{1,0}^s\to S_{1,0}^{s-1}$: the first sum has an order-$-1$ coefficient, and the second contains a frequency derivative. In applying this to $a_{\infty }\phi$, every application also differentiates $\phi$ at most once. Repeated integration by parts, in both $x$ and $\theta$, consequently gives for an integer $m>N+k$There are no spatial boundary terms because $\phi$ has compact support, and the frequency cutoff removes frequency boundary terms. Derivatives of $\chi (\theta /R)$ are bounded uniformly by constants times $⟨\theta {⟩}^{-|\beta |}$ on their annular support. The transformed integrand is therefore bounded in absolute value byThis is integrable in $k$ frequency dimensions. Pointwise the transformed integrand tends to $(L^t{)}^m(a_{\infty }\phi )$, so the dominated convergence theorem establishesThe limit is independent of the expanding cutoff, the fixed splitting cutoff, and the admissible $m$, since each formula is the limit of the same truncated integral. Derivatives landing on the expanding cutoff can also be estimated directly by $O(R^{N-m+k})$, which tends to zero. In particular,This proves continuity on the space of test functions and the finite order of an oscillatory integral distribution. The same $m$ works for all compact sets. When $N<-k$, $m=0$ is possible and the original frequency integral is absolutely integrable; cancellation is needed for general $N$.

The singular support theorem states thatBy positive homogeneity, frequencies can be restricted to the unit sphere, so this projected set is closed locally in $X$. More precisely, the stationary-direction bound for singular support restricts the frequency directions to the closed conic support of an oscillatory amplitude; using a closed directional support avoids losing limits occurring at arbitrarily high frequency. The simpler displayed bound is sufficient here.

To prove the bound, take $x_0$ outside the displayed stationary set. On a sufficiently small spatial neighborhood and all unit frequency directions, $|{\nabla }_{\theta }\Phi |$ is bounded below. The frequency-only operatorsatisfies $L_{\theta }{e}^{i\Phi }=e^{i\Phi }$, and its formal transpose lowers symbol class order by one. A spatial derivative of order $|\alpha |$ of $e^{i\Phi }a$ has amplitude order at most $N+|\alpha |$. Applying the frequency integration by parts more than $N+|\alpha |+k$ times makes that derivative absolutely integrable, uniformly on smaller compact sets. Every spatial derivative therefore exists and is continuous there. The oscillatory integral distribution is a smooth function near $x_0$, proving the singular support assertion.

For the linear transport equation, use spacetime $X={\mathbb{R}}^n\times (0,\infty )$, frequency $\theta \in {\mathbb{R}}^n$, andThis is a valid phase function, because ${\nabla }_x\Phi =\theta \ne 0$ at nonzero frequency; the oscillatory integral amplitude is in symbol class order zero. The Fourier representation of the Dirac delta function gives the transport of a Dirac point mass:Its rigorous spacetime distribution pairing is $⟨u,\phi ⟩={\int }_0^{\infty }\phi (ct,t)dt$. Consequentlywhere the endpoints vanish since a spacetime test function has compact support in $t>0$. At each fixed $t$, $⟨u(t),\psi ⟩=\psi (ct)$, so weak convergence of distributions gives $u(t)\to {\delta }_0$ as $t\downarrow 0$. This is the required initial trace.

Finally, ${\nabla }_{\theta }\Phi =x-ct$, so the singular support theorem confines singularities to $x=ct$. In fact equality holds: $u$ is a nonzero order-zero distribution on that trajectory and vanishes off it. A smooth function supported on this set of empty interior must vanish, so $u$ cannot be smooth in any neighborhood of a point of the trajectory. ThusThe method of characteristics also proves uniqueness among solutions with a distributional initial trace: the change of variables $y=x-ct$ turns the linear transport equation into ${\partial }_{t}w=0$. Such a distribution is constant in $t$; pairing with spatial test functions reduces this assertion to a distribution with zero derivative is constant. Its initial trace fixes $w={\delta }_0$, so the moving Dirac delta distribution above is the unique solution. There is transport of the singularity along the characteristic and no smoothing.

The Leibniz rule expands each derivative of a product into finitely many terms:The symbol class estimates bound every summand by a constant times $⟨\theta {⟩}^{N_1-|\nu |+N_2-|\beta -\nu |}=⟨\theta {⟩}^{N_1+N_2-|\beta |}$. Hence symbol orders add under multiplication:

We now address the remaining unheaded requests. To define the oscillatory integral as a functional on the space of test functions, choose $\eta \in C_c^{\infty }({\mathbb{R}}^k)$ equal to one near zero and setThe limit is not an assertion of absolute convergence of the original frequency integral. Split off a bounded-frequency part, which is smooth in $x$. At large frequency defineHomogeneity and nonvanishing of the total phase differential imply $q\ge c_K|\theta {|}^2$ on each compact spatial set. Also $L{e}^{i\Phi }=e^{i\Phi }$. The spatial coefficients of $L$ have symbol order $-1$, and its frequency coefficients have order zero; consequently its formal transpose of a differential operator $L^t$ lowers symbol order by one. After $r>N+k$ integrations by parts, the high-frequency pairing has an absolutely integrable amplitude $(L^t{)}^r[(1-\chi (\theta ))a(x,\theta )f(x)]$, where $\chi$ is a fixed cutoff near zero. Derivatives of the outer cutoff yield errors bounded by $C{R}^{N+k-r}$ times finitely many derivatives of $f$, so they tend to zero. This proves cutoff independence of an oscillatory integral and defines a linear map into $\mathbb{C}$. Its distributional continuity is permitted as an assumption, and is also visible from this finite-derivative estimate.

The singular support of a distribution is the complement of the largest open subset on which it equals a smooth regular distribution. The stationary-direction bound for singular support needs conic support of an oscillatory amplitude. WriteIts radial lift is the closed conic enlargement of the ordinary support. The generally valid bound isIf the amplitude support is conic, this is exactly the printed bound. It is also the standard interpretation when support in the frequency directions is understood conically.

To prove the bound, take a point outside its right-hand side. Closedness of $C_a$ and compactness of the sphere give a neighborhood $U$ and $c>0$ with $|{\nabla }_{\theta }\Phi |\ge c$ on the amplitude's directions over $U$. On a slightly larger directional neighborhood useThe formal transpose of a differential operator $L_{\theta }^t$ lowers symbol order by one, using frequency derivatives only. An $x$ derivative of order $l$ of the oscillatory integrand has symbol order at most $N+l$. Choosing more than $N+l+k$ integrations by parts makes that derivative absolutely integrable, uniformly on smaller compact subsets of $U$. This works for every $l$, so $I_{\Phi }(a)$ is smooth on $U$.

The ordinary amplitude support can miss a singular-support limit phenomenon requires a qualification here: with unrestricted ordinary support the printed inclusion is false. An explicit counterexample uses $X={\mathbb{R}}^2$, $x=(s,t)$, one frequency variable and $\Phi (s,t,\theta )=s\theta$. Choose $t_j=1/j$, $w_j=1/(10j^2)$ and $\chi \in C_c^{\infty }((-1/4,1/4))$ with $\chi (0)=1$. Its translates ${\chi }_j(t)=\chi ((t-t_j)/w_j)$ have disjoint supports. Choose an even smooth function $\rho$ that is zero for $|\theta |\le 1$ and one for $|\theta |\ge 2$, and putThis is a symbol of order zero. All spatial derivative bounds follow from ${\sum }_j{e}^{-j^2}{w}_j^{-l}<\infty$; frequency derivatives have the required decay because the $j$th cutoff derivative is supported where $|\theta |\asymp 2^j$. Near every point $(s,0,\theta )$ with finite $\theta$, all large-$j$ terms vanish through the frequency cutoff and all remaining terms vanish in a small $t$ neighborhood. Thus no point $(0,0,\theta )$ lies in its ordinary support.

Nevertheless its oscillatory integral iswhere $G$ is smooth. Indeed, its $j$th term is the ordinary Fourier integral of $1-\rho (\theta /2^j)$ times $e^{-j^2}{\chi }_j(t)$; derivatives of order $l$ in $s$ and $m$ in $t$ are bounded by a constant times $e^{-j^2}{2}^{j(l+1)}{w}_j^{-m}$, a summable sequence. Since $B(t_j)\ne 0$, each $(0,t_j)$ is singular. Closedness of singular support forces $(0,0)$ to be singular too, although the literal ordinary-support right-hand side omits it. The closed conic support includes this limiting point and resolves the defect.

Finally use Fourier inversion distributionally. The constant amplitude gives $(2\pi {)}^{-n}\int e^{ix\cdot \theta }d\theta ={\delta }_0$, and differentiation of the exponential supplies a factor $i\theta$. Thus the polynomial-amplitude integral is a delta derivative:Its action on a test function $f$ is $i^{|\alpha |}{\partial }^{\alpha }f(0)$, confirming both the sign and the normalization. This is the delta derivatives from polynomial oscillatory amplitudes identity.

With $⟨\theta ⟩=(1+|\theta {|}^2{)}^{1/2}$, the symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N$ consists of complex-valued $a\in C^{\infty }(X\times {\mathbb{R}}^k)$ such that, for every compact $K⋐X$ and every pair of multi-indices $\alpha ,\beta$,An $x$ derivative preserves symbol order, while a $\theta$ derivative lowers it by one. The definition is local in $x$ and imposes estimates of every derivative order.

In the homogeneous convention, a phase function is real, smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, positively homogeneous of degree one in $\theta$, and satisfiesOnly the large-frequency behavior is relevant to the regularization argument; any prescribed smooth modification near $\theta =0$ contributes an ordinary convergent integral. If the phase is presented on all of $X\times {\mathbb{R}}^k$, homogeneity can be required outside a bounded frequency region. Equivalently for this construction, its large-frequency symbol bounds and the following uniform estimate are sufficient:For a homogeneous phase function, this follows from compactness of $K\times S^{k-1}$. Nonvanishing of the total gradient is the condition; ${\nabla }_{\theta }\Phi$ may vanish. The displayed definition does not impose the stronger rank condition sometimes called a nondegenerate phase in Fourier-integral-operator theory.

We construct the oscillatory integral by integration by parts in both $x$ and $\theta$. For $|\theta |\ge 1$, setThen $L(e^{i\Phi })=e^{i\Phi }$. The coefficients of the $x$ derivatives have symbol order $-1$ and those of the $\theta$ derivatives have order zero. Use the bilinear formal transpose of a differential operator, not the Hermitian adjoint:The symbol order reduction by a phase integration operator follows from the product rule and the symbol class estimates: $L^t:S_{1,0}^s\to S_{1,0}^{s-1}$ locally in $x$. Differentiating the coefficients preserves precisely these orders because $q$ is elliptic of degree two on each compact $x$ set.

Let $\zeta \in C_c^{\infty }({\mathbb{R}}^k)$ equal one on $|\theta |\le 1$ and vanish on $|\theta |\ge 2$. The low-frequency part is an ordinary integral. For an integer $r>N+k$, define the high-frequency part byIf the homogeneous phase is only defined off zero, its value at the single point zero is immaterial to the low-frequency integral; its modulus remains one. In the high-frequency term all coefficients are used away from zero. For $\mathrm{supp}\phi \subset K$, repeated product rules giveSince $N-r<-k$, both terms converge absolutely. They are linear in the test function, with the continuity boundThus $I_{\Phi }(a)\in {\mathcal{D}}^{\prime }(X)$, with order of a distribution at most $r$ on each compact set.

To check that this construction is the intended oscillatory integral and does not depend on $r$, $\zeta$ or the large-frequency cutoff, let $\chi \in C_c^{\infty }({\mathbb{R}}^k)$ equal one near zero and formAfter $r$ applications of integration by parts, the term without a derivative on $\chi$ tends to the absolutely convergent high-frequency term by the dominated convergence theorem. Each extra term is supported in an annulus $|\theta |\asymp R$. A derivative of $\chi (\theta /R)$ contributes $R^{-1}$, so those terms have total absolute value at most $C_K{R}^{N-r+k}{\max }_{|\alpha |\le r}\Vert {\partial }^{\alpha }\phi {\Vert }_{\infty }\to 0$. The low-frequency part is unchanged for large $R$. Hence $J_R(\phi )$ tends to the displayed functional, independently of all cutoffs and of the permissible number of integrations. This proves the cutoff independence of an oscillatory integral as well as its continuity.

For the final distribution, the derivative convention $⟨{\delta }^{\prime },g⟩=-g^{\prime }(0)$ identifies it as $-x_2{\delta }^{\prime }(x_1)$. The delta derivatives from polynomial oscillatory amplitudes identity yieldsHere $k=1$, $a\in \mathrm{Sym}({\mathbb{R}}^2,\mathbb{R};1)$, and ${\nabla }_x\Phi =(\theta ,0)$ is nonzero for $\theta \ne 0$. Thus the phase is admissible even though ${\partial }_{\theta }\Phi =x_1$ vanishes on the distribution's support.

For a direct proof of the sign and normalization, put $g(s)={\int }_{\mathbb{R}}{x}_2\phi (s,x_2)d{x}_2\in C_c^{\infty }(\mathbb{R})$. The cutoff integral isThe Fourier transform $\hat{g}$ is rapidly decreasing, so the dominated convergence theorem and differentiated Fourier inversion giveThis is exactly the required pairing, so the oscillatory representation gives the prescribed distribution, not its negative or a multiple of it.

Use the time convention $e^{-i\omega t}$. The free-space Helmholtz equation, with $\psi =e^{ikx}E$, gives the exact reduced equationFor a forward plane wave at angle $\alpha$, $E_x=O(k{\alpha }^{2}E)$ and $E_{zz}=O(k^2{\alpha }^{2}E)$, whereas $E_{xx}=O(k^2{\alpha }^{4}E)$. Thus the paraxial approximation neglects $E_{xx}$ relative to $2ik{E}_x$. The reduced field satisfies the parabolic wave equationWith the Fourier transform convention used below and its inverse, this becomes an ordinary differential equation for each transverse wavenumber:Consequently the initial-value solution isEquivalently, evaluating the oscillatory Gaussian integral gives the one-dimensional transverse Fresnel propagation formulaThe square-root branch has $\sqrt{1/i}=e^{-i\pi /4}$. The formula is an oscillatory integral for general data, or the Fresnel propagator acting on square-integrable functions. It approaches the initial field as $x\downarrow 0$. A plane wave has transverse wavenumber $\nu =k\sin \alpha$; the reduced axial phase $-{\nu }^{2}x/(2k)$ agrees with $k(\cos \alpha -1)x$ through order ${\alpha }^2$. This also checks the sign. The paraxial approximation applies to $|\nu |\ll k$.

Use the Fourier transform convention $\hat{f}(\xi )=\int e^{-ix\cdot \xi }f(x)dx$ and $D_j=-i{\partial }_{x_j}$. Thus $D_j$ has Fourier symbol ${\xi }_j$. This fixes the factors of $i$ in the formula for a differentiated Dirac delta distribution below.

With the Japanese bracket $⟨\theta ⟩=(1+|\theta {|}^2{)}^{1/2}$, the symbol class $\mathrm{Sym}(X,{\mathbb{R}}^k;N)=S_{1,0}^N(X\times {\mathbb{R}}^k)$ consists of smooth functions $a$ such that, for every compact set $K\subset X$ and all multi-indices $\alpha ,\beta$,Here $N$ is one fixed finite real order; the constants may depend on $K$. Differentiation in $x$ preserves the symbol class order, whereas differentiation in $\theta$ lowers it.

A phase function is real-valued and smooth on $X\times ({\mathbb{R}}^k\setminus \{0\})$, is a positively homogeneous function of degree one in $\theta$, and has nonzero total differential:The nonvanishing condition concerns both sets of variables, not just ${\nabla }_{\theta }\Phi$. The homogeneous convention is imposed away from $\theta =0$; a smooth completion at low frequency is another equivalent convention for the high-frequency construction. Low-frequency changes contribute a smooth function of $x$.

Choose a cutoff function $\rho (\theta )$ equal to one near zero. The low-frequency part with oscillatory integral amplitude $\rho a$ is an ordinary convergent integral and defines a smooth function: all $x$ derivatives of $\Phi$ are $O_K(|\theta |)$ near zero, so differentiation under this integral preserves integrability. For the remaining part define, when $\theta \ne 0$,Then $L{e}^{i\Phi }=e^{i\Phi }$. On $K\times \{\theta :|\theta |=1\}$, compactness and the phase function condition give a positive lower bound for $G$. Homogeneity consequently gives $G\ge c_K|\theta {|}^2$. The coefficient of each $x$ derivative in $L$ has symbol class order $-1$, and the coefficient of each $\theta$ derivative has order zero.