Free Access
Issue
A&A
Volume 529, May 2011
Article Number A18
Number of page(s) 13
Section Astronomical instrumentation
DOI https://doi.org/10.1051/0004-6361/201116430
Published online 22 March 2011

© ESO, 2011

1. Introduction

Optical modules for X-ray telescopes consist of a number of grazing incidence mirror shells with a common axis and focus. In a widespread design, the Wolter’s, the mirrors comprise two consecutive segments, a paraboloid and a hyperboloid, in order to concentrate X-rays by means of a double reflection. The two reflections occur at the same incidence angle for X-rays coming from an on-axis source at astronomical distance (Van Speybroeck & Chase 1972). The optical design of the module is primarily dependent on the required effective area, which determines the telescope’s sensitivity.

The most representative indicator of the concentration power of an optical module is assumed in general to be the effective area on-axis. However, the effective area in general decreases as the X-ray source moves off-axis. This in turn diminishes the telescope’s sensitivity, hence may represent a severe limitation to its field of view. For this reason, the prediction of the on- and off-axis effective area modules is a very important task in the development of new X-ray telescopes such as NuSTAR (Hailey et al. 2010), NHXM (Basso et al. 2010), ASTRO-H (Kunieda et al. 2010), and IXO (Bookbinder 2010).

While the effective area on-axis is in general easy to calculate for a typical Wolter-I mirror module configuration, the computation becomes more difficult off-axis, because of the variable incidence angles over the two reflecting surfaces and the variable fraction of singly-reflected X-rays that are not focused and contribute to the stray light (see e.g., Cusumano et al. 2007). To make things worse, the assembled mirrors can shade each other if they are not spaced enough, which contributes to an even steeper decrease in the collecting area off-axis. A widespread method for performing the calculation, accounting for all these factors, has made use of accurate ray-tracing routines (see, e.g., Mangus & Underwood 1969; Zhao et al. 2004), which reconstruct the paths of a selection of X-rays impinging a mirror module. These numerical codes are in general accurate, because they simulate the real incidence of rays on the optical system. However, they are complex and time consuming, especially whenever the simulation includes wide-band multilayer coatings to extend the reflectivity beyond 10 keV (Joensen et al. 1995; Tawara et al. 1998). The reason is that the multilayer reflectivity computation, which has to be performed for every ray traced, is a complex procedure especially for wideband multilayers, which comprise many (~200) couples of layers.

For this reason, although the ray-tracing approach should not be disregarded, it is interesting to derive analytical formulae for the effective area off-axis. In attempting to achieve this, Van Speybroeck & Chase (1972) discovered, by analyzing the results of ray-tracing simulations, that the geometric collecting area of a Wolter-I mirror decreases with the off-axis angle θ of the source, according to the formula A(θ)=A(0)(12θ3α0),\begin{equation} A_{\infty}(\theta) = A_{\infty}(0)\,\left(1-\frac{2\theta}{3\alpha_0}\right), \label{eq:SC_formula} \end{equation}(1)where A(0) is the on-axis geometric area and α0 is the incidence angle for an astronomical source on-axis.

More recently, a semi-analytical method for computing the off-axis effective area has been applied to solve the problem of optimizing multilayer recipes to the telescope’s field of view (Mao et al. 1999; Mao et al. 2000; Madsen et al. 2009). In this approach, the effective area is computed by integrating the mirror reflectivity over the incidence angles, after weighting it over an appropriate function, Winc, derived from a ray-tracing.

In a previous paper (Spiga et al. 2009), we already derived a completely analytical method for computing the off-axis effective area of a Wolter-I mirror with any reflective coating. In a subsequent article (Spiga & Cotroneo 2010), we developed this formalism, deriving the analytical expression of the aforementioned Winc function, which represents the distribution of the off-axis effective area over the incidence angles of rays, and we also used it to face the problem of multilayer optimization (Cotroneo et al. 2010). The results were accurately verified as well, by comparison with the outcomes of a ray-tracing program. However, in these previous works, we did not consider the mutual shading (also known as vignetting) of mirrors, which may occur in mirror modules when shells are densely nested together. While in general the mirror module is designed to avoid any vignetting on-axis, the problem may arise for sources off-axis and cause a further loss of effective area. Consequently, the results could be applied only to single mirror shells, or to mirror assemblies for which the mutual obstruction is known to be negligible over all the field of view.

In this work, we overcome these limitations and extend the developed formalism to the general case of obstructed Wolter-I mirrors in X-ray optical modules. We still assume that the mirror profile can be approximated by a double cone as far as the sole effective area is concerned, a condition in general fulfilled by optics with large f-numbers. In Sect. 2, we briefly review the results obtained for unobstructed single mirror shells, and in addition we show how the analytical formalism can be inverted to derive the product of the two reflectivities from the desired effective area variation with the off-axis angle. In Sect. 3, we describe the geometrical parameters driving the nested mirror obstructions, derive the expression of the vignetting coefficients, and obtain an analytical formula for the off-axis effective area of an obstructed mirror (Eq. (38)). We then derive in Sect. 4 some analytical expressions for the obstructed geometric area, and we apply the results to the geometrical optimization of the module. In Sect. 5, we prove the validity of the analytical formulae by means of a ray-tracing routine. The results are briefly summarized in Sect. 6.

2. Analytical formulae for unobstructed Wolter-I mirrors

We consider a Wolter-I mirror shell (Fig. 1) with focal length f. We denote with R0 the radius of the reflective surface at the intersection plane and with α0 the incidence angle at the intersection plane for an astronomical source on-axis. They are related by the well-known relation R0=ftan(4α0).\begin{equation} R_0 = f \tan\left(4\alpha_0\right). \label{eq:Rf} \end{equation}(2)We hereafter assume that all incidence angles are shallow, therefore tan(4α0) ≈ 4α0. The optical axis of the mirror is oriented in the z-axis direction. We define L1 and L2 to be the lengths of the parabolic and hyperbolic segments (hereafter named primary and secondary), which are supposed to be undeformed and very smooth. We denote with RM and Rm the mirror radii at the entrance and exit pupil, respectively. We define rλ(α) to be the coating reflectivity for a generic incidence angle α, at the X-ray wavelength λ, and assume that the geometric optics can be applied, i.e., that λ is small enough to avoid aperture diffraction effects but large enough for the scattering due to surface roughness to be negligible (Raimondi & Spiga 2010). The source is assumed to be at a finite, although very large, distance D ≫ L1,2. Finally, θ ≥ 0 denotes the off-axis angle, i.e., the angle formed by the source direction with the optical axis. The x-axis of this reference frame is oriented such that the source lies in the xz plane, at x > 0 and z > 0.

thumbnail Fig. 1

Sketch of a grazing-incidence Wolter-I mirror shell with an off-axis X-ray source. We also show a ray direction vector before reflection (k0)\hbox{$(\underline{k}_0)$}, after the first reflection (k1)\hbox{$(\underline{k}_1)$}, and after two reflections (k2)\hbox{$(\underline{k}_2)$}.

2.1. A review of previous results

We briefly review the results about the computation of the off-axis effective area of unobstructed Wolter-I mirror shells, derived in previous works (Spiga et al. 2009; Spiga & Cotroneo 2010).

  • 1.

    Owing to the shallow incidence angles, the double coneapproximation is in general applicable to compute the effectivearea of Wolter-I mirrors. The error in the effective area estimation,by replacing a Wolter-I profile with a double cone, is on the orderof L/f or smaller, i.e., a few percent in real cases.

  • 2.

    For a source off-axis, the incidence angles on the primary and secondary segment, α1 and α2, essentially depend on the polar angles ϕ1 and ϕ2 of the points where the ray is reflected (Fig. 1). As long as θ is small, ϕ1 ≈ ϕ2. We thus denote their nearly-common value with ϕ. If the polar angle ϕ is measured from the x-axis, α1 and α2 have the expressions α1=α0+δθcosϕ,α2=α0δ+θcosϕ,\begin{eqnarray} \label{eq:angle1}\alpha_1 &=& \alpha_0+\delta-\theta\cos\varphi,\\ \label{eq:angle2}\alpha_2 &=& \alpha_0-\delta+\theta\cos\varphi , \end{eqnarray}where δ = R0/D is the X-ray semi-divergence due to the distance of the source: for an astronomical source, δ ≃ 0. Equations (3) and (4) are valid if α1 ≥ 0, α2 ≥ 0. We note that α1+α2=2α0.\begin{equation} \alpha_1+\alpha_2=2\alpha_0. \label{eq:anglesum} \end{equation}(5)

  • 3.

    The geometric ratio of rays that undergo the double reflection to those impinging the primary segment is expressed by the vignetting factor V(ϕ)=L2α2L1α1,\begin{equation} V(\varphi) = \frac{L_2\alpha_2}{L_1\alpha_1}, \label{eq:2reflvign} \end{equation}(6)with the constraint that 0 < V(ϕ) < 1.

  • 4.

    Using Eq. (6), it is easy to derive an integral formula for the effective area of the mirror shell AD(λ,θ)=2R00π()minrλ(α1)rλ(α2)dϕ,\begin{equation} A_D(\lambda, \theta) = 2R_0 \int_0^{\pi}\!\!(L\alpha)_{\mathrm{min}} \,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right) \,\mbox{d}\varphi, \label{eq:Aeff_fin_offaxis} \end{equation}(7)where ()min=max{0,min[L1α1(ϕ),L2α2(ϕ)]}.\begin{equation} (L\alpha)_{\mathrm{min}} = \max\left\{0, \min\left[L_1 \alpha_1(\varphi), L_2 \alpha_2(\varphi)\right]\right\}. \label{eq:Lalpha_min} \end{equation}(8)

  • 5.

    In the frequent case L1 = L2   (= L), Eq. (7) reduces to AD(λ,θ)=2R0L0παminrλ(α1)rλ(α2)dϕ,\begin{equation} A_{D}(\lambda, \theta) = 2R_0 L\int_0^{\pi}\!\!\alpha_{\mathrm{min}}\,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right)\, \mbox{d}\varphi, \label{eq:Aeff_fin} \end{equation}(9)where αmin = min [α1(ϕ),α2(ϕ) ]  if positive, and zero otherwise. The expression of the αmin angle can be written in a compact form αmin=max(0,α0|δθcosϕ|),\begin{equation} \alpha_{\mathrm{min}} = \max\left(0,\,\alpha_0-|\delta-\theta\cos\varphi|\right), \label{eq:theta_min} \end{equation}(10)

  • 6.

    In the ideal case of a perfectly-reflecting mirror, i.e., setting r = 1 for any α and λ, the integral in Eq. (9) can be solved explicitly. For instance, if δ = 0 and θ = 0, we obtain A(0)=2πR0Lα0,\begin{equation} A_{\infty}(0) =2\pi R_0 L \,\alpha_0, \label{eq:Ageom_1} \end{equation}(11)whereas, if δ = 0 and 0 < θ < α0, A(θ)=A(0)(12θπα0),\begin{equation} A_{\infty}(\theta) = A_{\infty}(0)\, \left(1- \frac{2\theta}{\pi\alpha_0}\right), \label{eq:Ageom_2} \end{equation}(12)where – and heretofore – we drop the dependence of the A’s on λ to denote the geometric mirror areas. Equation (12), after approximating π ≈ 3, becomes Eq. (1), the expression found by Van Speybroeck & Chase (1972). For θ > α0, the geometric area has the more complicated expression A(θ)=A(0)[12π(θα0θ2α021+arccosα0θ)].\begin{equation} A_{\infty}(\theta) = A_{\infty}(0)\, \left[1-\frac{2}{\pi}\left(\frac{\theta}{\alpha_0}-\sqrt{\frac{\theta^2}{\alpha_0^2}-1}+\arccos\frac{\alpha_0}{\theta}\right)\right]. \label{eq:Ageom_3} \end{equation}(13)

  • 7.

    If δ > 0, more analytical expressions can be obtained by solving the integral in Eq. (9). The predictions are in very good agreement with the findings of detailed ray-tracing routines. For more details, we refer to Spiga et al. (2009).

  • 8.

    If δ = 0 and 0 < θ < α0, we obtain an alternative form of Eq. (9), by changing the integration variable from ϕ to α1, A(λ,θ)=4R0Lα0θα0α1rλ(α1)rλ(α2)θ2(α0α1)2dα1.\begin{equation} A_{\infty}(\lambda, \theta) =4 R_0 L \int_{\alpha_0-\theta}^{\alpha_0}\!\frac{\alpha_1\,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right)}{\sqrt{\theta^2-\left(\alpha_0-\alpha_1\right)^2}}\, \mbox{d}\alpha_1. \label{eq:Aeff_alpha} \end{equation}(14)This equation can be used to derive the effective area as a function of θ, at a fixed λ.

2.2. Inverse computation: from A(λ,θ)\hbox{$\mathsf{_{\infty}(\lambda, \theta)}$} to the mirror reflectivity

As a further development of the equations reported in the previous section for unobstructed mirrors, we hereafter provide the inverse formula of Eq. (14). For a given λ, from the effective area variation with θ in  [0,α0] , for δ = 0, we derive the product of the reflectivities of the two segments Tλ(α1)=defrλ(α1)rλ(α2),\begin{equation} T_{\lambda}\left(\alpha_1\right) \stackrel{\mathrm{def}}{=}r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right), \label{eq:prod_def} \end{equation}(15)where α2 = 2α0 − α1 (Eq. (5)). Owing to the symmetry of α1 and α2 with respect to the y-axis when δ = 0, it is sufficient to compute Tλ(α1) for 0 < α1 < α0. To this end, we define Eλ(θ) to be the ratio of the effective area to the on-axis geometric area Eλ(θ)=A(λ,θ)A(0),\begin{equation} E_{\lambda}(\theta) = \frac{A_{\infty}(\lambda, \theta)}{A_{\infty}(0)}, \label{eq:EA_norm} \end{equation}(16)and then it is easy to demonstrate (see Appendix A) that the Tλ(α1) function, for α1 ∈ (0,α0), can be computed from the integral equation Tλ(α1)=α0α10π/2sint[ddθ(θEλ(θ))]θ=θ(t)dt,\begin{equation} T_{\lambda}\left(\alpha_1\right) = \frac{\alpha_0}{\alpha_1}\int_{0}^{\pi/2}\!\!\sin t \left[\frac{\mbox{d}}{\mbox{d}\theta}\left(\theta E_{\lambda}(\theta)\right)\right]_{\theta = \theta(t)}\!\!\!\!\!\!\!\mbox{d}t, \label{eq:refl_prod} \end{equation}(17)where the expression in the [ ] brackets, after computing the derivative, is to be evaluated at θ(t)=(α0α1)sint,\begin{equation} \theta(t) = \left(\alpha_0-\alpha_1\right)\sin t, \label{eq:tdef} \end{equation}(18)where 0 ≤ t ≤ π/2 is a dummy integration variable. The condition 0 < α1 < α0 then implies that 0 < θ < α0: in other words, the Eλ(θ) function has to be known in the entire interval (0,α0). By computing the integral in Eq. (17) for any normalized effective area function Eλ(θ), one obtains the corresponding reflectivity product at λ (Eq. (15)), for α1 taking on values in the same interval.

As a first example, we put in adimensional form the geometric vignetting for a source at infinity, Eq. (12), E(θ)=12θπα0,\begin{equation} E(\theta) =1- \frac{2\theta}{\pi\alpha_0}, \label{eq:refl_prod_geom0} \end{equation}(19)and substitute this into Eq. (17). We obtain Tλ(α1)=α0α10π/2(sint4(α0α1)πα0sin2t)dt,\begin{equation} T_{\lambda}\left(\alpha_1\right) = \frac{\alpha_0}{\alpha_1}\int_{0}^{\pi/2}\!\left(\sin t - \frac{4\left(\alpha_0-\alpha_1\right)}{\pi\alpha_0}\sin^2t \right)\,\mbox{d}t, \label{eq:refl_prod_geom1} \end{equation}(20)which can be immediately solved, yielding Tλ(α1)=1.\begin{equation} T_{\lambda}\left(\alpha_1\right) = 1. \label{eq:refl_prod_geom2} \end{equation}(21)This means that rλ(α) = 1 for any incidence angle and wavelength, as expected.

As a second example, we show the computation of a reflectivity product from the effective area of a mirror shell with a multilayer coating. The product of the two reflectivities at 30 keV, directly computed from the multilayer recipe, is shown in Fig. 2, as a function of α1. The effective area of the mirror shell at the same energy (Fig. 3) was computed in (0,α0), using Eq. (14). Finally, we used Eq. (17) to re-derive the Tλ(α1) function from the effective area curve. The result of the inverse computation (Fig. 2, dots) closely matches the original reflectance product (Fig. 2, line).

thumbnail Fig. 2

The product of the reflectivities of the two segments of a Wolter-I shell with α0 = 0.212 deg, at 30 keV. The reflective surface is a multilayer coating, a Pt/C graded stack with 200 couples of layers. The d-spacing in the multilayer follows the power-law model d(k) = a   (b + k) − c (Joensen et al. 1995), with a = 77.4 Å, b =  −0.94, c = 0.223, a Pt thickness ratio Γ = 0.42, and a 4 Å roughness rms. The solid line is directly computed from the multilayer structure, whilst the dots are the result of the inverse computation from the effective area in Fig. 3.

thumbnail Fig. 3

The effective area at 30 keV (λ = 0.4 Å) of a Wolter-I mirror shell with R0 = 148.5 mm, L = 300 mm, f = 10 m, α0 = 0.212 deg, and a reflectivity product as shown in Fig. 2 (solid line). The effective area was used to re-derive the Tλ(α1) function (Fig. 2, dots).

3. Obstructions in nested mirror modules

We now deal with a quantification of the obstructions that reduce the off-axis effective area of a mirror shell, when nested in a mirror module. Even if this effect is essentially a geometric shadowing, the effective area reduction depends on λ, because the incidence angles, and consequently the mirror reflectivity rλ(α), vary over the mirror surface when the source is off-axis. This is especially true when multilayer coatings, which exhibit a complicated rλ(α) function, are adopted. A reduction of the effective area is then relevant for the X-ray wavelengths that were reflected at the obstructed regions.

thumbnail Fig. 4

Obstruction in an assembly of mirror shells. Rays impinging a mirror shell can be blocked by the adjacent shell with smaller diameter in only three ways, as listed in the text. The impact points are highlighted. Other mirror shells are not shown. The obscured regions of the primary mirror are indicated in gray: angles and the mirror spacings are greatly exaggerated.

The clear aperture of a mirror shell can be obstructed by several factors: the dense packing of shells, the structures for mechanical support of mirrors, and the presence of pre-collimators designed to reduce the stray light. However, we hereafter limit the discussion to the first point, i.e., the mutual obstruction of nested mirror shells, assuming that they are all co-axial, co-focal, and all have the same intersection plane at z = 0 (Fig. 4).

3.1. Obstruction parameters

We now draw our attention to a particular mirror shell in the mirror module (Fig. 4), and adopt the same notation presented in Sect. 2. The obstruction of this shell (which we refer to as a reflective shell) can be assumed to be caused solely by the shadow cast by the adjacent shell with smaller radius, which we refer to as blocking shell. The radius of its outer (i.e., non-reflective) surface at z = 0 is denoted with R0\hbox{$R_0^*$}, which is forcedly smaller than R0. We then define RM\hbox{$R_{\mathrm M}^*$}, Rm\hbox{$R_{\mathrm m}^*$}, L1\hbox{$L_1^*$}, and L2\hbox{$L_2^*$} to be respectively the maximum radius, the minimum radius, the primary segment length, and the secondary length of the outer surface of the blocking shell. We admit that, in general, L1L1\hbox{$L_1 \neq L_1^*$} and L2L2\hbox{$L_2 \neq L_2^*$}. If L1=L2\hbox{$L_1^* = L_2^*$}, we denote their common value with L.

thumbnail Fig. 5

Geometrical meaning of the obstruction parameters, Φ, Ψ, and Σ, for a pair of Wolter-I nested shells.

A simple geometrical construction shows (Fig. 4) that an obstruction of the reflective shell can occur in three ways:

  • 1.

    If rays intersect the blocking shell’s primary segment beforeimpinging the reflective shell;

  • 2.

    If, after the first reflection, rays impact onto the outer surface of the blocking shell’s primary segment;

  • 3.

    If, after the second reflection, rays are blocked by the outer surface of the blocking shell’s secondary segment.

The obstruction amount depends on the closeness of the two mirrors. This is often expressed along with FF, the filling factor FF=RMR0·\begin{equation} FF = \frac{R_{\mathrm M}^*}{R_0}\cdot \label{eq:ff} \end{equation}(22)The configuration with FF = 1 causes each shell to exactly fit the clear area section of the adjacent shell, hence it maximizes the effective area on-axis at the expense of the off-axis area. If FF > 1, the mirror assembly is self-obstructed even on-axis. For this reason and to enlarge the field of view of the optics, solutions with FF < 1 are in general envisaged.

Nevertheless, it is convenient to adopt other parameters in the estimation of obstructions. The first of these is Φ=R0RML1+α0,\begin{equation} \Phi = \frac{R_0-R_{\mathrm M}^*}{L_1^*}+\alpha_0, \label{eq:phi} \end{equation}(23)which represents the maximum angle visible from the reflective shell at the intersection plane (see Fig. 5), through the entrance pupil. Clearly, Φ > α0 whenever FF < 1, and Φ < α0 if FF > 1. As we later see, the Φ parameter drives the obstruction of the first kind, i.e., at the entrance pupil.

The obstruction of the second kind, i.e., at the intersection plane, is chiefly determined by another parameter, Ψ=R0R0L1,\begin{equation} \Psi = \frac{R_0-R_0^*}{L_1}, \label{eq:psi} \end{equation}(24)which denotes the clear angular aperture at the intersection plane, as seen from the maximum diameter of the reflective mirror shell (see Fig. 5). The third obstruction parameter, which drives the third species of obstruction, is Σ=R0RmL23α0,\begin{equation} \Sigma = \frac{R_0-R_{\mathrm m}^*}{L_2^*}-3\alpha_0, \label{eq:sigma} \end{equation}(25)which represents the angular aperture visible from the reflective shell at the intersection plane (see Fig. 5), through the exit pupil. The importance of these angles becomes clearer in Sect. 3.3, when we derive the general formula for the obstructed mirror effective area (Eq. (38)), which depends on Φ, Ψ, and Σ as parameters.

The obstruction parameters can be related to each other Φ=L1L1Ψ+(α0α0),Σ=L1L2Ψ3(α0α0),\begin{eqnarray} \label{eq:rel1}\Phi &=& \frac{L_1}{L_1^*}\Psi + \left(\alpha_0 -\alpha^*_0\right), \\[1.5mm] \label{eq:rel2}\Sigma &=& \frac{L_1}{L_2^*}\Psi -3\left(\alpha_0-\alpha^*_0\right), \end{eqnarray}where α0\hbox{$\alpha_0^*$} is the on-axis incidence angle on the blocking shell. Using Eq. (2), these relations can also be written as Φ=L1Ψ(1L1+14f),Σ=L1Ψ(1L234f)·\begin{eqnarray} \label{eq:rel1a}\Phi &=& L_1\Psi\left(\frac{1}{L_1^*}+\frac{1}{4f}\right), \\[1.5mm] \label{eq:rel2a}\Sigma &=& L_1\Psi\left(\frac{1}{L_2^*} -\frac{3}{4f}\right)\cdot \end{eqnarray}From Eqs. (28) and (29) it follows, in particular, that

  • if L1L1L2\hbox{$L_1^* \leq L_1 \leq L_2^*$}, then Σ < Ψ < Φ;

  • if L1=L1=L2\hbox{$L_1^* = L_1 = L_2^*$} and f ≫ L1, then Σ ≈ Ψ ≈ Φ.

3.2. Vignetting coefficients

We now provide analytical expressions for the obstructions caused by the mirror nesting. To this end, we introduce a vignetting coefficient, Vn(ϕ), with n = 1,2,3, for every kind of obstruction as listed in Sect. 3.1. For an infinitesimal mirror sector located at the polar angle ϕ, the nth vignetting coefficient is defined as the fraction of primary segment’s length left clear by the nth obstruction. This definition is analogous to that of V(ϕ), the self-vignetting factor for double reflection (Eq. (6)), which we already treated in detail (Spiga et al. 2009). We then have, in addition to V(ϕ), three vignetting coefficients:

  • V1(ϕ), due to the shadow cast by the blocking shell’s primary segment before the first reflection;

  • V2(ϕ), due to the shadow cast by the blocking shell’s primary segment after the first reflection;

  • V3(ϕ), due to the shadow cast by the blocking shell’s secondary segment after the second reflection.

Firstly, we consider all vignetting factors to be independent of each other: the respective obstructed (i.e., lost) fractions of primary mirror length are denoted with Qn = 1 − Vn.

thumbnail Fig. 6

Vignetting coefficients as a function of the polar angle for L2 = 310 mm, L1 = 300 mm, L2\hbox{$L_2^*$} = 290 mm, L1\hbox{$L_1^*$} = 280 mm, α0 = 0.3 deg, δ = 0 deg, θ = 0.25 deg, Φ = 0.43 deg, Ψ = 0.4 deg, and Σ = 0.408 deg. We also plot V, the vignetting factor for double reflection (Eq. (6)).

thumbnail Fig. 7

Vignetting coefficients as a function of the polar angle for the same configuration as in Fig. 6, but this time the source is at a finite distance (D = 127 m, δ = 0.095 deg). We note that the obstruction at the exit pupil is less severe than for a source at infinity.

As we anticipated in the previous section, it is the Φ, Ψ, and Σ parameters that primarily determine the vignetting coefficients. More precisely, they are functions of the incidence angles α1 and α2, and they can be computed via the following formulae V1(α1)=1+L1(Φα1)L1α1,V2(α1)=Ψα1,V3(α1)=1+L2(Σα2)L1α1,\begin{eqnarray} \label{eq:V1_}V_1\left(\alpha_1\right) &=& 1+ \frac{L_1^*\left(\Phi-\alpha_1\right)}{L_1\alpha_1}, \\[1.5mm] \label{eq:V2_}V_2\left(\alpha_1\right) &=& \frac{\Psi}{\alpha_1},\\[1.5mm] \label{eq:V3_} V_3\left(\alpha_1\right) &=& 1+\frac{L_2^*\left(\Sigma-\alpha_2\right)}{L_1\alpha_1}, \end{eqnarray}where α1 + α2 = 2α0, and with the usual constraint 0 ≤ Vn ≤    1 for each n, otherwise we set Vn to 0 or 1, respectively. The derivation of Eqs. (30) to (32) is not difficult, but quite lengthy, so it has been postponed to Appendix B.

The dependence of the coefficients on ϕ is obtained by substituting the expressions of α1 and α2 (Eqs. (3) and (4)). We note that the obstructed region of the primary segment in the first and third kind of vignetting is located near the intersection plane, whereas V2 results from an obscuration of the primary mirror near z =  + L1 (see Fig. 4), in a precisely similar way to V.

Example plots of V1, V2, V3, and V, as functions of ϕ, are drawn in Figs. 6 and 7. The adopted values correspond to the case of two mirror shells with R0 = 210 mm, R0=207.9\hbox{$R_0^*=207.9$} mm, and f = 10 m: the other parameter values are reported in the figure caption, including the lengths of the mirrors, which have been chosen to fulfill the relations L1<L2<L1<L2\hbox{$L_1^* < L_2^* < L_1 <L_2$}. As we later see (Sect. 4.2), this choice is not accidental: it represents a compromise to minimize all obstructions.

Equations (30) to (32) and Figs. 6, and 7 show that:

  • if the conditionsα1 < Φ, α1 < Ψ, and α2 < Σ are fulfilled at a given ϕ, there is no obstruction at that polar angle;

  • if α1 > Ψ, but V < V2 as in Figs. 6 and 7, the obstruction at the intersection plane is ineffective, because all rays that were blocked would have missed the secondary segment;

  • the obstructions related to V1 and V2 are maximum at polar angles close to π;

  • the effect of V3 is larger at ϕ ≈ 0, where α1 is shallower and α2 is larger;

  • with the source at infinity, the obstruction related to V3 can be very large, especially near ϕ ≈ 0. In contrast, the effect of V3 is mitigated if the source is at a finite distance, because α2 becomes smaller;

  • if L1L1\hbox{$L_1 \geq L_1^*$}, then V1 > V2 for all ϕ;

  • if L1=L1\hbox{$L_1 = L_1^*$} and f ≫ L1, then Φ ≈ Ψ (Eq. (28)) and also V1 ≈ V2 for all ϕ.

3.3. General formula for the effective area of obstructed mirror shells

We already pointed out that the obstructions described by V and V2 are located close to the maximum diameter of a Wolter-I mirror shell, whereas the ones related to V1 and V3 concern mainly the region close to its intersection plane. Therefore, V and V2 are in competition for the obstruction near z =  + L1, while V1 and V3 do the same near z = 0 (Fig. 4). The total obstruction at the generic polar angle ϕ is then Qtot=max(Q,Q2,0)+max(Q1,Q3,0),\begin{equation} Q_{\mathrm{tot}} = \max\left(Q, Q_2, 0\right) + \max\left(Q_1, Q_3, 0\right), \label{eq:Q_total} \end{equation}(33)where the “0” value has been added to avoid negative obstructions. The corresponding total vignetting factor is 1 − Qtot, i.e., Vtot=min(V,V2,1)max(1V1,1V3,0).\begin{eqnarray} V_{\mathrm{tot}} = \min\left(V, V_2, 1\right) - \max\left(1-V_1, 1-V_3, 0\right). \label{eq:V_total} \end{eqnarray}(34)Using the expressions of V (Eq. (6)) and V2 (Eq. (31)), we can write the first term of Eq. (34) as min(V,V2,1)=min(L2α2L1α1,Ψα1,1).\begin{equation} \min\left(V, V_2, 1\right) = \min\left(\frac{L_2\alpha_2}{L_1\alpha_1}, \frac{\Psi}{\alpha_1}, 1\right). \label{eq:V_total1} \end{equation}(35)Using Eqs. (30) and (32), the second term reads max(1V1,1V3,0)=max(L1(α1Φ)L1α1,L2(α2Σ)L1α1,0).\begin{equation} \max\left(1-V_1, 1-V_3, 0\right) = \max\left(\frac{L_1^* \left(\alpha_1-\Phi\right)}{L_1\alpha_1}, \frac{L_2^* \left(\alpha_2-\Sigma \right)}{L_1\alpha_1},0 \right). \label{eq:V_total2} \end{equation}(36)If positive and not larger than 1, the total vignetting factor (Eq. (34)), can be used to compute the obstructed mirror effective area. An infinitesimal sector of the primary mirror with polar aperture Δϕ, as seen from the off-axis source, has a vignetted geometric area R0VtotL1α1Δϕ, hence the total effective area is AD(λ,θ)=2R00πVtotL1α1rλ(α1)rλ(α2)dϕ.\begin{equation} A_D(\lambda, \theta) = 2R_0 \int_0^{\pi}\!\!V_{\mathrm{tot}}L_1\alpha_1 \,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda} \left(\alpha_2\right) \,\mbox{d}\varphi. \label{eq:area_total} \end{equation}(37)By substituting Eq. (34) into Eq. (37), it is straightforward to derive the final expression for the off-axis effective area of the obstructed mirror shell AD(λ,θ)=2R00π[()minQmax]0rλ(α1)rλ(α2)dϕ,\begin{equation} A_D(\lambda, \theta) = 2R_0 \int_0^{\pi}\!\!\left[(L\alpha)_{\mathrm{min}}- Q_{\mathrm{max}}\right]_{\ge0} \,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right) \,\mbox{d}\varphi, \label{eq:area_total_fin} \end{equation}(38) where ()min=min(L1α1,L2α2,L1Ψ),Qmax=max[L1(α1Φ),L2(α2Σ),0],\begin{eqnarray} \label{eq:Lamin} (L\alpha)_{\mathrm{min}} & = & \min\left(L_1\alpha_1, L_2\alpha_2, L_1\Psi\right),\\[1.5mm] \label{eq:Qmax} Q_{\mathrm{max}} & = & \max\left[L_1^*\left(\alpha_1-\Phi\right), L_2^*\left(\alpha_2-\Sigma\right),0 \right], \end{eqnarray}provided that α1 ≥ 0, α2 ≥ 0, as usual. The  [    ]  ≥ 0 brackets in Eq. (38) mean that the enclosed expression ()min − Qmax, if negative, must be set to zero. For this reason, the integration cannot be carried out independently for the two terms. Moreover, because of the presence of Ψ in Eq. (39) and the different values of Φ and Σ in Eq. (40), the computation is asymmetric with respect to the y-axis, also when δ = 0. Therefore, the integral in Eq. (38) is not equivalent to twice the same integral over  [0/2 ] , unless all these conditions are fulfilled (as in Sect. 4.1): i) δ = 0, ii) Φ ≃ Σ, and iii) there is no vignetting of the second kind.

3.4. Conditions for obstruction-free mirror shells

We now discuss the conditions for the mirror shell to be obstruction-free. It can be seen, from Eqs. (39) and (40), that there is no obstruction if, and only if, for all ϕ{α1<Φ,α2<Σ,min(L1α1,L2α2)<L1Ψ.\begin{equation} \left\{\begin{array}{l} \alpha_1 < \Phi, \\ \alpha_2 < \Sigma, \\ \min\left(L_1\alpha_1, L_2\alpha_2\right) < L_1\Psi. \end{array}\right. \label{eq:obsfree} \end{equation}(41)If these inequalities are fulfilled, Eq. (38) correctly reduces to Eq. (7). The first two conditions are simply met if the maximum values of α1 and α2 do not exceed Φ and Σ, respectively. The third condition requires that either α1 < Ψ or α2<L1L2Ψ\hbox{$\alpha_2 < \frac{L_1}{L_2}\Psi$}, which in turn are equivalent1 to α1+α2<(1+L1L2)Ψ,\begin{equation} \alpha_1 +\alpha_2 < \left(1+\frac{L_1}{L_2}\right)\Psi, \label{eq:noV2_cond} \end{equation}(42)i.e., using Eq. (5), 2α0L2L1+L2<Ψ,\begin{equation} \frac{2\alpha_0L_2}{L_1+L_2} < \Psi, \label{eq:noV2_cond2} \end{equation}(43)which reduces to α0 < Ψ if L1 = L2. The definition of Ψ then allows us to write Eq. (43) as α02L1L2L1+L2<R0R0.\begin{equation} \alpha_0\frac{2L_1L_2}{L_1+L_2} <R_0-R_0^*. \label{eq:noV2_cond3} \end{equation}(44)This result has an immediate geometric interpretation, by noting that left-hand of Eq. (44) is exactly the maximum possible distance from the reflective shell of a ray reflected twice (Fig. 8). All other rays undergoing a double reflection cannot exceed this distance, which therefore represents the minimum separation for the two shells at z = 0.

thumbnail Fig. 8

There is no vignetting at the intersection plane if the distance of the ray reflected at the maximum and the minimum diameter is smaller than the spacing of the two mirror shells.

Finally, we define ˜L\hbox{$\tilde{L}$} to be the equivalent length of the mirror 2˜L=1L1+1L2,\begin{eqnarray} \frac{2}{\tilde{L}} = \frac{1}{L_1}+\frac{1}{L_2}, \label{eq:eqlen} \end{eqnarray}(45)and using Eqs. (3) and (4), we can express the obstruction avoidance conditions as α0+δ+θ<Φ,α0δ+θ<Σ,α0˜L<R0R0.\begin{eqnarray} \label{eq:noobs1}& &\alpha_0+\delta +\theta < \Phi,\\ \label{eq:noobs2}& &\alpha_0-\delta +\theta < \Sigma,\\ \label{eq:noobs3}& &\alpha_0\tilde{L} < R_0-R_0^*. \end{eqnarray}We note that if L1 = L2 = L, then also ˜L=L\hbox{$\tilde{L} = L$}, and that the last condition solely depends on the mirror pair geometry, not on the off-axis angle.

4. Some applications

4.1. Analytical expressions for the geometric area

As a first application, we derive some expressions for the geometric area (i.e., in the ideal case rλ(α) = 1 for all α) of an obstructed mirror as a function of θ. For simplicity, we only consider the case that δ = 0 and L1=L2=L1=L2\hbox{$L_1 = L_2 = L_1^* = L_2^*$}, and we suppose the mirror to be unobstructed on-axis, i.e., that Φ ≥ α0 (Eq. (23)). Finally, we reasonably assume that f ≫ L, so that Φ ≈ Ψ ≈ Σ (Sect. 3.1): consequently, Ψ ≥ α0 and there is no vignetting at the intersection plane (as in Fig. 8). Equation (39) then becomes, as in the unobstructed case, ()min=L·min(α1,α2).\begin{equation} (L\alpha)_{\mathrm{min}} = L\,\cdot \min\left(\alpha_1, \alpha_2\right). \label{eq:La_geom} \end{equation}(49)Adopting Φ as a unique obstruction parameter, and since min(α1,α2) + max(α1,α2) = 2α0 (Eq. (5)), Eq. (40) turns into QmaxL[2α0min(α1,α2)Φ],\begin{eqnarray} Q_{\mathrm{max}} \simeq L\left[2\alpha_0-\min\left(\alpha_1, \alpha_2\right)-\Phi\right], \label{eq:Qmax_geom} \end{eqnarray}(50)on the condition that both terms are non-negative. Substituting this expression into Eq. (38) and using Eq. (10) with δ = 0, we obtain AD(θ)=2R0L0π[αminβmax]0dϕ,\begin{equation} A_D(\theta)=2R_0 L \int_0^{\pi}\!\!\left[\alpha_{\mathrm{min}}- \beta_{\mathrm{max}}\right]_{\ge0}\,\mbox{d}\varphi, \label{eq:area_geom} \end{equation}(51)where we defined αmin=[α0θ|cosϕ|]0,βmax=[α0+θ|cosϕ|Φ]0,\begin{eqnarray} \label{eq:alphamin} \alpha_{\mathrm{min}} &=& \left[\alpha_0-\theta\,|\cos\varphi|\right]_{\ge0},\\[1.5mm] \label{eq:betamax} \beta_{\mathrm{max}} &=& \left[\alpha_0+\theta\,|\cos\varphi|-\Phi\right]_{\ge0}, \end{eqnarray}and where the  [    ]  ≥ 0 brackets have the same meaning as those appearing in Eq. (38).

We firstly consider the case α0 < Φ < 2α0, which implies that Φ − α0 < α0. As long as θ < Φ − α0, αmin > 0 but βmax = 0 for all ϕ, i.e., the mirror is not obstructed and Eq. (12) remains valid.

We now increase θ. Since Φ − α0 < Φ/2 < α0 by hypothesis, we are allowed to consider the case Φ − α0 < θ < Φ/2. For all ϕ, we still have αmin > 0, but this time βmax is positive for θ    | cosϕ |  > Φ − α0: for these values of ϕ, the mirror starts to be obstructed. Then Eq. (51) changes into A(θ)=4R0LarccosΦα0θπ/2(α0θcosϕ)dϕ+4R0L0arccosΦα0θ(Φ2θcosϕ)dϕ,\begin{eqnarray} A_{\infty}(\theta) &=& 4R_0 L \int_{\arccos\frac{\Phi-\alpha_0}{\theta}}^{\pi /2}\!\!\!\left(\alpha_0-\theta\cos\varphi\right)\,\mbox{d}\varphi\,\nonumber \\[1.5mm] \label{eq:obs_geom} && + 4R_0 L\int_0^{\arccos\frac{\Phi-\alpha_0}{\theta}}\!\!\!\! (\Phi-2\theta\cos\varphi)\,\mbox{d}\varphi, \end{eqnarray}(54)where the two terms are the unobstructed and the obstructed part of the area, respectively. We note that the two integrands are positive in the respective integration ranges, hence the integration can be carried out and we obtain A(θ)A(0)=12θπα0[1+S(Φα0θ)],\begin{equation} \frac{A_{\infty}(\theta)}{A_{\infty}(0)} = 1-\frac{2\theta}{\pi\alpha_0}\left[1+\, S\!\left(\frac{\Phi-\alpha_0}{\theta}\right)\right], \label{eq:obs_case1} \end{equation}(55)where we normalized the result to A(0), the on-axis geometric area (Eq. (11)), and defined the non-negative S(x) function, S(x)=1x2xarccosx,with0x1.\begin{equation} S(x) = \sqrt{1-x^2}-x\arccos x, \,\, \mbox{with} \,\, 0 \le x \le 1. \label{eq:Sfunction} \end{equation}(56)We recognize in the two first terms of Eq. (55) the usual unobstructed geometric area trend (Eq. (12)), whereas the term proportional to the S function represents the area lost because of the obstruction. Since S(1) = 0, Eq. (55) converges to the unobstructed area trend at θ = Φ − α0.

thumbnail Fig. 9

Normalized geometric area, A(θ)/A(0) of an obstructed mirror shell with α0 = 0.2 deg, δ = 0, and L1 = L2, as a function of the off-axis angle, for different values of the obstruction parameter Φ, with Φ ≈ Ψ ≈ Σ and L = L ∗ . The curves are traced using the analytical formulae reported in Sect. 4.1. The curve for the unobstructed mirror (the solid line) is also valid for any obstructed mirror with Φ > 2α0.

We now consider the case θ > Φ/2, still for Φ < 2α0 and δ = 0. The first term of Eq. (54) remains unchanged, even if the integration range shrinks because of the larger values of θ. In the second integral, indeed, we have to replace the lower integration limit with arccosΦ2θ\hbox{$\arccos\!\frac{\Phi}{2\theta}$} to keep the integrand non-negative. This integration limit also guarantees that αmin > 0 even for θ > α0. Because of this change, the resulting expression A(θ)A(0)=12θπα0[1+S(Φα0θ)2S(Φ2θ)],\begin{equation} \frac{A_{\infty}(\theta)}{A_{\infty}(0)} = 1-\frac{2\theta}{\pi\alpha_0}\left[1+ \,S\!\left(\frac{\Phi-\alpha_0}{\theta}\right)- 2\, S\!\left(\frac{\Phi}{2\theta}\right)\right], \label{eq:obs_case2} \end{equation}(57)has an additional term with respect to Eq. (55).

Finally, we assume that Φ ≥ 2α0. This time α0 ≤ Φ − α0, hence the geometric area decrease deviates from linearity (for θ > α0) before the mirror begins to be obstructed (θ > Φ − α0). The condition αmin > 0 then implies that the lower integration limit in the first term of Eq. (54) has to be replaced with arccosα0θ\hbox{$\arccos\frac{\alpha_0}{\theta}$}: consequently, αmin = 0 for 0<ϕ<arccosα0θ\hbox{$0<\varphi<\arccos\frac{\alpha_0}{\theta}$}, and the obstructed term is zero. We conclude that for Φ > 2α0 and δ = 0 the geometric area trend equals the unobstructed one, i.e., Eq. (12) for θ < α0, and Eq. (13) for θ > α0. As expected, in the limit Φ = 2α0 Eq. (57) becomes identical to Eq. (13).

As an example, we trace in Fig. 9 some geometric area curves of a mirror shell for different values of the obstruction parameter Φ ≥ α0, using Eqs. (12), (55), and (57) in the respective intervals of validity. We also plotted the unobstructed geometric area of the mirror (Eqs. (12) and (13)), which is also valid for Φ > 2α0. Some curves are also validated in Sect. 5 by means of an accurate ray-tracing computation.

4.2. Design of the mirror module

The obtained results can also be applied to the problem of designing a mirror module. In general, the radius and the length of the outermost shell of the module is assigned on the basis of the allocable space for the optics payload. Starting from this one, shells with decreasing radii are added, leaving a sufficient spacing to minimize the mirror vignetting for off-axis angles within the field of view of the optic. It is then convenient to reduce the obstruction, by increasing Φ and Σ (Sect. 3.1) for every pair of consecutive shells. For any choice of R0, R0\hbox{$R_0^*$}, and L1, this can be obtained by designing each mirror pair with L1<L1\hbox{$L_1^* < L_1$} (Eq. (26)) and L2<L1\hbox{$L_2^*< L_1$} (Eq. (27)). On the other side, mitigation of the mirror self-vignetting for double reflection (Eq. (6)) requires that L2 > L1 and L2>L1\hbox{$L_2^* >L_1^*$}. Hence, if we label the mirror shells with k = 1, 2, ...from large to small diameters, a performing module design might consist of mirror shells with ···L1,k+1L2,k+1L1,kL2,k···,\begin{equation} \cdots \leq L_{1,k+1} \leq L_{2,k+1} \leq L_{1,k} \leq L_{2,k} \leq \cdots, \label{eq:lengths} \end{equation}(58)i.e., decreasing lengths as the diameter is reduced, and with secondary segments longer than the respective primary but shorter than the primary of the adjacent mirror shell with larger diameter. A module design with mirror lengths scaled to the diameter has already been studied by Conconi et al. (2010) to minimize the defocusing due to the field curvature in the WFXT telescope.

Therefore, after choosing a sequence of mirror lengths according Eq. (58), we apply the obstruction-free conditions (Sect. 3.4). Using Eqs. (24) and (29), we rewrite Eq. (47) as α0,k+θ(R0,kR0,k+1τk+1)(1L2,k+134f),\begin{equation} \alpha_{0,k}+\theta \leq \left(R_{0,k}-R_{0,k+1}-\tau_{k+1}\right)\left(\frac{1}{L_{2,k+1}}-\frac{3}{4f}\right), \label{eq:cond1} \end{equation}(59)where τk + 1 is the thickness of the (k + 1)th shell, in general chosen to be proportional to the radius to maintain a constant mirror stiffness throughout the module. If we now define 1/L2f,k + 1 = 1/L2,k + 1 − 3/4f, we can write the last equation as R0,k+1+τk+1R0,k(α0,k+θ)L2f,k+1.\begin{equation} R_{0, k+1} +\tau_{k+1} \leq R_{0, k}- \left(\alpha_{0,k}+\theta\right)L_{2f, k+1}. \label{eq:des1} \end{equation}(60)In reality, Eq. (58) implies that Σ < Φ for all pairs of mirrors, hence Eq. (47) also fulfills Eq. (46). The last relation to be satisfied is Eq. (48), which simply reads R0,k+1+τk+1R0,kα0,k˜Lk+1.\begin{equation} R_{0, k+1} +\tau_{k+1} \leq R_{0, k}- \alpha_{0,k}\tilde{L}_{ k+1}. \label{eq:des2} \end{equation}(61)We then derive from Eqs. (60) and (61) the condition to be fulfilled by the kth couple of shells, in order to avoid obstructions up to an off-axis angle θR0,k+1+τk+1R0,kmax[α0,k˜Lk+1,(α0,k+θ)L2f,k+1].\begin{equation} R_{0, k+1} +\tau_{k+1} \leq R_{0, k}- \max\left[\alpha_{0,k}\tilde{L}_{ k+1}, \left(\alpha_{0,k}+\theta\right) L_{2f, k+1}\right]. \label{eq:des} \end{equation}(62)If mirror lengths are chosen according to Eq. (58), the last relation enables the computation of the maximum possible value of R0,k + 1 from R0,k, L1,k + 1, L2,k + 1, and the relation τk = τ(R0,k). When applied recursively from the outermost radius inwards, it provides us with the optimal population of mirror shells in the optical module.

thumbnail Fig. 10

Initial positions and destinations of 40 000 rays at the entrance pupil (points) for a Wolter-I mirror shell with L = L ∗ , and the angular parameters reported in the legend (Φ ≈ Ψ ≈ Σ). The radial scale has been expanded. Only rays that would have struck the primary mirror were traced. The limits of the regions of different vignetting (dashes) are computed from the vignetting coefficients.

5. Validation with ray-tracing results

We hereafter validate the formulae found in the previous sections by comparing them with the results of a detailed ray-tracing. As a first example, we validate the expressions of the vignetting coefficients (Sect. 3.2). In Fig. 10, we show the entrance section of an obstructed Wolter-I mirror shell with R0 = 139.6 mm, α0 = 0.2 deg, f = 10 m, and L1 = L2 = 300 mm. Rays coming from a source at infinite distance, off-axis by an angle θ = 0.15 deg, have been traced by simulating their reflection on the mirror: the positions at the entrance pupil of rays that impinged the primary segment are drawn, disregarding the others. The points then fill the geometric section of the primary, as viewed from the direction of the source. The blocking shell has the same focal and length of the reflective one, but a different radius R0=138.3\hbox{$R_0^* = 138.3$} mm.

Depending on their initial coordinates, the traced rays can obstructed in various ways, highlighted with different colors in Fig. 10, or reach the focal plane if they fall in the red region of the mirror aperture. We note that, in agreement with the discussion in Sect. 3.2, the obstructions at the entrance (yellow) and the exit pupil (blue) are mainly located close to the intersection plane, i.e., the inner contour of the colored area, while the obstruction after the first reflection (orange), and the vignetting for double reflection (green), mainly occur for rays firstly reflected at locations far from the intersection plane. Moreover, in this case (Ψ ≃ Φ > α0), the obstruction after the first reflection is ineffective because the orange region, which encloses the rays blocked after the first reflection, is completely surrounded by the green region, which comprises the rays that missed the second reflection. This is in more than qualitative agreement with the analytical tractation: in the same figure we have also traced the analytical expressions of the vignetting coefficients (dashed lines), after translating them into expressions of the radial coordinate (see Appendix B) and projecting them at the mirror’s entrance plane in the direction of the incident rays. We note that the lines perfectly follow the boundaries of the vignetted regions: the agreement shows that the expressions of the vignetting coefficients (Eqs. (30) through (32)) are correct.

thumbnail Fig. 11

Comparison between some analytical curves of Fig. 9 (lines), and the results of an accurate ray-tracing (symbols) for an obstructed Wolter-I mirror with the same α0 = 0.2 deg, R0 = 139.6 mm, L = L ∗  = 300 mm. The outer radius of blocking shell, R0\hbox{$R_0^*$}, takes on the values 137.24, 138.29, and 138.55 mm, for Φ to match the values used to draw the analytical curves (but the unobstructed one). The accord between the formulae and the tracing is excellent within the statistical error.

As a second example, we compare the analytical expressions of the geometric area found in Sect. 4.1 with the findings of the ray-tracing. Figure 11 reports some of the normalized geometric area curves (lines) already shown in Fig. 9, i.e., for no obstruction and for two values of Φ in the interval  [α0,2α0). We have also plotted as symbols the geometric area values obtained by applying a ray-tracing routine for a Wolter-I, reflective mirror shell and a co-focal blocking shell with the same length but variable values of R0<R0\hbox{$R_0^* < R_0$}. Two values of R0\hbox{$R_0^*$} (138.55 and 138.29 mm) were chosen to match the two values of Φ (0.2 and 0.25 deg) we used to trace the analytical curves, while the smallest one (137.24 mm) was selected to return Φ = 0.45 deg, i.e., a value larger than 2α0.

Inspection of Fig. 11 shows that the findings of the two methods are in excellent agreement, within the error bars of the ray-tracing. On the other hand, the ray-tracing routine used is a quite complex code, and the relative computation required a few hours time to reach a few percent accuracy, whilst the analytical curves can be traced immediately and without being affected by statistical errors. We also note that the results of the ray-tracing for Φ = 0.45 deg are perfectly reproduced by the unobstructed mirror analytical curve, in agreement with the conclusion in Sect. 4.1 that there is no obstruction at any off-axis angle if Φ > 2α0. This happens because, with such a loose mirror nesting, θ must be larger than α0 for the mirror to start being obstructed. In these conditions, the obstruction is ineffective because it is completely superseded at all polar angles by the vignetting for double reflection.

Finally, we directly compare in Fig. 12 some obstructed effective area off-axis curves, as a function of the X-ray energy, as computed from a ray-tracing and using Eq. (38). The reflective mirror shell has a fixed radius and a graded multilayer coating (see e.g., Joensen 1995) to enhance its reflectivity in hard X-rays up to 50 keV and beyond. The characteristics of the reflective shell and the description of the multilayer stack are reported in the figure caption. The off-axis angle θ = 6 arcmin is the same for all curves, while the outer radius of the blocking shell, R0\hbox{$R_0^*$}, has been varied. As expected, the effective area decreases as the tightness of the nesting increases (i.e., as R0\hbox{$R_0^*$} approaches R0), and moreover the findings of the ray-tracing (symbols) are in perfect agreement with those of the analytical computation (lines).

thumbnail Fig. 12

Effective area of an obstructed mirror shell with R0 = 148.5 mm, f = 10 m, α0 = 0.212 deg, L = L = 300 mm, and the X-ray source at infinite distance, 6 arcmin off-axis. The curves are computed for three different values of R0\hbox{$R_0^*$}. The reflective shell has a Pt/C graded multilayer, whose d-spacing variation in depth is a power law, d(k) = a   (b + k) − c, where k = 1 ... 200, a = 115.5 Å, b = 0.9, c = 0.27. We moreover assumed a Pt thickness-to-d-spacing ratio Γ = 0.35, a roughness rms σ = 4 Å, and a stack that ends with Pt on top.

We note that the high-energy part of the effective area is less affected by obstruction effects, because it results from reflection at polar angles ϕ ≈  ± π/2, where max(α1,α2) takes on the smallest value, i.e., α0, and the reflectivity is higher: this is also the angular region where the obstruction is lower (see Fig. 10).

6. Conclusions

We have developed the analytical formalism for the off-axis effective area of Wolter-I mirror shells, in double cone approximation, which we began to describe in a previous paper (Spiga et al. 2009).

We have shown that the analytical expression of the effective area off-axis can be inverted to derive the product of the reflectivity of the two segments (Sect. 2.2). This might be useful to future developments for computing a suitable multilayer recipe to return the desired effective area trend off-axis.

We have found analytical expressions for the vignetting coefficients (Sect. 3.2), for the three possible sources of obstruction in nested mirror modules, as a function of the azimuthal coordinate of the mirror surface.

Using the vignetting coefficients, we have derived an integral formula (Eq. (38)) for the obstructed effective area of a Wolter-I X-ray mirror in double cone approximation, with any reflective coating, including multilayers. The computation only requires the standard routines for the reflectivity of the coating, and an integration over the azimuthal coordinate of the mirror shell.

We have obtained analytical expressions of the obstructed geometric area (Sect. 4.1) for the case of a source at infinite distance, and applied the result to the problem of designing an optical module that does not suffer from the mutual obstructions of mirrors (Sect. 4.2).

Finally, the results have been validated by means of a comparison with the findings of a detailed ray-tracing (Sect. 5).

As a final application, we note that each vignetting coefficient can be adapted to estimate the unwanted vignetting caused by collimators aimed at reducing the stray light in mirror modules (see, e.g., Cusumano et al. 2007). For example, V1 would quantify the vignetting of the baffle at the entrance pupil, if RM\hbox{$R^*_{\mathrm M}$} is interpreted as the outer radius of the collimator ring and L1\hbox{$L_1^*$} as its distance from the intersection plane. With analogous substitutions, V3 would represent the vignetting of a baffle located at the exit pupil, and V2 would express the vignetting of the baffle at the intersection plane, even though this kind of baffle can be designed to avoid any obstruction of focused rays (Sect. 3.4).


1

From Eq. (42), it might alternatively follow that α1<L1L2Ψ\hbox{$\alpha_1 < \frac{L_1}{L_2}\Psi$} or α2 < Ψ. However, if L1 < L2 this would still imply that α1 < Ψ, whereas if L1 > L2 we would still have α2<Ψ<L1L2Ψ\hbox{$\alpha_2 < \Psi < \frac{L_1}{L_2}\Psi$}. In both cases, Eq. (42) remains valid.

2

From the expression of the k1\appendix \setcounter{section}{2} \hbox{$\underline{k}_1$} vector, it is easy to derive the stray light pattern on a detector at a generic distance d, i.e., at z = -d. The final coordinates of rays generated close to the intersection plane are x(ϕ)=R0cosϕk1,xk1,zd,y(ϕ)=R0sinϕk1,yk1,zd,\appendix \setcounter{section}{2} \begin{eqnarray} x(\varphi) &=& R_0\cos\varphi-\frac{k_{1,x}}{k_{1,z}}d\nonumber,\\ y(\varphi) &=& R_0\sin\varphi-\frac{k_{1,y}}{k_{1,z}}d\nonumber, \end{eqnarray}and, substituting the expression of the k1\appendix \setcounter{section}{2} \hbox{$\underline{k}_1$} vector components, we obtain the parametric equation of the pattern x(ϕ)=[R0(2α0+δ)d]cosϕ+θdcos2ϕ,y(ϕ)=[R0(2α0+δ)d]sinϕ+θdsin2ϕ.\appendix \setcounter{section}{2} \begin{eqnarray} x(\varphi) &=& \left[R_0-(2\alpha_0+\delta)d\right]\cos\varphi+\theta d \cos2\varphi,\nonumber\\ y(\varphi) &=& \left[R_0-\left(2\alpha_0+\delta\right)d\right]\sin\varphi+\theta d \sin2\varphi.\nonumber \\[-1.03cm] &&\nonumber \end{eqnarray}

3

From Eq. (B.22) we obtain that all rays reflected twice at locations close to the intersection plane (r = R0) converge to a single point F = (− θf′,0, − f′), as expected, with f=R04α0δ,\appendix \setcounter{section}{2} \begin{displaymath} f' = \frac{R_0}{4\alpha_0-\delta}, \end{displaymath}which, recalling Eq. (2) and that δ = R0/D, becomes the usual conjugate points formula, 1f=1f+1D·\appendix \setcounter{section}{2} \begin{displaymath} \frac{1}{f} = \frac{1}{f'}+\frac{1}{D}\cdot \end{displaymath} Focusing at that point, indeed, does not exactly occur for points more distant from the intersection plane. This is a well-known limit of the double cone approximation, even if it does not affect our computation of the effective area within large limits (Sect. 2.1).

Acknowledgments

This research is funded by ASI (Italian Space Agency, contract I/069/09/0). V. Cotroneo (Harvard-Smithsonian CfA, Boston, USA) is acknowledged for useful discussions and paper editing.

References

  1. Basso, S., Pareschi, G., Citterio, O., et al. 2010, in Proc. SPIE, 7732, 773218 [Google Scholar]
  2. Bookbinder, J. A. 2010, in Proc. SPIE, 7732, 77321 [Google Scholar]
  3. Conconi, P., Campana, S., Tagliaferri, G., et al. 2010, MNRAS, 405, 877 [NASA ADS] [Google Scholar]
  4. Cotroneo, V., Pareschi, G., Spiga, D., et al. 2010, in Proc. SPIE, 7732, 77322 [Google Scholar]
  5. Cusumano, G., Artale, M. A., Mineo, T., et al. 2007, in Proc. SPIE, 6688, 66880 [Google Scholar]
  6. Hailey, J. C., An, H., Blaedel, K. L., et al. 2010, in Proc. SPIE, 7732, 77320 [Google Scholar]
  7. Joensen, K. D., Voutov, P., Szentgyorgyi, A., et al. 1995, Appl. Opt., 34, 7935 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  8. Kunieda, H., Awaki, H., Furuzawa, A., et al. 2010, in proc. SPIE, 7732, 773214 [Google Scholar]
  9. Madsen, K. K., Harrison, F. A., Mao, P. H., et al. 2009, in Proc. SPIE, 7437, 743716 [Google Scholar]
  10. Mangus, J. D., & Underwood, J. H. 1969, Appl. Opt., 8, 95 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  11. Mao, P. H., Harrison, F. A., Windt, D. L., et al. 1999, Appl. Opt., 38, 4766 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  12. Mao, P. H., Bellan, L. M., Harrison, F. A., et al. 2000, in Proc. SPIE, 4138, 126 [Google Scholar]
  13. Raimondi, L., & Spiga, D. 2010, in Proc. SPIE, 7732, 77322Q [Google Scholar]
  14. Spiga, D., & Cotroneo, V. 2010, in Proc. SPIE, 7732, 77322K [Google Scholar]
  15. Spiga, D., Cotroneo, V., Basso, S, et al. 2009, A&A, 505, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Stover, J. C., 1995, Optical Scattering: measurement and analysis, SPIE Optical Engineering Press [Google Scholar]
  17. Tawara, Y., Yamashita, K., Kunieda, H., et al. 1998, in Proc. SPIE, 3444, 569 [Google Scholar]
  18. Van Speybroeck, L. P., & Chase, R. C. 1972, Appl. Opt., 11, 440 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  19. Zhao, P., Jerius, D. H., Edgar, R. J., et al. 2004, in Proc. SPIE, 5165, 482 [Google Scholar]

Appendix A: Inversion of the effective area integral

We hereafter report the derivation of the inverse integral equation (Eq. (17)), which allows us to derive the reflectivity product from the off-axis effective area (Sect. 2.2). We start from the expression of the unobstructed effective area for L1 = L2, δ = 0, and 0 < θ ≤ α0 (Eq. (14)), which we rewrite in terms of the adimensional ratio Eλ(θ) (Eq. (16)), Eλ(θ)=2πα0α0θα0α1rλ(α1)rλ(α2)θ2(α0α1)2dα1,\appendix \setcounter{section}{1} \begin{eqnarray} E_{\lambda}(\theta) =\frac{2}{\pi \alpha_0} \int_{\alpha_0-\theta}^{\alpha_0}\!\frac{\alpha_1\,r_{\lambda}\left(\alpha_1\right)\,r_{\lambda}\left(\alpha_2\right)}{\sqrt{\theta^2-\left(\alpha_0-\alpha_1\right)^2}}\, \mbox{d}\alpha_1, \label{eq:Aeff_alpha1} \end{eqnarray}(A.1)where 0 < θ ≤ α0. By setting the positive variable ω = α0 − α1, the integral becomes Eλ(θ)=2πα00θ(α0ω)Tλ(ω)θ2ω2dω.\appendix \setcounter{section}{1} \begin{eqnarray} E_{\lambda}(\theta) = \frac{2}{\pi\alpha_0} \int^{\,\theta}_{0}\!\frac{\left(\alpha_0-\omega\right)\, T_{\lambda}(\omega)}{\sqrt{\theta^2-\omega^2}}\, \mbox{d}\omega. \label{eq:Aeff_omega} \end{eqnarray}(A.2)We now set Kλ(ω) = ω2(α0 − ω)Tλ(ω), u = 1/ω and w = 1/θ. Equation (A.2) can be rewritten as 2πα0Eλ(w)w=4w+uKλ(u)u2w2du.\appendix \setcounter{section}{1} \begin{equation} 2\pi\alpha_0\frac{E_{\lambda}(w)}{w} = 4 \int_{w}^{+\infty}\!\frac{u K_{\lambda}(u)}{\sqrt{u^2-w^2}}\, \mbox{d}u. \label{eq:Aeff_omega2} \end{equation}(A.3)The right-hand side of Eq. (A.3) is the well-known Abel integral. It can be thereby solved for Kλ (see e.g., Stover 1995) Kλ(u)=α0u+1w2u2ddw(Eλ(w)w)dw.\appendix \setcounter{section}{1} \begin{equation} K_{\lambda}(u) = -\,\alpha_0\int_{u}^{+\infty}\!\!\frac{1}{\sqrt{w^2-u^2}}\,\frac{\mbox{d}}{\mbox{d}w}\!\left(\frac{E_{\lambda}(w)}{w}\right)\,\mbox{d}w. \label{eq:Aeff_omega3} \end{equation}(A.4)Then, restoring the T function and the θ and ω variables, we obtain Tλ(ω)=α0ω(α0ω)0ωθω2θ2ddθ[θEλ(θ)]dθ.\appendix \setcounter{section}{1} \begin{equation} T_{\lambda}(\omega) = \frac{\alpha_0}{\omega\left(\alpha_0-\omega\right)}\int_{0}^{\omega}\!\!\frac{\theta}{\sqrt{\omega^2-\theta^2}}\,\frac{\mbox{d}}{\mbox{d}\theta}\left[\theta E_{\lambda}(\theta)\right]\:\mbox{d}\theta. \label{eq:Aeff_omega4} \end{equation}(A.5)Finally, we substitute ω = α0 − α1 in Eq. (A.5) with 0 ≤ α1 < α0, and change the integration variable from θ to t, where θ = (α0 − α1)sint. Hence, we obtain the final result Tλ(α1)=α0α10π/2sint[ddθ(θEλ(θ))]θ=θ(t)dt,\appendix \setcounter{section}{1} \begin{equation} T_{\lambda}\left(\alpha_1\right) = \frac{\alpha_0}{\alpha_1}\int_{0}^{\pi/2}\!\!\sin t \,\left[\frac{\mbox{d}}{\mbox{d}\theta}\left(\theta E_{\lambda}(\theta)\right)\right]_{\theta = \theta(t)}\!\!\!\!\!\mbox{d}t, \label{eq:refl_product} \end{equation}(A.6)where the [ ] brackets mean that the enclosed expression has to be evaluated at (α0 − α1)sint. We have so obtained Eq. (17).

Appendix B: Derivation of the vignetting coefficients

We hereafter report the detailed derivation of the vignetting coefficients (Sect. 3.2). The geometry of mirror shell obstructions and the meaning of symbols are explained in Figs. 1, 4, and 5.

B.1. Vignetting at the entrance pupil, V1

This obstruction is caused by the shadow cast by the primary segment of the blocking shell onto the primary segment of the reflective shell. For this to occur, the shadow of the edge of the blocking shell at r0=(RMcosϕ0,RMsinϕ0,L1)\hbox{${\underline r}_0 = (R_{\mathrm M}^*\cos\varphi_0, R_{\mathrm M}^*\sin\varphi_0, L_1^*)$} has to intersect the reflective shell at some r1=(r1cosϕ1,r1sinϕ1,Z1)\hbox{$\underline{r}_1 = (r_1\cos\varphi_1, r_1\sin\varphi_1, Z_1)$}, with 0 < Z1 < L1. The shaded part of the mirror length is then Z1/L1.

To find Z1(ϕ1), we consider the generic ray emerging from the source, which is located at S=(Dsinθ,0,Dcosθ)\hbox{${\underline S} = (D\sin\theta, 0, D\cos\theta)$}, and passing by r0\hbox{${\underline r}_0$}. Now, r1\hbox{${\underline r}_1$}, the ray intersection point on the reflecting primary, whose equation is z1=r1R0α0,withr1>R0,\appendix \setcounter{section}{2} \begin{eqnarray} \label{eq:cone1} z_1 = \frac{r_1-R_0}{\alpha_0}, & \mbox{with}& r_1>R_0, \end{eqnarray}(B.1)in conical approximation, is identified via the vector equation (Sr0)×(r1r0)=0,\appendix \setcounter{section}{2} \begin{equation} \left({\underline S}-{\underline r}_0\right) \times \left({\underline r}_1-{\underline r}_0\right) =0, \label{eq:V1_eq} \end{equation}(B.2)where  ×  denotes a cross product. This returns two independent scalar equations {(1L1D)r1cosϕ1+δcosϕ0θα0(r1R0)=RMcosϕ0L1θ(1L1D)r1sinϕ1+δsinϕ0α0(r1R0)=RMsinϕ0,\appendix \setcounter{section}{2} \begin{equation} \left\{ \begin{array}{l} \left(1-\frac{L_1^*}{D}\right)r_1\cos\varphi_1 + \frac{\delta^*\cos\varphi_0-\theta}{\alpha_0}\left(r_1-R_0\right) = R_{\mathrm M}^*\cos\varphi_0-L_1^*\theta\\[2.5mm] \left(1-\frac{L_1^*}{D}\right)r_1\sin\varphi_1 +\frac{\delta^*\sin\varphi_0}{\alpha_0}\left(r_1-R_0\right) = R_{\mathrm M}^*\sin\varphi_0, \end{array}\right. \label{eq:V1_eqs} \end{equation}(B.3)where we defined δ=RM/D\appendix \setcounter{section}{2} \hbox{$\delta^* = R_{\mathrm M}^*/D$} to be the beam divergence at the blocking mirror shell.

The impact position is obtained by solving the previous equations for r1 and ϕ1. If θ = 0, the solution is expected to be independent of ϕ0, and the first points to be shaded are near the intersection plane, i.e., r1 ≈ R0. To treat the general case of θ ≥ 0, we search for a perturbative solution in the approximation of small θ, so we set r1 = R0 + ξ and ϕ1 = ϕ0 + ε, with 0 < ξ ≪ R0 and  | ε |  ≪ ϕ0. By substituting in Eqs. (B.3), approximating to the first order, and solving, we obtain a linear system whose solutions are ξε\appendix \setcounter{section}{2} \begin{eqnarray} \xi & \simeq & -\frac{\Phi-\left(\alpha_0+\delta-\theta\cos\varphi_0\right)}{\alpha_0+\delta^*-\theta\cos\varphi_0}L_1^*\alpha_0\label{eq:xi0},\\[1.5mm] \varepsilon& \simeq& \frac{L_1^*}{R_0}\,\frac{\delta^*+\Phi}{\left(\alpha_0+\delta^*-\theta\cos\varphi_0\right)}\,\theta\sin\varphi_0\,\label{eq:eps0}, \end{eqnarray}where we have used the definition of Φ (Eq. (23)) and neglected terms proportional to α0δ.

We note that ε is of the order of θ or less, i.e., it is negligible with respect to ϕ0 itself. We can then assume that ϕ1 ≈ ϕ0. Moreover, we know from Sect. 2.1 also that ϕ1 ≈ ϕ2, thus we denote with ϕ the nearly-common value of all these polar angles. We then rewrite Eq. (B.4), using Eq. (B.1), as Z1(ϕ)L1Φ(α0+δθcosϕ)α0+δθcosϕ.\appendix \setcounter{section}{2} \begin{equation} Z_1(\varphi) \simeq -L_1^*\frac{\Phi-\left(\alpha_0+\delta-\theta\cos\varphi\right)}{\alpha_0+\delta^*-\theta\cos\varphi}. \label{eq:z1_0} \end{equation}(B.6)Clearly, Z1 > 0 only if ξ > 0. If the spacing of the two mirrors is small enough, we may also approximate δ ≈ δ. The first vignetting factor is thereby V1=1Z1L1\appendix \setcounter{section}{2} \hbox{$V_1 = 1-\frac{Z_1}{L_1}$}, i.e., recalling the definition of α1 (Eq. (3)), V1(α1)1+L1(Φα1)L1α1,\appendix \setcounter{section}{2} \begin{equation} V_1\left(\alpha_1\right) \simeq 1+ \frac{L_1^*\left(\Phi-\alpha_1\right)}{L_1\alpha_1}, \label{eq:V1} \end{equation}(B.7)if positive and less than 1. We have so obtained Eq. (30). As usual, if this expression returns a negative value at some ϕ′, then V1(ϕ′) = 0, or, if larger than one, V1(ϕ′) = 1. If L1=L1\appendix \setcounter{section}{2} \hbox{$L_1 = L_1^*$} we find that V1 takes the simple form of V1(α1)Φα1,\appendix \setcounter{section}{2} \begin{equation} V_1\left(\alpha_1\right) \simeq \frac{\Phi}{\alpha_1}, \label{eq:V1_altfin} \end{equation}(B.8)where in this case we also note that, if V1 were positive and smaller than one for all ϕ, the geometric area of the obstructed primary segment would become A1(θ)=2R0L10πα1V1(α1)dϕ,\appendix \setcounter{section}{2} \begin{equation} A_1(\theta) = 2R_0L_1\int_0^{\pi}\!\alpha_1V_1\left(\alpha_1\right)\,\mbox{d}\varphi, \label{eq:A1g_obs} \end{equation}(B.9)which – as expected – returns the area of the corona delimited by the largest radii of the two shells A1(θ)=2πR0L1Φ=2πR0(RMRM),\appendix \setcounter{section}{2} \begin{equation} A_1(\theta) = 2\pi R_0L_1\Phi =2\pi R_0\left(R_{\mathrm M}-R_{\mathrm M}^*\right), \label{eq:A1g_obs2} \end{equation}(B.10)where we have used the relation RM ≃ R0 + α0L1.

B.2. Vignetting at the intersection plane, V2

In this case, vignetting occurs after the first reflection from the primary segment of the inner shell, i.e., a point of the blocking shell’s outer surface at z = 0, with coordinates r0=(R0cosϕ0,R0sinϕ0,0)\appendix \setcounter{section}{2} \hbox{${\underline r}_0 = (R_0^*\cos\varphi_0, R_0^*\sin\varphi_0, 0)$}, may intercept a ray reflected by the primary segment of the reflective shell at r1=(r1cosϕ1,r1sinϕ1,Z1)\appendix \setcounter{section}{2} \hbox{$\underline{r}_1=(r_1\cos\varphi_1, r_1\sin\varphi_1, Z_1)$}, with 0 < Z1 < L1. If it does, obstruction occurs at Z1(ϕ1) < z < L1. The coordinate Z1 is located via the vector equation k1×(r1r0)=0,\appendix \setcounter{section}{2} \begin{eqnarray} {\underline k}_1 \times \left({\underline r}_1-{\underline r}_0\right) =0, \label{eq:V2_eq} \end{eqnarray}(B.11)where k1\appendix \setcounter{section}{2} \hbox{$\underline{k}_1$} is the ray direction after the first reflection. This unit vector has the following expression2 (Spiga et al. 2009) k1(θ+(δ2α1)cosϕ1(δ2α1)sinϕ1-1).\appendix \setcounter{section}{2} \begin{equation} \underline{k}_1 \simeq \left( \begin{array}{c} -\theta+\left(\delta-2\alpha_1 \right)\cos\varphi_1\\ \left(\delta-2\alpha_1\right)\sin\varphi_1\\ -1 \end{array}\right). \label{eq:exit_dir1} \end{equation}(B.12)After some manipulation of Eq. (B.11), we obtain the two scalar equations {(α0cosϕ1+k1x)r1=k1xR0+α0R0cosϕ0(α0sinϕ1+k1y)r1=k1yR0+α0R0sinϕ0,\appendix \setcounter{section}{2} \begin{equation} \left\{\begin{array}{ccc} \Big(\alpha_0\cos\varphi_1+k_{1x}\Big)r_1 &=& k_{1x}R_0+\alpha_0R_0^*\cos\varphi_0 \\[1.5mm] \Big(\alpha_0\sin\varphi_1+k_{1y}\Big)r_1 &=& k_{1y}R_0+\alpha_0R_0^*\sin\varphi_0, \end{array}\right. \label{eq:V2_eqs} \end{equation}(B.13)and, using Eq. (B.12), we find that the solution of Eqs. (B.13) for θ = 0 is ϕ0 = ϕ1 and r1(0)=R0+α0α0+δ(R0R0).\appendix \setcounter{section}{2} \begin{equation} r_1^{(0)} = R_0 + \frac{\alpha_0}{\alpha_0+\delta}\left(R_0-R_0^*\right). \label{eq:soluz0} \end{equation}(B.14)Since R0>R0\appendix \setcounter{section}{2} \hbox{$R_0 > R_0^*$}, we have r1(0)>R0\appendix \setcounter{section}{2} \hbox{$r_1^{(0)} > R_0$}, hence Z1(0)>0\appendix \setcounter{section}{2} \hbox{$Z_1^{(0)} > 0$} (Eq. (B.1)). If θ > 0, we set r1=r1(0)+ξ\appendix \setcounter{section}{2} \hbox{$r_1 = r_1^{(0)}+\xi$}, ϕ0 = ϕ1 + ε, and proceed as in Sect. B.1. The solution to a first order approximation is ξα0(α0+δα1)α1(α0+δ)(R0R0).\appendix \setcounter{section}{2} \begin{equation} \xi \simeq \frac{\alpha_0\left(\alpha_0+\delta-\alpha_1\right)}{\alpha_1\left(\alpha_0+\delta\right)}\left(R_0-R_0^*\right). \label{eq:xi1} \end{equation}(B.15)We are not interested in the exact expression of ε, which is of the order of θ, so we can assume again that ϕ0 ≈ ϕ1 and neglect the ϕ’s subscript. In contrast, from ξ and Eq. (B.14) we can derive an expression for r1, and using Eq. (B.1) we obtain Z1Z1(ϕ)R0R0α1,\appendix \setcounter{section}{2} \begin{equation} Z_1(\varphi) \simeq \frac{R_0-R_0^*}{\alpha_1}, \label{eq:z1_1} \end{equation}(B.16)which is always non-negative. All points with z > Z1 are then obstructed. The resulting vignetting coefficient is V2 = Z1/L1, i.e., using the definition of Ψ (Eq. (24)), V2(α1)Ψα1·\appendix \setcounter{section}{2} \begin{equation} V_2\left(\alpha_1\right) \simeq \frac{\Psi}{\alpha_1}\cdot \label{eq:V2} \end{equation}(B.17)We have so obtained Eq. (31).

B.3. Vignetting at the exit pupil, V3

In this case, the obscuration occurs after the second reflection, on the secondary segment of the blocking shell. A generic point of the blocking shell at z=L2\appendix \setcounter{section}{2} \hbox{$z = -L_2^*$}, i.e., r0=(Rmcosϕ0,Rmsinϕ0,L2)\appendix \setcounter{section}{2} \hbox{${\underline r}_0 = (R_{\mathrm m}^*\cos\varphi_0, R_{\mathrm m}^*\sin\varphi_0, -L_2^*)$}, may intercept a ray after it was reflected at r2=(r2cosϕ2,r2sinϕ2,Z2)\appendix \setcounter{section}{2} \hbox{$\underline{r}_2 = (r_2\cos\varphi_2, r_2\sin\varphi_2, Z_2)$}, with  − L2 < Z2 < 0. If it does, the point r2\appendix \setcounter{section}{2} \hbox{$\underline{r}_2$} is located along with the usual vector equation k2×(r2r0)=0,\appendix \setcounter{section}{2} \begin{equation} {\underline k}_2 \times \left({\underline r}_2-{\underline r}_0\right) =0, \label{eq:V3_eq} \end{equation}(B.18)where the z-coordinate of the secondary segment is given by z2=r2R03α0,withr2<R0,\appendix \setcounter{section}{2} \begin{eqnarray} \label{eq:cone2} z_2 = \frac{r_2-R_0}{3\alpha_0}, & \mbox{with}& r_2<R_0, \end{eqnarray}(B.19)and k2\appendix \setcounter{section}{2} \hbox{$\underline{k}_2$} is the direction of the ray after the second reflection, k2=k12(k1·n2)n2,\appendix \setcounter{section}{2} \begin{equation} \underline{k}_2 = \underline{k}_1 -2\left(\underline{k}_1\cdot \underline{n}_2\right)\, \underline{n}_2, \label{eq:k2_def} \end{equation}(B.20)because the tangential component is conserved, while the normal component to the surface reverses its sign in the reflection. To the small angle approximation, since ϕ2 ≃ ϕ1 (Sect. 2.1), n2(cosϕ1sinϕ13α0),\appendix \setcounter{section}{2} \begin{equation} \underline{n}_2 \approx \left( \begin{array}{c} -\cos{\varphi_1}\\ -\sin{\varphi_1}\\ 3\alpha_0 \end{array} \right), \label{eq:n2} \end{equation}(B.21)and k1\appendix \setcounter{section}{2} \hbox{$\underline{k}_1$} is given by Eq. (B.12). The final ray direction3 is then k2(θ+(δ4α0)cosϕ1(δ4α0)sinϕ1-1).\appendix \setcounter{section}{2} \begin{equation} \underline{k}_2 \simeq \left( \begin{array}{c} -\theta+\left(\delta-4\alpha_0\right)\cos\varphi_1\\ \left(\delta-4\alpha_0\right)\sin\varphi_1\\ -1 \end{array} \right). \label{eq:exit_dir2} \end{equation}(B.22)Setting ϕ2 = ϕ1 in Eq. (B.18), we obtain the two independent equations {(3α0cosϕ1+k2x)r2=k2xRm+3α0Rmcosϕ0(3α0sinϕ1+k2y)r2=k2yRm+3α0Rmsinϕ0,\appendix \setcounter{section}{2} \begin{equation} \left\{\begin{array}{ccc} \Big(3\alpha_0\cos\varphi_1+k_{2x}\Big)r_2 &=& k_{2x}R_{\mathrm m}'+3\alpha_0R_{\mathrm m}^*\cos\varphi_0\\[1.5mm] \Big(3\alpha_0\sin\varphi_1+k_{2y}\Big)r_2 &=& k_{2y}R_{\mathrm m}'+3\alpha_0R_{\mathrm m}^*\sin\varphi_0, \end{array}\right. \label{eq:V3_eqs} \end{equation}(B.23)where we set RmR03α0L2\appendix \setcounter{section}{2} \hbox{$R_{\mathrm m}' \simeq R_0-3\alpha_0L_2^*$}. The solution for θ = 0 is ϕ0 = ϕ1 and r2(0)=Rm+3α0α0δ(RmRm).\appendix \setcounter{section}{2} \begin{equation} r_2^{(0)}= R_{\mathrm m}'+\frac{3\alpha_0}{\alpha_0-\delta}\left(R_{\mathrm m}'-R_{\mathrm m}^*\right). \label{eq:soluz0_2} \end{equation}(B.24)If θ > 0, we set r2=r2(0)+ξ\appendix \setcounter{section}{2} \hbox{$r_2 = r_2^{(0)}+\xi$}, ϕ0 = ϕ1 + ε as in Sects. B.1 and B.2. Substituting this into Eq. (B.23), approximating at the first order and solving the linear system, we obtain ϕ0 ≈ ϕ1 and the expression of ξ. Adding this to r2(0)\appendix \setcounter{section}{2} \hbox{$r_2^{(0)}$} we obtain r2=Rm+3α0α2(RmRm),\appendix \setcounter{section}{2} \begin{equation} r_2= R_{\mathrm m}'+\frac{3\alpha_0}{\alpha_2}\left(R_{\mathrm m}'-R_{\mathrm m}^*\right), \label{eq:soluz1_2} \end{equation}(B.25)where α2 is given by Eq. (4). The secondary mirror region at ϕ1 with radial coordinate between r2 and R0 is then obstructed: this is equivalent to an obscuration on the primary mirror between z = 0 and some z = Z1 < L1 (Fig. 4). The correspondence between r1 and r2 is given by the equation (Spiga et al. 2009) r22(α0+α1)R03α1r1α2·\appendix \setcounter{section}{2} \begin{equation} r_2 \simeq \frac{2\left(\alpha_0+\alpha_1\right)R_0-3\alpha_1r_1}{\alpha_2}\cdot \label{eq:r2} \end{equation}(B.26)By comparing Eqs. (B.25) and (B.26), we find r1, and the corresponding value of Z1 using Eq. (B.1), Z1(α2+3α0)L2(R0Rm)α1,\appendix \setcounter{section}{2} \begin{equation} Z_1 \simeq \frac{\left(\alpha_2+3\alpha_0\right)L_2^*-\left(R_0-R_{\mathrm m}^*\right)}{\alpha_1}, \label{eq:z1_2} \end{equation}(B.27)from which, recalling that α2 = 2α0 − α1 (Eq. (5)), we obtain the vignetting coefficient, V3 = 1 − Z1/L1, V3(α1)1+L2(Σα2)L1α1,\appendix \setcounter{section}{2} \begin{equation} V_3\left(\alpha_1\right) \simeq 1+\frac{L_2^*\left(\Sigma-\alpha_2\right)}{L_1\alpha_1}, \label{eq:V3} \end{equation}(B.28)where we have used the definition of Σ (Eq. (25)). We have thereby found Eq. (32).

All Figures

thumbnail Fig. 1

Sketch of a grazing-incidence Wolter-I mirror shell with an off-axis X-ray source. We also show a ray direction vector before reflection (k0)\hbox{$(\underline{k}_0)$}, after the first reflection (k1)\hbox{$(\underline{k}_1)$}, and after two reflections (k2)\hbox{$(\underline{k}_2)$}.

In the text
thumbnail Fig. 2

The product of the reflectivities of the two segments of a Wolter-I shell with α0 = 0.212 deg, at 30 keV. The reflective surface is a multilayer coating, a Pt/C graded stack with 200 couples of layers. The d-spacing in the multilayer follows the power-law model d(k) = a   (b + k) − c (Joensen et al. 1995), with a = 77.4 Å, b =  −0.94, c = 0.223, a Pt thickness ratio Γ = 0.42, and a 4 Å roughness rms. The solid line is directly computed from the multilayer structure, whilst the dots are the result of the inverse computation from the effective area in Fig. 3.

In the text
thumbnail Fig. 3

The effective area at 30 keV (λ = 0.4 Å) of a Wolter-I mirror shell with R0 = 148.5 mm, L = 300 mm, f = 10 m, α0 = 0.212 deg, and a reflectivity product as shown in Fig. 2 (solid line). The effective area was used to re-derive the Tλ(α1) function (Fig. 2, dots).

In the text
thumbnail Fig. 4

Obstruction in an assembly of mirror shells. Rays impinging a mirror shell can be blocked by the adjacent shell with smaller diameter in only three ways, as listed in the text. The impact points are highlighted. Other mirror shells are not shown. The obscured regions of the primary mirror are indicated in gray: angles and the mirror spacings are greatly exaggerated.

In the text
thumbnail Fig. 5

Geometrical meaning of the obstruction parameters, Φ, Ψ, and Σ, for a pair of Wolter-I nested shells.

In the text
thumbnail Fig. 6

Vignetting coefficients as a function of the polar angle for L2 = 310 mm, L1 = 300 mm, L2\hbox{$L_2^*$} = 290 mm, L1\hbox{$L_1^*$} = 280 mm, α0 = 0.3 deg, δ = 0 deg, θ = 0.25 deg, Φ = 0.43 deg, Ψ = 0.4 deg, and Σ = 0.408 deg. We also plot V, the vignetting factor for double reflection (Eq. (6)).

In the text
thumbnail Fig. 7

Vignetting coefficients as a function of the polar angle for the same configuration as in Fig. 6, but this time the source is at a finite distance (D = 127 m, δ = 0.095 deg). We note that the obstruction at the exit pupil is less severe than for a source at infinity.

In the text
thumbnail Fig. 8

There is no vignetting at the intersection plane if the distance of the ray reflected at the maximum and the minimum diameter is smaller than the spacing of the two mirror shells.

In the text
thumbnail Fig. 9

Normalized geometric area, A(θ)/A(0) of an obstructed mirror shell with α0 = 0.2 deg, δ = 0, and L1 = L2, as a function of the off-axis angle, for different values of the obstruction parameter Φ, with Φ ≈ Ψ ≈ Σ and L = L ∗ . The curves are traced using the analytical formulae reported in Sect. 4.1. The curve for the unobstructed mirror (the solid line) is also valid for any obstructed mirror with Φ > 2α0.

In the text
thumbnail Fig. 10

Initial positions and destinations of 40 000 rays at the entrance pupil (points) for a Wolter-I mirror shell with L = L ∗ , and the angular parameters reported in the legend (Φ ≈ Ψ ≈ Σ). The radial scale has been expanded. Only rays that would have struck the primary mirror were traced. The limits of the regions of different vignetting (dashes) are computed from the vignetting coefficients.

In the text
thumbnail Fig. 11

Comparison between some analytical curves of Fig. 9 (lines), and the results of an accurate ray-tracing (symbols) for an obstructed Wolter-I mirror with the same α0 = 0.2 deg, R0 = 139.6 mm, L = L ∗  = 300 mm. The outer radius of blocking shell, R0\hbox{$R_0^*$}, takes on the values 137.24, 138.29, and 138.55 mm, for Φ to match the values used to draw the analytical curves (but the unobstructed one). The accord between the formulae and the tracing is excellent within the statistical error.

In the text
thumbnail Fig. 12

Effective area of an obstructed mirror shell with R0 = 148.5 mm, f = 10 m, α0 = 0.212 deg, L = L = 300 mm, and the X-ray source at infinite distance, 6 arcmin off-axis. The curves are computed for three different values of R0\hbox{$R_0^*$}. The reflective shell has a Pt/C graded multilayer, whose d-spacing variation in depth is a power law, d(k) = a   (b + k) − c, where k = 1 ... 200, a = 115.5 Å, b = 0.9, c = 0.27. We moreover assumed a Pt thickness-to-d-spacing ratio Γ = 0.35, a roughness rms σ = 4 Å, and a stack that ends with Pt on top.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.