Assume initially that are integrable, so Fubini's theorem justifies interchanging the integrals. Transform the convolution and put :
This is the convolution theorem. No factor appears in this direction because the forward transform in the question is unnormalized.
Let be the indicator of . The finite integral is exactly the convolution . Its additional transform factor is
Applying the convolution theorem and part (d),
Equivalently this is . The value at zero also follows from .
By the definition of convolution and Fourier inversion,
The second line uses . The factor comes from the inverse-transform normalization. This is the product form of the convolution theorem, initially justified for sufficiently decaying functions and extended when the indicated products and convolutions exist.
For integrable functions , Fubini's theorem justifies exchanging the two integrals in the Fourier transform of the convolution. Set to obtain
This proves the convolution theorem with the stated normalization:
Direct evaluation gives the first Fourier transform:
For the even triangular function, the odd sine part integrates to zero. Two elementary integrations by parts or one evaluation of the cosine integral give
The values at zero are the respective areas, agreeing with the removable limits. Therefore . By the convolution theorem and uniqueness of the Fourier transform,
There is also a direct interpretation: the convolution is the length of overlap of the intervals and , giving exactly that triangular profile.
The Fourier transform of a derivative gives . For ,
This selects the solution decaying at both time infinities. Indeed, the inverse transform of is the one-dimensional modified Helmholtz Green function , and the convolution theorem gives . The printed equation has homogeneous solutions , so these two decay conditions matter; no causal initial-value condition was specified.
For , evaluate Fourier inversion by closing the contour in the lower half-plane, where decays. The contour is clockwise. The enclosed simple poles and residues of are
Thus inversion equals times the sum of the residues of , giving
As a check, applying gives . The convolution also gives for ; the two pieces have matching values and matching first derivatives at zero.
Use the Fourier transform convention
The convolution is . For integrable , Fubini's theorem and the substitution give
Hence the convolution theorem is
The convolution theorem gives the causal response
Define the integrated impulse response, for , by
Direct differentiation gives and . Thus
In particular, after switch-off the answer is
The response and its first derivative are continuous at both switching times, and the post-pulse motion is a freely decaying damped harmonic oscillator. The expression described as valid for in the paper is understood for ; before the pulse the causal output is zero.
Assume first that and are integrable functions, so that Fubini's theorem justifies changing the order of the two integrals. In the defining Fourier transform of the convolution, put :
This proves the convolution theorem with the stated normalization:
The convolution is the length of the overlap of with . This overlap has length when and zero otherwise, so . The convolution theorem now gives
The value at zero is also the area of the triangular function, providing a normalization check on its Fourier transform.
The Translation property of the Fourier transform gives
For the remaining integral, the same Translation property of the Fourier transform identifies as the Fourier transform of . By the product form of the convolution theorem, with inverse normalization ,
But is supported in , whereas the four translated triangular functions making up are all supported in . Their product is identically zero. Therefore the requested integral is for every real .
Extend to be zero for . Define for and
Then the full-line equation is . Use the Fourier transform convention
Under the corresponding growth bounds, the positive-support transform is analytic above a horizontal line and the negative-support transform below a horizontal line. If both functions have suitable subexponential growth, the basic domains are respectively the upper and lower half-planes; exponential decay can enlarge these to overlapping strips. The convolution theorem gives .
The Wiener-Hopf method factors into functions analytic and nonzero in the respective domains. After multiplicative or additive separation, the two sides continue to one entire function, whose growth determines a polynomial by Liouville's theorem. Real-axis zeros must be assigned consistently with the inversion contour, rather than divided out without specifying a side.
For , direct integration gives , so
Rearrange the Wiener-Hopf factorization as
The left side has the positive-support analytic continuation and the right side the negative-support one. Their growth, from the supplied bounds, makes the common continued function a constant . Consequently
For this kernel, on the negative axis is a multiple of , so the latter transform is analytic for . The common strip can therefore be chosen as . In particular the inversion line must pass above both real poles of .
Set . Then
This is the transform, on , of
Equivalently, Fourier inversion closed below the inversion line for picks up the two real poles with clockwise orientation.
As a direct check, let . Since , the equation gives . At zero, differentiation of the integral gives , so ; this sets the sine and cosine coefficients equal. Direct substitution, or the transforms above, then verifies the full integral equation, not merely its differentiated form.
For uniform illumination, the aperture is for and zero otherwise. With ,
Zeros occur at , ; the central maximum extends between . A monochromatic laser expanded into a collimated beam can illuminate the slit. A lens with the screen in its back focal plane produces the Fraunhofer pattern; alternatively use a sufficiently distant screen, . Align the slit and beam and avoid detector saturation.
For equally spaced slits, write the aperture as a convolution
The convolution theorem converts this into a single-slit transform times the finite-array transform. Therefore
Principal maxima are at for normal incidence, inside the single-slit envelope. Their height scales as times the envelope, and their first zeros are separated from a peak by .
The envelope-to-order-spacing ratio is : increasing it admits more orders within the central envelope and changes the relative order strengths. In particular, an order is missing when is a nonzero integer. Statements about absolute angular spacing need a fixed dimension: at fixed , increasing reduces order spacing; at fixed , decreasing broadens the envelope.
For fixed incident angle, differentiating the grating equation gives an order displacement . Equating this to the peak-to-first-zero width gives
Here is the illuminated grating width; this is the ideal Rayleigh limit, before slit width or aberrations broaden the spectral image.
With , the required count is . Thus 932 illuminated slits, or , suffice ideally. Checking diffraction-order feasibility, normal incidence cannot produce this second order because ; an oblique spectrograph geometry is required. A feasible reflection-Littrow example has , or .
Figure 1.
Original slit, finite-array convolution, diffraction pattern and Fourier-plane experimental layout
. Details.
Use Fourier transform convention . The convolution theorem gives
A finite array of very narrow slits at , , has transmission proportional to . It can be constructed by multiplying an infinite Dirac comb by a finite window. The infinite comb transforms to a reciprocal comb; the finite window convolves each reciprocal peak with a narrow sinc-shaped amplitude. Thus increasing the illuminated length narrows each Fraunhofer diffraction peak.
The exact finite sum makes the argument quantitative:
Its normalized intensity is . Principal maxima occur at ; their first zeros are displaced by . Replacing the ideal narrow slit by a finite slit function convolves the array in real space, multiplying this field by the single-slit Fourier transform. Field amplitudes are transformed and combined before their squared modulus is taken.
Figure 1.
A finite narrow-slit array and its Fourier intensity, showing how the principal peaks narrow with slit count
. Details.
The convolution of two aperture functions in one dimension is
With Fourier transform convention , the convolution theorem states
The first identity follows by substituting the convolution integral, putting , and separating the resulting double integral into times . For an absolutely integrable pair the interchange is justified by absolute convergence; finite optical apertures satisfy suitable integrability assumptions.
In Fraunhofer diffraction, the observed field is the Fourier transform of the transmitted complex aperture field. For two identical slits, let be one slit and use a pair of Dirac delta distributions to locate the centres:
Consequently . For a rectangular slit, , giving a single-slit diffraction envelope times the double-slit interference factor in intensity. A phase reversal replaces the plus between the deltas by a minus, giving instead, as in B7.
More generally, an -slit array is a single-aperture function convolved with a finite comb. Its field contains the factor
This separates the narrow diffraction grating maxima from the individual-aperture envelope. The same reasoning holds in two transverse dimensions. Multiplying an existing aperture by an additional finite mask instead convolves their field transforms and can broaden features. These are operations on complex fields: intensity is the squared modulus of the final field and cannot generally be substituted into the field-level convolution theorem.
The far-field complex amplitude is the two-dimensional Fourier transform of the aperture amplitude transmittance, evaluated at transverse spatial frequencies , . The light intensity is proportional to the squared magnitude of this transform.
Use linearity to add the transforms of simpler aperture components, preserving their relative phases. A translated component acquires a complex exponential factor. Repeated identical openings can be expressed as the convolution of one opening with a sum of Dirac delta functions; the convolution theorem then factors the field into a single-opening transform and an array factor. Multiplying aperture functions instead convolves their transforms. These operations simplify complex apertures without repeatedly evaluating the full diffraction integral. It is the complex amplitudes, rather than their wave intensities, that are added for coherent illumination.
Convolution describes the response of a linear time-invariant system to an arbitrary input. If is its impulse response and its input,
This is a weighted sum of responses to shifted impulses. For integrable functions the operation is commutative and associative. A causal system has for , restricting the integral to past inputs. With the Fourier convention , the convolution theorem gives .
For example, an RC low-pass frequency response has and , where is the Heaviside step function. Fast input variations are smoothed. The same mathematics describes optical blurring by a point spread function and many other linear response measurements.
Deconvolution tries to recover from measured and known . In noise-free data one formally sets where . Zeros of lose corresponding information; small values amplify measurement noise. If the measured spectrum is , the direct estimate has error , which can be very large at frequencies strongly suppressed by the instrument. An inverse is therefore not automatically unique or stable.
One regularized spectral deconvolution reconstruction is
This minimizes at each frequency and limits noise amplification, at the cost of bias and lost resolution. Choosing the regularization requires assumptions about noise and the desired signal; deconvolution cannot restore components that the data do not determine without such additional information.
There are three independent nonfactorized momentum integrals, namely a tadpole integral, the three-line sunset diagram integral, and the five-line banana integral. The six displayed Feynman diagrams involve six complete integral expressions:
The distinction between independent building blocks and their products removes an ambiguity in counting “different momentum integrals”. For the strict divergent-wave-function renormalization interpretation of ii, the graphs in c reduce to , comprising three full integral expressions and two independent nonfactorized building blocks.
Make a Wick rotation and use to separate ultraviolet divergences from infrared divergences. Write , and let be the Euclidean versions of the integrals above, with their harmless overall factors of removed. A spherical ultraviolet cutoff gives
Thus the one-vertex quartic graph is linearly divergent, and the one-vertex sextic graph is quadratically divergent as .
For the other two building blocks, the spatial Fourier transform is especially transparent. The three-dimensional Euclidean scalar propagator is
The convolution theorem converts the parallel lines into . Angular integration gives, with a short-distance cutoff ,
This is a different convenient cutoff regularization for the connected block; power-divergent coefficients depend on the cutoff convention, whereas the conclusions about powers and the presence of a kinetic logarithm do not. The sinc function is understood to be at .
For , the integrand starts as , so
Its divergent part is independent of external momentum. In particular,
The finite difference and its finite derivative show why the pure three-line sunset diagram does not need a divergent wave-function renormalization in three dimensions.
For , the small-distance Taylor expansion is
Consequently
In particular, the nonzero term requires a genuine divergent wave-function renormalization, in addition to the mass counterterm. The coefficients here refer to the displayed position-space cutoff and have not included Feynman-diagram symmetry factors or bare couplings.
The comparison with the superficial degree of divergence is
Thus the actual leading power agrees with in every case. The additional logarithms in products come from ultraviolet subdivergences, which uniform power counting alone does not resolve.
For and , differentiation with respect to removes the logarithmic divergence of , but leaves divergent factors or multiplying a finite, nonpolynomial function of . These are subdivergent contributions, not new nonlocal counterterms: subtracting each tadpole diagram subdivergence and the momentum-independent subdivergence removes them. For example, in a factorized subtraction scheme , and likewise for . The remaining kinetic logarithm explains the strict interpretation in ii. A dimensional regularization scheme can suppress power divergences; it does not invalidate this cutoff-based comparison with .
For integrable functions , the convolution theorem in the chosen convention is
Let be the Fourier transform of in the horizontal coordinate. Transforming Laplace's equation gives
For , exclude the growing exponential and obtain . At , continuity in frequency gives ; a constant zero-frequency transform is compatible with spatial decay of .
By the transform pair in part (a), the inverse transform of is the Poisson kernel for the upper half-plane, . Thus
The Poisson extension of exponential boundary data is harmonic, approaches the continuous boundary data and decays at spatial infinity. For uniqueness, apply the maximum principle for harmonic functions to the difference of two solutions on increasingly large half-discs: it is zero on the real boundary and uniformly small on the distant semicircle.
The displayed Poisson integral is the bounded harmonic function on the upper half-plane with continuous boundary values . Taking a Fourier transform in reduces Laplace's equation to , and the bounded solution has multiplier . The convolution theorem gives the displayed expression. At frequency zero the transformed function is constant in , not zero: its value is the total boundary integral, two.
The function decays uniformly as . For large , use and the boundary L1 norm two. When is bounded and is large, split the convolution into and its complement: the kernel is small on the first part, and the exponential boundary value is uniformly small on the second. This also distinguishes spatial decay of the harmonic function from decay of its zero-frequency Fourier coefficient.