Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition.

Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.

Comprehensive Studies of Magnetic Transitions and Spin-Phonon Couplings in the Tetrahedral Cobalt Complex Co(AsPh3)2I2.

A combination of inelastic neutron scattering (INS), far-IR magneto-spectroscopy (FIRMS), and Raman magneto-spectroscopy (RaMS) has been used to comprehensively probe magnetic excitations in Co(AsPh3)2I2 (1), a reported single-molecule magnet (SMM). With applied field, the magnetic zero-field splitting (ZFS) peak (2D') shifts to higher energies in each spectroscopy. INS placed the ZFS peak at 54 cm-1, as revealed by both variable-temperature (VT) and variable-magnetic-field data, giving results that agree well with those from both far-IR and Raman studies. Both FIRMS and RaMS also reveal the presence of multiple spin-phonon couplings as avoided crossings with neighboring phonons. Here, phonons refer to both intramolecular and lattice vibrations. The results constitute a rare case in which the spin-phonon couplings are observed with both Raman-active (g modes) and far-IR-active phonons (u modes; space group P21/c, no. 14, Z = 4 for 1). These couplings are fit using a simple avoided crossing model with coupling constants of ca. 1-2 cm-1. The combined spectroscopies accurately determine the magnetic excited level and the interaction of the magnetic excitation with phonon modes. Density functional theory (DFT) phonon calculations compare well with INS, allowing for the assignment of the modes and their symmetries. Electronic calculations elucidate the nature of ZFS in the complex. Features of different techniques to determine ZFS and other spin-Hamiltonian parameters in transition-metal complexes are summarized.

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories.

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. The results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

Calculation of the Thermal Neutron Scattering Cross-Section of Solids Using OCLIMAX.

The thermal neutron scattering cross-section of a solid depends on the energy (or wavelength) of the incident neutrons. Devising a method to calculate the energy dependence from first principles, without the approximations built in the scattering theory, has been a major undertaking in nuclear engineering. Here, we demonstrate such a calculation method using the program OCLIMAX. Our approach eliminates various approximations and limitations involved in a regular calculation with the LEAPR module of NJOY code, and the results are compared with available experimental and theoretical data. It is also demonstrated how additional insight can be obtained from the calculated full dynamical structure factor. The results reported here show the great potential and excellent platform provided by OCLIMAX for future development in the study of neutron thermalization in solid materials for different applications.

Simulation of Inelastic Neutron Scattering Spectra Using OCLIMAX.

Studying the vibration of atoms is of fundamental importance and can provide critical insight for the understanding of materials behavior, such as structure and phase transition, thermodynamics, and chemical reactions. The atomic vibration can be probed using vibrational spectroscopy with various incident particles such as photons, neutrons, or electrons. A major challenge when applying these techniques is often how to interpret the vibrational spectra and how to make connections to the theory. To this end, methods that can simulate the spectra from atomistic models are highly desired. In this paper, we present a program developed for the simulation of inelastic neutron scattering spectra. It has many new and useful features that were not previously available and will greatly facilitate the analysis and understanding of inelastic neutron scattering data.

Effect of magnetic fields on the methyl rotation in a paramagnetic cobalt(ii) complex. Quasielastic neutron scattering studies.

Molecular dynamics is a fundamental property of metal complexes. These dynamic processes, especially for paramagnetic complexes under external magnetic fields, are in general not well understood. Quasielastic neutron scattering (QENS) in 0-4 T magnetic fields has been used to study the dynamics of Co(acac)2(D2O)2 (1-d4, acac = acetylacetonate). At 80-100 K, rotation of the methyl groups on the acac ligands is the dominant dynamical process. This rotation is slowed down by the magnetic field increase. Rotation times at 80 K are 5.6(3) × 10-10 s at 0 T and 2.04(10) × 10-9 s at 4 T. The QENS studies suggest that methyl groups in these paramagnetic Co(ii) molecules do not behave as isolated units, which is consistent with results from earlier magnetic susceptibility studies indicating the presence of intermolecular interactions. DFT calculations show that unpaired electron spin density in 1 is dispersed to the atoms of both acac and H2O ligands. Methyl torsions in 1-d4 have also been observed at 5-100 K in inelastic neutron spectroscopy (INS). The QENS and INS results here help understand the dynamics of the compound in the solid state.

Spin-phonon couplings in transition metal complexes with slow magnetic relaxation.

Spin-phonon coupling plays an important role in single-molecule magnets and molecular qubits. However, there have been few detailed studies of its nature. Here, we show for the first time distinct couplings of g phonons of CoII(acac)2(H2O)2 (acac = acetylacetonate) and its deuterated analogs with zero-field-split, excited magnetic/spin levels (Kramers doublet (KD)) of the S = 3/2 electronic ground state. The couplings are observed as avoided crossings in magnetic-field-dependent Raman spectra with coupling constants of 1-2 cm-1. Far-IR spectra reveal the magnetic-dipole-allowed, inter-KD transition, shifting to higher energy with increasing field. Density functional theory calculations are used to rationalize energies and symmetries of the phonons. A vibronic coupling model, supported by electronic structure calculations, is proposed to rationalize the behavior of the coupled Raman peaks. This work spectroscopically reveals and quantitates the spin-phonon couplings in typical transition metal complexes and sheds light on the origin of the spin-phonon entanglement.

Capture of heavy hydrogen isotopes in a metal-organic framework with active Cu(I) sites.

The production of pure deuterium and the removal of tritium from nuclear waste are the key challenges in separation of light isotopes. Presently, the technological methods are extremely energy- and cost-intensive. Here we report the capture of heavy hydrogen isotopes from hydrogen gas by selective adsorption at Cu(I) sites in a metal-organic framework. At the strongly binding Cu(I) sites (32 kJ mol-1) nuclear quantum effects result in higher adsorption enthalpies of heavier isotopes. The capture mechanism takes place most efficiently at temperatures above 80 K, when an isotope exchange allows the preferential adsorption of heavy isotopologues from the gas phase. Large difference in adsorption enthalpy of 2.5 kJ mol-1 between D2 and H2 results in D2-over-H2 selectivity of 11 at 100 K, to the best of our knowledge the largest value known to date. Combination of thermal desorption spectroscopy, Raman measurements, inelastic neutron scattering and first principles calculations for H2/D2 mixtures allows the prediction of selectivities for tritium-containing isotopologues.

Hydrogen release reactions of Al-based complex hydrides enhanced by vibrational dynamics and valences of metal cations.

Hydrogen release from Al-based complex hydrides composed of metal cation(s) and [AlH4]- was investigated using inelastic neutron scattering viewed from vibrational dynamics. The hydrogen release followed the softening of translational and [AlH4]- librational modes, which was enhanced by vibrational dynamics and the valence(s) of the metal cation(s).

Hydrogen adsorption on two catalysts for the ortho- to parahydrogen conversion: Cr-doped silica and ferric oxide gel.

Molecular hydrogen exists in two spin-rotation coupled states: parahydrogen and orthohydrogen. Due to the variation of energy with rotational level, the occupation of ortho- and parahydrogen states is temperature dependent, with parahydrogen being the dominant species at low temperatures. The equilibrium at 20 K (99.8% parahydrogen) can be reached by natural conversion only after a lengthy process. With the use of a suitable catalyst, this process can be shortened significantly. Two types of commercial catalysts currently being used for ortho- to parahydrogen conversion are: iron(iii) oxide (Fe2O3, IONEX®), and chromium(ii) oxide doped silica catalyst (CrO·SiO2, OXISORB®). We investigate the interaction of ortho- and parahydrogen with the surfaces of these ortho-para conversion catalysts using neutron vibrational spectroscopy. The catalytic surfaces have been characterized using X-ray absorption fine structure (XAFS) and X-ray/neutron pair distribution function measurements.

Gate-opening effect in ZIF-8: the first experimental proof using inelastic neutron scattering.

The gate-opening phenomenon in ZIFs is of paramount importance to understand their behavior in industrial molecular separations. Here we show for the first time using in situ inelastic neutron scattering (INS) the swinging of the -CH3 groups and the imidazolate linkers in the prototypical ZIF-8 and ZIF-8@AC hybrid materials upon exposure to mild N2 pressure.

The free-energy barrier to hydride transfer across a dipalladium complex.

We use density-functional theory molecular dynamics (DFT-MD) simulations to determine the hydride transfer coordinate between palladium centres of the crystallographically observed terminal hydride locations, Pd-Pd-H, originally postulated for the solution dynamics of the complex bis-NHC dipalladium hydride [{(MesIm)2CH2}2Pd2H][PF6], and then calculate the free-energy along this coordinate. We estimate the transfer barrier-height to be about 20 kcal mol(-1) with a hydride transfer rate in the order of seconds at room temperature. We validate our DFT-MD modelling using inelastic neutron scattering which reveals anharmonicity of the hydride environment that is so pronounced that there is complete failure of the harmonic model for the hydride ligand. The simulations are extended to high temperature to bring the H-transfer to a rate that is accessible to the simulation technique.

Low energy structural dynamics and constrained libration of Li(NH3)4, the lowest melting point metal.

The lattice and molecular dynamics for the solid phases of the lowest melting-point metal, Li(NH3)4, are determined by incoherent inelastic neutron scattering. Measurements of internal molecular displacements and distortions of the Li(NH3)4 units have been modelled and assigned using density functional theory calculations for the solid and molecular system. Inelastic neutron scattering measurement allow for the first determination of NH3 librational transitions.

Proton vibrational dynamics in lithium imide investigated through incoherent inelastic and Compton neutron scattering.

In the present study we report neutron spectroscopic measurements on polycrystalline lithium imide, namely, incoherent inelastic neutron scattering at 20 K, and neutron Compton scattering from 10 K up to room temperature. From the former technique the H-projected density of phonon states up to 100 meV is derived, while the latter works out the spherically averaged single-particle (i.e., H, Li, and N) momentum distributions and, from this, the mean kinetic energies. Only for H at the lowest investigated temperature, non-gaussian components of its momentum distribution are detected. However, these components do not seem directly connected to the system anharmonicity, being fully compatible with the simple N-H bond anisotropy. Neutron data are also complemented by ab initio lattice dynamics simulations, both harmonic and, at room temperature, carried out in the framework of the so-called "quantum colored noise thermostat" method. The single-particle mean kinetic energies in lithium imide as a function of temperature show a quite peculiar behavior at the moment not reproduced by ab initio lattice dynamics methods, at least as far as H and Li are concerned. As matter of fact, neither their low temperature values nor their temperature trends can be precisely explained in terms of standard phonon calculations.

Vibrational properties of CaAlH5 and α-AlH3 with different AlH6 networks studied by inelastic neutron scattering.

We performed a combined study using inelastic neutron scattering (INS) and first-principles calculations of the vibrational properties of CaAlH(5) and α-AlH(3) with different AlH(6) networks, a zigzag one-dimensional AlH(6) network for CaAlH(5), and a three-dimensional AlH(6) network for α-AlH(3). Both materials showed qualitatively similar INS spectra, in which CaAlH(5)/α-AlH(3) was mainly divided into three regions: (i) the translational modes (318/316 cm(-1)), (ii) the librational modes of the octahedral AlH(6) units (external molecular motion) in the lower frequency range and H-Al-H bond-bending modes (intra molecular motion) at a higher frequency (420-1157/513-1038 cm(-1)), and (iii) the Al-H bond-stretching modes (1238-1750/1486-1942 cm(-1)). In region ii, the appearance of both librational and bond-bending modes was determined by the networked nature of the octahedral AlH(6) units. In addition, the librational modes of AlH(6) on α-AlH(3) exhibit higher frequencies than CaAlH(5) due to the tighter bonding between the octahedral AlH(6) units. With regard to average frequencies for the Al-H, ω(S), bond-stretching modes, and average Al-H bond distances on the aluminum-based hydrides including CaAlH(5) and α-AlH(3), ω(S) showed lower frequencies that correlate with lengthening of the Al-H bond distances.

A combined experimental inelastic neutron scattering, Raman and ab initio lattice dynamics study of α-lithium amidoborane.

A combination of inelastic neutron scattering (INS) spectroscopy and Raman spectroscopy with periodic density functional theory calculations is used to provide a complete assignment of the vibrational spectra of α-lithium amidoborane (α-LiNH(2)BH(3)). The Born charge density and the atomic motion up to the decomposition temperature have been modelled. These models not only explain the nature of bonding in α-LiNH(2)BH(3) but also provide an insight into the atomic mechanisms of its decomposition. The (INS) measurements were performed in the range of 0-4000 cm(-1) on the high-resolution time-of-flight TOSCA INS spectrometer at the ISIS Spallation Neutron Source at the Rutherford Appleton Laboratory.

Raman and inelastic neutron scattering study on a melt-infiltrated composite of NaAlH4 and nanoporous carbon.

Inelastic neutron scattering and Raman scattering spectra of a melt-infiltrated composite of NaAlH(4) and active carbon fibers have been measured at low temperature for two sample conditions: as prepared and subjected to hydrogen desorption-absorption cycling. After a careful data analysis, the present experimental results have been compared to the corresponding spectroscopic data taken from bulk NaAlH(4) and Na(3)AlH(6). Evident signatures induced by infiltration process onto the NaAlH(4) phonon bands have been detected, showing up as a strong peak broadening and smoothing together with, in some cases, an energy shift. Traces of Na(3)AlH(6), appearing as an extra intensity between 130 and 200 meV, seem also confirmed. A substantial agreement between neutron and Raman measurements has been found for the as-prepared melt-infiltrated sample, while for the cycled sample the two techniques produced rather dissimilar results. However, this apparent discrepancy can be explained by considering the different penetration depths of the two spectroscopic probes. Further work, both experimental and based on ab initio simulations, is surely needed in order to rationalize the finding of the present measurements.

Neutron Compton scattering investigation of sodium hydride: from bulk material to encapsulated nanoparticulates in amorphous silica gel.

In this study we utilize neutron Compton scattering (NCS) to determine differences in nuclear momentum distributions in NaH, both as bulk material and encapsulated as nanoscale particles (from 20 to 50 nm in diameter) within an amorphous silica-gel matrix (SiGNaH). In addition, elemental Na dispersed in such a matrix is also studied (SiGNa). Data treatment and fitting of experimental spectra yields comparison of the nuclear Compton profiles and radial momentum distributions for the proton in both bulk NaH and nanoscale SiGNaH, with resultant proton kinetic energies being in agreement with previous inelastic neutron studies of bulk NaH. Slight differences in proton radial momentum distributions for bulk and nanoscale systems are witnessed and discussed. The technique of stoichiometric-fixing is applied to the backscattering spectra of each system in order to examine changes in the Na profile width, and NCS is shown to be sensitive to the chemical environment change of this heavier nucleus. Examination of the Si and O profile widths in the gel samples also supports this method.

Mobility and dynamics in the complex hydrides LiAlH4 and LiBH4.

The dynamics and bonding of the complex hydrides LiBH4 and LiAlH4 have been investigated by vibrational spectroscopy. The combination of infrared, Raman, and inelastic neutron scattering (INS) spectroscopies on hydrided and deuterided samples reveals a complete picture of the dynamics of the BH4- and AlH4 anions respectively as well as the lattice. The straightforward interpretation of isotope effects facilitates tracer diffusion experiments revealing the diffusion coefficients of hydrogen containing species in LiBH4, and LiAlH4. LiBH4 exchanges atomic hydrogen starting at 200 degrees C. Despite having an iso-electronic structure, the mobility of hydrogen in LiAlH4 is different from that of LiBH4. Upon ball-milling of LiAlH4 and LiAlD4, hydrogen is exchanged with deuterium even at room temperature. However, the exchange reaction competes with the decomposition of the compound. The diffusion coefficients of the alanate and borohydride have been found to be D approximately equal 7 x 10(-14) m2 s(-1) at 473 K and D approximately equal 5 x 10(-16) m2 s(-1) at 348 K, respectively. The BH4 ion is easily exchanged by other ions such as I- or by NH2-. This opens the possibility of tailoring physical properties such as the temperature of the phase transition linked to the Li-ion conductivity in LiBH4 as measured by nuclear magnetic resonance and Raman spectroscopy. Temperature dependent Raman measurements on diffusion gradient samples Li(BH4)1-cIc demonstrate that increasing temperature has a similar impact to increasing the iodide concentration c: the system is driven towards the high-temperature phase of LiBH4. The influence of anion exchange on the hydrogen sorption properties is limited, though. For example, Li4(BH4)(NH2)3 does not exchange hydrogen easily even in the melt.

High resolution Raman and neutron investigation of Mg(BH4)2 in an extensive temperature range.

Raman spectra of Mg(BH(4))(2) have been measured in an extensive temperature range, from 15 to 473 K. Taking into account the high temperature conversion from the alpha to the beta phase, we have observed evident signatures of this phase transition and determined the Raman vibrational spectrum of each phase. The neutron scattering spectra of the beta phase sample were also recorded. The present experimental results have been compared to the density functional theory calculations available in the literature, and a substantial agreement has been found.