Fast Fourier transform-based phase screen simulations give accurate results only when the screen size (G) is much larger than the outer scale parameter (L_{0}). Otherwise, they fall short in correctly predicting both the low and high frequency behaviors of turbulence-induced phase distortions. Subharmonic compensation is a commonly used technique that aids in low-frequency correction but does not solve the problem for all values of screen size to outer scale parameter ratios (G / L_{0}). A subharmonics-based approach will lead to unequal sampling or weights calculation for subharmonics addition at the low-frequency range and patch normalization factor. We have modified the subharmonics-based approach by introducing a Gaussian phase autocorrelation matrix that compensates for these shortfalls. We show that the maximum relative error in structure function with respect to theoretical value is as small as 0.5% to 3% for (G / L_{0}) ratio of 1/1000 even for screen sizes up to 100 m diameter. |
1.IntroductionAccurately simulating the atmospheric turbulence behavior is well recognized as very challenging. For a variety of purposes such as the design and development of adaptive optics systems, speckle imaging techniques, and atmospheric propagation studies, it is essential to simulate good atmospheric phase screen models. Methods based on Zernike polynomial expansions,^{1} fast Fourier transform (FFT)-based methods,^{2}^{–}^{8} and low-frequency optimization method^{9} have been in use for this purpose. The Zernike polynomial method, which is widely in use, has a limitation due to the maximum number of coefficients needed for accurate compensation. The optimization method, which compensates accurately for low-frequency part of the spectrum using unequal sampling and unequal weight in low-frequency region, does not cover high-frequency deficiencies. Among these, FFT-based methods are computer memory size friendly and widely accepted. But, FFT operators assume uniform sampling for the non-uniformly distributed phase power spectrum, which can lead to underestimation in the low- and high-frequency out of band regions, as shown in Fig. 1. Thus, it has limitations in recreating the true phase power spectrum. To compensate for low-frequency components, Johansson and Gavel^{3} suggested employing the modified subharmonics equation (an adaptation from Lane et al.^{10}), which works well up to an infinite outer scale length. Sedmak^{6} later compared the performance of this method with that of Lane et al.^{10} by actually calculating the phase structure function from the simulated screen. He improved upon Lane et al.^{10} by employing different fine-tuned subharmonic weights for different $G/{L}_{0}$ ratios. Results from his analysis show that these FFT-based simulations are accurate for large screen size ($G$) to outer scale parameter (${L}_{0}$) ratios. For a screen size of $G=200\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ and outer scale of ${L}_{0}=25\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, the maximum relative error in the simulation approaches 1%. Our simulations demonstrate that the errors from low-frequency components start shooting up once we move to smaller $G/{L}_{0}$ ratios, even after compensating with modified subharmonics. In Fig. 1, we show^{11} a typical situation where the simulation band $(\frac{1}{G}-\frac{1}{\mathrm{\Delta}})$ is actually smaller than full band $(\frac{1}{{L}_{0}}-\frac{1}{{l}_{0}})$, where $\mathrm{\Delta}$ is the sampling size defined as the ratio of screen size $G$ to sampling number $N$ and ${l}_{0}$ is the inner scale parameter. In practice, the simulations are often curtailed at the low-frequency end, to a few times the optical beam size (say as determined by the telescope or laser beam diameter), whereas at the high frequency end, they often extend to only a few times that determined by the Fried parameter ${r}_{0}$. Clearly, the larger the simulation band to full band ratio, the more accurate the simulated results will be. On the one hand, the apertures of upcoming and future astronomical telescopes are often of the same order or even larger than the typical median outer scale sizes of about 20 to 25 m.^{12} On the other hand, wavefront sensing and compensation technologies are fast progressing that Nyquist sampling at ${r}_{0}$ scales even for large aperture telescopes are becoming quite possible. Thus, atmospheric turbulence simulations have to deal with a wide range in a multi-dimensional parameter space. For working with very small apertures relative to the outer scale, it may appear that we need to simulate only a relatively small screen size. But cutting out small apertures from a larger screen introduces deviation from phase structure function due to misrepresentation of low-frequency components present in the small screen power spectrum. In this paper, we present an approach and a corresponding algorithm to deal with phase screen simulations for a wide range of $G/{L}_{0}$ ratios, using the FFT-based method.^{13} Our technique builds upon the modified subharmonic approach of Johansson and Gavel^{3} and is inspired by Jingsong Xiang’s work.^{7} It works well for space- and time-invariant, zero intermittency atmospheric turbulence. Section 2 explains how to obtain phase autocorrealtion matrix using phase power spectrum, Sec. 3 presents the algorithm part to compensate for the remaining error in phase structure function calculation, Sec. 4 steps through the implementation of the algorithm with the help of a flow chart, Sec. 5 covers the validation of the technique using results from simulated phase screens, and Sec. 6 provides the concluding remarks. 2.Obtaining Phase Autocorrelation Matrix Using Phase Power SpectrumThe 2D phase structure function and phase autocorrelation matrix are related as follows:^{14} where ${B}_{\varphi}(m,n)$ is the phase autocorrelation matrix and $(m,n)$ is the coordinates along $x$ and $y$ axis. The 2D phase autocorrelation matrices for the FFT-based phase screen and the modified subharmonic method by Johansson and Gavel^{3} are represented as follows.Eq. (2)$${B}_{\varphi}^{FFT}(m,n)=\sum _{{m}^{\prime}=-{N}_{x}/2}^{{N}_{x}/2-1}\sum _{{n}^{\prime}=-{N}_{y}/2}^{{N}_{y}/2-1}{f}_{FFT}^{2}({m}^{\prime},{n}^{\prime}){e}^{i2\pi (\frac{{m}^{\prime}m}{{N}_{x}}+\frac{{n}^{\prime}n}{{N}_{y}})},$$Eq. (3)$${B}_{\varphi}^{SUB}(m,n)=\sum _{p=1}^{{N}_{p}}\sum _{{m}^{\prime}=-3}^{2}\sum _{{n}^{\prime}=-3}^{2}{f}_{SUB}^{2}({m}^{\prime},{n}^{\prime}){e}^{i2\pi {3}^{-p}(\frac{({m}^{\prime}+0.5)m}{{N}_{x}}+\frac{({n}^{\prime}+0.5)n}{{N}_{y}})},$$The 2D phase autocorrelation matrix after compensating with subharmonics is represented as 3.Algorithm to Compensate for Residual Error in Phase Structure FunctionTo calculate the remaining error in the final ${B}_{\varphi}(m,n)$, Eq. (4) is converted to phase structure matrix ${D}_{\varphi}(m,n)$ with the help of Eq. (1) with the assumption that ${B}_{\varphi}^{FFT}(\mathrm{0,0})$ and ${B}_{\varphi}^{SUB}(\mathrm{0,0})$ are zero because we are not concerned about the piston component. This gives the following equation where ${D}_{\text{theory}}(m,n)$ is the well-known theoretical von-Kármán phase structure matrix,^{3} given as follows:Eq. (6)$${D}_{\text{theory}}(r)=6.16{r}_{0}^{5/3}[0.6{({L}_{0}/2\pi )}^{5/3}-\frac{{(r{L}_{0}/4\pi )}^{5/6}}{\gamma (11/6)}{K}_{5/6}(2\pi r/{L}_{0})],$$We need to compensate ${D}_{\varphi}$ so that ${D}_{\text{error}}$ is minimized. However, simply adding error correction terms in the ${D}_{\varphi}$ matrix directly would only introduce further error into the system while taking the Fourier transform. This is because any matrix or curve in general will have higher order moments. Thus, if we take the Fourier transform of the adjusted equation, the resultant curve will have completely different moments and hence power spectrum. This is because the transition between two steps in the error matrix will not be smooth, which introduces additional errors due to Gibb’s phenomena such as overshoots. Just curve fitting with any function does not satisfy the additional requirement of leaving the power spectrum unaffected by the process. What we really need is to introduce a smoothening operator such as a Gaussian function in the phase autocorrelation matrix, which exactly compensates for ${D}_{\text{error}}$. For that we have developed an iterative algorithm (see the flow chart shown in Fig. 4) and implemented it in Matlab. The algorithm looks for the perfect Gaussian curve that minimizes the ${D}_{\text{error}}$ matrix. We use Matlab cftool to initially determine the correct 1D Gaussian matrix and later convert it into a 2D matrix by exploiting the fact that ${B}_{\varphi}(r)$, ${B}_{\text{theory}}(r)$, and ${B}_{\mathrm{SUB}}(r)$ all are dependent on $r$ only and hence are center symmetric functions. We call the fitted Gaussian phase structure matrix ${D}_{\text{gauss}}$ and the corresponding Gaussian phase autocorrelaiton matrix ${B}_{\text{gauss}}$ [using Eq. (1)]. The final equation for ${B}_{\text{tot}}$ can then be written as Eq. (7)$${B}_{\text{tot}}(m,n)={B}_{\varphi}^{FFT}(m,n)+{B}_{\varphi}^{SUB}(m,n)-{D}_{\text{gauss}}(m,n)/2.$$Here, we have used ${B}_{\text{gauss}}=-{D}_{\text{gauss}}/2$ from Eq. (1). A look at the power spectrum of ${B}_{\mathrm{tot}}(m,n)$ in Fig. 2 shows that it contains negative terms^{7} for the case of $G/{L}_{0}<1$. Directly putting those frequency terms equal to zero leads to a loss in the energy spectrum. Hence, ${B}_{\mathrm{tot}}(m,n)$ matrix needs to be preprocessed to eliminate most of these negative values in the power spectrum. Over small frequencies, piston and tip/tilt components account for most of these high magnitude negative elements. Therefore, we first extract the piston and tip/tilt components from the phase autocorrelation matrix ${B}_{\text{tot}}$. The tip/tilt component from phase autocorrelation matrix is given as^{7} where ${\sigma}_{\text{tilt}}^{2}$ is the variance of the random tilt angle in the $x$ or $y$ directions and given as follows:^{7}Eq. (9)$${\sigma}_{\text{tilt}}^{2}=\frac{{B}_{\mathrm{tot}}(G/2+\mathrm{\Delta})-{B}_{\mathrm{tot}}(G/2)}{\mathrm{\Delta}(G-\mathrm{\Delta})/2}.$$After setting, ${B}_{\text{tilt}}(0)=0$, the remaining phase autocorrelation matrix is given as follows: The power spectra ${f}_{\text{high}}^{2}$ and ${f}_{\text{tilt}}^{2}$ of the phase autocorrelation matrices ${B}_{\text{high}}(r)$ and ${B}_{\text{tilt}}(r)$ are obtained by standard Fourier transformation. Figure 3, shows the remaining negative power elements present in the power spectrum of ${B}_{\text{high}}$ matrix. In comparison to Fig. 2, the largest negative power contributions fall by factor of three order of magnitude. Now, we set the negative values in ${f}_{\text{high}}^{2}$ equal to zero by hand. The new error matrix is given as where ${B}_{\text{high}}^{\prime}(r)$ is the phase autocorrelation matrix obtained after setting the negative elements in ${f}_{\text{high}}^{2}$ to zero. The residual error that is present in the high-frequency region can then be reduced with the help of a Gaussian smoothing operator, using Matlab fmincon tool. The high-frequency compensated matrix is given asEq. (12)$${B}_{\text{high}}^{\mathrm{comp}}(r)={B}_{\text{high}}(r)-{H}_{\text{high}}^{\mathrm{comp}}(r){B}_{\text{high}}^{\mathrm{err}}(r),$$4.Implementation of the Compensation AlgorithmIn this section, we explain the error compensation algorithm with the help of the flow chart shown in Fig. 4. Brief explanations of each of the steps from ${L}_{1}$ to ${L}_{12}$ are given below.
fmincon: High Frequency error Optimization
Table 1 shows the result of curve fitting using cftool for different cases of $G/{L}_{0}$ and N, which demonstrates that the Gaussian error matrix can compensate for a wide range of $G/{L}_{0}$ ratios and under different sampling constraints. Table 1Result of curve fitting against Gaussian function for different cases of G/L0 and N in terms of MRE for fixed r0=0.2 m.
4.1.ExampleTo illustrate the robustness of the above algorithm, we have taken an example with $G=80\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, say for a large future telescope, $N=256$, and median value of ${L}_{0}=20\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. The output from the above algorithm corresponding to minimum error entry as in (step ${L}_{12}$), has been plotted against $EF=5$ and ${N}_{p}=8$. Figure 5 shows a 3D rendering of ${D}_{\text{error}}$ matrix with a maximum separation of up to 40 m, corresponding to Eq. (4). Figure 6 shows 1D ${D}_{\text{error}}$ matrix (radial section from 3D ${D}_{\text{error}}$ matrix) along with 1D fitted curve ${D}_{\text{gauss}}$ including the extrapolated part, for a maximum separation of up to 40 m. Lastly, MRE values are stored against 500 entries corresponding to ${N}_{p}$ ranging from 1 to 10, $EF$ ranging from 1 to 10 and $SoG$ ranging from 2 to 6 after performing cftool fitting. This has been arranged in descending order and shown in Fig. 7, which illustrates a large set of iterations where errors are $<1\%$ and entry with minimum MRE has been picked up. Typical time required to perform each iteration for this case is $\approx 4.9\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{s}$ on 2.3 GHz quad-core Intel Core i5 Macbook Pro 2018 model. 5.Validation via Phase Structure Function Calculated from Simulated Phase ScreenTo obtain the phase screen $\varphi (m,n)$ from the power spectrum, the following relation is used:^{7} Eq. (13)$$\varphi (m\mathrm{\Delta},n\mathrm{\Delta})=\sum _{{m}^{\prime}=-N/2}^{N/2-1}\sum _{{n}^{\prime}=-N/2}^{N/2-1}[{R}_{a}({m}^{\prime},{n}^{\prime})+\mathrm{i}{R}_{b}({m}^{\prime},{n}^{\prime})]f({m}^{\prime}{\mathrm{\Delta}}^{\prime},{n}^{\prime}{\mathrm{\Delta}}^{\prime})\mathrm{exp}[\mathrm{i}2\pi ({m}^{\prime}m+{n}^{\prime}n)/N],$$For validation, we consider scenarios of apertures up to 40 m, i.e., $G=80\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, at a median ${L}_{0}=20\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ for two different sampling levels $N=256$ and 512. The phase structure function, defined as an ensemble average of differences of phases at various separation,^{14} has been averaged over 100 K independent frames. The relative error in phase structure function is calculated as follows: Eq. (14)$$err(r)=\frac{{D}_{\varphi}^{\text{sim}}(r)-{D}_{\text{theory}}(r)}{{D}_{\text{theory}}(r)}.$$Here, ${D}_{\varphi}^{\text{sim}}(r)$ is the phase structure function from the simulated phase screen. The magnitude of the peak relative error $\mathrm{max}(|err(r)|)$ is $<1.6\%$ for $N=256$ and $<0.5\%$ for $N=512$ as shown in Fig. 8. We also illustrate the performance (shown in Fig. 9) with parameters $G=1\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, ${L}_{0}=100\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, and 1000 m, $N=128$, ${r}_{0}=0.2\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, which cover the extreme cases (very low $G/{L}_{0}$ ratios) which leads to the maximum error in the simulation. The magnitude of the peak relative error $\mathrm{max}(|\mathrm{err}(r)|)$ is $<1.6\%$ for ${L}_{0}=100\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ and $<1.8\%$ for ${L}_{0}=1000\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. Figure 10 shows one realization of the corresponding phase screen plots for ${L}_{0}=100\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. Figure 11 contains results of magnitude of the peak relative error in ${D}_{\varphi}^{\text{sim}}(r)$ for the case of different sampling points $N=128/256/512/1024$, for ${L}_{0}$ ranges up to 1024 m, ${r}_{0}=0.2\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$ and $G=2\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. Similarly, Fig. 12 contains results of magnitude of the peak relative error in ${D}_{\varphi}^{\text{sim}}(r)$ for the case of different sampling points $N=128/256/512/1024$, for $G$ ranges up to 100 m, ${r}_{0}=0.2\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$, and ${L}_{0}=25\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{m}$. There are some outliers that have a high residual error as shown in Figs. 11 and 12, because we have not set the phase autocorrelation matrix to zero for $r>G/2$. The reason for this stems from the non-zero value of ${B}_{\text{high}}(r)$ at $r>G/2$, where ${B}_{\text{high}}(r)$ is formed from the removal of piston and tilt from ${B}_{\mathrm{tot}}(r)$. This can be resolved using a better smoothening operator, which we can multiply with ${B}_{\text{high}}$ so that it falls to zero progressively and not sharply. This can provide further improvement in the compensation. 6.ConclusionIn this paper, we put forward a new method to compensate for the residual error in both the low and/or high-frequency region of FFT simulated phase screens that remain even after compensating with the modified subharmonic method. This method provides accurate phase screen structure for even $G/{L}_{0}$ ratios as small as 1/1000 plus screen sizes as large as 100 m. No patch normalization factor is needed, no need to calculate subharmonic weight coefficient^{10} and weights to compensate for high-frequency components, as done by Sedmak.^{6} While adequately large $G/{L}_{0}$ ratios may be the natural choice for modern large telescopes, simulations that deal with applications such as laser beam propagation through turbulent atmospheres would tend to have very small $G/{L}_{0}$ ratios. The method we propose is independent of the $G/{L}_{0}$ ratio choice. However, we emphasize that properly sampling ${r}_{0}$ and the high-frequency phase spectrum forces $N$ to be at least larger than $(2G/{r}_{0})$ and preferably up to the inner scale limit $(2G/{l}_{0})$. Currently, we have demonstrated this technique for only circular screens. We have used a GPU processor with total number of 128 cores, such that each iteration runs independently on each core. We have fixed the number of iterations to 500, although increasing this will lead to improvement of errors in some cases. Each core takes about $\sim 0.06$, $\sim 0.1$, $\sim 0.36$, and $\sim 1.1\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{min}$ for sampling sizes of 128, 256, 512, and 1024. On the above GPU system, this translates to total computing times for error minimization of about $\sim 0.2$, $\sim 0.5$, $\sim 1.25$, and $\sim 4\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{min}$ for sampling sizes of 128, 256, 512, and 1024, respectively. Once the coefficients are determined, generating multiple phase screen realizations from the corresponding power spectrum takes a few milliseconds at most on 2.3 GHz quad-core Intel Core i5 Macbook Pro 2018 model. Then it takes $<1\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{min}$ to $\sim 10\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{min}$ for averaging over 100k phase screens, for sampling sizes ranging from 128 to 1024. The uniqueness of our approach is its ability to deal with any $G/{L}_{0}$ ratio within a very broad range, in an automated iterative process with little human intervention needed for tuning of parameters. Any standard FFT-based approach (say Sedmak’s^{6} compensated approach) for a given computer platform is computationally fast, only if we already have determined proper measures of the various compensating components such as the patch normalization factor, subharmonic weights, and high-frequency weights. Typically, determining these compensations is where the difficulty is due to mathematical complexity, algorithmic limitations, and/or computational power requirements. Our algorithm accomplishes the determination of the required compensation in very little time, with fairly reasonable computational power while at the same time keeping the residual errors competitively low using an appropriate compensator. Other existing FFT based approaches have limitations in their operable $G/{L}_{0}$ range. For example, Xiang et al.^{7} offer a very computationally fast approach but does not apply subharmonic compensation. Zhang et al.^{9} does not consider compensation for the high-frequency error, thus leaving a residual error of more than 100% in the high-frequency region. Sedmak’s^{6} approach needs the determination of accurate subharmonic weights for different $G/{L}_{0}$ ratios. The accuracy of our method from low-frequency to the high-frequency range is between 0.5% and 3% for $G/{L}_{0}$ as low as 1/1000 and screen size up to 100 m in diameter. AcknowledgmentsWe would like to thank Sedmak for providing insights into the nature of atmospheric phase power spectrum through private communication. We also thank Xiang for sharing his MATLAB code that calculates the phase structure function quickly for a large number of phase screens. We acknowledge usage of IUCAA’s Pegasus cluster computer for running multiple independent iterations in parallel. ReferencesN. A. Roddier,
“Atmospheric wavefront simulation and Zernike polynomials,”
Proc. SPIE, 1237 668
–679
(1990). https://doi.org/10.1117/12.19346 PSISDG 0277-786X Google Scholar
B. J. Herman and L. A. Strugala,
“Method for inclusion of low-frequency contributions in numerical representation of atmospheric turbulence,”
Proc. SPIE, 1221 183
–192
(1990). https://doi.org/10.1117/12.18342 PSISDG 0277-786X Google Scholar
E. M. Johansson and D. T. Gavel,
“Simulation of stellar speckle imaging,”
Proc. SPIE, 2200 372
–383
(1994). https://doi.org/10.1117/12.177254 PSISDG 0277-786X Google Scholar
G. Sedmak,
“Performance analysis of and compensation for aspect-ratio effects of fast-Fourier-transform-based simulations of large atmospheric wave fronts,”
Appl. Opt., 37
(21), 4605
–4613
(1998). https://doi.org/10.1364/AO.37.004605 APOPAI 0003-6935 Google Scholar
B. L. McGlamery,
“Computer simulation studies of compensation of turbulence degraded images,”
Proc. SPIE, 0074 225
–233
(1976). https://doi.org/10.1117/12.954724 PSISDG 0277-786X Google Scholar
G. Sedmak,
“Implementation of fast-Fourier-transform-based simulations of extra-large atmospheric phase and scintillation screens,”
Appl. Opt., 43
(23), 4527
–4538
(2004). https://doi.org/10.1364/AO.43.004527 APOPAI 0003-6935 Google Scholar
J. Xiang,
“Fast and accurate simulation of the turbulent phase screen using fast Fourier transform,”
Opt. Eng., 53
(1), 016110
(2014). https://doi.org/10.1117/1.OE.53.1.016110 Google Scholar
J. Xiang,
“Accurate compensation of the low-frequency components for the FFT-based turbulent phase screen,”
Opt. Express, 20
(1), 681
–687
(2012). https://doi.org/10.1364/OE.20.000681 OPEXFF 1094-4087 Google Scholar
D. Zhang et al.,
“Accurate simulation of turbulent phase screen using optimization method,”
Optik, 178 1023
–1028
(2019). https://doi.org/10.1016/j.ijleo.2018.10.083 OTIKAJ 0030-4026 Google Scholar
R. Lane et al.,
“Simulation of a Kolmogorov phase screen,”
Waves Random Media, 2
(3), 209
–224
(1992). https://doi.org/10.1088/0959-7174/2/3/003 WRMEEV 0959-7174 Google Scholar
G. Sedmak,
(2019). Google Scholar
A. Ziad,
“Review of the outer scale of the atmospheric turbulence,”
Proc. SPIE, 9909 99091K
(2016). https://doi.org/10.1117/12.2231375 PSISDG 0277-786X Google Scholar
S. Chhabra et al.,
“Gaussian phase autocorrelation as an accurate compensator for FFT-based atmospheric phase screen simulations,”
Proc. SPIE, 11448 114487U
(2020). https://doi.org/10.1117/12.2575934 PSISDG 0277-786X Google Scholar
F. Roddier, Adaptive Optics in Astronomy, Cambridge University Press(1999). Google Scholar
BiographySorabh Chhabra is a PhD student at Inter University Center for Astronomy and Astrophysics (IUCAA), Pune. He received his B. Tech degree in electronics and communication from Delhi Technological University (formally Delhi College of Engineering) in 2016 and joined for his PhD in the same year in the Instrumentation Department at IUCAA. |