Molecular and particulate contamination of space optical surfaces can be extremely detrimental to optical system performance. In order to determine whether optical system is unusable, we must do a lot of tests. For reducing experimental effort, computational criterion for unusable space optical system is put forward, which is based on experimentation and calculation by 2D-FDTD (finite-difference time domain). And then factors making optical system unusable are studied when distribution of light intensity is consistent with the criterion. Results show that small particulate is more detrimental to optical system than big one for particle with high permittivity is. The distribution of particles making optical system unusable is also obtained.
1. Introduction
There are two types of contamination that affect optical performance: molecular and particulate contamination. Particulate contamination can occur through two phenomenons: obscuration and scattering effects. Particulate obscuration occurs when surface particles prevent light from reaching the underlying surface [1]. For particles with high permittivity, particulate obscuration plays a main role. Here, particulate with high permittivity is research subject.
Molecular and particulate contamination of spacecraft optical surfaces can be extremely detrimental to optical system performance [2]. Surface-deposited particles induce not only decreasing transmittance of optical system, but also reducing the image quality of optical system [3]. For system of space optical communication, contamination of optical surface can change the light intensity distribution [4].
As particles on optical surface accumulate, the space optical system will be unusable at last. Under what conditions the space optical system is unusable is concerned. As we know, particle size, the space between particles and absorbency of particles can influence quality of optical system. Then we should do a lot of experiments to ascertain whether optical system is unusable. In order to save money and time, a computational criterion for unusable space optical system is useful. In this paper, the criterion is put up. And then, for particle with high permittivity, factors making optical system unusable are studied.
2. The Criterion for Unusable Optical System
Kahnert [5] has analyzed some of the most widely used numerical techniques in detail. In these models, FDTD can resolve scattering by inhomogeneous, anisotropic, and arbitrarily shaped charged or uncharged particles [6], so that it has been applied to such diverse applications as antenna scattering [7], analysis of microstrip circuits [8], absorption in tissue [9], optical absorption by biology cellula [10–13] or scattering by atmospheric ice crystals [14], irregular aggregates [15], multifaceted concave and convex particles [16], spheres with nonspherical inclusions [17], hexagonal columns with air bubble and soot inclusions [18] or clusters of particles with or without inclusions [19], consistent and rigorous analysis of infinite graphene layers is developed [20], electrically thin dispersive layers [21], calculation of reflection coefficient for the plane wave incident on planar periodic structures [22], light scattering and absorption by nonspherical particles embedded in an absorbing dielectric medium [23], an arbitrary beam’s interaction with an arbitrary dielectric surface [24], to name but a few examples. So FDTD is chosen. Because this calculation is only used to reduce experimental effort, rather than for exact simulation, 2D-FDTD is appropriate. And 2D-FDTD can also decrease calculation workload.
According to resolution test, when the diameter (2r) of particulates is 0.3 mm, the maximal space (dmax) between particles, resulting in optical system being unusable, reaches 0.4 mm. In this case, far-field light intensity distribution calculated by 2D-FDTD is shown in Figure 1, and the relevant near-field distribution of light intensity is constant. Furthermore, light intensity distributions are also calculated when space between particles is less than 0.4 mm and diameter of particulates is more than 0.3 mm. Results show that both far-field light intensity distributionand near-field distribution are in accord with that of 2r=0.3 mm and d=0.4 mm [25]. Then a criterion for unusable optical system is put forward. That is to say, optical system can be thought unusable once near-field distribution of light intensity is constant and far-field light intensity distribution is like Figure 1. In [25], the resolution power is obtained by software HYRes 3.1, which is written by Olympus Corporation.
Distribution of transmitted light intensity while 2r=0.3 mm and d=0.4 mm.
3. Factors Making Optical System Unusable
As FDTD equations display, ε and Cs are factors. ε and Cs are different for different contaminants. That is to say, the kind of contaminant is a factor.
For the kind of contaminant (e.g., ε=1.5+5.0i), both the size of particulates and the space between particles are factors. Here, when distribution of light intensity is consistent with the criterion, the relationship between the maximal particle space and radius of particles is research emphasis.
3.1. Relationship between Radius of Particles and the Maximal Particle Space
Here, relevant maximal space between particles (dmax) is obtained, firstly, for unusable optical system, when radius of particulates varies from 0.0325 mm to 1.25 mm (error is less than 1 μm). See Figure 2.
Relationship between radius of particles and maximal particle space.
From Figure 2, the relationship between radius of particles and maximal space between particles is linear on the whole, although there are some exceptions. That is to say, as radius of particulates increases, maximal space between particles increases linearly. This is reasonable because the diameter of the Airy disk varies inversely with the particle’s radius, the smaller the particles are, the larger the Airy disk is. Then the change of light intensity overlap is less susceptible to the particle size. dmax in Figure 2 has smaller change with r. As the particle size becomes bigger, diffraction image of one particle becomes sharper. Then dmax changes with r rapidly.
From the linear relationship between dmax and r, it may be deduced that total coverage of contaminant, which causes the optical system to be unusable, is equal for all size particles. Therefore, ratio of opaque part to all is calculated when radius of particles is different.
3.2. Relationships between r and 2r/dmax
In order to study how total coverage of contaminant influences optical imaging system, 2r/dmax is calculated when radius of particles, r, is different. See Figure 3.
Relationship between r and 2r/dmax.
From Figure 3, total coverage of contaminant is not constant, which means total coverage of contaminant is not so important to imaging. It can also be seen from Figure 3 that the less radius of particles is, the less the quantity of contaminants which cause optical system to be unusable is. After r reaches a certain size, 2r/dmax moves up but slowly with the increasing of r. So a smaller quantity of smaller than larger particulates on optical surface can induce optical system to become unusable. In a word, small particulates are more detrimental to optical system than large particulates.
So the small particle’s contamination prevention should be the most important of particle’s contamination prevention.
3.3. Relationship between r and dmax – 2r
Because of particle’s extinction, dmax – 2r may be an important parameter for transmissivity of optical system. And then, dmax – 2r is calculated correspondingly, when distribution of light intensity is consistent with the criterion. Results are plotted in Figure 4.
Relationship between r and dmax – 2r.
Result of Figure 4 is consistent with that of 3.1. Furthermore, there is an interesting fact in Figure 4, that is, dmax – 2r is always less than 0.36 mm. This means that according to the criterion put up in this paper, it is possible to obtain an exact value of resolution power from HYRes 3.1 when the space between particles is more than 0.36 mm, although resolution power of optical system will decrease somewhat. However, it is not quite clear why this value is 0.36 mm.
4. Conclusions
In this paper, a criterion is put up based on test and calculation and it is used to determine whether optical system is usable or not in calculation. Then factors which may induce optical system unusable are studied. For particle with ε=1.5+5.0i, the following results can be deduced when optical system is unusable: (1) maximum space between particles is linear with radius of particles approximately, (2) the smaller the radius of particles is, the smaller is the minimal coverage of particle, (3) small particles are more detrimental to optical system than large particles are, and (4) dmax – 2r is always less than 0.36 mm.
In conclusion, when distribution of light intensity is consistent with the criterion, the smaller the size of particulates is, the smaller the maximum space between particles is. But this is not to say that optical system will be unusable so long as the coverage of particulates reaches a value, whether the particles on optical surface are large or small. In fact, the lesser radius of particle is, the less minimal coverage of particles is, when optical system becomes unusable. That is, small particulates are more detrimental to optical system than large particulates. From the above study, if dmax – 2r is more than 0.36 mm, resolution power can be read out by HYRes 3.1.
The results are benefit to study relevant problem of similar structure, such as periodic arrays of left-handed metamaterials [26], monolayer and bilayer graphene deposited on a SiO2/Si substrate [27], and multiwall carbon nanotubes [28]. And they may be of guidance for designs of some components, for example, design of graphene magnetic field sensors [29].
The paper presents that if near-field light intensity from FDTD is equal to the fixed value, optical system is hardly usable. However, the criterion is possibly not the best.
HoustonK. D.Comparative mirror cleaning study: a study on removing particulate contamination6291Optical Systems Degradation, Contamination, and Stray Light: Effects, Measurements, and Control IIAugust 2006629107Proceedings of SPIE2-s2.0-3375059522010.1117/12.683231FolkmanM. A.DarntonL. A.SilottS. G.FrinkM. E.Optical scatter and contamination effects facility1329Optical System Contamination: Effects, Measurement, Control IIJuly 1990221232Proceedings of SPIE2-s2.0-0025551741WuD.ZhouY.GuoS.Grain effect on imaging of spatial optical system6595Fundamental Problems of Optoelectronics and Microelectronics IIISeptember 2006659529Proceedings of SPIE2-s2.0-3424900338410.1117/12.726458WuD.ZhouY.Relationship between antenna contamination and laser wavelength in optical communication20093635645682-s2.0-6334909189710.3788/CJL20093603.0564KahnertF. M.Numerical methods in electromagnetic scattering theory200379-807758242-s2.0-003761604510.1016/S0022-4073(02)00321-7SunW.2000AndrewW. V.BalanisC. A.TirkasP. A.PengJ.BirtcherC. R.Finite-difference time-domain analysis of HF antennas on helicopter airframes19973921001132-s2.0-0031147449SheenD. M.AliS. M.AbouzahraM. D.KongJ. A.Application of the three-dimensional finite-difference time-domain method to the analysis of planar microstrip circuits19903878498572-s2.0-002546414810.1109/22.55775SullivanD. M.BorupD. T.GandhiO. P.Use of the finite-difference time-domain method in calculating EM Absorption in human tissues19873421481572-s2.0-0023289454DunnA.Richards-KortumR.Three-dimensional computation of light scatteing from cells1996248989052-s2.0-0030360920DunnA.SmithpeterC.WelchA. J.Richards-KortumR.Finite-difference time-domain simulation of light scattering from single cells1997232622662-s2.0-0031185944DrezekR.DunnA.Richards-KortumR.Light scattering from cells: finite-difference time-domain simulations and goniometric measurements19993816365136612-s2.0-0001560546DrezekR.DunnA.Richards-KortumR.A pulsed finite-difference time-domain (FDTD) method for calculating light scattering from biological cells over broad wavelength ranges2000671471572-s2.0-0000939545YangP.LiouK. N.Finite-difference time domain method for light scattering by small ice crystals in three-dimensional space19961310207220852-s2.0-0030269856VideenG.SunW.FuQ.Light scattering from irregular tetrahedral aggregates19981561–3592-s2.0-0032208831YangP.LiouK. N.MishchenkoM. I.GaoB. C.Efficient finite-difference time-domain scheme for light scattering by dielectric particles: application to aerosols20003921372737372-s2.0-0001300317MishchenkoM. I.HovenierJ. W.TravisL. D.2000San Diego, Calif, USAAcademicYangP.LiouK. N.Light scattering by hexagonal ice crystals: comparison of finite-difference time domain and geometric optics models19951211621762-s2.0-0029236096SunW.FuQ.ChenZ.Finite-difference time-domain solution of light scattering by dielectric particles with a perfectly matched layer absorbing boundary condition19993815314131512-s2.0-0032606342BouzianasG. D.KantartzisN. V.TsiboukisT. D.Subcell dispersive finite-difference time-domain schemes for infinite graphene-based structures201264377386KärkkäinenM. K.Subcell FDTD modeling of electrically thin dispersive layers2003516177417802-s2.0-003861473410.1109/TMTT.2003.812584AminianA.Rahmat-SamiiY.Spectral FDTD: a novel technique for the analysis of oblique incident plane wave on periodic structures2006546181818252-s2.0-3374520510910.1109/TAP.2006.875484SunW.LoebN. G.FuQ.Finite-difference time-domain solution of light scattering and absorption by particles in an absorbing medium20024127572857432-s2.0-0037144798SunW.PanH.VideenG.General finite-difference time-domain solution of an arbitrary electromagnetic source interaction with an arbitrary dielectric surface20094831601560252-s2.0-7044964706510.1364/AO.48.006015DakunW.YanpingZ.Performance test and analysis for imaging system of contaminated optical surface20103024114152-s2.0-7794990666210.3788/AOS20103002.0411BandlowB.SchuhmannR.LubkowskiG.WeilandT.Analysis of single-cell modeling of periodic metamaterial structures2008446166216652-s2.0-4404908712410.1109/TMAG.2007.916037AbergelD. S. L.RussellA.Fal'KoV. I.Visibility of graphene flakes on a dielectric substrate20079162-s2.0-3454784254210.1063/1.2768625063125MaffucciA.MianoG.VilloneF.ZamboniW.Analysis of multiwall carbon nanotubes using a three-dimensional integral formulation and a fluid model2008446161416172-s2.0-4404910914810.1109/TMAG.2007.915295PisanaS.BragancaP. M.MarineroE. E.GurneyB. A.Graphene magnetic field sensors2010466191019132-s2.0-7795281457810.1109/TMAG.2010.2041048