A fully quantum mechanical calculation of the diffusivity of hydrogen in iron using the tight binding approximation and path integral theory

We present calculations of free energy barriers and diffusivities as functions of temperature for the diffusion of hydrogen in bcc-Fe. This is a fully quantum mechanical approach since the total energy landscape is computed using a new self consistent, transferable tight binding model for interstitial impurities in magnetic iron. Also the hydrogen nucleus is treated quantum mechanically and we compare here two approaches in the literature both based in the Feynman path integral formulation of statistical mechanics. We find that the quantum transition state theory which admits greater freedom for the proton to explore phase space gives result in better agreement with experiment than the alternative which is based on fixed centroid calculations of the free energy barrier. We also find results in better agreement compared to recent centroid molecular dynamics (CMD) calculations of the diffusivity which employed a classical interatomic potential rather than our quantum mechanical tight binding theory. In particular we find first that quantum effects persist to higher temperatures than previously thought, and conversely that the low temperature diffusivity is smaller than predicted in CMD calculations and larger than predicted by classical transition state theory. This will have impact on future modeling and simulation of hydrogen trapping and diffusion

Self-energies in itinerant magnets: A focus on Fe and Ni

We present a detailed study of local and non-local correlations in the electronic structure of elemental transition metals carried out by means of the Quasiparticle Self-consistent GW (QSGW ) and Dynamical Mean Field Theory (DMFT). Recent high resolution ARPES and Haas-van Alphen data of two typical transition metal systems (Fe and Ni) are used as case study. (i) We find that the properties of Fe are very well described by QSGW. Agreement with cyclotron and very clean ARPES measurements is excellent, provided that final-state scattering is taken into account. This establishes the exceptional reliability of QSGW also in metallic systems. (ii) Nonetheless QSGW alone is not able to provide an adequate description of the Ni ARPES data due to strong local spin fluctuations. We surmount this deficiency by combining nonlocal charge fluctuations in QSGW with local spin fluctuations in DMFT (QSGW + 'Magnetic DMFT'). (iii) Finally we show that the dynamics of the local fluctuations are actually not crucial. The addition of an external static field can lead to similarly good results if non-local correlations are included through QSGW

Theoretical evaluation of the role of crystal defects on local equilibrium and effective diffusivity of hydrogen in iron

Hydrogen diffusion and trapping in ferrite is evaluated by quantum mechanically informed kinetic Monte Carlo simulations in defective microstructures. We find that the lattice diffusivity is attenuated by two to four orders of magnitude due to the presence of dislocations. We also find that pipe diffusivity is vanishingly small along screw dislocations and demonstrate that dislocations do not provide fast diffusion pathways for hydrogen as is sometimes supposed. We make contact between our simulations and the predictions of Oriani's theory of ‘effective diffusivity’, and find that local equilibrium is maintained between lattice and trap sites. We also find that the predicted effective diffusivity is in agreement with our simulated results in cases where the distribution of traps is $\textit{spatially homogeneous}$; in the trapping of hydrogen by dislocations where this condition is not met, the Oriani effective diffusivity is in agreement with the simulations to within a factor of two.We are grateful to the European Commission for Funding under the Seventh Framework Programme, Grant No. 263335, MultiHy (multiscale modelling of hydrogen embrittlement in crystalline materials) and Engineering and Physical Sciences Research Council under the HEmS programme grant EP/L014742

Formation of ultracold LiCs molecules

We present the first observation of ultracold LiCs molecules. The molecules are formed in a two-species magneto-optical trap and detected by two-photon ionization and time-of-flight mass spectrometry. The production rate coefficient is found to be in the range 10^{-18}\unit{cm^3s^{-1}} to 10^{-16}\unit{cm^3s^{-1}}, at least an order of magnitude smaller than for other heteronuclear diatomic molecules directly formed in a magneto-optical trap.Comment: 8 pages, 2 figure

Ultracold dense gas of deeply bound heteronuclear molecules

Recently, the quest for an ultracold and dense ensemble of polar molecules has attracted strong interest. Polar molecules have bright prospects for novel quantum gases with long-range and anisotropic interactions, for quantum information science, and for precision measurements. However, high-density clouds of ultracold polar molecules have so far not been produced. Here, we report a key step towards this goal. Starting from an ultracold dense gas of heteronuclear 40K-87Rb Feshbach molecules with typical binding energies of a few hundred kHz and a negligible dipole moment, we coherently transfer these molecules into a vibrational level of the ground-state molecular potential bound by >10 GHz. We thereby increase the binding energy and the expected dipole moment of the 40K-87Rb molecules by more than four orders of magnitude in a single transfer step. Starting with a single initial state prepared with Feshbach association, we achieve a transfer efficiency of 84%. While dipolar effects are not yet observable, the presented technique can be extended to access much more deeply bound vibrational levels and ultimately those exhibiting a significant dipole moment. The preparation of an ultracold quantum gas of polar molecules might therefore come within experimental reach.Comment: 5 pages, 5 figure

Questaal: A package of electronic structure methods based on the linear muffin-tin orbital technique

This paper summarises the theory and functionality behind Questaal, an open-source suite of codes for calculating the electronic structure and related properties of materials from first principles. The formalism of the linearised muffin-tin orbital (LMTO) method is revisited in detail and developed further by the introduction of short-ranged tight-binding basis functions for full-potential calculations. The LMTO method is presented in both Green’s function and wave function formulations for bulk and layered systems. The suite’s full-potential LMTO code uses a sophisticated basis and augmentation method that allows an efficient and precise solution to the band problem at different levels of theory, most importantly density functional theory, LDA +U, quasi-particle self-consistent GW and combinations of these with dynamical mean field theory. This paper details the technical and theoretical bases of these methods, their implementation in Questaal, and provides an overview of the code’s design and capabilities

Questaal: A package of electronic structure methods based on the linear muffin-tin orbital technique

This paper summarises the theory and functionality behind Questaal, an open-source suite of codes for calculating the electronic structure and related properties of materials from first principles. The formalism of the linearised muffin-tin orbital (LMTO) method is revisited in detail and developed further by the introduction of short-ranged tight-binding basis functions for full-potential calculations. The LMTO method is presented in both Green’s function and wave function formulations for bulk and layered systems. The suite’s full-potential LMTO code uses a sophisticated basis and augmentation method that allows an efficient and precise solution to the band problem at different levels of theory, most importantly density functional theory, LDA +U, quasi-particle self-consistent GW and combinations of these with dynamical mean field theory. This paper details the technical and theoretical bases of these methods, their implementation in Questaal, and provides an overview of the code’s design and capabilities