Thermal neutron scattering law calculations using ab initio molecular dynamics

Abstract

Abstract. In recent years, methods for the calculation of the thermal scattering law (i.e. S(α, β), where α and β are dimensionless momentum and energy transfer variables, respectively) were developed based on ab initio lattice dynamics (AILD) and/or classical molecular dynamics (CMD). While these methods are now mature and efficient, further advancement in the application of such atomistic techniques is possible using ab initio molecular dynamics (AIMD) methods. In this case, temperature effects are inherently included in the calculation, e.g. phonon density of states (DOS), while using ab initio force fields that eliminate the need for parameterized semi-empirical force fields. In this work, AIMD simulations were performed to predict the phonon spectra as a function of temperature for beryllium and graphite, which are representative nuclear reactor moderator and reflector materials. Subsequently, the calculated phonon spectra were utilized to predict S(α, β) using the LEAPR module of the NJOY code. The AIMD models of beryllium and graphite were 5 × 5 × 5 crystal unit cells (250 atoms and 500 atoms respectively). Electronic structure calculations for the prediction of Hellman-Feynman forces were performed using density functional theory with a GGA exchange correlation functional and corresponding core electron pseudopotentials. AIMD simulations of 1000-10,000 time-steps were performed with the canonical ensemble (NVT thermostat) for several temperatures between 300 K and 900 K. The phonon DOS were calculated as the power spectrum of the AIMD predicted velocity autocorrelation functions. The resulting AIMD phonon DOS and corresponding inelastic thermal neutron scattering cross sections at 300 K, where anharmonic effects are expected to be small, were found to be in reasonable agreement with the results generated using traditional AILD. This illustrated the validity of the AIMD approach. However, since the impact of the temperature on the phonon DOS (e.g. broadening of spectral peaks) was observed in AIMD analysis, this technique may be envisioned as the approach for deriving the needed atomistic data for thermal scattering law calculations under realistic temperature and structural conditions for a given material.

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures With Diffusion Experiments

“…This method was selected because PNDA simulations with MC21 were previously found to be highly successful in validating a theoretically-generated TSL for ice Ih over a wide temperature range [60]. Three ENDF light-water TSL models for H bound in H2O (H-H2O) are tested in this work: ENDF/B-VII.1 [61], ENDF/B-VIII.0 [8], and a version under development by North Carolina State University (NCSU) based on MD simulations [62,63].…”

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

“…The forces describing the small displacements of zirconium atoms are expected to be well-characterized within the harmonic approximation; however, hydrogen-due to its small mass and large vibrational displacement-may have significant anharmonic contributions to its motion even at 0 K (i.e., zero-point motion). AIMD simulations, which inherently include anharmonicity, were therefore performed using methods developed in previous work [20,21] to generate the hydrogen DOS. Using a reduced 350 eV planewave cutoff and a 1 × 1 × 1 k-mesh, a 3 × 3 × 3 (270 atoms) δ-ZrH 1.5 and 4 × 4 × 3 (288 atoms) -ZrH 2 supercell were first allowed to thermodynamically equilibrate for 1 ps at 0.00025 ps time-steps under a 300 K NVT thermostat.…”

The correction of Inelastic Neutron Scattering data of organic samples using the Average Functional Group Approximation

“…For this reason, experimentally obtained and tabulated TCS are available for just few systems. In order to overcome the impossibility of direct measurements, it is possible to retrieve TCS by calculations from simplified models [6,7], from molecular dynamics [8], as well as from ab initio calculations [9][10][11][12][13].…”