MMonCa: A flexible and powerful new Kinetic Monte Carlo simulator

Kinetic Monte Carlo (KMC) is a widely used technique to simulate the evolution of radiation damage inside solids. Despite de fact that this technique was developed several decades ago, there is not an established and easy to access simulating tool for researchers interested in this field, unlike in the case of molecular dynamics or density functional theory calculations. In fact, scientists must develop their own tools or use unmaintained ones in order to perform these types of simulations. To fulfil this need, we have developed MMonCa, the Modular Monte Carlo simulator. MMonCa has been developed using professional C++ programming techniques and has been built on top of an interpreted language to allow having a powerful yet flexible, robust but customizable and easy to access modern simulator. Both non lattice and Lattice KMC modules have been developed. We will present in this conference, for the first time, the MMonCa simulator. Along with other (more detailed) contributions in this meeting, the versatility of MMonCa to study a number of problems in different materials (particularly, Fe and W) subject to a wide range of conditions will be shown. Regarding KMC simulations, we have studied neutron-generated cascade evolution in Fe (as a model material). Starting with a Frenkel pair distribution we have followed the defect evolution up to 450 K. Comparison with previous simulations and experiments shows excellent agreement. Furthermore, we have studied a more complex system (He-irradiated W:C) using a previous parametrization [1]. He-irradiation at 4 K followed by isochronal annealing steps up to 500 K has been simulated with MMonCa. The He energy was 400 eV or 3 keV. In the first case, no damage is associated to the He implantation, whereas in the second one, a significant Frenkel pair concentration (evolving into complex clusters) is associated to the He ions. We have been able to explain He desorption both in the absence and in the presence of Frenkel pairs and we have also applied MMonCa to high He doses and fluxes at elevated temperatures. He migration and trapping dominate the kinetics of He desorption. These processes will be discussed and compared to experimental results. [1] C.S. Becquart et al. J. Nucl. Mater. 403 (2010) 75

Effect of ion flux on helium retention in helium-irradiated tungsten

Helium retention in irradiated tungsten leads to swelling, pore formation, sample exfoliation and embrittlement with deleterious consequences in many applications. In particular, the use of tungsten in future nuclear fusion plants is proposed due to its good refractory properties. However, serious concerns about tungsten survivability stems from the fact that it must withstand severe irradiation conditions. In magnetic fusion as well as in inertial fusion (particularly with direct drive targets), tungsten components will be exposed to low and high energy ion (helium) irradiation, respectively. A common feature is that the most detrimental situations will take place in pulsed mode, i.e., high flux irradiation. There is increasing evidence on a correlation between a high helium flux and an enhancement of detrimental effects on tungsten. Nevertheless, the nature of these effects is not well understood due to the subtleties imposed by the exact temperature profile evolution, ion energy, pulse duration, existence of impurities and simultaneous irradiation with other species. Physically based Kinetic Monte Carlo is the technique of choice to simulate the evolution of radiation-induced damage inside solids in large temporal and space scales. We have used the recently developed code MMonCa (Modular Monte Carlo simulator), presented in this conference for the first time, to study He retention (and in general defect evolution) in tungsten samples irradiated with high intensity helium pulses. The code simulates the interactions among a large variety of defects and impurities (He and C) during the irradiation stage and the subsequent annealing steps. In addition, it allows us to vary the sample temperature to follow the severe thermo-mechanical effects of the pulses. In this work we will describe the helium kinetics for different irradiation conditions. A competition is established between fast helium cluster migration and trapping at large defects, being the temperature a determinant factor. In fact, high temperatures (induced by the pulses) are responsible for large vacancy cluster formation and subsequent additional trapping with respect to low flux irradiation.

The use of the ESS Bilbao installations for laser fusion studies

In this work the use of ESS-Bilbao fast neutron lines for irradiation of materials for nuclear fusion is studied. For the comparison of ESS-Bilbao with an inertial fusion facility a simplified model of HiPER chamber has been used. Several positions for irradiation at ESS-Bilbao have been also compared. The material chosen for the damage analysis is silica due to its importance on IFC optics. In this work a detailed comparison between the two facilities for silica irradiation is given. The comparison covers the neutron fluxes, doses, defect production and PKA spectra. This study is also intended as a methodological approach or guideline for future works on other materials.

Ion-beam damage induced in SiO2: Excitonic mechanisms under heavy electronic excitation

The Centro de Micro-Análisis de Materiales (CMAM) in the Universidad Autónoma de Madrid is carrying out an extensive research program on the processes induced by high energy heavy mass ions (SHI) on dielectric materials and their photonic applications [1?21]. A significant part of this activity constitutes a relevant contribution to the scientific program associated to the TECHNOFUSION project. It is performed in collaboration with the Instituto de Fusion Nuclear at the UPM, the CIEMAT, the Departamento de Física de Materiales at UAM and several other national institutions (INTA) and international laboratories (GANIL, France), Legnaro Italy, Grenoble?. The program has led to a large number of publications in reputed international journals.

Silica final lens performance in laser fusion facilities: HiPER and LIFE

Silica final lens performance in laser fusion facilities: HiPER and LIFE

Monitorización y seguimiento del esfuerzo realizado por los estudiantes y de su asistencia a actividades presenciales

Este artículo documenta el planteamiento, la metodología y los primeros resultados de un plan de monitorización detallada del esfuerzo y de asistencia a actividades presenciales por parte de los estudiantes de las titulaciones ofertadas por la Escuela Técnica Superior de Ingenieros Navales de la Universidad Politécnica de Madrid durante el segundo cuatrimestre del curso 2011-2012. Se ha establecido un sistema mecánico de recogida de datos de esfuerzo por parte de los estudiantes utilizando una hoja tipo test especialmente configurada al efecto. Se pasa una hoja en todas y cada una de las actividades presenciales realizadas y en la hoja se solicita información sobre el trabajo "fuera de clase". Se documenta en este artículo cómo se ha estructurado esa hoja, qué tipo de datos se recogen, cómo se tratan mediante una base de datos creada al efecto, qué tipo de análisis se puede realizar y qué resultados preliminares obtenemos de dichos análisis.

Influence of the composition of natural gas-hydrogen blends (HCNG) on performance an CO2 emissions of a spark ignition engine

Addition of hydrogen to natural gas could be a short-term alternative to nowadays fossil fuels as the emissions of greenhouse gases may be reduced. The aim of this study is to evaluate the performance and emissions of a park ignition engine fuelled with pure natural gas, pure hydrogen and different blends of hydrogen and natural gas (HCNG). The increase of the hydrogen fraction leads to variations in the cylinder pressure and CO2 emissions. In this work, a combustion model based on thermodynamic equations is used considering separated zones for the burned and unburned gases. The results show that the maximum cylinder pressure gets higher as the fraction of hydrogen in the blend increases. The presence of hydrogen in the blend leads to a drecrease in the CO2 emissions. Due to hydrogen properties, leaner fuel-air mixtures can be used along with the appropiate spark timing, leading to an engine emissions improvement without a performance worsening.

CO2 emissions from a spark ignition engine operatin on natural-hydrogen blends (HCNG)

The addition of hydrogen to natural gas could be a short-term alternative to today’s fossil fuels, as greenhouse gas emissions may be reduced. The aim of this study is to evaluate the emissions and performance of a spark ignition engine fuelled by pure natural gas, pure hydrogen, and different blends of hydrogen and natural gas (HCNG). Increasing the hydrogen fraction leads to variations in cylinder pressure and CO2 emissions. In this study, a combustion model based on thermodynamic equations is used, considering separate zones for burned and unburned gases. The results show that the maximum cylinder pressure rises as the fraction of hydrogen in the blend increases. The presence of hydrogen in the blend leads to a decrease in CO2 emissions. Due to the properties of hydrogen, leaner fuel–air mixtures can be used along with the appropriate spark timing, leading to an improvement in engine emissions with no loss of performance.

Desarrollo de un laboratorio virtual para la medida de ciclos indicados de un motor alternativo policombustible de encendido provocado

Se ha desarrollado un laboratorio virtual para la medida de ciclos indicados en un motor alternativo policombustible de encendido provocado. El laboratorio virtual desarrollado permite por una parte simular y obtener el ciclo indicado de un motor alternativo de cuatro tiempos, pudiendo modificar los parámetros más importantes de operación (régimen de giro, presiones de admisión y escape, temperaturas, etc.) y de diseño (parámetros geométricos del motor). Por otra parte, el laboratorio virtual permite simular el ensayo del motor en banco de pruebas y la medida del ciclo indicado. Los modelos matemáticos necesarios se han escrito en lenguaje Fortran, y estos interaccionan con un interfaz gráfico de usuario (GUI) programado en VEE®. El laboratorio virtual puede implementarse en una plataforma virtual de enseñanza de manera que los alumnos puedan tener acceso al sistema desde cualquier ordenador, en cualquier momento y desde cualquier lugar.

Microtermometría, direcciones de flujo hidrotermal y metalogenia del yacimiento de tungsteno de Pasto Bueno (Ancash, Perú)

El yacimiento de Pasto Bueno se localiza en el extremo nordeste del Batolito de la Cordillera Blanca, comprende diversas vetas, generalmente subverticales, asociadas al stock cuarzomonzonitico de Consuzo, datado como Terciario Superior, que intruye a las pizarras de la fm. Chicama y cuarcitas de la fm. Chimu. Las principales vetas discurren con direccion N-S cortando al stock, aunque tambien existen sistemas NE-SW asi como NW-SE encajados sobre las rocas metamorficas. La mineralogia de mena reconocida comprende wolframita (hubnerita), tetraedrita/tenantita, esfalerita y galena, en una ganga de cuarzo, fluorita, sericita, pirita y carbonatos, ademas de molibdenita, calcopirita, bornita, arsenopirita, enargita (luzonita), stolzita, scheelita, zinnwaldita, topacio, tungstita y arsenico nativo. Estudios previos han caracterizado Pasto Bueno como un yacimiento con una gran componente de greisen, con una evolucion de las vetas desde un episodio temprano esteril de 400 oC, depositando la mineralizacion economica en torno a los 220-250 oC y con un evento postumo de 175-220 oC rico en CO2. La precipitacion de la wolframita se produjo a partir de un fluido netamente hidrotermal, sin embargo, dicha precipitacion estuvo controlada por el aporte al sistema de aguas externas meteoricas y/o metamorficas. El trabajo llevado a cabo ha consistido en la realizacion de un estudio microtermometrico de las 3 principales estructuras del distrito: Consuelo, Alonso-Fenix y Chabuca, para caracterizar la evolucion del fluido mineralizador desde el stock (veta Consuelo) hacia las rocas metasedimentarias de las fm. Chicama y Chimu (manto Alonso-Fenix y veta Chabuca). Para ello se realizo un muestreo sobre el evento principal de mineralizacion. Dichas muestras se sometieron a un estudio petrografico de lamina gruesa para seleccionar las muestras optimas para el posterior estudio microtermometrico. Previamente a la obtencion de las medidas de temperatura de fusion del hielo (criotermometria) y de homogenizacion del fluido; se realizo un estudio de petrografia de inclusiones fluidas para caracterizarlas y seleccionar las representativas. La interpretacion de los resultados ha permitido confirmar la existencia de un episodio previo de alta temperatura, superior a 282 oC y un evento mineralizador con temperaturas en torno a los 200-240 oC. Sin embargo, las salinidades obtenidas son mucho menores que las previamente publicadas, en torno al 5 % peso eq. NaCl, frente a 11-17 % peso eq. NaCl. Tambien se ha observado un fluido postumo rico en CO2, pero de temperatura superior, en torno a los 270 oC. Los gradientes isotermicos muestran dos focos para dichos fluidos hidrotermales: el primero asociado al stock en la veta Consuelo, y el segundo en la veta Chabuca, asociado a la zona de cabalgamiento de las pizarras de la fm. Chicama sobre las cuarcitas de la fm. Chimu. Este segundo foco puede corresponder con los aportes externos de aguas metamorficas. Para finalizar, se dan una serie de pautas para guiar las futuras exploraciones en el yacimiento. ABSTRACT The Pasto Bueno deposit is located at the northeastern end of the Cordillera Blanca Batholith. It comprises several veins, generally subvertical, associated with the quartz-monzonite stock of Consuzo, dated as Tertiary, which intrudes the Chicama fm. slates and the Chimu fm. quartzites. The main veins trend N-S cutting the stock, although there are also NE-SW and NWSE systems, hosted by the metamorphic rocks. The ore mineralogy comprises wolframite (hubnerite), tetrahedrite/tennantite, sphalerite and galena in a gangue of quartz, fluorite, sericite, pyrite and carbonates, and minor molybdenite, chalcopyrite, bornite, arsenopyrite, enargite (luzonite), stolzite, scheelite, zinnwaldite, topaz, tungstite and native arsenic. Previous studies have characterized Pasto Bueno as a deposit with a large component of greisen, with an evolution of the veins from an early barren 400 oC event , followed by economic mineralization of about 220-250 °C and a late event of 175 -220 oC rich in CO2. Wolframite precipitation occurred from a purely hydrothermal fluid; however, this precipitation was controlled by an external flux of meteoric and/or metamorphic waters. Microthermometric studies of the 3 main structures of the district (Consuelo, Alonso-Fenix and Chabuca veins) have been carried out to depict the evolution of the mineralizing fluid coming from the stock (Consuelo vein) into the metasedimentary rocks of the Chimu and Chicama fm. (Alonso-Fenix and Chabuca veins). The sampling was performed over the main event of mineralization. These samples were subject to a quick plate petrography study in order to select the optimal samples for further microthermometry studies. Before the freezing/heating measures, a fluid inclusion petrography study was done to characterize and select the representative F.I. Interpretation of results has confirmed the existence of a previous episode of higher temperature, over 282 °C, and a mineralizing event with temperatures of about 200-240 °C. However, obtained salinities, about 5 wt% NaCl equivalents, are much lower than those previously reported, about 11-17 wt% NaCl equivalents. A last fluid, rich in CO2, but of higher temperature, about 270 oC, has been characterized. Isothermal gradients show two foci for the hydrothermal fluids: the first one associated to the Consuzo stock as shown in the Consuelo vein, and the second one related to the thrust fault which places the Chicama fm. slates over the Chimu fm. quartzites in the Chabuca vein area. This second focus may correspond to an external input of metamorphic waters. Finally, some guidelines have been given to guide future explorations.

Silica final lens performance in laser fusion facilities: HiPER and LIFE

Nowadays, the projects LIFE (Laser Inertial Fusion Energy) in USA and HiPER (High Power Laser Energy Research) in Europe are the most advanced ones to demonstrate laser fusion energy viability. One of the main points of concern to properly achieve ignition is the performance of the final optics (lenses) under the severe irradiation conditions that take place in fusion facilities. In this paper, we calculate the radiation fluxes and doses as well as the radiation-induced temperature enhancement and colour centre formation in final lenses assuming realistic geometrical configurations for HiPER and LIFE. On these bases, the mechanical stresses generated by the established temperature gradients are evaluated showing that from a mechanical point of view lenses only fulfil specifications if ions resulting from the imploding target are mitigated. The absorption coefficient of the lenses is calculated during reactor startup and steady-state operation. The obtained results reveal the necessity of new solutions to tackle ignition problems during the startup process for HiPER. Finally, we evaluate the effect of temperature gradients on focal length changes and lens surface deformations. In summary, we discuss the capabilities and weak points of silica lenses and propose alternatives to overcome predictable problems

Effect of ion flux on helium retention in helium-irradiated tungsten

Helium retention in irradiated tungsten leads to swelling, pore formation, sample exfoliation and embrittlement with deleterious consequences in many applications. In particular, the use of tungsten in future nuclear fusion plants is proposed due to its good refractory properties. However, serious concerns about tungsten survivability stems from the fact that it must withstand severe irradiation conditions. In magnetic fusion as well as in inertial fusion (particularly with direct drive targets), tungsten components will be exposed to low and high energy ion irradiation (helium), respectively. A common feature is that the most detrimental situations will take place in pulsed mode, i.e., high flux irradiation. There is increasing evidence of a correlation between a high helium flux and an enhancement of detrimental effects on tungsten. Nevertheless, the nature of these effects is not well understood due to the subtleties imposed by the exact temperature profile evolution, ion energy, pulse duration, existence of impurities and simultaneous irradiation with other species. Object Kinetic Monte Carlo is the technique of choice to simulate the evolution of radiation-induced damage inside solids in large temporal and space scales. We have used the recently developed code MMonCa (Modular Monte Carlo simulator), presented at COSIRES 2012 for the first time, to study He retention (and in general defect evolution) in tungsten samples irradiated with high intensity helium pulses. The code simulates the interactions among a large variety of defects and during the irradiation stage and the subsequent annealing steps. The results show that the pulsed mode leads to significantly higher He retention at temperatures higher than 700 K. In this paper we discuss the process of He retention in terms of trap evolution. In addition, we discuss the implications of these findings for inertial fusion.

Effect of Sb incorporation on the electronic structure of InAs quantum dots

On the basis of optical characterization experiments and an eight band kp model, we have studied the effect of Sb incorporation on the electronic structure of InAs quantum dots (QDs). We have found that Sb incorporation in InAs QDs shifts the hole wave function to the center of the QD from the edges of the QD where it is otherwise pinned down by the effects of shear stress. The observed changes in the ground-state energy cannot merely be explained by a composition change upon Sb exposure but can be accounted for when the change in lateral size is taken into consideration. The Sb distribution inside the QDs produces distinctive changes in the density of states, particularly, in the separation between excitation shells. We find a 50% increase in the thermal escape activation energy compared with reference InAs quantum dots as well as an increment of the fundamental transition decay time with Sb incorporation. Furthermore, we find that Sb incorporation into quantum dots is strongly nonlinear with coverage, saturating at low doses. This suggests the existence of a solubility limit of the Sb incorporation into the quantum dots during growth.

Tunable Fano resonance in symmetric multilayered gold nanoshells

We have studied the evolution of dipole–dipole all-plasmonic Fano resonances (FRs) in symmetric multilayered nanoshells as a function of their geometrical parameters. We demonstrate that symmetry breaking is not mandatory for controlling the Fano resonance in such multilayer structures. By carefully selecting the geometrical parameters, the position of the FR can be tuned between 600 and 950 nm and its intensity can be increased up to four fold with respect to the non-optimized structures. Generation of FRs in such symmetric nanostructures presents clear advantages over their asymmetric counterparts, as they are easier to fabricate and can be used in a wider range of technological applications.

MMonCa: An Object Kinetic Monte Carlo simulator for damage irradiation evolution and defect diffusion

In this work, we introduce the Object Kinetic Monte Carlo (OKMC) simulator MMonCa and simulate the defect evolution in three different materials. We start by explaining the theory of OKMC and showing some details of how such theory is implemented by creating generic structures and algorithms in the objects that we want to simulate. Then we successfully reproduce simulated results for defect evolution in iron, silicon and tungsten using our simulator and compare with available experimental data and similar simulations. The comparisons validate MMonCa showing that it is powerful and flexible enough to be customized and used to study the damage evolution of defects in a wide range of solid materials.