Inconsistencies in high-pressure H2 adsorption data and a lack of comparative experiment-theory studies have made the evaluation of both new and existing metal-organic frameworks (MOFs) challenging in the context of hydrogen storage applications. In this work, we performed grand canonical Monte Carlo (GCMC) simulations in nearly 500 experimentally refined MOF structures to examine the variance in simulation results because of the equation of state, H2 potential, and the effect of density functional theory structural optimization. We find that hydrogen capacity at 77 K and 100 bar, as well as hydrogen 100-to-5 bar deliverable capacity, is correlated more strongly with the MOF pore volume than with the MOF surface area (the latter correlation is known as the Chahine's rule). The tested methodologies provide consistent rankings of materials. In addition, four prototypical MOFs (MOF-74, CuBTC, ZIF-8, and MOF-5) with a range of surface areas, pore structures, and surface chemistries, representative of promising adsorbents for hydrogen storage, are evaluated in detail with both GCMC simulations and experimental measurements. Simulations with a three-site classical potential for H2 agree best with our experimental data except in the case of MOF-5, in which H2 adsorption is best replicated with a five-site potential. However, for the purpose of ranking materials, these two choices for H2 potential make little difference. More significantly, 100 bar loading estimates based on more accurate equations of state for the vapor-liquid equilibrium yield the best comparisons with the experiment.
TiCl3 and TiF3 additives are known to facilitate hydrogenation and dehydrogenation in a variety of hydrogen storage materials, yet the associated mechanism remains under debate. Here, experimental and computational studies are reported for the reactivity with hydrogen gas of bulk and ball-milled TiCl3 and TiF3 at the temperatures and pressures for which these additives are observed to accelerate reactions when added to hydrogen storage materials. TiCl3, in either the α or δ polymorphic forms and of varying crystallite size ranging from ∼5 to 95 nm, shows no detectable reaction with prolonged exposure to hydrogen gas at elevated pressures (∼120 bar) and temperatures (up to 200 °C). Similarly, TiF3 with varying crystallite size from ∼4 to 25 nm exhibits no detectable reaction with hydrogen gas. Post-exposure vibrational and electronic structure investigations using Fourier transform infrared spectroscopy and x-ray absorption spectroscopy confirm this analysis. Moreover, there is no significant promotion of H2 dissociation at either interior or exterior surfaces, as demonstrated by H2/D2 exchange studies on pure TiF3. The computed energy landscape confirms that dissociative adsorption of H2 on TiF3 surfaces is thermodynamically inhibited. However, Ti-based additives could potentially promote H2 dissociation at interfaces where structural and compositional varieties are expected, or else by way of subsequent chemical transformations. At interfaces, metallic states could be formed intrinsically or extrinsically, possibly enabling hydrogen-coupled electronic transfer by donating electrons.
Because of their extraordinary surface areas and tailorable porosity, metal-organic frameworks (MOFs) have the potential to be excellent sensors of gas-phase analytes. MOFs with open metal sites are particularly attractive for detecting Lewis basic atmospheric analytes, such as water. Here, we demonstrate that thin films of the MOF HKUST-1 can be used to quantitatively determine the relative humidity (RH) of air using a colorimetric approach. HKUST-1 thin films are spin-coated onto rigid or flexible substrates and are shown to quantitatively determine the RH within the range of 0.1-5% RH by either visual observation or a straightforward optical reflectivity measurement. At high humidity (>10% RH), a polymer/MOF bilayer is used to slow the transport of H2O to the MOF film, enabling quantitative determination of RH using time as the distinguishing metric. Finally, the sensor is combined with an inexpensive light-emitting diode light source and Si photodiode detector to demonstrate a quantitative humidity detector for low humidity environments.
The ordered monoclinic phase of the alkali-metal decahydro-closo-decaborate salt Rb2B10H10 was found to be stable from about 250 K all the way up to an order-disorder phase transition temperature of ≈762 K. The broad temperature range for this phase allowed for a detailed quasielastic neutron scattering (QENS) and nuclear magnetic resonance (NMR) study of the protypical B10H10 2- anion reorientational dynamics. The QENS and NMR combined results are consistent with an anion reorientational mechanism comprised of two types of rotational jumps expected from the anion geometry and lattice structure, namely, more rapid 90° jumps around the anion C4 symmetry axis (e.g., with correlation frequencies of ≈2.6 × 1010 s-1 at 530 K) combined with order of magnitude slower orthogonal 180° reorientational flips (e.g., ≈3.1 × 109 s-1 at 530 K) resulting in an exchange of the apical H (and apical B) positions. Each latter flip requires a concomitant 45° twist around the C4 symmetry axis to preserve the ordered Rb2B10H10 monoclinic structural symmetry. This result is consistent with previous NMR data for ordered monoclinic Na2B10H10, which also pointed to two types of anion reorientational motions. The QENS-derived reorientational activation energies are 197(2) and 288(3) meV for the C4 fourfold jumps and apical exchanges, respectively, between 400 and 680 K. Below this temperature range, NMR (and QENS) both indicate a shift to significantly larger reorientational barriers, for example, 485(8) meV for the apical exchanges. Finally, subambient diffraction measurements identify a subtle change in the Rb2B10H10 structure from monoclinic to triclinic symmetry as the temperature is decreased from around 250 to 210 K.
Solid-state hydrogen storage materials undergo complex phase transformations whose behavior are collectively determined by thermodynamic (e.g., Gibbs free energy), mechanical (e.g., lattice and elastic constants), and mass transport (e.g., diffusivity) properties. These properties depend on the reaction conditions and evolve continuously during (de)hydrogenation. Thus, they are difficult to measure in experiments. Because of this, past progress to improve solid-state hydrogen storage materials has been prolonged. Using PdHx as a representative example for interstitial metal hydride, we have recently applied molecular dynamics simulations to quantify hydrogen diffusion in the entire reaction space of temperature and composition. Here, we have further applied molecular dynamics simulations to obtain well-converged expressions for lattice constants, Gibbs free energies, and elastic constants of PdHx at various stages of the reaction. Our studies confirm significant dependence of elastic constants on temperature and composition. Specifically, a new dynamic effect of hydrogen diffusion on elastic constants is discovered and discussed.
Confining NaAlH4 in nanoporous carbon scaffolds is known to alter the sorption kinetics and/or pathways of the characteristic bulk hydride reactions through interaction with the framework at the interface, increased specific surface area of the resulting nanoparticles, decreased hydrogen diffusion distances, and prevention of phase segregation. Although the nanosize effects have been well studied, the influence of the carbon scaffold surface chemistry remains unclear. Here we compare the hydrogen sorption characteristics of NaAlH4 confined by melt infiltration in nitrogen-doped/undoped ordered nanoporous carbon of two different geometries. 23Na and 27Al MAS NMR, N2 sorption, and PXRD verify NaAlH4 was successfully confined and remains intact in the carbon nanopores after infiltration. Both the N-doped/undoped nanoconfined systems demonstrate improved reversibility in relation to the bulk hydride during hydrogen desorption/absorption cycling. Isothermal kinetic measurements indicate a lowering of the activation energy for H2 desorption by as much as 70 kJ/mol in N-doped frameworks, far larger than the reduction in carbon-only frameworks. Most interestingly, this dramatic lowering of the activation energy is accompanied by an unexpected and anomalously low NaAlH4 desorption rate in the N-doped frameworks. This suggests that the framework surface chemistry plays an important role in the desorption process and that the rate limiting step for desorption may be associated with interactions of the hydride and host surface. Our results indicate that functionalization of carbon scaffold surface chemistry with heteroatoms provides a powerful method of altering the characteristic hydrogen sorption properties of confined metal hydride systems. Furthermore, this technique may prove beneficial in the path to a viable metal hydride-based hydrogen storage system.
Solid-state hydrogen storage materials undergo complex phase transformations whose kinetics is often limited by hydrogen diffusion. Among metal hydrides, palladium hydride undergoes a diffusional phase transformation upon hydrogen uptake, during which the hydrogen diffusivity varies with hydrogen composition and temperature. Here we perform robust statistically-averaged molecular dynamics simulations to obtain a well-converged analytical expression for hydrogen diffusivity in bulk palladium that is valid throughout all stages of the reaction. Our studies confirm significant dependence of the diffusivity on composition and temperature that elucidate key trends in the available experimental measurements. Whereas at low hydrogen compositions, a single process dominates, at high hydrogen compositions, diffusion is found to exhibit behavior consistent with multiple hopping barriers. Further analysis, supported by nudged elastic band computations, suggests that the multi-barrier diffusion can be interpreted as two distinct mechanisms corresponding to hydrogen-rich and hydrogen-poor local environments.
Lightweight complex metal hydrides are of interest for use as energy-dense on-board vehicular hydrogen stores. One material of particular interest, magnesium borohydride (Mg(BH4)2), has very high hydrogen capacity, at 14.9 wt.% H, but suffers from slow kinetics and the need for extreme conditions for both dehydrogenation and rehydrogenation from magnesium diboride (MgB2). In order to establish methods to improve the kinetic properties of this system, a greater understanding of the nucleation and growth of various solid phases is essential.
Here, diammonium dodecahydro-closo-dodecaborate (NH4)2B12H12 is the ionic compound combining NH4+ cations and [B12H12]2– anions, both of which can exhibit high reorientational mobility. To study the dynamical properties of this material, we measured the proton NMR spectra and spin–lattice relaxation rates in (NH4)2B12H12 over the temperature range of 6–475 K. Two reorientational processes occurring at different frequency scales have been revealed. In the temperature range of 200–475 K, the proton spin–lattice relaxation data are governed by thermally activated reorientations of the icosahedral [B12H12]2– anions. This motional process is characterized by the activation energy of 486(8) meV, and the corresponding reorientational jump rate reaches ~108 s–1 near 410 K. Below 100 K, the relaxation data are governed by the extremely fast process of NH4+ reorientations which are not “frozen out” at the NMR frequency scale down to 6 K. The experimental results in this range are described in terms of a gradual transition from the regime of low-temperature quantum dynamics (rotational tunneling of NH4 groups) to the regime of classical jump reorientations of NH4 groups with an activation energy of 26.5 meV. Our study offers physical insights into the rich dynamical behavior of (NH4)2B12H12 on an atomic level, providing a link between the microscopic and thermodynamic properties of this compound.
Recent theoretical predictions indicate that functional groups and additives could have a favorable impact on the hydrogen adsorption characteristics of sorbents; however, no definite evidence has been obtained to date and little is known about the impact of such modifications on the thermodynamics of hydrogen uptake and overall capacity. In this work, we investigate the effect of two types of additives on the cryoadsorption of hydrogen to mesoporous silica. First, Lewis and Brønsted acid sites were evaluated by grafting aluminum to the surface of mesoporous silica (MCF-17) and characterizing the resulting silicate materials' surface area and the concentration of Brønsted and Lewis acid sites created. Heat of adsorption measurements found little influence of surface acidity on the enthalpy of hydrogen cryoadsorption. Secondly, platinum nanoparticles of 1.5 nm and 7.1 nm in diameter were loaded into MCF-17, and characterized by TEM. Hydrogen absorption measurements revealed that the addition of small amounts of metallic platinum nanoparticles increases by up to two-fold the amount of hydrogen adsorbed at liquid nitrogen temperature. Moreover, we found a direct correlation between the size of platinum particles and the amount of hydrogen stored, in favor of smaller particles.
Jensen, Steffen R.H.; Paskevicius, Mark; Hansen, Bjarne R.S.; Jakobsen, Anders S.; Moller, Kasper T.; White, James L.; Allendorf, Mark; Stavila, Vitalie; Skibsted, Jorgen; Jensen, Torben R.
The hydrogen absorption properties of metal closo-borate/metal hydride composites, M2B10H10-8MH and M2B12H12-10MH, M = Li or Na, are studied under high hydrogen pressures to understand the formation mechanism of metal borohydrides. The hydrogen storage properties of the composites have been investigated by in situ synchrotron radiation powder X-ray diffraction at p(H2) = 400 bar and by ex situ hydrogen absorption measurements at p(H2) = 526 to 998 bar. The in situ experiments reveal the formation of crystalline intermediates before metal borohydrides (MBH4) are formed. On the contrary, the M2B12H12-10MH (M = Li and Na) systems show no formation of the metal borohydride at T = 400 °C and p(H2) = 537 to 970 bar. 11B MAS NMR of the M2B10H10-8MH composites reveal that the molar ratio of LiBH4 or NaBH4 and the remaining B species is 1:0.63 and 1:0.21, respectively. Solution and solid-state 11B NMR spectra reveal new intermediates with a B:H ratio close to 1:1. Our results indicate that the M2B10H10 (M = Li, Na) salts display a higher reactivity towards hydrogen in the presence of metal hydrides compared to the corresponding [B12H12]2- composites, which represents an important step towards understanding the factors that determine the stability and reversibility of high hydrogen capacity metal borohydrides for hydrogen storage.
The salt compounds NaB11H14, Na-7-CB10H13, Li-7-CB10H13, Na-7,8-C2B9H12, and Na-7,9-C2B9H12 all contain geometrically similar, monocharged, nido-undeca(carba)borate anions (i.e., truncated icosohedral-shaped clusters constructed of only 11 instead of 12 {B-H} + {C-H} vertices and an additional number of compensating bridging and/or terminal H atoms). We used first-principles calculations, X-ray powder diffraction, differential scanning calorimetry, neutron vibrational spectroscopy, neutron elastic-scattering fixed-window scans, quasielastic neutron scattering, and electrochemical impedance measurements to investigate their structures, bonding potentials, phase-transition behaviors, anion orientational mobilities, and ionic conductivities compared to those of their closo-poly(carba)borate cousins. All exhibited order-disorder phase transitions somewhere between room temperature and 375 K. All disordered phases appear to possess highly reorientationally mobile anions (> ~1010 jumps s-1 above 300 K) and cation-vacancy-rich, close-packed or body-center-cubic-packed structures [like previously investigated closo-poly(carba)borates]. Moreover, all disordered phases display superionic conductivities but with generally somewhat lower values compared to those for the related sodium and lithium salts with similar monocharged 1-CB9H10- and CB11H12- closo-carbaborate anions. This study significantly expands the known toolkit of solid-state, poly(carba)borate-based salts capable of superionic conductivities and provides valuable insights into the effect of crystal lattice, unit cell volume, number of carbon atoms incorporated into the anion, and charge polarization on ionic conductivity.
Li2B12H12, Na2B12H12, and their closo-borate relatives exhibit unusually high ionic conductivity, making them attractive as a new class of candidate electrolytes in solid-state Li- and Na-ion batteries. However, further optimization of these materials requires a deeper understanding of the fundamental mechanisms underlying ultrafast ion conduction. To this end, we use ab initio molecular dynamics simulations and density-functional calculations to explore the motivations for cation diffusion. We find that superionic behavior in Li2B12H12 and Na2B12H12 results from a combination of key structural, chemical, and dynamical factors that introduce intrinsic frustration and disorder. A statistical metric is used to show that the structures exhibit a high density of accessible interstitial sites and site types, which correlates with the flatness of the energy landscape and the observed cation mobility. Furthermore, cations are found to dock to specific anion sites, leading to a competition between the geometric symmetry of the anion and the symmetry of the lattice itself, which can facilitate cation hopping. Finally, facile anion reorientations and other low-frequency thermal vibrations lead to fluctuations in the local potential that enhance cation mobility by creating a local driving force for hopping. We discuss the relevance of each factor for developing new ionic conductivity descriptors that can be used for discovery and optimization of closo-borate solid electrolytes, as well as superionic conductors more generally.
The atomic force microscope (AFM) offers a rich observation window on the nanoscale, yet many dynamic phenomena are too fast and too weak for direct AFM detection. Integrated cavity-optomechanics is revolutionizing micromechanical sensing; however, it has not yet impacted AFM. Here, we make a groundbreaking advance by fabricating picogram-scale probes integrated with photonic resonators to realize functional AFM detection that achieve high temporal resolution (<10 ns) and picometer vertical displacement uncertainty simultaneously. The ability to capture fast events with high precision is leveraged to measure the thermal conductivity (η), for the first time, concurrently with chemical composition at the nanoscale in photothermal induced resonance experiments. The intrinsic η of metal-organic-framework individual microcrystals, not measurable by macroscale techniques, is obtained with a small measurement uncertainty (8%). The improved sensitivity (50×) increases the measurement throughput 2500-fold and enables chemical composition measurement of molecular monolayer-thin samples. Our paradigm-shifting photonic readout for small probes breaks the common trade-off between AFM measurement precision and ability to capture transient events, thus transforming the ability to observe nanoscale dynamics in materials.
Metal organic frameworks (MOFs) are extended, nanoporous crystalline compounds consisting of metal ions interconnected by organic ligands. Their synthetic versatility suggest a disruptive class of opto - electronic materials with a high degree of electrical tunability and without the property - degrading disorder of organic conductors. In this project we determined the factors controlling charge and energy transport in MOFs and evaluated their potential for thermoelectric energy conversion. Two strategies for a chieving electronic conductivity in MOFs were explored: 1) using redox active 'guest' molecules introduced into the pores to dope the framework via charge - transfer coupling (Guest@MOF), 2) metal organic graphene analogs (MOGs) with dispersive band structur es arising from strong electronic overlap between the MOG metal ions and its coordinating linker groups. Inkjet deposition methods were developed to facilitate integration of the guest@MOF and MOG materials into practical devices.
The Sandia HyMARC team continued its development of new synthetic, modeling, and diagnostic tools that are providing new insights into all major classes of storage materials, ranging from relatively simple systems such as PdHx and MgH2, to exceptionally complex ones, such as the metal borohydrides, as well as materials thought to be very well-understood, such as Ti-doped NaAlH4. This unprecedented suite of capabilities, capable of probing all relevant length scales within storage materials, is already having a significant impact, as they are now being used by both Seedling projects and collaborators at other laboratories within HyMARC. We expect this impact to grow as new Seedling projects begin and through collaborations with other scientists outside HyMARC. In the coming year, Sandia efforts will focus on the highest impact problems, in coordination with the other HyMARC National Laboratory partners, to provide the foundational science necessary to accelerate the discovery of new hydrogen storage materials.
Metal-organic frameworks (MOFs) are highly ordered, functionally tunable supramolecular materials with the potential to improve dye-sensitized solar cells (DSSCs). Several recent reports have indicated that photocurrent can be generated in Grätzel-type DSSC devices when MOFs are used as the sensitizer. However, the specific role(s) of the incorporated MOFs and the potential influence of residual MOF precursor species on device performance are unclear. Herein, we describe the assembly and characterization of a simplified DSSC platform in which isolated MOF crystals are used as the sensitizer in a planar device architecture. We selected a pillared porphyrin framework (PPF) as the MOF sensitizer, taking particular care to avoid contamination from light-absorbing MOF precursors. Photovoltaic and electrochemical characterization under simulated 1-sun and wavelength-selective illumination revealed photocurrent generation that is clearly ascribable to the PPF MOF. Continued refinement of highly versatile MOF structure and chemistry holds promise for dramatic improvements in emerging photovoltaic technologies. (Figure Presented).
Wood, Brandon C.; Stavila, Vitalie; Poonyayant, Natchapol; Heo, Tae W.; Ray, Keith G.; Klebanoff, Leonard E.; Udovic, Terrence J.; Lee, Jonathan R.I.; Angboonpong, Natee; Pakawatpanurut, Pasit
Internal interfaces in the Li3N/[LiNH2 + 2LiH] solid-state hydrogen storage system alter the hydrogenation and dehydrogenation reaction pathways upon nanosizing, suppressing undesirable intermediate phases to dramatically improve kinetics and reversibility. Finally, the key role of solid interfaces in determining thermodynamics and kinetics suggests a new paradigm for optimizing complex hydrides for solid-state hydrogen storage by engineering internal microstructure.
Mg(BH4)2 is a promising solid-state hydrogen storage material, releasing 14.9 wt% hydrogen upon conversion to MgB2. Although several dehydrogenation pathways have been proposed, the hydrogenation process is less well understood. Here, we present a joint experimental-theoretical study that elucidates the key atomistic mechanisms associated with the initial stages of hydrogen uptake within MgB2. Fourier transform infrared, X-ray absorption, and X-ray emission spectroscopies are integrated with spectroscopic simulations to show that hydrogenation can initially proceed via direct conversion of MgB2 to Mg(BH4)2 complexes. The associated energy landscape is mapped by combining ab initio calculations with barriers extracted from the experimental uptake curves, from which a kinetic model is constructed. The results from the kinetic model suggest that initial hydrogenation takes place via a multi-step process: molecular H2 dissociation, likely at Mg-terminated MgB2 surfaces, is followed by migration of atomic hydrogen to defective boron sites, where the formation of stable B-H bonds ultimately leads to the direct creation of Mg(BH4)2 complexes without persistent BxHy intermediates. Implications for understanding the chemical, structural, and electronic changes upon hydrogenation of MgB2 are discussed.
Robust time-averaged molecular dynamics has been developed to calculate finiteerature elastic constants of a single crystal. We find that when the averaging time exceeds a certain threshold, the statistical errors in the calculated elastic constants become very small. We applied this method to compare the elastic constants of Pd and PdH0.6 at representative low (10 K) and high (500 K) temperatures. The values predicted for Pd match reasonably well with ultrasonic experimental data at both temperatures. In contrast, the predicted elastic constants for PdH0.6 only match well with ultrasonic data at 10 K; whereas, at 500 K, the predicted values are significantly lower. We hypothesize that at 500 K, the facile hydrogen diffusion in PdH0.6 alters the speed of sound, resulting in significantly reduced values of predicted elastic constants as compared to the ultrasonic experimental data. Literature mechanical testing experiments seem to support this hypothesis.
Alkali metal borohydrides can reversibly store hydrogen; however, the materials display poor cyclability, oftentimes linked to the occurrence of stable closo-polyborate intermediate species. In an effort to understand the role of such intermediates on the hydrogen storage properties of metal borohydrides, several alkali metal dodecahydro-closo-dodecaborate salts were isolated in anhydrous form and characterized by diffraction and spectroscopic techniques. Mixtures of Li2B12H12, Na2B12H12, and K2B12H12 with the corresponding alkali metal hydrides were subjected to hydrogenation conditions known to favor partial or full reversibility in metal borohydrides. The stoichiometric mixtures of MH and M2B12H12 salts form the corresponding metal borohydrides MBH4 (M = Li, Na, K) in almost quantitative yield at 100 MPa H2 and 500°C. In addition, stoichiometric mixtures of Li2B12H12 and MgH2 were found to form MgB2 at 500°C and above upon desorption in vacuum. The two destabilization strategies outlined above suggest that metal polyhydro-closo-polyborate species can be converted into the corresponding metal borohydrides or borides, albeit under rather harsh conditions of hydrogen pressure and temperature. (Chemical Equation Presented).
As the world transitions from fossil fuels to clean energy sources in the coming decades, many technological challenges will require chemists and material scientists to develop new materials for applications related to energy conversion, storage, and efficiency. Because of their unprecedented adaptability, metal-organic frameworks (MOFs) will factor strongly in this portfolio. By utilizing the broad synthetic toolkit provided by the fields of organic and inorganic chemistry, MOF pores can be customized to suit a particular application. Of particular importance is the ability to tune the strength of the interaction between the MOF pores and guest molecules. By cleverly controlling these MOF-guest interactions, the chemist may impart new function into the Guest@MOF materials otherwise lacking in vacant MOF. Herein, we highlight the concept of the Guest@MOF as it relates to our efforts to develop these materials for energy-related applicatons. Our work in the areas of H2 and noble gas storage, hydrogenolysis of biomass, light-harvesting, and conductive materials will be discussed. Of relevance to light-harvesting applications, we report for the first time a postsynthetic modification strategy for increasing the loading of a light-sensitive electron-donor molecule in the pores of a functionalized MIL-101 structure. Through the demonstrated versatility of these approaches, we show that, by treating guest molecules as integral design elements for new MOF constructs, MOF science can have a significant impact on the advancement of clean energy technologies.
This paper uses molecular dynamics simulations to study surface and interface properties of PdHx that are relevant to hydrogen storage applications. In particular, surface energies, interfacial energies, surface diffusivities, and surface segregations are all determined as a function of temperature and composition. During the course of the calculations, we demonstrated robust molecular dynamics methods that can result in highly converged finite temperature properties. Challenging examples include accurate calculations of hydrogen surface diffusivities that account for all possible atomic jump mechanisms, and constructions of surface segregation composition profiles that have negligible statistical errors. Our robust calculations reveal that the Arrhenius plots of hydrogen surface diffusion is ideally linear at low compositions, and becomes nonlinear at high compositions. The fundamental cause for this behavior has been identified. This nonlinear surface diffusion behavioe is also in good agreement with available experimental data for bulk diffusion. The implication of our calculated properties on hydrogen storage application discussed.
Hydrogen diffusion impacts the performance of solid-state hydrogen storage materials and contributes to the embrittlement of structural materials under hydrogen-containing environments. In atomistic simulations, the diffusion energy barriers are usually calculated using molecular statics simulations where a nudged elastic band method is used to constrain a path connecting the two end points of an atomic jump. This approach requires prior knowledge of the "end points". For alloy and defective systems, the number of possible atomic jumps with respect to local atomic configurations is tremendous. Even when these jumps can be exhaustively studied, it is still unclear how they can be combined to give an overall diffusion behavior seen in experiments. Here we describe the use of molecular dynamics simulations to determine the overall diffusion energy barrier from the Arrhenius equation. This method does not require information about atomic jumps, and it has additional advantages, such as the ability to incorporate finite temperature effects and to determine the pre-exponential factor. As a test case for a generic method, we focus on hydrogen diffusion in bulk aluminum. We find that the challenge of this method is the statistical variation of the results. However, highly converged energy barriers can be achieved by an appropriate set of temperatures, output time intervals (for tracking hydrogen positions), and a long total simulation time. Our results help elucidate the inconsistencies of the experimental diffusion data published in the literature. The robust approach developed here may also open up future molecular dynamics simulations to rapidly study diffusion properties of complex material systems in multidimensional spaces involving composition and defects.
The search for solid electrolytes with sufficiently high ionic conductivities and stabilities is underway to enable the commercial viability of all-solid-state rechargeable batteries. LiCB9H10 and NaCB9H10 compounds exhibit the most impressive superionic conductivities yet among complex-hydride-based materials, including this class of large-polyhedral-anion-based salts. The pseudoaromatic nature of the CB9H10 anions makes them relatively stable like their B12H122-, B10H102-, and CB11H122- cousins, rendering their salts prime candidates for incorporation into next-generation, all-solid-state devices. Preliminary cyclic voltammetry measurements indicate that only cathodic and anodic currents are observed near 0 v corresponding to Li/Na deposition on the Au electrode and Li/Na stripping, respectively, without signifi cant anodic currents, at least ≤ 5 v for both LiCB9H10 (363 K) and NaCB9H10 (303 K). The similar conductivity behaviors with temperature for LiCB9H10 and NaCB9H10 compared to those for LiCB11H12 and NaCB11H12 , and their order-of-magnitude enhancements over disordered NaCB9H10, irrespective of structural symmetries, further reinforces the notion that anion monovalency better facilitates high cation translational mobility in these large- polyhedral-anion-based systems.
We demonstrate that metal-organic frameworks (MOFs) can catalyze hydrogenolysis of aryl ether bonds under mild conditions. Mg-IRMOF-74(I) and Mg-IRMOF-74(II) are stable under reducing conditions and can cleave phenyl ethers containing β-O-4, α-O-4, and 4-O-5 linkages to the corresponding hydrocarbons and phenols. Reaction occurs at 10 bar H2 and 120 °C without added base. DFT-optimized structures and charge transfer analysis suggest that the MOF orients the substrate near Mg2+ ions on the pore walls. Ti and Ni doping further increase conversions to as high as 82% with 96% selectivity for hydrogenolysis versus ring hydrogenation. Repeated cycling induces no loss of activity, making this a promising route for mild aryl-ether bond scission.