Research development of terahertz imaging systems

The electromagnetic radiation at a terahertz frequencies (from 0.1 THz to 10 THz) is situated in the frequency band comprised between the optical band and the radio band. The interest of the scientific community in this frequency band has grown up due to its large capabilities to develop innovative imaging systems. The terahertz waves are able to generate extremely short pulses that achieve good spatial resolution, good penetration capabilities and allow to identify microscopic structures using spectral analysis. The work carried out during the period of the grant has been based on the developement of system working at the aforementioned frequency band. The main system is based on a total power radiometer working at 0.1 THz to perform security imaging. Moreover, the development of this system has been useful to gain knowledge in the behavior of the component systems at this frequency band. Moreover, a vectorial network analyzer has been used to characterize materials and perform active raster imaging. A materials measurement system has been designed and used to measure material properties as permittivity, losses and water concentration. Finally, the design of a terahertz time-domain spectrometer (THz-TDS) system has been started. This system will allow to perform tomographic measurement with very high penetration resolutions while allowing the spectral characterization of the sample material. The application range of this kind of system is very wide: from the identification of cancerous tissues of a skin to the characterization of the thickness of a painted surface of a car.

On the analysis of travelling waves to a nonlinear flux limited reaction-diffusion equation

In this paper we study the existence and qualitative properties of travelling waves associated to a nonlinear flux limited partial differential equation coupled to a Fisher-Kolmogorov-Petrovskii-Piskunov type reaction term. We prove the existence and uniqueness of finite speed moving fronts of C2 classical regularity, but also the existence of discontinuous entropy travelling wave solutions.

Watching television: a previously unrecognized powerful trigger of λ waves.

OBJECTIVE: To assess whether lambda waves are elicited by watching television (TV) and their association with demographical and EEG features. METHODS: We retrospectively compared lambda wave occurrence in prolonged EEG monitorings of outpatients who were allowed to watch TV and in standard EEGs recorded in TV-free rooms. All EEGs were interpreted by the same two electroencephalographers. RESULTS: Of 2,072 standard EEG reports, 36 (1.7 %) mentioned lambda waves versus 46 (32.2%) of 143 prolonged EEG monitoring reports (P < 0.001). Multivariable comparison of prolonged EEG monitorings and standard EEGs disclosed that recordings performed in rooms with a TV (odds ratio, 20.6; 95% confidence interval, 4.8-88.0) and normal EEGs (odds ratio, 3.03; 95% confidence interval, 1.5-6.25) were independently associated with lambda waves. In the prolonged EEG monitoring group, all recordings with lambda waves also had positive occipital sharp transients of sleep. CONCLUSIONS: Watching TV likely represents a powerful and previously unrecognized stimulus for lambda waves. Furthermore, this study confirms the benign nature of this EEG variant and its strong association with positive occipital sharp transients of sleep.

Comparison of three commercially available radio frequency coils for human brain imaging at 3 Tesla.

OBJECTIVE: To evaluate a transverse electromagnetic (TEM), a circularly polarized (CP) (birdcage), and a 12-channel phased array head coil at the clinical field strength of B0 = 3T in terms of signal-to-noise ratio (SNR), signal homogeneity, and maps of the effective flip angle alpha. MATERIALS AND METHODS: SNR measurements were performed on low flip angle gradient echo images. In addition, flip angle maps were generated for alpha(nominal) = 30 degrees using the double angle method. These evaluation steps were performed on phantom and human brain data acquired with each coil. Moreover, the signal intensity variation was computed for phantom data using five different regions of interest. RESULTS: In terms of SNR, the TEM coil performs slightly better than the CP coil, but is second to the smaller 12-channel coil for human data. As expected, both the TEM and the CP coils show superior image intensity homogeneity than the 12-channel coil, and achieve larger mean effective flip angles than the combination of body and 12-channel coil with reduced radio frequency power deposition. CONCLUSION: At 3T the benefits of TEM coil design over conventional lumped element(s) coil design start to emerge, though the phased array coil retains an advantage with respect to SNR performance.

Could radiotherapy effectiveness be enhanced by electromagnetic field treatment?

One of the main goals in radiobiology research is to enhance radiotherapy effectiveness without provoking any increase in toxicity. In this context, it has been proposed that electromagnetic fields (EMFs), known to be modulators of proliferation rate, enhancers of apoptosis and inductors of genotoxicity, might control tumor recruitment and, thus, provide therapeutic benefits. Scientific evidence shows that the effects of ionizing radiation on cellular compartments and functions are strengthened by EMF. Although little is known about the potential role of EMFs in radiotherapy (RT), the radiosensitizing effect of EMFs described in the literature could support their use to improve radiation effectiveness. Thus, we hypothesized that EMF exposure might enhance the ionizing radiation effect on tumor cells, improving the effects of RT. The aim of this paper is to review reports of the effects of EMFs in biological systems and their potential therapeutic benefits in radiotherapy.

Two waves of recombinase gene expression in developing thymocytes.

During T cell development in the thymus, T cell receptor (TCR) alpha, beta, gamma, and delta genes are rearranged and expressed. TCR rearrangement strictly depends upon the coordinate activity of two recombinase activating genes, Rag-1 and Rag-2. In this study we have followed the expression of these genes at different stages of intrathymic development. The results indicate that there are two periods of high Rag-1 and Rag-2 mRNA expression. The first wave peaks early at the CD25+CD4-CD8-CD3- stage of development and coincides with the initial appearance of transcripts derived from fully rearranged TCR beta, gamma, and delta genes, whereas the second wave occurs later at the CD4+CD8+ stage coincident with full-length TCR alpha mRNA expression. Active downregulation of Rag-1 and Rag-2 mRNA expression appears to occur in vivo between the two peaks of recombinase activity. This phenomenon can be mimicked in vitro in response to artificial stimuli such as phorbol myristate acetate and calcium ionophore. Collectively our data suggest that recombinase expression is actively regulated during early thymus development independently of cell surface expression of a mature heterodimeric TCR protein complex.

Reaction-diffusion waves of advance in the transition to agricultural economics

In a previous paper [J.Fort and V.Méndez, Phys. Rev. Lett. 82, 867 (1999)], the possible importance of higher-order terms in a human population wave of advance has been studied. However, only a few such terms were considered. Here we develop a theory including all higher-order terms. Results are in good agreement with the experimental evidence involving the expansion of agriculture in Europe

Investigation of Electromagnetic Gauges for Determining In-Place HMA Density, 2007

Density is an important component of hot-mix asphalt (HMA) pavement quality and long-term performance. Insufficient density of an in-place HMA pavement is the most frequently cited construction-related performance problem. This study evaluated the use of electromagnetic gauges to nondestructively determine densities. Field and laboratory measurements were taken with two electromagnetic gauges—a PaveTracker and a Pavement Quality Indicator (PQI). Test data were collected in the field during and after paving operations and also in a laboratory on field mixes compacted in the lab. This study revealed that several mix- and project-specific factors affect electromagnetic gauge readings. Consequently, the implementation of these gauges will likely need to be done utilizing a test strip on a project- and mix-specific basis to appropriately identify an adjustment factor for the specific electromagnetic gauge being used for quality control and quality assurance (QC/QA) testing. The substantial reduction in testing time that results from employing electromagnetic gauges rather than coring makes it possible for more readings to be used in the QC/QA process with real-time information without increasing the testing costs.

Endogenous spillovers, increased competition and re-organization waves

We consider an entrepreneur that is the sole producer of a costreducing skill, but the entrepreneur that hires a team to usethe skill cannot prevent collusive trade for the innovation related knowledge between employees and competitors. We showthat there are two types of diffusion avoiding strategies forthe entrepreneur to preempt collusive communication i) settingup a large productive capacity (the traditional firm) and ii)keeping a small team (the lean firm). The traditional firm ischaracterized by its many "marginal" employees that work shortdays, receive flat wages and are incompletely informed about the innovation. The lean firm is small in number of employees,engages in complete information sharing among members, that are paid with stock option schemes. We find that the lean firm is superior to the traditional firm when technological entry costsare low and when the sector is immature.

Feasibility and electromagnetic compatibility study of the ClearPEM front-end electronics for simultaneous PET-MR imaging

In this work we present a first feasibility study of the ClearPEM technology for simultaneous PET-MR imaging. The mutual electromagnetic interference (EMI) effects between both systems were evaluated on a 7 T magnet by characterizing the response behavior of the ClearPEM detectors and front-end electronics to pulsed RF power and switched magnetic field gradients; and by analyzing the MR system performance degradation from noise pickup into the RF receiver chain, and from magnetic susceptibility artifacts caused by PET front-end materials.

The geophysical electromagnetic prospecting in the spatial location of the travertine deposits of the Banyoles Depression (Girona)

The results of the application of the geophysical electromagnetic prospection methods in the resolution of the problems of the spatial location of the travertine quaternary formations of the Banyoles depression are presented

Electromagnetic radiation initiated by hadronic jets from microquasars in the ISM

Microquasars are potential candidates to produce a non-negligible fraction of the observed galactic cosmic rays. The protons accelerated at the jet termination shock interact with the interstellar medium and may produce detectable fluxes of extended emission at different energy bands: high-energy and very high-energy gamma-rays produced by neutral pion-decay, synchrotron and bremsstrahlung emission in a wide energy range generated by the secondary electrons produced by charged pion-decay. We discuss the association between this scenario and some of the unidentified EGRET sources in the galactic plane.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.