Finite-difference Time-domain Research Articles

Colorimetric Fabry-Pérot Sensor with Hetero-Structured Dielectric for Humidity Monitoring.

A full-color colorimetric humidity sensor with high brightness is proposed by using a hetero-structured dielectric film in a metal-insulator-metal (MIM) resonator. A humidity-responsive polymer is designed to graft on top of a metal-organic frameworks (MOFs) thin film (MOFs-Polymer) as insulator layer in the resonator. Programmable tuning of reflected color is achieved by controlling the polymer thicknesses, and finite difference time domain simulation of light-matter interactions at subwavelength scales proves the dependence of the reflected wavelength on dielectric layer thickness of the resonator. Vivid full-color changing is realized during tracking humidity process due to swelling of the stimuli-responsive polymer. Ultrafast response (≈0.75 s) is achieved for tracking trace H2O from H2O/methanol mixture, which is ≈104 faster than that of the pure polymer-based MIM resonator. Meanwhile, the study observes significant spectral redshift because the porous MOFs film facilitates the preconcentration of external stimulus and improves the detection sensitivity of the resonator. Further, double-channel anti-counterfeiting multiplexing imaging is devised on the MIM resonator by photomask technology. Patterned encoding for security label is achieved on the MIM resonator by engineering humidity-tunable pixels of Au/MOFs-Polymer/Au and humidity-invalid pixels of Au/MOFs/Au.

Research on the random generation model of 2D irregular aggregates in concrete based on fast grid region division method and its ground penetrating radar characteristics

Concrete can be regarded as a composite material with aggregates embedded in hardened slurries. The size, shape and distribution of aggregates directly affect multiple properties of concrete, including mechanical properties, ultrasonic field characteristics of concrete and ground penetrating radar (GPR) wave field characteristics. Therefore, conducting in-depth research on the random generation model of two-dimensional irregular aggregates is particularly important. This article introduces a fast grid region division method (FGRDM), which is an efficient reconstruction algorithm designed for two-dimensional concrete aggregate models. The algorithm takes into consideration the random micro-aggregate structure of concrete. A two-dimensional random distribution model of irregular aggregates in concrete was successfully established by generating concrete models with different particle size distributions and aggregate contents. Furthermore, the finite difference time domain method was used to analyze the GPR wave field of the random model. The significant influence of aggregate shape and distribution on the GPR wave field was discussed. Compared with existing two-dimensional irregular aggregate random generation models, the FGRDM proposed in this paper avoids computational redundancy and generates irregular aggregates at an increasingly fast speed. In the analysis of GPR data, electromagnetic waves are greatly affected by the shape and distribution of aggregates, highlighting the critical role of two-dimensional random distribution models in GPR data analysis. The paper uses FGRDM to generate random irregular aggregates to simulate the distribution of aggregates in actual concrete. Different diameter glass fiber reinforced polymer (GFRP) bars are inserted into the actual concrete and the numerical model to investigate the effect of aggregate size on the recognition of different diameter GFRP bars. The simulation results are generally consistent with the experimental results.

Electromagnetic Response of a Uniformly Moving Resistive Sheet

ABSTRACTBy using the finite‐difference time‐domain (FDTD) method and a two‐port network model, this letter investigates the response of a uniformly moving resistive sheet to illuminating electromagnetic plane waves, at normal incidence. The numerical Doppler frequency shifts and scattering matrix are matched with curve‐fitting expressions. The effects of motion on the impedance matrix, the chain matrix, and the continuity equations are analyzed. Moreover, the total energy of the transferred and reflected waves are studied for the cases of a resistive sheet moving toward or away from the plane wave source.

Enhancing Passive Radiative Cooling Films with Hollow Yttrium-Oxide Spheres Insights from FDTD Simulation.

Passive daytime radiative cooling (PDRC) presents a promising avenue for efficient thermal management without relying on electrical power. In this study, the potential of integrating Hollow Yttrium-Oxide Spheres (HYSs) within a Polydimethylsiloxane (PDMS) matrix to enhance PDRC is investigated. Through a combination of experimental characterization and computational analysis, the optical properties and radiative cooling performance of PDMS films embedded with HYSs are evaluated. These results demonstrate that HYSs significantly improve both solar reflectivity and long-wave infrared (LWIR) emissivity of the PDMS matrix. Finite-Difference Time-Domain (FDTD) simulations confirm the scattering efficiency of HYSs across various wavelength ranges, highlighting their effectiveness as additives for enhancing the radiative properties of passive cooling materials. Experimental validation reveals enhanced reflectivity and emissivity of PDMS films with embedded HYSs, resulting in superior cooling performance compared to non-HYS counterparts. Overall, this study underscores the potential of HYS-infused PDMS films as a promising solution for passive radiative cooling, with broad applicability in diverse domains requiring efficient thermal management solutions. Additionally, these research insights pave the way for establishing an AI database for passive radiative cooling research, offering new avenues for further exploration and application in this field.

Finite-difference time-domain simulation of THz light pulse with orbital angular momentum produced using a spiral phase plate

A silicon spiral phase plate (SPP) was proposed for a pulsed terahertz signal to produce orbital angular momentum with a topological charge of $l=1$. The resulting field distribution after the signal passes through the SPP was simulated using finite-difference time-domain method and observed at different propagation distances. It shows that the helical spatial structure and optical vortex were best observed at the farthest distance of 7.5 cm. Additionally, the cross-section field distribution at different times as the signal passes through the SPP indicates a helical wavefront structure. Moreover, the output of a source with different central frequencies was observed when using the same SPP and showed a hollow distribution indicating a vortex despite not being located in the center.

On the suitability of rigorous coupled-wave analysis for fast optical force simulations

Abstract Optical force responses underpin nanophotonic actuator design, which requires a large number of force simulations to optimize structures. Commonly used computation methods, such as the finite-difference time-domain (FDTD) method, and finite element methods (FEM), are resource intensive and require large amounts of calculation time when multiple structures need to be
compared during optimization. This research demonstrates that performing optical force calculations on 2D-periodic structures using the rigorous coupled-wave analysis (RCWA) method is typically on the order of 10 times faster than other approaches and with sufficient accuracy to suit optical design purposes. Moreover, this speed increase is available on consumer grade laptops, avoiding the need for a high performance computing resource.

Three-Dimensional Reconstruction of Retaining Structure Defects from Crosshole Ground Penetrating Radar Data Using a Generative Adversarial Network

Crosshole ground penetrating radar (GPR) is an efficient method for ensuring the quality of retaining structures without the need for excavation. However, interpreting crosshole GPR data is time-consuming and prone to inaccuracies. To address this challenge, we proposed a novel three-dimensional (3D) reconstruction method based on a generative adversarial network (GAN) to recover 3D permittivity distributions from crosshole GPR images. The established framework, named CGPR2VOX, integrates a fully connected layer, a residual network, and a specialized 3D decoder in the generator to effectively translate crosshole GPR data into 3D permittivity voxels. The discriminator was designed to enhance the generator’s performance by ensuring the physical plausibility and accuracy of the reconstructed models. This adversarial training mechanism enables the network to learn non-linear relationships between crosshole GPR data and subsurface permittivity distributions. CGPR2VOX was trained using a dataset generated through finite-difference time-domain (FDTD) simulations, achieving precision, recall and F1-score of 91.43%, 96.97% and 94.12%, respectively. Model experiments validate that the relative errors of the estimated positions of the defects were 1.67%, 1.65%, and 1.30% in the X-, Y-, and Z-direction, respectively. Meanwhile, the method exhibits noteworthy generalization capabilities under complex conditions, including condition variations, heterogeneous materials and electromagnetic noise, highlighting its reliability and effectiveness for practical quality assurance of retaining structures.

Inverse design of polaritonic devices

Polaritons, arising from the strong coupling between excitons and photons within microcavities, hold promise for optoelectronic and all-optical devices. They have found applications in various domains, including low-threshold lasers and quantum information processing. To realize complex functionalities, non-intuitive designs for polaritonic devices are required. In this contribution, we use finite-difference time-domain simulations of the dissipative Gross–Pitaevskii equation, written in a differentiable manner, and combine it with an adjoint formulation. Such a method allows us to use topology optimization to engineer the potential landscape experienced by polariton condensates to tailor its characteristics on demand. The potential directly translates to a blueprint for a functional device, and various fabrication and optical control techniques can experimentally realize it. We inverse-design a selection of polaritonic devices, i.e., a structure that spatially shapes the polaritons into a flat-top distribution, a metalens that focuses a polariton, and a nonlinearly activated isolator. The functionalities are preserved when employing realistic fabrication constraints such as minimum feature size and discretization of the potential. Our results demonstrate the utility of inverse design techniques for polaritonic devices, providing a stepping stone toward future research in optimizing systems with complex light–matter interactions.

Xenon plasma-focused ion beam milling for fabrication of high-purity, bright single-photon sources operating in the C-band

Electron beam lithography is a standard method for fabricating photonic micro and nanostructures around semiconductor quantum dots (QDs), which are crucial for efficient single and indistinguishable photon sources in quantum information processing. However, this technique is difficult for direct 3D control of the structure shape, complicating the design and enlarging the 2D footprint to suppress in-plane photon leakage while directing photons into the collecting lens aperture. Here, we present an alternative approach to employ xenon plasma-focused ion beam (Xe-PFIB) technology as a reliable method for the 3D shaping of photonic structures containing low-density self-assembled InAs/InP quantum dots emitting in the C-band range of the 3rd telecommunication window. The method is optimized to minimize the possible ion-beam-induced material degradation, which allows exploration of both non-deterministic and deterministic fabrication approaches, resulting in photonic structures naturally shaped as truncated cones. As a demonstration, we fabricate mesas using a heterogeneously integrated structure with a QD membrane atop an aluminum mirror and silicon substrate. Finite-difference time-domain simulations show that the angled sidewalls significantly increase the emission collection efficiency to approx. 0.9 for NA = 0.65. We demonstrate experimentally a high purity of pulsed single-photon emission (∼99%) and a superior extraction efficiency value reported in the C-band of η = 24 ± 4%.

Evaluation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using a high figure-of-merit plasmonic multimode refractive index optical sensor

In recent years, following the outbreak of the COVID-19 pandemic, there has been a significant increase in cases of SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) and related deaths worldwide. Despite the pandemic nearing its end due to the introduction of mass-produced vaccines against SARS-CoV-2, early detection and diagnosis of the virus remain crucial in preventing disease progression. This article explores the rapid identification of SARS-CoV-2 by implementing a multimode plasmonic refractive index (MMRI) optical sensor, developed based on the split ring resonator (SRR) design. The Finite Difference Time Domain (FDTD) numerical solution method simulates the sensor. The studied sensor demonstrates three resonance modes within the reflection spectrum ranging from 800 nm to 1400 nm. Its material composition and dimensional parameters are optimized to enhance the sensor’s performance. The research indicates that all three resonance modes exhibit strong performance with high sensitivity and figures of merit. Notably, the first mode achieves an exceptional sensitivity of 557 nm/RIU, while the third mode exhibits a commendable sensitivity of 453 nm/RIU and a Figure of Merit (FOM) of 45 RIU−1. These findings suggest that the developed MMRI optical sensor holds significant potential for the early and accurate detection of SARS-CoV-2, contributing to improved disease management and control efforts.

Advances in Design and Fabrication of Micro-Structured Solid Targets for High-Power Laser-Matter Interaction

Accelerated particles have multiple applications in materials research, medicine, and the space industry. In contrast to classical particle accelerators, laser-driven acceleration at intensities greater than 1018 W/cm2, currently achieved at TW and PW laser facilities, allow for much larger electric field gradients at the laser focus point, several orders of magnitude higher than those found in conventional kilometer-sized accelerators. It has been demonstrated that target design becomes an important factor to consider in ultra-intense laser experiments. The energetic and spatial distribution of the accelerated particles strongly depends on the target configuration. Therefore, target engineering is one of the key approaches to optimizing energy transfer from the laser to the accelerated particles. This paper provides an overview of recent progress in 2D and 3D micro-structured solid targets, with an emphasis on fabrication procedures based on laser material processing. Recently, 3D laser lithography, which involves Two-Photon Absorption (TPA) effects in photopolymers, has been proposed as a technique for the high-resolution fabrication of 3D micro-structured targets. Additionally, laser surface nano-patterning followed by the replication of the patterns through molding, has been proposed and could become a cost-effective and reliable solution for intense laser experiments at high repetition rates. Recent works on numerical simulations have also been presented. Using particle-in-cell (PIC) simulation software, the importance of structured micro-target design in the energy absorption process of intense laser pulses—producing localized extreme temperatures and pressures—was demonstrated. Besides PIC simulations, the Finite-Difference Time-Domain (FDTD) numerical method offers the possibility to generate the specific data necessary for defining solid target material properties and designing their optical geometries with high accuracy. The prospects for the design and technological fabrication of 3D targets for ultra-intense laser facilities are also highlighted.

Terahertz Non-Destructive Testing of Porosity in Multi-Layer Thermal Barrier Coatings Based on Small-Sample Data

Accurately characterizing the internal porosity rate of thermal barrier coatings (TBCs) was essential for prolonging their service life. This work concentrated on atmospheric plasma spray (APS)-prepared TBCs and proposed the utilization of terahertz non-destructive detection technology to evaluate their internal porosity rate. The internal porosity rates were ascertained through a metallographic analysis and scanning electron microscopy (SEM), followed by the reconstruction of the TBC model using a four-parameter method. Terahertz time-domain simulation data corresponding to various porosity rates were generated employing the time-domain finite difference method. In simulating actual test signals, white noise with a signal-to-noise ratio of 10 dB was introduced, and various wavelet transforms were utilized for denoising purposes. The effectiveness of different signal processing techniques in mitigating noise was compared to extract key features associated with porosity. To address dimensionality challenges and further enhance model performance, kernel principal component analysis (kPCA) was employed for data processing. To tackle issues related to limited sample sizes, this work proposed to use the Siamese neural network (SNN) and generative adversarial network (GAN) algorithms to solve this challenge in order to improve the generalization ability and detection accuracy of the model. The efficacy of the constructed model was assessed using multiple evaluation metrics; the results indicate that the novel hybrid WT-kPCA-GAN model achieves a prediction accuracy exceeding 0.9 while demonstrating lower error rates and superior predictive performance overall. Ultimately, this work presented an innovative, convenient, non-destructive online approach that was safe and highly precise for measuring the porosity rate of TBCs, particularly in scenarios involving small sample sizes facilitating assessments regarding their service life.

Compact THz Polarization Splitter with Large Bandwidth Based on Cascade Hexagonal Porous Two-Core PCF

ABSTRACT A novel cyclic olefin copolymer (COC)-based terahertz (THz) polarization beam splitter (PBS) employing photonic crystal fiber (PCF) with cascade hexagonal porous two-core structure is proposed. The porous two-core is achieved by replacing the original two large air holes with three layers of small air holes. The mode coupling characteristics between the even and odd modes in the x- and y-polarization directions for the proposed THz PBS are explored using the time-domain finite difference (FDTD) method. The influence of PBS structural parameters on its coupling length is investigated. Numerical analysis results indicate that the bandwidths of 227.8 GHz (ranging from 1.8970 to 2.1248 THz) and 147 GHz (ranging from 1.928 to 2.075 THz) are acquired for the x- and y-polarization fundamental modes, respectively. The highest extinction ratios (ERs) are 114.5 dB and 72.3 dB, and the minimum insertion losses (ILs) are 0.0142 dB and 0.0369 dB for the x- and y-polarization fundamental modes, respectively. Moreover, the length of the proposed PBS is merely 11.655 mm. The proposed PBS will find potential applications in polarization-diverse THz systems, including broadband line-of-sight communication, real-time security checks, and high-resolution imaging.

Free Electron Density Gradients Enhanced Biosensor for Ultrasensitive and Accurate Affinity Assessment of the Immunotherapy Drugs.

Accurate affinity assessments play an important role in drug discovery, screening, and efficacy evaluation. Label-free affinity biosensors are recognized as a dependable and standard technology for addressing this challenge. This study constructs a free electron density gradient-enhanced meta-surface plasmon resonance (FED-MSPR) biosensor through a finite-difference time-domain simulation model, the biosensor demonstrates superior detection performance in accurately determining affinity and kinetic rate constants. By controlling the dielectric properties of the metal on the surface of the nanocup arrays, the plasmon resonance effects are easily tuned without changing the nanostructure design. Compared with the single-layer gold chip, the triple-layer FED-MSPR chip demonstrated a four-fold improvement in resolution at the optimal resonance peak. Additionally, the sensitivity and figure of merit (FOM) of the multi-layer chip increased by 3.5 and 7.99 times, respectively. Following modification with high- and low-staggered carboxylation, the noise-signal ratio and baseline stability of the real-time kinetic curves based on these chips are significantly enhanced. The developed carboxylation FED-MSPR platform is successfully used to perform affinity assays for Adalimumab and TNF-α protein, resulting in favorable dynamic curves. These findings validate the proposed FED-MSPR biosensor platform as cost-effective, rapid, sensitive, and label-free, facilitating real-time quality control in drug development.

Application of the FDTD Method for Multivariate Analysis of the Influence of Conductivity and the Arrangement of Hollows Inside Bricks on the Values of Electric Field Intensity

The article contains a numerical analysis of the effects of electromagnetic wave propagation in an area containing a non-ideal, non-uniform, and absorbing dielectric. The analysis concerned the influence of the structure of the building material and its electrical parameters on the electric field intensity. The analysis took into account the variability of the number of hollows in the brick, the width of hollows, as well as the arrangement of these hollows relative to each other using the example of two types of bricks. The article also provides the most commonly used values of electrical parameters for building materials (brick, plaster). For this reason, the article includes results for different values of conductivity (0–0.2 S/m). The FDTD (Finite Difference Time Domain) method was used for multivariate analysis. The aim was to verify the correctness of the numerical assumptions adopted. Using the example of the most commonly used wall structure in construction, the results obtained using the FDTD method were compared with values obtained using another numerical method, the finite element method (FEM). The influence of an additional layer of plaster on the considered wall on the electric field was also checked. The analysis showed that a symmetrical arrangement of bricks results in higher values of the electric field by an average of 20%. Of course, this depends on the length of the hollows and the number of holes. The highest field values occur at low conductivities (0–0.04 S/m). A brick wall with a larger number of hollows and a symmetrical brick arrangement shows the highest electric field intensity, especially for hollow sizes (0.009–0.015 m).

Construction of Ag-modified ZnO/g-C3N4 heterostructure for enhanced photocatalysis performance.

ZnO/g-C3N4 heterojunction modified with Ag nanoparticles (ZnO/CN/Ag) was synthesized by depositing ZnO nanorods/Ag nanoparticles onto g-C3N4 nanosheets. Under xenon lamp irradiation, 99% of Rhodamine B (RhB) was degraded by ZnO/CN/Ag-5% composite within 30min, which was much higher than the degradation efficiency of ZnO and ZnO/CN. The synergistic effect of g-C3N4 and ZnO, along with the localized surface plasmon resonance effect of Ag NPs, contributes to the improvement of photocatalytic performance. Ag nanoparticle provides another charge transfer path from g-C3N4 to ZnO, which speeds up the separation of electron-hole pairs. Meanwhile, the catalyst had good stability and recyclability. Finite-difference time-domain method and the density functional theory were used to obtain the charge transfer process. The photodegradation process has been studied in depth.

High-Resolution Distance Dependence Interrogation of Scanning Ion Conductance Microscopic Tip-Enhanced Raman Spectroscopy Enabled by Two-Dimensional Molybdenum Disulfide Substrates.

Scanning ion conductance microscopy (SICM) is a powerful surface imaging tool used in the electrolytic environment. Tip-enhanced Raman spectroscopy (TERS) can give more information in addition to the morphology provided by the SICM by utilizing label-free Raman spectroscopy aided by the localized plasmonic enhancement from the metal-coated probes. In this study, the integration of SICM with TERS is demonstrated through employing a silver-coated plasmonic nanopipette. Leveraging a two-dimensional (2D) molybdenum disulfide (MoS2) as a model system, the SICM-TERS enhancement factor was estimated to be ∼105, supported by finite-difference time-domain (FDTD) simulation. Moreover, the subnanometer distance dependence SICM-TERS study reveals the tensile stress and structural changes caused by the nanopipette. These findings illustrate the potential of SICM-TERS for providing comprehensive morphological and chemical insights into electrolytic environments, paving the way for future investigations of electrocatalytic and biological systems.

Microstructured Reflective Coatings on Commodity Textiles for Passive Personal Cooling.

As the effects of climate change become more severe and widespread, maintaining personal thermal homeostasis becomes necessary for survival. In principle, advanced textiles and garments have the ability to leverage light absorption, transmission and/or reflection, in addition to straightforward convection, to heat or cool bodies in extreme temperature conditions. For cooling, in particular, surfaces adept at selectively reflecting or refracting high-energy incident light (200 nm-2.5 mm) from the sun while transmitting or emitting infrared light (8-13 mm) from radiant body heat boast the ability to maintain cooler body temperatures, even when exposed to direct sunlight and the open sky. Here, we present a strategy to transform common clothing into implements for passive personal cooling. As confirmed by Mie scattering calculations, cheap and biocompatible calcium carbonate and barium sulfate micro/nanoparticles are found to serve as suitable reflectors for radiative cooling. Finite-difference time domain simulations reveal, surprisingly, that higher reflectance is achieved with surface coatings containing these materials, as compared to extruded metamaterial fibers containing CaCO3 and BaSO4 particles embedded within a polymer matrix. A stepwise process involving photoinitiated chemical vapor deposition and ion-exchange driven crystal growth is used to create a lamellar composite coating comprised of alternating CaCO3 and BaSO4 nano/microparticle layers directly on the surface of common fabrics. A polyester poplin fabric coated in this manner shows a cooling ability of up to 8 °C compared to an uncoated sample, achieving a maximum cooling of 6 °C below ambient temperature. Wash and durability testing of the lamellar coating reveal no mechanical degradation and no evident attenuation in the material's performance, affirming its resilience and long-term effectiveness as a functional textile coating for personal cooling. We also assess the performance of our coated fabrics in multiple outdoor environments to conclude that we can achieve up to 3.4 °C of sub-ambient cooling in optically complex built environments.

Multifunctional electrically switchable metalens

Here, we propose three all-solid-state, electrically switchable, transmissive, and multifunctional metalens arrays comprising tunable metal-oxide material BaTiO3 (BTO) nanopillars with different structural parameters. To produce the required phase profile for each metalens, the rectangular nanopillars, as the unit cell of the BTO structure, are developed to support arbitrary combinations of two independent phase shifts (0−2π) with bias voltage states of 0 V and 60 V, at first. Second, the structure parameters of different functional metalens arrays are generated efficiently and accurately based on the dual-phase modulation characteristics of a single structure. Finally, we use the finite-difference time-domain method to simulate the three kinds of switchable metalenses. The results show that the three metalenses can realize the function of adjustable focus position, switchable beam focusing and beam deflection, and switchable beam focusing and beam splitting.

Application of the FDTD Method to Analyze the Influence of Brick Complexity on Electromagnetic Wave Propagation

This article presents a numerical analysis of the effects related to the propagation of electromagnetic waves in an area containing a non-ideal, non-uniform, and absorbing dielectric. The analysis concerns the influence of electrical parameters, the structure of the building material, and the layering of the wall on the values of the electric field intensity. A multivariate analysis was carried out with different conductivity values. Homogeneous materials (e.g., solid brick) can be analyzed using the analytical method. In the case of complex materials containing, e.g., hollows (brick with hollows, hollow block), it is necessary to use the numerical method. The FDTD (finite difference time domain) method was used to assess the dependence of the electric field intensity on the layering, the length of hollows in bricks, and the material loss. In order to check the correctness of the adopted numerical assumptions, a series of tests related to the discretization of the model was carried out. The article also presents the influence of changing the length of hollows in bricks on the values of the electric field intensity at a frequency of 2.4 GHz. The instantaneous field distributions and maximum values of the electric field intensity are presented. In the model with a two-layer wall, regardless of the conductivity, the field values were the same for the two models, where the difference in the percentage of ceramic mass in the brick was 8%. A 12% decrease in the percentage of ceramic mass in the brick resulted in a 15% increase in the value of the area between a single-layer and a double-layer wall made of clinker bricks. At a conductivity of 0.04 S/m for a single-layer wall, the field values were similar for all brick variants.