Relationship between the beam propagation method and linear canonical and fractional Fourier transforms.

The beam propagation method (BPM) can be viewed as a chain of alternating convolutions and multiplications, as filtering operations alternately in the space and frequency domains or as multiplication operations sandwiched between linear canonical or fractional Fourier transforms. These structures provide alternative models of inhomogeneous media and potentially allow mathematical tools and algorithms associated with these transforms to be applied to the BPM. As an example, in the case where quadratic approximation is possible, it is shown that the BPM can be represented as a single LCT system, leading to significantly faster computation of the output field.

Linear algebraic theory of partial coherence: continuous fields and measures of partial coherence.

This work presents a linear algebraic theory of partial coherence for optical fields of continuous variables. This approach facilitates use of linear algebraic techniques and makes it possible to precisely define the concepts of incoherence and coherence in a mathematical way. We have proposed five scalar measures for the degree of partial coherence. These measures are zero for incoherent fields, unity for fully coherent fields, and between zero and one for partially coherent fields.

Beyond Nyquist sampling: a cost-based approach.

A sampling-based framework for finding the optimal representation of a finite energy optical field using a finite number of bits is presented. For a given bit budget, we determine the optimum number and spacing of the samples in order to represent the field with as low error as possible. We present the associated performance bounds as trade-off curves between the error and the cost budget. In contrast to common practice, which often treats sampling and quantization separately, we explicitly focus on the interplay between limited spatial resolution and limited amplitude accuracy, such as whether it is better to take more samples with lower amplitude accuracy or fewer samples with higher accuracy. We illustrate that in certain cases sampling at rates different from the Nyquist rate is more efficient.

Phase-space window and degrees of freedom of optical systems with multiple apertures.

We show how to explicitly determine the space-frequency window (phase-space window) for optical systems consisting of an arbitrary sequence of lenses and apertures separated by arbitrary lengths of free space. If the space-frequency support of a signal lies completely within this window, the signal passes without information loss. When it does not, the parts that lie within the window pass and the parts that lie outside of the window are blocked, a result that is valid to a good degree of approximation for many systems of practical interest. Also, the maximum number of degrees of freedom that can pass through the system is given by the area of its space-frequency window. These intuitive results provide insight and guidance into the behavior and design of systems involving multiple apertures and can help minimize information loss.

Teaching science, technology, and society to engineering students: a sixteen year journey.

The course Science, Technology, and Society is taken by about 500 engineering students each year at Bilkent University, Ankara. Aiming to complement the highly technical engineering programs, it deals with the ethical, social, cultural, political, economic, legal, environment and sustainability, health and safety, reliability dimensions of science, technology, and engineering in a multidisciplinary fashion. The teaching philosophy and experiences of the instructor are reviewed. Community research projects have been an important feature of the course. Analysis of teaching style based on a multi-dimensional model is given. Results of outcome measurements performed for ABET assessment are provided. Challenges and solutions related to teaching a large class are discussed.

Representation of optical fields using finite numbers of bits.

We consider the problem of representation of a finite-energy optical field, with a finite number of bits. The optical field is represented with a finite number of uniformly spaced finite-accuracy samples (there is a finite number of amplitude levels that can be reliably distinguished for each sample). The total number of bits required to encode all samples constitutes the cost of the representation. We investigate the optimal number and spacing of these samples under a total cost budget. Our framework reveals the trade-off between the number, spacing, and accuracy of the samples. When we vary the cost budget, we obtain trade-off curves between the representation error and the cost budget. We also discuss the effect of degree of coherence of the field.

Fundamental structure of Fresnel diffraction: longitudinal uniformity with respect to fractional Fourier order.

Fresnel integrals corresponding to different distances can be interpreted as scaled fractional Fourier transformations observed on spherical reference surfaces. Transverse samples can be taken on these surfaces with separation that increases with propagation distance. Here, we are concerned with the separation of the spherical reference surfaces along the longitudinal direction. We show that these surfaces should be equally spaced with respect to the fractional Fourier transform order, rather than being equally spaced with respect to the distance of propagation along the optical axis. The spacing should be of the order of the reciprocal of the space-bandwidth product of the signals. The space-dependent longitudinal and transverse spacings define a grid that reflects the structure of Fresnel diffraction.

Fundamental structure of Fresnel diffraction: natural sampling grid and the fractional Fourier transform.

Fresnel integrals corresponding to different distances can be interpreted as scaled fractional Fourier transformations observed on spherical reference surfaces. We show that by judiciously choosing sample points on these curved reference surfaces, it is possible to represent the diffracted signals in a nonredundant manner. The change in sample spacing with distance reflects the structure of Fresnel diffraction. This sampling grid also provides a simple and robust basis for accurate and efficient computation, which naturally handles the challenges of sampling chirplike kernels.

Synthesis of three-dimensional light fields with binary spatial light modulators.

Computation of a binary spatial light modulator (SLM) pattern that generates a desired light field is a challenging quantization problem for which several algorithms have been proposed, mainly for far-field or Fourier plane reconstructions. We study this problem assuming that the desired light field is synthesized within a volumetric region in the non-far-field range after free space propagation from the SLM plane. We use Fresnel and Rayleigh-Sommerfeld scalar diffraction theories for propagation of light. We show that, when the desired field is confined to a sufficiently narrow region of space, the ideal gray-level complex-valued SLM pattern generating it becomes sufficiently low pass (oversampled) so it can be successfully halftoned into a binary SLM pattern by solving two decoupled real-valued constrained halftoning problems. Our simulation results indicate that, when the synthesis region is considered, the binary SLM is indistinguishable from a lower resolution full complex gray-level SLM. In our approach, free space propagation related computations are done only once at the beginning, and the rest of the computation time is spent on carrying out standard image halftoning.

Full-complex amplitude modulation with binary spatial light modulators.

Imperfections and nonrobust behavior of practical multilevel spatial light modulators (SLMs) degrade the performance of many proposed full-complex amplitude modulation schemes. We consider the use of more robust binary SLMs for this purpose. We propose a generic method, by which, out of K binary (or 1 bit) SLMs of size M×N, we effectively create a new 2(K)-level (or K bit) SLM of size M×N. The method is a generalization of the well-known concepts of bit plane representation and decomposition for ordinary gray scale digital images and relies on forming a properly weighted superposition of binary SLMs. When K is sufficiently large, the effective SLM can be regarded as a full-complex one. Our method is as efficient as possible from an information theoretical perspective. A 4f system is discussed as a possible optical implementation. This 4f system also provides a means for eliminating the undesirable higher diffraction orders. The components of the 4f system can easily be customized for different production technologies.

Equivalence of linear canonical transform domains to fractional Fourier domains and the bicanonical width product: a generalization of the space-bandwidth product.

Linear canonical transforms (LCTs) form a three-parameter family of integral transforms with wide application in optics. We show that LCT domains correspond to scaled fractional Fourier domains and thus to scaled oblique axes in the space-frequency plane. This allows LCT domains to be labeled and ordered by the corresponding fractional order parameter and provides insight into the evolution of light through an optical system modeled by LCTs. If a set of signals is highly confined to finite intervals in two arbitrary LCT domains, the space-frequency (phase space) support is a parallelogram. The number of degrees of freedom of this set of signals is given by the area of this parallelogram, which is equal to the bicanonical width product but usually smaller than the conventional space-bandwidth product. The bicanonical width product, which is a generalization of the space-bandwidth product, can provide a tighter measure of the actual number of degrees of freedom, and allows us to represent and process signals with fewer samples.

Fast and accurate algorithm for the computation of complex linear canonical transforms.

A fast and accurate algorithm is developed for the numerical computation of the family of complex linear canonical transforms (CLCTs), which represent the input-output relationship of complex quadratic-phase systems. Allowing the linear canonical transform parameters to be complex numbers makes it possible to represent paraxial optical systems that involve complex parameters. These include lossy systems such as Gaussian apertures, Gaussian ducts, or complex graded-index media, as well as lossless thin lenses and sections of free space and any arbitrary combinations of them. Complex-ordered fractional Fourier transforms (CFRTs) are a special case of CLCTs, and therefore a fast and accurate algorithm to compute CFRTs is included as a special case of the presented algorithm. The algorithm is based on decomposition of an arbitrary CLCT matrix into real and complex chirp multiplications and Fourier transforms. The samples of the output are obtained from the samples of the input in approximately N log N time, where N is the number of input samples. A space-bandwidth product tracking formalism is developed to ensure that the number of samples is information-theoretically sufficient to reconstruct the continuous transform, but not unnecessarily redundant.

Fast and accurate computation of two-dimensional non-separable quadratic-phase integrals.

We report a fast and accurate algorithm for numerical computation of two-dimensional non-separable linear canonical transforms (2D-NS-LCTs). Also known as quadratic-phase integrals, this class of integral transforms represents a broad class of optical systems including Fresnel propagation in free space, propagation in graded-index media, passage through thin lenses, and arbitrary concatenations of any number of these, including anamorphic/astigmatic/non-orthogonal cases. The general two-dimensional non-separable case poses several challenges which do not exist in the one-dimensional case and the separable two-dimensional case. The algorithm takes approximately N log N time, where N is the two-dimensional space-bandwidth product of the signal. Our method properly tracks and controls the space-bandwidth products in two dimensions, in order to achieve information theoretically sufficient, but not wastefully redundant, sampling required for the reconstruction of the underlying continuous functions at any stage of the algorithm. Additionally, we provide an alternative definition of general 2D-NS-LCTs that shows its kernel explicitly in terms of its ten parameters, and relate these parameters bidirectionally to conventional ABCD matrix parameters.

Efficient computation of quadratic-phase integrals in optics.

We present a fast NlogN time algorithm for computing quadratic-phase integrals. This three-parameter class of integrals models propagation in free space in the Fresnel approximation, passage through thin lenses, and propagation in quadratic graded-index media as well as any combination of any number of these and is therefore of importance in optics. By carefully managing the sampling rate, one need not choose N much larger than the space-bandwidth product of the signals, despite the highly oscillatory integral kernel. The only deviation from exactness arises from the approximation of a continuous Fourier transform with the discrete Fourier transform. Thus the algorithm computes quadratic-phase integrals with a performance similar to that of the fast-Fourier-transform algorithm in computing the Fourier transform, in terms of both speed and accuracy.

Complex signal recovery from multiple fractional Fourier-transform intensities.

The problem of recovering a complex signal from the magnitudes of any number of its fractional Fourier transforms at any set of fractional orders is addressed. This problem corresponds to the problem of phase retrieval from the transverse intensity profiles of an optical field at arbitrary locations in an optical system involving arbitrary concatenations of lenses and sections of free space. The dependence of the results on the number of orders, their spread, and the noise is investigated. Generally, increasing the number of orders improves the results, but with diminishing return beyond a certain point. Selecting the measurement planes such that their fractional orders are well separated or spread as much as possible also leads to better results.

Fractional free space, fractional lenses, and fractional imaging systems.

Continuum extensions of common dual pairs of operators are presented and consolidated, based on the fractional Fourier transform. In particular, the fractional chirp multiplication, fractional chirp convolution, and fractional scaling operators are defined and expressed in terms of their common nonfractional special cases, revealing precisely how they are interpolations of their conventional counterparts. Optical realizations of these operators are possible with use of common physical components. These three operators can be interpreted as fractional lenses, fractional free space, and fractional imaging systems, respectively. Any optical system consisting of an arbitrary concatenation of sections of free space and thin lenses can be interpreted as a fractional imaging system with spherical reference surfaces. As a special case, a system departing from the classical single-lens imaging condition can be interpreted as a fractional imaging system.

Linear algebraic theory of partial coherence: discrete fields and measures of partial coherence.

A linear algebraic theory of partial coherence is presented that allows precise mathematical definitions of concepts such as coherence and incoherence. This not only provides new perspectives and insights but also allows us to employ the conceptual and algebraic tools of linear algebra in applications. We define several scalar measures of the degree of partial coherence of an optical field that are zero for full incoherence and unity for full coherence. The mathematical definitions are related to our physical understanding of the corresponding concepts by considering them in the context of Young's experiment.