## Abstract

The scattering of a partially-coherent wave from a statistically rough material surface is investigated via derivation of the scattered field cross-spectral density function. Two forms of the cross-spectral density are derived using the physical optics approximation. The first is applicable to smooth-to-moderately rough surfaces and is a complicated expression of source and surface parameters. Physical insight is gleaned from its analytical form and presented in this work. The second form of the cross-spectral density function is applicable to very rough surfaces and is remarkably physical. Its form is discussed at length and closed-form expressions are derived for the angular spectral degree of coherence and spectral density radii. Furthermore, it is found that, under certain circumstances, the cross-spectral density function maintains a Gaussian Schell-model form. This is consistent with published results applicable only in the paraxial regime. Lastly, the closed-form cross-spectral density functions derived here are rigorously validated with scatterometer measurements and full-wave electromagnetic and physical optics simulations. Good agreement is noted between the analytical predictions and the measured and simulated results.

© 2013 OSA

## 1. Introduction

Rough surface scattering has been an active area of research for a half century. The published research on the subject can generally be divided into two main groups. The first deals with rough surface scattering research performed by the RF/microwave community for synthetic aperture radar and remote sensing applications. Predominately concerned with fully-coherent, monochromatic plane-wave scattering from rough surfaces (with some exceptions), the common approaches employed by the RF/microwave community to the rough surface scattering problem are physical optics (PO) [1–7], perturbation [4, 8], and computational/full-wave [9–14] methods. The interested reader is referred to Elfouhaily and Guérin [15], Warnick and Chew [16], Maradudin [17], Ogilvy [18], and Nieto-Vesperinas [19] for excellent summaries of these techniques.

The second of the two main groups deals with rough surface scattering research performed by the optics community. Initially, this work mainly dealt with how incoherent light interacted with surfaces, modeled via the bidirectional reflectance distribution function (BRDF), for applications in passive visible/near-IR remote sensing and computer graphics [20–24]. More recently, with the proliferation of laser-based systems (LIDAR/LADAR and directed energy), the interaction of coherent laser light with rough surfaces, in particular the statistical behaviors of the resulting speckle patterns, gained considerable interest in such fields as metrology and remote sensing. Some of the early notable research in this area was performed by Dainty [25], Goodman [26], Parry [27], Fujii and Asakura [28,29], Pedersen [30], and Yoshimura *et al.*[31].

Since the presence of speckle is typically detrimental in applications involving coherent light, techniques for suppressing speckle naturally followed. While there are many such techniques [32], due predominately to the work of Wolf [33–35], the use of partially-coherent light, instead of laser light, in active illumination systems is becoming tremendously popular for speckle suppression. In particular, much literature is dedicated to the properties of partially-coherent light whose cross-spectral density (CSD) function possesses a Gaussian Schell-model (GSM) form [33, 35–39]. In regards to the scattering of partially-coherent light, most of the current literature deals with scattering from low-contrast scatterers, i.e., for scatterers in which the Born approximation is valid [33, 35, 36, 40–45]. Far less work has been performed analyzing the scattering of partially-coherent light from rough metallic surfaces. Of the work that has been published, the following approaches are common: the phase-screen model [32, 46–48], ABCD matrices [47–51], and the coherent-mode representation [52].

The purpose of this rough surface scattering work is to extend the traditional, fully-coherent approaches to cases involving partially-coherent illumination. Previous work by the authors derived analytical forms for the scattering of a partially-coherent beam from a statistically rough perfectly-reflecting surface using the PO approximation [53]. In this work, the previously derived scattering solutions are generalized to material surfaces and rigorously verified via experiment and simulation.

Two forms of the scattered field CSD function are derived and discussed. The first, applicable to surfaces of smooth-to-moderate roughness, is expressed in terms of an infinite series. While its behavior depends in a complex way on source and surface parameters, its analytical form is physically intuitive. The second form of the scattered field CSD function is applicable to very rough surfaces. This form of the CSD function is incredibly physical and, under certain circumstances, maintains a GSM form, which is in agreement with the literature cited above valid in the paraxial regime. As such, closed-form expressions are derived for the angular spectral degree of coherence (SDoC) and spectral density (SD) radii. These expressions are, in general, complicated functions of both the source and surface parameters. It is demonstrated that for many scenarios of interest, the SDoC radius can be safely approximated as a function of just the source parameters and the SD radius can be simplified to a function of just the surface parameters.

To verify the theoretical analysis, experimental and Monte Carlo simulation results using a full-wave electromagnetic technique (all the physics of the wave/surface interaction are included) and one based on the PO approximation are presented and compared to the predictions of the analytical models. This paper is concluded with a summary of the work and contributions presented.

## 2. Methodology

The scattering geometry utilized in this analysis is shown in Fig. 1. The surface height is described by the function *h*(*x*) with mean, standard deviation, and correlation length equal to 0, *σ _{h}*, and

*ℓ*, respectively. The surface (of length 2

_{h}*L*) is illuminated by a partially-coherent beam (parameters

*w*and

_{s}*ℓ*defined below) emanating from the source plane specified by the coordinates (

_{s}*u*,

*v*,

*z*). Its origin, relative to the rough surface coordinate system, is located at the point (−

*x*,

_{s}*y*). The observation vector

_{s}**=**

*ρ***x̂**

*x*+

**ŷ**

*y*points from the rough surface origin to the observation point. The medium below the rough interface is electrically described by the permittivity

*ε*and permeability

*μ*; the medium above the rough interface is vacuum (

*ε*

_{0},

*μ*

_{0}).

A GSM form for the incident field CSD is used to model the partially-coherent illumination, viz.,

*w*and

_{s}*ℓ*are the source radius and source coherence radius, respectively and the functional dependence of

_{s}*W*,

^{i}*E*

_{0},

*w*, and

_{s}*ℓ*on the radian frequency

_{s}*ω*is omitted for brevity [33, 35]. Both s-pol and p-pol scattering solutions are derived below. Note that the rough surface and the source excitation are invariant in the

*z*direction resulting in a two-dimensional scattering problem. As mentioned by Johnson [17], light scattered from one-dimensional surfaces demonstrates the same physical behaviors as light scattered from two-dimensional surfaces. The mechanisms which are not captured in two-dimensional scattering problems are cross-polarized scattering and out-of-plane scattering [17]. In this way, the expressions derived in this paper are equivalent to scalar-wave scattering solutions.

#### 2.1. PO expression for the scattered field

For the sake of brevity, the analysis to follow focuses on perpendicular polarization. The solution for parallel polarization is provided at the end of the section. For perpendicular polarization, the incident field in the source plane takes the form

Utilizing the plane-wave spectrum representation of electromagnetic fields [54], the incident electric and magnetic fields become

*ρ**is a vector that points from the source plane origin to the rough surface origin, ${\mathbf{k}}^{i}=\widehat{\mathbf{u}}{k}_{u}^{i}+\widehat{\mathbf{v}}{k}_{v}^{i}$, and*

_{s}**k**

*| =*

^{i}*k*

_{0}= 2

*π*/

*λ*, thus ${k}_{v}^{i}=\sqrt{{k}_{0}^{2}-{\left({k}_{u}^{i}\right)}^{2}}$.

The scattered field in the far zone can be found by utilizing the transverse components of the far-field vector potentials **N** and **L**[55] and the equivalent electric **J** and magnetic **M** currents induced on the rough surface by the incident field, namely,

**′ =**

*ρ***x̂**

*x*′ +

**ŷ**

*h*(

*x*′). Note that the integral is over the parameterized surface contour

*C*with d

*c*equal to the arc length, i.e., $\text{d}c=\text{d}{x}^{\prime}\sqrt{1+{\left[{h}^{\prime}\left({x}^{\prime}\right)\right]}^{2}}$.

Using the PO approximation, the equivalent currents induced by the incident field on the material surface are

**n̂**is the unit outward normal to the surface given by

*r*

_{⊥}is the s-pol Fresnel complex amplitude reflection coefficient. Substitution of Eqs. (6) and (7) into Eq. (5) and subsequent simplification yields the following expression for ${E}_{z}^{s}$:

**=**

*ν**k*

_{0}

**̂ −**

*ρ***k**

*and*

^{i}*r*

_{⊥}is primed because it is a function of surface coordinates (

*x*′,

*h*(

*x*′)).

No analytical solution has yet been found for the integral over the rough surface coordinate *x*′ in Eq. (8) because of *ℋ*_{⊥}, which is a complicated function of *h*(*x*′) and its derivative *h _{x}*

_{′}[1]. This integral is typically evaluated using the stationary-phase approximation in which

*h*

_{x}_{′}≈ −

*ν*/

_{x}*ν*[1,19]. This approximation physically dictates that reflection from the rough surface is locally specular, i.e., local diffraction effects are excluded [1, 17–19, 23]. Applying this approximation to Eq. (8) and rearranging the integrals results in

_{y}#### 2.2. Scattered field cross-spectral density function

Applying the definition of the autocorrelation of a random process to Eq. (10) yields

*x*′

_{1},

*ρ*_{1}, ${k}_{u1}^{i}$ or

*x*′

_{2},

*ρ*_{2}, ${k}_{u2}^{i}$ (whichever is applicable), respectively, and the expectations are computed over the ensembles of rough surface and incident field realizations. In arriving at the above expression, it has been assumed that the incident field plane-wave spectrum is statistically independent of the rough surface; this assumption is physically intuitive. Note that the autocorrelation of the incident field plane-wave spectrum is equivalent to the Fourier transform of the CSD function given in Eq. (1). The second autocorrelation expression in Eq. (11), after some simple manipulation, is equivalent to the joint characteristic function of the random variables

*h*(

*x*′

_{1}) and

*h*(

*x*′

_{2}) [23]. The incident field term and the rough surface term are represented as $\mathrm{\Phi}\left({k}_{u1}^{i},{k}_{u2}^{i}\right)$ and $\chi \left({k}_{u1}^{i},{k}_{u2}^{i};{x}_{1}^{\prime},{x}_{2}^{\prime}\right)$ in the analysis to follow.

The integrals over the incident field plane-wave spectrum can be approximated using the method of stationary phase [33]. Using the method of stationary phase to evaluate the plane-wave spectrum integrals has two implications. The first is that
${k}_{v}^{i}\gg {k}_{u}^{i}$ and is therefore approximated as
${k}_{v}^{i}\approx {k}_{0}-{\left({k}_{u}^{i}\right)}^{2}/\left(2{k}_{0}\right)$ for all “phase” terms (i.e., exp[j(⋯)]) and as
${k}_{v}^{i}\approx {k}_{0}$ for all amplitude terms. This physically implies that the incident field is highly directional being predominately directed along the *v* direction in Fig. 1. The second implication is that the distance from the rough surface to the source plane *ρ _{s}* must be much greater than

*L*—this is generally the case. To provide some idea of how much greater

*ρ*must be than

_{s}*L*,

*ρ*= 10

_{s}*L*,

*ρ*= 20

_{s}*L*, and

*ρ*= 50

_{s}*L*are approximately 95%, 97.5%, and 99% accurate, respectively. Note that these two implications mean that the ${k}_{u}^{i}\widehat{\mathbf{v}}$ term in

*ℋ*

_{⊥}, see Eq. (9), diminishes at the rate of 1/

*ρ*and therefore can be neglected. Recall that this term corresponded to the

_{s}**v̂**component of the incident magnetic field implying that the field incident on the rough surface is transverse electromagnetic (in particular, TEM

*) when the method of stationary phase is utilized. Applying the method of stationary phase and substituting in $\mathrm{\Phi}\left({k}_{u1}^{i},{k}_{u2}^{i}\right)$ produces*

^{v}To simplify Eq. (12) further, a form for the rough surface characteristic function *χ* must be chosen. A very common choice for the statistical distribution of the rough surface is to assume that the surface heights are Gaussian distributed and Gaussian correlated, i.e.,

*χ*into Eq. (12), performing the common variable transformation

*x*=

_{d}*x*′

_{1}−

*x*′

_{2}and

*x*=

_{a}*x*′

_{1}+

*x*′

_{2}, and subsequent simplification yields

At this point it is necessary to handle the exponential term containing the surface autocorrelation function. Classically, this has been done in one of two ways. The first is to expand the exponential term in a Taylor series and proceed with the evaluation of the integrals [23]. Mathematically, this approach is applicable to all surfaces; however, computationally (because the series is slowly convergent), this approach is limited to smooth-to-moderately rough surfaces and is investigated in the next section. The other approach involves expanding the surface autocorrelation function in a Taylor series and only retaining the first two terms, i.e., $\text{exp}\left(-{x}_{d}^{2}/{\ell}_{h}^{2}\right)\approx 1-{x}_{d}^{2}/{\ell}_{h}^{2}$ [3, 23]. This treatment is applicable to very rough surfaces and is discussed in Section 2.2.2.

### 2.2.1. Smooth-to-moderately rough surfaces

As briefly discussed in the previous section, further progress can be made on Eq. (15) by expanding the exponential term containing the surface autocorrelation function in a Taylor series, specifically,

*x*results in

_{d}In order to evaluate the remaining integral and arrive at a closed-form solution, the complex error function (erf) must be properly handled. While it is possible to determine the order of the error function argument by using physical insight and examining the associated exponential functions in Eq. (17), it is much simpler to examine the integrand in Eq. (15) and determine the condition at which the *x _{d}* integration limits can be extended from [|

*x*| − 2

_{a}*L*, 2

*L*− |

*x*|] to (−∞, ∞). A necessary condition for this approximation occurs when

_{a}*δ*is a user-defined parameter and denotes the point at which the exponential function no longer has significant value. Taking note that the minimum of the exponential function’s argument occurs when

*m*= 0 and that the maximum value

*x*takes on in the integration is 2

_{d}*L*, the following condition can be derived:

*x*integration limits in Eq. (15) [|

_{d}*x*| − 2

_{a}*L*, 2

*L*− |

*x*|] ≈ (−∞, ∞). How well the incident beam “fits” is determined by

_{a}*δ*—the smaller the

*δ*, the more accurate the approximation. If the projected fully-coherent incident beam size does not “fit” on the rough surface, the erf ≠ 1 and the remaining expression must be evaluated numerically.

Assuming that the above condition holds, the remaining integral in Eq. (17) can be evaluated, namely,

*α*=

*ℓ*/

_{s}*w*. This condition physically means that the projected partially-coherent beam size must “fit” on the rough surface in order to safely assume that the sum of the erf’s ≈ 2. This condition is more stringent than Eq. (19).

_{s}The physical interpretation of Eq. (20) is obscured somewhat by the summation; however, the exponential term on the second line of the expression generally drives the angular extent of the scattered SD. The first exponential term on the third line generally determines the angular extent over which the scattered field is correlated, i.e., the SDoC radius. The SD *S* and SDoC *μ* definitions can be found in Refs. [33, 35, 36]. When one considers a fully-coherent incident Gaussian beam (*ℓ _{s}* → ∞) and quasi-monochromatic light [58], the coherent scattered irradiance, i.e.,
${I}_{c}^{s}={\left|\u3008{E}_{z}^{s}\u3009\right|}^{2}$, can be determined quite easily from Eq. (20) by considering just the

*m*= 0 term of the series and setting

*ρ*_{1}=

*ρ*_{2}=

**. With some minor algebraic manipulation, one arrives at the two-dimensional form of the coherent scattered irradiance expression derived by Wang**

*ρ**et al.*[2]. One can also derive from Eq. (20) the incoherent scattered irradiance ${I}_{i}^{s}$ for a fully-coherent incident Gaussian beam and quasi-monochromatic light. Note that ${I}_{i}^{s}$ is equivalent to the variance of the scattered field, i.e., ${I}_{i}^{s}={I}^{s}-{I}_{c}^{s}$, where ${I}^{s}=\u3008{E}_{z}^{s}{E}_{z}^{s*}\u3009$ is the average scattered irradiance. This expression also simplifies to the form derived by Wang

*et al.*[2].

### 2.2.2. Very rough surfaces

Consider the form of the joint characteristic function:

*χ*only has significant value when ${\vartheta}_{y1}/{\vartheta}_{y2}+{\vartheta}_{y2}/{\vartheta}_{y1}-2\text{exp}\left(-{x}_{d}^{2}/{\ell}_{h}^{2}\right)\approx 0$. Considering that

*ϑ*

_{y}_{1}/

*ϑ*

_{y}_{2}≈

*ϑ*

_{y}_{2}/

*ϑ*

_{y}_{1}, ${\vartheta}_{y1}+{\vartheta}_{y2}+{\vartheta}_{y2}/{\vartheta}_{y1}-2\text{exp}\left(-{x}_{d}^{2}/{\ell}_{h}^{2}\right)\approx 0$ is only possible for small

*x*. Thus, it makes sense to expand $\text{exp}\left(-{x}_{d}^{2}/{\ell}_{h}^{2}\right)$ and retain the first two terms, i.e., $\text{exp}\left(-{x}_{d}^{2}/{\ell}_{h}^{2}\right)\approx 1-{x}_{d}^{2}/{\ell}_{h}^{2}$. Substituting this into Eq. (15) and carrying out the integrations, one arrives at the desired result for the scattered field CSD function:

_{d}The above expression for the scattered field CSD function is remarkably physical. The first exponential term on the second line of Eq. (23) is predominately responsible for the angular extent of the scattered SD. This exponential term is a function of the sum of the squares of the projected observation angles *ϑ*_{x1} and *ϑ*_{x2}. The exponential term on the third line of Eq. (23) determines the SDoC radius. Note that this term is a function of the difference of the projected observation angles, i.e.,
$\left|{\vartheta}_{x1}-{\vartheta}_{x2}\right|=\left|\text{sin}{\theta}_{1}^{r}-\text{sin}{\theta}_{2}^{r}\right|$. Thus, it could be stated that the scattered field CSD function maintains its GSM form with respect to *ϑ*_{x1} and *ϑ*_{x2} if Eq. (21) is satisfied.

Because of the magnitude of the argument of the “correlation” exponential (specifically, the ${k}_{0}^{2}{\rho}_{s}^{4}{\sigma}_{h}^{2}{\vartheta}_{y1}{\vartheta}_{y2}$ term), $\text{sin}{\theta}_{1}^{r}\approx \text{sin}{\theta}_{2}^{r}$, implying that ${\theta}_{1}^{r}\approx {\theta}_{2}^{r}$, for the exponential to have a significant value. This implies that the correlation exponential is approximately a function of $\mathrm{\Delta}{\theta}^{r}={\theta}_{1}^{r}-{\theta}_{2}^{r}$. Since this term is predominately responsible for the behavior of the SDoC, the SDoC is also a function (approximately) of $\mathrm{\Delta}{\theta}^{r}={\theta}_{1}^{r}-{\theta}_{2}^{r}$. This is in agreement with the findings of previous studies which were restricted to the paraxial regime.

#### 2.3. Angular spectral degree of coherence radius

In this section, a theoretical form for the angular SDoC radius is derived from Eq. (23). Note that because of the summation in Eq. (20), it is not possible to derive an expression for the angular SDoC radius using that form of the scattered field CSD function.

Assuming Eq. (21) is satisfied and setting the exponential term on line three of Eq. (23) equal to 1/*e* results in the following expression:

*α*=

*ℓ*/

_{s}*w*(previously defined), Ω =

_{s}*w*/

_{s}*ρ*is the half-angle subtended by the source viewed from the rough surface, and ${\sigma}_{{h}^{\prime}}=\sqrt{2}{\sigma}_{h}/{\ell}_{h}$ is the surface slope standard deviation [22, 23]. Note that the same procedure applied to

_{s}*ϑ*

_{x1}and

*ϑ*

_{x2}(discussed above) has been applied to

*ϑ*

_{y}_{1}and

*ϑ*

_{y}_{2}in deriving the above expression.

The term involving Ω and *σ _{h}*

_{′}can be neglected when one considers that typical values for these parameters range from approximately 10

^{−4}rad and 0.05–0.5 rad, respectively. Note that the first term under the radical contains only source parameters, while the second term contains mostly rough surface parameters. In many scenarios of interest the source term is much greater than the rough surface term; thus, factoring out the source term and expanding the resulting radical in a Taylor series yields

The first term is entirely composed of source parameters and therefore represents the source contribution to the angular SDoC radius. Because of the presence of *σ _{h}*

_{′}, the second term provides a small correction to the angular SDoC radius due to the roughness of the surface. In most cases, the angular SDoC radius can be safely approximated by utilizing only the first term making the expression dependent only on the properties of the source illumination. This is consistent with the classic narrow-band, fully-coherent result derived by Goodman [26]. Note that when

*θ*≈

^{r}*θ*, i.e., when observation is in the specular direction, the dependence of the expression on observation and incident angles disappears.

^{i}#### 2.4. Angular spectral density radius

As is the case with the angular SDoC radius, an expression for the angular SD radius can only be derived using the very rough surface form of the scattered field CSD function in Eq. (23). Furthermore, the expression can only be derived for near-normal incidence.

Assuming Eq. (21) is satisfied, letting *θ ^{i}* = 0 and
${\theta}_{1}^{r}={\theta}_{2}^{r}={\theta}_{1/e}^{r}$, and setting the first exponential term on line two of Eq. (23) equal to 1/

*e*yields the following expression after some simplification:

*ψ*root makes physical sense), and utilizing the quadratic formula, yields the following roots for

*ψ*:

In Eq. (30), the first term in the argument of the inverse cosine depends solely on the characteristics of the rough surface and it physically represents the rough surface contribution to the angular extent of the SD. The second argument contains terms dealing with the source parameters and provides a small correction to the angular SD radius due to the size and coherence of the source. In most cases, the angular SD radius can be approximated using only the first term making the expression solely dependent on the roughness of the surface.

#### 2.5. Parallel-polarization cross-spectral density function

Considering that the rough surface’s features are large compared to wavelength, intuition dictates that the p-pol scattered field should have the same general behavior as the s-pol scattered field. This is the case. The only difference between the s-pol and p-pol CSD functions is the parameter *𝒮*. Setting *𝒮*^{⊥} = *𝒮*^{‖} in Eqs. (20) and (23) yields the p-pol scattered field CSD functions. Expressions for *𝒮*^{‖} and *𝒮*^{⊥} can be found in the appendix.

## 3. Validation of analytical solutions

In this section, the analytical solutions for the scattered field CSD functions derived above are validated via experiments, full-wave electromagnetic simulations, and PO simulations. For the experimental validation, the Complete Angle Scatter Instrument (CASI) [60] at the Air Force Institute of Technology was used. Measurements of the scattered SD *S ^{s}* versus

*θ*from a LabSphere Infragold [61] 5.08 cm × 5.08 cm coupon were made at

^{r}*λ*= 3.39

*μ*m (HeNe MWIR laser) and

*θ*= 20

^{i}*°*, 40

*°*, and 60

*°*. Before the scattering measurements were made, the surface statistics of the Infragold coupon were determined using a KLA Tencor Alpha-Step IQ surface profiler [62]. Four 1 cm scans (step size 0.2

*μ*m), performed along different directions, were taken and analyzed. The measured surface height standard deviation

*σ*, surface slope standard deviation

_{h}*σ*

_{h}_{′}, and correlation length

*ℓ*were 11.09

_{h}*μ*m, 0.2441 rad, and 116.9

*μ*m, respectively.

For the full-wave electromagnetic simulations, the Method of Moments (MoM) [63] was used to solve the following impedance-boundary-condition-based s-pol electric field integral equation for the unknown electric current **J**:

**A**and

**F**are the magnetic and electric vector potentials, respectively [64]. Once

**J**had been found, the scattered field was determined via convolution with the far-field form of the free-space Green’s function. A narrow-band, fully-coherent Gaussian beam was used for the incident field ${E}_{z}^{i}$.

*S*versus

^{s}*θ*and the modulus of the SDoC |

^{r}*μ*| versus

^{s}*Δθ*were calculated using the scattered field predicted from 500 rough surface realizations at

^{r}*λ*= 3.39

*μ*m, the Infragold surface statistics stated above, and a complex index of refraction

*n*− j

*κ*for gold of 1.995 − j20.95 [65]. The 500 rough surface realizations were generated using the method outlined by Yura and Hanson [66].

Both the experimental and full-wave simulations used narrow-band, fully-coherent incident fields. To verify the predictions of the analytical solutions for partially-coherent incident fields, PO Monte Carlo simulations were also performed. To model the random nature of the incident field, 500 source plane field instances were generated using the phase screen method described by Xiao and Voelz [67]. For the random surface, 500 rough surface instances were generated using the technique detailed in Ref. [66]. The 500 source plane field instances were propagated to and evaluated at each rough surface instance. The scattered fields were then computed numerically using the relations in Eq. (5). The end result of this process was a set of 500^{2} scattered fields, from which |*μ ^{s}*| and

*S*were computed. These simulations were performed at

^{s}*λ*= 3.39

*μ*m using hypothetical rough aluminum surfaces

*n*− j

*κ*= 5.366 − j33.33 [65] with varying degrees of surface roughness. In addition, the properties of the incident field (size and coherence) were also varied to examine their effects.

#### 3.1. Results

Figure 2 shows the normalized *S ^{s}* CASI measurement, MoM simulation, PO analytical solution, and PO simulation results versus

*θ*. The simulations were setup to best match the CASI experimental scenario, i.e.,

^{r}*ρ*= 185 cm and

_{s}*w*≈ 0.5 mm. For MoM computational reasons,

_{s}*L*= 2.5 mm, which resulted in an incident beam which over-illuminated the rough surface by less than 2%. While all the traces possess a Gaussian-like shape, there is a significant difference between the measured and MoM trace widths versus the PO results. As briefly discussed in Section 2.2, previous work had found that the SE-SE model more accurately represented surfaces roughened by random industrial processes [57]. If the Infragold sample measured and analyzed here adhered to the G-G model of rough surfaces, ${\sigma}_{{h}^{\prime}}=\sqrt{2}{\sigma}_{h}/{\ell}_{h}\approx 0.1342$ rad. The measured

*σ*

_{h}_{′}was roughly twice that value. Thus, in this case, the G-G model is not an ideal choice. Note that in these MoM simulation results, SE-SE surfaces were used, hence the good agreement with the CASI results. SE-SE surfaces could have been used in the PO simulations; however, as discussed previously, it would not have been possible to present analytical results. These results are presented for two main reasons. The first is to verify the analytical solutions which, based on the excellent agreement between the PO simulation and PO analytical results, is clearly accomplished. The other is to demonstrate the validity of the ubiquitous G-G model of rough surfaces. As can be seen, the model tends to underpredict the width of the scattered SD. Nevertheless, considering that analytical scattering solutions are possible using the G-G model, this is likely an acceptable drawback.

To truly validate the PO analytical solutions, MoM simulations were also performed using hypothetical Infragold G-G surfaces. These results, along with PO simulation and PO analytical solution results, are shown in Fig. 3. The source distance *ρ _{s}*, source radius

*w*, and surface 1/2 length

_{s}*L*were the same as above. The figure reports the normalized

*S*versus

^{s}*θ*in Fig. 3(a) and

^{r}**|**

*μ*| versus Δ

^{s}*θ*in Figs. 3(b) and 3(c). In Figs. 3(a) and 3(b), the simulated Infragold surface possessed the same

^{r}*σ*and

_{h}*ℓ*values reported above. Since this hypothetical surface qualifies as a very rough surface, the theoretical angular SD and SDoC radii, given in Eqs. (30) and (26), respectively, are also plotted as vertical dashed line in Figs. 3(a) and 3(b). Note the excellent agreement between the MoM simulation, PO simulation, and analytical PO solution results. In Fig. 3(c), |

_{h}*μ*| is shown for hypothetical Infragold G-G surfaces which possessed

^{s}*σ*= 0.05

_{h}*λ*, 0.1

*λ*, and 0.25

*λ*with

*ℓ*held constant at 8

_{h}*λ*. Surfaces with these parameters do not qualify as very rough surfaces and therefore, the angular SDoC radii are not shown. Again, note the very good agreement among the results.

Having validated the PO analytical solutions (as well as the PO simulation procedure) via experiment and full-wave electromagnetic simulations which utilized narrow-band, fully-coherent incident fields, attention can now be turned to validating the PO solutions via PO simulations using partially-coherent incident fields. These results are reported in Figs. 4 and 5. Figure 4 shows the theoretical (solid curves) and simulated (circles) *S ^{s}* for

*σ*

_{h}_{′}= 0.01 rad, 0.05 rad, 0.1 rad, and 0.25 rad (with

*ℓ*= 20

_{h}*λ*, ${\sigma}_{h}={\ell}_{h}{\sigma}_{{h}^{\prime}}/\sqrt{2}$, and

*L*= 0.35 m) versus

*θ*at normal incidence. Figures 4(a)–4(d) depict these curves for

^{r}*α*= 2, 1, 0.5, and 0.25 (with

*w*= 5 mm), respectively, i.e., a coherent source to a relatively incoherent source. The vertical dashed lines in the figures mark the locations of ${\theta}_{1/e}^{r}$ given by Eq. (30). Note the excellent agreement between the simulated and theoretical predictions. Also note that although the coherence properties of the source are very different for the curves plotted in the figures, the

_{s}*S*for the very rough surface curves

^{s}*σ*

_{h}_{′}= 0.05 rad, 0.1 rad, and 0.25 rad are nearly identical. This physically implies that

*S*(or the angular spread of the spectral density) for very rough surfaces is driven by surface properties not source parameters. Recall that this was theoretically predicted.

^{s}Figure 5 shows the theoretical (solid traces) and simulated (circles) |*μ ^{s}*| for

*α*= 2, 1, 0.5, and 0.25 (with the same

*w*as above) versus

_{s}*Δθ*. Figures 5(a)–5(d) show these curves for

^{r}*σ*

_{h}_{′}= 0.01 rad, 0.05 rad, 0.1 rad, and 0.25 rad (with the same

*ℓ*,

_{h}*σ*, and

_{h}*L*as above), respectively, i.e., a smooth surface to a very rough surface. The vertical dashed lines in the figures mark the locations of |

*Δθ*|

^{r}_{1/e}given by Eq. (26). Recall that Eq. (26) is applicable to very rough surfaces, a condition not met by the surfaces whose results are shown in Fig. 5(a). Note the excellent agreement between the Monte Carlo simulation results and the theoretical predictions. Also note that |

*Δθ*|

^{r}_{1/e}depend almost exclusively on the coherence properties of the source

*α*, as predicted.

## 4. Conclusion

The scattering of a partially-coherent field from a statistically rough surface was investigated. This work significantly extended previous efforts by incorporating the effects of material parameters (*ε* and *μ*) in the analysis and by rigorously validating the derived solutions via experiment and extensive simulation. Two forms of the scattered-field CSD function were derived in this work using the PO approximation. The first, applicable to smooth-to-moderately rough surfaces, was represented as an infinite series. While being a rather complicated expression dependent on both source and surface parameters, physical insight was gleaned from its analytical form. The second form of the scattered field CSD function was applicable to very rough surfaces. This expression was examined at length to include derivations of the angular SDoC and SD radii. It was noted that under certain circumstances, this form of the CSD function maintained a GSM form in agreement with published results valid only in the paraxial regime. Lastly, the closed-form expressions for the scattered field CSD functions were rigorously validated with scatterometer measurements, as well as, MoM and PO Monte Carlo simulations. The analytical predictions were found to be in good agreement with the measurement and simulation results.

## Appendix

The symbols first used in Eq. (12) are

*r*

_{‖}is the p-pol Fresnel complex amplitude reflection coefficient. Note that

*r*

_{⊥,‖}= −1 and 1 for a perfect electric conductor (i.e., a perfect reflector) and a perfect magnetic conductor, respectively.

## Acknowledgments

This research was supported in part by an appointment to the Postgraduate Research Participation Program at the Air Force Institute of Technology administered by the Oak Ridge Institute for Science and Education through an interagency agreement between the U.S. Department of Energy and AFIT.

The views expressed in this paper are those of the authors and do not reflect the official policy or position of the U.S. Air Force, the Department of Defense, or the U.S. Government.

## References and links

**1. **F. T. Ulaby, R. K. Moore, and A. K. Fung, *Microwave Remote Sensing: Active and Passive* (Artech House, 1986, vol. 2).

**2. **M.-J. Wang, Z.-S. Wu, and Y.-L. Li, “Investigation on the scattering characteristics of Gaussian beam from two dimensional dielectric rough surfaces based on the Kirchhoff approximation,” Prog. Electromagn. Res. B **4**, 223–235 (2008) [CrossRef] .

**3. **A. Ishimaru, C. Le, Y. Kuga, L. A. Sengers, and T. K. Chan, “Polarimetric scattering theory for high slope rough surfaces,” Prog. Electromagn. Res. **14**, 1–36 (1996).

**4. **A. Ishimaru, *Wave Propagation and Scattering in Random Media* (IEEE, 1997).

**5. **G. Guo, S. Li, and Q. Tan, “Statistical properties of laser speckles generated from far rough surfaces,” Int. J. Infrared Millimeter Waves **22**, 1177–1191 (2001) [CrossRef] .

**6. **V. I. Mandrosov, *Coherent Fields and Images in Remote Sensing* (SPIE, 2004) [CrossRef] .

**7. **M. E. Knotts, T. R. Michel, and K. A. O’Donnell, “Angular correlation functions of polarized intensities scattered from a one-dimensionally rough surface,” J. Opt. Soc. Am. A **9**, 1822–1831 (1992) [CrossRef] .

**8. **S. O. Rice, “Reflection of electromagnetic wave by slightly rough surfaces,” Commun. Pure Appl. Math. **4**, 351–378 (1951) [CrossRef] .

**9. **A. K. Fung and K. S. Chen, *Microwave Scattering and Emission Models for Users* (Artech House, 2010).

**10. **R. Axline and A. Fung, “Numerical computation of scattering from a perfectly conducting random surface,” IEEE Trans. Antennas Propag. **26**, 482–488 (1978) [CrossRef] .

**11. **R. E. Collin, “Full wave theories for rough surface scattering: an updated assessment,” Radio Sci. **29**, 1237–1254 (1994) [CrossRef] .

**12. **R. Collin, “Scattering of an incident Gaussian beam by a perfectly conducting rough surface,” IEEE Trans. Antennas Propag. **42**, 70–74 (1994) [CrossRef] .

**13. **E. I. Thorsos, “The validity of the Kirchhoff approximation for rough surface scattering using a Gaussian roughness spectrum,” J. Acoust. Soc. Am. **83**, 78–92 (1988) [CrossRef] .

**14. **E. Bahar, “Full-wave solutions for the depolarization of the scattered radiation fields by rough surfaces of arbitrary slope,” IEEE Trans. Antennas Propag. **29**, 443–454 (1981) [CrossRef] .

**15. **T. M. Elfouhaily and C.-A. Guérin, “A critical survey of approximate scattering wave theories from random rough surfaces,” Waves Random Media **14**, R1–R40 (2004) [CrossRef] .

**16. **K. F. Warnick and W. C. Chew, “Numerical simulation methods for rough surface scattering,” Waves Random Media **11**, R1–R30 (2001) [CrossRef] .

**17. **A. A. Maradudin, ed., *Light Scattering and Nanoscale Surface Roughness* (Springer, 2007) [CrossRef] .

**18. **J. A. Ogilvy, *Theory of Wave Scattering from Random Rough Surfaces* (IOP Publishing, 1991).

**19. **M. Nieto-Vesperinas, *Scattering and Diffraction in Physical Optics*, 2 ed. (World Scientific, 2006) [CrossRef] .

**20. **K. E. Torrance and E. M. Sparrow, “Theory for off-specular reflection from roughened surfaces,” J. Opt. Soc. Am. **57**, 1105–1112 (1967) [CrossRef] .

**21. **R. G. Priest and S. R. Meier, “Polarimetric microfacet scattering theory with applications to absorptive and reflective surfaces,” Opt. Eng. **41**, 988–993 (2002) [CrossRef] .

**22. **Y. Sun, “Statistical ray method for deriving reflection models of rough surfaces,” J. Opt. Soc. Am. A **24**, 724–744 (2007) [CrossRef] .

**23. **P. Beckmann and A. Spizzichino, *The Scattering of Electromagnetic Waves from Rough Surfaces* (Pergamon, 1963).

**24. **X. D. He, K. E. Torrance, F. X. Sillion, and D. P. Greenberg, “A comprehensive physical model for light reflection,” SIGGRAPH Computer Graphics **25**, 175–186 (1991) [CrossRef] .

**25. **J. Dainty, “Some statistical properties of random speckle patterns in coherent and partially coherent illumination,” Opt. Acta **17**, 761–772 (1970) [CrossRef] .

**26. **J. Goodman, “Statistical properties of laser speckle patterns,” in “*Laser Speckle and Related Phenomena*,” vol. 9 of *Topics in Applied Physics*, J. C. Dainty, ed. (Springer-Verlag, 1975, pp. 9–75) [CrossRef] .

**27. **G. Parry, “Speckle patterns in partially coherent light,” in “*Laser Speckle and Related Phenomena*,” vol. 9 of *Topics in Applied Physics*, J. C. Dainty, ed. (Springer-Verlag, 1975, pp. 77–121) [CrossRef] .

**28. **H. Fujii and T. Asakura, “A contrast variation of image speckle intensity under illumination of partially coherent light,” Opt. Commun. **12**, 32–38 (1974) [CrossRef] .

**29. **H. Fujii and T. Asakura, “Statistical properties of image speckle patterns in partially coherent light,” Nouv. Rev. Opt. **6**, 5–14 (1975) [CrossRef] .

**30. **H. Pedersen, “The roughness dependence of partially developed, monochromatic speckle patterns,” Opt. Commun. **12**, 156–159 (1974) [CrossRef] .

**31. **T. Yoshimura, K. Kato, and K. Nakagawa, “Surface-roughness dependence of the intensity correlation function under speckle-pattern illumination,” J. Opt. Soc. Am. A **7**, 2254–2259 (1990) [CrossRef] .

**32. **J. W. Goodman, *Speckle Phenomena in Optics: Theory and Applications* (Ben Roberts & Company, 2007).

**33. **L. Mandel and E. Wolf, *Optical Coherence and Quantum Optics* (Cambridge University, 1995).

**34. **E. Wolf, “Unified theory of coherence and polarization of random electromagnetic beams,” Phys. Lett. A **312**, 263–267 (2003) [CrossRef] .

**35. **E. Wolf, *Introduction to the Theory of Coherence and Polarization of Light* (Cambridge University, 2007).

**36. **G. Gbur and T. Visser, “The structure of partially coherent fields,” Progress in Optics **55**, 285–341 (2010) [CrossRef] .

**37. **O. Korotkova, M. Salem, and E. Wolf, “The far-zone behavior of the degree of polarization of electromagnetic beams propagating through atmospheric turbulence,” Opt. Commun. **233**, 225–230 (2004) [CrossRef] .

**38. **O. Korotkova, T. Visser, and E. Wolf, “Polarization properties of stochastic electromagnetic beams,” Opt. Commun. **281**, 515–520 (2008) [CrossRef] .

**39. **T. Leskova, A. Maradudin, and J. Munoz-Lopez, “The design of one-dimensional randomly rough surfaces that act as Collett-Wolf sources,” Opt. Commun. **242**, 123–133 (2004) [CrossRef] .

**40. **J. Li, Y. Chen, S. Xu, Y. Wang, M. Zhou, Q. Zhao, Y. Xin, and F. Chen, “Condition for invariant spectral degree of coherence of an electromagnetic plane wave on scattering,” Opt. Commun. **284**, 724–728 (2011) [CrossRef] .

**41. **T. Wang and D. Zhao, “Stokes parameters of an electromagnetic light wave on scattering,” Opt. Commun. **285**, 893–895 (2012) [CrossRef] .

**42. **T. Shirai and T. Asakura, “Multiple light scattering from spatially random media under the second-order Born approximation,” Opt. Commun. **123**, 234–249 (1996) [CrossRef] .

**43. **G. Gbur and E. Wolf, “Determination of density correlation functions from scattering of polychromatic light,” Opt. Commun. **168**, 39–45 (1999) [CrossRef] .

**44. **F. Gori, C. Palma, and M. Santarsiero, “A scattering experiment with partially coherent light,” Opt. Commun. **74**, 353–356 (1990) [CrossRef] .

**45. **E. Wolf, J. T. Foley, and F. Gori, “Frequency shifts of spectral lines produced by scattering from spatially random media,” J. Opt. Soc. Am. A **6**, 1142–1149 (1989) [CrossRef] .

**46. **B. G. Hoover and V. L. Gamiz, “Coherence solution for bidirectional reflectance distributions of surfaces with wavelength-scale statistics,” J. Opt. Soc. Am. A **23**, 314–328 (2006) [CrossRef] .

**47. **L. C. Andrews and R. L. Phillips, *Laser Beam Propagation through Random Media*, 2 ed. (SPIE, 2005) [CrossRef] .

**48. **O. Korotkova, *Partially Coherent Beam Propagation in Turbulent Atmosphere with Applications* (VDM Verlag Dr. Müller, 2009).

**49. **R. S. Hansen, H. T. Yura, and S. G. Hanson, “First-order speckle statistics: an analytic analysis using ABCD matrices,” J. Opt. Soc. Am. A **14**, 3093–3098 (1997) [CrossRef] .

**50. **H. Yura and S. Hanson, “Variance of intensity for Gaussian statistics and partially developed speckle in complex ABCD optical systems,” Opt. Commun. **228**, 263–270 (2003) [CrossRef] .

**51. **S. Sahin, Z. Tong, and O. Korotkova, “Sensing of semi-rough targets embedded in atmospheric turbulence by means of stochastic electromagnetic beams,” Opt. Commun. **283**, 4512–4518 (2010) [CrossRef] .

**52. **J. Huttunen, A. T. Friberg, and J. Turunen, “Scattering of partially coherent electromagnetic fields by microstructured media,” Phys. Rev. E **52**, 3081–3092 (1995) [CrossRef] .

**53. **M. W. Hyde, S. Basu, S. J. Cusumano, and M. F. Spencer, “Scalar wave solution for the scattering of a partially coherent beam from a statistically rough metallic surface,” Proc. SPIE **8550** (2012) [CrossRef] .

**54. **T. Hansen and A. Yaghjian, *Plane-Wave Theory of Time-Domain Fields: Near-Field Scanning Applications* (IEEE, 1999) [CrossRef] .

**55. **S. Ramo, J. Whinnery, and T. Van Duzer, *Fields and Waves in Communication Electronics*, 3 ed. (John Wiley & Sons, 1994).

**56. **C. A. Mack, “Analytic form for the power spectral density in one, two, and three dimensions,” J. Micro/Nanolithogr. MEMS MOEMS **10**, 040501 (2011).

**57. **S. Basu, S. J. Cusumano, M. W. Hyde, M. A. Marciniak, and S. T. Fiorino, “Validity of using Gaussian Schell model for extended beacon studies,” Proc. SPIE **8380** (2012).

**58. **Under quasi-monochromatic conditions (*Δω* ≪ *ω̄*), *S* (** ρ**,

*ω*)

*≈*Γ(

**,**

*ρ***, 0)**

*ρ**δ*(

*ω*−

*ω̄*) and

*μ*(

*ρ*_{1},

*ρ*_{2},

*ω*) ≈

*γ*(

*ρ*_{1},

*ρ*_{2},

*τ*) exp (j

*ω̄τ*). Here, Γ is the mutual coherence function and

*γ*is the complex degree of coherence [59].

**59. **A. T. Friberg and E. Wolf, “Relationships between the complex degrees of coherence in the space-time and in the space-frequency domains,” Opt. Lett. **20**, 623–625 (1995) [CrossRef] [PubMed] .

**60. **B. Balling, “A comparative study of the bidirectional reflectance distribution function of several surfaces as a mid-wave infrared diffuse reflectance standard,” Master’s thesis, Graduate School of Engineering and Management, Air Force Institute of Technology (AETC), Wright-Patterson AFB OH (2009). http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA495933.

**61. **LabSphere Inc., “A guide to reflectance coatings and materials” (LabSphere, Inc., 2012). http://www.labsphere.com/uploads/technical-guides/a-guide-to-reflectance-materials-and-coatings.pdf.

**62. **KLA-Tencor Corporation, “Alpha-Step IQ Surface Profiler” (KLA-Tencor Corporation, 2012). http://www.kla-tencor.com/surface-profiling/alpha-step-iq.html.

**63. **R. Harrington, *Field Computation by Moment Methods* (IEEE, 1993) [CrossRef] .

**64. **A. F. Peterson, S. L. Ray, and R. Mittra, *Computational Methods for Electromagnetics* (IEEE, 1998).

**65. **RefractiveIndex.INFO, “RefractiveIndex.INFO Refractive index database” (RefractiveIndex.INFO, 2012). http://refractiveindex.info.

**66. **H. T. Yura and S. G. Hanson, “Digital simulation of an arbitrary stationary stochastic process by spectral representation,” J. Opt. Soc. Am. A **28**, 675–685 (2011) [CrossRef] .

**67. **X. Xiao and D. Voelz, “Wave optics simulation approach for partial spatially coherent beams,” Opt. Express **14**, 6986–6992 (2006) [CrossRef] [PubMed] .