Search NASA⌕ Search

Engineering topics

Machida, Masahiko

Publications and source records attributed to Machida, Masahiko.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Correlations for the specific heat capacity of ( U x Pu 1 - x ) 1 - y Gd y O 2 - z derived from molecular dynamics

We report UO 2 is the primary conventional fuel used in most nuclear reactors with Gd 2 O 3 commonly added as a burnable absorber to produce a more level power distribution in the reactor core at the beginning of operation. It can also be mixed with other actinide oxides to produce mixed oxide (MOx) fuel. In this study, molecular dynamics simulations were used to predict the specific heat capacity of Gd-doped PuO 2 , UO 2 and (U, Pu)O 2 MOx accommodating Gd 3+ substituted at cation sites via two charge compensation mechanisms - oxygen vacancy formation and the oxidation of U 4+ to U 5+ . The specific heat capacity values for PuO 2 and UO 2 are in good agreement with other studies showing a distinct peak at high temperatures - above 1800 K. As Gd 3+ is added, the peak height reduces for each composition considered. An analytical fit was applied to the data where Gd 3+ was fully charge compensated by either oxygen vacancies or U 5+ . The expression was then validated by predicting the specific heat capacity for three compositions of (Ux Pu 1-x ) 1-y Gd y O 2-z containing both oxygen vacancies and U 5+ , and compared to molecular dynamics data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Radiocesium interaction with clay minerals: Theory and simulation advances Post–Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Transmutation effects on long-term Cs retention in phyllosilicate minerals from first principles

The accidental release and incorporation of radiocesium into soil minerals represents a massive environmental, technical and social challenge. Accurately forecasting the evolving distribution and fate of long- and medium-lived isotopes such as 137Cs and 134Cs over decadal time scales is essential. The cesium cation has long been modeled as a strongly and selectively sorbed species into clay mineral interlayers; however, because of the time scales involved by the radioisotopes half-lives, the effects of radioactive decay on Cs retention have been unknown. We report density functional theory (DFT) simulations of transmutation effects of radiocesium on long-term Cs retention in phlogopite. The calculations show that the progressive appearance of daughter product Ba2+ is accompanied by a proportional increase in thermodynamic driving force to preferentially discharge remaining Cs, both radioactive and stable, back into aqueous solution. Based on thermodynamic analysis, the findings indicate that radiocesium transmutation provides a mean to weaken the binding of Cs in phyllosilicate minerals, therefore potentially involving a premature re-release of Cs back into the environment. In the case where radiogenic Ba2+ ions accumulate in the mineral, collateral effects would ultimately be an increase in the overall interlayer binding energy and a lower resorption capacity.

Sassi, Michel JPC↗
DOE OSTI2017

Origin of 6-fold coordinated aluminum at (010)-type pyrophyllite edges

To better understand the aqueous chemical reactivity of clay mineral edges we explored the relationships between hydration and the structure of (010)-type edges of pyrophyllite. In particular, we used density functional theory and the quantum theory of atoms in molecules to evaluate the stability of 6-fold coordinated hydrated aluminum at the edge in terms of the electron density distribution. Geometrical optimization revealed an intra-edge hydrogen bond network between aluminol hydroxyls and water ligands completing the aluminum coordination shell. From the electron density isosurfaces one water ligand is not covalently bonded to aluminum. Bader charge analysis revealed that there is a small charge transfer between this water ligand and the aluminum that stabilizes its location, but which derives from the formation of the hydrogen bond network and not from Al-OH 2 interaction per se. To clarify which edge species make the most significant hydrogen bond to this water ligand, we evaluate Bader charges in edge systems including 5-fold coordinated aluminums with and without geometrical optimization. From this comparison the charge transfer is a combined result of charge transfer from the aluminum to aluminol hydroxyls and then from aluminol hydroxyls to OH 2 ligands. This suggests that the weakly bound water ligand in question, and more generally 6-fold hydrated edge Al coordination, is stabilized primarily by the hydrogen bond network which in turn leads to weak ionic attraction to the aluminum center itself. The finding highlights the importance of cooperative effects between solvent structure and the coordination of metal cations exposed at clay mineral edges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Molecular dynamics simulations of cesium adsorption on illite nanoparticles

Here we present that the charged surfaces of micaceous minerals, especially illite, regulate the mobility of the major radioisotopes of Cs ( 134 Cs, 135 Cs, 137 Cs) in the geosphere. Despite the long history of Cs adsorption studies, the nature of the illite surface sites remains incompletely understood. To address this problem, we present atomistic simulations of Cs competition with Na for three candidate illite adsorption sites – edge, basal plane, and interlayer. Our simulation results are broadly consistent with affinities and selectivities that have been inferred from surface complexation models. Cation exchange on the basal planes is thermodynamically ideal, but exchange on edge surfaces and within interlayers shows complex, thermodynamically non-ideal behavior. The basal planes are weakly Cs-selective, while edges and interlayers have much higher affinity for Cs. The dynamics of Na-Cs exchange are rapid for both cations on the basal planes, but considerably slower for Cs localized on edge surfaces. Lastly, in addition to new insights into Cs adsorption and exchange with Na on illite, we report the development of a methodology capable of simulating fully-flexible clay mineral nanoparticles with stable edge surfaces using a well-tested interatomic potential model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Molecular Simulation of Cesium Adsorption at the Basal Surface of Phyllosilicate Minerals

A better understanding of the thermodynamics of radioactive cesium uptake at the surfaces of phyllosilicate minerals is needed to understand mechanisms of its selective adsorption and help guide the development of practical and inexpensive decontamination techniques. In this work, molecular dynamics simulations were carried out to determine the thermodynamics of adsorption of Cs + at the basal surface of six 2:1 phyllosilicate minerals, namely pyrophyllite, illite, muscovite, phlogopite, celadonite, and margarite. These minerals were selected to isolate the effects of the magnitude of the permanent layer charge (≤ 2), its location (tetrahedral versus octahedral sheet), and the structure of the octahedral sheet (dioctahedral versus trioctahedral). Good agreement was obtained with experiment in terms of the hydration free energy of Cs + and the structure and thermodynamics of Cs + adsorption at the muscovite basal surface, for which published data were available for comparison. With the exception of pyrophyllite, which did not exhibit an inner-sphere free energy minimum, all phyllosilicate minerals showed similar behavior with respect to Cs + adsorption; notably, Cs + adsorption was predominantly inner-sphere whereas outer-sphere adsorption was very weak with the simulations predicting the formation of an extended outer-sphere complex. For a given location of the layer charge, the free energy of adsorption as an inner-sphere complex was found to vary linearly with the magnitude of the layer charge. For a given location and magnitude of the layer charge, adsorption at phlogopite (trioctahedral sheet structure) was much less favorable than at muscovite (dioctahedral sheet structure) due to the electrostatic repulsion between the adsorbed Cs + and the hydrogen atom of the hydroxyl group directly below the six-membered siloxane ring cavity. For a given magnitude of the layer charge and structure of the octahedral sheet, adsorption at celadonite (layer charge located in the octahedral sheet) was favored over muscovite (layer charge located in the tetrahedral sheet) due to the increased distance with surface potassium ions.

(001) surface↗
DOE OSTI2016

Radiation-Damage Resistance in Phyllosilicate Minerals From First Principles and Implications for Radiocesium and Strontium Retention in Soils

Accidental discharges of the hazardous nuclear fission products 137Cs+ and 90Sr2+ into the environment, such as during the Fukushima Dai-ichi nuclear accident, have repeatedly occurred throughout the nuclear age. Numerous studies of their fate and transport in soils and sediments have demonstrated their strong and selective binding to phyllosilicates such as clay minerals, primarily via cation exchange into interlayer sites. The locally concentrated amounts of these radioactive beta-emitters that can be found in these host minerals raises important questions regarding the long-term interplay and durability of radioisotope-clay association, which is not well known. This study goes beyond the usual short-term focus to address the permanence of radioisotope retention in clay minerals, by developing a general theoretical understanding of their resistance to defect creation. We report ab initio molecular dynamics (AIMD) calculations of the threshold displacement energy (TDE) of each symmetry-unique atomic specie comprising the unit cell of model vermiculite. The determined TDE values are material specific and radiation independent. We use them to estimate the probability of Frenkel pair creation via direct electron-ion collision, as could be induced by the passage of a high energy electron emitted during the beta-decay of 137Cs, 90Sr and daughter 90Y. For 137Cs and 90Sr, we found that the probability is about 36%, while for 90Y the probability is much higher with about 89%. The long-term retention picture that emerges is that decay will progressively alter the clay interlayer structure and charge, likely leading to delamination of the clay and re-release of residual parent isotopes. Future work examining the effect of Frenkel defect accumulation on the binding energy of parent and daughter radionuclides in the interlayer is thus justifiable, and potentially important for accurate long-term forecasting of radionuclide transport in the environment.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2011

Spectral properties of trapped one-dimensional ultracold fermions loaded on optical lattices

We examine spectral properties on one-dimensional (1D) trapped ultracold Fermi atoms loaded on optical lattices by using the dynamical density-matrix renormalization group (DDMRG) method. We find that spectra are very rich due to the interplay of free or repulsive interaction and the harmonic-trap potential. One of the rich examples is that a 1D Tomonaga-Luttinger (TL) liquid spectrum emerges beneath multiple bound-state flat levels in the Mott-plateau phase and the TL spectrum gradually disappears with a partial breakdown of the central Mott plateau by increasing the trap strength. A more striking example is the growth of a distinct low-energy broadened band from the bound-state flat levels together with the breakdown. It demonstrates a carrier doping effect into the Mott phase in optical lattice systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Collective modes in multiband superfluids and superconductors: Multiple dynamical classes

One important way to characterize the states having a gauge symmetry spontaneously broken over multibands is to look at their collective excitation modes. We find that a three-band system has multiple Leggett modes with significantly different masses, which can be classified into different dynamical classes according to whether multiple interband Josephson currents add or cancel. This provides a way to dynamically characterize multiband superconductivity while the pairing symmetry is a static property.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗