Source author record

Fernando L. Teixeira

Fernando L. Teixeira appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2019arXiv

Performance Analysis and Dynamic Evolution of Deep Convolutional Neural Network for Nonlinear Inverse Scattering

The solution of nonlinear electromagnetic (EM) inverse scattering problems is typically hindered by several challenges such as ill-posedness, strong nonlinearity, and high computational costs. Recently, deep learning has been demonstrated to be a promising tool in addressing these challenges. In particular, it is possible to establish a connection between a deep convolutional neural network (CNN) and iterative solution methods of nonlinear EM inverse scattering. This has led to the development of an efficient CNN-based solution to nonlinear EM inverse problems, termed DeepNIS. It has been shown that DeepNIS can outperform conventional nonlinear inverse scattering methods in terms of both image quality and computational time. In this work, we quantitatively evaluate the performance of DeepNIS as a function of the number of layers using structure similarity measure (SSIM) and mean-square error (MSE) metrics. In addition, we probe the dynamic evolution behavior of DeepNIS by examining its near-isometry property. It is shown that after a proper training stage the proposed CNN is near optimal in terms of the stability and generalization ability.

preprint2019arXiv

Quantum Information Propagation Preserving Computational Electromagnetics

We propose a new methodology, called numerical canonical quantization, to solve quantum Maxwell's equations useful for mathematical modeling of quantum optics physics, and numerical experiments on arbitrary passive and lossless quantum-optical systems. It is based on: (1) the macroscopic (phenomenological) electromagnetic theory on quantum electrodynamics (QED), and (2) concepts borrowed from computational electromagnetics. It was shown that canonical quantization in inhomogeneous dielectric media required definite and proper normal modes. Here, instead of ad-hoc analytic normal modes, we numerically construct complete and time-reversible normal modes in the form of traveling waves to diagonalize the Hamiltonian. Specifically, we directly solve the Helmholtz wave equations for a general linear, reciprocal, isotropic, non-dispersive, and inhomogeneous dielectric media by using either finite-element or finite-difference methods. To convert a scattering problem with infinite number of modes into one with a finite number of modes, we impose Bloch-periodic boundary conditions. This will sparsely sample the normal modes with numerical Bloch-Floquet-like normal modes. Subsequent procedure of numerical canonical quantization is straightforward using linear algebra. We provide relevant numerical recipes in detail and show an important numerical example of indistinguishable two-photon interference in quantum beam splitters, exhibiting Hong-Ou-Mandel effect, which is purely a quantum effect. Also, the present methodology provides a way of numerically investigating existing or new macroscopic QED theories. It will eventually allow quantum-optical numerical experiments of high fidelity to replace many real experiments as in classical electromagnetics.

preprint2016arXiv

Exact Charge-Conserving Scatter-Gather Algorithm for Particle-in-Cell Simulations on Unstructured Grids: A Geometric Perspective

We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This algorithm recovers, at a fundamental level, the scatter-gather algorithms presented recently by Campos-Pinto et al. [1] (to first-order) and by Squire et al. [2], but it is derived here in a streamlined fashion from a geometric viewpoint. Some ingredients reflecting this viewpoint are (1) the use of (discrete) differential forms of various degrees to represent fields, currents, and charged particles and provide localization rules for the degrees of freedom thereof on the various grid elements (nodes, edges, facets), (2) use of Whitney forms as basic interpolants from discrete differential forms to continuum space, and (3) use of a Galerkin formula for the discrete Hodge star operators (i.e., "mass matrices" incorporating the metric datum of the grid) applicable to generally irregular, unstructured grids. The expressions obtained for the scatter charges and scatter currents are very concise and do not involve numerical quadrature rules. Appropriate fractional areas within each grid element are identified that represent scatter charges and scatter currents within the element, and a simple geometric representation for the (exact) charge conservation mechanism is obtained by such identification. The field update is based on the coupled first-order Maxwell's curl equations to avoid spurious modes with secular growth (otherwise present in formulations that discretize the second-order wave equation). Examples are provided to verify preservation of discrete Gauss' law for all times.

preprint2016arXiv

Interface-Flattening Transform for EM Field Modeling in Tilted, Cylindrically-Stratified Geophysical Media

We propose and investigate an "interface-flattening" transformation, hinging upon Transformation Optics (T.O.) techniques, to facilitate the rigorous analysis of electromagnetic (EM) fields radiated by sources embedded in tilted, cylindrically-layered geophysical media. Our method addresses the major challenge in such problems of appropriately approximating the domain boundaries in the computational model while, in a full-wave manner, predicting the effects of tilting in the layers. When incorporated into standard pseudo-analytical algorithms, moreover, the proposed method is quite robust, as it is not limited by absorption, anisotropy, and/or eccentering profile of the cylindrical geophysical formations, nor is it limited by the radiation frequency. These attributes of the proposed method are in contrast to past analysis methods for tilted-layer media that often place limitations on the source and medium characteristics. Through analytical derivations as well as a preliminary numerical investigation, we analyze and discuss the method's strengths and limitations.

preprint2016arXiv

Perfectly Reflectionless Omnidirectional Absorbers and Electromagnetic Horizons

We demonstrate the existence of metamaterial blueprints describing, and fundamental limitations concerning, perfectly reflectionless omnidirectional electromagnetic absorbers (PR-OEMA). Previous attempts to define PR-OEMA blueprints have led to active (gain), rather than passive, media. We explain this fact and unveil new, distinct limitations of true PR-OEMA devices, including the appearance of an "electromagnetic horizon" on physical solutions. As practical alternatives we introduce alternative OEMA blueprints corresponding to media that, while not reflectionless, are nonetheless effective in absorbing incident waves in a manner robust to incident wave diversity.

preprint2016arXiv

Space-Time Block Diagonalization for Frequency-Selective MIMO Broadcast Channels

The most relevant linear precoding method for frequency-flat MIMO broadcast channels is block diagonalization (BD) which, under certain conditions, attains the same nonlinear dirty paper coding channel capacity. However, BD is not easily translated to frequency-selective channels, since space-time information is required for transceiver design. In this paper, we demonstrate that BD is feasible in frequency-selective MIMO broadcast channels to eliminate inter-user interference (IUI) if the transmit block length is sufficiently large, and if the number of transmit antennas is greater than the number of users. We also propose three different approaches to mitigate/eliminate inter-symbol interference (ISI) in block transmissions: i) time-reversal-based BD (TRBD) which maximizes spatial focusing around the receivers using transmitter processing only, ii) equalized BD (EBD) which minimizes the ISI using transmitter processing only, and iii) joint processing BD (JPBD), which uses linear processing at the transmitter and the receiver to suppress ISI. We analyze the theoretical diversity and multiplexing gains of these techniques, and we demonstrate that JPBD approximates full multiplexing gain for a sufficiently large transmit block length. Extensive numerical simulations show that the achievable rate and probability of error performance of all the proposed methods improve those of conventional time-reversal beamforming. Moreover, JPBD provides the highest achievable rate region for frequency-selective MIMO broadcast channels.

preprint2016arXiv

Spectral-Domain Computation of Fields Radiated by Sources in Non-Birefringent Anisotropic Media

We derive the key expressions to robustly address the eigenfunction expansion-based analysis of electromagnetic (EM) fields produced by current sources within planar non-birefringent anisotropic medium (NBAM) layers. In NBAM, the highly symmetric permeability and permittivity tensors can induce directionally-dependent, but polarization independent, propagation properties supporting "degenerate" characteristic polarizations, i.e. four linearly-independent eigenvectors associated with only two (rather than four) unique, non-defective eigenvalues. We first explain problems that can arise when the source(s) specifically reside within NBAM planar layers when using canonical field expressions. To remedy these problems, we exhibit alternative spectral-domain field expressions, immune to such problems, that form the foundation for a robust eigenfunction expansion-based analysis of time-harmonic EM radiation and scattering within such type of planar-layered media. Numerical results demonstrate the high accuracy and stability achievable using this algorithm.

preprint2016arXiv

Spectral-Domain-Based Scattering Analysis of Fields Radiated by Distributed Sources in Planar-Stratified Environments with Arbitrarily Anisotropic Layers

We discuss the numerically stable, spectral-domain computation and extraction of the scattered electromagnetic field excited by distributed sources embedded in planar-layered environments, where each layer may exhibit arbitrary and independent electrical and magnetic anisotropic response and loss profiles. This stands in contrast to many standard spectral-domain algorithms that are restricted to computing the fields radiated by Hertzian dipole sources in planar-layered environments where the media possess azimuthal-symmetric material tensors (i.e., isotropic, and certain classes of uniaxial, media). Although computing the scattered field, particularly when due to distributed sources, appears (from the analytical perspective, at least) relatively straightforward, different procedures within the computation chain, if not treated carefully, are inherently susceptible to numerical instabilities and (or) accuracy limitations due to the potential manifestation of numerically overflown and (or) numerically unbalanced terms entering the chain. Therefore, primary emphasis herein is given to effecting these tasks in a numerically stable and robust manner for all ranges of physical parameters. After discussing the causes behind, and means to mitigate, these sources of numerical instability, we validate the algorithm's performance against closed-form solutions. Finally, we validate and illustrate the applicability of the proposed algorithm in case studies concerning active remote sensing of marine hydrocarbon reserves embedded deep within lossy, planar-layered media.

preprint2016arXiv

Tensor Green's Function Evaluation in Arbitrarily Anisotropic, Layered Media using Complex-Plane Gauss-Laguerre Quadrature

We discuss the application of Complex-Plane Gauss-Laguerre Quadrature (CGLQ) to efficiently evaluate two-dimensional Fourier integrals arising as the solution to electromagnetic fields radiated by elementary dipole antennas embedded within planar-layered media with arbitrary material parameters. More specifically, we apply CGLQ to the long-standing problem of rapidly and efficiently evaluating the semi-infinite length "tails" of the Fourier integral path while simultaneously and robustly guaranteeing absolute, exponential convergence of the field solution despite diversity in the doubly anisotropic layer parameters, source type (i.e., electric or equivalent magnetic dipole), source orientation, observed field type (magnetic or electric), (non-zero) frequency, and (non-zero) source-observer separation geometry. The proposed algorithm exhibits robustness despite unique challenges arising for the fast evaluation of such two-dimensional integrals. Herein, we develop the mathematical treatment to rigorously evaluate the tail integrals using CGLQ and discuss and address the specific issues posed to the CGLQ method when anisotropic, layered media are present. To empirically demonstrate the CGLQ algorithm's computational efficiency, versatility, and accuracy, we perform a convergence analysis along with two case studies related to modeling of electromagnetic resistivity tools employed in geophysical prospection of layered, anisotropic Earth media and validating the ability of isoimpedance substrates to enhance the radiation performance of planar antennas placed in close proximity to metallic ground planes.

preprint2015arXiv

Full-Wave Algorithm to Model Effects of Bedding Slopes on the Response of Subsurface Electromagnetic Geophysical Sensors near Unconformities

We propose a full-wave pseudo-analytical numerical electromagnetic (EM) algorithm to model subsurface induction sensors, traversing planar-layered geological formations of arbitrary EM material anisotropy and loss, which are used, for example, in the exploration of hydrocarbon reserves. Unlike past pseudo-analytical planar-layered modeling algorithms that impose parallelism between the formation's bed junctions however, our method involves judicious employment of Transformation Optics techniques to address challenges related to modeling relative slope (i.e., tilting) between said junctions (including arbitrary azimuth orientation of each junction). The algorithm exhibits this flexibility, both with respect to loss and anisotropy in the formation layers as well as junction tilting, via employing special planar slabs that coat each "flattened" (i.e., originally tilted) planar interface, locally redirecting the incident wave within the coating slabs to cause wave fronts to interact with the flattened interfaces as if they were still tilted with a specific, user-defined orientation. Moreover, since the coating layers are homogeneous rather than exhibiting continuous material variation, a minimal number of these layers must be inserted and hence reduces added simulation time and computational expense. As said coating layers are not reflectionless however, they do induce artificial field scattering that corrupts legitimate field signatures due to the (effective) interface tilting. Numerical results, for two half-spaces separated by a tilted interface, quantify error trends versus material and sensor characteristics. We finally exhibit responses of sensors traversing three-layered media, where we vary the anisotropy, loss, and relative tilting of the formations and explore the sensitivity of the sensor's complex-valued measurements.

preprint2015arXiv

Interference-Nulling Time-Reversal Beamforming for mm-Wave Massive MIMO in Multi-User Frequency-Selective Indoor Channels

Millimeter wave (mm-wave) and massive MIMO have been proposed for next generation wireless systems. However, there are many open problems for the implementation of those technologies. In particular, beamforming is necessary in mm-wave systems in order to counter high propagation losses. However, conventional beamsteering is not always appropriate in rich scattering multipath channels with frequency selective fading, such as those found in indoor environments. In this context, time-reversal (TR) is considered a promising beamforming technique for such mm-wave massive MIMO systems. In this paper, we analyze a baseband TR beamforming system for mm-wave multi-user massive MIMO. We verify that, as the number of antennas increases, TR yields good equalization and interference mitigation properties, but inter-user interference (IUI) remains a main impairment. Thus, we propose a novel technique called interference-nulling TR (INTR) to minimize IUI. We evaluate numerically the performance of INTR and compare it with conventional TR and equalized TR beamforming. We use a 60 GHz MIMO channel model with spatial correlation based on the IEEE 802.11ad SISO NLoS model. We demonstrate that INTR outperforms conventional TR with respect to average BER per user and achievable sum rate under diverse conditions, providing both diversity and multiplexing gains simultaneously.

preprint2014arXiv

Computation of potentials from current electrodes in cylindrically stratified media: A stable, rescaled semi-analytical formulation

We present an efficient and robust semi-analytical formulation to compute the electric potential due to arbitrary-located point electrodes in three-dimensional cylindrically stratified media, where the radial thickness and the medium resistivity of each cylindrical layer can vary by many orders of magnitude. A basic roadblock for robust potential computations in such scenarios is the poor scaling of modified-Bessel functions used for computation of the semi-analytical solution, for extreme arguments and/or orders. To accommodate this, we construct a set of rescaled versions of modified-Bessel functions, which avoids underflows and overflows in finite precision arithmetic, and minimizes round-off errors. In addition, several extrapolation methods are applied and compared to expedite the numerical evaluation of the (otherwise slowly convergent) associated Sommerfeld-type integrals. The proposed algorithm is verified in a number of scenarios relevant to geophysical exploration, but the general formulation presented is also applicable to other problems governed by Poisson equation such as Newtonian gravity, heat flow, and potential flow in fluid mechanics, involving cylindrically stratified environments.