Inverse scattering of 2d photonic structures by layer-stripping

Design and reconstruction of 2d and 3d photonic structures are usually carried out by forward simulations combined with optimization or intuition. Reconstruction by means of layer-stripping has been applied in seismic processing as well as in design and characterization of 1d photonic structures such as fiber Bragg gratings. Layer-stripping is based on causality, where the earliest scattered light is used to recover the structure layer-by-layer. Our set-up is a 2d layered nonmagnetic structure probed by plane polarized harmonic waves entering normal to the layers. It is assumed that the dielectric permittivity in each layer only varies orthogonal to the polarization. Based on obtained reflectance data covering a suitable frequency interval, time-localized pulse data are synthesized and applied to reconstruct the refractive index profile in the leftmost layer by identifying the local, time-domain Fresnel reflection at each point. Once the first layer is known, its impact on the reflectance data is stripped off, and the procedure repeated for the next layer. Through numerical simulations it will be demonstrated that it is possible to reconstruct structures consisting of several layers. The impact of evanescent modes and limited bandwidth is discussed

IBM Watson Supporting Space Radiation

The NASA Human Research Program (HRP) Space Radiation (SR) Program Element has been working with IBM Watson Explorer (WEX) to create a tool that allows researchers to search the NASA SR-funded research corpus to help streamline research and maximize efficiency. The entire corpus of publications from research funded by the NASA SR Program Element has been ingested into WEX to allow for examination of: synergies across funded research areas, gaps in research, and collaboration opportunities. This information will be valuable to both scientists and managers as it will allow analysis related to specific scientific questions, inform key decisions and support cross validation of study results. NASA will also evaluate the potential to make WEX publicly available in order to facilitate proposal generation and to enhance collaborations within and across disciplines

Blind deconvolution of sources from the transmission responses of one-dimensional inhomogeneous continuous and discrete layered absorbing media

Non-destructive testing and reconstruction of lossy non-uniform transmission lines require solving inverse scattering problems for inhomogeneous one-dimensional (1D) lossy media. Often the source waveform used to probe the medium is unknown, since it cannot be measured separately. This paper shows that the blind deconvolution problem of reconstructing both an unknown source and an unknown lossy 1D medium can be solved, provided the medium absorption is sufficiently large. We generalize the minimum-phase property of the transmission response to a more severe property for lossy media. We also show how partial knowledge of the source can be used to reconstruct low-loss and lossless 1D media. Both discrete and continuous 1D layered media are considered. Numerical examples are presented.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49102/2/ip7417.pd

Connections between three-dimensional inverse scattering and linear least-squares estimation of random fields

The three-dimensional Schrödinger equation inverse scattering problem with a nonspherically-symmetric potential is related to the filtering problem of computing the linear leastsquares estimate of the three-dimensional random field on the surface of a sphere from noisy observations inside the sphere. The relation consists of associating an estimation problem with the inverse scattering problem, and vice-versa. This association allows equations and quantities for one problem to be given interpretations in terms of the other problem. A new fast algorithm is obtained for the estimation of random fields using this association. The present work is an extension of the connections between estimation and inverse scattering already known to exist for stationary random processes and one-dimensional scattering potentials, and isotropic random fields and radial scattering protentials.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/41642/1/10440_2004_Article_BF00046966.pd

Differential and integral methods for three-dimensional inverse scattering problems with a non-local potential

The inverse scattering problem for the three-dimensional Schrodinger equation with a non-local potential diagonal in the radius variable is considered. Generalised Gel'fand-Levitan and Marchenko integral equations are derived using Gram-Schmidt orthogonalisations of free-space wave solutions. A generalised miracle equation is then used to recover the non-local potential. This is the first unified presentation of both the generalised Gel'fand-Levitan and Marchenko equations for this problem. New differential fast algorithms for solving these integral equations and recovering the potential are presented. An application to the problem of linear least-squares estimation of a random field is noted.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49089/2/ipv4i2p549.pd

Discrete matrix Schrodinger equation equivalents of discrete two-component matrix wave systems

Multichannel discrete lossless two-component wave systems describe physical scattering for multiple coupled transmission lines, for elastic media, and for electromagnetic wave propagation in layered media. Such a wave system is transformed into an equivalent discrete matrix Schrodinger equation that includes the effects of transmission losses and transmission scattering. These effects are not included in the continuous Schrodinger equation, so that discretisation of the latter omits significant effects. Applications include fast algorithms for synthesis and inversion of the reflection responses of multichannel scattering media. These algorithms require roughly half as many matrix multiplications as previous algorithms that employ the multichannel two-component wave system.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49091/2/ipv5i3p425.pd

Inverse scattering for time-varying one-dimensional layered media: algorithms and applications

We formulate and present algorithms and applications for the inverse scattering problem for a one-dimensional layered medium whose reflection coefficients vary with time. We show that given the reflection responses to impulsive plane waves incident on the medium at different times, this problem can be solved either by solving a series of nested non-Toeplitz systems of equations, or by propagating a set of coupled layer-stripping algorithms with reflection coefficients varying in time. All multiple scattering effects are included here. We show that our results reduce to previous results in the special case of time-invariant media. Applications include the two-dimensional (2D) inverse resistivity problem of reconstructing the 2D resistivity of a 2D medium from surface measurements, and the reconstruction of rapidly curing layered media. We also present a time-varying discrete Miura transform linking the time-varying discrete wavesystem inverse scattering problem to a time-varying discrete Schrödinger equation, and a new important feasibility condition on the reflection responses, which seems to be new in the context of the 2D inverse resistivity problem. A numerical example is provided.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49101/2/ip7318.pd

Discrete Gel'fand-Levitan and Marchenko matrix equations and layer stripping algorithms for the discrete two-dimensional Schrödinger equation inverse scattering problem with a nonlocal potential

We develop discrete counterparts to the Gel'fand-Levitan and Marchenko integral equations for the two-dimensional (2D) discrete inverse scattering problem in polar coordinates with a nonlocal potential. We also develop fast layer stripping algorithms that solve these systems of equations exactly. The significance of these results is: (1) they are the first numerical implementation of Newton's multidimensional inverse scattering theory; (2) they show that the result will almost always be a nonlocal potential, unless the data are `miraculous'; (3) they show that layer stripping algorithms implement fast `split' signal processing fast algorithms; (4) they link 2D discrete inverse scattering with 2D discrete random field linear least-squares estimation; and (5) they formulate and solve 2D discrete Schrödinger equation inverse scattering problems in polar coordinates.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49103/2/ip8322.pd

Inversion of the Bloch transform in magnetic resonance imaging using asymmetric two-component inverse scattering

In magnetic resonance imaging, the relation between the radio-frequency modulation of the magnetic field and the desired final magnetisation state is called the Bloch transform. Selective excitation then amounts to inverting this transform, which is highly nonlinear. Previous attempts to formulate this problem as an inverse scattering problem have restricted attention to solutions using reflectionless potentials. The author uses fast numerical algorithms for inverse scattering problems to obtain a much larger set of solutions. Numerical examples are included.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49094/2/ipv6i1p133.pd

COMMENT: Reply to a comment by R Newton

Following Newton's comment (see ibid., vol.5, p.437 (1989)), the author points out that: (1) the major contribution of his previous paper (see ibid., vol.4, p.549 (1988)) is the generalised fast Cholesky algorithm, not the integral equations; (2) the Marchenko integral equation is correct as written, but does not determine the scattered field; and (3) it is implicitly assumed that the regular solution exists. The last two points suggest that the generalised fast Cholesky algorithm may be the only viable solution to the non-local problem.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/49092/2/ipv5i3p439.pd