Simulation and Analysis of Co-phasing Errors of the Segmented Primary Mirror Tiled by Petal-shaped Segments in LOT*

2022-12-12 08:24SHENShidongCUIXiangqunZHANGYong
天文学报 2022年6期

SHEN Shi-dong CUI Xiang-qun ZHANG Yong

(1 National Astronomical Observatories/Nanjing Institute of Astronomical Optics & Technology,Chinese Academy of Sciences,Nanjing 210042)(2 CAS Key Laboratory of Astronomical Optics & Technology,Nanjing Institute of Astronomical Optics & Technology,Nanjing 210042)(3 University of Chinese Academy of Sciences,Beijing 100049)

ABSTRACT In order to fulfill the requirement for astronomical observation in the opticalinfrared band using large universal aperture optical telescopes,China is planning to construct a 12-meter large aperture optical-infrared telescope(LOT).LOT is of particular significance to Chinese astronomy and would be co-phased in the near-infrared waveband.The primary mirror of LOT will be constructed from either hexagonal segments or petal-shaped segments.Analysis of co-phasing errors using a segmented primary mirror tiled by petal-shaped segments for LOT has been presented.The proposed primary mirror of LOT has 72 petal-shaped segments and a 4.5-meter circular mirror in the center.Detailed theoretical derivation of co-phasing errors and simulation about out-plane errors separately and comprehensively have been provided.SR(Strehl Ratio)of LOT is greater than 0.8 when the RMS(Root Mean Square)value of normally distributed tip-tilt errors is less than 0.016′′or the RMS value of normally distributed piston error is less than 42.5 nm at the co-phasing wavelength 1 μm.The simulation results are of great importance for the segmented primary mirror with active optics in LOT.

Key words telescopes,methods:analytical,methods:data analysis,techniques:high angular resolution

1 Introduction

With the rapid development of astronomical science and technology,and the urgent need for large telescopes with higher light collecting capacity and resolution,many large aperture opticalinfrared telescopes have been built or are under construction both in China and abroad.The main approaches for the design of large aperture telescopes include:(1)the thin mirror active optics approach used in the VLT(Very Large Telescope)of the European Southern Observatory[1];(2)the segmented mirror active optics approach used in the KECK(Keck telescope)of the United States[2];and(3)the combined thin mirror and segmented mirror active optics approach pioneered by LAMOST(Large Sky Area Multi-Object Fiber Spectroscopic Telescope)of China[3–5].However,due to mirror size limitations,the current aperture limit of a single mirror is 8.4 m,which was successfully constructed by the Mirror Laboratory at University of Arizona.In order to construct larger or giant aperture telescopes,segmented mirror active optics has become the approach of choice taking into consideration factors such as technical difficulty,cost,and construction time.

To fulfill the requirements for astronomical observation in the optical-infrared band with large universal aperture optical telescopes,CFGT(Chinese Future Giant Telescope)was proposed in 2000,using a 30-meter primary mirror design with fan-shaped segments[6–11].Japan subsequently proposed JELT(Japan Extremely Large Telescope),the primary mirror of which would also have been 30 m in size[12].Japan has since built a middle aperture telescope,namely the Seimei telescope.The Seimei telescope is a 3.8-meter opticalinfrared telescope,the primary mirror of which consists of 18 petal-shaped segments[13–14].

Chinese astronomers proposed LOT(12-meter large aperture optical-infrared telescope)in 2016 and have considered the use of a segmented primary mirror[15–17].However,there are two proposed segmentation schemes:one type is a segmented primary mirror tiled into 84 hexagonal segments which is commonly used in segmented mirror telescopes,and the other is an innovative design using a 4.5-meter circular center mirror and 72 petal-shaped segments.Fig.1 depicts two types of segmentation schemes.

Fig.1 Two types of segmentation schemes

The segmented primary mirror tiled into petal-shaped segments has numerous advantages compared to the segmented primary mirror tiled into hexagonal segments.For example,the primary mirror can have a circular profile,and the type of segments are consistent in every ring,regardless of physical size or profile.The number of petalshaped segments required is also fewer than would be required using hexagonal segments for the same size of primary mirror,even though the area of an individual hexagonal segment and a petal-shaped segment is almost the same[18–20].This paper will therefore focus on developing an innovative scheme as described in Table 1,in which the size parameters and quantity of each segment are introduced in detail.We can also consider adopting all fanshaped segments for LOT.However,changing the inner 4.5-meter diameter to fan-shaped segments will not affect the research described in this paper from a practical point of view.Furthermore,a perfect co-phasing wavefront in visible imaging can be obtained by using this 4.5-meter mirror;thin mirror active optics technology can also be used in the central mirror to achieve a combination of thin mirror active optics and fan-shaped segmented mirror active optics if necessary.

Table 1 The primary mirror of LOT tiled with fan-shaped segments

The development of LOT is of great significance to Chinese astronomy,and the segmented mirror active optics approach is the preferred choice for LOT taking into consideration co-focus in the visible band and co-phase in the infrared band[21–23].Theoretical derivation and simulation of co-phasing errors are presented in this paper when the primary mirror of LOT reaches co-phase at 1 μm[24–27].

Error analysis and simulation are the prerequisites and necessary conditions for error allocation during the feasibility evaluation for any large aperture segmented telescope.Chanan et al.[28]got SR(Strehl Ratio)for KECK as a function of the RMS(Root Mean Square)segment co-phasing error.Besides,Troy et al.[29]analyzed the influence of gaps,secondary mirror obstruction and reflectivity of mirrors for TMT(Thirty Meter Telescope).In addition,Yaitskova et al.[30–32]derived the analytical expression for formulas describing highly segmented telescopes,established an optical calculation model of PSF(Point Spread Function)based on FFT(Fast Fourier Transform),and simulated the influence of piston and tip-tilt on SR of the telescope,providing data support for the construction of ELT(European Extremely Large Telescope).

In order to evaluate our proposed scheme and study the sensitivity of co-phasing errors,in this study we undertake theoretical derivation and simulation analysis of co-phasing errors in LOT in which the primary mirror is tiled by petal-shaped segments.The results provide a reference and basis for the selection and construction of the primary mirror shape scheme of the LOT.

2 Analysis of co-phasing errors

Fig.2 introduces the global coordinate system of the primary mirror and the local coordinate system of the segment.TheXCYplane is the projection plane of the primary mirror in Fig.2,and theZ-axis is the direction of the optical axis perpendicular to theXCYplane.PointCis the origin of the global coordinate systemXY Z.The origin point,o,of the local coordinate system(xyz)of each segment is located at the center of circumcircle of the segment.They-direction is the radial direction of the primary mirror and thez-direction is the normal direction of a segment,wherexis perpendicular toyandz.

Fig.2 XCY plane in the primary mirror global coordinate system and yoz plane in the segment local coordinate system

Wis assumed to be the wavefront of the primary mirror,andΔWis the wavefront error of the primary mirror.ωiis the wavefront of theith segment,and(Δω)iis the wavefront error of theith segment.The formula to describe the approximate relationship between the Strehl Ratio and wavefront is as follows[33]:

whereσwfis the RMS value of wavefront error(in wavelength units,λ)in Eq.(1).

In addition,σwfcan be expressed as[34]:

where

Note that in Eq.(3),Arepresents the integral region of the primary mirror.

In Eq.(4)and Eq.(5),Nrepresents the number of segments,and dzirepresents the integral region of the segment.

The radius of the circumcircle of a petalshaped segment isa,the center of which is the center of its circumcircle.If the polar coordinate system(ρ,θ)is considered with regard to the local coordinate system whereρ=r/a,the surface of segment in the local coordinates can be expressed as a expansion[35]:

αmnis the local coordinate expansion coefficient,bothmandnare integers andm≥n≥0 andm-n=even in Eq.(6).

The radius of curvature of the primary mirror,Rc,is-38400 mm.prepresents an off-axis distance from the center of a segment to the center of the primary mirror,andp=(wherexoandyoare the coordinates of pointoin Fig.2).The conic constant of the primary,k,has a value of-0.9837843.α20is the coefficient of the first term in Eq.(6),represents defocus.The coefficient of the second term in Eq.(6)that represents astigmatism,α22,iska2p2/4.The coefficient of the third term in Eq.(6),representing coma,α31,iska3p/2.In this study,we primarily focus on the effects of astigmatism as this dominates in terms of wavefront aberrations[36].

2.1 The influence of piston

If there is some random piston errorαiwhich obeys normal distributionNof theith segment,the standard deviationσαis the RMS value of piston error,we can obtain the wavefront of theith segment:

The wavefront of the segmented primary mirror can also be obtained:

wherezirepresents the integral region of theith segment and segirepresents area of theith segment.By substituting this result in Eq.(2),it can be concluded that:

At this time,the RMS value of piston error is 0.0375λwhen SR is 0.8 orσwfis 0.075λ.The RMS value of piston should be less than 37.5 nm when the co-phasing wavelength is 1 μm.

2.2 The influence of tip-tilt

We use the notationai,γi,βito represent the radius of circumcircle of theith petal-shaped segment,the amount of tilt in thex-direction of theith petal-shaped segment(inλunits),and the amount of tilt in they-direction of theith petal-shaped segment(inλunits)respectively.Ifγi~Nandβi~N,the standard deviationσγis the RMS value of tip and the standard deviationσβis the RMS value of tilt,the wavefront of theith petal-shaped segment can be obtained:

Assuming that the standard deviation of the normal distribution of tip-tilt is the same,we can conclude that:

The RMS value of tip-tilt errorσis 0.065λwhen SR is 0.8 orσwfis 0.075λ.The RMS value of tiptilt error should be less than 0.016′′(σwf/ai)when the co-phasing wavelength is 1 μm.The theoretical co-phasing errors of the two segmented mirror schemes using hexagonal segments[37]and petalshaped segments are compared in Table 2.

Table 2 Comparison of RMS value ofco-phasing errors between two segment shapes

2.3 The influence of in-plane errors

Decenter inxaxis can be ignored among the three in-plane errors,based on the rotational symmetry of the quadric surface.Here,the radial error decenter inyaxis and clocking error are mainly analyzed.We assume that the radial error of radial distancepof theith fan-shaped segment can be expressed as(Δp)i(inλunits),and,σΔpis the RMS value of radial translation error.The radial distance of fanshaped segments in the first ring isp1,with a value of 2.919352 m.The radial distance of fan-shaped segments in the second ring isp2,4.160594 m and the radial distance of fan-shaped segments in the third ring isp3,5.404607 m,then it can be calculated as follows:

wherepiis the radial distance of theith segment.

In this case,theσwfvalue is 0.075λwhen SR is 0.8.The RMS value of radial translation errorσΔpis 6162.8λ.The RMS value of radial translation error of each fan-shaped segment should be less than 6.162 mm when the co-phasing wavelength is 1 μm.

Assuming that the clocking error of theith fan-shaped segment can be expressed as(Δθ)i(radian units),and(Δθ)i~N,the wavefront error of theith fan-shaped segment can be obtained:

The house on Beechwood Avenue had a big Front yard with shade trees and fruite trees. In the middle of the yard was a water pump where the four little girls pumped water for cooking, cleaning and watering the garden. On one side of the yard, Grandma and Grandpa planted tomatoes, beans, squash2(), cucumbers, peppers and strawberries to feed their growing family. They planted roses geraniums lilacs() and irises3 on the other side of the yard,around the statue of the Blessed Mother.

In this case,theσwfvalue is 0.075λwhen SR is 0.8.The RMS value of radial translation errorσΔθis 2.93′when the co-phasing wavelength is 1 μm.

3 Simulation results of cophasing errors

As described above,a number of approximate calculations are required in the theoretical derivation of the co-phasing errors of the primary mirror which need to be implemented in the simulation software.The LOT model was established to simulate real co-phasing errors,which were used for the final reference and comparative verification of the simulation analysis.Nasmyth focal system of LOT will be built as a priority in order to meet the requirements of domestic astronomers for fine detail observation.The parameters are shown in Table 3.

Table 3 Main parameters of Nasmyth focal system of LOT(Focal ratio F number/14.4)

Fig.3 is a ring of petal-shaped segments in the aperture.The inner radius and outer radius of a segment are given byr,Rrespectively and the space between adjacent segments is 2din Fig.3.A1 and B1 represent two vertices of a segment respectively.

Fig.3 Diagram of calculating the aperture of petal-shaped segment

The four vertices of a segment are calculated,and a program was written to generate the segment aperture,which would then be imported into ZEMAX software for the next stage of the simulation.

The fan-shaped segments in each ring can be solved by calculating values for one of the fanshaped segments in a given position,and then rotated to obtain the positions of the other segments in the same ring.For example,the coordinates A1 and B1 of a segment in the second ring are calculated according to

and

The segment aperture program is used to generate the corresponding aperture,which would be imported into ZEMAX software for the next simulation.

The loss of light-gathering energy is 0.41% which is very small compared to the light loss caused by the secondary mirror and the central hole of the M4 mirror in LOT,assuming the space between the fan-shaped segments is 3 mm.It is further calculated that the energy loss values are 1.3% and 2.6% when the gap between fan-shaped segments is 10 mm and 20 mm,respectively.The SR value of LOT is 1 in both cases,and it can therefore be concluded that the gap between the segments has no effect on the image quality except for the shadow effect on the final diffraction pattern.

Each petal-shaped segment has six degrees of freedom and different reflective surface properties,and all segments are independent in the nonsequence mode of ZEMAX software.This method can simulate the actual segmented mirror system as shown in Fig.4.Since three in-plane degrees of freedom of a segment,such as radial displacement and rotation about the center of a segment,can be guaranteed by the machining accuracy,the influence of the three out-of-plane degrees of freedom will be mainly simulated.

Fig.4 3D layout of LOT

3.1 Comprehensive simulation analysis of tip-tilt

The Monte Carlo method is used to analyze the sensitivity of tip-tilt errors.The tilt angles of 72 segments in ZEMAX software obey the same normal distribution.The tolerance width is±2σ(whereσis the standard deviation),where the two ends of the width are equal to the maximum and minimum values of the tolerance operands,and RMS of tilt angles isσif wavefront error RMS is used as the error evaluation standard.If the tilt angles of all segments are uniformly distributed,the tolerance width is±ν,and both ends of the width are equal to the maximum and minimum values of the tolerance operands,the RMS value of tilt angles is given byν/

Firstly,a comprehensive simulation analysis of tip-tilt is carried out.It is assumed that the two types of errors obey the same normal distribution with a tolerance range of±0.032′′.The simulation results are shown in Table 4,and the wavefront of the primary mirror is better than 0.0755λ(λis 1 μm)for 90% of the processes.The calculated SR of telescope will be larger than 0.8 when the tip-tilt RMS value is less than 0.016′′,which is consistent with the theoretically derived value.A simulation result is shown in Fig.5,in which the RMS value of tip-tilt is 0.015′′,the RMS value of the wavefront of LOT is 0.068λ(λis 1 μm)and the SR value of LOT is 0.84.

Table 4 Simulation results of the influence of tip-tilt on SR(subject to normal distribution)

Fig.5 Wavefront of LOT(RMS value of tip-tilt is 0.015′′)

In addition,it is assumed that tip-tilt is uniformly distributed at the same time and the tolerance range is±0.025′′.The simulation results are as follows in Table 5.

Table 5 Simulation results of the influence of tip-tilt on SR(subject to uniform distribution)

The simulation results show that the wavefront of the telescope in 90% of the processes is better than 0.0758λ,indicating that the SR value of the telescope is greater than 0.8 when the RMS value of the uniformly distributed tip-tilt errors is less than 0.0145′′.One of the simulation results is shown in Fig.6,where the RMS value of tip-tilt is 0.013′′,the RMS value of the wavefront of LOT is 0.0588λ(co-phasing wavelengthλis 1 μm)and the SR of LOT is 0.847.

Fig.6 Wavefront of LOT in a simulation of uniformly distributed tip-tilt(RMS value of tip-tilt is 0.013′′)

3.2 Simulation analysis of piston

Monte Carlo analysis for the piston of all segments is undertaken,as the piston of each segment can be assumed to follow a normal distribution.The average value of piston is 0,and the tolerance range is±85 nm.90% of processes will be better than 0.075λ(λis 1 μm),which is taken as the evaluation criterion,and simulation results are shown in Table 6.The RMS value of the wavefront of the telescope in 90% of the processes is less than 0.072λ,indicating that the SR of telescope is greater than 0.8 when RMS of piston is less than 42.5 nm.One of the simulation results is shown in Fig.7,in which the RMS value of piston is 30.5 nm,the RMS value of the wavefront of LOT is 0.0596λ(co-phasing wavelengthλis 1 μm)and the SR of LOT is 0.87.

Fig.7 Wavefront of LOT(RMS value of piston is 30.5 nm)

Table 6 Simulation results of the influence of piston on SR(subject to normal distribution)

In addition,the simulation results indicate the RMS value of the wavefront in 90% of processes in Monte Carlo analysis is superior to 0.075λ,assuming that piston of segments are uniformly distributed with a tolerance range of±71 nm.Simulation results are shown in Table 7.The SR of the telescope is greater than 0.8 when the RMS value of uniformly distributed piston is less than 41 nm.One of the results of the simulation is shown in Fig.8,in which the RMS value of piston is 36.1 nm,the RMS value of the wavefront of LOT is 0.07λ(cophasing wavelengthλis 1 μm)and the SR of LOT is 0.822.

Fig.8 Wavefront of LOT(RMS value of piston is 36.1 nm)

3.3 Comprehensive simulation analysis of tip-tilt and piston

Table 8 shows theoretical values of tip-tilt and piston that are normally distributed with the same weighting.Monte Carlo analysis of the out-ofplane errors is presented below,aiming for 90% of the simulation processes to be better than 0.075λ(λis 1 μm).The piston of all segments is normally distributed and has a range of±56 nm and an average value of 0.In addition,tip-tilt values of all segments follow a normal distribution,with a range of±0.024′′and an average value of 0.The simulation results are shown in Table 9.

Table 7 Simulation results of the influence of piston on SR(subject to uniform distribution)

Table 8 Theoretical values of tip-tilt and piston with equal weighting effect

Table 9 Simulation results of out-of-plane errors(subject to normal distribution,equal weight)

The simulation results show that the RMS value of the wavefront for 90% of processes is better than 0.07235λ.Therefore,the SR of the telescope is greater than 0.8 when the RMS value of tip-tilt with a normal distribution is less than 0.012′′and the RMS value of piston is less than 28 nm.The RMS value of the wavefront of the telescope is 0.06125λand SR=0.867,given one of the integrated simulation processes,seen in Fig.9,where the RMS value of tip-tilt is 0.009′′and the RMS value of piston is 23.4 nm.

Fig.9 Wavefront of LOT(RMS value of tip-tilt is 0.009′′and RMS value of piston is 23.4 nm).

The requirement for the tip-tilt error of segments can be relaxed and the RMS value of piston can be reduced to 20 nm or 15 nm in practical engineering.The comprehensive simulation results are shown in Table 10 and Table 11.

RMS of the wavefront of LOT is better than 0.074λduring 90% of the processes in Table 10,indicating that the SR of the telescope is better than 0.8 when the RMS value of piston is less than 20 nm and the RMS value of tip-tilt is less than 0.014′′.The RMS of the wavefront of LOT is better than 0.0746λduring 90% of the processes in Table 11,indicating that the SR of the telescope is greater than 0.8 when the RMS of piston is less than 15 nm and the RMS of tip-tilt is less than 0.015′′.

Table 10 Monte Carlo simulation results of out-of-plane errors(RMS value of piston is 20 nm and RMS value of tip-tilt is 0.014′′)

Table 11 Monte Carlo simulation results of out-of-plane errors(RMS value of piston is 15 nm and RMS value of tip-tilt is 0.015′′)

4 Conclusion

The Chinese optical infrared astronomy community has strong potential to develop into an international-class community;however,it is challenging to grow this potential without a large telescope like LOT as a platform.Therefore,there is a strong need of a general-purpose optical telescope for China,which could be equipped with different scientific instruments.

We study tiling schemes using fan-shaped segments,and consider the influence of various tiling errors from the segmented primary mirror on image quality,especially the Strehl Ratio,and complete the theoretical derivation for the whole segmented mirror active optics simulation procedure.The segment errors,including three out-plane segment errors,are simulated using both numerical and synthetic approaches,with a lot of simulation results successfully achieved.The SR of LOT is greater than 0.8 when the RMS value of normally distributed tip-tilt errors is less than 0.016′′or the RMS value of normally distributed piston error is less than 42.5 nm with the co-phasing wavelength 1 μm.The comprehensive simulation analysis of the three out-of-plane errors shows that the RMS value of tip-tilt should be less than 0.012′′and the RMS value of piston should be less than 28 nm when tip-tilt and piston are normally distributed with the same weighting.In addition,the RMS value of normally distributed tip-tilt should be less than 0.014′′when the RMS value of normally distributed piston is 20 nm.The RMS value of normally distributed tip-tilt should be less than 0.015′′when the RMS value of normally distributed piston is 15 nm.

The fan-shaped segmented mirror scheme is similar to the hexagonal segmented mirror scheme in terms of co-phasing errors,but the fan-shaped segmented mirror scheme is a better choice based on the processing cost and time.The results can be used as a reference for LOT and will provide LOT with the right to choose the scheme.Our results have important implications for the future development of large diameter optical-infrared telescopes.