Development of Novel Time-Domain Electromagnetic Methods for Offshore Groundwater Studies: A Data Application from Bat Yam, Israel

Recent marine long-offset transient electromagnetic (LOTEM) measurements yielded the offshore delineation of a fresh groundwater body beneath the seafloor in the region of Bat Yam, Israel. The LOTEM application was effective in detecting this freshwater body underneath the Mediterranean Sea and allowed an estimation of its seaward extent. However, the measured data set was insufficient to understand the hydrogeological configuration and mechanism controlling the occurrence of this fresh groundwater discovery. Especially the lateral geometry of the freshwater boundary, important for the hydrogeological modelling, could not be resolved. Without such an understanding, a rational management of this unexploited groundwater reservoir is not possible. Two new high-resolution marine time-domain electromagnetic methods are theoretically developed to derive the hydrogeological structure of the western aquifer boundary. The first is called Circular Electric Dipole (CED). It is the land-based analogous of the Vertical Electric Dipole (VED), which is commonly applied to detect resistive structures in the subsurface. Although the CED shows exceptional detectability characteristics in the step-off signal towards the sub-seafloor freshwater body, an actual application was not carried out in the extent of this study. It was found that the method suffers from an insufficient signal strength to adequately delineate the resistive aquifer under realistic noise conditions. Moreover, modelling studies demonstrated that severe signal distortions are caused by the slightest geometrical inaccuracies. As a result, a successful application of CED in Israel proved to be rather doubtful. A second method called Differential Electric Dipole (DED) is developed as an alternative to the intended CED method. Compared to the conventional marine time-domain electromagnetic system that commonly applies a horizontal electric dipole transmitter, the DED is composed of two horizontal electric dipoles in an in-line configuration that share a common central electrode. Theoretically, DED has similar detectability/resolution characteristics compared to the conventional LOTEM system. However, the superior lateral resolution towards multi-dimensional resistivity structures make an application desirable. Furthermore, the method is less susceptible towards geometrical errors making an application in Israel feasible. In the extent of this thesis, the novel marine DED method is substantiated using several one-dimensional (1D) and multi-dimensional (2D/3D) modelling studies. The main emphasis lies on the application in Israel. Preliminary resistivity models are derived from the previous marine LOTEM measurement and tested for a DED application. The DED method is effective in locating the two-dimensional resistivity structure at the western aquifer boundary. Moreover, a prediction regarding the hydrogeological boundary conditions are feasible, provided a brackish water zone exists at the head of the interface. A seafloor-based DED transmitter/receiver system is designed and built at the Institute of Geophysics and Meteorology at the University of Cologne. The first DED measurements were carried out in Israel in April 2016. The acquired data set is the first of its kind. The measured data is processed and subsequently interpreted using 1D inversion. The intended aim of interpreting both step-on and step-off signals failed, due to the insufficient data quality of the latter. Yet, the 1D inversion models of the DED step-on signals clearly detect the freshwater body for receivers located close to the Israeli coast. Additionally, a lateral resistivity contrast is observable in the 1D inversion models that allow to constrain the seaward extent of this freshwater body. A large-scale 2D modelling study followed the 1D interpretation. In total, 425 600 forward calculations are conducted to find a sub-seafloor resistivity distribution that adequately explains the measured data. The results indicate that the western aquifer boundary is located at 3600 m - 3700 m before the coast. Moreover, a brackish water zone of 3 Omega*m to 5 Omega*m with a lateral extent of less than 300 m is likely located at the head of the freshwater aquifer. Based on these results, it is predicted that the sub-seafloor freshwater body is indeed open to the sea and may be vulnerable to seawater intrusion.

Numerical and experimental investigation of wedge tip radius effect on wedge plasmons

We report numerical analysis and experimental observation of strongly localized plasmons guided by triangular metal wedges and pay special attention to the effect of smooth (nonzero radius) tips. Dispersion, dissipation, and field structure of such wedge plasmons are analyzed using the compact two-dimensional finite-difference time-domain algorithm. Experimental observation is conducted by the end-fire excitation and near-field scanning optical microscope detection of the predicted plasmons on 40°silver nanowedges with the wedge tip radii of 20, 85, and 125 nm that were fabricated by the focused-ion beam method. The effect of smoothing wedge tips is shown to be similar to that of increasing wedge angle. Increasing wedge angle or wedge tip radius results in increasing propagation distance at the same time as decreasing field localization (decreasing wave number). Quantitative differences between the theoretical and experimental propagation distances are suggested to be due to a contribution of scattered bulk and surface waves near the excitation region as well as the addition of losses due to surface roughness. The theoretical and measured propagation distances are several plasmon wavelengths and are useful for a range of nano-optical applications

Numerical modelling of industrial microwave heating

The numerical modelling of electromagnetic waves has been the focus of many research areas in the past. Some specific applications of electromagnetic wave scattering are in the fields of Microwave Heating and Radar Communication Systems. The equations that govern the fundamental behaviour of electromagnetic wave propagation in waveguides and cavities are the Maxwell's equations. In the literature, a number of methods have been employed to solve these equations. Of these methods, the classical Finite-Difference Time-Domain scheme, which uses a staggered time and space discretisation, is the most well known and widely used. However, it is complicated to implement this method on an irregular computational domain using an unstructured mesh. In this work, a coupled method is introduced for the solution of Maxwell's equations. It is proposed that the free-space component of the solution is computed in the time domain, whilst the load is resolved using the frequency dependent electric field Helmholtz equation. This methodology results in a timefrequency domain hybrid scheme. For the Helmholtz equation, boundary conditions are generated from the time dependent free-space solutions. The boundary information is mapped into the frequency domain using the Discrete Fourier Transform. The solution for the electric field components is obtained by solving a sparse-complex system of linear equations. The hybrid method has been tested for both waveguide and cavity configurations. Numerical tests performed on waveguides and cavities for inhomogeneous lossy materials highlight the accuracy and computational efficiency of the newly proposed hybrid computational electromagnetic strategy.

A compact interferometric sensor design using three waveguide coupling

The use of metal stripes for the guiding of plasmons is a well established technique for the infrared regime and has resulted in the development of a myriad of passive optical components and sensing devices. However, the plasmons suffer from large losses around sharp bends, making the compact design of nanoscale sensors and circuits problematic. A compact alternative would be to use evanescent coupling between two sufficiently close stripes, and thus we propose a compact interferometer design using evanescent coupling. The sensitivity of the design is compared with that achieved using a hand-held sensor based on the Kretschmann style surface plasmon resonance technique. Modeling of the new interferometric sensor is performed for various structural parameters using finite-difference time-domain and COMSOL Multiphysics. The physical mechanisms behind the coupling and propagation of plasmons in this structure are explained in terms of the allowed modes in each section of the device.

Optical waveguide and cavity effects on whispering-gallery mode resonances in a ZnO nanonail

Spatially resolved cathodoluminescence (CL) study of a ZnO nanonail, having thin shank, tapered neck, and hexagonal head sections, is reported. Monochromatic imaging and line scan profiling indicate that the wave guiding and leaking from growth imperfections in addition to the oxygen deficiency variation determine the spatial contrast of CL emissions. Occurrence of resonance peaks at identical wavelengths regardless of CL-excitation spots is inconsistent with the whispering-gallery mode (WGM) resonances of a two-dimensional cavity in the finite difference time domain simulation. However, three dimensioanl cavity simulation produced WGM peaks that are consistent with the experimental spectra, including transverse-electric resonances that are comparable to transverse-magnetic ones.

A modified dual-band microstrip monopole antenna

A modified microstrip-fed planar monopole antenna with open circuited coupled line is presented in this paper. The operational bandwidth of the proposed antenna covers the 2.4 GHz ISM band (2.42-2.48 GHz) and the 5 GHz WLAN band (5 GHz to 6 GHz). The radiating elements occupy a small area of 23×8 mm2. The Finite Difference Time Domain method is used to predict the input impedance of the antenna. The calculated return loss shows very good agreement with measured data. Reasonable antenna gain is observed across the operating band. The measured radiation patterns are similar to those of a simple monopole antenna.

Self assembly of plasmonic core-satellite nano-assemblies mediated by hyperbranched polymer linkers

The morphology of plasmonic nano-assemblies has a direct influence on optical properties, such as localised surface plasmon resonance (LSPR) and surface enhanced Raman scattering (SERS) intensity. Assemblies with core-satellite morphologies are of particular interest, because this morphology has a high density of hot-spots, while constraining the overall size. Herein, a simple method is reported for the self-assembly of gold NPs nano-assemblies with a core-satellite morphology, which was mediated by hyperbranched polymer (HBP) linkers. The HBP linkers have repeat units that do not interact strongly with gold NPs, but have multiple end-groups that specifically interact with the gold NPs and act as anchoring points resulting in nano-assemblies with a large (~48 nm) core surrounded by smaller (~15 nm) satellites. It was possible to control the number of satellites in an assembly which allowed optical parameters such as SPR maxima and the SERS intensity to be tuned. These results were found to be consistent with finite-difference time domain (FDTD) simulations. Furthermore, the multiplexing of the nano-assemblies with a series of Raman tag molecules was demonstrated, without an observable signal arising from the HBP linker after tagging. Such plasmonic nano-assemblies could potentially serve as efficient SERS based diagnostics or biomedical imaging agents in nanomedicine.

Plasmonic effect of annealed gold islands for improving efficiency of organic solar cells

Embedding metallic nanoparticles in organic solar cells can enhance the photoabsorption through light trapping processes. This paper investigates how gold islands obtained by annealing 1–5 nm thick Au layers affect the photoabsorption. Using finite-difference time-domain simulations, the cell efficiency for various island geometries and thicknesses are analyzed and the properties of the islands for maximal photocurrent are discussed. It is shown that a careful choice of size and concentration of gold islands could contribute to enhance the power conversion efficiencies when compared to standard organic solar cell devices. The conclusions are then compared to experimental data for thermally annealed gold islands in bulk heterojunction solar cells. The results of this paper will contribute to the optimization of plasmonic organic solar cell systems and will pave the way for the development of highly efficient organic solar cell devices.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.

Lightning-Induced Current and Voltage on a Rocket in the Presence of Its Trailing Exhaust Plume

This paper presents time-domain characteristics of induced current and voltage on a rocket in the presence of its exhaust plume when an electromagnetic (EM) wave generated by a nearby lightning discharge is incident on it. For the EM-field interaction with the rocket, the finite-difference time-domain technique has been used. The distributed electrical parameters, such as capacitance and inductance of the rocket and its exhaust plume, are computed using the method of moments technique. For the electrical characterization of the exhaust plume, the computational fluid dynamics technique has been used. The computed peak value of the electrical conductivity of the exhaust plume is 0.12 S/m near the exit plane and it reduces to 0.02 S/m at the downstream end. The relative permittivity varies from 0.91 to 0.99. The exhaust plume behaves as a good conductor for EM fields with frequencies less than 2.285 GHz. It has been observed that the peak value of the induced current on the rocket gets enhanced significantly in the presence of the conducting exhaust plume for the rocket and exhaust plume dimensions and parameters studied. The magnitude of the time-varying induced current at the tail is much more than that of any other section of the rocket.

Lightning Induced Voltages on a Rocket and its Exhaust Plume

A preliminary attempt has been made to study the time domain characteristics of the induced voltage and current on the rocket and its exhaust plume (ionized trail) when it is coupled with the transient electromagnetic field generated by a nearby lightning discharge. For the computation, finite difference time domain (FDTD) technique has been used where the object is assumed to be a finite vertical nonuniform transmission line above a perfectly conducting ground. It is seen that the amplitude of the first peak of the induced voltage and current at the mid point of the object is 23.5 kV and 4.9 kA respectively.

Tailoring local density of optical states to control emission intensity and anisotropy of quantum dots in hybrid photonic-plasmonic templates

We report results of controlled tuning of the local density of states (LDOS) in versatile, flexible, and hierarchical self assembled plasmonic templates. Using 5 nm diameter gold (Au) spherical nanoantenna within a polymer template randomly dispersed with quantum dots, we show how the photoluminescence intensity and lifetime anisotropy of these dots can be significantly enhanced through LDOS tuning. Finite difference time domain simulations corroborate the experimental observations and extend the regime of enhancement to a wider range of geometric and spectral parameters bringing out the versatility of these functional plasmonic templates. It is also demonstrated how the templates act as plasmonic resonators for effectively engineer giant enhancement of the scattering efficiency of these nano antenna embedded in the templates. Our work provides an alternative method to achieve spontaneous emission intensity and anisotropy enhancement with true nanoscale plasmon resonators. (C) 2015 AIP Publishing LLC.

On the superposition principle in interference experiments

The superposition principle is usually incorrectly applied in interference experiments. This has recently been investigated through numerics based on Finite Difference Time Domain (FDTD) methods as well as the Feynman path integral formalism. In the current work, we have derived an analytic formula for the Sorkin parameter which can be used to determine the deviation from the application of the principle. We have found excellent agreement between the analytic distribution and those that have been earlier estimated by numerical integration as well as resource intensive FDTD simulations. The analytic handle would be useful for comparing theory with future experiments. It is applicable both to physics based on classical wave equations as well as the non-relativistic Schrodinger equation.