Comparison among MCNP-based depletion codes applied to burnup calculations of pebble-bed HTR lattices

The double-heterogeneity characterising pebble-bed high temperature reactors (HTRs) makes Monte Carlo based calculation tools the most suitable for detailed core analyses. These codes can be successfully used to predict the isotopic evolution during irradiation of the fuel of this kind of cores. At the moment, there are many computational systems based on MCNP that are available for performing depletion calculation. All these systems use MCNP to supply problem dependent fluxes and/or microscopic cross sections to the depletion module. This latter then calculates the isotopic evolution of the fuel resolving Bateman's equations. In this paper, a comparative analysis of three different MCNP-based depletion codes is performed: Montburns2.0, MCNPX2.6.0 and BGCore. Monteburns code can be considered as the reference code for HTR calculations, since it has been already verified during HTR-N and HTR-N1 EU project. All calculations have been performed on a reference model representing an infinite lattice of thorium-plutonium fuelled pebbles. The evolution of k-inf as a function of burnup has been compared, as well as the inventory of the important actinides. The k-inf comparison among the codes shows a good agreement during the entire burnup history with the maximum difference lower than 1%. The actinide inventory prediction agrees well. However significant discrepancy in Am and Cm concentrations calculated by MCNPX as compared to those of Monteburns and BGCore has been observed. This is mainly due to different Am-241 (n,γ) branching ratio utilized by the codes. The important advantage of BGCore is its significantly lower execution time required to perform considered depletion calculations. While providing reasonably accurate results BGCore runs depletion problem about two times faster than Monteburns and two to five times faster than MCNPX. © 2009 Elsevier B.V. All rights reserved.

Efficient generation of one-group cross sections for coupled Monte Carlo depletion calculations

Coupled Monte Carlo depletion systems provide a versatile and an accurate tool for analyzing advanced thermal and fast reactor designs for a variety of fuel compositions and geometries. The main drawback of Monte Carlo-based systems is a long calculation time imposing significant restrictions on the complexity and amount of design-oriented calculations. This paper presents an alternative approach to interfacing the Monte Carlo and depletion modules aimed at addressing this problem. The main idea is to calculate the one-group cross sections for all relevant isotopes required by the depletion module in a separate module external to Monte Carlo calculations. Thus, the Monte Carlo module will produce the criticality and neutron spectrum only, without tallying of the individual isotope reaction rates. The onegroup cross section for all isotopes will be generated in a separate module by collapsing a universal multigroup (MG) cross-section library using the Monte Carlo calculated flux. Here, the term "universal" means that a single MG cross-section set will be applicable for all reactor systems and is independent of reactor characteristics such as a neutron spectrum; fuel composition; and fuel cell, assembly, and core geometries. This approach was originally proposed by Haeck et al. and implemented in the ALEPH code. Implementation of the proposed approach to Monte Carlo burnup interfacing was carried out through the BGCORE system. One-group cross sections generated by the BGCORE system were compared with those tallied directly by the MCNP code. Analysis of this comparison was carried out and led to the conclusion that in order to achieve the accuracy required for a reliable core and fuel cycle analysis, accounting for the background cross section (σ0) in the unresolved resonance energy region is essential. An extension of the one-group cross-section generation model was implemented and tested by tabulating and interpolating by a simplified σ0 model. A significant improvement of the one-group cross-section accuracy was demonstrated.

The stochastic implicit Euler method - A stable coupling scheme for Monte Carlo burnup calculations

Existing Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the coupling schemes of the existing Monte Carlo burnup codes can be numerically unstable. Here we develop the Stochastic Implicit Euler method - a stable and efficient new coupling scheme. The implicit solution is obtained by the stochastic approximation at each time step. Our test calculations demonstrate that the Stochastic Implicit Euler method can provide an accurate solution to problems where the methods in the existing Monte Carlo burnup codes fail. © 2013 Elsevier Ltd. All rights reserved.

Defects in gallium nitride nanowires: First principles calculations

Atomic configurations and formation energies of native defects in an unsaturated GaN nanowire grown along the [001] direction and with (100) lateral facets are studied using large-scale ab initio calculation. Cation and anion vacancies, antisites, and interstitials in the neutral charge state are all considered. The configurations of these defects in the core region and outermost surface region of the nanowire are different. The atomic configurations of the defects in the core region are same as those in the bulk GaN, and the formation energy is large. The defects at the surface show different atomic configurations with low formation energy. Starting from a Ga vacancy at the edge of the side plane of the nanowire, a N-N split interstitial is formed after relaxation. As a N site is replaced by a Ga atom in the suboutermost layer, the Ga atom will be expelled out of the outermost layers and leaves a vacancy at the original N site. The Ga interstitial at the outmost surface will diffuse out by interstitialcy mechanism. For all the tested cases N-N split interstitials are easily formed with low formation energy in the nanowires, indicating N-2 molecular will appear in the GaN nanowire, which agrees well with experimental findings.

Magnetic coupling properties of mn-doped ZnO nanowires: First-principles calculations

Based on the density functional theory, we study the magnetic coupling properties of Mn-doped ZnO nanowires. For the nanowires with passivated surfaces, the antiferromagnetic state is found and the Mn atoms have a clustering tendency. When the distance between two Mn atoms is large, the system energetically favors the paramagnetic or spin-glass state. For the nanowires with unpassivated surfaces, the ferromagnetic (FM) coupling states appear between the two nearest Mn atoms, and the zinc vacancies can further stabilize the FM states between them. The electrons with enough concentration possibly mediate the FM coupling due to the negative exchange splitting of conduction band minimum induced by the s-d coupling, which could be useful in nanomaterial design for spintronics. (C) 2008 American Institute of Physics.

Electronic structures and mechanical properties of uranium monocarbide from first-principles LDA plus U and GGA plus U calculations

The electronic structure, elastic constants, Poisson's ratio, and phonon dispersion curves of UC have been systematically investigated from the first-principles calculations by the projector-augmented-wave (PAW) method. In order to describe precisely the strong on-site Coulomb repulsion among the localized U 5f electrons, we adopt the local density approximation (LDA) + U and generalized gradient approximation (GGA) + U formalisms for the exchange correlation term. We systematically study how the electronic properties and elastic constants of UC are affected by the different choice of U as well as the exchange-correlation potential. We show that by choosing an appropriate Hubbard U parameter within the GGA + U approach, most of our calculated results are in good agreement with the experimental data. Therefore. the results obtained by the GGA + U with effective Hubbard parameter U chosen around 3 eV for UC are considered to be reasonable. (C) 2009 Elsevier B.V. All rights reserved.

Magnetic coupling properties of rare-earth metals (Gd, Nd) doped ZnO: First-principles calculations

The electronic structure and magnetic coupling properties of rare-earth metals (Gd, Nd) doped ZnO have been investigated using first-principles methods. We show that the magnetic coupling between Gd or Nd ions in the nearest neighbor sites is ferromagnetic. The stability of the ferromagnetic coupling between Gd ions can be enhanced by appropriate electron doping into ZnO Gd system and the room-temperature ferromagnetism can be achieved. However, for ZnO Nd system, the ferromagnetism between Nd ions can be enhanced by appropriate holes doping into the sample. The room-temperature ferromagnetism can also be achieved in the n-conducting ZnO Nd sample. Our calculated results are in good agreement with the conclusions of the recent experiments. The effect of native defects (V-Zn, V-O) on the ferromagnetism is also discussed. (C) 2009 American Institute of Physics. [DOI 10.1063/1.3176490]

Deep levels in high resistivity GaN detected by thermally stimulated luminescence and first-principles calculations

Thermally stimulated luminescence spectroscopy has been applied to study the deep centres in unintentionally doped high resistivity GaN epilayers grown by the metal organic chemical vapour deposition method on c-sapphire substrates. Two trap states with activation energies of 0.12 and 0.62 eV are evaluated from two luminescence peaks at 141.9 and 294.7 K in the luminescence curve. Our spectroscopy measurement, in combination with more accurate first-principles studies, provided insights into the microscopic origin of these levels. Our investigations suggest that the lower level at 0.12 eV might originate from C-N, which behaves as a hole trap state; the deeper level at 0.62 eV can be correlated with V-Ga that corresponds to the yellow luminescence band observed in low-temperature photoluminescence spectra.

Variational calculations of neutral bound excitons in GaAs quantum-well wires

The binding energy of an exciton bound to a neutral donor (D-0,X) in GaAs quantum-well wires is calculated variationally as a function of the wire width for different positions of the impurity inside the wire by using a two-parameter wavefunction. There is no artificial parameter added in our calculation. The results we have obtained show that the binding energies are closely correlated to the sizes of the wire, the impurity position, and also that their magnitudes are greater than those in the two-dimensional quantum wells compared. In addition, we also calculate the average interparticle distance as a function of the wire width. The results are discussed in detail.

Effects of magnetic field on the electronic structure of wurtzite quantum dots: Calculations using effective-mass envelope function theory

The Hamiltonian of the wurtzite quantum dots in the presence of an external homogeneous magnetic field is given. The electronic structure and optical properties are studied in the framework of effective-mass envelope function theory. The energy levels have new characteristics, such as parabolic property, antisymmtric splitting, and so on, different from the Zeeman splitting. With the crystal field splitting energy Delta(c)=25 meV, the dark excitons appear when the radius is smaller than 25.85 A in the absence of external magnetic field. This result is more consistent with the experimental results reported by Efros [Phys. Rev. B 54, 4843 (1996)]. It is found that dark excitons become bright under appropriate magnetic field depending on the radius of dots. The circular polarization factors of the optical transitions of randomly oriented dots are zero in the absence of external magnetic field and increase with the increase of magnetic field, in agreement with the experimental results. The circular polarization factors of single dots change from nearly 0 to about 1 as the orientation of the magnetic field changes from the x axis of the crystal structure to the z axis, which can be used to determine the orientation of the z axis of the crystal structure of individual dots. The antisymmetric Hamiltonian is very important to the effects of magnetic field on the circular polarization of the optical transition of quantum dots.

Giant Zeeman splitting in manganese-doped ZnSe quantum spheres: Eight-band effective-mass calculations

We deduce the eight-band effective-mass Hamiltonian model for a manganese-doped ZnSe quantum sphere in the presence of the magnetic field, including the interaction between the conduction and valence bands, the spin-orbit coupling within the valence bands, the intrinsic spin Zeeman splitting, and the sp-d exchange interaction between the carriers and magnetic ion in the mean-field approximation. The size dependence of the electron and hole energy levels as well as the giant Zeeman splitting energies are studied theoretically. We find that the hole giant Zeeman splitting energies decrease with the increasing radius, smaller than that in the bulk material, and are different for different J(z) states, which are caused by the quantum confinement effect. Because the quantum sphere restrains the excited Landau states and exciton states, in the experiments we can observe directly the Zeeman splitting of basic states. At low magnetic field, the total Zeeman splitting energy increases linearly with the increasing magnetic field and saturates at modest field which is in agreement with recent experimental results. Comparing to the undoped case, the Zeeman splitting energy is 445 times larger which provides us with wide freedom to tailor the electronic structure of DMS nanocrystals for technological applications.