Journal of Physical Chemistry C

Singlet-to-Triplet Spin Transitions Facilitate Selective 1-Butene Formation during Ethylene Dimerization in Ni(II)-MFU-4l

J. L. Mancuso, C. A. Gaggioli, L. Gagliardi, C. H. Hendon
Journal of Physical Chemistry C
Volume 125, Issue 40, Page 22036-22043
2021
A1
Published while none of the authors were employed at the CMM

Abstract 

Catalytic ethylene dimerization is an important chemical reaction that suffers from a lack of selectivity for the desired product, 1-butene. Metal-organic frameworks (MOFs) bearing Ni-based catalytic sites have been shown to yield record selectivity for 1-butene. Early efforts to understand this selectivity revealed that chain propagation and 2-butene formation are competitive with 1-butene, seemingly at odds with experimental evidence that these products are disfavored. Here, we present an alternative mechanism for selective 1-butene formation in the highest performing MOF, Ni(II)-MFU-4l. Our study reveals competing electronic spin configuration pathways that intersect along the reaction coordinate. Intersystem crossing provides an explanation for the selective formation of 1-butene in the MOF. Furthermore, we explore intersystem crossing as a unique design principle for MOF catalyst design and highlight a departure from conventional molecular catalyst design paradigms.

Elucidating the Vibrational Fingerprint of the Flexible Metal-Organic Framework MIL-53(Al) Using a Combined Experimental/Computational Approach

A.E.J. Hoffman, L. Vanduyfhuys, I. Nevjestic, J. Wieme, S.M.J. Rogge, H. Depauw, P. Van der Voort, H. Vrielinck, V. Van Speybroeck
Journal of Physical Chemistry C
122, 5, 2734-2746
2018
A1

Abstract 

In this work mid-infrared (mid-IR), far-IR, and Raman spectra are presented for the distinct (meta)stable phases of the flexible metal-organic framework MIL-53(Al). Static density functional theory (DFT) simulations are performed allowing for the identification of all IR active modes, which is unprecedented in the low-frequency region. A unique vibrational fingerprint is revealed, resulting from aluminum-oxide backbone stretching modes, which can be used to clearly distinguish the IR spectra of the closed- and large-pore phases. Furthermore, molecular dynamics simulations based on a DFT description of the potential energy surface enable to determine the theoretical Raman spectrum of the closed- and large-pore phases for the first time. An excellent correspondence between theory and experiment is observed. Both the low-frequency IR and Raman spectra show major differences in vibrational modes between the closed- and large-pore phases indicating changes in lattice dynamics between the two structures. In addition, several collective modes related to the breathing mechanism in MIL-53(Al) are identified. In particular, we rationalize the importance of the trampoline-like motion of the linker for the phase transition.

Open Access version available at UGent repository
Gold Open Access

Methane Adsorption in Zr-Based MOFs: Comparison and Critical Evaluation of Force Fields

S. Vandenbrande, T. Verstraelen, J. J. Gutierrez-Sevillano, M. Waroquier, V. Van Speybroeck
Journal of Physical Chemistry C
121 (45), 25309-25322
2017
A1

Abstract 

The search for nanoporous materials that are highly performing for gas storage and separation is one of the contemporary challenges in material design. The computational tools to aid these experimental efforts are widely available and adsorption isotherms are routinely computed for huge sets of (hypothetical) frameworks. Clearly the computational results depend on the interactions between the adsorbed species and the adsorbent, which are commonly described using force fields. In this paper, an extensive comparison and in-depth investigation of several force fields from literature is reported for the case of methane adsorption in the Zr-based Metal-Organic Frameworks UiO-66, UiO-67, DUT-52, NU-1000 and MOF-808. Significant quantitative differences in the computed uptake are observed when comparing different force fields, but most qualitative features are common which suggests some predictive power of the simulations when it comes to these properties. More insight into to the host-guest interactions is obtained by benchmarking the force fields with an extensive number of ab initio computed single molecule interaction energies. This analysis at the molecular level reveals that especially ab initio derived force fields perform well in reproducing the ab initio interaction energies. Finally, the high sensitivity of uptake predictions on the underlying potential energy surface is explored.

Open Access version available at UGent repository
Gold Open Access

Precision of Electric-Field Gradient Predictions by Density Functional Theory and Implications for the Nuclear Quadrupole Moment and Its Error Bar of the 111Cd 245 keV 5/2+ Level

L.A. Errico, K. Lejaeghere, J. Runco, S.N. Mishra, M. Renteria, S. Cottenier
Journal of Physical Chemistry C
120 (40), 23111-23120
2016
A1

Abstract 

We present ab initio calculated electric-field gradient tensors at Cd sites in a set of simple. yet diverse noncubic metals. By combining these predictions with carefully selected published experimental data, the nuclear quadrupole moment of the 245 keV 5/2(+) level of Cd-111 is determined to be 0.76(2) b. Knowing this quadrupole moment is important for time-differential perturbed angular correlation spectroscopy: decades of experimentally obtained nuclear quadrupole coupling constants for solids can now be more reliably converted into electronic structure information. For nuclear physics systematics, this is a rare opportunity to have reliable quadrupole moment information for a short-lived, level, that is not accessible to regular experimental methods. Much effort is spent on the determination of a meaningful error bar, which is an-aspect that gained only recently more attention in the context of density-functional theory predictions. This required assessing the numerical uncertainty in density functional theory predictions for electric-field gradient-tensors in solids. In contrast to quantum chemistry methods, these density functional theory predictions cannot detect systematic errors. By comparing our quadrupole moment value with an independent value Obtained from quantum chemistry calculations and experiment, we show that systematic errors are small for the systems studied here. Yet; there are indications that-density functional theory underestimates by a few percent the electric-field gradient, and therefore overestimates the quadrupole Moment by the saute amount. We point out which future work needs' to be done to characterize the possible deviations inherent to density functional theory.

Exploring the flexibility of MIL-47(V)-type materials using force field molecular dynamics simulations

J. Wieme, L. Vanduyfhuys, S.M.J. Rogge, M. Waroquier, V. Van Speybroeck
Journal of Physical Chemistry C
120 (27), 14934-14947
2016
A1

Abstract 

The flexibility of three MIL-47(V)-type materials (MIL-47, COMOC-2 and COMOC-3) has been explored by constructing the pressure-versus-volume and free energy-versus-volume curves at various temperatures ranging from 100 K to 400 K. This is done with first-principles based force fields using the recently proposed QuickFF parameterization protocol. Specific terms were added for the materials at hand to describe the asymmetry of the one-dimensional vanadium-oxide chain and to account for the flexibility of the organic linkers. The force fields are used in a series of molecular dynamics simulations at fixed volumes, but varying unit cell shapes. The three materials show a distinct pressure-versus-volume behavior, which underlines the ability to tune the mechanical properties by varying the linkers towards different applications such as nanosprings, dampers and shock absorbers.

Open Access version available at UGent repository

Determination of the Nature of the Cu coordination complexes formed in the presence of NO and NH3 within SSZ-13

I. Lezcano-Gonzalez, D. Wragg, W. A. Slawinski, K. Hemelsoet, A. Van Yperen-De Deyne, M. Waroquier, V. Van Speybroeck, A. Beale
Journal of Physical Chemistry C
119 (43), 24393-24403
2015
A1

Abstract 

Ammonia-selective catalytic reduction (NH3-SCR) using Cu zeolites is a well-established strategy for the abatement of NOx gases. Recent studies have demonstrated that Cu is particularly active when exchanged into the SSZ-13 zeolite, and its location in either the 6r or 8r renders it an excellent model system for fundamental studies. In this work, we examine the interaction of NH3-SCR relevant gases (NO and NH3) with the Cu2+ centers within the SSZ-13 structure, coupling powder diffraction (PD), X-ray absorption spectroscopy (XAFS), and density functional theory (DFT). This combined approach revealed that, upon calcination, cooling and gas exposure Cu ions tend to locate in the 8r window. After NO introduction, Cu ions are seen to coordinate to two framework oxygens and one NO molecule, resulting in a bent Cu–nitrosyl complex with a Cu–N–O bond angle of ∼150°. Whilst Cu seems to be partially reduced/changed in coordination state, NO is partially oxidized. On exposure to NH3 while the PD data suggest the Cu2+ ion occupies a similar position, simulation and XAFS pointed toward the formation of a Jahn–Teller distorted hexaamine complex [Cu(NH3)6]2+ in the center of the cha cage. These results have important implications in terms of uptake and storage of these reactive gases and potentially for the mechanisms involved in the NH3-SCR process.

Mechanical properties from periodic plane wave QM codes: the challenge of the flexible nanoporous MIL-47 (V) framework

D.E.P. Vanpoucke, K. Lejaeghere, V. Van Speybroeck, M. Waroquier, A. Ghysels
Journal of Physical Chemistry C
119, 23752-23766
2015
A1

Abstract 

Modeling the flexibility of metal–organic frameworks (MOFs) requires the computation of mechanical properties from first principles, e.g., for screening of materials in a database, for gaining insight into structural transformations, and for force field development. However, this paper shows that computations with periodic density functional theory are challenged by the flexibility of these materials: guidelines from experience with standard solid-state calculations cannot be simply transferred to flexible porous frameworks. Our test case, the MIL-47(V) material, has a large-pore and a narrow-pore shape. The effect of Pulay stress (cf. Pulay forces) leads to drastic errors for a simple structure optimization of the flexible MIL-47(V) material. Pulay stress is an artificial stress that tends to lower the volume and is caused by the finite size of the plane wave basis set. We have investigated the importance of this Pulay stress, of symmetry breaking, and of k-point sampling on (a) the structure optimization and (b) mechanical properties such as elastic constants and bulk modulus, of both the large-pore and narrow-pore structure of MIL-47(V). We found that, in the structure optimization, Pulay effects should be avoided by using a fitting procedure, in which an equation of state E(V) (EOS) is fit to a series of energy versus volume points. Manual symmetry breaking could successfully lower the energy of MIL-47(V) by distorting the vanadium–oxide distances in the vanadyl chains and by rotating the benzene linkers. For the mechanical properties, the curvature of the EOS curve was compared with the Reuss bulk modulus, derived from the elastic tensor in the harmonic approximation. Errors induced by anharmonicity, the eggbox effect, and Pulay effects propagate into the Reuss modulus. The strong coupling of the unit cell axes when the unit cell deforms expresses itself in numerical instability of the Reuss modulus. For a flexible material, it is therefore advisible to resort to the EOS fit procedure.

Open Access version available at UGent repository

Shape-selective diffusion of olefins in 8-ring solid acid microporous zeolites

A. Ghysels, S.L. Moors, K. Hemelsoet, K. De Wispelaere, M. Waroquier, G. Sastre, V. Van Speybroeck
Journal of Physical Chemistry C
119, 41, 23721-23734
2015
A1

Abstract 

The diffusion of olefins through 8-ring solid acid microporous zeolites is investigated using molecular dynamics simulations techniques and using a newly developed flexible force field. Within the context of the Methanol to Olefin (MTO) process and the observed product distribution, knowledge on the diffusion paths is essential to obtain molecular level control over the process conditions. Eight-ring zeotype materials are favorably used for the MTO process as they give a selective product distribution towards low carbon olefins. To investigate how composition, acidity and flexibility influence the diffusion paths of ethene and propene, a series of isostructural aluminosilicates (zeolites) and silicoaluminophosphates (AlPOs and SAPOs) are investigated with and without randomly distributed acidic sites. Distinct variations in diffusion of ethene are observed in terms of temperature, composition, acidity, and topology (AEI, CHA, AFX). In general, diffusion of ethene is an activated process for which free energy barriers for individual rings may be determined. We observe ring dependent diffusion behavior which can not solely be described in terms of the composition and topology of the rings. A new descriptor had to be introduced namely the accessible window area (AWA), inspired by implicit solvation models of proteins and small molecules. The AWA may be determined throughout the molecular dynamics trajectories and correlates well with the number of ring crossings at the molecular level and the free energy barriers for ring crossings from one cage to the other. The overall observed diffusivity is determined by molecular characteristics of individual rings for which AWA is a proper descriptor. Temperature-induced changes in framework dynamics and diffusivity may be captured by following the new descriptor throughout the simulations.

Open Access version available at UGent repository
Green Open Access

Reactivity of CO on carbon covered cobalt surfaces in Fischer-Tropsch Synthesis

L. Joos, I. Filot, S. Cottenier, E. Hensen, M. Waroquier, V. Van Speybroeck, R.A. van Santen
Journal of Physical Chemistry C
118 (10), 5317–5327
2014
A1

Abstract 

Fischer–Tropsch synthesis is an attractive process to convert alternative carbon sources, such as biomass, natural gas, or coal, to fuels and chemicals. Deactivation of the catalyst is obviously undesirable, and for a commercial plant it is of high importance to keep the catalyst active as long as possible during operating conditions. In this study, the reactivity of CO on carbon-covered cobalt surfaces has been investigated by means of density functional theory (DFT). An attempt is made to provide insight into the role of carbon deposition on the deactivation of two cobalt surfaces: the closed-packed Co(0001) surface and the corrugated Co(112̅1) surface. We also analyzed the adsorption and diffusion of carbon atoms on both surfaces and compared the mobility. Finally, the results for Co(0001) and Co(112̅1) are compared, and the influence of the surface topology is assessed.

Pages

Subscribe to RSS - Journal of Physical Chemistry C