A robust and low frequency stable time domain PMCHWT equation

The time domain PMCHWT equation models transient scattering by piecewise homogeneous dielectrics. After discretization, it can be solved using the marching-on-in-time algorithm. Unfortunately, the PMCHWT equation suffers from DC instability: it supports constant in time regime solutions. Upon discretization, the corresponding poles of the system response function shift into the unstable region of the complex plane, rendering the MOT algorithm unstable. Furthermore, the discrete system becomes ill-conditioned when a large time step is used. This phenomenon is termed low frequency breakdown. In this contribution, the quasi Helmholtz components of the PMCHWT equation are separated using projector operators. Judicially integrating or differentiating these components of the basis and testing functions leads to an algorithm that (i) does not suffer from unstable modes even in the presence of moderate numerical errors, (ii) remains well-conditioned for large time steps, and (iii) can be applied effectively to both simply and multiply connected geometries

How do we understand and visualize uncertainty?

Geophysicists are often concerned with reconstructing subsurface properties using observations collected at or near the surface. For example, in seismic migration, we attempt to reconstruct subsurface geometry from surface seismic recordings, and in potential field inversion, observations are used to map electrical conductivity or density variations in geologic layers. The procedure of inferring information from indirect observations is called an inverse problem by mathematicians, and such problems are common in many areas of the physical sciences. The inverse problem of inferring the subsurface using surface observations has a corresponding forward problem, which consists of determining the data that would be recorded for a given subsurface configuration. In the seismic case, forward modeling involves a method for calculating a synthetic seismogram, for gravity data it consists of a computer code to compute gravity fields from an assumed subsurface density model. Note that forward modeling often involves assumptions about the appropriate physical relationship between unknowns (at depth) and observations on the surface, and all attempts to solve the problem at hand are limited by the accuracy of those assumptions. In the broadest sense then, exploration geophysicists have been engaged in inversion since the dawn of the profession and indeed algorithms often applied in processing centers can all be viewed as procedures to invert geophysical data

A space-time mixed Galerkin marching-on-in-time scheme for the time-domain combined field integral equation

The time domain combined field integral equation (TD-CFIE), which is constructed from a weighted sum of the time domain electric and magnetic field integral equations (TD-EFIE and TD-MFIE) for analyzing transient scattering from closed perfect electrically conducting bodies, is free from spurious resonances. The standard marching-on-in-time technique for discretizing the TD-CFIE uses Galerkin and collocation schemes in space and time, respectively. Unfortunately, the standard scheme is theoretically not well understood: stability and convergence have been proven for only one class of space-time Galerkin discretizations. Moreover, existing discretization schemes are nonconforming, i.e., the TD-MFIE contribution is tested with divergence conforming functions instead of curl conforming functions. We therefore introduce a novel space-time mixed Galerkin discretization for the TD-CFIE. A family of temporal basis and testing functions with arbitrary order is introduced. It is explained how the corresponding interactions can be computed efficiently by existing collocation-in-time codes. The spatial mixed discretization is made fully conforming and consistent by leveraging both Rao-Wilton-Glisson and Buffa-Christiansen basis functions and by applying the appropriate bi-orthogonalization procedures. The combination of both techniques is essential when high accuracy over a broad frequency band is required

A Stable Higher Order Space-Time Galerkin Scheme for Time Domain Integral Equations

Stability of time domain integral equation (TDIE) solvers has remained an elusive goal for many years. Advancement of this research has largely progressed on four fronts: (1) Exact integration, (2) Lubich quadrature, (3) smooth temporal basis functions, and (4) Space-time separation of convolutions with the retarded potential. The latter method was explored in [Pray et al. IEEE TAP 2012]. This method's efficacy in stabilizing solutions to the time domain electric field integral equation (TD-EFIE) was demonstrated on first order surface descriptions (flat elements) in tandem with 0th order functions as the temporal basis. In this work, we develop the methodology necessary to extend to higher order surface descriptions as well as to enable its use with higher order temporal basis functions. These higher order temporal basis functions are used in a Galerkin framework. A number of results that demonstrate convergence, stability, and applicability are presented.Comment: 8 pages, 12 figure

Accurate temporal discretization of time domain boundary integral equations

In this contribution, a novel temporal discretization scheme for time domain boundary integral equations is introduced. It distinguishes itself by (i) a new approach to the construction of higher order temporal basis functions, and (ii) the use of temporal Petrov-Galerkin testing as opposed to the widespread collocation in time method. The retarded potential boundary integral equation and its classic collocation in time discretization will be revisited. Next, the new temporal basis and testing functions will be introduced. The space-time Petrov-Galerkin discretization using these functions will be elucidated. Finally, numerical results are presented testifying to the improved accuracy of the novel scheme

Calderon multiplicative preconditioner for the PMCHWT equation applied to chiral media

In this contribution, a Calderon preconditioned algorithm for the modeling of scattering of time harmonic electromagnetic waves by a chiral body is introduced. The construction of the PMCHWT in the presence of chiral media is revisited. Since this equation reduces to the classic PMCHWT equation when the chirality parameter tends to zero, it shares its spectral properties. More in particular, it suffers from dense grid breakdown. Based on the work in [1], [2], a regularized version of the PMCHWT equation is introduced. A discretization scheme is described. Finally, the validity and spectral properties are studied numerically. More in particular, it is proven that linear systems arising in the novel scheme can be solved in a small number of iterations, regardless the mesh parameter

Nanobody technology : a versatile toolkit for microscopic imaging, protein-protein interaction analysis, and protein function exploration

Over the last two decades, nanobodies or single-domain antibodies have found their way in research, diagnostics, and therapy. These antigen-binding fragments, derived from Camelid heavy chain only antibodies, possess remarkable characteristics that favor their use over conventional antibodies or fragments thereof, in selected areas of research. In this review, we assess the current status of nanobodies as research tools in diverse aspects of fundamental research. We discuss the use of nanobodies as detection reagents in fluorescence microscopy and focus on recent advances in super-resolution microscopy. Second, application of nanobody technology in investigating protein-protein interactions is reviewed, with emphasis on possible uses in mass spectrometry. Finally, we discuss the potential value of nanobodies in studying protein function, and we focus on their recently reported application in targeted protein degradation. Throughout the review, we highlight state-of-the-art engineering strategies that could expand nanobody versatility and we suggest future applications of the technology in the selected areas of fundamental research

On the Hierarchical Preconditioning of the PMCHWT Integral Equation on Simply and Multiply Connected Geometries

We present a hierarchical basis preconditioning strategy for the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) integral equation considering both simply and multiply connected geometries.To this end, we first consider the direct application of hierarchical basis preconditioners, developed for the Electric Field Integral Equation (EFIE), to the PMCHWT. It is notably found that, whereas for the EFIE a diagonal preconditioner can be used for obtaining the hierarchical basis scaling factors, this strategy is catastrophic in the case of the PMCHWT since it leads to a severly ill-conditioned PMCHWT system in the case of multiply connected geometries. We then proceed to a theoretical analysis of the effect of hierarchical bases on the PMCHWT operator for which we obtain the correct scaling factors and a provably effective preconditioner for both low frequencies and mesh refinements. Numerical results will corroborate the theory and show the effectiveness of our approach