Conceptual design of the liquid metal laboratory of the TECHNOFUSION facility

The application of liquid metal technology in fusion devices requires R&D related to many phenomena: interaction between liquid metals and structural material as corrosion, erosion and passivation techniques; magneto-hydrodynamics; free surface fluid-dynamics and any other physical aspect that will be needed for their safe reliable operation. In particular, there is a significant shortage of experimental facilities dedicated to the development of the lithium technology. In the framework of the TECHNOFUSION project, an experimental laboratory devoted to the lithium technology development is proposed, in order to shed some light in the path to IFMIF and the design of chamber's first wall and divertors. The conceptual design foresee a development in two stages, the first one consisting on a material testing loop. The second stage proposes the construction of a mock-up of the IFMIF target that will allow to assess the behaviour of a free-surface lithium target under vacuum conditions. In this paper, such conceptual design is addressed.

Modelización y validación experimental del transporte de tritio en metal líquido y otros materiales para reactores de fusión nuclear magnética e inercial

La presente tesis se centra en el estudio de los fenómenos de transporte de los isótopos de hidrógeno, y más concretamente del tritio, en materiales de interés para los reactores de fusión nuclear. Los futuros reactores de fusión nuclear necesitarán una Planta de Tritio, con una envoltura regeneradora (breeding blanket) y unos sistemas auxiliares claves para su diseño. Por lo tanto su desarrollo y cualificación son cruciales para demostrar que los reactores de fusión son una opción viable como futura fuente de energía. Se han resaltado los diferentes retos de la difusión y retención de estas especies ligeras para cada sistema de la Planta de Tritio, y se han identificado las necesidades experimentales y paramétricas para abordar las simulaciones de difusión, como factores de transporte como la difusividad, absorción/desorción, solubilidad y atrapamiento. Se han estudiado los fenómenos de transporte y parámetros del T en el metal líquido LiPb, componente del breeding blanket tanto para una planta de fusión magnética como inercial. Para ello se han utilizado dos experimentos con características diversas, uno de ellos se ha llevado a cabo en un reactor de alto flujo, y por lo tanto, en condiciones de irradiación, y el otro sin irradiación. Los métodos de simulación numérica aplicados se han adaptado a los experimentos para las mediciones y para estudiar el régimen de transporte. En el estudio de estos experimentos se ha obtenido un valor para algunos de los parámetros claves en el transporte y gestión del tritio en el reactor. Finalmente se realiza un cálculo de la acumulación y difusión de tritio en una primera pared de tungsteno para un reactor de fusión inercial. En concreto para el proyecto de fusión por láser europeo, HiPER (para sus fases 4a y 4b). Se ha estudiado: la implantación de los isótopos de H y He en la pared de W tras una reacción de fusión por iluminación directa con un láser de 48MJ; el efecto en el transporte de T de los picos de temperatura en el W debido a la frecuencia de los eventos de fusión; el régimen de transporte en la primera pared. Se han identificado la naturaleza de las trampas más importantes para el T y se ha propuesto un modelo avanzado para la difusión con trampas. ABSTRACT The present thesis focuses into study the transport phenomenons of hydrogen isotopes, more specifically tritium, in materials of interest for nuclear fusion reactors. The future nuclear reactors will be provided of a Tritium Plant, with its breeding blanket and its auxiliary systems, all of them essential components for the plant. Therefore a reliable development and coalification are key issues to prove the viability of the nuclear fusion reactors as an energy source. The currently challenges for the diffusion and accumulation of these light species for each system of the TP has been studied. Experimental and theoretical needs have been identified and analyzed, specially from the viewpoint of the parameters. To achieve reliable simulations of tritium transport, parameters as diffusivity, absorption/desorption, solubility and trapping must be reliables. Transport phenomenon and parameters of T in liquid metal have been studied. Lead lithium is a key component of the breeding blanket, either in magnetic or inertial fusion confinement. Having this aim in mind, two experiments with different characteristics have been used; one of them has been realized in a high flux reactor, and hence, under irradiation conditions. The other one has been realized without radiation. The mathematical methods for the simulation have been adapted to the experiments, for the measures and also to study the transport behavior. A value for some key parameters for tritium management has been obtained in these studies. Finally, tritium accumulation and diffusion in a W first wall of an inertial nuclear fusion reactor has been assessed. A diffusion model of the implanted H, D, T and He species for the two initial phases of the proposed European laser fusion Project HiPER (namely, phase 4a and phase 4b) has been implemented using Tritium Migration Analysis Program, TMAP7. The effect of the prompt and working temperatures and the operational pulsing modes on the diffusion are studied. The nature of tritium traps in W and their performance has been analyzed and discussed.

Monte Carlo Neutronics and Thermal-hydraulics coupling applied to Fast Reactors

Un escenario habitualmente considerado para el uso sostenible y prolongado de la energía nuclear contempla un parque de reactores rápidos refrigerados por metales líquidos (LMFR) dedicados al reciclado de Pu y la transmutación de actínidos minoritarios (MA). Otra opción es combinar dichos reactores con algunos sistemas subcríticos asistidos por acelerador (ADS), exclusivamente destinados a la eliminación de MA. El diseño y licenciamiento de estos reactores innovadores requiere herramientas computacionales prácticas y precisas, que incorporen el conocimiento obtenido en la investigación experimental de nuevas configuraciones de reactores, materiales y sistemas. A pesar de que se han construido y operado un cierto número de reactores rápidos a nivel mundial, la experiencia operacional es todavía reducida y no todos los transitorios se han podido entender completamente. Por tanto, los análisis de seguridad de nuevos LMFR están basados fundamentalmente en métodos deterministas, al contrario que las aproximaciones modernas para reactores de agua ligera (LWR), que se benefician también de los métodos probabilistas. La aproximación más usada en los estudios de seguridad de LMFR es utilizar una variedad de códigos, desarrollados a base de distintas teorías, en busca de soluciones integrales para los transitorios e incluyendo incertidumbres. En este marco, los nuevos códigos para cálculos de mejor estimación ("best estimate") que no incluyen aproximaciones conservadoras, son de una importancia primordial para analizar estacionarios y transitorios en reactores rápidos. Esta tesis se centra en el desarrollo de un código acoplado para realizar análisis realistas en reactores rápidos críticos aplicando el método de Monte Carlo. Hoy en día, dado el mayor potencial de recursos computacionales, los códigos de transporte neutrónico por Monte Carlo se pueden usar de manera práctica para realizar cálculos detallados de núcleos completos, incluso de elevada heterogeneidad material. Además, los códigos de Monte Carlo se toman normalmente como referencia para los códigos deterministas de difusión en multigrupos en aplicaciones con reactores rápidos, porque usan secciones eficaces punto a punto, un modelo geométrico exacto y tienen en cuenta intrínsecamente la dependencia angular de flujo. En esta tesis se presenta una metodología de acoplamiento entre el conocido código MCNP, que calcula la generación de potencia en el reactor, y el código de termohidráulica de subcanal COBRA-IV, que obtiene las distribuciones de temperatura y densidad en el sistema. COBRA-IV es un código apropiado para aplicaciones en reactores rápidos ya que ha sido validado con resultados experimentales en haces de barras con sodio, incluyendo las correlaciones más apropiadas para metales líquidos. En una primera fase de la tesis, ambos códigos se han acoplado en estado estacionario utilizando un método iterativo con intercambio de archivos externos. El principal problema en el acoplamiento neutrónico y termohidráulico en estacionario con códigos de Monte Carlo es la manipulación de las secciones eficaces para tener en cuenta el ensanchamiento Doppler cuando la temperatura del combustible aumenta. Entre todas las opciones disponibles, en esta tesis se ha escogido la aproximación de pseudo materiales, y se ha comprobado que proporciona resultados aceptables en su aplicación con reactores rápidos. Por otro lado, los cambios geométricos originados por grandes gradientes de temperatura en el núcleo de reactores rápidos resultan importantes para la neutrónica como consecuencia del elevado recorrido libre medio del neutrón en estos sistemas. Por tanto, se ha desarrollado un módulo adicional que simula la geometría del reactor en caliente y permite estimar la reactividad debido a la expansión del núcleo en un transitorio. éste módulo calcula automáticamente la longitud del combustible, el radio de la vaina, la separación de los elementos de combustible y el radio de la placa soporte en función de la temperatura. éste efecto es muy relevante en transitorios sin inserción de bancos de parada. También relacionado con los cambios geométricos, se ha implementado una herramienta que, automatiza el movimiento de las barras de control en busca d la criticidad del reactor, o bien calcula el valor de inserción axial las barras de control. Una segunda fase en la plataforma de cálculo que se ha desarrollado es la simulació dinámica. Puesto que MCNP sólo realiza cálculos estacionarios para sistemas críticos o supercríticos, la solución más directa que se propone sin modificar el código fuente de MCNP es usar la aproximación de factorización de flujo, que resuelve por separado la forma del flujo y la amplitud. En este caso se han estudiado en profundidad dos aproximaciones: adiabática y quasiestática. El método adiabático usa un esquema de acoplamiento que alterna en el tiempo los cálculos neutrónicos y termohidráulicos. MCNP calcula el modo fundamental de la distribución de neutrones y la reactividad al final de cada paso de tiempo, y COBRA-IV calcula las propiedades térmicas en el punto intermedio de los pasos de tiempo. La evolución de la amplitud de flujo se calcula resolviendo las ecuaciones de cinética puntual. Este método calcula la reactividad estática en cada paso de tiempo que, en general, difiere de la reactividad dinámica que se obtendría con la distribución de flujo exacta y dependiente de tiempo. No obstante, para entornos no excesivamente alejados de la criticidad ambas reactividades son similares y el método conduce a resultados prácticos aceptables. Siguiendo esta línea, se ha desarrollado después un método mejorado para intentar tener en cuenta el efecto de la fuente de neutrones retardados en la evolución de la forma del flujo durante el transitorio. El esquema consiste en realizar un cálculo cuasiestacionario por cada paso de tiempo con MCNP. La simulación cuasiestacionaria se basa EN la aproximación de fuente constante de neutrones retardados, y consiste en dar un determinado peso o importancia a cada ciclo computacial del cálculo de criticidad con MCNP para la estimación del flujo final. Ambos métodos se han verificado tomando como referencia los resultados del código de difusión COBAYA3 frente a un ejercicio común y suficientemente significativo. Finalmente, con objeto de demostrar la posibilidad de uso práctico del código, se ha simulado un transitorio en el concepto de reactor crítico en fase de diseño MYRRHA/FASTEF, de 100 MW de potencia térmica y refrigerado por plomo-bismuto. ABSTRACT Long term sustainable nuclear energy scenarios envisage a fleet of Liquid Metal Fast Reactors (LMFR) for the Pu recycling and minor actinides (MAs) transmutation or combined with some accelerator driven systems (ADS) just for MAs elimination. Design and licensing of these innovative reactor concepts require accurate computational tools, implementing the knowledge obtained in experimental research for new reactor configurations, materials and associated systems. Although a number of fast reactor systems have already been built, the operational experience is still reduced, especially for lead reactors, and not all the transients are fully understood. The safety analysis approach for LMFR is therefore based only on deterministic methods, different from modern approach for Light Water Reactors (LWR) which also benefit from probabilistic methods. Usually, the approach adopted in LMFR safety assessments is to employ a variety of codes, somewhat different for the each other, to analyze transients looking for a comprehensive solution and including uncertainties. In this frame, new best estimate simulation codes are of prime importance in order to analyze fast reactors steady state and transients. This thesis is focused on the development of a coupled code system for best estimate analysis in fast critical reactor. Currently due to the increase in the computational resources, Monte Carlo methods for neutrons transport can be used for detailed full core calculations. Furthermore, Monte Carlo codes are usually taken as reference for deterministic diffusion multigroups codes in fast reactors applications because they employ point-wise cross sections in an exact geometry model and intrinsically account for directional dependence of the ux. The coupling methodology presented here uses MCNP to calculate the power deposition within the reactor. The subchannel code COBRA-IV calculates the temperature and density distribution within the reactor. COBRA-IV is suitable for fast reactors applications because it has been validated against experimental results in sodium rod bundles. The proper correlations for liquid metal applications have been added to the thermal-hydraulics program. Both codes are coupled at steady state using an iterative method and external files exchange. The main issue in the Monte Carlo/thermal-hydraulics steady state coupling is the cross section handling to take into account Doppler broadening when temperature rises. Among every available options, the pseudo materials approach has been chosen in this thesis. This approach obtains reasonable results in fast reactor applications. Furthermore, geometrical changes caused by large temperature gradients in the core, are of major importance in fast reactor due to the large neutron mean free path. An additional module has therefore been included in order to simulate the reactor geometry in hot state or to estimate the reactivity due to core expansion in a transient. The module automatically calculates the fuel length, cladding radius, fuel assembly pitch and diagrid radius with the temperature. This effect will be crucial in some unprotected transients. Also related to geometrical changes, an automatic control rod movement feature has been implemented in order to achieve a just critical reactor or to calculate control rod worth. A step forward in the coupling platform is the dynamic simulation. Since MCNP performs only steady state calculations for critical systems, the more straight forward option without modifying MCNP source code, is to use the flux factorization approach solving separately the flux shape and amplitude. In this thesis two options have been studied to tackle time dependent neutronic simulations using a Monte Carlo code: adiabatic and quasistatic methods. The adiabatic methods uses a staggered time coupling scheme for the time advance of neutronics and the thermal-hydraulics calculations. MCNP computes the fundamental mode of the neutron flux distribution and the reactivity at the end of each time step and COBRA-IV the thermal properties at half of the the time steps. To calculate the flux amplitude evolution a solver of the point kinetics equations is used. This method calculates the static reactivity in each time step that in general is different from the dynamic reactivity calculated with the exact flux distribution. Nevertheless, for close to critical situations, both reactivities are similar and the method leads to acceptable practical results. In this line, an improved method as an attempt to take into account the effect of delayed neutron source in the transient flux shape evolutions is developed. The scheme performs a quasistationary calculation per time step with MCNP. This quasistationary simulations is based con the constant delayed source approach, taking into account the importance of each criticality cycle in the final flux estimation. Both adiabatic and quasistatic methods have been verified against the diffusion code COBAYA3, using a theoretical kinetic exercise. Finally, a transient in a critical 100 MWth lead-bismuth-eutectic reactor concept is analyzed using the adiabatic method as an application example in a real system.

Silica final lenses under HiPER laser fusion reactor operation conditions

We have studied the thermo-mechanical response and atomistic degradation of final lenses in HiPER project. Final silica lenses are squares of 75 × 75 cm2 with a thickness of 5 cm. There are two scenarios where lenses are located at 8 m from the centre: •HiPER 4a, bunches of 100 shots (maximum 5 DT shots <48 MJ at ≈0.1 Hz). No blanket in chamber geometry. •HiPER 4b, continuous mode with shots ≈50 MJ at 10 Hz to generate 0.5 GW. Liquid metal blanket in chamber design.

Diffusion of Zn into GaAs and AlGaAs from isothermal Liquid-phase epitaxy solutions

In this work we present results of zinc diffusion in GaAs using the liquid phase epitaxy technique from liquid solutions of Ga‐As‐Zn and Ga‐As‐Al‐Zn. Using silicon‐doped n‐GaAs substrates, working at a diffusion temperature of 850 °C, and introducing a dopant concentration ranging 1018–1019 cm−3, the most important findings regarding the diffusion properties are as follows: (a) zinc concentration in the solid depends on the square root of zinc atomic fraction in the liquid; (b) the diffusion is dominated by the interstitial‐substitutional process; (c) the diffusivity D varies as about C3 in the form D=2.9×10−67C3.05; (d) aluminum plays the role of the catalyst of the diffusion process, if it is introduced in the liquid solution, since it is found that D varies as (γAsXlAs)−1; (e) the zinc interstitial is mainly doubly ionized (Zn++i); (f) the zinc diffusion coefficient in Al0.85 Ga0.15 As is about four times greater than in GaAs; (g) by means of all these results, it is possible to control zinc diffusion processes in order to obtain optimized depth junctions and doping levels in semiconductor device fabrication.

Liquid lead and Lithium. A Comprehensive molecular dynamics study

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. One of the main issues in these programs is the problem of liquid metals breeder blanket behavior. Structural material of the blanket should meet high requirements because of extreme operating conditions. Therefore the knowledge of eutectic properties like optimal composition, physical and thermodynamic behavior or diffusion coefficients of Tritium are extremely necessary for current designs. In particular, the knowledge of the function linking the tritium concentration dissolved in liquid materials with the tritium partial pressure at a liquid/gas interface in equilibrium, CT=f(PT), is of basic importance because it directly impacts all functional properties of a blanket determining: tritium inventory, tritium permeation rate and tritium extraction efficiency. Nowadays, understanding the structure and behavior of this compound is a real goal in fusion engineering and materials science. Simulations of liquids can provide much information to the community; not only supplementing experimental data, but providing new tests of theories and ideas, making specific predictions that require experimental tests, and ultimately helping to lead to the deeper understanding and better predictive behavior.

Simulación atomística de metales líquidos y aleaciones para reactores de fusión nuclear = Atomistic simulation of liquid metals and alloys for nuclear fusion reactors

La presente tesis comprende un estudio de metales líquidos, Li, Pb y eutéctico Li17Pb en el ámbito de la tecnología de fusión nuclear. Uno de los problemas fundamentales en futuros reactores de fusión es la producción y extracción de tritio en la denominada envoltura regeneradora (blanket en inglés). Dicho blanket tendrá dos propósitos, la extracción del calor generado por las reacciones de fusión para su posterior conversión en energía eléctrica así como la producción de tritio para realimentar el proceso. Dicha producción se realizará mediante el “splitting” del Li con los neutrones provenientes de la fusión. Esta reacción produce tritio y helio por lo que la interacción del T y el He con el metal líquido, con los materiales estructurales así como con el He es un problema fundamental aun no bien entendido y de gran importancia para futuros diseños. Los capítulos 1 2 y 3 presentan una introducción a dichos problemas. El capítulo 1 introduce al lector en la tecnología de fusión nuclear. El segundo capítulo explica en mayor detalle el uso de metales líquidos en reactores de fusión, no solo en blankets sino también como primera pared, divertor etc, lo que se denomina en general “plasma facing materials”. Por último se ofrece una breve introducción a las técnicas de dinámica molecular clásica (CMD) y un breve resumen de los potenciales más usados. El estudio se ha llevado a cabo utilizando simulación atomística mediante potenciales semi-empíricos del tipo átomo embebido (EAM). La Tesis consta de 4 partes bien definidas. En primer lugar se verificó la idoneidad de los potenciales tipo EAM para simular las propiedades de los metales Li y Pb en fase líquida. Dicho estudio se detalla en el Capítulo 4 y en su extensión, el Apéndice 1, en el que se estudia los límites de validez de esta aproximación. Los resultados de dicho estudio han sido publicados y presentados en diversos congresos internacionales. Un resumen de la metodología seguida fue publicado como capítulo de libro en Technofusión 2011. Los resultados se presentaron en diversos congresos internacionales, entre ellos ICENES 2011, (Artículo en ICENES Proceedings) ICOPS-SOFE 2011, en una presentación oral etc. El trabajo ha sido aceptado recientemente en Journal of Nuclear Materiales (Fraile et al 2012). La segunda parte y más importante comprende el desarrollo de un potencial para el estudio de la mezcla de ambos metales. Éste es el trabajo más novedoso e importante dado que no existía en la literatura un potencial semejante. Se estudiaron dos aproximaciones distintas al problema, un potencial tipo EAM/cd y un potencial EAM/alloy. Ambos potenciales dan resultados satisfactorios para la simulación del eutéctico (y concentraciones de Li menores que el 17%). Sin embargo el sistema LiPb en todas las concentraciones es un sistema que se aparta enormemente de una solución ideal y dicho potencial no daba buenos resultados para mezclas PbLi con concentraciones de Li grandes. Este problema fue solventado mediante el desarrollo de un segundo potencial, esta vez tipo EAM/alloy (segunda parte del Capítulo 5). Dicho trabajo será enviado a Physical Review Letters o a Phys. Rev. B, y una extensión junto con un estudio detallado de las propiedades del eutéctico de acuerdo con nuestras simulaciones se enviará a continuación a Phys. Rev. B. En tercer lugar se estudió el problema de la difusividad del H en metales líquidos aprovechando distintos potenciales existentes en la literatura. El problema del H en metales líquidos es importante en metalurgia. En dicho capítulo se estudió la difusividad del H en Pd, Ni y Al con potenciales tipo EAM, y también con un potencial más sofisticado que tiene en cuenta la dependencia angular de las interacciones (ADP por sus siglas en inglés). De este modo disponemos de un estudio detallado del problema con diferentes modelos y diferentes metales. La conclusión apunta a que si se compara con los resultados experimentales (muy escasos) los resultados obtenidos mediante CMD dan valores bajos de la difusividad del H. Las razones de dicho desacuerdo entre simulación y experimentos se detallan en el Capítulo 6. Este trabajo ha sido presentado en una presentación oral en el reciente congreso internacional “Trends on Nanothecnology” TNT 2012 celebrado en Madrid. El trabajo será publicado en un futuro próximo. Por último, como se dijo anteriormente, el estudio del He, la formación de burbujas en metales líquidos, su difusión nucleación y cavitación es otro problema deseable de ser estudiado mediante técnicas atomísticas. Para ello es necesario el desarrollo de diversos potenciales, He-Li, He-Pb y un potencial ternario Pb-Li-He. Para ello se han realizado simulación ab initio de los sistemas Pb+He y Li+He. Dicho estudio pretende calcular las fuerzas entre los átomos del metal (Pb o Li) con intersticiales de He. De este modo aplicaremos el “force matching method” (FMM) para el desarrollo de dichos potenciales. En el Capítulo 7 se detallan los resultados obtenidos referidos a las posiciones más favorables de las impurezas de He dentro de redes cristalinas de Pb y Li así como el efecto de tener en cuenta el acoplo spin-orbita (SOC en inglés). El análisis de los resultados en términos de transferencia de carga y análisis de las densidades electrónicas, así como la creación de los potenciales mencionados está en progreso. En conjunto la tesis presenta un estudio de los diversos problemas relacionados con el uso de metales líquidos en reactores de fusión y representa un primer paso en la determinación de parámetros de gran importancia para el diseño de blankets y sistemas de primera pared. Con la simulación MD de dichos problemas mediante, importante, potenciales realistas, valores de difusión, solubilidad etc de especies ligeras, H (o sus isotopos) y He en metales líquidos podrá ser calculada complementando así la base de datos que presenta enormes incertidumbres.

Estudios de especiación de arsénico y acumulación de metales en muestras de interés medioambiental : Arsenic speciation and metal accumulation studies in environmental samples

Se ha estudiado la determinación de especies de arsénico y de contenidos totales de arsénico y metales pesados, específicamente cadmio, cromo, cobre, níquel, plomo y cinc, en muestras de interés medioambiental por su elevada capacidad acumuladora de metales, concretamente algas marinas comestibles y plantas terrestres procedentes de suelos contaminados por la actividad minera. La determinación de contenidos totales se ha llevado a cabo mediante espectrometría de emisión atómica con plasma de acoplamiento inductivo (ICP‐AES), así como por espectrometría de fluorescencia atómica con generación de hidruros (HG‐AFS), para bajos contenidos de arsénico. Las muestras fueron mineralizadas en medio ácido y calentamiento en horno de microondas. Los métodos fueron validados a través de su aplicación a materiales de referencia de matriz similar a la de las muestras, certificados en contenidos totales de los elementos seleccionados. Los resultados obtenidos mostraron su elevada capacidad de bioabsorción, especialmente en relación a los elevados contenidos de arsénico encontrados en algunas especies de algas pardas (Phaeophytas). En las plantas, se calcularon los factores de translocación, acumulación y biodisponibilidad de los elementos estudiados, permitiendo identificar a la especie Corrigiola telephiifolia como posible acumuladora de plomo e hiperacumuladora de arsénico. La determinación de especies de arsénico hidrosolubles en las muestras objeto de estudio, se llevó a cabo por cromatografía líquida de alta eficacia (HPLC) acoplado a ICP‐AES, HG‐ICP‐AES y HG‐AFS, incluyendo una etapa previa de foto‐oxidación. Los métodos desarrollados, mediante intercambio aniónico y catiónico, permitieron la diferenciación de hasta once especies de arsénico. Para el análisis de las muestras, fue necesaria la optimización de métodos de extracción, seleccionándose la extracción asistida por microondas (MAE) con agua desionizada. Asimismo, se realizaron estudios de estabilidad de arsénico total y de las especies hidrosolubles presentes en las algas, tanto sobre la muestra sólida como en sus extractos acuosos, evaluando las condiciones de almacenamiento adecuadas. En el caso de las plantas, la aplicación del diseño factorial de experimentos permitió optimizar el método de extracción y diferenciar entre las especies de arsénico presentes en forma de iones sencillos de mayor movilidad y el arsénico más fuertemente enlazado a componentes estructurales. Los resultados obtenidos permitieron identificar la presencia de arseniato (As(V)) y arsenito (As(III)) en las plantas, así como de ácido monometilarsónico (MMA) y óxido de trimetilarsina (TMAO) en algunas especies. En la mayoría de las algas se encontraron especies tóxicas, tanto mayoritarias (arseniato) como minoritarias (ácido dimetilarsínico (DMA)), así como hasta cuatro arsenoazúcares. Los resultados obtenidos y su estudio a través de la legislación vigente, mostraron la necesidad de desarrollar una reglamentación específica para el control de este tipo de alimentos. La determinación de especies de arsénico liposolubles en las muestras de algas se llevó a cabo mediante HPLC, en modo fase inversa, acoplado a espectrometría de masas con plasma de acoplamiento inductivo (ICP‐MS) y con ionización por electrospray (ESI‐MS), permitiendo la elucidación estructural de estos compuestos a través de la determinación de sus masas moleculares. Para ello, fue necesaria la puesta a punto de métodos extracción y purificación de los extractos. La metodología desarrollada permitió identificar hasta catorce especies de arsénico liposolubles en las algas, tres de ellas correspondientes a hidrocarburos que contienen arsénico, y once a arsenofosfolípidos, además de dos especies desconocidas. Las masas moleculares de las especies identificadas fueron confirmadas mediante cromatografía de gases acoplada a espectrometría de masas (GC‐MS) y espectrometría de masas de alta resolución (HR‐MS). ABSTRACT The determination of arsenic species and total arsenic and heavy metal contents (cadmium, chromium, cooper, nickel, lead and zinc) in environmental samples, with high metal accumulator capacity, has been studied. The samples studied were edible marine algae and terrestrial plants from soils polluted by mining activities. The determination of total element contents was performed by inductively coupled plasma atomic emission spectrometry (ICP‐AES), as well as by hydride generation atomic fluorescence spectrometry (HG‐AFS) for low arsenic contents. The samples studied were digested in an acidic medium by heating in a microwave oven. The digestion methods were validated against reference materials, with matrix similar to sample matrix and certified in total contents of the elements studied. The results showed the high biosorption capacity of the samples studied, especially regarding the high arsenic contents in some species of brown algae (Phaeophyta division). In terrestrial plants, the translocation, accumulation and bioavailability factors of the elements studied were calculated. Thus, the plant species Corrigiola telephiifolia was identified as possible lead accumulator and arsenic hyperaccumulator. The determination of water‐soluble arsenic species in the samples studied was carried out by high performance liquid chromatography (HPLC) coupled to ICP‐AES, HG‐ICP‐AES and HG‐AFS, including a prior photo‐oxidation step. The chromatographic methods developed, by anion and cation exchange, allowed us to differentiate up to eleven arsenic species. The sample analysis required the optimization of extraction methods, choosing the microwave assisted extraction (MAE) with deionized water. On the other hand, the stability of total arsenic and water‐soluble arsenic species in algae, both in the solid samples and in the water extracts, was studied, assessing the suitable storage conditions. In the case of plant samples, the application of a multivariate experimental design allowed us to optimize the extraction method and differentiate between the arsenic species present as simple ions of higher mobility and the arsenic more strongly bound to structural components. The presence of arsenite (As(III)) and arsenate (As(V)) was identified in plant samples, as well as monomethylarsonic acid (MMA) and trimethylarsine oxide (TMAO) in some cases. Regarding algae, toxic arsenic species were found in most of them, both As(V) and dimethylarsinic acid (DMA), as well as up to four arsenosugars. These results were discussed according to the current legislation, showing the need to develop specific regulations to control this kind of food products. The determination of lipid‐soluble arsenic species in alga samples was performed by reversed‐phase HPLC coupled to inductively coupled plasma and electrospray mass spectrometry (ICP‐MS and ESI‐MS), in order to establish the structure of these compounds by determining the corresponding molecular masses. For this purpose, it was necessary to develop an extraction method, as well as a clean‐up method of the extracts. The method developed permitted the identification of fourteen lipid‐soluble arsenic compounds in algae, corresponding to three arsenic‐hydrocarbons and eleven arsenosugarphospholipids, as well as two unknown compounds. Accurate mass measurements of the identified compounds were performed by gas chromatography coupled to mass spectrometry (GC‐MS) and high resolution mass spectrometry (HR‐MS).

Engineering metal precipitate size distributions to enhance gettering in multicrystalline silicon

The extraction of metal impurities during phosphorus diffusion gettering (PDG) is one of the crucial process steps when fabricating high-efficiency solar cells using low-cost, lower-purity silicon wafers. In this work, we show that for a given metal concentration, the size and density of metal silicide precipitates strongly influences the gettering efficacy. Different precipitate size distributions can be already found in silicon wafers grown by different techniques. In our experiment, however, the as-grown distribution of precipitated metals in multicrystalline Si sister wafers is engineered through different annealing treatments in order to control for the concentration and distribution of other defects. A high density of small precipitates is formed during a homogenization step, and a lower density of larger precipitates is formed during extended annealing at 740º C. After PDG, homogenized samples show a decreased interstitial iron concentration compared to as-grown and ripened samples, in agreement with simulations.

Degradation in seawater of structural adhesives for hybrid fibre-metal laminated materials

The adhesives used for applications in marine environments are subject to particular chemical conditions, which are mainly characterised by an elevated chlorine ion content and intermittent wetting/drying cycles, among others.These conditions can limit the use of adhesives due to the degradation processes that they experience. In this work, the chemical degradation of two different polymers, polyurethane and vinylester, was studied in natural seawater under immersion for different periods of time.The diffusion coefficients and concentration profiles of water throughout the thickness of the adhesiveswere obtained.Microstructural changes in the polymer due to the action of water were observed by SEM, and the chemical degradation of the polymer was monitored with the Fourier transform infrared spectroscopy (FTIR) and differential scanning calorimetry (DSC). The degradation of the mechanical properties of the adhesive was determined by creep tests withMixed Cantilever Beam (MCB) specimens at different temperatures. After 180 days of immersion of the specimens, it was concluded that the J-integral value (depending on the strain) implies a loss of stiffness of 51% and a decrease in the failure load of 59% for the adhesive tested.

Understanding phosphorus diffusion into silicon in a MOVPE environment for III-V on silicon solar cells

Dual-junction solar cells formed by a GaAsP or GaInP top cell and a silicon bottom cell seem to be attractive candidates to materialize the long sought-for integration of III?V materials on silicon for photovoltaic applications. When manufacturing a multi-junction solar cell on silicon, one of the first processes to be addressed is the development of the bottom subcell and, in particular, the formation of its emitter. In this study, we analyze, both experimentally and by simulations, the formation of the emitter as a result of phosphorus diffusion that takes place during the first stages of the epitaxial growth of the solar cell. Different conditions for the Metal-Organic Vapor Phase Epitaxy (MOVPE) process have been evaluated to understand the impact of each parameter, namely, temperature, phosphine partial pressure, time exposure and memory effects in the final diffusion profiles obtained. A model based on SSupremIV process simulator has been developed and validated against experimental profiles measured by ECV and SIMS to calculate P diffusion profiles in silicon formed in a MOVPE environment taking in consideration all these factors.

Functional and expression analysis of the metal-inducible dmeRF system from Rhizobium legumionosarum bv. viciae

A gene encoding a homolog to the cation diffusion facilitator protein DmeF from Cupriavidus metallidurans has been identified in the genome of Rhizobium leguminosarum UPM791. The R. leguminosarum dmeF gene is located downstream of an open reading frame (designated dmeR) encoding a protein homologous to the nickel- and cobalt-responsive transcriptional regulator RcnR from Escherichia coli. Analysis of gene expression showed that the R. leguminosarum dmeRF genes are organized as a transcriptional unit whose expression is strongly induced by nickel and cobalt ions, likely by alleviating the repressor activity of DmeR on dmeRF transcription. An R. leguminosarum dmeRF mutant strain displayed increased sensitivity to Co(II) and Ni(II), whereas no alterations of its resistance to Cd(II), Cu(II), or Zn(II) were observed. A decrease of symbiotic performance was observed when pea plants inoculated with an R. leguminosarum dmeRF deletion mutant strain were grown in the presence of high concentrations of nickel and cobalt. The same mutant induced significantly lower activity levels of NiFe hydrogenase in microaerobic cultures. These results indicate that the R. leguminosarum DmeRF system is a metal-responsive efflux mechanism acting as a key element for metal homeostasis in R. leguminosarum under free-living and symbiotic conditions. The presence of similar dmeRF gene clusters in other Rhizobiaceae suggests that the dmeRF system is a conserved mechanism for metal tolerance in legume endosymbiotic bacteria.

Functional characterization of Rhizobium leguminosarum bv. viciae DmeRF, a cation diffusion facilitator system involved in nickel and cobalt resistance

In prokaryotes, nickel is an essential element participating in the structure of enzymes involved in multiple cellular processes. Nickel transport is a challenge for microorganisms since, although essential, high levels of this metal inside the cell are toxic. For this reason, bacteria have developed high-affinity nickel transporters as well as nickel-specific detoxification systems. Ultramafic soils, and soils contaminated with heavy metals are excellent sources of nickel resistant bacteria. Molecular analysis of strains isolated in the habitats has revealed novel genetic systems involved in adaptation to such hostile conditions.

Impact of metal-organic vapor phase epitaxy environment on silicon bulk lifetime for III–V-on-Si multijunction solar cells

With the final goal of integrating III-V materials on silicon substrates for tandem solar cells, the influence of the Metal-Organic Vapor Phase Epitaxy (MOVPE) environment on the minority carrier properties of silicon wafers has been evaluated. These properties will essentially determine the photovoltaic performance of the bottom cell in a III-V-on-Si tandem solar cell. A comparison of the base minority carrier lifetimes obtained for different thermal processes carried out in a MOVPE reactor on Czochralski silicon wafers has been carried out. An important degradation of minority carrier lifetime during the surface preparation (i.e. H2 anneal) has been observed. Three different mechanisms have been proposed for explaining this behavior: 1) the introduction of extrinsic impurities coming from the reactor; 2) the activation of intrinsic lifetime killing impurities coming from the wafer itself; and finally, 3) the formation of crystal defects, which eventually become recombination centers. The effect of the emitter formation by phosphorus diffusion has also been evaluated. In this sense, it has been reported that lifetime can be recovered during the emitter formation either by the effect of the P on extracting impurities, or by the role of the atomic hydrogen on passivating the defects.

Optimizing phosphorus diffusion for photovoltaic applications: peak doping, inactive phosphorus,gettering, and contact formation

The phosphosilicate glass (PSG), fabricated by tube furnace diffusion using a POCl3 source, is widely used as a dopant source in the manufacturing of crystalline silicon solar cells. Although it has been a widely addressed research topic for a long time, there is still lack of a comprehensive understanding of aspects such as the growth, the chemical composition, possible phosphorus depletion, the resulting in-diffused phosphorus profiles, the gettering behavior in silicon, and finally the metal-contact formation. This paper addresses these different aspects simultaneously to further optimize process conditions for photovoltaic applications. To do so, a wide range of experimental data is used and combined with device and process simulations, leading to a more comprehensive interpretation. The results show that slight changes in the PSG process conditions can produce high-quality emitters. It is predicted that PSG processes at 860 °C for 60 min in combination with an etch-back and laser doping from PSG layer results in high-quality emitters with a peak dopant density Npeak = 8.0 × 1018 cm−3 and a junction depth dj = 0.4 μm, resulting in a sheet resistivityρsh = 380 Ω/sq and a saturation current-density J0 below 10 fA/cm2. With these properties, the POCl3 process can compete with ion implantation or doped oxide approaches.