scispace - formally typeset
Search or ask a question

Showing papers in "Journal of The Optical Society of America A-optics Image Science and Vision in 2012"


Journal ArticleDOI
TL;DR: It is shown that an alternative form of multi-measurement imaging, ptychography, can be extended to three dimensions and can successfully recover images in the presence of multiple scattering and when the projection approximation is not applicable.
Abstract: Generally, methods of three-dimensional imaging such as confocal microscopy and computed tomography rely on two essentials: multiple measurements (at a range of focus positions or rotations) and a weakly scattering specimen (to avoid distortion of the focal spot in the confocal microscope or to satisfy the projection approximation in tomography). Here we show that an alternative form of multi-measurement imaging, ptychography, can be extended to three dimensions and can successfully recover images in the presence of multiple scattering and when the projection approximation is not applicable. We demonstrate our technique experimentally using visible light, where it has applications in imaging thick samples such as biological tissues; however the results also have important implications for x ray and electron imaging.

319 citations


Journal ArticleDOI
TL;DR: It is proposed that founder events and genetic drift, rather than natural selection, are the cause of differences in the prevalence of inherited red-green color deficiency in Europeans and Asians.
Abstract: Literature that describes the prevalence of inherited red-green color deficiency in different populations is reviewed. Large random population surveys show that the prevalence of deficiency in European Caucasians is about 8% in men and about 0.4% in women and between 4% and 6.5% in men of Chinese and Japanese ethnicity. However, the male: female prevalence ratio is markedly different in Europeans and Asians. Recent surveys suggest that the prevalence is rising in men of African ethnicity and in geographic areas that have been settled by incoming migrants. It is proposed that founder events and genetic drift, rather than natural selection, are the cause of these differences.

185 citations


Journal ArticleDOI
TL;DR: The behavior of the spectral density and the state of coherence of beamlike fields generated by planar stochastic sources on propagation in free space and linear isotropic random media is examined.
Abstract: In a recent publication [Opt. Lett.37, 2970 (2012)10.1364/OL.37.002970], a novel class of planar stochastic sources, generating far fields with flat intensity profiles, was introduced. In this paper we examine the behavior of the spectral density and the state of coherence of beamlike fields generated by such sources on propagation in free space and linear isotropic random media. In particular, we find that at sufficiently large distances from the source, the medium destroys the flat intensity profile, even if it remains such for intermediate distances from the source.

177 citations


Journal ArticleDOI
TL;DR: This work reports on an architecture that acquires the two-dimensional spatial Fourier transform of the target object and determines its image signature, resolution, and signal-to-noise ratio in the presence of practical constraints such as atmospheric turbulence, background radiation, and photodetector noise.
Abstract: Computational ghost imaging is a structured-illumination active imager coupled with a single-pixel detector that has potential applications in remote sensing. Here we report on an architecture that acquires the two-dimensional spatial Fourier transform of the target object (which can be inverted to obtain a conventional image). We determine its image signature, resolution, and signal-to-noise ratio in the presence of practical constraints such as atmospheric turbulence, background radiation, and photodetector noise. We consider a bistatic imaging geometry and quantify the resolution impact of nonuniform Kolmogorov-spectrum turbulence along the propagation paths. We show that, in some cases, short-exposure intensity averaging can mitigate atmospheric-turbulence-induced resolution loss. Our analysis reveals some key performance differences between computational ghost imaging and conventional active imaging, and identifies scenarios in which theory predicts that the former will perform better than the latter.

167 citations


Journal ArticleDOI
TL;DR: The computational ghost imaging with a phase spatial light modulator (SLM) for wave field coding is considered and an approximative Gaussian distribution with an invariant variance results in the algorithm that is efficient for Poissonian observations.
Abstract: The computational ghost imaging with a phase spatial light modulator (SLM) for wave field coding is considered. A transmission-mask amplitude object is reconstructed from multiple intensity observations. Compressive techniques are used in order to gain a successful image reconstruction with a number of observations (measurement experiments), which is smaller than the image size. Maximum likelihood style algorithms are developed, respectively, for Poissonian and approximate Gaussian modeling of random observations. A sparse and overcomplete modeling of the object enables the advanced high accuracy and sharp imaging. Numerical experiments demonstrate that an approximative Gaussian distribution with an invariant variance results in the algorithm that is efficient for Poissonian observations.

132 citations


Journal ArticleDOI
TL;DR: Numerical analysis of the theoretical sensitivity limit of the localized surface plasmon resonance (LSPR) to the surrounding dielectric environment shows that, for thin (≤10 nm) analyte layers, sensitivity of the LSPR is comparable with a traditional surface plAsmon resonance sensor and L SPR has the potential to be significantly less sensitive to temperature fluctuations.
Abstract: In this paper, the theoretical sensitivity limit of the localized surface plasmon resonance (LSPR) to the surrounding dielectric environment is discussed. The presented theoretical analysis of the LSPR phenomenon is based on perturbation theory. Derived results can be further simplified assuming quasistatic limit. The developed theory shows that LSPR has a detection capability limit independent of the particle shape or arrangement. For a given structure, sensitivity is directly proportional to the resonance wavelength and depends on the fraction of the electromagnetic energy confined within the sensing volume. This fraction is always less than unity; therefore, one should not expect to find an optimized nanofeature geometry with a dramatic increase in sensitivity at a given wavelength. All theoretical results are supported by finite-difference time-domain calculations for gold nanoparticles of different geometries (rings, split rings, paired rings, and ring sandwiches). Numerical sensitivity calculations based on the shift of the extinction peak are in good agreement with values estimated by perturbation theory. Numerical analysis shows that, for thin (≤10 nm) analyte layers, sensitivity of the LSPR is comparable with a traditional surface plasmon resonance sensor and LSPR has the potential to be significantly less sensitive to temperature fluctuations.

91 citations


Journal ArticleDOI
TL;DR: The reduced-support constraint approach is shown to be effective in escaping stagnation caused by the twin-image problem and arises when the retrieved Fourier-domain phase is divided into sets of regions, some of which reconstruct the object while others the twin.
Abstract: The twin-image problem in phase retrieval is characterized by the simultaneous occurrence of features from the original object and its inversion about the origin (twin image). This problem can occur in reconstructions for which the object support is centrosymmetric or loose, and in severe cases it can greatly hinder image quality. In this paper we examine this problem and find that it arises when the retrieved Fourier-domain phase is divided into sets of regions, some of which reconstruct the object while others the twin. We examine sample reconstructions that present the twin-image problem to different extents and find that, even when the twin-image problem is not visually evident, it can exist in small regions of the retrieved Fourier phase. The reduced-support constraint approach is shown to be effective in escaping stagnation caused by the twin-image problem.

84 citations


Journal ArticleDOI
TL;DR: Attenuation by the presence of dust degrades the IR channel but exhibits almost no measurable impact on the THz signal, and numerical simulations of THz attenuation with different dust concentrations agree with the measured results.
Abstract: In order to study and compare propagation features of terahertz (THz) links with infrared (IR) links under different weather conditions such as turbulence, fog, and dust particles, THz and IR free space communication links at 625 GHz carrier frequency and 1.5 μm wavelength, respectively, with a maximum data rate of 2.5 Gb/s have been developed. After propagating through the same channel perturbation caused by dust, attenuation of the carrier frequencies by dust as well as scintillation effects on both channels are analyzed by measuring the power attenuation and bit error rates. Attenuation by the presence of dust degrades the IR channel but exhibits almost no measurable impact on the THz signal. Numerical simulations of THz attenuation with different dust concentrations are presented and agree with the measured results.

82 citations


Journal ArticleDOI
TL;DR: It is shown that the new sources are capable of producing beams with polarization properties that evolve on propagation in a manner much more complex compared to the well-known electromagnetic Gaussian Schell-model beams.
Abstract: A class of electromagnetic sources with nonuniformly distributed field correlations is introduced. The conditions on source parameters guaranteeing that the source generates a physical beam are derived. It is shown that the new sources are capable of producing beams with polarization properties that evolve on propagation in a manner much more complex compared to the well-known electromagnetic Gaussian Schell-model beams.

81 citations


Journal ArticleDOI
TL;DR: A lab setup for analyzing impairments of terahertz (THz) and infrared (IR) free space links caused by local refraction index changes in the signal's propagation paths that could be induced by turbulence, particles, humidity, etc.
Abstract: We describe a lab setup for analyzing impairments of terahertz (THz) and infrared (IR) free space links caused by local refraction index changes in the signal’s propagation paths that could be induced by turbulence, particles, humidity, etc. A THz signal comprising a 2.5 Gb/s data load modulated on a carrier at 625 GHz, is launched through a weather emulating chamber, detected, and its performance analyzed. An IR beam at 1.5 um wavelength carrying the same data load is superposed with the THz beam, propagating through the same weather conditions and also performance analyzed. We modulate the IR channel with a usual non-return-to-zero (NRZ) format but use duobinary coding for driving our THz source, which enables signaling at high data rate and higher output power. As both beams pass through the same channel perturbations and as their degradations are recorded simultaneously we can simultaneously compare the weather impact on both. We investigate scintillation and fog attenuation effects for the THz and IR signals by measuring bit error rates (BER), signal power, and phase front distortions.

78 citations


Journal ArticleDOI
TL;DR: This paper shows that this adaptation of the Cumulative Reconstructor gives the same reconstruction quality as the original algorithm and leads to a significant improvement with respect to noise propagation.
Abstract: The Cumulative Reconstructor is an accurate, extremely fast reconstruction algorithm for Shack–Hartmann wavefront sensor data. But it has shown an unacceptable high noise propagation for large apertures. Therefore, in this paper we describe a domain decomposition approach to deal with this drawback. We show that this adaptation of the algorithm gives the same reconstruction quality as the original algorithm and leads to a significant improvement with respect to noise propagation. The method is combined with an integral control and compared to the classical matrix vector multiplication algorithm on an end-to-end simulation of a single conjugate adaptive optics system. The reconstruction time is 20n (number of subapertures), and the method is parallelizable.

Journal ArticleDOI
TL;DR: This paper presents a detailed numerical study on the performance of the standard phasing algorithms with random phase illumination (RPI), and it is shown that RPI with σ=2 is sufficient for phasing complex-valued images under a sector condition and ρ=1 forphasing nonnegative images.
Abstract: This paper presents a detailed numerical study on the performance of the standard phasing algorithms with random phase illumination (RPI). Phasing with high resolution RPI and the oversampling ratio σ=4 determines a unique phasing solution up to a global phase factor. Under this condition, the standard phasing algorithms converge rapidly to the true solution without stagnation. Excellent approximation is achieved after a small number of iterations, not just with high resolution but also low resolution RPI in the presence of additive as well multiplicative noises. It is shown that RPI with σ=2 is sufficient for phasing complex-valued images under a sector condition and σ=1 for phasing nonnegative images. The error-reduction algorithm with RPI is proved to converge to the true solution under proper conditions.

Journal ArticleDOI
TL;DR: Initial benign deterioration of color discrimination in the 40+ decade is revealed as an incremental loss of discrimination along the Deutan axis (Trivector test), and in the 50+ decade, as an elongation of the major axes of all three ellipses (Ellipses test).
Abstract: Color discrimination was estimated using the Cambridge Colour Test (CCT) in 160 normal trichromats of four life decades, 20-59 years of age. For each age cohort, medians and tolerance limits of the CCT parameters are tabulated. Compared across the age cohorts (Kruskal-Wallis test), the Trivector test showed increases in the three vectors, Protan, Deutan, and Tritan, with advancing age; the Ellipses test revealed significant elongation of the major axes of all three ellipses but no changes in either the axis ratio or the angle of the ellipse major axis. Multiple comparisons (Mann-Whitney test) between the cohorts of four age decades (20+,…,50+) revealed initial benign deterioration of color discrimination in the 40+ decade, as an incremental loss of discrimination along the Deutan axis (Trivector test), and in the 50+ decade, as an elongation of the major axes of all three ellipses (Ellipses test).

Journal ArticleDOI
TL;DR: Algorithms are presented to construct synthetic images in which local image statistics--including luminance distributions, pair-wise correlations, and higher-order correlations--are explicitly specified and all other statistics are determined implicitly by maximum-entropy, to measure the sensitivity of the human visual system to local imageStatistics and to sample their interactions.
Abstract: The space of visual signals is high-dimensional and natural visual images have a highly complex statistical structure. While many studies suggest that only a limited number of image statistics are used for perceptual judgments, a full understanding of visual function requires analysis not only of the impact of individual image statistics, but also, how they interact. In natural images, these statistical elements (luminance distributions, correlations of low and high order, edges, occlusions, etc.) are intermixed, and their effects are difficult to disentangle. Thus, there is a need for construction of stimuli in which one or more statistical elements are introduced in a controlled fashion, so that their individual and joint contributions can be analyzed. With this as motivation, we present algorithms to construct synthetic images in which local image statistics—including luminance distributions, pair-wise correlations, and higher-order correlations—are explicitly specified and all other statistics are determined implicitly by maximum-entropy. We then apply this approach to measure the sensitivity of the human visual system to local image statistics and to sample their interactions.

Journal ArticleDOI
TL;DR: In this paper, a pixelated gradient-based optical proximity correction (OPC) and phase-shifting mask (PSM) optimization method is proposed to enhance the resolution in immersion lithography systems.
Abstract: Recently, a set of gradient-based optical proximity correction (OPC) and phase-shifting mask (PSM) optimization methods has been developed to solve for the inverse lithography problem under scalar imaging models, which are only accurate for numerical apertures (NAs) of less than approximately 0.4. However, as lithography technology enters the 45 nm realm, immersion lithography systems with hyper-NA (NA>1) are now extensively used in the semiconductor industry. For the hyper-NA lithography systems, the vector nature of the electromagnetic field must be taken into account, leading to the vector imaging models. Thus, the OPC and PSM optimization approaches developed under the scalar imaging models are inadequate to enhance the resolution in immersion lithography systems. This paper focuses on developing pixelated gradient-based OPC and PSM optimization algorithms under a vector imaging model. We first formulate the mask optimization framework, in which the imaging process of the optical lithography system is represented by an integrative and analytic vector imaging model. A gradient-based algorithm is then used to optimize the mask iteratively. Subsequently, a generalized wavelet penalty is proposed to keep a balance between the mask complexity and convergence errors. Finally, a set of methods is exploited to speed up the proposed algorithms.

Journal ArticleDOI
TL;DR: This review paper addresses typical mistakes and omissions that involve theoretical research and modeling of optical propagation through atmospheric turbulence and some misconceptions regarding short-exposure imaging, propagation of polarized waves, and calculations of the scintillation index of the beam waves.
Abstract: This review paper addresses typical mistakes and omissions that involve theoretical research and modeling of optical propagation through atmospheric turbulence. We discuss the disregard of some general properties of narrow-angle propagation in refractive random media, the careless use of simplified models of turbulence, and omissions in the calculations of the second moment of the propagating wave. We also review some misconceptions regarding short-exposure imaging, propagation of polarized waves, and calculations of the scintillation index of the beam waves.

Journal ArticleDOI
TL;DR: It is shown that arrangements of multiple convex and concave axicons may be implemented to optimize the depth of focus in a miniaturized OCT system, using a telescopic optical arrangement of considerably shorter optical system length than that achievable with classical micro-optics.
Abstract: Employing Bessel beams in imaging takes advantage of their self-reconstructing properties to achieve small focal points while maintaining a large depth of focus. Bessel beams are efficiently generated using axicons, and their utility in scanning imaging systems, such as optical coherence tomography (OCT), has been demonstrated. As these systems are miniaturized to allow, for example, endoscopic implementations, micro-axicons are required to assure the maintenance of a large depth of focus. We demonstrate here the design, fabrication, and application of molded micro-axicons for use in silicon-based micro-optical benches. It is shown that arrangements of multiple convex and concave axicons may be implemented to optimize the depth of focus in a miniaturized OCT system, using a telescopic optical arrangement of considerably shorter optical system length than that achievable with classical micro-optics.

Journal ArticleDOI
TL;DR: A multimodal adaptive optics system for high-resolution multifunctional use in a variety of research and clinical applications and for efficient location and orientation of retinal regions of interest is developed.
Abstract: Optical coherence tomography (OCT) and scanning laser ophthalmoscopy (SLO) are complementary imaging modalities, the combination of which can provide clinicians with a wealth of information to detect retinal diseases, monitor disease progression, or assess new therapies. Adaptive optics (AO) is a tool that enables correction of wavefront distortions from ocular aberrations. We have developed a multimodal adaptive optics system (MAOS) for high-resolution multifunctional use in a variety of research and clinical applications. The system integrates both OCT and SLO imaging channels into an AO beam path. The optics and hardware were designed with specific features for simultaneous SLO/OCT output, for high-fidelity AO correction, for use in humans, primates, and small animals, and for efficient location and orientation of retinal regions of interest. The MAOS system was tested on human subjects and rodents. The design, performance characterization, and initial representative results from the human and animal studies are presented and discussed.

Journal ArticleDOI
TL;DR: The semianalytic approach allows efficient and accurate calculation of the absorption of an array with a complex unit cell and gives direct physical insight into the absorption mechanism in such structures, which can be used to enhance the absorption.
Abstract: A finite element-based modal formulation of diffraction of a plane wave by an absorbing photonic crystal slab of arbitrary geometry is developed for photovoltaic applications. The semianalytic approach allows efficient and accurate calculation of the absorption of an array with a complex unit cell. This approach gives direct physical insight into the absorption mechanism in such structures, which can be used to enhance the absorption. The verification and validation of this approach is applied to a silicon nanowire array, and the efficiency and accuracy of the method is demonstrated. The method is ideally suited to studying the manner in which spectral properties (e.g., absorption) vary with the thickness of the array, and we demonstrate this with efficient calculations that can identify an optimal geometry.

Journal ArticleDOI
TL;DR: The Stokes vector lidar equation is introduced, which is a full description of polarization in lidar from laser output to detector, and this theoretical description is interpreted in the context of forward polar decomposition of Mueller matrices where distinct polarization attributes of diattenuation, retardance, and depolarization are elucidated.
Abstract: Polarization measurements have become nearly indispensible in lidar cloud and aerosol studies. Despite polarization’s widespread use in lidar, its theoretical description has been widely varying in accuracy and completeness. Incomplete polarization lidar descriptions invariably result in poor accountability for scatterer properties and instrument effects, reducing data accuracy and disallowing the intercomparison of polarization lidar data between different systems. We introduce here the Stokes vector lidar equation, which is a full description of polarization in lidar from laser output to detector. We then interpret this theoretical description in the context of forward polar decomposition of Mueller matrices where distinct polarization attributes of diattenuation, retardance, and depolarization are elucidated. This decomposition can be applied to scattering matrices, where volumes consisting of randomly oriented particles are strictly depolarizing, while oriented ice crystals can be diattenuating, retarding, and depolarizing. For instrument effects we provide a description of how different polarization attributes will impact lidar measurements. This includes coupling effects due to retarding and depolarization attributes of the receiver, which have no description in scalar representations of polarization lidar. We also describe how the effects of polarizance in the receiver can result in nonorthogonal polarization detection channels. This violates one of the most common assumptions in polarization lidar operation.

Journal ArticleDOI
TL;DR: Within the Rayleigh approximation, an individual ellipsoidal metal nanorod that is optically confined in three dimensions using a single focused laser beam is investigated, revealing a significantly different behavior if the trapping wavelength is longer or shorter than the wavelength corresponding to the longitudinal plasmon resonance mode.
Abstract: Within the Rayleigh approximation, we investigate the behavior of an individual ellipsoidal metal nanorod that is optically confined in three dimensions using a single focused laser beam. We focus on the description of the optical torque and optical force acting upon the nanorod placed into a linearly polarized Gaussian beam (scalar description of the electric field) or a strongly focused beam (vector field description). The study comprises the influence of the trapping laser wavelength, the angular aperture of focusing optics, the orientation of the ellipsoidal nanorod, and the aspect ratio of its principal axes. The results reveal a significantly different behavior of the nanorod if the trapping wavelength is longer or shorter than the wavelength corresponding to the longitudinal plasmon resonance mode. Published experimental observations are compared with our theoretical predictions with satisfactory results.

Journal ArticleDOI
TL;DR: Different beam combinations for stimulated emission depletion microscopy are compared, finding that azimuthal polarization combination consistently produces spots 15%-30% smaller than the commonly used, circularly polarized light combination.
Abstract: We compare different beam combinations for stimulated emission depletion microscopy. The four considered copolarized, mutually symmetric, but complementary write + erase beam combinations are circularly polarized beam + circularly polarized vortex with charge +1 or −1, azimuthally polarized with a vortex + azimuthally polarized, and radially polarized beam + radially polarized with a vortex. The resulting fluorescent spot was calculated for plane incident pump and erase beams, for plane waves with added high NA annular ring apertures, and when both incident beams were optimized with amplitude–phase masks. For all three incident wave cases, the azimuthal polarization combination consistently produces spots 15%–30% smaller than the commonly used, circularly polarized light combination (the first from above). The two other polarization combinations produce even smaller, of the order of nanometers/0.003λ, fluorescent spots with a caveat of having nonzero erase beam intensity in the center. Nevertheless, these combinations can be advantageous when exploiting PF, i.e., using molecules that respond solely to the longitudinal (or only to transversal) component of the illuminating field.

Journal ArticleDOI
TL;DR: The three-dimensional radiative transfer equation is solved for modeling the light propagation in anisotropically scattering semi-infinite media such as biological tissue, considering the effect of internal reflection at the interfaces using the two-dimensional Fourier transform and the modified spherical harmonics method.
Abstract: The three-dimensional radiative transfer equation is solved for modeling the light propagation in anisotropically scattering semi-infinite media such as biological tissue, considering the effect of internal reflection at the interfaces. The two-dimensional Fourier transform and the modified spherical harmonics method are applied to derive the general solution to the associated homogeneous problem in terms of analytical functions. The obtained solution is used for solving boundary-value problems, which are important for applications in the biomedical optics field. The derived equations are successfully verified by comparisons with Monte Carlo simulations.

Journal ArticleDOI
TL;DR: Comparing spatial frequency tuning curves for luminance gratings and gratings isolating cone input to the receptive field center reveals indications in some cells of a more complex receptive structure than a simple difference of Gaussians model.
Abstract: Receptive fields of midget ganglion cells and parvocellular lateral geniculate nucleus (LGN) neurons show color-opponent responses because they receive antagonistic input from the middle- and long-wavelength sensitive cones. It has been controversial as to whether this opponency can derive from random connectivity; if receptive field centers of cells near the fovea are cone-specific due to midget morphology, this would confer some degree of color opponency even with random cone input to the surround. A simple test of this mixed surround hypothesis is to compare spatial frequency tuning curves for luminance gratings and gratings isolating cone input to the receptive field center. If tuning curves for luminance gratings were bandpass, then with the mixed surround hypothesis tuning curves for gratings isolating the receptive field center cone class should also be bandpass, but to a lesser extent than for luminance. Tuning curves for luminance, chromatic, and cone-isolating gratings were measured in macaque retinal ganglion cells and LGN cells. We defined and measured a bandpass index to compare luminance and center cone-isolating tuning curves. Midget retinal ganglion cells and parvocellular LGN cells had bandpass indices between 0.1 and 1 with luminance gratings, but the index was usually near 1 (meaning low-pass tuning) when the receptive field center cone class alone was modulated. This is strong evidence for a considerable degree of cone-specific input to the surround. A fraction of midget and parvocellular cells showed evidence of incomplete specificity. Fitting the data with receptive field models revealed considerable intercell variability, with indications in some cells of a more complex receptive structure than a simple difference of Gaussians model.

Journal ArticleDOI
TL;DR: Findings suggest that V1, especially double-opponent cells, may function to extract spatial information from color boundaries correlated with scene-structure cues, such as shadows lit by ambient blue sky juxtaposed with surfaces reflecting sunshine.
Abstract: Colors defined by the two intermediate directions in color space, "orange-cyan" and "lime-magenta," elicit the same spatiotemporal average response from the two cardinal chromatic channels in the lateral geniculate nucleus (LGN). While we found LGN functional magnetic resonance imaging (fMRI) responses to these pairs of colors were statistically indistinguishable, primary visual cortex (V1) fMRI responses were stronger to orange-cyan. Moreover, linear combinations of single-cell responses to cone-isolating stimuli of V1 cone-opponent cells also yielded stronger predicted responses to orange-cyan over lime-magenta, suggesting these neurons underlie the fMRI result. These observations are consistent with the hypothesis that V1 recombines LGN signals into "higher-order" mechanisms tuned to noncardinal color directions. In light of work showing that natural images and daylight samples are biased toward orange-cyan, our findings further suggest that V1 is adapted to daylight. V1, especially double-opponent cells, may function to extract spatial information from color boundaries correlated with scene-structure cues, such as shadows lit by ambient blue sky juxtaposed with surfaces reflecting sunshine.

Journal ArticleDOI
TL;DR: Four different Monte Carlo methods, widely used in the field of tissue optics, that are based on four different ways to build photons' trajectories are considered, providing both theoretical arguments and numerical results showing the statistical equivalence of the four methods.
Abstract: In the field of photon migration in turbid media, different Monte Carlo methods are usually employed to solve the radiative transfer equation. We consider four different Monte Carlo methods, widely used in the field of tissue optics, that are based on four different ways to build photons’ trajectories. We provide both theoretical arguments and numerical results showing the statistical equivalence of the four methods. In the numerical results we compare the temporal point spread functions calculated by the four methods for a wide range of the optical properties in the slab and semi-infinite medium geometry. The convergence of the methods is also briefly discussed.

Journal ArticleDOI
TL;DR: The wavefront sensor in active and adaptive telescopes is usually not in the optical path toward the scientific detector, which may generate additional wavefront aberrations, which have to be separated from the errors due to the telescope optics by a series of measurements taken in the center of the field.
Abstract: The wavefront sensor in active and adaptive telescopes is usually not in the optical path toward the scientific detector. It may generate additional wavefront aberrations, which have to be separated from the errors due to the telescope optics. The aberrations that are not rotationally symmetric can be disentangled from the telescope aberrations by a series of measurements taken in the center of the field, with the wavefront sensor at different orientation angles with respect to the focal plane. This method has been applied at the VLT Survey Telescope on the ESO Paranal observatory.

Journal ArticleDOI
TL;DR: Using azimuthal polarization, it is shown that an amplitude-phase filter allows generation of a subdiffractive dark spot in a prescribed finite area.
Abstract: We compare generation of a dark spot using focusing of beams with azimuthal polarizion, radial polarization with a vortex, and a circular polarization with either a first or second order vortex. By optimization of the amplitude-phase pupil, it is ascertained that azimuthal polarization is the most suitable one to obtain the diffraction bounded dark spot per se whose scalar approximation limit has FWHM=0.29λ. Consequently, for dark spot generation, this polarization plays the role of the radial polarization in creation of the diffraction-limited bright spot. Using azimuthal polarization, it is shown that an amplitude-phase filter allows generation of a subdiffractive dark spot in a prescribed finite area.

Journal ArticleDOI
TL;DR: This paper compares the most relevant solutions proposed hitherto to a novel multirate algorithm using the linear-quadratic-Gaussian (LQG) approach capable of upsampling the correction to further reduce the impact of vibrations.
Abstract: Vibration suppression in astronomical adaptive optics (AO) systems has gathered great attention in the context of next-generation instrumentation for current telescopes and future Extremely Large Telescopes. Laser tomographic AO systems require natural guide stars to measure the low-order modes such as tip–tilt (TT) and TT-anisoplanatism. To increase the sky coverage, the guide stars are often faint, thus requiring lower temporal sampling frequencies to work on a more favorable signal-to-noise regime. Such sampling frequencies can be of the order of, or even lower than, the range of frequencies where vibrations are likely to appear. Ideally, vibrations affecting the low-order modes could be corrected at the higher laser loop frame rate using an upsampling procedure. This paper compares the most relevant solutions proposed hitherto to a novel multirate algorithm using the linear-quadratic-Gaussian (LQG) approach capable of upsampling the correction to further reduce the impact of vibrations. Results from numerical Monte Carlo simulations span a large range of parameters from pure sinusoids to relatively broad peak vibrations, covering the likely-to-be signals in a realistic AO system. The improvement is shown at sampling frequencies from 20 to 800 Hz, including below the vibration itself, in the example of 29.5 Hz on a Thirty Meter Telescope-like scenario. The multirate LQG ensures the least residual for both faint and bright stars for all the peak widths considered based on telemetry from the Keck Observatory.

Journal ArticleDOI
TL;DR: A multipole expansion, based on spherical harmonics, provides an efficient method for calculating the field in the focal region of a lens for radially polarized illumination, or other illumination polarization and phase distributions, including vortex beams.
Abstract: A multipole expansion, based on spherical harmonics, provides an efficient method for calculating the field in the focal region of a lens for radially polarized illumination, or other illumination polarization and phase distributions, including vortex beams. The multipole approach also has the benefit of providing a simple measure of the purity of the longitudinal field mode. The method is also convenient for calculation of fields scattered by particles and calculation of optical trapping forces.