Trends in Siting of Metals in Heterometallic Nd-Yb Metal-Organic Frameworks and Molecular Crystals

Several studies suggest that metal ordering within metal-organic frameworks (MOFs) is important for understanding how MOFs behave in relevant applications; however, these siting trends can be difficult to determine experimentally. To garner insight into the energetic driving forces that may lead to nonrandom ordering within heterometallic MOFs, we employ density functional theory (DFT) calculations on several bimetallic metal-organic crystals composed of Nd and Yb metal atoms. We also investigate the metal siting trends for a newly synthesized MOF. Our DFT-based energy of mixing results suggest that Nd will likely occupy sites with greater access to electronegative atoms and that local homometallic domains within a mixed-metal Nd-Yb system are favored. We also explore the use of less computationally extensive methods such as classical force fields and cluster expansion models to understand their feasibility for large system sizes. This study highlights the impact of metal ordering on the energetic stability of heterometallic MOFs and crystal structures.

Tuned Hydrogen Bonding in Rare-Earth Metal-Organic Frameworks for Design of Optical and Electronic Properties: An Exemplar Study of Y-2,5-Dihydroxyterephthalic Acid

Organic linkers in metal-organic framework (MOF) materials exhibit differences in hydrogen bonding (H-bonding), which can alter the geometric, electronic, and optical properties of the MOF. Density functional theory (DFT) simulations were performed on a photoluminescent Y-2,5-dihydroxyterephthalic acid (DOBDC) MOF with H-bonding concentrations between 0 and 100%; the H-bonds were located on both bidentate- and monodentate-bound DOBDC linkers. At 0% H-bond concentration in the framework, the lattice parameters contracted, the density increased, and simulated X-ray diffraction patterns shifted. Comparison with published experimental data identified that Y-DOBDC MOF structures must have a degree of H-bond concentration. The concentration of H-bonds in the system shifted the calculated band gap energy from 2.25 eV at 100% to 3.00 eV at 0%. The band gap energies also indicate a distinction of H-bonds formed on bidentate-coordinated linkers compared to those on monodentate linkers. Additionally, when the calculated optical spectra are compared with experimental data, the ligand-to-ligand charge-transfer luminescence in Y-DOBDC MOFs is expected to result from an average of 20-40% H-bonding with at least 50% of the bidentate linkers containing H-bonding. Therefore, the type of H-bonding within the DOBDC linker determines the electronic structure and the optical absorption of the MOF framework structure. Tuning of the H-bonding in rare-earth MOFs provides an opportunity to control the specific optical and adsorption properties of the MOF framework on the basis of reactions between the linker and the environment.

NOx Adsorption and Optical Detection in Rare Earth Metal-Organic Frameworks

Acid gases (e.g., NO_x and SO_x), commonly found in complex chemical and petrochemical streams, require material development for their selective adsorption and removal. Here, we report the NO_x adsorption properties in a family of rare earth (RE) metal-organic frameworks (MOFs) materials. Fundamental understanding of the structure-property relationship of NO_x adsorption in the RE-DOBDC materials platform was sought via a combined experimental and molecular modeling study. No structural change was noted following humid NO_x exposure. Density functional theory (DFT) simulations indicated that H₂O has a stronger affinity to bind with the metal center than NO₂, while NO₂ preferentially binds with the DOBDC ligands. Further modeling results indicate no change in binding energy across the RE elements investigated. Also, stabilization of the NO₂ and H₂O molecules following adsorption was noted, predicted to be due to hydrogen bonding between the framework ligands and the molecules and nanoconfinement within the MOF structure. This interaction also caused distinct changes in emission spectra, identified experimentally. Calculations indicated that this is due to the adsorption of NO₂ molecules onto the DOBDC ligand altering the electronic transitions and the resulting photoluminescent properties, a feature that has potential applications in future sensing technologies.

Structure and electronic properties of rare earth DOBDC metal-organic-frameworks

Here, we apply density functional theory (DFT) to investigate rare-earth metal organic frameworks (RE-MOFs), RE12(μ3-OH)16(C8O6H4)8(C8O6H5)4 (RE = Y, Eu, Tb, Yb), and characterize the level of theory needed to accurately predict structural and electronic properties in MOF materials with 4f-electrons. A two-step calculation approach of geometry optimization with spin-restricted DFT and large core potential (LCPs), and detailed electronic structures with spin-unrestricted DFT with a full valence potential + Hubbard U correction is investigated. Spin-restricted DFT with LCPs resulted in good agreement between experimental lattice parameters and optimized geometries, while a full valence potential is necessary for accurate representation of the electronic structure. The electronic structure of Eu-DOBDC MOF indicated a strong dependence on the treatment of highly localized 4f-electrons and spin polarization, as well as variation within a range of Hubbard corrections (U = 1-9 eV). For Hubbard corrected spin-unrestricted calculations, a U value of 1-4 eV maintains the non-metallic character of the band gap with slight deviations in f-orbital energetics. When compared with experimentally reported results, the importance of the full valence calculation and the Hubbard correction in correctly predicting the electronic structure is highlighted.

Theoretical insight into the single-atom catalytic mechanism of CeO2-supported Ag catalysts in CO oxidation

Revealing the accurate active center structure and the functional mechanism of CeO2-supported Ag catalysts during catalysis is extremely important for their accurate synthesis. In this work, a series of AgnCeO2 (n = 1, 2, 3, 4 and 10) model catalysts was constructed, and a DFT investigation of the reaction mechanism of CO oxidation, as a probe reaction on those catalysts, was carried out. It was found that the entire catalytic reaction was completed coordinately by Ag, lattice O and O vacancies, which could be considered as the active centers. Noticeably, the mobility of Ag atoms played an important role in the reaction process, leading to the observation of a single-atom catalytic mechanism, wherein a series of single Ag atomic species was formed during the reaction, which was beneficial to CO oxidation. With the completion of some elementary reactions, the single Ag formed during the migration of CO-Ag could return to the Ag cluster again. As expected, the single-AgCeO2 catalyst exhibited extremely high activity due to the absence of the binding effect of Ag-Ag. Nevertheless, the AgnCeO2 (n > 1) catalysts showed similar catalytic activity, which was slightly worse than that of single AgCeO2, indicating that the size effect of the Ag cluster was not obvious. These results provide the theoretical basis for further understanding the functional mechanism of the AgnCeO2 catalyst and are helpful for designing various catalysts with tailored functionalities.

A comprehensive assessment of the low-temperature thermal properties and thermodynamic functions of CeO2

Reported is an experimental and computational investigation of the low temperature heat capacity, thermodynamic functions, and thermal conductivity of stoichiometric, polycrystalline CeO₂. The experimentally measured heat capacity at T < 15 K provides an important correction to the historically accepted experimental values, and the low temperature thermal conductivity serves as the most comprehensive data set at T < 400 K available. Below 10 K, the heat capacity is observed to obey the Debye T³ law, with a Debye temperature of Θ_D = 455 K. The entropy, enthalpy, and Gibbs free energy functions are obtained from the experimental heat capacity and compared with predictions from Hubbard-corrected density functional perturbation theory calculations using the Perdew, Burke, and Ernzerhof parameterization revised for solids. The thermal conductivity is determined using the Maldonado continuous measurement technique, along with laser flash analysis, and analyzed according to the Klemens-Callaway model.

Density Functional Theory Calculations of Structural, Electronic, and Magnetic Properties of the 3d Metal Trifluorides MF3 (M = Ti-Ni) in the Solid State

We employ density functional theory with Hubbard U correction or hybrid functionals to study the series of magnetic 3d metal trifluorides MF₃ (M = Ti-Ni). Experimental lattice parameters are reproduced with an error margin of 0.5%-4.3%. Cooperative Jahn-Teller distortions are reproduced for MnF₃ , but also found in TiF₃ and CoF₃ at smaller levels compared to MnF₃ . Trends in electronic structure with respect to positions of the d bands are linked to the magnetic properties where M = Ti-Cr are weak magnetic Mott-Hubbard insulators, M = Fe-Ni are strong magnetic charge-transfer insulators and MnF₃ falls in between. Our work contributes to the characterization of the relatively unknown NiF₃ , since FeF₃ and CoF₃ have similar electronic and magnetic properties. However, NiF₃ does not show a Jahn-Teller distortion as present in CoF₃ . © 2019 Wiley Periodicals, Inc.

First-Principles Structural, Mechanical, and Thermodynamic Calculations of the Negative Thermal Expansion Compound Zr2(WO4)(PO4)2

The negative thermal expansion (NTE) material Zr₂(WO₄)(PO₄)₂ has been investigated for the first time within the framework of the density functional perturbation theory (DFPT). The structural, mechanical, and thermodynamic properties of this material have been predicted using the Perdew, Burke and Ernzerhof for solid (PBEsol) exchange-correlation functional, which showed superior accuracy over standard functionals in previous computational studies of the NTE material α-ZrW₂O₈. The bulk modulus calculated for Zr₂(WO₄)(PO₄)₂ using the Vinet equation of state at room temperature is K ₀ = 63.6 GPa, which is in close agreement with the experimental estimate of 61.3(8) at T = 296 K. The computed mean linear coefficient of thermal expansion is -3.1 × 10^-6 K^-1 in the temperature range ∼0-70 K, in line with the X-ray diffraction measurements. The mean Grüneisen parameter controlling the thermal expansion of Zr₂(WO₄)(PO₄)₂ is negative below 205 K, with a minimum of -2.1 at 10 K. The calculated standard molar heat capacity and entropy are C _P ⁰ = 287.6 and S ⁰ = 321.9 J·mol^-1·K^-1, respectively. The results reported in this study demonstrate the accuracy of DFPT/PBEsol for assessing or predicting the relationship between structural and thermomechanical properties of NTE materials.

Model representations of kerogen structures: An insight from density functional theory calculations and spectroscopic measurements

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.