Introduction

This simulation investigates the influence of the FFTS performance on the quality of the reconstructions. This means that the AO will give a strehl of 0.30 and the FFTS will give a residual OPD error.

The setup of the simulation is described in table 1. The input data for the test uses a strehl value of 0.30 for the target and the calibrator which has a magnitude of 14. In addition the target and calibrator images are influenced by a gaussian distributed phase error with the same standard deviation, resulting in a fringe contrast loss. This error is modelled as a phase error on the second pupil and for each phase screen (the properties of the atmosphere will result in about 600 phase screens per position angle) a random phase error is introduced. This phase error is specified as a gaussian error with a standard deviation of the OPD between 0.0 (perfect FFTS) and 0.5 lambda (worst case, no FFTS available).

ParameterDescription
Atmosphere/AOThe AO delivers a strehl of about 30 percent for the target and the psf star. The psfs are separately simulated and therefore not equal!
FFTSThe FFTS will leave a residual phase error which is given as a standard deviation of 0.0, 0.05, 0.1, 0.2, and 0.5 lambda.
ObservationThe observation target is NGC4151 with an overlayed dust tori (simulated) at 10.1mag integral brightness. The brightness of the calibrator is set to 14mag. For each position angle (108, 144, 180, 216, 252 degree) an image equivalent to 60s exposure time was created. The wavelength was 1.2 micron (whole J-Band) and 2.2 micron (whole K-band).
Table 1: Setup for the simulation in J-Band and K-Band, for the common test setup see table 1.

The raw data were generated according to the common scheme described in section Overview. In addition the phase screens for the psfs are influenced by a phase error on the second pupil, and the images for the target and psf-star are generated separately which means, that the calibrator image is not influenced by the target (no halo, etc.).

All simulated input data for the LN DRS pipeline are available as a tar-file (ex3_j_input.tar.gz (18MB) and ex3_k_input.tar.gz (20MB)). The corresponding results are also available as tar files (ex3_j_results.tar.gz (2.7MB) and ex3_k_results.tar.gz (2.8MB)).

The generated calibrator raw images are shown in section Calibrator. The raw images and deconvolution results are shown in section Raw Frames and Results for an AGN for the AGN and section Raw Frames and Results for a Star Cluster for the star cluster.

Calibrator

The simulated raw LBT interferograms used for the deconvolution are ideal images, they are not influenced by detector effects like different pixel gain or bad pixels. In figure 1 the calibrators for the different phase errors and one position angle are shown.

J-Band, calibrator, phase error 0.0J-Band, calibrator, phase error 0.1J-Band, calibrator, phase error 0.2J-Band, calibrator, phase error 0.5
K-Band, calibrator, phase error 0.0K-Band, calibrator, phase error 0.1K-Band, calibrator, phase error 0.2K-Band, calibrator, phase error 0.5
Figure 1: Central 128x128 pixels of the generated J-Band (top row) and K-Band (bottpon row) images of a calibrator with phase errors of (from top left to bottom right) 0.0, 0.1, 0.2, and 0.5 lambda.

Raw Frames and Results for an AGN

The basis of the simulation is an image of NGC4151 with an overlayed dust torus (see figure 1). The simulated raw images for a position angle of 108 degree are shown in figure 2.

J-Band raw image for a phase error of 0.0J-Band raw image for a phase error of 0.1J-Band raw image for a phase error of 0.2J-Band raw image for a phase error of 0.5
K-Band raw image for a phase error of 0.0K-Band raw image for a phase error of 0.1K-Band raw image for a phase error of 0.2K-Band raw image for a phase error of 0.5
Figure 2: The simulated J-Band (top row) and K-Band (bottom row) raw images for a position angle of 108 degree of NGC4151 at 10.1 mag including sky background with a phase error of 0.0, 0.1, 0.2, and 0.5 lambda (from left to right).

A comparison of the reconstructions depending on the phase error is presented in figure 3 for the J-Band. For the K-Band the results are presented in figure 4.

J-Band coadded raw image, phase error 0.0J-Band coadded raw image, phase error 0.1J-Band coadded raw image, phase error 0.2J-Band coadded raw image, phase error 0.3
J-Band Richardson-Lucy reconstruction, phase error 0.0J-Band Richardson-Lucy reconstruction, phase error 0.1J-Band Richardson-Lucy reconstruction, phase error 0.2J-Band Richardson-Lucy reconstruction, phase error 0.5
J-Band Building-Block reconstruction, phase error 0.0J-Band Building-Block reconstruction, phase error 0.1J-Band Building-Block reconstruction, phase error 0.2J-Band Building-Block reconstruction, phase error 0.5
Figure 3: The top row shows the coadded J-Band raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The phase error is 0.0, 0.1, 0.2, and 0.5 lambda (from left to right). In the third row, the Building-Block algorithm was used
K-Band coadded raw image, phase error 0.0K-Band coadded raw image, phase error 0.1K-Band coadded raw image, phase error 0.2K-Band coadded raw image, phase error 0.3
K-Band Richardson-Lucy reconstruction, phase error 0.0K-Band Richardson-Lucy reconstruction, phase error 0.1K-Band Richardson-Lucy reconstruction, phase error 0.2K-Band Richardson-Lucy reconstruction, phase error 0.5
K-Band Building-Block reconstruction, phase error 0.0K-Band Building-Block reconstruction, phase error 0.1K-Band Building-Block reconstruction, phase error 0.2K-Band Building-Block reconstruction, phase error 0.5
Figure 4: The top row shows the coadded K-Band raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The phase error is 0.0, 0.1, 0.2, and 0.5 lambda (from left to right). In the third row, the Building-Block algorithm was used.

In table 2 (J-Band) and table 3 (K-Band) the errors for the star cluster depending on the OPD error are presented. In figure 5 the image errors depending on the phase error jitter are shown.

TestcaseRichardson-LucyBuilding-Block
Phase ErrorImage ErrorImage Error
0.0 0.00.012 (0.012)0.007 (0.007)
0.0 0.050.013 (0.013)0.008 (0.008)
0.0 0.10.016 (0.016)0.008 (0.008)
0.0 0.20.031 (0.031)0.012 (0.012)
0.0 0.30.062 (0.062)0.019 (0.019)
0.0 0.40.081 (0.081)0.036 (0.036)
0.0 0.50.088 (0.088)0.038 (0.038)
Table 2: AGN: Errors depending on the OPD error (J-Band).
TestcaseRichardson-LucyBuilding-Block
Phase ErrorImage ErrorImage Error
0.0 0.00.020 (0.020)0.011 (0.011)
0.0 0.050.027 (0.027)0.016 (0.016)
0.0 0.10.028 (0.028)0.014 (0.014)
0.0 0.20.053 (0.053)0.022 (0.022)
0.0 0.30.108 (0.108)0.035 (0.035)
0.0 0.40.145 (0.144)0.057 (0.057)
0.0 0.50.149 (0.147)0.082 (0.082)
Table 3: AGN: Errors depending on the OPD error (K-Band).
Image errors depending on the phase error jitter
Figure 5: In this figure, the image errors depending on the phase error jitter are shown.

In figure 6 the image error describing the difference between the reconstruction and the reference image is shown for each iteration (in steps of 100 iterations). The effect of the OPD error on the reconstruction error and the optimal iteration number is clearly visible.

Richardson-Lucy reconstruction error (J-Band)Richardson-Lucy reconstruction error (K-Band)
Building-Block reconstruction error (J-Band)Building-Block reconstruction error (K-Band)
Figure 6: Reconstruction error depending on the iteration number for the J-Band (left side) and K-Band (right side). For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.

In figure 7 (J-Band) and figure 8 (K-Band) a horizontal cut slightly above the intensity maximum of the reconstructions of NHC4151 compared with the ideal reference image is shown.

Richardson-Lucy: Profile through the reconstructed galaxy (J-Band)Richardson-Lucy: Central part of the profile through the reconstructed galaxy (J-Band)
Building-Block: Profile through the reconstructed galaxy (J-Band)Building-Block: Central part of the profile through the reconstructed galaxy (J-Band)
Figure 7: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (J-Band). On the right side only the central part is shown. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.
Richardson-Lucy: Profile through the reconstructed galaxy (K-Band)Richardson-Lucy: Central part of the profile through the reconstructed galaxy (K-Band)
Building-Block: Profile through the reconstructed galaxy (K-Band)Building-Block: Central part of the profile through the reconstructed galaxy (K-Band)
Figure 8: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (K-Band). On the right side only the central part is shown. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.

Raw Frames and Results for a Star Cluster

The basis of the simulation is an image of a star cluster (see figure 2). The simulated raw images for a position angle of 108 degree are shown in figure 9.

J-Band raw image for a phase error of 0.0J-Band raw image for a phase error of 0.1J-Band raw image for a phase error of 0.2J-Band raw image for a phase error of 0.5
K-Band raw image for a phase error of 0.0K-Band raw image for a phase error of 0.1K-Band raw image for a phase error of 0.2K-Band raw image for a phase error of 0.5
Figure 9: The simulated J-Band (top row) and K-Band (bottom row) raw images for a position angle of 108 degree of the star cluster at 10.1 mag including sky background with a phase error of 0.0, 0.1, 0.2, and 0.5 lambda (from left to right).

A comparison of the reconstructions depending on the phase error is presented in figure 10 for the J-Band. For the K-Band the results are presented in figure 11.

J-Band coadded raw image, phase error 0.0J-Band coadded raw image, phase error 0.1J-Band coadded raw image, phase error 0.2J-Band coadded raw image, phase error 0.3
J-Band Richardson-Lucy reconstruction, phase error 0.0J-Band Richardson-Lucy reconstruction, phase error 0.1J-Band Richardson-Lucy reconstruction, phase error 0.2J-Band Richardson-Lucy reconstruction, phase error 0.5
J-Band Building-Block reconstruction, phase error 0.0J-Band Building-Block reconstruction, phase error 0.1J-Band Building-Block reconstruction, phase error 0.2J-Band Building-Block reconstruction, phase error 0.5
Figure 10: The top row shows the coadded J-Band raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The phase error is 0.0, 0.1, 0.2, and 0.5 lambda (from left to right). In the third row, the Building-Block algorithm was used
K-Band coadded raw image, phase error 0.0K-Band coadded raw image, phase error 0.1K-Band coadded raw image, phase error 0.2K-Band coadded raw image, phase error 0.3
K-Band Richardson-Lucy reconstruction, phase error 0.0K-Band Richardson-Lucy reconstruction, phase error 0.1K-Band Richardson-Lucy reconstruction, phase error 0.2K-Band Richardson-Lucy reconstruction, phase error 0.5
K-Band Building-Block reconstruction, phase error 0.0K-Band Building-Block reconstruction, phase error 0.1K-Band Building-Block reconstruction, phase error 0.2K-Band Building-Block reconstruction, phase error 0.5
Figure 11: The top row shows the coadded K-Band raw images. The second row show the reconstructions using the Richardson-Lucy algorithm. The phase error is 0.0, 0.1, 0.2, and 0.5 lambda (from left to right). In the third row, the Building-Block algorithm was used

In table 4 (J-Band) and table 5 (K-Band) the errors for the star cluster depending on the OPD error are presented. The image errors are also summarized in figure 12 and the photometric errors in figure 13.

TestcaseRichardson-LucyBuilding-Block
Phase ErrorImage ErrorMagnitude ErrorPosition ErrorImage ErrorMagnitude ErrorPosition Error
0.0 0.00.013 (0.013)0.02, 0.01, 0.02, 0.19, 0.000.00, 0.00, 0.07, 0.08, 0.210.013 (0.013)0.02, 0.00, 0.10, 0.29, 0.380.00, 0.00, 0.05, 0.01, 0.04
0.0 0.050.016 (0.016)0.01, 0.00, 0.04, 0.06, 0.090.01, 0.02, 0.27, 0.05, 0.020.015 (0.015)0.01, -0.00, 0.16, 0.17, 0.440.01, 0.01, 0.09, 0.01, 0.05
0.0 0.10.014 (0.014)0.01, 0.01, 0.02, 0.11, 0.340.01, 0.01, 0.21, 0.06, 0.070.012 (0.012)0.01, 0.00, 0.13, 0.22, 0.850.00, 0.00, 0.07, 0.01, 0.10
0.0 0.20.024 (0.024)0.01, 0.01, -0.00, 0.10, 0.240.00, 0.00, 0.17, 0.14, 0.040.020 (0.020)0.01, 0.00, 0.06, 0.18, 0.590.00, 0.01, 0.06, 0.10, 0.06
0.0 0.30.055 (0.055)0.01, 0.00, -0.00, 0.04, 0.190.01, 0.00, 0.26, 0.17, 0.490.025 (0.025)0.01, 0.00, 0.09, 0.17, 0.600.00, 0.00, 0.23, 0.08, 0.35
0.0 0.40.098 (0.098)0.01, 0.01, 0.04, 0.07, 0.280.01, 0.01, 0.15, 0.20, 0.770.036 (0.036)-0.00, 0.01, 0.06, 0.11, 0.430.01, 0.02, 0.12, 0.17, 0.77
0.0 0.50.112 (0.112)0.01, 0.01, -0.01, 0.05, 0.030.01, 0.02, 0.06, 0.02, 0.560.054 (0.053)0.01, 0.01, 0.04, 0.12, 0.230.01, 0.02, 0.09, 0.00, 0.56
Table 4: Star cluster: Errors depending on the OPD error (J-Band).
TestcaseRichardson-LucyBuilding-Block
Phase ErrorImage ErrorMagnitude ErrorPosition ErrorImage ErrorMagnitude ErrorPosition Error
0.0 0.00.043 (0.043)0.01, 0.00, 0.08, 0.20, 1.850.07, 0.01, 0.69, 1.52, 1.180.029 (0.029)0.02, 0.01, 0.13, 0.29, 1.840.05, 0.00, 0.77, 1.55, 0.77
0.0 0.050.057 (0.057)0.05, 0.02, -0.04, 0.43, 2.050.08, 0.03, 1.16, 0.16, 0.480.042 (0.042)0.06, 0.02, 0.00, 0.52, 1.920.08, 0.03, 1.11, 0.11, 0.48
0.0 0.10.055 (0.055)0.04, 0.01, 0.31, 0.57, 1.590.05, 0.03, 0.75, 1.47, 0.530.034 (0.034)0.04, 0.01, 0.36, 0.65, 1.730.03, 0.04, 0.51, 1.59, 0.17
0.0 0.20.110 (0.110)0.04, 0.01, 0.29, 0.32, 1.980.07, 0.11, 0.26, 0.95, 0.680.065 (0.065)0.03, 0.00, 0.28, 0.40, 1.910.07, 0.10, 0.28, 0.94, 0.57
0.0 0.30.155 (0.155)0.08, 0.02, 0.56, 0.90, 2.080.12, 0.05, 1.77, 1.08, 0.400.070 (0.070)0.04, -0.01, 0.47, 1.02, 1.930.12, 0.04, 2.24, 1.58, 0.44
0.0 0.40.186 (0.186)0.09, 0.05, 0.28, 0.08, 2.090.09, 0.02, 1.20, 0.86, 0.400.128 (0.128)0.04, 0.02, 0.21, 0.08, 1.370.10, 0.03, 1.37, 0.90, 1.33
0.0 0.50.204 (0.204)0.11, 0.05, 0.34, 0.54, 0.220.06, 0.02, 2.20, 0.95, 1.530.124 (0.124)0.05, 0.01, 0.31, 0.57, 0.310.07, 0.02, 2.77, 1.00, 1.61
Table 5: Star cluster: Errors depending on the OPD error(K-Band).
Image errors depending on the phase error jitter
Figure 12: In this figure, the image errors depending on the phase error jitter are shown.
Photometric errors depending on the phase error jitter using Richardson-Lucy (J-Band)Photometric errors depending on the phase error jitter using Richardson-Lucy (K-Band)
Photometric errors depending on the phase error jitter using Building-Block (J-Band)Photometric errors depending on the phase error jitter using Building-Block (K-Band)
Figure 13: Photometric errors depending on the phase error jitter in the J-Band (left side) and K-Band (right side). The reconstruction in the first row uses the Richardson-Lucy algorithm, in the second row, the Building-Block method was used.

In figure 14 the image error describing the difference between the reconstruction and the reference image is shown for each iteration (in steps of 100 iterations). The effect of the OPD error on the reconstruction error and the optimal iteration number is clearly visible.

Richardson-Lucy reconstruction error (J-Band)Richardson-Lucy reconstruction error (K-Band)
Building-Block reconstruction error (J-Band)Building-Block reconstruction error (K-Band)
Figure 14: Reconstruction error depending on the iteration number for the J-Band (left side) and K-Band (right side). For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.

In the figures figure 15, figure 16, figure 17, and figure 18 the image and photometric errors depending on the calibrator strehl and the iteration number is shown. The overall shape is the same, but the optimum number of iterations is difficult to select, because not all photometric errors have their minimum at the same iteration number.

Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, phase error 0.0)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, phase error 0.0)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, phase error 0.0)Building-Block: Image and photometric errors depending on the iteration number (K-Band, phase error 0.0)
Figure 15: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The phase error is 0.0. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, phase error 0.1)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, phase error 0.1)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, phase error 0.1)Building-Block: Image and photometric errors depending on the iteration number (K-Band, phase error 0.1)
Figure 16: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The phase error is 0.1. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, phase error 0.2)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, phase error 0.2)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, phase error 0.2)Building-Block: Image and photometric errors depending on the iteration number (K-Band, phase error 0.2)
Figure 17: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The phase error is 0.2. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, phase error 0.5)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, phase error 0.5)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, phase error 0.5)Building-Block: Image and photometric errors depending on the iteration number (K-Band, phase error 0.5)
Figure 18: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The phase error is 0.5. For the reconstructions in the top row, Richardson-Lucy was used. In the second row, the Building-Block algorithm was used.