An Algorithm for Fresnel Diffraction Computing Based on Fractional Fourier Transform

The fractional Fourier transform (FrFT) is used for the solution of the diffraction integral in optics. A scanning approach is proposed for finding the optimal FrFT order. In this way, the process of diffraction computing is speeded up. The basic algorithm and the intermediate results at each stage are demonstrated.

Formation of soliton trains in Bose-Einstein condensates as a nonlinear Fresnel diffraction of matter waves

The problem of generation of atomic soliton trains in elongated Bose-Einstein condensates is considered in framework of Whitham theory of modulations of nonlinear waves. Complete analytical solution is presented for the case when the initial density distribution has sharp enough boundaries. In this case the process of soliton train formation can be viewed as a nonlinear Fresnel diffraction of matter waves. Theoretical predictions are compared with results of numerical simulations of one- and three-dimensional Gross-Pitaevskii equation and with experimental data on formation of Bose-Einstein bright solitons in cigar-shaped traps. (C) 2003 Elsevier B.V. All rights reserved.

Point diffraction interferometer with a liquid crystal monopixel

In this work a novel point diffraction interferometer based on a variable liquid crystal wave plate (LCWP) has been implemented. The LCWP consists of a 3x3 cm2 monopixel cell with parallel alignment. The monopixel cell was manufactured such that the electrode covers the entire surface except in a centered circular area of 50 μm of diameter. This circle acts as a point perturbation which diffracts the incident wave front giving rise to a spherical reference wave. By applying a voltage to the LCWP we can change the phase of the wave front that passes through the monopixel, except at the center. Phase shifting techniques are used in order to calculate the amplitude and phase distribution of the object wave front. The system allows a digital hologram to be obtained, and by using the Fresnel diffraction integral it is possible to digitally reconstruct the different planes that constitute the three dimensional object.

In-situ measurement and characterization of cloud particles using digital in-line holography

Satellite measurement validations, climate models, atmospheric radiative transfer models and cloud models, all depend on accurate measurements of cloud particle size distributions, number densities, spatial distributions, and other parameters relevant to cloud microphysical processes. And many airborne instruments designed to measure size distributions and concentrations of cloud particles have large uncertainties in measuring number densities and size distributions of small ice crystals. HOLODEC (Holographic Detector for Clouds) is a new instrument that does not have many of these uncertainties and makes possible measurements that other probes have never made. The advantages of HOLODEC are inherent to the holographic method. In this dissertation, I describe HOLODEC, its in-situ measurements of cloud particles, and the results of its test flights. I present a hologram reconstruction algorithm that has a sample spacing that does not vary with reconstruction distance. This reconstruction algorithm accurately reconstructs the field to all distances inside a typical holographic measurement volume as proven by comparison with analytical solutions to the Huygens-Fresnel diffraction integral. It is fast to compute, and has diffraction limited resolution. Further, described herein is an algorithm that can find the position along the optical axis of small particles as well as large complex-shaped particles. I explain an implementation of these algorithms that is an efficient, robust, automated program that allows us to process holograms on a computer cluster in a reasonable time. I show size distributions and number densities of cloud particles, and show that they are within the uncertainty of independent measurements made with another measurement method. The feasibility of another cloud particle instrument that has advantages over new standard instruments is proven. These advantages include a unique ability to detect shattered particles using three-dimensional positions, and a sample volume size that does not vary with particle size or airspeed. It also is able to yield two-dimensional particle profiles using the same measurements.

Orthogonal chirp division multiplexing for coherent optical fiber communications

In this paper, we propose an orthogonal chirp division multiplexing (OCDM) technique for coherent optical communication. OCDM is the principle of orthogonally multiplexing a group of linear chirped waveforms for high-speed data communication, achieving the maximum spectral efficiency (SE) for chirp spread spectrum, in a similar way as the orthogonal frequency division multiplexing (OFDM) does for frequency division multiplexing. In the coherent optical (CO)-OCDM, Fresnel transform formulates the synthesis of the orthogonal chirps; discrete Fresnel transform (DFnT) realizes the CO-OCDM in the digital domain. As both the Fresnel and Fourier transforms are trigonometric transforms, the CO-OCDM can be easily integrated into the existing CO-OFDM systems. Analyses and numerical results are provided to investigate the transmission of CO-OCDM signals over optical fibers. Moreover, experiments of 36-Gbit/s CO-OCDM signal are carried out to validate the feasibility and confirm the analyses. It is shown that the CO-OCDM can effectively compensate the dispersion and is more resilient to fading and noise impairment than OFDM.

Modeling diffraction in free-space optical interconnects by the mode expansion method

Free-space optical interconnects (FSOIs), made up of dense arrays of vertical-cavity surface-emitting lasers, photodetectors and microlenses can be used for implementing high-speed and high-density communication links, and hence replace the inferior electrical interconnects. A major concern in the design of FSOIs is minimization of the optical channel cross talk arising from laser beam diffraction. In this article we introduce modifications to the mode expansion method of Tanaka et al. [IEEE Trans. Microwave Theory Tech. MTT-20, 749 (1972)] to make it an efficient tool for modelling and design of FSOIs in the presence of diffraction. We demonstrate that our modified mode expansion method has accuracy similar to the exact solution of the Huygens-Kirchhoff diffraction integral in cases of both weak and strong beam clipping, and that it is much more accurate than the existing approximations. The strength of the method is twofold: first, it is applicable in the region of pronounced diffraction (strong beam clipping) where all other approximations fail and, second, unlike the exact-solution method, it can be efficiently used for modelling diffraction on multiple apertures. These features make the mode expansion method useful for design and optimization of free-space architectures containing multiple optical elements inclusive of optical interconnects and optical clock distribution systems. (C) 2003 Optical Society of America.

Simultaneous observations of ionospheric quasiperiodic scintillations from short and long meridional baselines using VHF transmissions from transit satellites

For the first time it was possible to observe regular quasiperiodic scintillations (QPS) in VHF radio-satellite transmissions from orbiting satellites simultaneously at short (2.1 km) and long (121 km) meridional baselines in the vicinity of a typical mid-latitude station (Brisbane; 27.5degreesS and 152.9degreesE geog. and 35.6degrees invar.lat.), using three sites (St. Lucia-S, Taringa-T in Brisbane and Boreen Pt.-B, north of Brisbane). A few pronounced quasiperiodic (QP) events were recorded showing unambiguous regular structures at the sites which made it possible to deduce a time displacement of the regular fading minimum at S, T and B. The QP structure is highly dependent on the geometry of the ray-path from a satellite to the observer which is manifested as a change of a QP event from symmetrical to non-symmetrical for stations separated by 2.1 km, and to a radical change in the structure of the event over a distance of 121 km. It is suggested the short-duration intense QP events are due to a Fresnel diffraction (or a reflection mechanism) of radio-satellite signals by a single ionospheric irregularity in a form of an ellipsoid with a large ionization gradient along the major axis. The structure of a QP event depends on the angle of viewing of the irregular blob from a radio-satellite. In view of this it is suggested that the reported variety of the ionization formation, responsible for different types of QPS, is only apparent but not real. (C) 2003 Elsevier Science Ltd. All rights reserved.

Dynamics of bright matter wave solitons in a Bose-Einstein condensate

Recent experimental and theoretical advances in the creation and description of bright matter wave solitons are reviewed. Several aspects are taken into account, including the physics of soliton train formation as the nonlinear Fresnel diffraction, soliton-soliton interactions, and propagation in the presence of inhomogeneities. The generation of stable bright solitons by means of Feshbach resonance techniques is also discussed. © World Scientific Publishing Company.

Diffraction theory of Fresnel lenses encoded in low-resolution devices

A mathematical model that describes the behavior of low-resolution Fresnel lenses encoded in any low-resolution device (e.g., a spatial light modulator) is developed. The effects of low-resolution codification, such the appearance of new secondary lenses, are studied for a general case. General expressions for the phase of these lenses are developed, showing that each lens behaves as if it were encoded through all pixels of the low-resolution device. Simple expressions for the light distribution in the focal plane and its dependence on the encoded focal length are developed and commented on in detail. For a given codification device an optimum focal length is found for best lens performance. An optimization method for codification of a single lens with a short focal length is proposed.

Diffraction theory of optimized low-resolution Fresnel encoded lenses

A mathematical model describing the behavior of low-resolution Fresnel encoded lenses (LRFEL's) encoded in any low-resolution device (e.g., a spatial light modulator) has recently been developed. From this model, an LRFEL with a short focal length was optimized by our imposing the maximum intensity of light onto the optical axis. With this model, analytical expressions for the light-amplitude distribution, the diffraction efficiency, and the frequency response of the optimized LRFEL's are derived.

On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

A Nyström method for a class of integral equations on the real line with applications to scattering by diffraction gratings and rough surfaces

We propose a Nystr¨om/product integration method for a class of second kind integral equations on the real line which arise in problems of two-dimensional scalar and elastic wave scattering by unbounded surfaces. Stability and convergence of the method is established with convergence rates dependent on the smoothness of components of the kernel. The method is applied to the problem of acoustic scattering by a sound soft one-dimensional surface which is the graph of a function f, and superalgebraic convergence is established in the case when f is infinitely smooth. Numerical results are presented illustrating this behavior for the case when f is periodic (the diffraction grating case). The Nystr¨om method for this problem is stable and convergent uniformly with respect to the period of the grating, in contrast to standard integral equation methods for diffraction gratings which fail at a countable set of grating periods.

On the solvability of second kind integral equations on the real line

This paper considers general second kind integral equations of the form(in operator form φ − kφ = ψ), where the functions k and ψ are assumed known, with ψ ∈ Y, the space of bounded continuous functions on R, and k such that the mapping s → k(s, · ), from R to L1(R), is bounded and continuous. The function φ ∈ Y is the solution to be determined. Conditions on a set W ⊂ BC(R, L1(R)) are obtained such that a generalised Fredholm alternative holds: If W satisfies these conditions and I − k is injective for all k ∈ W then I − k is also surjective for all k ∈ W and, moreover, the inverse operators (I − k) − 1 on Y are uniformly bounded for k ∈ W. The approximation of the kernel in the integral equation by a sequence (kn) converging in a weak sense to k is also considered and results on stability and convergence are obtained. These general theorems are used to establish results for two special classes of kernels: k(s, t) = κ(s − t)z(t) and k(s, t) = κ(s − t)λ(s − t, t), where κ ∈ L1(R), z ∈ L∞(R), and λ ∈ BC((R\{0}) × R). Kernels of both classes arise in problems of time harmonic wave scattering by unbounded surfaces. The general integral equation results are here applied to prove the existence of a solution for a boundary integral equation formulation of scattering by an infinite rough surface and to consider the stability and convergence of approximation of the rough surface problem by a sequence of diffraction grating problems of increasingly large period.

Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade

O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente.

Migração Kirchhoff pré-empilhamento em profundidade modificada usando o operador de feixes gaussianos

A teoria dos feixes gaussianos foi introduzida na literatura sísmica no início dos anos 80 por pesquisadores russos e tchecos, e foi originalmente utilizada no cálculo do campo de ondas eletromagnéticas, baseado na teoria escalar da difração. Na teoria dos feixes gaussianos, o campo de ondas sísmicas é obtido por uma integral, cujo o integrando é constituído de duas partes, a saber: (1) as amplitudes dos campos das ondas na vizinhança do ponto de observação e (2) a função fase de cada um desses campos de ondas, que neste caso é representada por um tempo de trânsito paraxial complexo. Como ferramenta de imageamento, mais precisamente como operador de migração, os primeiros trabalhos usando feixes gaussianos datam do final da década de 80 e início dos anos 90. A regularidade dos campos de ondas descritos pelos feixes gaussianos, além de sua alta precisão em regiões singulares do modelo de velocidades, tornaram o uso de feixes gaussianos como uma alternativa híbrida viável para a migração. Nesse trabalho, unimos a flexibilidade da migração tipo Kirchhoff em profundidade em verdadeira amplitude com a regularidade da descrição do campo de ondas, representado pela sobreposição de feixes gaussianos. Como forma de controlar de forma estável quantidades usadas na construção de feixes gaussianos, utilizamos informações advindas do volume de Fresnel, mais precisamente a zona de Fresnel ao redor do ponto de reflexão e a zona de Fresnel projetada, localizada ao redor do ponto de registro do sismograma e cuja a informação se encontra nas curvas de reflexão de dados sísmico. Nosso processo de migração pode ser chamado como uma migração Kirchhoff em verdadeira amplitude usando um operador de feixes gaussianos.