3D Wavefront Research Articles

Refraction/reflection reversal in two-dimensional acoustic metagratings

Unlike acoustic metasurfaces that rely solely on phase gradients, acoustic metagratings (AMs) operate based on both phase gradients and grating diffraction, thus further extending the generalized Snell's law (GSL). In particular, AMs can achieve reversal of refraction and reflection based on the parity of the number of wave propagations inside the AMs. So far, discussions of this GSL extension have largely been applied to one-dimensional periodic AMs, while the designs of two-dimensional (2D) periodic AMs and their performance in three-dimensional (3D) space have been quite limited. Here, we study the GSL extension in 3D space and experimentally demonstrate a series of functional 2D periodic AMs. The designed AMs can achieve sound refraction/reflection under any incidence angle in 3D space, without restrictions to certain critical ranges; adjusting incident angles only enables the reversal of refraction and reflection. Additionally, we demonstrate two types of dual-layer sound lenses based on two AMs, whose reversal of refraction and reflection can be realized by simply attaching or separating the two AMs. Our work paves the way to complex 3D wavefront manipulation of AMs, which may find potential use in practical acoustic devices.

Spin‐Selective Trifunctional Metasurfaces for Deforming Versatile Nondiffractive Beams along the Optical Trajectory

AbstractExploring and taming the diffraction phenomena and divergence of light are foundational to enhancing comprehension of nature and developing photonic technologies. Despite the numerous types of nondiffraction beam generation technologies, the 3D deformation and intricate wavefront shaping of structures during propagation have yet to be studied through the lens of nanophotonic devices. Herein, the dynamic conversion of a circular Airy beam (CAB) to a Bessel beam with a single‐layer spin‐selective metasurface is demonstrated. This spatial transformation arises from the interplay of 1D local and 2D global phases, facilitating the 3D control of non‐diffractive light fields. An additional overall phase gradient and orbital angular momentum are introduced, which effectively altering the propagation direction and transverse fields of complex amplitude beams along the optical path. The manifested samples exhibit superior defect resistance, laying a crucial application in micro/nanolithography technologies. This approach expands the in‐plane spin‐selective mechanism and leverages the out‐of‐plane propagation dimension, allowing for integrated high‐resolution imaging, on‐chip optical micromanipulation, and micro/nanofabrication within a versatile nanophotonic platform.

Wavefront Modulation and Beam Shaping into Arbitrary Three-Dimensional Intensity Distributions

In this study the methods of three-dimensional (3D) wavefront intensity modulation by employing contrast-inverted holography, previously introduced as Gabor inverted holography, are further investigated. The present study provides the recipes for creating 3D wavefront intensity modulations using phase-only and amplitude-only modulators and compares the results. The 3D wavefront modulation using spherical waves is also demonstrated, and the miniaturization of 3D intensity beams is discussed; it is shown that both the resolution and the size of the created 3D structures are ultimately given by the wavelength of the employed radiation. The manuscript also addresses the quality of the formed 3D intensity curves and determines the parameters that provide the best smooth appearance of the 3D curves. The presented methods of 3D intensity wavefront modulation can be realized for all kinds of waves: light, X-ray, electron, etc, provided the modulator can be manufactured for the corresponding wavelength. The methods of 3D intensity wavefront modulation can be applied in various techniques: lithography, micro-robotics, particle trapping, etc.

Moiré-Ronchigram analysis applied in the characterization of aberration surfaces and optical surface parameters from 3D wavefronts.

In this work, we propose the construction of the Moiré-Ronchigram as a pattern obtained from simple Ronchigrams, superposed and with a slight degree of rotation, to obtain a 3D wavefront. The aberrations in the obtained wavefront are fitted to the aberration Zernike polynomial of degree 12 in order to obtain a polynomial expression of a high degree of efficiency, which is subsequently compared with the general equation of quadrics (considering the coefficients accordingly) to obtain essential parameters (as focus, eccentricity, or F number) that allow us to understand and manipulate the optical elements under test optimally. Finally, we compare our results analyzing optical surfaces to make a weak statistic of our proposed method.

Reliability of data-driven wavefront attributes in laterally heterogeneous media

3D wavefront attributes play a major role in many processing steps, such as prestack data enhancement, diffraction separation, and wavefront tomography. For the determination of the 3D wavefront attributes, various stacking operators can be used by adopting semblance optimization. These operators are derived for laterally homogeneous media. In praxis, however, they are applied in real geologic environments with even strong lateral velocity variations such as salt structures. This leads to the question of the quality of the 3D wavefront attributes using these operators when determined in the presence of strong lateral velocity changes. We compared the 3D wavefront attributes determined by 3D common-reflection-surface (CRS) operator (called data-driven wavefront attributes) with the 3D wavefront attributes computed by 3D kinematic and dynamic ray tracing (called model-driven wavefront attributes). For the determination of the 3D CRS wavefront attributes, we have developed a global optimization scheme based on differential evolution. Reflection seismic data of the laterally heterogeneous 3D SEG C3WA salt model are considered, and the model-driven wavefront attributes are computed for a smoothed version of the 3D SEG salt model. The comparison reveals that the wavefront attributes for the normal-incidence-point ray indicate a very good match not only in areas of mild lateral velocity variation but even in regions with strong lateral velocity variations. Approximately 80%–90% of the total picks indicate the good match with a relative error of less than 10% when a semblance threshold of 0.1 is considered in the automatic picking process. This confirms the validity of the determined wavefront attributes even in the presence of strong lateral velocity changes. Using a higher semblance threshold in the automatic picking leads to fewer picks but with an even better match between model- and data-driven wavefront attributes.

3D wavefront attribute determination and conflicting dip processing

The knowledge of 3D wavefront attributes allows many important applications, such as stacking, 5D interpolation, 3D diffraction separation and imaging, and 3D wavefront tomography, just to name a few. For the determination of wavefront attributes, we use the common-reflection-surface (CRS) operator. We adopt a simultaneous search for the determination of wavefront attributes and combine it with conflicting dip processing. For the simultaneous search, we compare three heuristic global optimization algorithms such as particle swarm optimization (PSO), genetic algorithm (GA), and differential evolution (DE). For conflicting dip processing, a dip angle decomposition method for the probed sample is introduced and the simultaneous search is independently performed in specified dip ranges to individually obtain attributes and semblance for each range. Results for the laterally heterogeneous 3D SEG C3WA data indicate that DE has superior performance to determine the 3D wavefront attributes when compared with PSO, GA, and the conventional pragmatic approach because a higher semblance and an improved set of wavefront attributes are achieved. A comparison of the data-driven wavefront attributes obtained from the DE with the model-driven wavefront attributes computed by kinematic and dynamic ray tracing reveals the validity of the data-driven wavefront attributes. Combining the simultaneous search with conflicting dip processing for the 3D CRS stack further improved reflected energy and diffraction details when compared with results without simultaneous search and/or conflicting dip processing.

3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching

ABSTRACTSubsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration.

Measurement of three-dimensional wavefronts using the Ichikawa-Lohmann-Takeda solution to the irradiance transport equation.

In this paper, we use the irradiance transport equation and the Fourier transform-based experimental solution given by Ichikawa-Lohmann-Takeda. We analyze experimental factors such as the digital filter, the introduced error for the rotation and period of the Ronchi ruling, and a new method is demonstrated for the measurement of 3D wavefront information.

Tunable liquid crystal cylindrical micro-optical array for aberration compensation

A tunable aberration compensation device for rectangular micro-optical systems is proposed and demonstrated. This device, which is based in nematic liquid crystal and a micro-electrode structure, forms gradients in the index of refraction as a function of voltage. We have developed a fringe skeletonizing application in order to extract the 3D wavefront from an interference pattern. This software tool obtains the optical aberrations using Chebyshev polynomials. By using phase shifted electrical signals the aberrations can be controlled independently. A complete independent control over the spherical and coma aberration has been demonstrated. Also, an independent control over the astigmatism aberration has been demonstrated in a broad range. This device has promising applications where aberration compensation is required. The independent compensation achieved for some coefficients, such as astigmatism for example, is more than 2.4 waves.

Analysis of modal tomography for three-dimensional wavefront sensing of atmosphere turbulence

Three-dimensional (3D) wavefront sensing of atmosphere turbulence is a key step to realize multi-conjugate adaptive optics technology. In this article, model tomography, which is one of the most important algorithms for 3D wavefront sensing, is analyzed in theory, the principle limitation of model tomography is proposed. Based on this view the reason for tomography error is discussed, and the simulation results for different errors are produced finally. The analysis shows that part Zernike model basis is used as a new basis of polynomial decomposition in model tomography, which leads to modal coupling error and aliasing error. The un-correlation of part Zernike model basis is the prerequisite to avoid modal aliasing error, but modal coupling error cannot be removed and we can only restrain its influence. Combined with simulation result a method with large field of view (FOV) sensing and small FOV reconstruction is proposed and gives a good performance to control modal coupling error.

3D imaging and wavefront sensing with a plenoptic objective

3D imaging and wavefront sensing with a plenoptic objective

Extracting Intramural Wavefront Orientation from Optical Upstroke Shapes in Whole Hearts

Information about intramural propagation of electrical excitation is crucial to understanding arrhythmia mechanisms in thick ventricular muscle. There is currently a controversy over whether it is possible to extract such information from the shape of the upstroke in optical mapping recordings. We show that even in the complex geometry of a whole guinea pig heart, optical upstroke morphology reveals the 3D wavefront orientation near the surface. To characterize the upstroke morphology, we use V F ∗ , the fractional level at which voltage-sensitive fluorescence, V F, has maximal time derivative. Low values of V F * ( ∼ 0.2 ) indicate a wavefront moving away from the surface, high values of V F * ( ∼ 0.4 ) a wavefront moving toward the surface, and intermediate values of V F * ( ∼ 0.6 ) a wavefront moving parallel to the surface. We further performed computer simulations using Luo-Rudy II electrophysiology and a simplified 3D geometry. The simulated V * maps for free wall and apical stimulations as well as for sinus rhythm are in good quantitative agreement with the averaged experimental results. Furthermore, computer simulations show that the effect of the curvature of the heart on wave propagation is negligible.

Advancements in the Interferometric Measurements of Real Time Finishing Birefringent Filter's Crystal Plates

The finishing of birefringent plates consists of two processes: polishing and evaluation of a surface, which have been performed separately till now. The purpose of this work is achieving of high accuracy of the evaluation and machining of the plane-parallel plates from birefringent crystals, in particular of crystal plates of birefringent filters during their finishing. The developed process combines evaluation and polishing in an interactive way. We have found modes of treatment, shape of polisher, have designed interferometer, with a mirror arranged in polisher. Visual checking of optical thickness comparatively with reference plate was carried out using the interference fringes of equal birefringence, and checking of an optical wedge - by interference rings of an equal inclination. The automated processing of TV camera interference fringes was impossible, because of gaps of interference fringes on polishing cells above the mirror. Therefore a special software was developed for processing of a complex fringe pattern interferogram. Software FastInterf uses furrier analysis technique which allows to process an interferogram with multiply gaps. Interferograms are registered by a high resolution TV camera (1280 ×1024). Automatic processing of a fringe interferogram using FastInterf software takes less then one second. The influence of gaps is excluded, and the flat field is taken into account. Software provides full 3D surface and wavefront maps. Aberration analysis of a wavefront gives information on thickness of a plate comparatively with a reference one, optical wedge of plate and azimuth of an inclination of wave front. Moreover, software provides a control of surface quality. The measuring device, features of the software are described and process of interferometric evaluation during polishing is illustrated.

3D imaging: wave front sensing utilizing a birefringent crystal

Utilizing a rotatable thin birefringent crystal and two polarizers, the three-dimensional (3D) wave front of any light wave can be measured conveniently. In first experiments phase-front distortions as small as 15 &mgr;m are detected with a dynamic range of 3 mm and a spatial resolution of 50 &mgr;m. Such a dynamic range and spatial resolution exceed the performance of conventional wave front sensors of, e.g., the Shack-Hartmann type, significantly. Furthermore, the new system is rather insensitive against mechanical instabilities in opposition to interferometric and holographic techniques.