A method to obtain the trapping energy and trapping range between hydrogen and defects at finite temperature

Abstract

The binding energy between hydrogen (H) and defects in solid phase materials has been widely studied, which is of vital importance to understand the H retention effects and defect growth mechanisms. However, present studies of binding energy through density functional theory (DFT) or the molecular statics (MS) method were usually performed at 0 K, which could not take the influence of entropy into consideration. In this work, a thermodynamic method has been proposed to obtain the trapping energy between H and defects at finite temperatures. The method is based on the rate theory, which uses trapping energy (V) and trapping range (δ) to describe the trapping properties of defects. Ultimately, a parameterized H spatial cumulative distribution function at thermodynamic equilibrium state could be given. The trapping energy and trapping range parameters in the function can be determined by contrast with the results obtained from molecular dynamics or other methods. This method has been applied to calculate the trapping energies and trapping ranges of H to helium bubble and grain boundary, respectively. Further discussion has been made on the discrepancy between trapping energies obtained by this method and the conventional DFT/MS method.