First light on an adaptive optics system using a non-modulation pyramid wavefront sensor for a 1.8 m telescope

Similar documents
Sky demonstration of potential for ground layer adaptive optics correction

SOLAR ADAPTIVE OPTICS SYSTEM FOR 1-M NEW VACUUM SOLAR TELESCOPE

Optimized calibration of the adaptive optics system on the LAM Pyramid bench

arxiv: v1 [astro-ph.im] 12 Jul 2018

A STUDY OF PYRAMID WFS BEHAVIOUR UNDER IMPERFECT ILLUMINATION

Shack-Hartmann wavefront sensor sensitivity loss factor estimation in partial correction regime

Adaptive Optics for the Giant Magellan Telescope. Marcos van Dam Flat Wavefronts, Christchurch, New Zealand

Linearity of the pyramid wavefront sensor

Germany; of Arizona, Tucson, AZ USA ABSTRACT

Development of Field of View for Ground-based Optical Telescopes in Adaptive Optics Xiaochun Zhong 1, 2, a, Shujuan Wang 2, b, Zhiliang Huang 3, c

A novel laser guide star: Projected Pupil Plane Pattern

MAORY (Multi conjugate Adaptive Optics RelaY) for E-ELT. Paolo Ciliegi. On behalf of the MAORY Consortium

NA LASER GUIDE STAR AO WITH DYNAMICAL REFOCUS

1. Give short answers to the following questions. a. What limits the size of a corrected field of view in AO?

1. INTRODUCTION ABSTRACT

Analysis of the Sequence Of Phase Correction in Multiconjugate Adaptive Optics

Polarization Shearing Interferometer (PSI) Based Wavefront Sensor for Adaptive Optics Application. A.K.Saxena and J.P.Lancelot

PROCEEDINGS OF SPIE. Large Binocular Telescope Adaptive Optics System: new achievements and perspectives in adaptive optics

Layer Oriented Multi-Conjugate Adaptive Optics systems: Performance analysis by numerical simulations

Measuring tilt and focus for sodium beacon adaptive optics on the Starfire 3.5 meter telescope -- Conference Proceedings (Preprint)

Sky Projected Shack-Hartmann Laser Guide Star

AOL Spring Wavefront Sensing. Figure 1: Principle of operation of the Shack-Hartmann wavefront sensor

Speckles and adaptive optics

SPATIO-TEMPORAL PREDICTION FOR ADAPTIVE OPTICS WAVEFRONT RECONSTRUCTORS

Error Budgets, and Introduction to Class Projects. Lecture 6, ASTR 289

Weston, ACT 2611, Australia; ABSTRACT

Pupil matching of Zernike aberrations

Astronomy 203 practice final examination

Measuring Segment Piston with a Non-Redundant Pupil Mask on the Giant Magellan Telescope

The MAORY Multi-Conjugate Adaptive Optics module Emiliano Diolaiti Istituto Nazionale di Astrofisica

Satellite and debris characterisation in LEO and GEO using adaptive optics

An Introduction to. Adaptive Optics. Presented by. Julian C. Christou Gemini Observatory

Telescopes & Adaptive Optics. Roberto Ragazzoni INAF Astronomical Observatory of Padova

The IPIE Adaptive Optical System Application For LEO Observations

Stray light analysis of an on-axis three-reflection space optical system

International Journal of Scientific & Engineering Research, Volume 4, Issue 7, July ISSN

Adaptive Optics for Satellite and Debris Imaging in LEO and GEO

Wavefront Sensing using Polarization Shearing Interferometer. A report on the work done for my Ph.D. J.P.Lancelot

THE SUBARU CORONAGRAPHIC EXTREME AO HIGH SENSITIVITY VISIBLE WAVEFRONT SENSORS

Astronomie et astrophysique pour physiciens CUSO 2015

Wavefront Correction of Model-based Sensorless Adaptive Optics System

The Flexible Adaptive Optics Concept

Long- and short-term average intensity for multi-gaussian beam with a common axis in turbulence

Control of the Keck and CELT Telescopes. Douglas G. MacMartin Control & Dynamical Systems California Institute of Technology

Double-distance propagation of Gaussian beams passing through a tilted cat-eye optical lens in a turbulent atmosphere

Techniques for direct imaging of exoplanets

Astronomical Seeing. Northeast Astro-Imaging Conference. Dr. Gaston Baudat Innovations Foresight, LLC. April 7 & 8, Innovations Foresight

Wavefront errors due to atmospheric turbulence Claire Max

Astronomy. Optics and Telescopes

Dispersed fringe sensor for the Giant Magellan Telescope

Laboratory Experiments of Laser Tomographic Adaptive Optics at Visible Wavelengths on a 10-meter Telescope

Modelling the multi-conjugate adaptive optics system of the European Extremely Large Telescope

Field Tests of elongated Sodium LGS wave-front sensing for the E-ELT

Ground-Layer Adaptive Optics Christoph Baranec (IfA, U. Hawai`i)

Wavefront Sensing in Astronomy

Adaptive secondary mirror for LBT and its capacitive sensors: how can we calibrate them?

Wavefront reconstruction for adaptive optics. Marcos van Dam and Richard Clare W.M. Keck Observatory

The High Order Test Bench: Evaluating High Contrast Imaging Concepts for SPHERE and EPICS

The measurement of the Interaction Matrices for the LINC-NIRVANA MCAO

Adaptive Optics Overview Phil Hinz What (Good) is Adaptive Optics?

September 9, Wednesday 3. Tools for Solar Observations-I

A Digital Holographic Approach for Co-Phasing of Segmented Telescopes: Proof of Concept Using Numerical Simulations

Wave-front generation of Zernike polynomial modes with a micromachined membrane deformable mirror

Spatial Light Modulator for wavefront correction 1. INTRODUCTION

Deformable mirror fitting error by correcting the segmented wavefronts

Residual phase variance in partial correction: application to the estimate of the light intensity statistics

Introduction to Adaptive Optics. Tim Morris

Designing a Computer Generated Hologram for Testing an Aspheric Surface

arxiv: v1 [astro-ph.im] 4 Jul 2018 ABSTRACT

arxiv:astro-ph/ v1 5 Nov 1999

Sensing and control of segmented mirrors with a pyramid wavefront sensor in the presence of spiders

NB: from now on we concentrate on seeing, as scintillation for large telescopes is unimportant

Layer-oriented adaptive optics for solar telescopes

Exoplanets Direct imaging. Direct method of exoplanet detection. Direct imaging: observational challenges

Atmospheric dispersion correction for the Subaru AO system

McMath-Pierce Adaptive Optics Overview. Christoph Keller National Solar Observatory, Tucson

Telescopes. Astronomy 320 Wednesday, February 14, 2018

Concept, modeling, and performance prediction of a low-cost, large deformable mirror

Atmospheric Turbulence and its Influence on Adaptive Optics. Mike Campbell 23rd March 2009

Point spread function reconstruction at W.M. Keck Observatory : progress and beyond

Analysis of second-harmonic generation microscopy under refractive index mismatch

Why Use a Telescope?

Imaging through Kolmogorov model of atmospheric turbulence for shearing interferometer wavefront sensor

NGAO trade study report: LOWFS architecture KAON 487 (WBS )

Phase Retrieval for the Hubble Space Telescope and other Applications Abstract: Introduction: Theory:

An experimental indoor phasing system based on active optics using dispersed Hartmann sensing technology in the visible waveband

solar telescopes Solar Physics course lecture 5 Feb Frans Snik BBL 707

Present and future Chinese large telescope projects

Exoplanet High Contrast Imaging Technologies Ground

Effect of adaptive telescope mirror dynamics on the residual of atmospheric turbulence correction

Adaptive Optics with Laser Guide Stars - The ALFA system

High contrast imaging at 3-5 microns. Philip M. Hinz University of Arizona Matt Kenworthy, Ari Heinze, John Codona, Roger Angel

More Optical Telescopes

End-to-end model for the Polychromatic Laser Guide Star project (ELP-OA)

대기외란보정을위한단일연속변형거울에관한연구 A Study on a Single Continuous Deformable Mirror for Correction of Atmospheric Disturbance

Measurement of Atmospheric Turbulence with a Shack Hartmann Wavefront Sensor at the new MMT s Prime Focus

High Dynamic Range and the Search for Planets

Telescopes: Portals of Discovery

Chapter 6 Telescopes: Portals of Discovery

Chapter 6 Lecture. The Cosmic Perspective. Telescopes Portals of Discovery Pearson Education, Inc.

Transcription:

First light on an adaptive optics system using a non-modulation pyramid wavefront sensor for a 1.8 m telescope Shengqian Wang ( 王胜千 ) 1,2, **, Kai Wei ( 魏凯 ) 1,2, Wenjia Zheng ( 郑文佳 ) 1,2, and Changhui Rao ( 饶长辉 ) 1,2, * 1 Key Laboratory on Adaptive Optics, Chinese Academy of Sciences, Chengdu 610209, China 2 Institute of Optics and Electronics, Chinese Academy of Sciences, Chengdu 610209, China *Corresponding author: chrao@ioe.ac.cn; **corresponding author: wsq_ioe@126.com Received April 13, 2016; accepted August 19, 2016; posted online September 21, 2016 Our adaptive optics system based on a non-modulation pyramid wavefront sensor is integrated into a 1.8 m astronomical telescope installed at the Yunnan Observatory in LiJiang, and the first light with high-resolution imaging of an astronomical star is successfully achieved. In this Letter, the structure and performance of this system are introduced briefly, and then the observation results of star imaging are reported to show that the angular resolution of an adaptive optics system using a non-modulation pyramid wavefront sensor can approach the diffraction limit quality of a 1.8 m telescope. OCIS codes: 010.1080, 010.7350. doi: 10.3788/COL201614.100101. Optical images of astronomical objects viewed through ground-based telescopes are blurred by the atmosphere [1]. Adaptive optics (AO) attempts to correct, in real time, the distortion induced by the turbulence in the incoming wavefront from the astronomical object of interest, thereby enabling the telescope to reach a diffractionlimited image quality [2]. The AO system uses a wavefront sensor to make measurements of the perturbations introduced by the atmosphere [3], so the wavefront sensor is a key element of the AO system, and its performance has a large influence on the effectiveness of wavefront correction of the AO system [4]. The pyramid wavefront sensor (PWFS) was first proposed by Ragazzoni in 1996 for astronomical AO [5], and it has been widely investigated for applications due to its superiority in adjustable gain and flexible sampling sub-apertures, which makes the PWFS increase the number of accessible scientific targets by more efficiently using guiding star photons. Thus, the PWFS is an attractive option for the next large telescope AO system compared with the Shack Hartmann wavefront sensor [6,7]. The PWFS has been successfully used in the Telescopio Nazionale Galileo telescope and the Large Binocular Telescope, and encouraging on-sky results have been achieved [8,9]. This sensor is also investigated for implementation on the giant Magellan telescope and the Very Large Telescope Optics Facility [10,11]. At present, there are no relevant reports of the PWFS being used in large telescopes in China. In 2014, the PWFS was successfully applied to the 1.8 m telescope at the Yunnan Observatory, and first light on the AO system based on this sensor is reported in this Letter. The 1.8 m telescope at the Yunnan Observatory is a Cassegrain optical structure, and light from a star is reflected by the primary mirror, the secondary mirror, the third mirror, and the other five reflective mirrors in the Coude room. The 127-element AO system installed on the 1.8 m telescope works based on the a Shack Hartmann wavefront sensor with a 13 13 sub-aperture microlens array, and the first light on high-resolution imaging for stars was achieved in 2009 [12]. In 2014, the new AO system, based on the PWFS and installed on the telescope, provided an alternate wavefront sensor for the existing AO system with a Shack Hartmann wavefront sensor. The AO system based on the PWFS was installed in the Coude room (Fig. 1). First, the light rays originating from the natural star reach the Coude focal plane. After propagating through the reflective mirror M1 and the off-axis parabolic mirror PR1_1, the parallel light beam matching the pupil of the tracking system was obtained, and then, through a beam splitter BS1, the light was separated into the tracking system and the reflective mirror M2. The tracking system was used to provide the preliminary correction of the wavefront tilt, which ensures that the beam of light can enter the detection field of the PWFS. With the off-axis parabolic mirrors PR1_2 and PR1_3, the light beam reflected by BS1 was extended to match the pupil of the tilt mirror and the 127-element deformable mirror (DM). Finally, the beam was reduced using the off-axis parabolic mirrors PR2_1 and PR2_2 to match the pupil of the PWFS and the far-field imaging system. The transmission beam from the beam splitter BS2 enters the PWFS, and the reflective beam was transmitted into the far-field imaging system. The detector of the PWFS was an OCAM2 camera, and the corresponding frame rate of the PWFS is 1500 Hz. A Cascade camera was used as the far-field imaging detector, and the imaging band was 700 to 900 nm. 1671-7694/2016/100101(5) 100101-1 2016 Chinese Optics Letters

Fig. 1. Optical layout of the AO system based on the PWFS for 1.8 m telescope. M: reflective mirror, PR: parabolic mirror, TM: tilt mirror, BS: beam splitter mirror. The spectral bandwidth of the imaging system is 700 to 900 nm, and the telescope diameter is 1760 mm. If the center wavelength 800 nm is used as the calculation parameter, the diffraction-limited angular resolution can be calculated according to the following formula: λ D ¼ 0.8 μm 1760 mm ¼ 0.4546 urad ¼ 0.093800 : (1) The optical layout of the PWFS is illustrated in Fig. 2. The light rays through a focusing lens are focused on the top of a refractive pyramid with four sides, which splits the incoming beam into four separate beams, and then a relay lens re-images the four sub-pupils corresponding to each beam onto the detector camera. If the wavefront is not aberrated, the light is equally split by the refractive pyramid and the four sub-pupils are equally and uniformly illuminated. If wavefront is aberrated by turbulence, the intensity difference of these four sub-pupil images contains essential information about the optical aberrations. The intensity distributions of these four sub-pupil images on the detector camera are expressed, respectively, as I 1, I 2, I 3, and I 4, and the detection signals of the PWFS in the x and y directions are calculated with the following formulas [5] : S x ðx; yþ ¼ I 1ðx; yþþi 4 ðx; yþ I 2 ðx; yþ I 3 ðx; yþ I 1 ðx; yþþi 2 ðx; yþþi 3 ðx; yþþi 4 ðx; yþ ; (2) S y ðx; yþ ¼ I 1ðx; yþþi 2 ðx; yþ I 3 ðx; yþ I 4 ðx; yþ I 1 ðx; yþþi 2 ðx; yþþi 3 ðx; yþþi 4 ðx; yþ : (3) Fig. 2. Schematic diagram of the PWFS. When the PWFS is used with modulation by oscillating the optical component, the linear measurement range is increased. Let us consider that the beam is modulated with a circular path with an angular amplitude θ. If the gradient of the aberration is smaller than the modulation θ in the linear measurement range, the PWFS signals relate to the wavefront slope by the formulas [13] W x θ S xðx; yþ; (4) W y θ S yðx; yþ: (5) This linear modulation scheme is used for the calibration of the AO loop closing. In an astronomical AO system like the one presented here, atmospheric turbulence provides for the modulation of the beam around the PWFS [14]. The non-modulation mode is thus selected here to simplify the AO design. If no modulation is applied, the formula between the wavefront ϕðx; yþ and the PWFS signal S x ðx; yþ is given by the following expression (the same formula as the PWFS signal S y ðx; yþ) [15] : S x ðx; yþ ¼ 1 π Z BðyÞ BðyÞ p:v: sin½ϕðx0 ; yþ ϕðx; yþš x 0 dx 0 : (6) x In Eq. (6), the parameter BðyÞ is the x value of the pupil edge at y. Equation (6) shows that the PWFS signal is generated by the phase variations on the entire pupil across the considered direction, and these phase variations contribute to the signal weighted by the inverse of the distance between the integration point ðx 0 ; yþ and the measured point ðx; yþ. The linear relation between the PWFS signal and the phase tilt cannot be obtained using this formula, but we have demonstrated that the correction orientation for the non-modulation PWFS is right. The AO system can be closed with a successive iteration process [16], and the feasibility of this non-modulation method has been demonstrated [17]. Consequently, this AO system 100101-2

for the 1.8 m telescope uses the PWFS without any modulation. In order to analyze the characteristics of the nonmodulation PWFS using in the 1.8 m telescope, a closed-loop AO simulation module is established. The simulation parameters are as follows (consistent with the actual 1.8 m telescope): Telescope diameter: 1760 mm Central obscuration: 280 mm Number of actuators of DM: 127 Number of sub-apertures of PWFS: 13 13 Working mode of PWFS: non-modulation Fried parameter r 0 :135 mm First, the theoretical limit fitting ability of the 127- element DM to the Zernike aberration is analyzed and simulated. The relative error variance (defined as parameter β) is used to describe the correction ability [18]. The parameter β is given by the following expression: β ¼ RMS residual RMS initial : (7) In Eq. (7), RMS initial is the initial RMS error of the wavefront aberration, and RMS residual is the residual RMS error after correction. The smaller the value is, the smaller the residual wavefront error will be, and the stronger the aberration correction ability. The green histograms in Fig. 3 represent the relative error variance for the top 35 orders of Zernike aberrations, that is to say, the correction ability of the AO system using this 127-element DM is not less than the theoretical fitting ability. Second, the closed-loop AO simulation module based on the non-modulation PWFS and the 127-element DM is used to correct the top 35 orders of Zernike aberrations, and the blue histograms in Fig. 3 show the correction effect. Due to the limited space sampling of the PWFS with 13 13 sub-apertures, the AO system cannot achieve the theoretical correction ability of the 127-element DM. From the simulation figure, the PWFS correction result is 0.04 larger than that of the theoretical correction for the 28th mode, and the difference is less than 0.015 for the other modes. Consequently, the 127-element AO system using the non-modulation PWFS can be used, as expected, to correct the top 35 orders of Zernike aberrations. Third, the correction ability for the Kolmogorov phase screen composed of 3 to 104 Zernike modes (according to Fried parameter r 0 ¼ 135 mm) is also simulated with the 127-element AO system based on the non-modulation PWFS (with 13 13 sampling sub-apertures). The initial and corrected far-field images are shown in Fig. 4, and the correction effect reaches the diffraction-limit level (the Strehl ratio is about 0.97 for the corrected far-field spot). The initial and corrected PWFS detection images are illustrated in Fig. 5; when the closed loop is stabilized, the four sub-pupil images have the same intensity distribution (each of the four sub-pupil images has a 13 13 pixels region). The full width at half-maximum (FWHM) for the twodimensional section of the closed-loop spot is illustrated in Fig. 6. If the diffraction-limited FWHM is τ 0, it can be calculated that the x-direction FWHM is 1.02τ 0 and the y-direction FWHM is also 1.02τ 0. The above simulation results verify the feasibility of the non-modulation mode for the PWFS and also lay a good theoretical foundation for the development of the field experiment. The detection CCD of the PWFS for the 1.8 m telescope is an OCAM2 camera with a 240 240 pixelarray. Four sub-pupil images of the PWFS are distributed on the CCD, and each corresponds to the 65 65 pixels region (see Fig 7(a)). With the 5 5 pixel-rebin processing, the sampling of each sub-pupil changs to 13 13 pixels (see Fig. 7(b)), which is equivalent to the original AO system based on a Shack Hartmann wavefront sensor Fig. 4. Simulation results of the initial and corrected far-field images for the turbulence phase screen (D ¼ 1760 mm, r 0 ¼ 135 mm). Fig. 3. Relative error variance of the top 35 orders of Zernike aberrations. Fig. 5. Simulation results of the initial and corrected PWFS images for the turbulence phase screen (D ¼ 1760 mm, r 0 ¼ 135 mm). 100101-3

Fig. 9. Space arrangement of the 127 actuators (the outer circle represents the mirror edge of the DM, and the gray area represents 1.8 m telescope pupil plane). Fig. 6. FWHM of the diffraction-limited spot, and the twodimensional FWHM of the simulated closed-loop spot. Fig. 10. Open-loop and closed-loop star images (Name: HIP113881, Magnitude: 2.44). Fig. 7. (a) Original image (each sub-pupil has 65 65 pixels) and (b) rebin image (each sub-pupil has 13 13 pixels). In late October 2014, we performed the final integration, alignment, and calibration of the AO system based on the PWFS, and high-resolution imaging for the first light of the natural stars was obtained. The open-loop and closed-loop star images of the AO system are shown in Fig. 10. The name of this imaging stellar star was HIP113881, and the magnitude of this star was with 13 13 sampling sub-apertures. Figure 7 is a calibration image that is used to determine the position of the four PWFS sub-pupil images, and the following Fig. 8 illustrates the open-loop and closed-loop images for the actual field experiment. The optical path length aberrations of a beam are corrected by the mechanical displacement of a DM s surface. There are three main DM families: discrete actuator mirrors using a continuous faceplate; bimorph mirrors, in which the actuator is combined with the faceplate; and segmented mirrors [2]. The DM used in the AO system of the 1.8 m telescope is a 127-discrete actuators mirror, and the arrangement of these 127 actuators in the DM is shown in Fig. 9. Fig. 8. Open-loop and closed-loop experiment PWFS images. Fig. 11. Two-dimensional section of the closed-loop image. 100101-4

2.44. The open-loop imaging spot was scattered on the detection surface of the camera due to the atmospheric turbulence, and after closed-loop correction with the AO system, the spot energy was concentrated. The FWHM of the imaging spot is a criterion method of the actual angular resolution. In order to make a quantitative study on the imaging performance of the AO system, the FWHM for the two-dimensional section of the closed-loop star image is illustrated in Fig. 11. It can be calculated that the x-direction FWHM is 0.1109 00, and the y-direction FWHM is 0.1015 00. The x-direction FWHM is about 1.18 times the diffractionlimited value (0.1109 0.0938 1.18), and the y-direction FWHM is about 1.08 times the diffraction-limited value (0.1015 0.0938 1.08). Compared with the above simulation results, the x-direction FWHM reaches about 1.16 times the simulation value (1.18 1.02 1.16), and the y-direction FWHM reaches about 1.06 times the simulation value (1.08 1.02 1.06). In conclusion, the first light of an AO system based on a non-modulation PWFS on a 1.8 m telescope is reported. The results show that the AO system using the nonmodulation PWFS can work satisfactorily and the diffraction limit performance can be obtained. A significant advantage of the PWFS is that the number of subapertures can be adjusted according to the change of the turbulence, so that the PWFS has better adaptability compared with a Shack Hartmann wavefront sensor. The corresponding experiment has been carried out and will be reported in the future. I would like to acknowledge and thank Professor Wenhan Jiang for his guidance and support in the course of this study. This work was supported by the National Natural Science Foundation of China under Grant No. 61008038. References 1. R. J. Noll, J. Opt. Soc. Am. 66, 207 (1976). 2. J. W. Hardy, Adaptive Optics for Astronomical Telescopes (Oxford University, 1998). 3. W. H. Southwell, J. Opt. Soc. Am. 70, 998 (1980). 4. L. Yang, L. Hu, D. Li, Z. Cao, Q. Mu, J. Ma, and L. Xuan, Chin. Opt. Lett. 13, 120801 (2015). 5. R. Ragazzoni, J. Mod. Opt. 43, 289 (1996). 6. R. Ragazzoni and J. Farinato, Astron Astrophys. 350, 23 (1999). 7. T. Y. Chew, R. M. Clare, and R. G. Lane, Opt. Commun. 268, 189 (2006). 8. R. Ragazzoni, A. Ghedina, A. Baruffolo, E. Marchetti, J. Farinato, T. Niero, G. Crimi, and M. Ghigo, Proc. SPIE 4007, 423 (2000). 9. S. Esposito, A. Riccardi, L. Fini, A. T. Puglisi, E. Pinna, M. Xompero, R. Briguglio, F. Quirós-Pacheco, P. Stefanini, J. C. Guerra, L. Busoni, A. Tozzi, F. Pieralli, G. Agapito, G. Brusa- Zappellini, R. Demers, J. Brynnel, C. Arcidiacono, and P. Salinari, Proc. SPIE 7736, 773609 (2010). 10. A. H. Bouche, D. S. Acton, R. Biasi, R. Conan, B. Espeland, S. Esposito, J. Filgueira, D. Gallieni, B. A. McLeod, E. Pinna, F. Santoro, G. Trancho, and M. A. Van Dam, Proc. SPIE 9148, 91480W (2014). 11. F. Quirós-Pacheco, G. Agapito, A. Riccardi, S. Esposito, M. Le Louarn, and E. Marchetti, Proc. SPIE 8447, 84475L (2012). 12. K. Wei, X. Zhang, H. Xian, W. Ma, A. Zhang, L. Zhou, C. Guan, M. Li, D. Chen, S. Chen, Z. Liao, C. Rao, and Y. Zhang, Chin. Opt. Lett. 8, 1019 (2010). 13. A. Burvall, E. Daly, S. R. Chamot, and C. Dainty, Opt. Express 14, 11925 (2006). 14. J. B. Costa, Appl. Opt. 44, 60 (2005). 15. C. Vérinaud, Opt. Commun. 233, 27 (2004). 16. S. Wang, C. Rao, H. Xian, J. Zhang, J. Wang, and Z. Liu, Opt. Express 19, 8135 (2011). 17. S. Wang, C. Rao, A. Zhang, X. Zhang, K. Wei, Y. Tian, Z. Liao, C. Zhang, H. Xian, X. Zhang, and L. Wei, Proc. SPIE 8447, 84476K (2012). 18. H. Duan, E. Li, H. Wang, Z. Yang, and Y. Zhang, Acta Opt. Sin. 23, 1143 (2003). 100101-5