On the monitoring of surface displacement in connection with volcano reactivation in Tenerife, Canary Islands, using space techniques

Geodetic volcano monitoring in Tenerife has mainly focused on the Las Cañadas Caldera, where a geodetic micronetwork and a levelling profile are located. A sensitivity test of this geodetic network showed that it should be extended to cover the whole island for volcano monitoring purposes. Furthermore, InSAR allowed detecting two unexpected movements that were beyond the scope of the traditional geodetic network. These two facts prompted us to design and observe a GPS network covering the whole of Tenerife that was monitored in August 2000. The results obtained were accurate to one centimetre, and confirm one of the deformations, although they were not definitive enough to confirm the second one. Furthermore, new cases of possible subsidence have been detected in areas where InSAR could not be used to measure deformation due to low coherence. A first modelling attempt has been made using a very simple model and its results seem to indicate that the deformation observed and the groundwater level variation in the island may be related. Future observations will be necessary for further validation and to study the time evolution of the displacements, carry out interpretation work using different types of data (gravity, gases, etc) and develop models that represent the island more closely. The results obtained are important because they might affect the geodetic volcano monitoring on the island, which will only be really useful if it is capable of distinguishing between displacements that might be linked to volcanic activity and those produced by other causes. One important result in this work is that a new geodetic monitoring system based on two complementary techniques, InSAR and GPS, has been set up on Tenerife island. This the first time that the whole surface of any of the volcanic Canary Islands has been covered with a single network for this purpose. This research has displayed the need for further similar studies in the Canary Islands, at least on the islands which pose a greater risk of volcanic reactivation, such as Lanzarote and La Palma, where InSAR techniques have been used already.

Modelización de un proceso de biorremediación de suelos contaminados con gasoil

El objeto de esta tesis doctoral es encontrar, mediante herramientas estadísticas, un modelo matemático que prediga la evolución temporal de las concentraciones de los hidrocarburos totales del petróleo (TPH) en un proceso de biorremediación de suelos contaminados con gasoleo. Obtenido el modelo matemático realizar las comparaciones pertinentes para determinar la eficiencia de distintos tratamientos, la influencia tanto del suelo como del nivel de concentración. Para cumplir este objetivo se realizó un diseño de experimentos que tomaba en cuenta dos tipos de suelos, dos niveles de concentración de hidrocarburo y seis tratamientos. Se realizaron en total 1824 ensayos en laboratorio repartidos en 8 campañas de campo durante un año de experiencia tanto en parcelas de 4 x 4 m como en microcosmos. Los resultados obtenidos muestran que el proceso de biorremediación se puede ajustar a un comportamiento exponencial, que es una ecuación de cinética de primer orden, y que las comparaciones realizadas han mostrado resultados satisfactorios de la eficiencia de algunos tratamientos. ABSTRACT The purpose of this thesis is using statistical tools find a mathematical model to predict the time evolution of the total petroleum hydrocarbons (TPH) concentrations in the bioremediation process of diesel contaminated soil. Retrieved mathematical model make relevant comparisons to determine the efficiency of different treatments and the influence of both soil and concentration levels. To achieve this goal a design of experiments was developed, it took into account two types of soil, two concentration levels of hydrocarbon and six treatments. There were a total of 1824 laboratory tests spread over 8 field campaigns during a year of experience in both plots of 4 x 4 m in microcosm. The results obtained show that the bioremediation process can be adjusted to an exponential model, it is an equation of kinetic of first order and that the comparisons have shown satisfactory results of the efficiency of some treatments.

Diabatic quantum gates

At present, several models for quantum computation have been proposed. Adiabatic quantum computation scheme particularly offers this possibility and is based on a slow enough time evolution of the system, where no transitions take place. In this work, a new strategy for quantum computation is provided from the opposite point of view. The objective is to control the non-adiabatic transitions between some states in order to produce the desired exit states after the evolution. The model is introduced by means of an analogy between the adiabatic quantum computation and an inelastic atomic collision. By means of a simple two-state model, several quantum gates are reproduced, concluding the possibility of diabatic universal faulttolerant quantum computation. Going a step further, a new quantum diabatic computation model is glimpsed, where a carefully chosen Hamiltonian could carry out a non-adiabatic transition between the initial and the sought final state.

Simulation of academic performance rates, estimation of indicators of subjects and of a complete degree

The purpose of this paper is to present a program written in Matlab-Octave for the simulation of the time evolution of student curricula, i.e, how students pass their subjects along time until graduation. The program computes, from the simulations, the academic performance rates for the subjects of the study plan for each semester as well as the overall rates, which are a) the efficiency rate defined as the ratio of the number of students passing the exam to the number of students who registered for it and b) the success rate, defined as the ratio of the number of students passing the exam to the number of students who not only registered for it but also actually took it. Additionally, we compute the rates for the bachelor academic degree which are established for Spain by the National Quality Evaluation and Accreditation Agency (ANECA) and which are the graduation rate (measured as the percentage of students who finish as scheduled in the plan or taking an extra year) and the efficiency rate (measured as the percentage of credits which a student who graduated has really taken). The simulation is done in terms of the probabilities of passing all the subjects in their study plan. The application of the simulator to Polytech students in Madrid, where requirements for passing are specially stiff in first and second year subjects, is particularly relevant to analyze student cohorts and the probabilities of students finishing in the minimum of four years, or taking and extra year or two extra years, and so forth. It is a very useful tool when designing new study plans. The calculation of the probability distribution of the random variable "number of semesters a student has taken to complete the curricula and graduate" is difficult or even unfeasible to obtain analytically, and this is even truer when we incorporate uncertainty in parameter estimation. This is why we apply Monte Carlo simulation which not only provides illustration of the stochastic process but also a method for computation. The stochastic simulator is proving to be a useful tool for identification of the subjects most critical in the distribution of the number of semesters for curriculum vitae (CV) completion and subsequently for a decision making process in terms of CV planning and passing standards in the University. Simulations are performed through a graphical interface where also the results are presented in appropriate figures. The Project has been funded by the Call for Innovation in Education Projects of Universidad Politécnica de Madrid (UPM) through a Project of its school Escuela Técnica Superior de Ingenieros Industriales ETSII during the period September 2010-September 2011.

A molecular dynamics study of swift heavy ion irradiation of amorphous silica: the role of thermal effects

Irradiation with swift heavy ions (SHI), roughly defined as those having atomic masses larger than 15 and energies exceeding 1 MeV/amu, may lead to significant modification of the irradiated material in a nanometric region around the (straight) ion trajectory (latent tracks). In the case of amorphous silica, SHI irradiation originates nano-tracks of higher density than the virgin material (densification). As a result, the refractive index is increased with respect to that of the surroundings. Moreover, track overlapping leads to continuous amorphous layers that present a significant contrast with respect to the pristine substrate. We have recently demonstrated that SHI irradiation produces a large number of point defects, easily detectable by a number of experimental techniques (work presented in the parallel conference ICDIM). The mechanisms of energy transfer from SHI to the target material have their origin in the high electronic excitation induced in the solid. A number of phenomenological approaches have been employed to describe these mechanisms: coulomb explosion, thermal spike, non-radiative exciton decay, bond weakening. However, a detailed microscopic description is missing due to the difficulty of modeling the time evolution of the electronic excitation. In this work we have employed molecular dynamics (MD) calculations to determine whether the irradiation effects are related to the thermal phenomena described by MD (in the ps domain) or to electronic phenomena (sub-ps domain), e.g., exciton localization. We have carried out simulations of up to 100 ps with large boxes (30x30x8 nm3) using a home-modified version of MDCASK that allows us to define a central hot cylinder (ion track) from which heat flows to the surrounding cold bath (unirradiated sample). We observed that once the cylinder has cooled down, the Si and O coordination numbers are 4 and 2, respectively, as in virgin silica. On the other hand, the density of the (cold) cylinder increases with respect to that of silica and, furthermore, the silica network ring size decreases. Both effects are in agreement with the observed densification. In conclusion, purely thermal effects do not explain the generation of point defects upon irradiation, but they do account for the silica densification.

Collective Organization of Complex Networks: Emergent Scales, Vulnerability, and Targeting of Desired Dynamics

Cuando una colectividad de sistemas dinámicos acoplados mediante una estructura irregular de interacciones evoluciona, se observan dinámicas de gran complejidad y fenómenos emergentes imposibles de predecir a partir de las propiedades de los sistemas individuales. El objetivo principal de esta tesis es precisamente avanzar en nuestra comprensión de la relación existente entre la topología de interacciones y las dinámicas colectivas que una red compleja es capaz de mantener. Siendo este un tema amplio que se puede abordar desde distintos puntos de vista, en esta tesis se han estudiado tres problemas importantes dentro del mismo que están relacionados entre sí. Por un lado, en numerosos sistemas naturales y artificiales que se pueden describir mediante una red compleja la topología no es estática, sino que depende de la dinámica que se desarrolla en la red: un ejemplo son las redes de neuronas del cerebro. En estas redes adaptativas la propia topología emerge como consecuencia de una autoorganización del sistema. Para conocer mejor cómo pueden emerger espontáneamente las propiedades comúnmente observadas en redes reales, hemos estudiado el comportamiento de sistemas que evolucionan según reglas adaptativas locales con base empírica. Nuestros resultados numéricos y analíticos muestran que la autoorganización del sistema da lugar a dos de las propiedades más universales de las redes complejas: a escala mesoscópica, la aparición de una estructura de comunidades, y, a escala macroscópica, la existencia de una ley de potencias en la distribución de las interacciones en la red. El hecho de que estas propiedades aparecen en dos modelos con leyes de evolución cuantitativamente distintas que siguen unos mismos principios adaptativos sugiere que estamos ante un fenómeno que puede ser muy general, y estar en el origen de estas propiedades en sistemas reales. En segundo lugar, proponemos una medida que permite clasificar los elementos de una red compleja en función de su relevancia para el mantenimiento de dinámicas colectivas. En concreto, estudiamos la vulnerabilidad de los distintos elementos de una red frente a perturbaciones o grandes fluctuaciones, entendida como una medida del impacto que estos acontecimientos externos tienen en la interrupción de una dinámica colectiva. Los resultados que se obtienen indican que la vulnerabilidad dinámica es sobre todo dependiente de propiedades locales, por tanto nuestras conclusiones abarcan diferentes topologías, y muestran la existencia de una dependencia no trivial entre la vulnerabilidad y la conectividad de los elementos de una red. Finalmente, proponemos una estrategia de imposición de una dinámica objetivo genérica en una red dada e investigamos su validez en redes con diversas topologías que mantienen regímenes dinámicos turbulentos. Se obtiene como resultado que las redes heterogéneas (y la amplia mayora de las redes reales estudiadas lo son) son las más adecuadas para nuestra estrategia de targeting de dinámicas deseadas, siendo la estrategia muy efectiva incluso en caso de disponer de un conocimiento muy imperfecto de la topología de la red. Aparte de la relevancia teórica para la comprensión de fenómenos colectivos en sistemas complejos, los métodos y resultados propuestos podrán dar lugar a aplicaciones en sistemas experimentales y tecnológicos, como por ejemplo los sistemas neuronales in vitro, el sistema nervioso central (en el estudio de actividades síncronas de carácter patológico), las redes eléctricas o los sistemas de comunicaciones. ABSTRACT The time evolution of an ensemble of dynamical systems coupled through an irregular interaction scheme gives rise to dynamics of great of complexity and emergent phenomena that cannot be predicted from the properties of the individual systems. The main objective of this thesis is precisely to increase our understanding of the interplay between the interaction topology and the collective dynamics that a complex network can support. This is a very broad subject, so in this thesis we will limit ourselves to the study of three relevant problems that have strong connections among them. First, it is a well-known fact that in many natural and manmade systems that can be represented as complex networks the topology is not static; rather, it depends on the dynamics taking place on the network (as it happens, for instance, in the neuronal networks in the brain). In these adaptive networks the topology itself emerges from the self-organization in the system. To better understand how the properties that are commonly observed in real networks spontaneously emerge, we have studied the behavior of systems that evolve according to local adaptive rules that are empirically motivated. Our numerical and analytical results show that self-organization brings about two of the most universally found properties in complex networks: at the mesoscopic scale, the appearance of a community structure, and, at the macroscopic scale, the existence of a power law in the weight distribution of the network interactions. The fact that these properties show up in two models with quantitatively different mechanisms that follow the same general adaptive principles suggests that our results may be generalized to other systems as well, and they may be behind the origin of these properties in some real systems. We also propose a new measure that provides a ranking of the elements in a network in terms of their relevance for the maintenance of collective dynamics. Specifically, we study the vulnerability of the elements under perturbations or large fluctuations, interpreted as a measure of the impact these external events have on the disruption of collective motion. Our results suggest that the dynamic vulnerability measure depends largely on local properties (our conclusions thus being valid for different topologies) and they show a non-trivial dependence of the vulnerability on the connectivity of the network elements. Finally, we propose a strategy for the imposition of generic goal dynamics on a given network, and we explore its performance in networks with different topologies that support turbulent dynamical regimes. It turns out that heterogeneous networks (and most real networks that have been studied belong in this category) are the most suitable for our strategy for the targeting of desired dynamics, the strategy being very effective even when the knowledge on the network topology is far from accurate. Aside from their theoretical relevance for the understanding of collective phenomena in complex systems, the methods and results here discussed might lead to applications in experimental and technological systems, such as in vitro neuronal systems, the central nervous system (where pathological synchronous activity sometimes occurs), communication systems or power grids.

Local effects of depth-dependent water content of ice and snow and firn layers temperature on a conjectured subglacial lake below Amundsenisen Icefield (Svalbard)

It is known that the Amundsenisen Icefield in Southern Spitzbergen (Svalbard achipelago) is temperate with an upper layer of snow and firn. It is an accumulation area and, though ice/water mass balance is clearly subject to time evolution, observation data on the long-term elevation changes over the past 40 years (Nuth et al., 2010) allow to assume constant icefield surface. Within our study of the plausibility of a subglacial lake (Glowacki et al., 2007), here, we focus on the sensitivity of the system to the thermal effect of the firn and snow layers.

A Special Perturbation Method in Orbital Dynamics

Bead models are used in dynamical simulation of tethers. These models discretize a cable using beads distributed along its length. The time evolution is obtained nu- merically. Typically the number of particles ranges between 5 and 50, depending on the required accuracy. Sometimes the simulation is extended over long periods (several years). The complex interactions between the cable and its spatial environment require to optimize the propagators —both in runtime and precisión that constitute the central core of the process. The special perturbation method treated on this article conjugates simpleness of computer implementation, speediness and precision, and is capable to propagate the orbit of whichever material particle. The paper describes the evolution of some orbital elements, which are constants in a non-perturbed problem, but which evolve in the time scale imposed by the perturbation. It can be used with any kind of orbit and it is free of sin- gularities related to small inclination and/or small eccentricity. The use of Euler parameters makes it robust.

Calculation of dielectrophoretic potential generated under illumination patterns on x- and z-cut LiNbO3 : application to photovoltaic particle trapping

The dielectrophoretic potential generated near the surface of a z-cut LiNbO3 by photovoltaic charge transport has been calculated for first time. The procedure and results are compared with the ones corresponding to x-cut. Diferences in the position, sharpness and time evolution are reported, and their implication on particle trapping are discussed.

A procedure to track vibration modes under changing external factors: application to a pedestrian bridge

A low-cost vibration monitoring system has been developed and installed on an urban steel- plated stress-ribbon footbridge. The system continuously measures: the acceleration (using 18 triaxial MEMS accelerometers distributed along the structure), the ambient temperature and the wind velocity and direction. Automated output-only modal parameter estimation based on the Stochastic Subspace Identification (SSI) is carried out in order to extract the modal parameters, i.e., the natural frequencies, damping ratios and modal shapes. Thus, this paper analyzes the time evolution of the modal parameters over a whole-year data monitoring. Firstly, for similar environmental/operational factors, the uncertainties associated to the time window size used are studied and quantified. Secondly, a methodology to track the vibration modes has been established since several of them with closely-spaced natural frequencies are identified. Thirdly, the modal parameters have been correlated against external factors. It has been shown that this stress-ribbon structure is highly sensitive to temperature variation (frequency changes of more than 20%) with strongly seasonal and daily trends

Spectral diffusion and the energy landscape of a protein

We present a novel type of spectral diffusion experiment in the millikelvin range to characterize the energy landscape of a protein as compared with that of a glass. We measure the time evolution of spectral holes for more than 300 hr after well-defined initial nonequilibrium conditions. We show that the model of noninteracting two-level systems can describe spectral diffusion in the glass, but fails for the protein. Our results further demonstrate that randomness in the energy landscape of a protein shows features of organization. There are “deep minimum” states separated by barriers, the heights of which we are able to estimate. The energy landscape of a glass is featureless by comparison.

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the ≈10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

Photoconductivity and Electro-Optical Properties of Wide Band Gap Metal Oxides

The PhD activity described in this Thesis was focused on the study of metal-oxide wide-bandgap materials, aiming at fabricating new optoelectronic devices such as solar-blind UV photodetectors, high power electronics, and gas sensors. Photocurrent spectroscopy and DC photocurrent time evolution were used to investigate the performance of prototypes under different atmospheres, temperatures and excitation wavelengths (or dark conditions). Cathodoluminescence, absorption spectroscopy, XRD and SEM were used to assess structural, morphologic, electrical and optical properties of materials. This thesis is divided into two main sections, each describing the work done on a different metal-oxide semiconductor. 1) MOVPE-grown Ga2O3 thin films for UV solar-blind photodetectors and high power devices The semiconducting oxides, among them Ga2O3, have been employed for several decades as transparent conducting oxide (TCO) electrodes for fabrication of solar cells, displays, electronic, and opto-electronic devices. The interest was mainly confined to such applications, as these materials tend to grow intrinsically n-type, and attempts to get an effective p-type doping has consistently failed. The key requirements of TCO electrodes are indeed high electrical conductivity and good transparency, while crystallographic perfection is a minor issue. Furthermore, for a long period no high-quality substrates and epi-layers were available, which in turn impeded the development of a truly full-oxide electronics. Recently, Ga2O3 has attracted renewed interest, as large single crystals and high-quality homo- and hetero-epitaxial layers became available, which paved the way to novel application areas. Our research group spent the last two years in developing a low temperature (500-700°C) MOVPE growth procedure to obtain thin films of Ga2O3 on different substrates (Dept. of Physics and IMEM-CNR at UNIPR). We obtained a significant result growing on oriented sapphire epitaxial films of high crystalline, undoped, pure phase -Ga2O3 (hexagonal). The crystallographic properties of this phase were investigated by XRD, in order to clarify the lattice parameters of the hexagonal cell. First design and development of solar blind UV photodetectors based on -phase was carried out and the optoelectronic performance is evaluated by means of photocurrent spectroscopy. The UV-response is adequately fast and reliable to render this unusual phase a subject of great interest for future applications. The availability of a hexagonal phase of Ga2O3 stable up to 700°C, belonging to the same space group of gallium nitride, with high crystallinity and tunable electrical properties, is intriguing in view of the development of nitride-based devices, by taking advantage of the more favorable symmetry and epitaxial relationships with respect to the monoclinic β-phase. In addition, annealing at temperatures higher than 700°C demonstrate that the hexagonal phase converts totally in the monoclinic one. 2) ZnO nano-tetrapods: charge transport mechanisms and time-response in optoelectronic devices and sensors Size and morphology of ZnO at the nanometer scale play a key role in tailoring its physical and chemical properties. Thanks to the possibility of growing zinc oxide in a variety of different nanostructures, there is a great variety of applications, among which gas sensors, light emitting diodes, transparent conducting oxides, solar cells. Even if the operation of ZnO nanostructure-based devices has been recently demonstrated, the mechanisms of charge transport in these assembly is still under debate. The candidate performed an accurate investigation by photocurrent spectroscopy and DC-photocurrent time evolution of electrical response of both single-tetrapod and tetrapod-assembly devices. During the research done for this thesis, a thermal activation energy enables the performance of samples at high temperatures (above about 300°C). The energy barrier is related to the leg-to-leg interconnection in the assembly of nanotetrapods. Percolation mechanisms are responsible for both the very slow photo-response (minutes to hours or days) and the significant persistent photocurrent. Below the bandgap energy, electronic states were investigated but their contribution to the photocurrent are two-three order of magnitude lower than the band edge. Such devices are suitable for employ in photodetectors as well as in gas sensors, provided that the mechanism by which the photo-current is generated and gas adsorption on the surface modify the conductivity of the material are known.

Estabilidade de vórtices em condensados de Bose-Einstein

Neste trabalho de mestrado é estudada a estabilidade de vórtices em condensados de Bose-Einstein com interação atrativa entre os átomos através da solução numérica da equação de Gross-Pitaevskii. Inicialmente são reproduzidos resultados da literatura, nos quais são estudados vórtices em condensados bidimensionais atrativos com potencial interatômico homogêneo em todo o condensado. A estabilidade de tais sistemas é inferida através da solução numérica das equações de Bogoliubov-de Gennes e da evolução temporal dos vórtices. Demonstra-se que esses vórtices são estáveis, até um certo número crítico de átomos, apenas para valores de vorticidade S=1. Em seguida foi proposto um modelo no qual a interação entre os átomos é espacialmente modulada. Neste caso é possível demonstrar que vórtices com valores de vorticidade de até S=6, pelo menos, são estáveis. Finalmente é estudada a estabilidade de vórtices em condensados tridimensionais atrativos, novamente com potencial interatômico homogêneo em todo o condensado. Assim como no caso bidimensional mostra-se que tais vórtices são estáveis para valores de vorticidade de S=1. Espera-se em breve estudar a estabilidade de vórtices em condesados tridimensionais com potencial de interação espacialmente modulado.