Photo-acoustic tomography in a rotating setting

Photo-acoustic tomography is a coupled-physics (hybrid) medical imaging modality that aims to reconstruct optical parameters in biological tissues from ultrasound measurements. As propagating light gets partially absorbed, the resulting thermal expansion generates minute ultrasonic signals (the photo-acoustic effect) that are measured at the boundary of a domain of interest. Standard inversion procedures first reconstruct the source of radiation by an inverse ultrasound (boundary) problem and second describe the optical parameters from internal information obtained in the first step. This paper considers the rotating experimental setting. Light emission and ultrasound measurements are fixed on a rotating gantry, resulting in a rotation-dependent source of ultrasound. The two-step procedure we just mentioned does not apply. Instead, we propose an inversion that directly aims to reconstruct the optical parameters quantitatively. The mapping from the unknown (absorption and diffusion) coefficients to the ultrasound measurement via the unknown ultrasound source is modeled as a composition of a pseudo-differential operator and a Fourier integral operator. We show that for appropriate choices of optical illuminations, the above composition is an elliptic Fourier integral operator. Under the assumption that the coefficients are unknown on a sufficiently small domain, we derive from this a (global) injectivity result (measurements uniquely characterize our coefficients) combined with an optimal stability estimate. The latter is the same as that obtained in the standard (non-rotating experimental) setting

Simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection

Through the assimilation of present-day mantle seismic structure, adjoint methods can be used to constrain the structure of the mantle at earlier times, i.e., mantle initial conditions. However, the application to geophysical problems is restricted through both the high computational expense from repeated iteration between forward and adjoint models and the need to know mantle properties (such as viscosity and the absolute magnitude of temperature or density) a priori. We propose that an optimal first guess to the initial condition can be obtained through a simple backward integration (SBI) of the governing equations, thus lessening the computational expense. Given a model with known mantle properties, we show that a solution based on an SBI-generated first guess has smaller residuals than arbitrary guesses. Mantle viscosity and the effective Rayleigh number are crucial for mantle convection models, neither of which is exactly known. We place additional constraints on these basic mantle properties when the convection-induced dynamic topography on Earth's surface is considered within an adjoint inverse method. Besides assimilating present-day seismic structure as a constraint, we use dynamic topography and its rate of change in an inverse method that allows simultaneous inversion of the absolute upper and lower mantle viscosities, scaling between seismic velocity and thermal anomalies, and initial condition. The theory is derived from the governing equations of mantle convection and validated by synthetic experiments for both one-layer viscosity and two-layer viscosity regionally bounded spherical shells. For the one-layer model, at any instant of time, the magnitude of dynamic topography is controlled by the temperature scaling while the rate of change of topography is controlled by the absolute value of viscosity. For the two-layer case, the rate of change of topography constrains upper mantle viscosity while the magnitude of dynamic topography determines the temperature scaling (lower mantle viscosity) when upper-mantle (lower-mantle) density anomaly dominates the flow field; this two-stage scheme minimizes the tradeoff between temperature and lower mantle viscosity. For both cases, we show that the theory can constrain mantle properties with errors arising through the adjoint recovery of the initial condition; for the two-layer model, this error is manifest as a tradeoff between the temperature scaling and lower mantle viscosity

Monitoring thermal ablation via microwave tomography. An ex vivo experimental assessment

Thermal ablation treatments are gaining a lot of attention in the clinics thanks to their reduced invasiveness and their capability of treating non-surgical patients. The effectiveness of these treatments and their impact in the hospital's routine would significantly increase if paired with a monitoring technique able to control the evolution of the treated area in real-time. This is particularly relevant in microwave thermal ablation, wherein the capability of treating larger tumors in a shorter time needs proper monitoring. Current diagnostic imaging techniques do not provide effective solutions to this issue for a number of reasons, including economical sustainability and safety. Hence, the development of alternative modalities is of interest. Microwave tomography, which aims at imaging the electromagnetic properties of a target under test, has been recently proposed for this scope, given the significant temperature-dependent changes of the dielectric properties of human tissues induced by thermal ablation. In this paper, the outcomes of the first ex vivo experimental study, performed to assess the expected potentialities of microwave tomography, are presented. The paper describes the validation study dealing with the imaging of the changes occurring in thermal ablation treatments. The experimental test was carried out on two ex vivo bovine liver samples and the reported results show the capability of microwave tomography of imaging the transition between ablated and untreated tissue. Moreover, the discussion section provides some guidelines to follow in order to improve the achievable performances

Cloaking via mapping for the heat equation

This paper explores the concept of near-cloaking in the context of time-dependent heat propagation. We show that after the lapse of a certain threshold time instance, the boundary measurements for the homogeneous heat equation are close to the cloaked heat problem in a certain Sobolev space norm irrespective of the density-conductivity pair in the cloaked region. A regularised transformation media theory is employed to arrive at our results. Our proof relies on the study of the long time behaviour of solutions to the parabolic problems with high contrast in density and conductivity coefficients. It further relies on the study of boundary measurement estimates in the presence of small defects in the context of steady conduction problem. We then present some numerical examples to illustrate our theoretical results.Comment: 31 pages, 7 figures, 1 tabl

Quantum state tomography with non-instantaneous measurements, imperfections and decoherence

Tomography of a quantum state is usually based on positive operator-valued measure (POVM) and on their experimental statistics. Among the available reconstructions, the maximum-likelihood (MaxLike) technique is an efficient one. We propose an extension of this technique when the measurement process cannot be simply described by an instantaneous POVM. Instead, the tomography relies on a set of quantum trajectories and their measurement records. This model includes the fact that, in practice, each measurement could be corrupted by imperfections and decoherence, and could also be associated with the record of continuous-time signals over a finite amount of time. The goal is then to retrieve the quantum state that was present at the start of this measurement process. The proposed extension relies on an explicit expression of the likelihood function via the effective matrices appearing in quantum smoothing and solutions of the adjoint quantum filter. It allows to retrieve the initial quantum state as in standard MaxLike tomography, but where the traditional POVM operators are replaced by more general ones that depend on the measurement record of each trajectory. It also provides, aside the MaxLike estimate of the quantum state, confidence intervals for any observable. Such confidence intervals are derived, as the MaxLike estimate, from an asymptotic expansion of multi-dimensional Laplace integrals appearing in Bayesian Mean estimation. A validation is performed on two sets of experimental data: photon(s) trapped in a microwave cavity subject to quantum non-demolition measurements relying on Rydberg atoms; heterodyne fluorescence measurements of a superconducting qubit.Comment: 11 pages, 4 figures, submitte

Quantitative photoacoustic tomography with piecewise constant material parameters

The goal of quantitative photoacoustic tomography is to determine optical and acoustical material properties from initial pressure maps as obtained, for instance, from photoacoustic imaging. The most relevant parameters are absorption, diffusion and Grueneisen coefficients, all of which can be heterogeneous. Recent work by Bal and Ren shows that in general, unique reconstruction of all three parameters is impossible, even if multiple measurements of the initial pressure (corresponding to different laser excitation directions at a single wavelength) are available. Here, we propose a restriction to piecewise constant material parameters. We show that in the diffusion approximation of light transfer, piecewise constant absorption, diffusion and Gr\"uneisen coefficients can be recovered uniquely from photoacoustic measurements at a single wavelength. In addition, we implemented our ideas numerically and tested them on simulated three-dimensional data

In-situ steel solidification imaging in continuous casting using magnetic induction tomography

: Solidification process in continuous casting is a critical part of steel production. The speed and quality of the solidification process determines the quality of final product. Computational fluid dynamics (CFD) simulations are often used to describe the process and design of its control system, but so far, there is no any tool that provides an on-line measurement of the solidification front of hot steel during the continuous casting process. This paper presents a new tool based on magnetic induction tomography (MIT) for real time monitoring of this process. The new MIT system was installed at the end of the secondary cooling chamber of a casting unit and tested during several days in a real production process. MIT is able to create an internal map of electrical conductivity of hot steel deep inside the billet. The image of electrical conductivity is then converted to temperature profile that allows the measurement of the solid, mushy and liquid layers. In this study, such a conversion is done by synchronizing in one time step the MIT measurement and the thermal map generated with the actual process parameters available at that time. The MIT results were then compared with the results obtained of the CFD and thermal modelling of the industrial process. This is the first in-situ monitoring of the interior structure during a real continuous casting.The SHELL-THICK project has received funding from EU Research Fund for Coal and Steel under grant number 709830. This study reflects only the author's views and the European Commission is not responsible for any use that may be made of the information contained therein

Slabs in the lower mantle and their modulation of plume formation

Numerical mantle convection models indicate that subducting slabs can reach the core-mantle boundary (CMB) for a wide range of assumed material properties and plate tectonic histories. An increase in lower mantle viscosity, a phase transition at 660 km depth, depth-dependent thermal expansivity, and depth-dependent thermal diffusivity do not preclude model slabs from reaching the CMB. We find that ancient slabs could be associated with lateral temperature anomalies ~500°C cooler than ambient mantle. Plausible increases of thermal conductivity with depth will not cause slabs to diffuse away. Regional spherical models with actual plate evolutionary models show that slabs are unlikely to be continuous from the upper mantle to the CMB, even for radially simple mantle structures. The observation from tomography showing only a few continuous slab-like features from the surface to the CMB may be a result of complex plate kinematics, not mantle layering. There are important consequences of deeply penetrating slabs. Our models show that plumes preferentially develop on the edge of slabs. In areas on the CMB free of slabs, plume formation and eruption are expected to be frequent while the basal thermal boundary layer would be thin. However, in areas beneath slabs, the basal thermal boundary layer would be thicker and plume formation infrequent. Beneath slabs, a substantial amount of hot mantle can be trapped over long periods of time, leading to “mega-plume” formation. We predict that patches of low seismic velocity may be found beneath large-scale high seismic velocity structures at the core-mantle boundary. We find that the location, buoyancy, and geochemistry of mega-plumes will differ from those plumes forming at the edge of slabs. Various geophysical and geochemical implications of this finding are discussed