Probleemstelling:

Crystalline, nanoporous materials are inherently present in several day-to-day applications such as energy storage, gas separation, and heterogeneous catalysis. Their unique combination of structural order and tunable pores makes these materials tractable for a variety of applications relying on high surface areas, high porosities, and chemical versatility. Covalent organic frameworks (COFs) form a recently proposed class of these materials consisting solely of organic moieties [1]. In general, these materials can be subdivided into 3D COFs, which are completely connected by strong covalent bonds, and 2D COFs, which form layered structures whose stacking order is determined by much weaker forces, such as dispersion and Coulomb interactions [2]. This latter category of COFs is capable of forming 1D channels, similar to graphene (see Figure 1), thereby forming extended electron clouds that promote charge migration and make these materials exceptionally interesting for photocatalysis and use in photoactive devices. However, to ensure their suitability for these applications, it is essential to limit the extent of shearing between adjacent layers of these 2D COFs so that the materials remain sufficiently stable. As the interlayer interactions are relatively weak, it is therefore crucial to understand how one can optimize the rigidity of the material to prevent interlayer shearing by tuning the components of these 2D COFs.

In recent years, the Center for Molecular Modeling (CMM) has developed various computational tools to study the behavior of nanoporous materials under external stimuli such as temperature and mechanical pressure, and we have proven their reliability through comparison with experiment [3]. While it is already possible to investigate the response of nanoporous materials upon hydrostatic pressures, the performance of the algorithm to predict the response of these materials upon more general stresses such as shear stresses has not yet been proven. Until now, one has to rely on determining the relevant mechanical moduli near equilibrium (see Figure 2), which does not provide all relevant information for the material’s resilience against shear stresses which are crucial to induce interlayer shearing. The development of a reliable protocol to predict the response of 2D COFs on more general stresses would therefore enable the in silico screening of stacking in 2D COFs, paving the way for their use in several industrial applications.

Doelstelling:

The aim of this thesis is twofold. In a first step, the student will become familiar with the basics of statistical physics and ensemble theory, obtaining the necessary background on the requirements for the derivation of a constant-stress ensemble. Here, the student will benefit from the expertise present at the CMM obtained during the derivation of constant-temperature and constant-pressure ensembles [4]. With this knowledge, it is the aim to thoroughly test the current implementations of the two constant-stress ensembles available in our in-house software code Yaff [5, 6], potentially extending them with new constant-stress algorithms from literature if deemed necessary. In a second step, the algorithm will be applied to study the shear stress necessary to induce layer shearing in 2D COFs and identify the weakest directions along which this shearing motion can be induced. While we will start with calculating the necessary shear stress for existing COFs, it is the aim to investigate how the effect of interlayer shearing can be counteracted by the inclusion of judiciously chosen, alternative building blocks in the material [2]. Moreover, the results obtained in this thesis will be compared with results based on calculating the mechanical moduli, which are easier to determine but only provide information on the material near equilibrium.

To reach the final goal of limiting the interlayer shearing in 2D COFs, different pathways can be chosen, such as the inclusion of extended linkers, linkers with bulky substituents, linkers with multiple aromatic moieties, or linkers containing atoms with large differences in electronegativity. By choosing building blocks with a particular shape or with large partial charges, resulting in a lock-and-key-like molecular conformation, one could fix the alignment of adjacent layers due to steric hindrance or electronic confinement [7]. These effects can be efficiently characterized through the free energy surface (FES) which dictates the behavior of the system. For instance, for COF-5, a conventional COF lacking an intrinsic mechanism to counteract interlayer shearing, the FES in the top pane of Figure 3 is characterized by a broadly distributed minimum. This indicates that adjacent layers have a large conformational freedom and can shear easily, thereby jeopardizing the channel size in the material. However, by introducing electronic confinement or steric hindrance in the material, as in the bottom panes of Figure 3, the spatial extent of the minimum in the FES is strongly reduced, fixing the relative position of adjacent layers. Following these pathways, newly in silico designed COFs will be optimized in this thesis to alter the dispersion and/or Coulomb interactions between the different layers, providing the necessary resistance to interlayer shearing. Subsequently, the effect of introducing guest atoms on the shearing will be probed, corresponding to the molecular environment of the layers during the use of these COFs for practical applications. The results obtained in this thesis can then be critically compared with preliminary trends that we extracted earlier based on the stacking profiles for these 2D COFs, as visualized in Figure 1.

The student will be actively coached to make him/her acquainted with the advanced simulations techniques early in the thesis year, and to transfer necessary programming skills needed to perform the research. A strong interest in programming is, however, a prerequisite.