skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00PM ET on Friday, December 15 until 2:00 AM ET on Saturday, December 16 due to maintenance. We apologize for the inconvenience.

Title: A Computational Design Framework for Efficient, Fabrication Error-Tolerant, Planar THz Diffractive Optical Elements

We demonstrate ultra-thin (1.5-3λ0), fabrication-error tolerant efficient diffractive terahertz (THz) optical elements designed using a computer-aided optimization-based search algorithm. The basic operation of these components is modeled using scalar diffraction of electromagnetic waves through a pixelated multi-level 3D-printed polymer structure. Through the proposed design framework, we demonstrate the design of various ultrathin planar THz optical elements, namely (i) a high Numerical Aperture (N.A.), broadband aberration rectified spherical lens (0.1 THz–0.3 THz), (ii) a spectral splitter (0.3 THz–0.6 THz) and (iii) an on-axis broadband transmissive hologram (0.3 THz–0.5 THz). Such an all-dielectric computational design-based approach is advantageous against metallic or dielectric metasurfaces from the perspective that it incorporates all the inherent structural advantages associated with a scalar diffraction based approach, such as (i) ease of modeling, (ii) substrate-less facile manufacturing, (iii) planar geometry, (iv) high efficiency along with(v)broadband operation, (vi) area scalability and (vii) fabrication error-tolerance. With scalability and error tolerance being two major bottlenecks of previous design strategies. This work is therefore, a significant step towards the design of THz optical elements by bridging the gap between structural and computational design i.e. through a hybrid design-based approach enabling considerably less computational resources than the previous state of the art. Furthermore, the approach used herein can be expanded to a myriad of optical elements at any wavelength regime.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. This work discusses the design and fabrication of a dual-plane terahertz (THz) hologram and an extended-depth-of-focus THz diffractive lens. The dual-plane THz hologram consists of 50 × 50 diffractive optical elements with identical element pixel size 1×1 mm, and the extended-depth-of-focus THz diffractive lens is designed with 25 concentric rings with identical ring width of 1 mm, resulting in same device dimension 50 mm × 50 mm. The height of the hologram pixels and concentric rings of the diffractive lens are optimized by nonlinear optimization algorithms with scalar diffraction theory based on Ray-Sommerfeld diffraction equation. Finite-Difference Time-Domain (FDTD) simulation results agree with optimization results obtained from the scalar diffraction theory for both the THz hologram and the THz diffractive lens. The demonstrated experimental results show that the proposed THz hologram and THz diffractive lens can generate the desired diffraction patterns. These diffractive structures have the potential to be applied in areas such as THz imaging, data storage, and displays.

    more » « less
  2. null (Ed.)
    Abstract Aerosol jet printing (AJP) is a direct-write additive manufacturing technique, which has emerged as a high-resolution method for the fabrication of a broad spectrum of electronic devices. Despite the advantages and critical applications of AJP in the printed-electronics industry, AJP process is intrinsically unstable, complex, and prone to unexpected gradual drifts, which adversely affect the morphology and consequently the functional performance of a printed electronic device. Therefore, in situ process monitoring and control in AJP is an inevitable need. In this respect, in addition to experimental characterization of the AJP process, physical models would be required to explain the underlying aerodynamic phenomena in AJP. The goal of this research work is to establish a physics-based computational platform for prediction of aerosol flow regimes and ultimately, physics-driven control of the AJP process. In pursuit of this goal, the objective is to forward a three-dimensional (3D) compressible, turbulent, multiphase computational fluid dynamics (CFD) model to investigate the aerodynamics behind: (i) aerosol generation, (ii) aerosol transport, and (iii) aerosol deposition on a moving free surface in the AJP process. The complex geometries of the deposition head as well as the pneumatic atomizer were modeled in the ansys-fluent environment, based on patented designs in addition to accurate measurements, obtained from 3D X-ray micro-computed tomography (μ-CT) imaging. The entire volume of the constructed geometries was subsequently meshed using a mixture of smooth and soft quadrilateral elements, with consideration of layers of inflation to obtain an accurate solution near the walls. A combined approach, based on the density-based and pressure-based Navier–Stokes formation, was adopted to obtain steady-state solutions and to bring the conservation imbalances below a specified linearization tolerance (i.e., 10−6). Turbulence was modeled using the realizable k-ε viscous model with scalable wall functions. A coupled two-phase flow model was, in addition, set up to track a large number of injected particles. The boundary conditions of the CFD model were defined based on experimental sensor data, recorded from the AJP control system. The accuracy of the model was validated using a factorial experiment, composed of AJ-deposition of a silver nanoparticle ink on a polyimide substrate. The outcomes of this study pave the way for the implementation of physics-driven in situ monitoring and control of AJP. 
    more » « less
  3. We present two prescriptions for broadband (), millimeter-wave antireflection coatings for cryogenic, sintered polycrystalline aluminum oxide optics: one for large-format (700 mm diameter) planar and plano–convex elements, the other for densely packed arrays of quasi-optical elements—in our case, 5 mm diameter half-spheres (called “lenslets”). The coatings comprise three layers of commercially available, polytetrafluoroethylene-based, dielectric sheet material. The lenslet coating is molded to fit the 150 mm diameter arrays directly, while the large-diameter lenses are coated using a tiled approach. We review the fabrication processes for both prescriptions, then discuss laboratory measurements of their transmittance and reflectance. In addition, we present the inferred refractive indices and loss tangents for the coating materials and the aluminum oxide substrate. We find that at 150 GHz and 300 K the large-format coating sample achievestransmittance, and the lenslet coating sample achievestransmittance.

    more » « less
  4. null (Ed.)
    The dielectric waveguide (WG) is an important building block of high-speed and high-bandwidth optical and opto-electronic interconnect networks that operate in the THz frequency regime. At the interface of Si/SiO 2 dielectric waveguides with width above w = 2.5 μm and anisotropic surface roughness, transverse electric (TE) mode surface wave propagation can experience a loss of approximately a = 2 dB/cm; however, propagation losses increase rapidly to near a = 44 dB/cm as the width decreases to w = 500 nm, due to increased interaction of surface waves and sidewall surface roughness that exhibits random distribution inherent to the manufacturing process. Previous works have developed analytic expressions for computing propagation loss in a single dielectric waveguide exhibiting random roughness. More recent works report a = 0.4 dB/cm noting the non-trivial estimation errors in previous theoretical formulations which relied on planar approximations, and highlight the discrepancy in planar approximations vs. the 3-D Volume Current Method. A challenge that remains in the path of designing nanoscale optical interconnects is the dearth of efficient 3-D stochastic computational electromagnetic (CEM) models for multiple tightly coupled optical dielectric waveguides that characterize propagation loss due to random surface roughness in waveguide sidewalls. Through a series of theoretical and numerical experiments developed in the method of finite-difference time-domain (FDTD), we aim to develop stochastic CEM models to quantify propagation loss and facilitate signal & power integrity modeling & simulation of arbitrary configurations of multiple tightly-coupled waveguides, and to gain further insights into loss mechanisms due to random surface roughness in optical interconnects. 
    more » « less
  5. In Part I [Appl. Opt.58,6067(2019)APOPAI003-693510.1364/AO.58.006067], we used a coupled optoelectronic model to optimize a thin-filmCuIn1−<#comment/>ξ<#comment/>Gaξ<#comment/>Se2(CIGS) solar cell with a graded-bandgap photon-absorbing layer and a periodically corrugated backreflector. The increase in efficiency due to the periodic corrugation was found to be tiny and that, too, only for very thin CIGS layers. Also, it was predicted that linear bandgap-grading enhances the efficiency of the CIGS solar cells. However, a significant improvement in solar cell efficiency was found using a nonlinearly (sinusoidally) graded-bandgap CIGS photon-absorbing layer. The optoelectronic model comprised two submodels: optical and electrical. The electrical submodel applied the hybridizable discontinuous Galerkin (HDG) scheme directly to equations for the drift and diffusion of charge carriers. As our HDG scheme sometimes fails due to negative carrier densities arising during the solution process, we devised a new, to the best of our knowledge, computational scheme using the finite-difference method, which also reduces the overall computational cost of optimization. An unfortunate normalization error in the electrical submodel in Part I came to light. This normalization error did not change the overall conclusions reported in Part I; however, some specifics did change. The new algorithm for the electrical submodel is reported here along with updated numerical results. We re-optimized the solar cells containing a CIGS photon-absorbing layer with either (i) a homogeneous bandgap, (ii) a linearly graded bandgap, or (iii) a nonlinearly graded bandgap. Considering the meager increase in efficiency with the periodic corrugation and additional complexity in the fabrication process, we opted for a flat backreflector. The new algorithm is significantly faster than the previous algorithm. Our new results confirm efficiency enhancement of 84% (resp. 63%) when the thickness of the CIGS layer is 600 nm (resp. 2200 nm), similarly to Part I. A hundredfold concentration of sunlight can increase the efficiency by an additional 27%. Finally, the currently used 110-nm-thick layer ofMgF2performs almost as well as optimal single- and double-layer antireflection coatings.

    more » « less