V. Van Speybroeck

Nuclear quantum effects on zeolite proton hopping kinetics explored with machine learning potentials and path integral molecular dynamics

M. Bocus, R. Goeminne, A. Lamaire, M. Cools-Ceuppens, T. Verstraelen, V. Van Speybroeck
Nature Communications
14, 1008
2023
A1

Abstract 

Proton hopping is a key reactive process within zeolite catalysis. However, the accurate determination of its kinetics poses major challenges both for theoreticians and experimentalists. Nuclear quantum effects (NQEs) are known to influence the structure and dynamics of protons, but their rigorous inclusion through the path integral molecular dynamics (PIMD) formalism was so far beyond reach for zeolite catalyzed processes due to the excessive computational cost of evaluating all forces and energies at the Density Functional Theory (DFT) level. Herein, we overcome this limitation by training first a reactive machine learning potential (MLP) that can reproduce with high fidelity the DFT potential energy surface of proton hopping around the first Al coordination sphere in the H-CHA zeolite. The MLP offers an immense computational speedup, enabling us to derive accurate reaction kinetics beyond standard transition state theory for the proton hopping reaction. Overall, more than 0.6 μs of simulation time was needed, which is far beyond reach of any standard DFT approach. NQEs are found to significantly impact the proton hopping kinetics up to ~473 K. Moreover, PIMD simulations with deuterium can be performed without any additional training to compute kinetic isotope effects over a broad range of temperatures.

Gold Open Access

Unfolding the terahertz spectrum of soft porous crystals: rigid unit modes and their impact on phase transitions

A.E.J. Hoffman, I. Senkovska, J. Wieme, A. Krylov, S. Kaskel, V. Van Speybroeck
Journal of Materials Chemistry A
10 (33), 17254-17266
2022
A1

Abstract 

Phase transitions in exible metal-organic frameworks or soft porous crystals are mediated by low-frequency phonons or rigid-unit modes. The alteration of specic building blocks may change the lattice dynamics of these frameworks, which can inuence the phase transition mechanism. In this work, the impact of building block substitution on the rigid-unit modes in exible MIL-53 analogs with a winerack topology will be investigated via ab initio lattice dynamics calculations. First, the accuracy of the theoretical simulations is veried via experimental Raman measurements, which provide unique ngerprint vibrations in the terahertz range to characterize the phase transition. Following analysis of the low-frequency vibrations shows that there exists a set of universal rigid-unit modes inducing translations and/or rotations of the building blocks. The theoretical results demonstrate that linker substitutions have a large eect on the rigid-unit mode frequencies, whereas this is less so for inorganic chain substitutions. These ndings may help to rationally tune the phonon frequencies in soft porous crystals.

Gold Open Access

How the Layer Alignment in Two-dimensional Nanoporous Covalent Organic Frameworks Impacts Its Electronic Properties

K. S. Rawat, S. Borgmans, T. Braeckevelt, C.V. Stevens, P. Van der Voort, V. Van Speybroeck
ACS Applied Nano Materials
5, 10, 14377-14387
2022
A1

Abstract 

Two-dimensional nanoporous covalent organic frame-works (2D COFs) have gathered significant interest due to their wide range of applications. Due to the lack of strong covalent interlayer interactions, their layers can be stacked in countless ways, each resulting in unique nanoscale characteristics impacting the structural, chemical, and electronic properties. To characterize and understand the layer stacking in 2D COFs and its effect on the structural and electronic properties, we carried out a detailed density functional theory investigation on four materials, CTF-1, COF-1, COF-5, and Pc-PBBA. This entailed an in-depth evaluation of the potential energy as a function of the interlayer distance and offset, the powder X-ray diffraction (PXRD) pattern, and the electronic properties. From the potential energy surfaces, the typical slipped AA-stacking configuration was confirmed as optimal for each of the 2D COFs, with a slight offset from a perfect alignment of the layers. The statically calculated PXRD patterns based on these optimized stacking configurations showed discrepancies when compared to experimental data. Instead, when properly accounting for dynamic fluctuations by calculating the average diffraction pattern over the course of a molecular dynamics simulation, a better agreement with the experiment is obtained. Different stacking configurations also profoundly affected the electronic band structure of COFs as the interlayer pi-pi interactions are significantly impacted by the layer offset. Evidently, with decreasing layer offsets, the pi-pi interactions increase due to the layer alignment, leading to a decrease in the band gap and an increase in interlayer charge mobility. Our study highlights the need for accurate modeling of the stacking configuration in 2D COFs as a small-scale deviation in the adjacent layer position can significantly affect the structural and electronic properties.

Accurately Determining the Phase Transition Temperature of CsPbI3 via Random-Phase Approximation Calculations and Phase-Transferable Machine Learning Potentials

T. Braeckevelt, R. Goeminne, S. Vandenhaute, S. Borgmans, T. Verstraelen, J.A. Steele, M. Roeffaers, J. Hofkens, S.M.J. Rogge, V. Van Speybroeck
Chemistry of Materials
34, 19, 8561–8576
2022
A1

Abstract 

While metal halide perovskites (MHPs) have shown great potential for various optoelectronic applications, their widespread adoption in commercial photovoltaic cells or photosensors is currently restricted, given that MHPs such as CsPbI3 and FAPbI3 spontaneously transition to an optically inactive nonperovskite phase at ambient conditions. Herein, we put forward an accurate first-principles procedure to obtain fundamental insight into this phase stability conundrum. To this end, we computationally predict the Helmholtz free energy, composed of the electronic ground state energy and thermal corrections, as this is the fundamental quantity describing the phase stability in polymorphic materials. By adopting the random phase approximation method as a wave function-based method that intrinsically accounts for many-body electron correlation effects as a benchmark for the ground state energy, we validate the performance of different exchange-correlation functionals and dispersion methods. The thermal corrections, accessed through the vibrational density of states, are accessed through molecular dynamics simulations, using a phase-transferable machine learning potential to accurately account for the MHPs’ anharmonicity and mitigate size effects. The here proposed procedure is critically validated on CsPbI3, which is a challenging material as its phase stability changes slowly with varying temperature. We demonstrate that our procedure is essential to reproduce the experimental transition temperature, as choosing an inadequate functional can easily miss the transition temperature by more than 100 K. These results demonstrate that the here validated methodology is ideally suited to understand how factors such as strain engineering, surface functionalization, or compositional engineering could help to phase-stabilize MHPs for targeted applications.

Open Access version available at UGent repository
Gold Open Access

How Reproducible are Surface Areas Calculated from the BET Equation?

J.W.M. Osterrieth, J. Rampersad, D. Madden, N. Rampal, L. Skoric, B. Connolly, M.D. Allendorf, V. Stavila, J.L Snider, R. Ameloot, J. Marreiros, C. Ania, D. Azevedo, E. Vilarrasa-Garcia, B.F. Santos, X.-H. Bu, Z. Chang, H. Bunzen, N.R. Champness, S.L. Griffin, B. Cheng, R.-B. Lin, B. Coasne, S. Cohen, J.C. Moreton, Y.J. Colón, L. Chen, R. Clowes, F.-X. Coudert, Y. Cui, B. Hou, D.M. D'Alessandro, P.W. Doheny, M. Dincă, C. Sun, C. Doonan, M.T. Huxley, J.D. Evans, P. Falcaro, R. Ricco, O. Farha, K.B. Idrees, T. Islamoglu, P. Feng, H. Yang, R.S. Forgan, D. Bara, S. Furukawa, E. Sanchez, J. Gascon, S. Telalović, S.K. Ghosh, S. Mukherjee, M.R. Hill, M.M. Sadiq, P. Horcajada, P. Salcedo-Abraira, K. Kaneko, R. Kukobat, J. Kenvin, S. Keskin, S. Kitagawa, K.-i. Otake, R.P. Lively, S.J.A. DeWitt, P.L. Llewellyn, B.V. Lotsch, S.T. Emmerling, A.M. Pütz, C. Martí-Gastaldo, N.M. Padial, J. García-Martínez, N. Linares, D. Maspoch, J.A. Suárez del Pino, P.Z. Moghadam, R. Oktavian, R.E. Morris, P.S. Wheatley, J. Navarro, C. Petit, D. Danaci, M.J. Rosseinsky, A.P. Katsoulidis, M. Schroeder, X. Han, S. Yang, C. Serre, G. Mouchaham, D.S. Sholl, R. Thyagarajan, D. Siderius, R.Q. Snurr, R.B. Goncalves, S. Telfer, S.J. Lee, V.P. Ting, J.L. Rowlandson, T. Uemura, T. Iiyuka, M.A. van der Veen, D. Rega, V. Van Speybroeck, S.M.J. Rogge, A. Lamaire, K.S. Walton, L.W. Bingel, S. Wuttke, J. Andreo, O. Yaghi, B. Zhang, C.T. Yavuz, T.S. Nguyen, F. Zamora, C. Montoro, H. Zhou, A. Kirchon, D. Fairen-Jimenez
Advanced Materials
34, 27, 2201502
2022
A1

Abstract 

Porosity and surface area analysis play a prominent role in modern materials science. At the heart of this sits the Brunauer–Emmett–Teller (BET) theory, which has been a remarkably successful contribution to the field of materials science. The BET method was developed in the 1930s for open surfaces but is now the most widely used metric for the estimation of surface areas of micro- and mesoporous materials. Despite its widespread use, the calculation of BET surface areas causes a spread in reported areas, resulting in reproducibility problems in both academia and industry. To prove this, for this analysis, 18 already-measured raw adsorption isotherms were provided to sixty-one labs, who were asked to calculate the corresponding BET areas. This round-robin exercise resulted in a wide range of values. Here, the reproducibility of BET area determination from identical isotherms is demonstrated to be a largely ignored issue, raising critical concerns over the reliability of reported BET areas. To solve this major issue, a new computational approach to accurately and systematically determine the BET area of nanoporous materials is developed. The software, called “BET surface identification” (BETSI), expands on the well-known Rouquerol criteria and makes an unambiguous BET area assignment possible.

Gold Open Access

Machine Learning Potentials for Metal-Organic Frameworks using an Incremental Learning Approach

S. Vandenhaute, M. Cools-Ceuppens, S. DeKeyser, T. Verstraelen, V. Van Speybroeck
npj Computational Materials
9, 1, 19
2023
A1

Abstract 

Computational modeling of physical processes in metal-organic frameworks (MOFs) is highly challenging due to the presence of spatial heterogeneities and complex operating conditions which affect their behavior. Density functional theory (DFT) may describe interatomic interactions at the quantum mechanical level, but is computationally too expensive for systems beyond the nanometer and picosecond range. Herein, we propose an incremental learning scheme to construct accurate and data-efficient machine learning potentials for MOFs. The scheme builds on the power of equivariant neural network potentials in combination with parallelized enhanced sampling and on-the-fly training to simultaneously explore and learn the phase space in an iterative manner. With only a few hundred single-point DFT evaluations per material, accurate and transferable potentials are obtained, even for flexible frameworks with multiple structurally different phases. The incremental learning scheme is universally applicable and may pave the way to model framework materials in larger spatiotemporal windows with higher accuracy.

 

A flexible and scalable implementation of the methodology is available in Psiflow.

Gold Open Access

Stable Amorphous Solid Dispersion of Flubendazole with High Drug Loading via Solvent Electrospinning

J. Becelaere, E. Van den Broeck, E. Schoolaert, V. Vanhoorne, J. F.R. Van Guyse, M. Vergaelen, S. Borgmans, K. Creemers, V. Van Speybroeck, C. Vervaet, R. Hoogenboom, K. De Clerck
Journal of controlled release
351, November 2022, Pages 123-126
2022
A1

Abstract 

In this work, an important step is taken towards the bioavailability improvement of poorly water-soluble drugs, such as flubendazole (Flu), posing a challenge in the current development of many novel oral-administrable therapeutics. Solvent electrospinning of a solution of the drug and poly(2-ethyl-2-oxazoline) is demonstrated to be a viable strategy to produce stable nanofibrous amorphous solid dispersions (ASDs) with ultrahigh drug-loadings (up to 55 wt% Flu) and long-term stability (at least one year). Importantly, at such high drug loadings, the concentration of the polymer in the electrospinning solution has to be lowered below the concentration where it can be spun in absence of the drug as the interactions between the polymer and the drug result in increased solution viscosity. A combination of experimental analysis and molecular dynamics simulations revealed that this formulation strategy provides strong, dominant and highly stable hydrogen bonds between the polymer and the drug, which is crucial to obtain the high drug-loadings and to preserve the long-term amorphous character of the ASDs upon storage. In vitro drug release studies confirm the remarkable potential of this electrospinning formulation strategy by significantly increased drug solubility values and dissolution rates (respectively tripled and quadrupled compared to the crystalline drug), even after storing the formulation for one year.

A general synthesis of azetidines by copper-catalysed photoinduced anti-Baldwin radical cyclization of ynamides

C. Jacob, H. Baguia, A. Dubart, S. Oger, P. Thilmany, J. Beaudelot, C. Deldaele, S. Peruško, Y. Landrain, B. Michelet, S. Neale, E. Romero, C. Moucheron, V. Van Speybroeck, C. Theunissen, G. Evano
Nature Communications
13, 560
2022
A1

Abstract 

A general anti-Baldwin radical 4-exo-dig cyclization from nitrogen-substituted alkynes is reported. Upon reaction with a heteroleptic copper complex in the presence of an amine and under visible light irradiation, a range of ynamides were shown to smoothly cyclize to the corresponding azetidines, useful building blocks in natural product synthesis and medicinal chemistry, with full control of the regioselectivity of the cyclization resulting from a unique and underrated radical 4-exo-dig pathway.

Gold Open Access

Mechanistic characterization of zeolite-catalyzed aromatic electrophilic substitution at realistic operating conditions

M. Bocus, L. Vanduyfhuys, F. De Proft, B.M. Weckhuysen, V. Van Speybroeck
JACS Au (Journal of the American Chemical Society)
2, 2, 502-514
2022
A1

Abstract 

Zeolite-catalyzed benzene ethylation is an important industrial reaction, as it is the first step in the production of styrene for polymer manufacturing. Furthermore, it is a prototypical example of aromatic electrophilic substitution, a key reaction in the synthesis of many bulk and fine chemicals. Despite extensive research, the reaction mechanism and the nature of elusive intermediates at realistic operating conditions is not properly understood. More in detail, the existence of the elusive arenium ion (better known as Wheland complex) formed upon electrophilic attack on the aromatic ring is still a matter of debate. Temperature effects and the presence of protic guest molecules such as water are expected to impact the reaction mechanism and lifetime of the reaction intermediates. Herein, we used enhanced sampling ab initio molecular dynamics simulations to investigate the complete mechanism of benzene ethylation with ethene and ethanol in the H-ZSM-5 zeolite. We show that both the stepwise and concerted mechanisms are active at reaction conditions and that the Wheland intermediate spontaneously appears as a shallow minimum in the free energy surface after the electrophilic attack on the benzene ring. Addition of water enhances the protonation kinetics by about 1 order of magnitude at coverages of one water molecule per Brønsted acidic site. In the fully solvated regime, an overstabilization of the BAS as hydronium ion occurs and the rate enhancement disappears. The obtained results give critical atomistic insights in the role of water to selectively tune the kinetics of protonation reactions in zeolites.

Gold Open Access

Acidity effect on benzene methylation kinetics over substituted H-MeAlPO-5 catalysts

M. Morten, T. Cordero-Lanzac, P. Cnudde, E. A. Redekop, S. Svelle, V. Van Speybroeck, U. Olsbye
Journal of Catalysis
404, 594-606
2021
A1

Abstract 

Methylation of aromatic compounds is a key reaction step in various industrial processes such as the aromatic cycle of methanol-to-hydrocarbons chemistry. The study of isolated methylation reactions and of the influence of catalyst acidity on their kinetics is a challenging task. Herein, we have studied unidirectional metal-substituted H-MeAlPO-5 materials to evaluate the effect of acid strength on the kinetics of benzene methylation with DME. First-principle simulations showed a direct correlation between the methylation barrier and acid site strength, which depends on the metal substituent. Three H-MeAlPO-5 catalysts with high (Me = Mg), moderate (Me = Si) and low acidity (Me = Zr) were experimentally tested, confirming a linear relationship between the methylation activation energy and acid strength. The effects of temperature and reactant partial pressure were evaluated, showing significant differences in the byproduct distribution between H-MgAlPO-5 and H-SAPO-5. Comparison with propene methylation suggested that the Mg substituted catalyst is also the most active for the selective methylation of alkenes.

Open Access version available at UGent repository
Gold Open Access

Pages

Subscribe to RSS - V. Van Speybroeck