Appendix C Closed-Form Solution for First-Order Diffraction

For first-order diffraction, the position of the diffraction point can be determined in closed form. Consider a source ๐ฌ, a target ๐ญ, and an edge defined by its direction ๐ž^. For convenience, we assume the edge passes through the origin of the coordinate system, denoted by ๐จ.

According to Fermatโ€™s principle, the diffraction point ๐ฏ on the edge minimizes the total path length ๐ฌโ†’๐ฏโ†’๐ญ. This corresponds to minimizing the function

โ„’โข(x)=โˆฅ๐ฌโˆ’๐ฏโˆฅ2+โˆฅ๐ฏโˆ’๐ญโˆฅ2 (183)

where

๐ฏ=xโข๐ž^. (184)

As shown inย [22], โ„’โข(โ‹…) is strictly convex and thus has a unique minimizer.

We introduce the set of vectors (๐ฎ^1,๐ฎ^2,๐ฎ^3), where

๐ฎ^1 =๐ฌโˆ’(๐ž^๐–ณโข๐ฌ)โข๐ž^โˆฅ๐ฌโˆ’(๐ž^๐–ณโข๐ฌ)โข๐ž^โˆฅ2 (185)
๐ฎ^2 =๐ญโˆ’(๐ž^๐–ณโข๐ฌ)โข๐ž^โˆฅ๐ญโˆ’(๐ž^๐–ณโข๐ฌ)โข๐ž^โˆฅ2 (186)
๐ฎ^3 ={๐ฎ^1ร—๐ฎ^2ย ifย โข๐ฎ^1ร—๐ฎ^2โ‰ ๐ŸŽ๐ž^ย ifย โข๐ฎ^1ร—๐ฎ^2=๐ŸŽ. (187)

Note that ๐ฎ^3=ยฑ๐ž^, and that the source and the edge lie within the plane (๐จ,๐ฎ^1,๐ž^). A key observation is that rotating the target ๐ญ about the edge does not alter the path length. Specifically, let ๐‘๐ฎ^3โข(ฯ•) denote the rotation matrix for an angle ฯ• around the edge (๐จ,๐ž^). Then,

โ„’โข(x) =โˆฅ๐ฌโˆ’๐ฏโˆฅ2+โˆฅ๐ฏโˆ’๐ญโˆฅ2 (188)
=โˆฅ๐ฌโˆ’๐ฏโˆฅ2+โˆฅ๐ฏโˆ’๐‘๐ฎ^3โข(ฯ•)โข๐ญโŸ๐ญโ€ฒโˆฅ2 (189)

holds for any ฯ•โˆˆโ„. Especially, by rotating the target around the axis (๐จ,๐ฎ^3) by the angle

ฯ•=ฯ€โˆ’arccosโก(๐ฎ^1๐–ณโข๐ฎ^2), (191)

the rotated target ๐ญโ€ฒ is placed in the same plane (๐จ,๐ฎ^1,๐ž^) as the source and the edge, but on the side opposite to the source with respect to the edge, as depicted in Figureย 36. The rotation matrix ๐‘๐ฎ^3โข(ฯ•) is given by Rodriguesโ€™ rotation formulaย [30].

Refer to caption
(a) View in the plane (๐จ,๐ฎ^1,๐ฎ^2), which is perpendicular to the edge ๐ž^.
Refer to caption
(b) View in the plane (๐จ,๐ฎ^1,๐ž^), containing ๐ฌ, the edge, and the rotated target ๐ญโ€ฒ.
Figure 36: Illustration of rotating the target ๐ญ about the edge ๐ž^ by the angle ฯ• to obtain the rotated target ๐ญโ€ฒ, which is coplanar with the source ๐ฌ and the edge. The diffraction point ๐ฏ is such that ๐ฌ, ๐ฏ, and ๐ญโ€ฒ are collinear.

Minimizing โ„’โข(x) with the rotated target ๐ญโ€ฒ is equivalent to ๐ฌ, ๐ฏ, and ๐ญโ€ฒ being collinear, i.e.

(๐ฏโˆ’๐ฌ)ร—๐ฌ๐ญโ€ฒโ†’ =๐ŸŽ (192)
โ‡”xโข(๐ž^ร—๐ฌ๐ญโ€ฒโ†’) =๐ฌร—๐ฌ๐ญโ€ฒโ†’ (193)

where ๐ฌ๐ญโ€ฒโ†’=๐ญโ€ฒโˆ’๐ฌ. Since ๐ž^ร—๐ฌ๐ญโ€ฒโ†’ and ๐ฌร—๐ฌ๐ญโ€ฒโ†’ are parallel, the solution for x is

x=signโก([๐ž^ร—๐ฌ๐ญโ€ฒโ†’]๐–ณโข[๐ฌร—๐ฌ๐ญโ€ฒโ†’])โขโˆฅ๐ฌร—๐ฌ๐ญโ€ฒโ†’โˆฅ2โˆฅ๐ž^ร—๐ฌ๐ญโ€ฒโ†’โˆฅ2. (194)