Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Composition, nanostructure, and optical properties of silver and silver-copper lusters

Lusters are composite thin layers of coinage metal nanoparticles in glass displaying peculiar optical properties and obtained by a process involving ionic exchange, diffusion, and crystallization. In particular, the origin of the high reflectance (golden-shine) shown by those layers has been subject of some discussion. It has been attributed to either the presence of larger particles, thinner multiple layers or higher volume fraction of nanoparticles. The object of this paper is to clarify this for which a set of laboratory designed lusters are analysed by Rutherford backscattering spectroscopy, transmission electron microscopy, x-ray diffraction, and ultraviolet-visible spectroscopy. Model calculations and numerical simulations using the finite difference time domain method were also performed to evaluate the optical properties. Finally, the correlation between synthesis conditions, nanostructure, and optical properties is obtained for these materials.

Development and Analysis of a Compact Dual-Band Coplanar Antenna

The coplanar wave guide is an attractive device in microwave integrated circuits due to its uniplanar nature, ease of fabrication and low production cost. Several attempts are already done to explore the radiating modes in coplanar wave guide transmission lines. Usually coplanar wave guides are excited by an SMA connector with its centre conductor connected to the exact middle of the centre strip and the outer ground conductor to the two ground strips. The mode excited on it is purely a bound mode. The E-field distribution in the two slots are out of phase and there for cancels at the far field. This thesis addresses an attempt to excite an in phase E-field distribution in the two slots of the co planar wave guide by employing a feed asymmetry, in order to get radiation from the two large slot discontinuities of the coplanar waveguide. The omni directional distribution of the radiating energy can be achieved by widening the centre strip.The first part of the thesis deals with the investigations on the resonance phenomena of conventional coplanar waveguides at higher frequency bands. Then an offset fed open circuited coplanar waveguide supporting resonance/radiation phenomena is analyzed. Finally, a novel compact co planar antenna geometry with dual band characteristics, suitable for mobile terminal applications is designed and characterized using the inferences from the above study.

Octagonal Microstrip Patch Antenna For Dual Band Applications

A dual port dual polarized octagonal microstrip patch antenna suitable for dual band applications is discussed theoretically and experimentally. The antenna exhibits good impedance bandwidth, gain and broad radiation patterns. Parameters predicted by the Conformal Finite Difference Time Domain algorithm show good agreement with the simulated results and experimental observations

Análise e síntese de antenas e superfícies seletivas de frequência utilizando computação evolucionária e inteligência de enxames

The frequency selective surfaces, or FSS (Frequency Selective Surfaces), are structures consisting of periodic arrays of conductive elements, called patches, which are usually very thin and they are printed on dielectric layers, or by openings perforated on very thin metallic surfaces, for applications in bands of microwave and millimeter waves. These structures are often used in aircraft, missiles, satellites, radomes, antennae reflector, high gain antennas and microwave ovens, for example. The use of these structures has as main objective filter frequency bands that can be broadcast or rejection, depending on the specificity of the required application. In turn, the modern communication systems such as GSM (Global System for Mobile Communications), RFID (Radio Frequency Identification), Bluetooth, Wi-Fi and WiMAX, whose services are highly demanded by society, have required the development of antennas having, as its main features, and low cost profile, and reduced dimensions and weight. In this context, the microstrip antenna is presented as an excellent choice for communications systems today, because (in addition to meeting the requirements mentioned intrinsically) planar structures are easy to manufacture and integration with other components in microwave circuits. Consequently, the analysis and synthesis of these devices mainly, due to the high possibility of shapes, size and frequency of its elements has been carried out by full-wave models, such as the finite element method, the method of moments and finite difference time domain. However, these methods require an accurate despite great computational effort. In this context, computational intelligence (CI) has been used successfully in the design and optimization of microwave planar structures, as an auxiliary tool and very appropriate, given the complexity of the geometry of the antennas and the FSS considered. The computational intelligence is inspired by natural phenomena such as learning, perception and decision, using techniques such as artificial neural networks, fuzzy logic, fractal geometry and evolutionary computation. This work makes a study of application of computational intelligence using meta-heuristics such as genetic algorithms and swarm intelligence optimization of antennas and frequency selective surfaces. Genetic algorithms are computational search methods based on the theory of natural selection proposed by Darwin and genetics used to solve complex problems, eg, problems where the search space grows with the size of the problem. The particle swarm optimization characteristics including the use of intelligence collectively being applied to optimization problems in many areas of research. The main objective of this work is the use of computational intelligence, the analysis and synthesis of antennas and FSS. We considered the structures of a microstrip planar monopole, ring type, and a cross-dipole FSS. We developed algorithms and optimization results obtained for optimized geometries of antennas and FSS considered. To validate results were designed, constructed and measured several prototypes. The measured results showed excellent agreement with the simulated. Moreover, the results obtained in this study were compared to those simulated using a commercial software has been also observed an excellent agreement. Specifically, the efficiency of techniques used were CI evidenced by simulated and measured, aiming at optimizing the bandwidth of an antenna for wideband operation or UWB (Ultra Wideband), using a genetic algorithm and optimizing the bandwidth, by specifying the length of the air gap between two frequency selective surfaces, using an optimization algorithm particle swarm

Detecção e localização de descargas parciais em transformadores de potência utilizando sinais acústicoas transitórios e o método GSPS

Neste trabalho, é desenvolvido um método de localização de descargas parciais, em transformadores de potência, baseado no algoritmo GPS (Global Positioning System). Para a análise da estrutura, foi desenvolvido um solftware, no qual as equações diferenciais que representam a propagação de ondas acústicas são resolvidas numericamente através do método Acoustic Finite Difference Time Domain (AFDTD), cujo domínio computacional é truncado através da técnica CPML (Convolutional Perfectly Matched Layer). Os resultados obtidos são comparados a estimativas produzidas utilizando-se sinais elétricos relativos às descargas.

Study of the thermal effects induced by magnetic resonance on endocardial leads: numerical models and experimental validation

Magnetic resonance imaging (MRI) is today precluded to patients bearing active implantable medical devices AIMDs). The great advantages related to this diagnostic modality, together with the increasing number of people benefiting from implantable devices, in particular pacemakers(PM)and carioverter/defibrillators (ICD), is prompting the scientific community the study the possibility to extend MRI also to implanted patients. The MRI induced specific absorption rate (SAR) and the consequent heating of biological tissues is one of the major concerns that makes patients bearing metallic structures contraindicated for MRI scans. To date, both in-vivo and in-vitro studies have demonstrated the potentially dangerous temperature increase caused by the radiofrequency (RF) field generated during MRI procedures in the tissues surrounding thin metallic implants. On the other side, the technical evolution of MRI scanners and of AIMDs together with published data on the lack of adverse events have reopened the interest in this field and suggest that, under given conditions, MRI can be safely performed also in implanted patients. With a better understanding of the hazards of performing MRI scans on implanted patients as well as the development of MRI safe devices, we may soon enter an era where the ability of this imaging modality may be more widely used to assist in the appropriate diagnosis of patients with devices. In this study both experimental measures and numerical analysis were performed. Aim of the study is to systematically investigate the effects of the MRI RF filed on implantable devices and to identify the elements that play a major role in the induced heating. Furthermore, we aimed at developing a realistic numerical model able to simulate the interactions between an RF coil for MRI and biological tissues implanted with a PM, and to predict the induced SAR as a function of the particular path of the PM lead. The methods developed and validated during the PhD program led to the design of an experimental framework for the accurate measure of PM lead heating induced by MRI systems. In addition, numerical models based on Finite-Differences Time-Domain (FDTD) simulations were validated to obtain a general tool for investigating the large number of parameters and factors involved in this complex phenomenon. The results obtained demonstrated that the MRI induced heating on metallic implants is a real risk that represents a contraindication in extending MRI scans also to patient bearing a PM, an ICD, or other thin metallic objects. On the other side, both experimental data and numerical results show that, under particular conditions, MRI procedures might be consider reasonably safe also for an implanted patient. The complexity and the large number of variables involved, make difficult to define a unique set of such conditions: when the benefits of a MRI investigation cannot be obtained using other imaging techniques, the possibility to perform the scan should not be immediately excluded, but some considerations are always needed.

Nanoporous alumina based hypersonic phononic hybrid structures studied by brillouin spectroscopy

Phononische Kristalle sind strukturierte Materialien mit sich periodisch ändernden elastischen Moduln auf der Wellenlängenskala. Die Interaktion zwischen Schallwellen und periodischer Struktur erzeugt interessante Interferenzphänomene, und phononische Kristalle erschließen neue Funktionalitäten, die in unstrukturierter Materie unzugänglich sind. Hypersonische phononische Kristalle im Speziellen, die bei GHz Frequenzen arbeiten, haben Periodizitäten in der Größenordnung der Wellenlänge sichtbaren Lichts und zeigen daher die Wege auf, gleichzeitig Licht- und Schallausbreitung und -lokalisation zu kontrollieren, und dadurch die Realisierung neuartiger akusto-optischer Anordnungen. Bisher bekannte hypersonische phononische Kristalle basieren auf thermoplastischen Polymeren oder Epoxiden und haben nur eingeschränkte thermische und mechanische Stabilität und mechanischen Kontrast. Phononische Kristalle, die aus mit Flüssigkeit gefüllten zylindrischen Kanälen in harter Matrix bestehen, zeigen einen sehr hohen elastischen Kontrast und sind bislang noch unerforscht. In dieser Dissertation wird die experimentelle Untersuchung zweidimensionaler hypersonischer phononischer Kristalle mit hexagonaler Anordnung zylindrischer Nanoporen basierend auf der Selbstorganisation anodischen Aluminiumoxids (AAO) beschrieben. Dazu wird die Technik der hochauflösenden inelastischen Brillouin Lichtstreuung (BLS) verwendet. AAO ist ein vielsetiges Modellsystem für die Untersuchung reicher phononischer Phänomene im GHz-Bereich, die eng mit den sich in den Nanoporen befindlichen Flüssigkeiten und deren Interaktion mit der Porenwand verknüpft sind. Gerichteter Fluss elastischer Energie parallel und orthogonal zu der Kanalachse, Lokalisierung von Phononen und Beeinflussung der phononischen Bandstruktur bei gleichzeitig präziser Kontrolle des Volumenbruchs der Kanäle (Porosität) werden erörtert. Außerdem ermöglicht die thermische Stabilität von AAO ein temperaturabhängiges Schalten phononischer Eigenschaften infolge temperaturinduzierter Phasenübergänge in den Nanoporen. In monokristallinen zweidimensionalen phononischen AAO Kristallen unterscheiden sich die Dispersionsrelationen empfindlich entlang zweier hoch symmetrischer Richtungen in der Brillouinzone, abhängig davon, ob die Poren leer oder gefüllt sind. Alle experimentellen Dispersionsrelationen werden unter Zuhilfenahme theoretische Ergebnisse durch finite Elemente Analyse (FDTD) gedeutet. Die Zuordnung der Verschiebungsfelder der elastischen Wellen erklärt die Natur aller phononischen Moden.

New enhancement strategies for plasmon-enhanced fluorescence biosensors

This thesis investigates metallic nanostructures exhibiting surface plasmon resonance for the amplification of fluorescence signal in sandwich immunoassays. In this approach, an analyte is captured by an antibody immobilized on a plasmonic structure and detected by a subsequently bound fluorophore labeled detection antibody. The highly confined field of surface plasmons originates from collective charge oscillations which are associated with high electromagnetic field enhancements at the metal surface and allow for greatly increased fluorescence signal from the attached fluorophores. This feature allows for improving the signal-to-noise ratio in fluorescence measurements and thus advancing the sensitivity of the sensor platform. In particular, the thesis presents two plasmonic nanostructures that amplify fluorescence signal in devices that rely on epifluorescence geometry, in which the fluorophore absorbs and emits light from the same direction perpendicular to the substrate surface.rnThe first is a crossed relief gold grating that supports propagating surface plasmon polaritons (SPPs) and second, gold nanoparticles embedded in refractive index symmetric environment exhibiting collective localized surface plasmons (cLSPs). Finite-difference time-domain simulations are performed in order to design structures for the optimum amplification of established Cy5 and Alexa Fluor 647 fluorophore labels with the absorption and emission wavelengths in the red region of spectrum. The design takes into account combined effect of surface plasmon-enhanced excitation rate, directional surface plasmon-driven emission and modified quantum yield for characteristic distances in immunoassays. Homebuilt optical instruments are developed for the experimental observation of the surface plasmon mode spectrum, measurements of the angular distribution of surface plasmon-coupled fluorescence light and a setup mimicking commercial fluorescence reading systems in epifluorescence geometry.rnCrossed relief grating structures are prepared by interference lithography and multiple copies are made by UV nanoimprint lithography. The fabricated crossed diffraction gratings were utilized for sandwich immunoassay-based detection of the clinically relevant inflammation marker interleukin 6 (IL-6). The enhancement factor of the crossed grating reached EF=100 when compared to a flat gold substrate. This result is comparable to the highest reported enhancements to date, for fluorophores with relatively high intrinsic quantum yield. The measured enhancement factor excellently agrees with the predictions of the simulations and the mechanisms of the enhancement are explained in detail. Main contributions were the high electric field intensity enhancement (30-fold increase) and the directional fluorescence emission at (4-fold increase) compared to a flat gold substrate.rnCollective localized surface plasmons (cLSPs) hold potential for even stronger fluorescence enhancement of EF=1000, due to higher electric field intensity confinement. cLSPs are established by diffractive coupling of the localized surface plasmon resonance (LSPR) of metallic nanoparticles and result in a narrow resonance. Due to the narrow resonance, it is hard to overlap the cLSPs mode with the absorption and emission bands of the used fluorophore, simultaneously. Therefore, a novel two resonance structure that supports SPP and cLSP modes was proposed. It consists of a 2D array of cylindrical gold nanoparticles above a low refractive index polymer and a silver film. A structure that supports the proposed SPP and cLSP modes was prepared by employing laser interference lithography and the measured mode spectrum was compared to simulation results.rn

Optical properties of a nanomatch-like plasmonic structure

The optical properties of a match-like plasmonic nanostructure are numerically investigated using full-wave finite-difference time-domain analysis in conjunction with dispersive material models. This work is mainly motivated by the developed technique enabling reproducible fabrication of nanomatch structures as well as the growing applications that utilize the localized field enhancement in plasmonic nanostructures. Our research revealed that due to the pronounced field enhancement and larger resonance tunabilities, some nanomatch topologies show potentials for various applications in the field of, e.g., sensing as well as a novel scheme for highly reproducible tips in scanning near field optical microscopy, among others. Despite the additional degrees of freedom that are offered by the composite nature of the proposed nanomatch topology, the paper also reflects on a fundamental complication intrinsic to the material interfaces especially in the nanoscale: stoichiometric mixing. We conclude that the specificity in material modeling will become a significant issue in future research on functionalized composite nanostructures.

Detecting Improvised Explosive Devices Via Forward Looking Ground Penetrating Radar

Forward-looking ground penetrating radar shows promise for detection of improvised explosive devices in active war zones. Because of certain insurmountable physical limitations, post-processing algorithm development is the most popular research topic in this field. One such investigative avenue explores the worthiness of frequency analysis during data post-processing. Using the finite difference time domain numerical method, simulations are run to test both mine and clutter frequency response. Mines are found to respond strongest at low frequencies and cause periodic changes in ground penetrating radar frequency results. These results are called into question, however, when clutter, a phenomenon generally known to be random, is also found to cause periodic frequency effects. Possible causes, including simulation inaccuracy, are considered. Although the clutter models used are found to be inadequately random, specular reflections of differing periodicity are found to return from both the mine and the ground. The presence of these specular reflections offers a potential alternative method of determining a mine’s presence.

Study of and design procedure for dual circularly polarized waveguide slot arrays

For the past sixty years, waveguide slot radiator arrays have played a critical role in microwave radar and communication systems. They feature a well-characterized antenna element capable of direct integration into a low-loss feed structure with highly developed and inexpensive manufacturing processes. Waveguide slot radiators comprise some of the highest performance—in terms of side-lobe-level, efficiency, etc. — antenna arrays ever constructed. A wealth of information is available in the open literature regarding design procedures for linearly polarized waveguide slots. By contrast, despite their presence in some of the earliest published reports, little has been presented to date on array designs for circularly polarized (CP) waveguide slots. Moreover, that which has been presented features a classic traveling wave, efficiency-reducing beam tilt. This work proposes a unique CP waveguide slot architecture which mitigates these problems and a thorough design procedure employing widely available, modern computational tools. The proposed array topology features simultaneous dual-CP operation with grating-lobe-free, broadside radiation, high aperture efficiency, and good return loss. A traditional X-Slot CP element is employed with the inclusion of a slow wave structure passive phase shifter to ensure broadside radiation without the need for performance-limiting dielectric loading. It is anticipated this technology will be advantageous for upcoming polarimetric radar and Ka-band SatCom systems. The presented design methodology represents a philosophical shift away from traditional waveguide slot radiator design practices. Rather than providing design curves and/or analytical expressions for equivalent circuit models, simple first-order design rules – generated via parametric studies — are presented with the understanding that device optimization and design will be carried out computationally. A unit-cell, S-parameter based approach provides a sufficient reduction of complexity to permit efficient, accurate device design with attention to realistic, application-specific mechanical tolerances. A transparent, start-to-finish example of the design procedure for a linear sub-array at X-Band is presented. Both unit cell and array performance is calculated via finite element method simulations. Results are confirmed via good agreement with finite difference, time domain calculations. Array performance exhibiting grating-lobe-free, broadside-scanned, dual-CP radiation with better than 20 dB return loss and over 75% aperture efficiency is presented.

Heat transfer in microwave heating

Heat transfer is considered as one of the most critical issues for design and implement of large-scale microwave heating systems, in which improvement of the microwave absorption of materials and suppression of uneven temperature distribution are the two main objectives. The present work focuses on the analysis of heat transfer in microwave heating for achieving highly efficient microwave assisted steelmaking through the investigations on the following aspects: (1) characterization of microwave dissipation using the derived equations, (2) quantification of magnetic loss, (3) determination of microwave absorption properties of materials, (4) modeling of microwave propagation, (5) simulation of heat transfer, and (6) improvement of microwave absorption and heating uniformity. Microwave heating is attributed to the heat generation in materials, which depends on the microwave dissipation. To theoretically characterize microwave heating, simplified equations for determining the transverse electromagnetic mode (TEM) power penetration depth, microwave field attenuation length, and half-power depth of microwaves in materials having both magnetic and dielectric responses were derived. It was followed by developing a simplified equation for quantifying magnetic loss in materials under microwave irradiation to demonstrate the importance of magnetic loss in microwave heating. The permittivity and permeability measurements of various materials, namely, hematite, magnetite concentrate, wüstite, and coal were performed. Microwave loss calculations for these materials were carried out. It is suggested that magnetic loss can play a major role in the heating of magnetic dielectrics. Microwave propagation in various media was predicted using the finite-difference time-domain method. For lossy magnetic dielectrics, the dissipation of microwaves in the medium is ascribed to the decay of both electric and magnetic fields. The heat transfer process in microwave heating of magnetite, which is a typical magnetic dielectric, was simulated by using an explicit finite-difference approach. It is demonstrated that the heat generation due to microwave irradiation dominates the initial temperature rise in the heating and the heat radiation heavily affects the temperature distribution, giving rise to a hot spot in the predicted temperature profile. Microwave heating at 915 MHz exhibits better heating homogeneity than that at 2450 MHz due to larger microwave penetration depth. To minimize/avoid temperature nonuniformity during microwave heating the optimization of object dimension should be considered. The calculated reflection loss over the temperature range of heating is found to be useful for obtaining a rapid optimization of absorber dimension, which increases microwave absorption and achieves relatively uniform heating. To further improve the heating effectiveness, a function for evaluating absorber impedance matching in microwave heating was proposed. It is found that the maximum absorption is associated with perfect impedance matching, which can be achieved by either selecting a reasonable sample dimension or modifying the microwave parameters of the sample.

Optimization of a label-free biosensor vertically characterized based on a periodic lattice of high aspect ratio SU-8 nano-pillars with a simplified 2D theoretical model

We have recently demonstrated a biosensor based on a lattice of SU8 pillars on a 1 μm SiO2/Si wafer by measuring vertically reflectivity as a function of wavelength. The biodetection has been proven with the combination of Bovine Serum Albumin (BSA) protein and its antibody (antiBSA). A BSA layer is attached to the pillars; the biorecognition of antiBSA involves a shift in the reflectivity curve, related with the concentration of antiBSA. A detection limit in the order of 2 ng/ml is achieved for a rhombic lattice of pillars with a lattice parameter (a) of 800 nm, a height (h) of 420 nm and a diameter(d) of 200 nm. These results correlate with calculations using 3D-finite difference time domain method. A 2D simplified model is proposed, consisting of a multilayer model where the pillars are turned into a 420 nm layer with an effective refractive index obtained by using Beam Propagation Method (BPM) algorithm. Results provided by this model are in good correlation with experimental data, reaching a reduction in time from one day to 15 minutes, giving a fast but accurate tool to optimize the design and maximizing sensitivity, and allows analyzing the influence of different variables (diameter, height and lattice parameter). Sensitivity is obtained for a variety of configurations, reaching a limit of detection under 1 ng/ml. Optimum design is not only chosen because of its sensitivity but also its feasibility, both from fabrication (limited by aspect ratio and proximity of the pillars) and fluidic point of view. (© 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)