Core physics and safety analysis of Generation-IV Sodium Fast Reactors using existing and newly developed computational tools

El futuro de la energía nuclear de fisión dependerá, entre otros factores, de la capacidad que las nuevas tecnologías demuestren para solventar los principales retos a largo plazo que se plantean. Los principales retos se pueden resumir en los siguientes aspectos: la capacidad de proporcionar una solución final, segura y fiable a los residuos radiactivos; así como dar solución a la limitación de recursos naturales necesarios para alimentar los reactores nucleares; y por último, una mejora robusta en la seguridad de las centrales que en definitiva evite cualquier daño potencial tanto en la población como en el medio ambiente como consecuencia de cualquier escenario imaginable o más allá de lo imaginable. Siguiendo estas motivaciones, la Generación IV de reactores nucleares surge con el compromiso de proporcionar electricidad de forma sostenible, segura, económica y evitando la proliferación de material fisible. Entre los sistemas conceptuales que se consideran para la Gen IV, los reactores rápidos destacan por su capacidad potencial de transmutar actínidos a la vez que permiten una utilización óptima de los recursos naturales. Entre los refrigerantes que se plantean, el sodio parece una de las soluciones más prometedoras. Como consecuencia, esta tesis surgió dentro del marco del proyecto europeo CP-ESFR con el principal objetivo de evaluar la física de núcleo y seguridad de los reactores rápidos refrigerados por sodio, al tiempo que se desarrollaron herramientas apropiadas para dichos análisis. Efectivamente, en una primera parte de la tesis, se abarca el estudio de la física del núcleo de un reactor rápido representativo, incluyendo el análisis detallado de la capacidad de transmutar actínidos minoritarios. Como resultado de dichos análisis, se publicó un artículo en la revista Annals of Nuclear Energy [96]. Por otra parte, a través de un análisis de un hipotético escenario nuclear español, se evalúo la disponibilidad de recursos naturales necesarios en el caso particular de España para alimentar una flota específica de reactores rápidos, siguiendo varios escenarios de demanda, y teniendo en cuenta la capacidad de reproducción de plutonio que tienen estos sistemas. Como resultado de este trabajo también surgió una publicación en otra revista científica de prestigio internacional como es Energy Conversion and Management [97]. Con objeto de realizar esos y otros análisis, se desarrollaron diversos modelos del núcleo del ESFR siguiendo varias configuraciones, y para diferentes códigos. Por otro lado, con objeto de poder realizar análisis de seguridad de reactores rápidos, son necesarias herramientas multidimensionales de alta fidelidad específicas para reactores rápidos. Dichas herramientas deben integrar fenómenos relacionados con la neutrónica y con la termo-hidráulica, entre otros, mediante una aproximación multi-física. Siguiendo este objetivo, se evalúo el código de difusión neutrónica ANDES para su aplicación a reactores rápidos. ANDES es un código de resolución nodal que se encuentra implementado dentro del sistema COBAYA3 y está basado en el método ACMFD. Por lo tanto, el método ACMFD fue sometido a una revisión en profundidad para evaluar su aptitud para la aplicación a reactores rápidos. Durante ese proceso, se identificaron determinadas limitaciones que se discutirán a lo largo de este trabajo, junto con los desarrollos que se han elaborado e implementado para la resolución de dichas dificultades. Por otra parte, se desarrolló satisfactoriamente el acomplamiento del código neutrónico ANDES con un código termo-hidráulico de subcanales llamado SUBCHANFLOW, desarrollado recientemente en el KIT. Como conclusión de esta parte, todos los desarrollos implementados son evaluados y verificados. En paralelo con esos desarrollos, se calcularon para el núcleo del ESFR las secciones eficaces en multigrupos homogeneizadas a nivel nodal, así como otros parámetros neutrónicos, mediante los códigos ERANOS, primero, y SERPENT, después. Dichos parámetros se utilizaron más adelante para realizar cálculos estacionarios con ANDES. Además, como consecuencia de la contribución de la UPM al paquete de seguridad del proyecto CP-ESFR, se calcularon mediante el código SERPENT los parámetros de cinética puntual que se necesitan introducir en los típicos códigos termo-hidráulicos de planta, para estudios de seguridad. En concreto, dichos parámetros sirvieron para el análisis del impacto que tienen los actínidos minoritarios en el comportamiento de transitorios. Concluyendo, la tesis presenta una aproximación sistemática y multidisciplinar aplicada al análisis de seguridad y comportamiento neutrónico de los reactores rápidos de sodio de la Gen-IV, usando herramientas de cálculo existentes y recién desarrolladas ad' hoc para tal aplicación. Se ha empleado una cantidad importante de tiempo en identificar limitaciones de los métodos nodales analíticos en su aplicación en multigrupos a reactores rápidos, y se proponen interesantes soluciones para abordarlas. ABSTRACT The future of nuclear reactors will depend, among other aspects, on the capability to solve the long-term challenges linked to this technology. These are the capability to provide a definite, safe and reliable solution to the nuclear wastes; the limitation of natural resources, needed to fuel the reactors; and last but not least, the improved safety, which would avoid any potential damage on the public and or environment as a consequence of any imaginable and beyond imaginable circumstance. Following these motivations, the IV Generation of nuclear reactors arises, with the aim to provide sustainable, safe, economic and proliferationresistant electricity. Among the systems considered for the Gen IV, fast reactors have a representative role thanks to their potential capacity to transmute actinides together with the optimal usage of natural resources, being the sodium fast reactors the most promising concept. As a consequence, this thesis was born in the framework of the CP-ESFR project with the generic aim of evaluating the core physics and safety of sodium fast reactors, as well as the development of the approppriated tools to perform such analyses. Indeed, in a first part of this thesis work, the main core physics of the representative sodium fast reactor are assessed, including a detailed analysis of the capability to transmute minor actinides. A part of the results obtained have been published in Annals of Nuclear Energy [96]. Moreover, by means of the analysis of a hypothetical Spanish nuclear scenario, the availability of natural resources required to deploy an specific fleet of fast reactor is assessed, taking into account the breeding properties of such systems. This work also led to a publication in Energy Conversion and Management [97]. In order to perform those and other analyses, several models of the ESFR core were created for different codes. On the other hand, in order to perform safety studies of sodium fast reactors, high fidelity multidimensional analysis tools for sodium fast reactors are required. Such tools should integrate neutronic and thermal-hydraulic phenomena in a multi-physics approach. Following this motivation, the neutron diffusion code ANDES is assessed for sodium fast reactor applications. ANDES is the nodal solver implemented inside the multigroup pin-by-pin diffusion COBAYA3 code, and is based on the analytical method ACMFD. Thus, the ACMFD was verified for SFR applications and while doing so, some limitations were encountered, which are discussed through this work. In order to solve those, some new developments are proposed and implemented in ANDES. Moreover, the code was satisfactorily coupled with the thermal-hydraulic code SUBCHANFLOW, recently developed at KIT. Finally, the different implementations are verified. In addition to those developments, the node homogenized multigroup cross sections and other neutron parameters were obtained for the ESFR core using ERANOS and SERPENT codes, and employed afterwards by ANDES to perform steady state calculations. Moreover, as a result of the UPM contribution to the safety package of the CP-ESFR project, the point kinetic parameters required by the typical plant thermal-hydraulic codes were computed for the ESFR core using SERPENT, which final aim was the assessment of the impact of minor actinides in transient behaviour. All in all, the thesis provides a systematic and multi-purpose approach applied to the assessment of safety and performance parameters of Generation-IV SFR, using existing and newly developed analytical tools. An important amount of time was employed in identifying the limitations that the analytical nodal diffusion methods present when applied to fast reactors following a multigroup approach, and interesting solutions are proposed in order to overcome them.

Third-generation synchrotron x-ray diffraction of 6-μm crystal of raite, ≈Na3Mn3Ti0.25Si8O20(OH)2⋅10H2O, opens up new chemistry and physics of low-temperature minerals

The crystal structure of raite was solved and refined from data collected at Beamline Insertion Device 13 at the European Synchrotron Radiation Facility, using a 3 × 3 × 65 μm single crystal. The refined lattice constants of the monoclinic unit cell are a = 15.1(1) Å; b = 17.6(1) Å; c = 5.290(4) Å; β = 100.5(2)°; space group C2/m. The structure, including all reflections, refined to a final R = 0.07. Raite occurs in hyperalkaline rocks from the Kola peninsula, Russia. The structure consists of alternating layers of a hexagonal chicken-wire pattern of 6-membered SiO4 rings. Tetrahedral apices of a chain of Si six-rings, parallel to the c-axis, alternate in pointing up and down. Two six-ring Si layers are connected by edge-sharing octahedral bands of Na+ and Mn3+ also parallel to c. The band consists of the alternation of finite Mn–Mn and Na–Mn–Na chains. As a consequence of the misfit between octahedral and tetrahedral elements, regions of the Si–O layers are arched and form one-dimensional channels bounded by 12 Si tetrahedra and 2 Na octahedra. The channels along the short c-axis in raite are filled by isolated Na(OH,H2O)6 octahedra. The distorted octahedrally coordinated Ti4+ also resides in the channel and provides the weak linkage of these isolated Na octahedra and the mixed octahedral tetrahedral framework. Raite is structurally related to intersilite, palygorskite, sepiolite, and amphibole.

EPCIS Event-Based traceability in pharmaceutical supply chains via automated generation of linked pedigrees

In this paper we show how event processing over semantically annotated streams of events can be exploited, for implementing tracing and tracking of products in supply chains through the automated generation of linked pedigrees. In our abstraction, events are encoded as spatially and temporally oriented named graphs, while linked pedigrees as RDF datasets are their specific compositions. We propose an algorithm that operates over streams of RDF annotated EPCIS events to generate linked pedigrees. We exemplify our approach using the pharmaceuticals supply chain and show how counterfeit detection is an implicit part of our pedigree generation. Our evaluation results show that for fast moving supply chains, smaller window sizes on event streams provide significantly higher efficiency in the generation of pedigrees as well as enable early counterfeit detection.

Machine learning "as a service" for high energy physics (MLaaS4HEP): evolution of a framework for ML-based physics analyses

The scientific success of the LHC experiments at CERN highly depends on the availability of computing resources which efficiently store, process, and analyse the amount of data collected every year. This is ensured by the Worldwide LHC Computing Grid infrastructure that connect computing centres distributed all over the world with high performance network. LHC has an ambitious experimental program for the coming years, which includes large investments and improvements both for the hardware of the detectors and for the software and computing systems, in order to deal with the huge increase in the event rate expected from the High Luminosity LHC (HL-LHC) phase and consequently with the huge amount of data that will be produced. Since few years the role of Artificial Intelligence has become relevant in the High Energy Physics (HEP) world. Machine Learning (ML) and Deep Learning algorithms have been successfully used in many areas of HEP, like online and offline reconstruction programs, detector simulation, object reconstruction, identification, Monte Carlo generation, and surely they will be crucial in the HL-LHC phase. This thesis aims at contributing to a CMS R&D project, regarding a ML "as a Service" solution for HEP needs (MLaaS4HEP). It consists in a data-service able to perform an entire ML pipeline (in terms of reading data, processing data, training ML models, serving predictions) in a completely model-agnostic fashion, directly using ROOT files of arbitrary size from local or distributed data sources. This framework has been updated adding new features in the data preprocessing phase, allowing more flexibility to the user. Since the MLaaS4HEP framework is experiment agnostic, the ATLAS Higgs Boson ML challenge has been chosen as physics use case, with the aim to test MLaaS4HEP and the contribution done with this work.

Populations III.1 and III.2 gamma-ray bursts: constraints on the event rate for future radio and X-ray surveys

Aims. We calculate the theoretical event rate of gamma-ray bursts (GRBs) from the collapse of massive first-generation (Population III; Pop III) stars. The Pop III GRBs could be super-energetic with the isotropic energy up to E(iso) greater than or similar to 10(55-57) erg, providing a unique probe of the high-redshift Universe. Methods. We consider both the so-called Pop III.1 stars (primordial) and Pop III.2 stars (primordial but affected by radiation from other stars). We employ a semi-analytical approach that considers inhomogeneous hydrogen reionization and chemical evolution of the intergalactic medium. Results. We show that Pop III.2 GRBs occur more than 100 times more frequently than Pop III.1 GRBs, and thus should be suitable targets for future GRB missions. Interestingly, our optimistic model predicts an event rate that is already constrained by the current radio transient searches. We expect similar to 10-10(4) radio afterglows above similar to 0.3 mJy on the sky with similar to 1 year variability and mostly without GRBs (orphans), which are detectable by ALMA, EVLA, LOFAR, and SKA, while we expect to observe maximum of N < 20 GRBs per year integrated over at z > 6 for Pop III.2 and N < 0.08 per year integrated over at z > 10 for Pop III.1 with EXIST, and N < 0.2 for Pop III.2 GRBs per year integrated over at z > 6 with Swift.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Lepton sector of a fourth generation

In extensions of the standard model with a heavy fourth generation, one important question is what makes the fourth-generation lepton sector, particularly the neutrinos, so different from the lighter three generations. We study this question in the context of models of electroweak symmetry breaking in warped extra dimensions, where the flavor hierarchy is generated by choosing the localization of the zero-mode fermions in the extra dimension. In this setup the Higgs sector is localized near the infrared brane, whereas the Majorana mass term is localized at the ultraviolet brane. As a result, light neutrinos are almost entirely Majorana particles, whereas the fourth-generation neutrino is mostly a Dirac fermion. We show that it is possible to obtain heavy fourth-generation leptons in regions of parameter space where the light neutrino masses and mixings are compatible with observation. We study the impact of these bounds, as well as the ones from lepton flavor violation, on the phenomenology of these models.

Directed flow at midrapidity in event-by-event hydrodynamics

Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to result in a new type of directed flow (v(1)) that, unlike the usual directed flow, is also present at midrapidity. We compute this new v(1) versus transverse momentum and centrality for Au-Au collisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of v(1) is correlated with the angle of the initial dipole of the distribution, as predicted, though with a large dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with results inferred from STAR correlation data.

Optical second harmonic generation in the centrosymmetric magnetic semiconductors EuTe and EuSe

The magnetic europium chalcogenide semiconductors EuTe and EuSe are investigated by the spectroscopy of second harmonic generation (SHG) in the vicinity of the optical band gap formed by transitions involving the 4f and 5d electronic orbitals of the magnetic Eu(2+) ions. In these materials with centrosymmetric crystal lattice the electric-dipole SHG process is symmetry forbidden so that no signal is observed in zero magnetic field. Signal appears, however, in applied magnetic field with the SHG intensity being proportional to the square of magnetization. The magnetic field and temperature dependencies of the induced SHG allow us to introduce a type of nonlinear optical susceptibility determined by the magnetic-dipole contribution in combination with a spontaneous or induced magnetization. The experimental results can be described qualitatively by a phenomenological model based on a symmetry analysis and are in good quantitative agreement with microscopic model calculations accounting for details of the electronic energy and spin structure.

Theory of near-gap second harmonic generation in centrosymmetric magnetic semiconductors: Europium chalcogenides

Second harmonic generation is strictly forbidden in centrosymmetric materials, within the electric dipole approximation. Recently, it was found that the centrosymmetric magnetic semiconductors EuTe and EuSe can generate near-gap second harmonics, if the system is submitted to an external magnetic field. Here, a theoretical model is presented, which well describes the observed phenomena. The model shows that second harmonic generation becomes efficient when the magnetic dipole oscillations between the band-edge excited states of the system, induced by the excitation light, enter the in-phase regime, which can be achieved by applying a magnetic field to the material.

Spin-Induced Optical Second Harmonic Generation in the Centrosymmetric Magnetic Semiconductors EuTe and EuSe

Spectroscopy of the centrosymmetric magnetic semiconductors EuTe and EuSe reveals spin-induced optical second harmonic generation (SHG) in the band gap vicinity at 2.1-2.4 eV. The magnetic field and temperature dependence demonstrates that the SHG arises from the bulk of the materials due to a novel type of nonlinear optical susceptibility caused by the magnetic dipole contribution combined with spontaneous or induced magnetization. This spin-induced susceptibility opens access to a wide class of centrosymmetric systems by harmonics generation spectroscopy.

Dynamical generation of mass in the D = (2+1) Wess-Zumino model

In this work we study the dynamical generation of mass in the massless N = 1 Wess-Zumino model in a three-dimensional spacetime. Using the tadpole method to compute the effective potential, we observe that supersymmetry is dynamically broken together with the discrete symmetry A(x) -> A(x). We show that this model, different from nonsupersymmetric scalar models, exhibits a consistent perturbative dynamical generation of mass after two-loop corrections to the effective potential.

Femtosecond third-harmonic generation in a glass ceramic containing sodium niobate nanocrystals

The third-harmonic optical susceptibility, chi((3))(3 omega; omega, omega, omega) of a silicate glass ceramic containing sodium niobate nanocrystals was measured for incident broadband light with central frequency omega corresponding to 1900nm. Absolute values of |chi((3))| and the dispersion of the refractive index from 600 to 1900nm were measured using the spectrally resolved femtosecond Maker fringes technique. The experiments show that |chi((3))| is 1 order of magnitude larger than silica, and it grows by similar to 50% when the volume fraction occupied by the nanocrystals increases up to 40%. (C) 2011 Optical Society of America

Generation of nonground-state Bose-Einstein condensates by modulating atomic interactions

A technique is proposed for creating nonground-state Bose-Einstein condensates in a trapping potential by means of the temporal modulation of atomic interactions. Applying a time-dependent spatially homogeneous magnetic field modifies the atomic scattering length. A modulation of the scattering length excites the condensate, which, under special conditions, can be transferred to an excited nonlinear coherent mode. It is shown that a phase-transition-like behavior occurs in the time-averaged population imbalance between the ground and excited states. The application of the technique is analyzed and it is shown that the considered effect can be realized for experimentally available condensates.

The influence of crystallographic orientation on the generation of multiple structural phases generation in Silicon by cyclic microindentation

A Raman scattering study on multiple phase generation in silicon submitted to successive Vickers microindentation cycles, in different crystallographic orientations, was performed. The microindentations were perfon-ned in a virgin single crystal (100)-oriented surface, in the [001] and [011] directions. The results indicated that the formation of multiple phases by cyclic microindentation may depend on the crystallographic direction and number of successive cycles: the onset of several different structural phases was detected after the third cycle for the [001] direction and only after 15 cycles for the [011] direction, indicating that there is a crystallographic orientation dependence for multiple phase generation. (C) 2007 Elsevier B.V. All rights reserved.