Publication: Research › Ph.d. thesis – Annual report year: 2012
The work presented here revolved around assisting in the data analysis by performing density functional theory calculations on the possible dynamical events. For the Mg(BH4)2, in good agreement with the experiments, C2-type rotations occur at lower temperature than C3-type rotations and approximately 15% of the BH4 units activate at a lower temperature than the rest. For the Ca(BH4)2, in addition to the rotational dynamics, an unidentified event was detected which, according to the calculations was most likely due to H2-interstitial defects. In good agreement with the experiments, C3-type rotations activate at lower temperature than C2-type rotations. In order to investigate the environment of reaction pathways, a method for finding the ridge between first order saddle points on a multidimensional surface was developed. Information about the ridge can be used to test the validity of the harmonic approximation to transition state theory for reaction rates, in particular to verify that second order saddle points - maxima along the ridge - are high enough compared to the first order saddle points. Furthermore, corrections to the harmonic approximation can be estimated by direct evaluation of the configuration integral along the ridge. New minima along the ridge can also be identified during the path optimisation, thereby revealing additional transition mechanisms. The method is based on modifying the gradient of a set of points along a path connecting the saddle points to iteratively converge to the ridge. At each iteration during the optimisation, the gradient is inverted along
an unstable eigenmode perpendicular to the path, locally mapping the ridge to a minimum energy path which can be located using various techniques. The method was applied to Al adatom diffusion on the Al(100) surface to find the ridge between 2-, 3- and 4-atom concerted displacement and hop mechanisms for diffusion. Significant corrections were offered for the 3- and 4-atom concerted
displacements. The method offers a simple-to-use way to check the validity of reaction rates but has the potential to offer more accurate rates on its own by representing the transition state with the ridge.
|Number of pages||126|
Note re. dissertation
The work was financially supported by DTU and CAMD, as well as the Nordic Center of Excellence on Hydrogen Storage Materials, the European Graduate School for Sustainable Energy Technologies - The Molecular Approach and Catalysis for Sustainable Energy (CASE).