Elastic (stress-strain) halo associated to ion-induced nano-tracks in lithium niobate: role of crystal anisotropy

The elastic strain/stress fields (halo) around a compressed amorphous nano-track (core) caused by a single high-energy ion impact on LiNbO3 are calculated. A method is developed to approximately account for the effects of crystal anisotropy of LiNbO3 (symmetry 3m) on the stress fields for tracks oriented along the crystal axes (X, Y or Z). It only considers the zero-order (axial) harmonic contribution to the displacement field in the perpendicular plane and uses effective Poisson moduli for each particular orientation. The anisotropy is relatively small; however, it accounts for some differential features obtained for irradiations along the crystallographic axes X, Y and Z. In particular, the irradiation-induced disorder (including halo) and the associated surface swelling appear to be higher for irradiations along the X- or Y-axis in comparison with those along the Z-axis. Other irradiation effects can be explained by the model, e.g. fracture patterns or the morphology of pores after chemical etching of tracks. Moreover, it offers interesting predictions on the effect of irradiation on lattice parameters

Study of the damage evolution of the concrete under freeze-thaw cycles using traditional and non-traditional techniques

Some experiments have been performed to investigate the cyclic freeze-thaw deterioration of concrete, using traditional and non-traditional techniques. Two concrete mixes, with different pore structure, were tested in order to compare the behavior of a freeze-thaw resistant concrete from one that is not. One of the concretes was air entrained, high content of cement and low w/c ratio, and the other one was a lower cement content and higher w/c ratio, without air-entraining agent. Concrete specimens were studied under cyclic freeze-thaw conditions according to UNE-CENT/TS 12390-9 test, using 3% NaCl solution as freezing medium (CDF test: Capillary Suction, De-icing agent and Freeze-thaw Test). The temperature and relative humidity were measured during the cycles inside the specimens using embedded sensors placed at different heights from the surface in contact with the de-icing agent solution. Strain gauges were used to measure the strain variations at the surface of the specimens. Also, measurements of ultrasonic pulse velocity through the concrete specimens were taken before, during, and after the freeze-thaw cycles. According to the CDF test, the failure of the non-air-entraining agent concrete was observed before 28 freeze-thaw cycles; contrariwise, the scaling of the air-entraining agent concrete was only 0.10 kg/m 2 after 28 cycles, versus 3.23 kg/m 2 in the deteriorated concrete, after 28 cycles. Similar behavior was observed on the strain measurements. The residual strain in the deteriorated concrete after 28 cycles was 1150 m versus 65 m, in the air-entraining agent concrete. By means of monitoring the changes of ultrasonic pulse velocity during the freeze-thaw cycles, the deterioration of the tested specimens were assessed

Elastoplastic consolidation at finite strain. Part 2: finite element implementation and numerical examples

A mathematical model for finite strain elastoplastic consolidation of fully saturated soil media is implemented into a finite element program. The algorithmic treatment of finite strain elastoplasticity for the solid phase is based on multiplicative decomposition and is coupled with the algorithm for fluid flow via the Kirchhoff pore water pressure. A two-field mixed finite element formulation is employed in which the nodal solid displacements and the nodal pore water pressures are coupled via the linear momentum and mass balance equations. The constitutive model for the solid phase is represented by modified Cam—Clay theory formulated in the Kirchhoff principal stress space, and return mapping is carried out in the strain space defined by the invariants of the elastic logarithmic principal stretches. The constitutive model for fluid flow is represented by a generalized Darcy's law formulated with respect to the current configuration. The finite element model is fully amenable to exact linearization. Numerical examples with and without finite deformation effects are presented to demonstrate the impact of geometric nonlinearity on the predicted responses. The paper concludes with an assessment of the performance of the finite element consolidation model with respect to accuracy and numerical stability.

Novel system for bite-force sensing and monitoring based on magnetic eear field communication

Intraoral devices for bite-force sensing have several applications in odontology and maxillofacial surgery, as bite-force measurements provide additional information to help understand the characteristics of bruxism disorders and can also be of help for the evaluation of post-surgical evolution and for comparison of alternative treatments. A new system for measuring human bite forces is proposed in this work. This system has future applications for the monitoring of bruxism events and as a complement for its conventional diagnosis. Bruxism is a pathology consisting of grinding or tight clenching of the upper and lower teeth, which leads to several problems such as lesions to the teeth, headaches, orofacial pain and important disorders of the temporomandibular joint. The prototype uses a magnetic field communication scheme similar to low-frequency radio frequency identification (RFID) technology (NFC). The reader generates a low-frequency magnetic field that is used as the information carrier and powers the sensor. The system is notable because it uses an intra-mouth passive sensor and an external interrogator, which remotely records and processes information regarding a patient?s dental activity. This permits a quantitative assessment of bite-force, without requiring intra-mouth batteries, and can provide supplementary information to polysomnographic recordings, current most adequate early diagnostic method, so as to initiate corrective actions before irreversible dental wear appears. In addition to describing the system?s operational principles and the manufacture of personalized prototypes, this report will also demonstrate the feasibility of the system and results from the first in vitro and in vivo trials.

The evolution of the laws of software evolution. A discussion based on a systematic literature review

After more than 40 years of life, software evolution should be considered as a mature field. However, despite such a long history, many research questions still remain open, and controversial studies about the validity of the laws of software evolution are common. During the first part of these 40 years the laws themselves evolved to adapt to changes in both the research and the software industry environments. This process of adaption to new paradigms, standards, and practices stopped about 15 years ago, when the laws were revised for the last time. However, most controversial studies have been raised during this latter period. Based on a systematic and comprehensive literature review, in this paper we describe how and when the laws, and the software evolution field, evolved. We also address the current state of affairs about the validity of the laws, how they are perceived by the research community, and the developments and challenges that are likely to occur in the coming years.

Resistance to malathion and lambda-cyhalothrin in Ceratitis capitata (Wiedemann) (Diptera: Tephritidae)

La mosca mediterránea de la fruta, Ceratitis capitata (Wiedemann, 1824) (Diptera: Tephritidae), es una de las plagas de mayor incidencia económica en cítricos y otros frutales a nivel mundial. En España las medidas de control de esta plaga en cítricos, desde mediados de los 90 hasta 2009, se basaron principalmente en el monitoreo de las poblaciones y en la aplicación de tratamientos aéreos y terrestres con malatión cebo. Sin embargo, desde la retirada en la Unión Europea en 2009 de los productos fitosanitarios que contienen malatión, los insecticidas más utilizados para el control de esta plaga han sido lambda-cihalotrina y spinosad. En 2004-2005 se detectaron poblaciones españolas de C. capitata resistentes a malatión. Esta resistencia se ha asociado a una mutación (G328A) en la acetilcolinesterasa (AChE), a una duplicación del gen de la AChE (Ccace2) (una de las copias lleva la mutación G328A) y a resistencia metabólica mediada por esterasas (posiblemente aliesterasas). Sin embargo, cuando se secuenció la aliesterasa CcE7 en individuos de una línea resistente a malatión, no se encontró ninguna de las mutaciones (G137D y/o W251L/S/G) asociadas a resistencia en otras especies, si bien se encontraron otras mutaciones al compararlos con individuos de una línea susceptible. Asimismo, mediante la selección en laboratorio de una línea resistente a malatión (W-4Km) con lambda-cihalotrina, se ha podido obtener una línea resistente a lambda-cihalotrina (W-1K). Finalmente, se ha demostrado la capacidad de esta especie para desarrollar resistencia a spinosad mediante selección en laboratorio. Los múltiples mecanismos de resistencia identificados evidencian el potencial de esta especie para desarrollar resistencia a insecticidas con diferentes modos de acción. Los objetivos de esta tesis doctoral son: 1) evaluar la susceptibilidad de poblaciones españolas de campo de C. capitata a lambda-cihalotrina y dilucidar los mecanismos de resistencia en la línea W-1Kλ; 2) comparar la herencia, el coste biológico y la estabilidad de la resistencia a malatión mediada por la mutación G328A y la duplicación del gen Ccace2 (una de las copias lleva la mutación G328A); y 3) investigar el papel de las mutaciones identificadas en la aliesterasa CcαE7 en la resistencia a malatión. Estos estudios son de utilidad para el desarrollo de estrategias de manejo de la resistencia que puedan prevenir o retrasar la aparición de resistencia y aumentar la sostenibilidad de los insecticidas disponibles para el control de esta plaga. Nuestros resultados indican que las poblaciones españolas de C. capitata analizadas han desarrollado resistencia a lambda-cihalotrina. Los valores de CL50 estimados para las poblaciones recogidas en la Comunidad Valenciana, Cataluña y Andalucía oscilaron entre 129 ppm y 287 ppm, igualando o sobrepasando la concentración recomendada para los tratamientos de campo (125 ppm). Estos resultados contrastan con los obtenidos con tres poblaciones de campo recogidas en Túnez, cuya susceptibilidad fue similar a la de la línea control (C). La línea resistente a lambda-cihalotrina W-1K se continuó seleccionando en el laboratorio alcanzándose unos niveles de resistencia de 205 veces con respecto a la línea C, siendo su CL50 (4224 ppm) más de 30 veces superior a la concentración recomendada para los tratamientos de campo. Esta línea resistente mostró altos niveles de resistencia cruzada a deltametrina (150 veces) y a etofenprox (240 veces), lo que sugiere que el desarrollo de resistencia a lambda-cihalotrina podría comprometer la eficacia de otros piretroides para el control de esta plaga. Hemos demostrado que la resistencia de la línea W-1K a lambda-cihalotrina fue casi completamente suprimida por el sinergista PBO, lo que indica que las enzimas P450 desempeñan un papel muy importante en la resistencia a este insecticida. Sin embargo, tanto las moscas de la línea susceptible C como las de la línea resistente W-1K perdieron inmediatamente la capacidad de caminar (efecto “knock-down”) al ser tratadas tópicamente con lambda-cihalotrina, lo que sugiere que la resistencia no está mediada por alteraciones en la molécula diana (resistencia tipo “kdr”). La resistencia metabólica mediada por P450 fue analizada comparando la expresión de 53 genes CYP (codifican enzimas P450) de las familias CYP4, CYP6, CYP9 y CYP12 en adultos de la línea resistente W-1K y de la línea susceptible C. Nuestros resultados muestran que el gen CYP6A51 (número de acceso GenBank XM_004534804) fue sobreexpresado (13-18 veces) en la línea W-1K. Por otra parte, la expresión del gen CYP6A51 fue inducida tanto en adultos de la línea W-1K como de la línea C al ser tratados con lambda-cihalotrina. Sin embargo, no se obtuvieron diferencias significativas entre la línea susceptible C y la línea resistente W-1K al comparar la cantidad de P450 y la actividad NADPH-citocromo c reductasa presente en fracciones microsomales obtenidas a partir de abdómenes. Asimismo, no hemos podido correlacionar el metabolismo de deltametrina, estimado in vitro mediante la incubación de este insecticida con fracciones microsomales, con el nivel de resistencia a este piretroide observado en los bioensayos con la línea W-1K. Por otro lado, no se encontró ninguna alteración en la región promotora 5'UTR del gen CYP6A51 (-500 pb desde el inicio de la traducción) que pudiera explicar su sobreexpresión en la línea W-1K. Los datos obtenidos sugieren que la resistencia a lambda-cihalotrina en la línea W-1K está mediada por P450 y que la sobreexpresión de CYP6A51 puede desempeñar un papel importante, aunque se necesitan más evidencias para establecer una asociación directa de la resistencia con este gen. Hemos estudiado la herencia, el coste biológico y la estabilidad de la resistencia a malatión mediada por la mutación G328A y la duplicación del gen Ccace2 (una de las copias lleva la mutación G328A). La línea susceptible C, donde no se encuentra la mutación G328A (genotipo S/S), se cruzó con dos isolíneas establecidas para representar genotipos únicos correspondientes a los dos mecanismos de resistencia asociados a la molécula diana: 1) la isolínea 267Y (genotipo R/R) establecida a partir de una pareja que portaba la mutación G328A en homocigosis; 2) la isolínea 306TY (genotipo RS/RS) establecida a partir de una pareja que portaba en homocigosis la duplicación del gen Ccace2. No se realizaron cruces recíprocos, ya que mediante experimentos de hibridación in situ en cromosomas politénicos se pudo comprobar que el locus de la AChE y la duplicación (probablemente en tándem) se localizan en el cromosoma autosómico 2L. La susceptibilidad al malatión de los parentales resistentes (R/R o RS/RS) y susceptibles (S/S), los cruces F1 (S/R, S/RS y R/RS) y los retrocruzamientos indican que la resistencia a malatión es semi-dominante en ambos casos. Sin embargo, nuestros resultados no fueron concluyentes con respecto a la naturaleza monogénica de la resistencia a malatión en estas isolíneas. Por lo tanto, no podemos descartar que otros genes que contribuyan a la resistencia, además de la mutación G328A (isolínea 267Y) y de la duplicación del gen Ccace2 (isolínea 306TY), puedan haber sido seleccionados durante el proceso de selección de 267Y y 306TY. Varios parámetros biológicos fueron evaluados para determinar si estos dos mecanismos de resistencia a malatión suponen un coste biológico para los genotipos resistentes. Individuos con genotipo R/R mostraron un retraso en el tiempo de desarrollo de huevo a pupa, un peso de pupa reducido y una menor longevidad de los adultos, en comparación con los individuos con genotipo S/S. Sin embargo, el peso de pupa de los individuos con genotipo RS/RS fue similar al de los individuos S/S, y su desarrollo de huevo a pupa intermedio entre S/S y R/R. Estas diferencias en el coste biológico pueden estar relacionadas con la reducción de la eficiencia catalítica de la AChE mutada en los individuos R/R, y al efecto compensatorio que la copia no mutada del gen tiene en los individuos RS/RS que portan la duplicación. La estabilidad de la resistencia a malatión mediada por la mutación G328A y la duplicación se analizó mediante el seguimiento de los caracteres de resistencia en la progenie de retrocruzamientos S/R x R/R y S/RS x RS/RS a lo largo de varias generaciones en ausencia de presión de selección con insecticidas. Nuestros resultados muestran que la frecuencia del alelo que porta la mutación G328A disminuyó desde 67,5% en la primera generación del retrocruzamiento S/R x R/R (75% esperado, asumiendo segregación mendeliana y que sólo hay dos alelos: uno mutado y otro no mutado) a 12% después de 10 generaciones. Por el contrario, la frecuencia de la duplicación sólo disminuyó desde 75% en en la primera generación del retrocruzamiento S/RS x RS/RS (75% esperado, asumiendo segregación Mendeliana y que la duplicación segrega como un único alelo) a 50% en el mismo período, lo que indica que la duplicación es más estable que la mutación. Asimismo, se analizó la presencia de la mutación y de la duplicación en poblaciones de campo recogidas en seis localidades en 2004-2007, cuando todavía se usaba el malatión, y se comparó con poblaciones recogidas en los mismos campos en 2010, un año después de la prohibición del malatión en la Unión Europea. La frecuencia media del genotipo susceptible (S/S) aumentó del 55,9% en el período 2004-2007 a 70,8% en 2010, mientras que la frecuencia de los genotipos portadores de la mutación en homocigosis o heterocigosis (R/R y S/R) disminuyó del 30,4 al 9,2%, los que llevan la duplicación en homocigosis o heterocigosis (RS/RS y S/RS) aumentaron levemente desde 12,8 hasta 13,3%, y los que llevan a la vez la mutación y la duplicación (R/RS) también aumentaron del 1 al 6,7%. Estos resultados son consistentes con que la duplicación del gen Ccace2 (con una copia con la mutación G328A y la otra copia no mutada) es más ventajosa que la mutación G328A por si sola, ya que la duplicación mantiene los niveles de resistencia a la vez que limita el coste biológico. Para investigar la asociación entre la resistencia a malatión y las mutaciones encontradas previamente en CcE7, hemos generado isolíneas con mutaciones específicas seleccionadas por su ubicación próxima a la entrada al centro activo de la enzima. La isolínea Sm2 (procedente de una hembra heterocigota para la mutación V96L y un macho homocigoto para el alelo no mutado) mantuvo altos niveles de resistencia a malatión, incluso después de 30 generaciones sin presión de selección. Por el contrario, la isolínea 267Y (compuesta por individuos homocigotos para la mutación L267Y) y la línea 306TY (compuesta por individuos homocigotos para la doble mutación R306T-N307Y) mostraron una reducción significativa en los niveles de resistencia. También hemos encontrado que la resistencia a malatión de la línea Sm2 fue parcialmente revertida por DEF y TPP, y que Sm2 mostró una reducción significativa en la actividad MTB, como se ha descrito en otras especies que muestran resistencia específica a malatión mediada por aliesterases. Además, fue posible asociar la presencia de la mutación V96L en individuos de la línea Sm2 con supervivencia a una concentración discriminante de malatión (5,000 ppm) y con una baja actividad MTB. Estos resultados sugieren una posible relación entre la mutación V96L en la aliesterasa CcE7 y la resistencia a malatión, aunque todavía no se puede concluir que la resistencia es causada por esta mutación, siendo necesarios más estudios para comprobar su contribución a la resistencia. En conclusión, se ha encontrado por primera vez resistencia a lambda-cihalotrina en poblaciones de campo de C. capitata, y nuestros resultados indican que las P450 son el principal mecanismo de resistencia en la línea W-1K. Esta situación se suma al caso previamente descrito de resistencia en campo a malatión asociada a la mutación G328A, a la duplicación del gen Ccace2 (una de las copias lleva la mutación G328A) y a resistencia metabólica mediada por esterasas. Nuestros resultados también indican que la alteración de la molécula diana AChE parece ser responsable de un cierto nivel de resistencia a malatión en C. capitata, que puede ser estimada como aproximadamente 25-40 veces para la mutación G328A y 40-60 veces para la duplicación; mientras que la resistencia mediada por esterasas y que ha sido asociada en este estudio con la mutación V96L en CcE7 puede conferir un efecto multiplicativo (por un factor de 5 a 10) aumentando la resistencia a malatión a 200-400 veces. Por otra parte, hemos demostrado que los insectos resistentes que llevan la duplicación tienen un coste biológico menor y muestran una estabilidad mayor que aquellos con la mutación G328A en ausencia de presión de selección con insecticidas. Esto representa un escenario en el que los genotipos con la duplicación permanecerán en el campo en frecuencias bajas a moderadas, pero podrían ser seleccionados rápidamente si se utilizan malatión u otros insecticidas que muestren resistencia cruzada. Estos resultados tienen importantes implicaciones para los programas de manejo de la resistencia, ya que el repertorio de insecticidas eficaces para el control de C. capitata es cada vez más limitado. Además, la coexistencia de múltiples mecanismos de resistencia en poblaciones de campo ofrece el potencial para desarrollar resistencia frente a otros insecticidas disponibles para el control de esta plaga. Estrategias para de manejo de la resistencia basadas en la alternancia de insecticidas con diferentes modos de acción, y su combinación con otros métodos de control, deben ser implementadas para evitar el desarrollo de resistencia en campo. ABSTRACT The Mediterranean fruit fly (Medfly), Ceratitis capitata (Wiedemann, 1824) (Diptera: Tephritidae), is one of the most economically damaging pests of citrus and other fruit crops worldwide. Control measures in citrus crops in Spain from the mid 90's to 2009 were mainly based on field monitoring of population levels and aerial and ground treatments with malathion bait sprays. However, since the withdrawal of phytosanitary products containing malathion in the European Union in 2009, lambda-cyhalothrin and spinosad have become the most widely used insecticides for the control of this pest. Resistance to malathion was found in Spanish field populations of C. capitata in 2004-2005. This resistance has been associated with a mutation G328A in the acetylcholinesterase (AChE), a duplication of the AChE gene (Ccace2) (one of the copies bearing the mutation G328A), and metabolic resistance mediated by esterases (probably aliesterases). However, when the gene of the aliesterase CcE7 was sequenced in individuals from a malathion resistant strain of C. capitata, none of the known G137D and/or W251L/S/G mutations associated to resistance in other species were found, though other mutations were detected when compared with individuals from a susceptible strain. Noteworthy, a lambda-cyhalothrin resistant strain (W-1K) was obtained by selecting a field-derived malathion resistant strain (W-4Km) with lambda-cyhalothrin. Moreover, it has also been demonstrated the capacity of this species to develop resistance to spinosad by laboratory selection. The multiple resistance mechanisms identified highlight the potential of this species to develop resistance to insecticides with different modes of action. The objectives of this PhD Thesis are: 1) to assess the susceptibility of Spanish field populations of C. capitata to lambda-cyhalothrin and to elucidate the resistance mechanisms in the W-1Kλ strain; 2) to compare the inheritance, fitness cost and stability of the malathion resistance mediated by the G328A mutation and the duplication of the Ccace2 gene (with one of the copies bearing the mutation G328A); and 3) to investigate the role of the aliesterase CcαE7 mutations in malathion resistance. All these studies will be of use for devising proactive resistance management strategies that could prevent or delay resistance development and would increase the sustainability of the insecticides available for Medfly control. Our results indicate that Spanish field populations of C. capitata have developed resistance to lambda-cyhalothrin. The LC50 values estimated for populations collected at Comunidad Valenciana, Cataluña and Andalucía ranged from 129 ppm to 287 ppm, equaling or overpassing the recommended concentration for field treatments (125 ppm). These results contrast with those obtained with three different Tunisian field populations, whose susceptibility was similar to that of the control (C) strain. The lambda-cyhalothrin resistant W-1K strain has been further selected to achieve a 205-fold resistance compared to the C strain, being its LC50 (4,224 ppm) more than 30 times higher than the recommended concentration for field applications. This resistant strain showed high levels of cross-resistance to deltamethrin (150-fold) and etofenprox (240-fold), suggesting that the development of resistance to lambda-cyhalothrin may compromise the effectiveness of other pyrethroids for the control of this species. We have shown that the resistance of the W-1K strain to lambda-cyhalothrin was almost completely suppressed by the synergist PBO, indicating that P450 enzymes play a very important role in resistance to this insecticide. However, both susceptible C and resistant W-1K flies were knocked down after topical treatment with lambda-cyhalothrin, suggesting that kdr resistance mediated by alterations of the target site is not playing a major role. Metabolic resistance mediated by P450 was further analyzed by comparing the expression of 53 genes of the families CYP4, CYP6, CYP9 and CYP12 in adults flies from the resistant W-1K and the susceptible C strains. We found that the gene CYP6A51 (GenBank accession number XM_004534804) was overexpressed (13-18-fold) in the W-1K strain. Moreover, the expression of the CYP6A51 gene was induced when adults of the W-1K and C strains were treated with lambda-cyhalothrin. However, no significant differences were obtained between susceptible C and resistant W-1K strains for the quantity of P450 and for the activity of NADPH- cytochrome c reductase measured in microsomal fractions obtained from abdomens. Moreover, we failed to correlate the metabolism of deltamethrin, analyzed in vitro by incubating this insecticide with microsomal fractions, with the resistance level against this pyrethroid observed in bioassays with W-1K. The sequencing of the 5´UTR region of the CYP6A51 gene failed in finding an alteration in the promoter region (-500 bp from translation start site) that could explain overexpression in the W-1K strain. All data obtained suggest that resistance to lambda-cyhalothrin in the W- 1K strain is mediated by P450 and that overexpression of CYP6A51 may play a major role, although further evidences are needed to establish a direct association of resistance with this gene. We have studied the inheritance, fitness cost and stability of the malathion resistance mediated by the G328A mutation and the duplication of the Ccace2 gene (with one of the copies bearing the mutation G328A). The malathion-susceptible C strain where the G328A mutation is not found (S/S genotype) was crossed with two isolines established to represent unique genotypes corresponding to the two target-site resistance mechanisms: 1) the 267Y isoline (genotype R/R) was established from a couple bearing the mutation G328A in homozygosis; and 2) the 306TY isoline (genotype RS/RS) was established from a couple being homozygous for the duplication of the Ccace2 gene. Reciprocal crosses have not been performed, since in situ hybridization on polythene chromosomes showed that the AChE locus and the duplication (most probably in tandem) are placed at the autosomal chromosome 2L. Mortality responses to malathion of resistant isolines (R/R or RS/RS) and susceptible (S/S) genotypes, F1 crosses (S/R, S/RS, and R/RS), and the back-crosses indicated that resistance to malathion is inherited as a semi-dominant trait in both cases. However, our results were not conclusive about the monogenic nature of the resistance to malathion in these isolines. Thus, we can not discard that other genes contributing to resistance, in addition to the mutation G328A (isoline 267Y) and the duplication of the Ccace2 gene (isoline 306TY), may have been selected during the selection process of 267Y and 306TY. Several biological parameters were evaluated to determine if these two malathion resistance mechanisms impose a fitness cost for resistant genotypes. Individuals with genotype R/R have a reduced fitness in terms of developmental time from egg to pupa, pupal weight and adult longevity, when compared to susceptible individuals (genotype S/S). Interestingly, the fitness cost was substantially diminished in individuals with genotype RS/RS. These differences in fitness may be related to the reduction of the catalytic efficiency of mutated AChE in individuals R/R, and the compensatory effect that the non-mutated copy of the gene has on individuals RS/RS bearing the duplication. The stability of malathion reistance associated with the mutation G328A or the duplication was analyzed by following these resistant traits in the progeny of the back-crosses S/RS x RS/RS and S/R x R/R over consecutive generations in the absence of insecticide selection pressure. Our results show that the frequency of the allele bearing the mutation G328A decreased from 67.5% at the first generation of the back-cross S/R x R/R (75% expected, assuming Mendelian segregation and that there are only two alleles: one mutated and the other non-mutated) to 12% after 10 generations. By contrast, the frequency of the duplication only declined from 75% at the first generation of the back-cross S/RS x RS/RS (75% expected, assuming Mendelian segregation and that the duplication segregates as an unique allele) to 50% in the same period, indicating that the duplication is more stable than the mutation. The presence of the mutation and the duplication was analyzed in field populations collected in six localities in 2004-2007, when malathion was still used, and compared to populations collected in the same fields in 2010, one year after the prohibition of malathion in the European Union. The average frequency of the susceptible genotype (S/S) increased from 55.9% in the period 2004-2007 to 70.8% in 2010, whereas the frequency of those genotypes carrying the mutation in homozygosis or heterozygosis (R/R and S/R) declined from 30.4 to 9.2%, those carrying the duplication in homozygosis or heterozygosis (RS/RS and S/RS) increased slightly from 12.8 to 13.3%, and those carrying both the mutation and the duplication (R/RS) also increased from 1 to 6.7%. These results are consistent with the duplication of the Ccace2 gene (with one of the copies bearing the mutation G328A and the other copy non-mutated) being more advantageous than the G328A mutation alone by maintaining resistance while restoring part of the fitness. In order to investigate the association of malathion resistance with mutations previously found in the aliesterase CcE7, we have generated isolines bearing specific mutations selected by their putative location near the upper part of the active site gorge of the enzyme. The isoline Sm2 (originating from a female heterozygous for the mutation V96L and a male homozygous for the non-mutated allele) kept high levels of resistance to malathion, even after 30 generations without selection pressure. On the contrary, the isoline 267Y (composed by individuals homozygous for the mutation L267Y) and the strain 306TY (composed by homozygous for the double mutation R306T-N307Y) showed a significant reduction in the levels of resistance. We have found also that resistance to malathion in the Sm2 isoline was partially reverted by DEF and TPP, and that Sm2 showed a significant reduction in MTB activity, as reported for other species showing malathion-specific resistance mediated by aliesterases. Besides, it was possible to associate the presence of the mutation V96L in individuals from the Sm2 isoline with both survival to a discriminating concentration of malathion (5,000 ppm) and low MTB activity. Our results point out to a possible connection betwen the mutation V96L in the aliesterase CcE7 and resistance to malathion, though we can not yet conclude that the resistance is caused by the mutation, being needed further work to understand its contribution to resistance. In conclusion, resistance to lambda-cyhalothrin has been found for the first time in field populations of C. capitata, and metabolic resistance mediated by P450 appears to be the main resistance mechanism in the resistant strain W-1K. These findings add to the previously reported case of field resistance to malathion, associated to the G328A mutation and the duplication of the Ccace2 gene (with one of the copies bearing the mutation G328A) and to metabolic resistance mediated by esterases. Our results also indicate that altered target site AChE appears to be responsible for a certain level of resistance to malathion in C. capitata, that can be estimated as about 25-40-fold for the mutation G328A and 40-60-fold for the duplication; whereas metabolic resistance mediated by esterases and associated in this study with the mutation V96L in CcE7 may confer a multiplicative effect (by a factor of 5 to10) increasing malathion resistance to 200-400-fold. Moreover, we have shown that resistant insects carrying the duplication have better fitness and exhibit a higher stability than those with the mutation G328A in the absence of insecticide pressure. This represents a scenario where genotypes with the duplication will remain in the field at low to moderate frequencies, but could be rapidly selected if malathion or other insecticides showing cross-resistance are used. These findings have important implications for resistance management programs, as the repertoire of effective insecticides for C. capitata control is becoming very limited. Besides, multiple resistance mechanisms coexisting in field populations provide the potential to develop resistance to other available insecticides for the control of this pest. Appropriate resistance management strategies based on the alternation of insecticides with different modes of action, and their combination with other control methods, must then be implemented to avoid the evolution of resistance in the field.

Evolution of microstructure, macrotexture and mechanical properties of commercially pure Ti during ECAP-conform processing and drawing

Long-length ultrafine-grained (UFG) Ti rods are produced by equal-channel angular pressing via the conform scheme (ECAP-C) at 200 °C, which is followed by drawing at 200 °C. The evolution of microstructure, macrotexture, and mechanical properties (yield strength, ultimate tensile strength, failure stress, uniform elongation, elongation to failure) of pure Ti during this thermo-mechanical processing is studied. Special attention is also paid to the effect of microstructure on the mechanical behavior of the material after macrolocalization of plastic flow. The number of ECAP-C passes varies in the range of 1–10. The microstructure is more refined with increasing number of ECAP-C passes. Formation of homogeneous microstructure with a grain/subgrain size of 200 nm and its saturation after 6 ECAP-C passes are observed. Strength properties increase with increasing number of ECAP passes and saturate after 6 ECAP-C passes to a yield strength of 973 MPa, an ultimate tensile strength of 1035 MPa, and a true failure stress of 1400 MPa (from 625, 750, and 1150 MPa in the as-received condition). The true strain at failure failure decreases after ECAP-C processing. The reduction of area and true strain to failure values do not decrease after ECAP-C processing. The sample after 6 ECAP-C passes is subjected to drawing at 200¯C resulting in reduction of a grain/subgrain size to 150 nm, formation of (10 1¯0) fiber texture with respect to the rod axis, and further increase of the yield strength up to 1190 MPa, the ultimate tensile strength up to 1230 MPa and the true failure stress up to 1600 MPa. It is demonstrated that UFG CP Ti has low resistance to macrolocalization of plastic deformation and high resistance to crack formation after necking.

Mechanical behavior and residual stress field of dissimilar laser welded joints

The present investigation addresses the mechanical behavior and residual stress field of dissimilar joints produced by laser welding. Microstructure characterization and residual strain scanning, carried out by neutron diffraction, were used to assess the joints features. It was found that the heat source position influences the base metals dilution and the residual stress field associated to the welding process. The tensile behavior of the joint, different zones achieved by using a video-image based system (VIC-2D) reveals that the residual stress field, together with the positive difference in yield between the weld metal and the base materials protects the joint from being plastically deformed.

Geomorphological analysis and evolution of an altered floodplain drainage system. The case of the Partido Stream (Spain)

The Partido Stream is a small torrential course that flows into the marsh of the Doñana National Park, an area that was declared a World Heritage Site in 1994. Before 1981, floods occurred, and the stream overflowed onto a floodplain. As an old alluvial fan, the floodplain has its singular orography and functionality. Fromthe floodplain, several drainage channels, locally called caño, discharged into themarsh. The Partido Streamhad themorphology of a caño and covered approximately 8 km from the old fan to the marsh. The stream was straightened and channelised in 1981 to cultivate the old fan. This resulted in floods that were concentrated between the banks in the following years, which caused the depth of water and the shear stress to increase, thus, scouring the river bed and river banks. In this case, the eroded materials were carried towards the marsh where a new alluvial fan evolved. Control measures on the old fan were implemented in 2006 to stop the development of the new alluvial fan downstream over the marsh. Thus, the stream would partially recover its original behaviour that it had before channelisation, moving forwards in a new, balanced state. The present study describes the geomorphological evolution that channelisation has caused since 1981 and the later slow process of recovery of the original hydraulic-sedimentation regime since 2006. Additionally, it deepens the understanding of the original hydraulic behaviour of the stream, combining field data and 2D simulations.

Caracterización dinámica del glaciar Hurd combinando observaciones de campo y simulaciones numéricas

El objetivo fundamental de esta tesis es la caracterización de la morfología y del estado de deformaciones y tensiones del Glaciar Hurd (Isla Livingston, Archipiélago de las Shetland del Sur, Antártida), mediante una combinación de observaciones de campo, registros de georradar y simulaciones numéricas. La morfología y el estado de deformaciones y tensiones actuales son la expresión de la evolución dinámica del glaciar desde tiempos pretéritos hasta recientes, y su análisis nos dará las pautas con las cuales ser capaces de predecir, con el apoyo de las simulaciones numéricas, su evolución futura. El primer aspecto que se aborda es el estudio de las estructuras que pueden observarse en la superficie del glaciar. Describimos las distintas técnicas utilizadas (medidas de campo, fotointerpretación de ortofotografías, análisis geoquímico de cenizas volcánicas, etc.) y presentamos el análisis e interpretación de los resultados morfo-estructurales, así como la correlación, mediante análisis geoquímicos (fluorescencia de rayos X), entre las cenizas volcánicas que extruyen en la superficie del Glaciar Hurd y las del volcán Decepción, origen de las cenizas. Esto nos permite realizar una datación de las mismas como Tefra 1, correspondiente a la erupción de 1970, Tefra 2, correspondiente a las erupciones pre-1829, y el conjunto Tefra 3, asociado a las erupciones más antiguas. En segundo lugar nos ocupamos de las estructuras presentes en el interior del glaciar, cuya herramienta de detección fundamental es el georradar. Identificadas estas estructuras internas, las vinculamos con las observadas en la superficie del glaciar. También hemos estudiado la estructura hidrotérmica del glaciar, obteniendo una serie de evidencias adicionales de su carácter politérmico. Entre éstas se contaban, hasta ahora, las basadas en el valor del parámetro de rigidez de la relación constitutiva del hielo determinada por ajuste de modelos dinámicos y observaciones realizados por Otero (2008) y las basadas en las velocidades de las ondas de radar en el hielo determinadas con el método de punto medio común por Navarro y otros (2009). Las evidencias adicionales que aportamos en esta tesis son: 1) la presencia de estructuras típicas de régimen compresivo en la zona terminal del glaciar y de cizalla en los márgenes del mismo, y 2) la presencia de un estrato superficial de hielo frío (por encima de otro templado) en la zona de ablación de los tres lóbulos del Glaciar Hurd –Sally Rocks, Argentina y Las Palmas–, que alcanzan espesores de 70, 50 y 40 m, respectivamente. Este estrato de hielo frío está probablemente congelado al lecho subglaciar en la zona terminal (Molina y otros, 2007; esta tesis). Por último, nos ocupamos de la simulación numérica de la dinámica glaciar. Presentamos el modelo físico-matemático utilizado, discutimos sus condiciones de contorno y cómo éstas se miden en los trabajos de campo, y describimos el procedimiento de resolución numérica del sistema de ecuaciones parciales del modelo. Presentamos los resultados para los campos de velocidades, deformaciones y tensiones, comparando estos resultados con las estructuras observadas. También incluimos el análisis de las elipses de deformación acumulativa, que proporcionan información sobre las estructuras a las que puede dar lugar la evolución del estado de deformaciones y tensiones a las que se ve sometido el hielo según avanza, lentamente, desde la cabecera hasta la zona terminal del glaciar, con tiempos de tránsito de hasta 1.250 años, recogiendo así la historia de deformaciones en el glaciar. Concluyendo, ponemos de manifiesto en esta tesis que las medidas de campo de las estructuras y niveles de cenizas, las medidas de georradar y las simulaciones numéricas de la dinámica glaciar, realizadas de forma combinada, permiten caracterizar el régimen actual de velocidades, deformaciones y tensiones del glaciar, entender su evolución en el pasado y predecir su evolución futura. ABSTRACT The main objective of this thesis is to characterize the morphology and the state of strains and stresses of Hurd Glacier (Livingston Island, South Shetland Islands archipelago, Antarctica) through a combination of field observations, ground-penetrating radar measurements and numerical simulations. The morphology and the current state of strain and stresses are the expression of the dynamic evolution of the glacier from the past to recent times, and their analysis gives us the guidelines to be able to predict, with the support of numerical simulations, its future evolution. The first subject addressed is the study of structures that can be observed on the glacier surface. We describe the different techniques used (field measurements, photointerpretation of orthophotos, geochemical analysis of volcanic ashes, etc.) and we present the analysis and interpretation of the morpho-structural results, as well as the correlation with geochemical analysis (XRF) between the volcanic ashes extruded to the surface of Hurd Glacier and those of Deception Island volcano, from which the ashes originate. This allows us dating the ashes as Tephra 1, corresponding to the 1970 eruption, Tephra 2, corresponding to the pre-1829 eruptions, and the Tephra 3 group, associated with older eruptions. Secondly we focus on the study of the structures present within the glacier, which are detected with the help of ground-penetrating radar. Once identified, we link these internal structures with those observed on the glacier surface. We also study the hydrothermal structure of the glacier, getting a series of additional evidences of its polythermal structure. Among the evidences available so far, we can mention those based on the value of the stiffness parameter of the constitutive relation of ice, determined by fitting dynamic models to observations, as done by Otero (2008), and those based on the velocity of propagation of the radar waves through the glacier ice, measured using the common midpoint method, as done by Navarro et al. (2009). The additional evidences that we provide in this thesis are: 1) the presence of structures typical of compressive regime in the terminal zone of the glacier, together with shear at its margins, and 2) the presence of a surface layer of cold ice (overlying a layer of temperate ice) in the ablation zone of the three lobes of Hurd Glacier –Sally Rocks, Argentina and Las Palmas–, reaching thicknesses of 70, 50 and 40 m, respectively. This cold layer is probably frozen to the subglacial bed in the terminal zone (Molina and others 2007; this thesis). Finally, we deal with the numerical simulation of glacier dynamics. We present the physical-mathematical model, discuss its boundary conditions and how they are measured in the field work, and describe the method of numerical solution of the model’s partial differential equations. We present the results for the velocity, strain and stress fields, comparing these results with the observed structures. We also include an analysis of the ellipses of cumulative deformation, which provide information about the structures that can result from the evolution of the strain and stress regime of the glacier ice as it moves slowly from the head to the snout of the glacier, with transit times of up to 1,250 years, so picking the history of deformation of the glacier. Summarizing, we show in this thesis that field measurements of structures and ash layers, ground-penetrating radar measurements and numerical simulations of glacier dynamics, performed in combination, allow us to characterize the current regime of velocities, strains and stresses of the glacier, to understand its past evolution and to predict its future evolution.

Application of PZT and FBG impedance sensors for structural health monitoring of reinforced concrete beams strengthened with FRP composite materials

Esta Tesis tiene como objetivo principal el desarrollo de métodos de identificación del daño que sean robustos y fiables, enfocados a sistemas estructurales experimentales, fundamentalmente a las estructuras de hormigón armado reforzadas externamente con bandas fibras de polímeros reforzados (FRP). El modo de fallo de este tipo de sistema estructural es crítico, pues generalmente es debido a un despegue repentino y frágil de la banda del refuerzo FRP originado en grietas intermedias causadas por la flexión. La detección de este despegue en su fase inicial es fundamental para prevenir fallos futuros, que pueden ser catastróficos. Inicialmente, se lleva a cabo una revisión del método de la Impedancia Electro-Mecánica (EMI), de cara a exponer sus capacidades para la detección de daño. Una vez la tecnología apropiada es seleccionada, lo que incluye un analizador de impedancias así como novedosos sensores PZT para monitorización inteligente, se ha diseñado un procedimiento automático basado en los registros de impedancias de distintas estructuras de laboratorio. Basándonos en el hecho de que las mediciones de impedancias son posibles gracias a una colocación adecuada de una red de sensores PZT, la estimación de la presencia de daño se realiza analizando los resultados de distintos indicadores de daño obtenidos de la literatura. Para que este proceso sea automático y que no sean necesarios conocimientos previos sobre el método EMI para realizar un experimento, se ha diseñado e implementado un Interfaz Gráfico de Usuario, transformando la medición de impedancias en un proceso fácil e intuitivo. Se evalúa entonces el daño a través de los correspondientes índices de daño, intentando estimar no sólo su severidad, sino también su localización aproximada. El desarrollo de estos experimentos en cualquier estructura genera grandes cantidades de datos que han de ser procesados, y algunas veces los índices de daño no son suficientes para una evaluación completa de la integridad de una estructura. En la mayoría de los casos se pueden encontrar patrones de daño en los datos, pero no se tiene información a priori del estado de la estructura. En este punto, se ha hecho una importante investigación en técnicas de reconocimiento de patrones particularmente en aprendizaje no supervisado, encontrando aplicaciones interesantes en el campo de la medicina. De ahí surge una idea creativa e innovadora: detectar y seguir la evolución del daño en distintas estructuras como si se tratase de un cáncer propagándose por el cuerpo humano. En ese sentido, las lecturas de impedancias se emplean como información intrínseca de la salud de la propia estructura, de forma que se pueden aplicar las mismas técnicas que las empleadas en la investigación del cáncer. En este caso, se ha aplicado un algoritmo de clasificación jerárquica dado que ilustra además la clasificación de los datos de forma gráfica, incluyendo información cualitativa y cuantitativa sobre el daño. Se ha investigado la efectividad de este procedimiento a través de tres estructuras de laboratorio, como son una viga de aluminio, una unión atornillada de aluminio y un bloque de hormigón reforzado con FRP. La primera ayuda a mostrar la efectividad del método en sencillos escenarios de daño simple y múltiple, de forma que las conclusiones extraídas se aplican sobre los otros dos, diseñados para simular condiciones de despegue en distintas estructuras. Demostrada la efectividad del método de clasificación jerárquica de lecturas de impedancias, se aplica el procedimiento sobre las estructuras de hormigón armado reforzadas con bandas de FRP objeto de esta tesis, detectando y clasificando cada estado de daño. Finalmente, y como alternativa al anterior procedimiento, se propone un método para la monitorización continua de la interfase FRP-Hormigón, a través de una red de sensores FBG permanentemente instalados en dicha interfase. De esta forma, se obtienen medidas de deformación de la interfase en condiciones de carga continua, para ser implementadas en un modelo de optimización multiobjetivo, cuya solución se haya por medio de una expansión multiobjetivo del método Particle Swarm Optimization (PSO). La fiabilidad de este último método de detección se investiga a través de sendos ejemplos tanto numéricos como experimentales. ABSTRACT This thesis aims to develop robust and reliable damage identification methods focused on experimental structural systems, in particular Reinforced Concrete (RC) structures externally strengthened with Fiber Reinforced Polymers (FRP) strips. The failure mode of this type of structural system is critical, since it is usually due to sudden and brittle debonding of the FRP reinforcement originating from intermediate flexural cracks. Detection of the debonding in its initial stage is essential thus to prevent future failure, which might be catastrophic. Initially, a revision of the Electro-Mechanical Impedance (EMI) method is carried out, in order to expose its capabilities for local damage detection. Once the appropriate technology is selected, which includes impedance analyzer as well as novel PZT sensors for smart monitoring, an automated procedure has been design based on the impedance signatures of several lab-scale structures. On the basis that capturing impedance measurements is possible thanks to an adequately deployed PZT sensor network, the estimation of damage presence is done by analyzing the results of different damage indices obtained from the literature. In order to make this process automatic so that it is not necessary a priori knowledge of the EMI method to carry out an experimental test, a Graphical User Interface has been designed, turning the impedance measurements into an easy and intuitive procedure. Damage is then assessed through the analysis of the corresponding damage indices, trying to estimate not only the damage severity, but also its approximate location. The development of these tests on any kind of structure generates large amounts of data to be processed, and sometimes the information provided by damage indices is not enough to achieve a complete analysis of the structural health condition. In most of the cases, some damage patterns can be found in the data, but none a priori knowledge of the health condition is given for any structure. At this point, an important research on pattern recognition techniques has been carried out, particularly on unsupervised learning techniques, finding interesting applications in the medicine field. From this investigation, a creative and innovative idea arose: to detect and track the evolution of damage in different structures, as if it were a cancer propagating through a human body. In that sense, the impedance signatures are used to give intrinsic information of the health condition of the structure, so that the same clustering algorithms applied in the cancer research can be applied to the problem addressed in this dissertation. Hierarchical clustering is then applied since it also provides a graphical display of the clustered data, including quantitative and qualitative information about damage. The performance of this approach is firstly investigated using three lab-scale structures, such as a simple aluminium beam, a bolt-jointed aluminium beam and an FRP-strengthened concrete specimen. The first one shows the performance of the method on simple single and multiple damage scenarios, so that the first conclusions can be extracted and applied to the other two experimental tests, which are designed to simulate a debonding condition on different structures. Once the performance of the impedance-based hierarchical clustering method is proven to be successful, it is then applied to the structural system studied in this dissertation, the RC structures externally strengthened with FRP strips, where the debonding failure in the interface between the FRP and the concrete is successfully detected and classified, proving thus the feasibility of this method. Finally, as an alternative to the previous approach, a continuous monitoring procedure of the FRP-Concrete interface is proposed, based on an FBGsensors Network permanently deployed within that interface. In this way, strain measurements can be obtained under controlled loading conditions, and then they are used in order to implement a multi-objective model updating method solved by a multi-objective expansion of the Particle Swarm Optimization (PSO) method. The feasibility of this last proposal is investigated and successfully proven on both numerical and experimental RC beams strengthened with FRP.

Atomistic modeling of Ge damage accumulation, amorphization and solid phase epitaxial regrowth

En los últimos años, el Ge ha ganado de nuevo atención con la finalidad de ser integrado en el seno de las existentes tecnologías de microelectrónica. Aunque no se le considera como un canddato capaz de reemplazar completamente al Si en el futuro próximo, probalemente servirá como un excelente complemento para aumentar las propiedades eléctricas en dispositivos futuros, especialmente debido a su alta movilidad de portadores. Esta integración requiere de un avance significativo del estado del arte en los procesos de fabricado. Técnicas de simulación, como los algoritmos de Monte Carlo cinético (KMC), proporcionan un ambiente atractivo para llevar a cabo investigación y desarrollo en este campo, especialmente en términos de costes en tiempo y financiación. En este estudio se han usado, por primera vez, técnicas de KMC con el fin entender el procesado “front-end” de Ge en su fabricación, específicamente la acumulación de dañado y amorfización producidas por implantación iónica y el crecimiento epitaxial en fase sólida (SPER) de las capas amorfizadas. Primero, simulaciones de aproximación de clisiones binarias (BCA) son usadas para calcular el dañado causado por cada ión. La evolución de este dañado en el tiempo se simula usando KMC sin red, o de objetos (OKMC) en el que sólamente se consideran los defectos. El SPER se simula a través de una aproximación KMC de red (LKMC), siendo capaz de seguir la evolución de los átomos de la red que forman la intercara amorfo/cristalina. Con el modelo de amorfización desarrollado a lo largo de este trabajo, implementado en un simulador multi-material, se pueden simular todos estos procesos. Ha sido posible entender la acumulación de dañado, desde la generación de defectos puntuales hasta la formación completa de capas amorfas. Esta acumulación ocurre en tres regímenes bien diferenciados, empezando con un ritmo lento de formación de regiones de dañado, seguido por una rápida relajación local de ciertas áreas en la fase amorfa donde ambas fases, amorfa y cristalina, coexisten, para terminar en la amorfización completa de capas extensas, donde satura el ritmo de acumulación. Dicha transición ocurre cuando la concentración de dañado supera cierto valor límite, el cual es independiente de las condiciones de implantación. Cuando se implantan los iones a temperaturas relativamente altas, el recocido dinámico cura el dañado previamente introducido y se establece una competición entre la generación de dañado y su disolución. Estos efectos se vuelven especialmente importantes para iones ligeros, como el B, el cual crea dañado más diluido, pequeño y distribuido de manera diferente que el causado por la implantación de iones más pesados, como el Ge. Esta descripción reproduce satisfactoriamente la cantidad de dañado y la extensión de las capas amorfas causadas por implantación iónica reportadas en la bibliografía. La velocidad de recristalización de la muestra previamente amorfizada depende fuertemente de la orientación del sustrato. El modelo LKMC presentado ha sido capaz de explicar estas diferencias entre orientaciones a través de un simple modelo, dominado por una única energía de activación y diferentes prefactores en las frecuencias de SPER dependiendo de las configuraciones de vecinos de los átomos que recristalizan. La formación de maclas aparece como una consecuencia de esta descripción, y es predominante en sustratos crecidos en la orientación (111)Ge. Este modelo es capaz de reproducir resultados experimentales para diferentes orientaciones, temperaturas y tiempos de evolución de la intercara amorfo/cristalina reportados por diferentes autores. Las parametrizaciones preliminares realizadas de los tensores de activación de tensiones son también capaces de proveer una buena correlación entre las simulaciones y los resultados experimentales de velocidad de SPER a diferentes temperaturas bajo una presión hidrostática aplicada. Los estudios presentados en esta tesis han ayudado a alcanzar un mejor entendimiento de los mecanismos de producción de dañado, su evolución, amorfización y SPER para Ge, además de servir como una útil herramienta para continuar el trabajo en este campo. In the recent years, Ge has regained attention to be integrated into existing microelectronic technologies. Even though it is not thought to be a feasible full replacement to Si in the near future, it will likely serve as an excellent complement to enhance electrical properties in future devices, specially due to its high carrier mobilities. This integration requires a significant upgrade of the state-of-the-art of regular manufacturing processes. Simulation techniques, such as kinetic Monte Carlo (KMC) algorithms, provide an appealing environment to research and innovation in the field, specially in terms of time and funding costs. In the present study, KMC techniques are used, for the first time, to understand Ge front-end processing, specifically damage accumulation and amorphization produced by ion implantation and Solid Phase Epitaxial Regrowth (SPER) of the amorphized layers. First, Binary Collision Approximation (BCA) simulations are used to calculate the damage caused by every ion. The evolution of this damage over time is simulated using non-lattice, or Object, KMC (OKMC) in which only defects are considered. SPER is simulated through a Lattice KMC (LKMC) approach, being able to follow the evolution of the lattice atoms forming the amorphous/crystalline interface. With the amorphization model developed in this work, implemented into a multi-material process simulator, all these processes can be simulated. It has been possible to understand damage accumulation, from point defect generation up to full amorphous layers formation. This accumulation occurs in three differentiated regimes, starting at a slow formation rate of the damage regions, followed by a fast local relaxation of areas into the amorphous phase where both crystalline and amorphous phases coexist, ending in full amorphization of extended layers, where the accumulation rate saturates. This transition occurs when the damage concentration overcomes a certain threshold value, which is independent of the implantation conditions. When implanting ions at relatively high temperatures, dynamic annealing takes place, healing the previously induced damage and establishing a competition between damage generation and its dissolution. These effects become specially important for light ions, as B, for which the created damage is more diluted, smaller and differently distributed than that caused by implanting heavier ions, as Ge. This description successfully reproduces damage quantity and extension of amorphous layers caused by means of ion implantation reported in the literature. Recrystallization velocity of the previously amorphized sample strongly depends on the substrate orientation. The presented LKMC model has been able to explain these differences between orientations through a simple model, dominated by one only activation energy and different prefactors for the SPER rates depending on the neighboring configuration of the recrystallizing atoms. Twin defects formation appears as a consequence of this description, and are predominant for (111)Ge oriented grown substrates. This model is able to reproduce experimental results for different orientations, temperatures and times of evolution of the amorphous/crystalline interface reported by different authors. Preliminary parameterizations for the activation strain tensors are able to also provide a good match between simulations and reported experimental results for SPER velocities at different temperatures under the appliance of hydrostatic pressure. The studies presented in this thesis have helped to achieve a greater understanding of damage generation, evolution, amorphization and SPER mechanisms in Ge, and also provide a useful tool to continue research in this field.

Coupled Nonlinear Ginzburg-Landau and Mechanics Model for Martensitic Transformations in Polycrystals

Las transformaciones martensíticas (MT) se definen como un cambio en la estructura del cristal para formar una fase coherente o estructuras de dominio multivariante, a partir de la fase inicial con la misma composición, debido a pequeños intercambios o movimientos atómicos cooperativos. En el siglo pasado se han descubierto MT en diferentes materiales partiendo desde los aceros hasta las aleaciones con memoria de forma, materiales cerámicos y materiales inteligentes. Todos muestran propiedades destacables como alta resistencia mecánica, memoria de forma, efectos de superelasticidad o funcionalidades ferroicas como la piezoelectricidad, electro y magneto-estricción etc. Varios modelos/teorías se han desarrollado en sinergia con el desarrollo de la física del estado sólido para entender por qué las MT generan microstructuras muy variadas y ricas que muestran propiedades muy interesantes. Entre las teorías mejor aceptadas se encuentra la Teoría Fenomenológica de la Cristalografía Martensítica (PTMC, por sus siglas en inglés) que predice el plano de hábito y las relaciones de orientación entre la austenita y la martensita. La reinterpretación de la teoría PTMC en un entorno de mecánica del continuo (CM-PTMC) explica la formación de los dominios de estructuras multivariantes, mientras que la teoría de Landau con dinámica de inercia desentraña los mecanismos físicos de los precursores y otros comportamientos dinámicos. La dinámica de red cristalina desvela la reducción de la dureza acústica de las ondas de tensión de red que da lugar a transformaciones débiles de primer orden en el desplazamiento. A pesar de las diferencias entre las teorías estáticas y dinámicas dado su origen en diversas ramas de la física (por ejemplo mecánica continua o dinámica de la red cristalina), estas teorías deben estar inherentemente conectadas entre sí y mostrar ciertos elementos en común en una perspectiva unificada de la física. No obstante las conexiones físicas y diferencias entre las teorías/modelos no se han tratado hasta la fecha, aun siendo de importancia crítica para la mejora de modelos de MT y para el desarrollo integrado de modelos de transformaciones acopladas de desplazamiento-difusión. Por lo tanto, esta tesis comenzó con dos objetivos claros. El primero fue encontrar las conexiones físicas y las diferencias entre los modelos de MT mediante un análisis teórico detallado y simulaciones numéricas. El segundo objetivo fue expandir el modelo de Landau para ser capaz de estudiar MT en policristales, en el caso de transformaciones acopladas de desplazamiento-difusión, y en presencia de dislocaciones. Comenzando con un resumen de los antecedente, en este trabajo se presentan las bases físicas de los modelos actuales de MT. Su capacidad para predecir MT se clarifica mediante el ansis teórico y las simulaciones de la evolución microstructural de MT de cúbicoatetragonal y cúbicoatrigonal en 3D. Este análisis revela que el modelo de Landau con representación irreducible de la deformación transformada es equivalente a la teoría CM-PTMC y al modelo de microelasticidad para predecir los rasgos estáticos durante la MT, pero proporciona una mejor interpretación de los comportamientos dinámicos. Sin embargo, las aplicaciones del modelo de Landau en materiales estructurales están limitadas por su complejidad. Por tanto, el primer resultado de esta tesis es el desarrollo del modelo de Landau nolineal con representación irreducible de deformaciones y de la dinámica de inercia para policristales. La simulación demuestra que el modelo propuesto es consistente fcamente con el CM-PTMC en la descripción estática, y también permite una predicción del diagrama de fases con la clásica forma ’en C’ de los modos de nucleación martensítica activados por la combinación de temperaturas de enfriamiento y las condiciones de tensión aplicada correlacionadas con la transformación de energía de Landau. Posteriomente, el modelo de Landau de MT es integrado con un modelo de transformación de difusión cuantitativa para elucidar la relajación atómica y la difusión de corto alcance de los elementos durante la MT en acero. El modelo de transformaciones de desplazamiento y difusión incluye los efectos de la relajación en borde de grano para la nucleación heterogenea y la evolución espacio-temporal de potenciales de difusión y movilidades químicas mediante el acoplamiento de herramientas de cálculo y bases de datos termo-cinéticos de tipo CALPHAD. El modelo se aplica para estudiar la evolución microstructural de aceros al carbono policristalinos procesados por enfriamiento y partición (Q&P) en 2D. La microstructura y la composición obtenida mediante la simulación se comparan con los datos experimentales disponibles. Los resultados muestran el importante papel jugado por las diferencias en movilidad de difusión entre la fase austenita y martensita en la distibución de carbono en las aceros. Finalmente, un modelo multi-campo es propuesto mediante la incorporación del modelo de dislocación en grano-grueso al modelo desarrollado de Landau para incluir las diferencias morfológicas entre aceros y aleaciones con memoria de forma con la misma ruptura de simetría. La nucleación de dislocaciones, la formación de la martensita ’butterfly’, y la redistribución del carbono después del revenido son bien representadas en las simulaciones 2D del estudio de la evolución de la microstructura en aceros representativos. Con dicha simulación demostramos que incluyendo las dislocaciones obtenemos para dichos aceros, una buena comparación frente a los datos experimentales de la morfología de los bordes de macla, la existencia de austenita retenida dentro de la martensita, etc. Por tanto, basado en un modelo integral y en el desarrollo de códigos durante esta tesis, se ha creado una herramienta de modelización multiescala y multi-campo. Dicha herramienta acopla la termodinámica y la mecánica del continuo en la macroescala con la cinética de difusión y los modelos de campo de fase/Landau en la mesoescala, y también incluye los principios de la cristalografía y de la dinámica de red cristalina en la microescala. ABSTRACT Martensitic transformation (MT), in a narrow sense, is defined as the change of the crystal structure to form a coherent phase, or multi-variant domain structures out from a parent phase with the same composition, by small shuffles or co-operative movements of atoms. Over the past century, MTs have been discovered in different materials from steels to shape memory alloys, ceramics, and smart materials. They lead to remarkable properties such as high strength, shape memory/superelasticity effects or ferroic functionalities including piezoelectricity, electro- and magneto-striction, etc. Various theories/models have been developed, in synergy with development of solid state physics, to understand why MT can generate these rich microstructures and give rise to intriguing properties. Among the well-established theories, the Phenomenological Theory of Martensitic Crystallography (PTMC) is able to predict the habit plane and the orientation relationship between austenite and martensite. The re-interpretation of the PTMC theory within a continuum mechanics framework (CM-PTMC) explains the formation of the multivariant domain structures, while the Landau theory with inertial dynamics unravels the physical origins of precursors and other dynamic behaviors. The crystal lattice dynamics unveils the acoustic softening of the lattice strain waves leading to the weak first-order displacive transformation, etc. Though differing in statics or dynamics due to their origins in different branches of physics (e.g. continuum mechanics or crystal lattice dynamics), these theories should be inherently connected with each other and show certain elements in common within a unified perspective of physics. However, the physical connections and distinctions among the theories/models have not been addressed yet, although they are critical to further improving the models of MTs and to develop integrated models for more complex displacivediffusive coupled transformations. Therefore, this thesis started with two objectives. The first one was to reveal the physical connections and distinctions among the models of MT by means of detailed theoretical analyses and numerical simulations. The second objective was to expand the Landau model to be able to study MTs in polycrystals, in the case of displacive-diffusive coupled transformations, and in the presence of the dislocations. Starting with a comprehensive review, the physical kernels of the current models of MTs are presented. Their ability to predict MTs is clarified by means of theoretical analyses and simulations of the microstructure evolution of cubic-to-tetragonal and cubic-to-trigonal MTs in 3D. This analysis reveals that the Landau model with irreducible representation of the transformed strain is equivalent to the CM-PTMC theory and microelasticity model to predict the static features during MTs but provides better interpretation of the dynamic behaviors. However, the applications of the Landau model in structural materials are limited due its the complexity. Thus, the first result of this thesis is the development of a nonlinear Landau model with irreducible representation of strains and the inertial dynamics for polycrystals. The simulation demonstrates that the updated model is physically consistent with the CM-PTMC in statics, and also permits a prediction of a classical ’C shaped’ phase diagram of martensitic nucleation modes activated by the combination of quenching temperature and applied stress conditions interplaying with Landau transformation energy. Next, the Landau model of MT is further integrated with a quantitative diffusional transformation model to elucidate atomic relaxation and short range diffusion of elements during the MT in steel. The model for displacive-diffusive transformations includes the effects of grain boundary relaxation for heterogeneous nucleation and the spatio-temporal evolution of diffusion potentials and chemical mobility by means of coupling with a CALPHAD-type thermo-kinetic calculation engine and database. The model is applied to study for the microstructure evolution of polycrystalline carbon steels processed by the Quenching and Partitioning (Q&P) process in 2D. The simulated mixed microstructure and composition distribution are compared with available experimental data. The results show that the important role played by the differences in diffusion mobility between austenite and martensite to the partitioning in carbon steels. Finally, a multi-field model is proposed by incorporating the coarse-grained dislocation model to the developed Landau model to account for the morphological difference between steels and shape memory alloys with same symmetry breaking. The dislocation nucleation, the formation of the ’butterfly’ martensite, and the redistribution of carbon after tempering are well represented in the 2D simulations for the microstructure evolution of the representative steels. With the simulation, we demonstrate that the dislocations account for the experimental observation of rough twin boundaries, retained austenite within martensite, etc. in steels. Thus, based on the integrated model and the in-house codes developed in thesis, a preliminary multi-field, multiscale modeling tool is built up. The new tool couples thermodynamics and continuum mechanics at the macroscale with diffusion kinetics and phase field/Landau model at the mesoscale, and also includes the essentials of crystallography and crystal lattice dynamics at microscale.

Education Planning Evolution for Forest Engineering in Spain

Forest engineering in Spain has a long tradition and active presence in the engineering field. It is also one of the first educational institutions that shaped the Spanish technological panorama in the mid nineteenth century. The actual situation of the forest systems in Spain is the result of 166 years of observation, research, education and the application of specific techniques and principles that forest engineers acquired with the successive study plans that were implanted in educational institutions. In this paper, the planning historical process of education in Forest engineer is analyzed, differentiating between four historical periods. The analysis of the stages focuses on the contents of the study plans, the orientation towards educational objectives, the duration of the studies and the causes for the modifications that had an impact on the evolution through time within the framework of the acquired experience and the technological advances

An Escherichia coli strain with all chromosomal rRNA operons inactivated: Complete exchange of rRNA genes between bacteria

Current global phylogenies are built predominantly on rRNA sequences. However, an experimental system for studying the evolution of rRNA is not readily available, mainly because the rRNA genes are highly repeated in most experimental organisms. We have constructed an Escherichia coli strain in which all seven chromosomal rRNA operons are inactivated by deletions spanning the 16S and 23S coding regions. A single E. coli rRNA operon carried by a multicopy plasmid supplies 16S and 23S rRNA to the cell. By using this strain we have succeeded in creating microorganisms that contain only a foreign rRNA operon derived from either Salmonella typhimurium or Proteus vulgaris, microorganisms that have diverged from E. coli about 120–350 million years ago. We also were able to replace the E. coli rRNA operon with an E. coli/yeast hybrid one in which the GTPase center of E. coli 23S rRNA had been substituted by the corresponding domain from Saccharomyces cerevisiae. These results suggest that, contrary to common belief, coevolution of rRNA with many other components in the translational machinery may not completely preclude the horizontal transfer of rRNA genes.