Introduction

This test investigates the dependency of the reconstruction error on the atmospheric dispersion (refraction). In addition, a very red target (0.5 mag brighter in the red) and a slightly red calibrator (0.2 mag brighter in the red) were used. For this test, the PSFs are generated using several (41) monochromatic PSFs. This simulation only tests the difficult cases where the target has a declination between -30 and +20 degree. For each target, three parallactic angles were used which corresponds to the largest angle which can be achieved going down to 75 degree zenith angle.

The setup of the simulation is described in table 1. The input data for the test uses a strehl value of 0.30 for target and calibrator. The psf star has a magnitude of 14. In addition, the target and psf star PSF are generated by using 11 (no atmospheric dispersion) or 41 (with atmospheric dispersion) monochromatic PSFs.

The multi-monochromatic PSFs are created by generating a monochromatic OPD screen for the shortest wavelength λ0. The OPD values are then rescaled for a given wavelength λ (we used 11 or 41 discrete wavelength to sample a spectral band) by λ0/λ. The result was put into an array which was enlarged by λ/λ0. Then the atmospheric dispersion was simulated by moving the center of the psf according to the assumed zenith angle and a model atmosphere. A PSF for a specific wavelength was then calculated by |FFT-1|2, weighted by the spectral intensity and all summed up to the final multi-monochromatic PSF.

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!. The model of the atmosphere is from Allan's Astrophysical Quantities, 1999, section 11.20 Atmospheric refraction and air path, page 262.
SpectraThe target and calibrator show a spectra which was sampled at 41 equally spaced wavelength over a part of the J-Band (1.156 - 1.250 μm: J shorter half) and K-Band (2.0 - 2.4 μm). The spectrum of the target is 0.5 mag brighter on the red side and the spectrum of the calibrator is 0.2 mag brighter at the red side.
ObservationThe observation target is NGC4151 with an overlayed dust tori (simulated) at 10.1mag integral brightness and a small star cluster. The brightness of the psf-star is set to 14mag. For each position angle an image equivalent to 60s exposure time was created. The central wavelength was 1.203 micron (shorter side of the J-Band) and 2.2 micron (whole K-band). This target was placed at -30, -20, -10, 0, 10, and 20 degree declination.
Table 1: Setup for the simulation in J-Band and K-Band, for the common test setup see table 1.

In table 2 all combinations of declination, culmination, smallest altitude, and parallactic angles are shown. In addition, each test case was also created without atmospheric dispersion.

Declination [deg]Culmination [deg]Smallest altitude [deg]PA [deg]
-3027.315.6144 180 216
-2037.316.5135 180 225
-1047.315.0128 180 232
057.315.1124 180 236
1067.328.9122 180 238
2077.355.7120 180 240
Table 2: Combinations of declination, culmination, smallest altitude, and achieved parallactic angles for target and calibrator.

The raw data were generated according to the common scheme described in section Overview. In addition, the target and psf star PSF are generated by using 41 monochromatic PSFs.

All simulated input data for the LN DRS pipeline are available as a tar-file (ex8_j_input.tar.gz (3.9MB) and ex8_k_input.tar.gz (4.3MB)). The corresponding results are also available as tar files (ex8_j_results.tar.gz (552KB) and ex8_k_results.tar.gz (546KB)).

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 and figure 2 the raw calibrator for different declination is shown.

J-Band, calibrator, decl -30, angle 144J-Band, calibrator, decl -20, angle 135J-Band, calibrator, decl 0, angle 124J-Band, calibrator, decl +20, angle 120
J-Band, calibrator, decl -30, angle 180J-Band, calibrator, decl -20, angle 180J-Band, calibrator, decl 0, angle 180J-Band, calibrator, decl +20, angle 180
Figure 1: Central 128x128 pixels of the generated J-Band calibrator images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the calibrator images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the calibrator images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).
K-Band, calibrator, decl -30, angle 144K-Band, calibrator, decl -20, angle 135K-Band, calibrator, decl 0, angle 124K-Band, calibrator, decl +20, angle 120
K-Band, calibrator, decl -30, angle 180K-Band, calibrator, decl -20, angle 180K-Band, calibrator, decl 0, angle 180K-Band, calibrator, decl +20, angle 180
Figure 2: Central 128x128 pixels of the generated K-Band calibrator images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the calibrator images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the calibrator images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).

Raw Frames and Results for an AGN

The basis of this simulation is an image of NGC4151 with an overlayed dust torus (see figure 1). The simulated raw images are shown in figure 3 (J-Band) and in figure 4 (K-Band).

J-Band, NGC4151, decl -30, angle 144J-Band, NGC4151, decl -20, angle 135J-Band, NGC4151, decl 0, angle 124J-Band, NGC4151, decl +20, angle 120
J-Band, NGC4151, decl -30, angle 180J-Band, NGC4151, decl -20, angle 180J-Band, NGC4151, decl 0, angle 180J-Band, NGC4151, decl +20, angle 180
Figure 3: Central 128x128 pixels of the generated J-Band AGN images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the target images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the AGN images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).
K-Band, NGC4151, decl -30, angle 144K-Band, NGC4151, decl -20, angle 135K-Band, NGC4151, decl 0, angle 124K-Band, NGC4151, decl +20, angle 120
K-Band, NGC4151, decl -30, angle 180K-Band, NGC4151, decl -20, angle 180K-Band, NGC4151, decl 0, angle 180K-Band, NGC4151, decl +20, angle 180
Figure 4: Central 128x128 pixels of the generated K-Band AGN images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the target images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the AGN images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).

A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 5 and figure 6 for the J-Band and K-Band.

J-Band, coadded, -30 degree declinationJ-Band, coadded, -20 degree declinationJ-Band, coadded, 0 degree declinationJ-Band, coadded, +20 degree declination
J-Band, Richardson-Lucy reconstruction, -30 degree declinationJ-Band, Richardson-Lucy reconstruction, -20 degree declinationJ-Band, Richardson-Lucy reconstruction, 0 degree declinationJ-Band, Richardson-Lucy reconstruction, +20 degree declination
J-Band, Buidling-Block reconstruction, -30 degree declinationJ-Band, Buidling-Block reconstruction, -20 degree declinationJ-Band, Buidling-Block reconstruction, 0 degree declinationJ-Band, Buidling-Block reconstruction, +20 degree declination
Figure 5: Central 128x128 pixels of the coadded and reconstructed J-Band NGC4151 images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the coadded raw target images. The second row contain the reconstructions created with the Richardson-Lucy algorithm. For the third row, the Building-Block method was used.
K-Band, coadded, -30 degree declinationK-Band, coadded, -20 degree declinationK-Band, coadded, 0 degree declinationK-Band, coadded, +20 degree declination
K-Band, Richardson-Lucy reconstruction, -30 degree declinationK-Band, Richardson-Lucy reconstruction, -20 degree declinationK-Band, Richardson-Lucy reconstruction, 0 degree declinationK-Band, Richardson-Lucy reconstruction, +20 degree declination
K-Band, Buidling-Block reconstruction, -30 degree declinationK-Band, Buidling-Block reconstruction, -20 degree declinationK-Band, Buidling-Block reconstruction, 0 degree declinationK-Band, Buidling-Block reconstruction, +20 degree declination
Figure 6: Central 128x128 pixels of the coadded and reconstructed K-Band NGC4151 images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the coadded raw target images. The second row contain the reconstructions created with the Richardson-Lucy algorithm. For the third row, the Building-Block method was used.

In figure 7 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 atmospheric dispersion on the reconstruction error and the optimal iteration number is clearly visible.

Richardson-Lucy reconstruction error (J-Band)Richardson-Lucy reconstruction error (K-Band)
Buidling-Block reconstruction error (J-Band)Buidling-Block reconstruction error (K-Band)
Figure 7: Reconstruction error depending on the iteration number for the J-Band (left side) and K-Band (right side). The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.

In figure 8 (J-Band) and figure 9 (K-Band) a horizontal cut slightly above the intensity maximum of the reconstructions of NGC4151 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 8: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (J-Band). On the left side only the central part is shown. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method 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 9: Horizontal profile through the center of the reconstructed galaxy compared to the ideal image (K-Band). On the left side only the central part is shown. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.

In table 3 (J-Band) and table 4 (K-Band) the errors for the AGN depending on the atmospheric dispersion are presented.

TestcaseRichardson-LucyBuilding-Block
DecPAAltitudesImage ErrorImage Error
-30144 180 21615.6 27.3 15.60.156 (0.155)0.120 (0.120)
-20135 180 22516.5 37.3 16.50.143 (0.142)0.103 (0.103)
-10128 180 23215.0 47.3 15.00.124 (0.123)0.088 (0.088)
0124 180 23615.1 57.3 15.10.091 (0.091)0.059 (0.059)
10122 180 23828.9 67.3 28.90.068 (0.068)0.047 (0.047)
20120 180 24055.7 77.3 55.70.037 (0.037)0.021 (0.021)
Table 3: NGC4151: Errors depending on the atmospheric dispersion (short side of the J-Band).
TestcaseRichardson-LucyBuilding-Block
DecPAAltitudesImage ErrorImage Error
-30144 180 21615.6 27.3 15.60.157 (0.157)0.127 (0.127)
-20135 180 22516.5 37.3 16.50.118 (0.117)0.079 (0.079)
-10128 180 23215.0 47.3 15.00.095 (0.094)0.046 (0.046)
0124 180 23615.1 57.3 15.10.075 (0.075)0.041 (0.041)
10122 180 23828.9 67.3 28.90.044 (0.044)0.022 (0.022)
20120 180 24055.7 77.3 55.70.038 (0.038)0.017 (0.017)
Table 4: NGC4151: Errors depending on the atmospheric dispersion (K-Band).

Raw Frames and Results for a Star Cluster

The basis of this simulation is a small star cluster (see figure 2). The simulated raw images are shown in figure 10 (J-Band) and in figure 11 (K-Band).

J-Band, starcluster, decl -30, angle 144J-Band, starcluster, decl -20, angle 135J-Band, starcluster, decl 0, angle 124J-Band, starcluster, decl +20, angle 120
J-Band, starcluster, decl -30, angle 180J-Band, starcluster, decl -20, angle 180J-Band, starcluster, decl 0, angle 180J-Band, starcluster, decl +20, angle 180
Figure 10: Central 128x128 pixels of the generated J-Band star cluster images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the calibrator images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the star cluster images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).
K-Band, starcluster, decl -30, angle 144K-Band, starcluster, decl -20, angle 135K-Band, starcluster, decl 0, angle 124K-Band, starcluster, decl +20, angle 120
K-Band, starcluster, decl -30, angle 180K-Band, starcluster, decl -20, angle 180K-Band, starcluster, decl 0, angle 180K-Band, starcluster, decl +20, angle 180
Figure 11: Central 128x128 pixels of the generated K-Band star cluster images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the calibrator images where the position angle is smaller than 180 degree and the altitude is 15.6 degree, 16.5 degree, 15.1 degree, and 55.7 degree. The second row contain the star cluster images at culmination (altitudes 27.3 degree, 37.3 degree, 57.3 degree, and 77.3 degree).

A comparison of the reconstructions depending on the atmospheric dispersion is presented in figure 12 (J-Band) and figure 13 (K-Band).

J-Band, coadded, -30 degree declinationJ-Band, coadded, -20 degree declinationJ-Band, coadded, 0 degree declinationJ-Band, coadded, +20 degree declination
J-Band, Richardson-Lucy reconstruction, -30 degree declinationJ-Band, Richardson-Lucy reconstruction, -20 degree declinationJ-Band, Richardson-Lucy reconstruction, 0 degree declinationJ-Band, Richardson-Lucy reconstruction, +20 degree declination
J-Band, Buidling-Block reconstruction, -30 degree declinationJ-Band, Buidling-Block reconstruction, -20 degree declinationJ-Band, Buidling-Block reconstruction, 0 degree declinationJ-Band, Buidling-Block reconstruction, +20 degree declination
Figure 12: Central 128x128 pixels of the coadded and reconstructed J-Band star cluster images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the coadded raw target images. The second row contain the reconstructions created with the Richardson-Lucy algorithm. For the third row, the Building-Block method was used.
K-Band, coadded, -30 degree declinationK-Band, coadded, -20 degree declinationK-Band, coadded, 0 degree declinationK-Band, coadded, +20 degree declination
K-Band, Richardson-Lucy reconstruction, -30 degree declinationK-Band, Richardson-Lucy reconstruction, -20 degree declinationK-Band, Richardson-Lucy reconstruction, 0 degree declinationK-Band, Richardson-Lucy reconstruction, +20 degree declination
K-Band, Buidling-Block reconstruction, -30 degree declinationK-Band, Buidling-Block reconstruction, -20 degree declinationK-Band, Buidling-Block reconstruction, 0 degree declinationK-Band, Buidling-Block reconstruction, +20 degree declination
Figure 13: Central 128x128 pixels of the coadded and reconstructed K-Band star cluster images at a declination of -30 degree (first column), -20 degree (second column), 0 degree (third column), and +20 degree (fourth column). The first row contains the coadded raw target images. The second row contain the reconstructions created with the Richardson-Lucy algorithm. For the third 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 atmospheric dispersion on the reconstruction error and the optimal iteration number is clearly visible.

Richardson-Lucy reconstruction error (J-Band)Richardson-Lucy reconstruction error (K-Band)
Buidling-Block reconstruction error (J-Band)Buidling-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). The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.

In the figures figure 15, figure 16, figure 17, and figure 18 the image and photometric errors depending on the atmospheric dispersion 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, -30 degree declination)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, -30 degree declination)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, -30 degree declination)Building-Block: Image and photometric errors depending on the iteration number (K-Band, -30 degree declination)
Figure 15: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The declination is -30 degree for the target and calibrator. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, -20 degree declination)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, -20 degree declination)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, -20 degree declination)Building-Block: Image and photometric errors depending on the iteration number (K-Band, -20 degree declination)
Figure 16: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The declination is -20 degree for the target and calibrator. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, 0 degree declination)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, 0 degree declination)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, 0 degree declination)Building-Block: Image and photometric errors depending on the iteration number (K-Band, 0 degree declination)
Figure 17: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The declination is 0 degree for the target and calibrator. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.
Richardson-Lucy: Image and photometric errors depending on the iteration number (J-Band, +20 degree declination)Richardson-Lucy: Image and photometric errors depending on the iteration number (K-Band, +20 degree declination)
Building-Block: Image and photometric errors depending on the iteration number (J-Band, +20 degree declination)Building-Block: Image and photometric errors depending on the iteration number (K-Band, +20 degree declination)
Figure 18: Image and photometric errors depending on the iteration number in the J-Band (left side) and K-Band (right side). The declination is +20 degree for the target and calibrator. The first row contain the reconstructions created with the Richardson-Lucy algorithm. For the second row, the Building-Block method was used.

In table 5 (J-Band) and table 6 (K-Band) the errors for the AGN depending on the atmospheric dispersion are presented.

TestcaseRichardson-LucyBuilding-Block
DecPAAltitudesImage ErrorMagnitude ErrorPosition ErrorImage ErrorMagnitude ErrorPosition Error
-30144 180 21615.6 27.3 15.60.100 (0.100)0.07, 0.04, 0.56, 0.10, 2.680.05, 0.04, 0.25, 0.57, 0.590.110 (0.110)0.07, 0.03, 0.85, 0.24, 2.430.05, 0.05, 0.49, 0.54, 0.60
-20135 180 22516.5 37.3 16.50.104 (0.104)0.04, 0.06, 0.53, 0.40, 0.980.06, 0.05, 0.18, 0.07, 0.570.108 (0.108)0.04, 0.05, 0.74, 0.54, 2.200.05, 0.06, 0.38, 0.01, 0.84
-10128 180 23215.0 47.3 15.00.070 (0.070)0.04, 0.04, 0.43, 0.40, 0.360.03, 0.02, 0.62, 0.05, 1.230.070 (0.070)0.03, 0.03, 0.70, 0.52, 0.960.03, 0.03, 0.36, 0.01, 1.17
0124 180 23615.1 57.3 15.10.049 (0.049)0.02, 0.03, 0.36, 0.26, 0.860.00, 0.02, 0.75, 0.60, 0.070.049 (0.049)0.02, 0.01, 0.59, 0.40, 2.210.00, 0.03, 0.32, 0.42, 0.50
10122 180 23828.9 67.3 28.90.072 (0.072)0.04, 0.03, 0.05, -0.02, 0.390.04, 0.04, 0.67, 0.36, 1.110.068 (0.068)0.04, 0.02, 0.15, 0.01, 0.680.04, 0.03, 0.60, 0.25, 0.92
20120 180 24055.7 77.3 55.70.030 (0.030)0.02, 0.01, 0.10, 0.24, 0.580.03, 0.03, 0.06, 0.14, 0.510.026 (0.026)0.02, 0.00, 0.19, 0.34, 1.150.01, 0.03, 0.07, 0.01, 0.47
Table 5: Star cluster: Errors depending on the atmospheric dispersion (short side of the J-Band).
TestcaseRichardson-LucyBuilding-Block
DecPAAltitudesImage ErrorMagnitude ErrorPosition ErrorImage ErrorMagnitude ErrorPosition Error
-30144 180 21615.6 27.3 15.60.167 (0.167)0.12, 0.03, 0.90, 0.55, 2.080.29, 0.15, 1.82, 1.19, 0.360.150 (0.150)0.13, 0.04, 0.81, 0.63, 1.940.27, 0.14, 1.62, 1.30, 0.36
-20135 180 22516.5 37.3 16.50.093 (0.093)0.03, 0.03, 0.88, 0.65, 2.070.08, 0.06, 0.75, 0.59, 0.390.065 (0.065)0.04, 0.03, 0.80, 0.76, 1.920.04, 0.05, 0.71, 0.48, 0.44
-10128 180 23215.0 47.3 15.00.070 (0.070)0.05, 0.03, 1.17, 0.23, 1.900.10, 0.05, 0.59, 1.39, 0.940.052 (0.052)0.05, 0.04, 1.06, 0.33, 1.510.09, 0.05, 0.24, 1.27, 1.72
0124 180 23615.1 57.3 15.10.075 (0.075)0.06, 0.03, 0.17, 0.09, 2.010.08, 0.03, 0.24, 1.29, 0.420.061 (0.061)0.07, 0.05, 0.19, 0.13, 1.880.08, 0.01, 0.19, 1.36, 0.45
10122 180 23828.9 67.3 28.90.067 (0.067)0.03, 0.02, 0.75, 0.70, 1.690.07, 0.05, 1.57, 2.26, 1.790.053 (0.053)0.03, 0.03, 0.68, 0.75, 1.750.06, 0.06, 1.89, 2.32, 1.07
20120 180 24055.7 77.3 55.70.057 (0.057)0.05, 0.03, 0.91, 0.98, 1.920.05, 0.02, 0.71, 0.60, 0.650.046 (0.046)0.05, 0.03, 0.87, 1.07, 1.900.03, 0.02, 0.61, 0.78, 0.45
Table 6: Star cluster: Errors depending on the atmospheric dispersion (K-Band).