No abstract

The design of new and interesting inorganic frameworks is an ongoing challenge in materials sciences. New structures containing double-four-ring (D4R) units have recently received particular attention. The present work focuses on the computational design of new three-dimensional frameworks made of D4R units exclusively. In a first step, our simulations explore the possible ways to assemble predefined D4R units in 3D space using a sophisticated cascade of simulated annealing/minimizations steps (autoassembly of secondary building units method). While the existing zeotype topologies were successfully generated, new topologies were predicted including very open frameworks containing new types of cages. In a second step, lattice energy minimizations were performed to estimate the viability of these hypothetical frameworks as silicate, aluminophosphaste, and gallophosphate candidates. When comparing the hypothetical structures to existing compounds, our results raise the challenging question of the appropriate chemical composition that should be aimed at for a given framework topology of interest.

A systematic computational study of gallophosphates was undertaken. First, lattice energy minimization calculations using a formal-charge shell model potential have been carried out on a series of hypothetical gallium phosphates derived from their metallogallophosphate, aluminophosphate, or aluminosilicate analogues through atomic substitution. The minimized structures show the typical features in terms of bond angles and distances as expected in zeolitic gallophosphates. Second, the crystal structures of several gallophosphates in their calcined forms have been predicted, using for each compound lattice energy minimization and an initial model derived from its as-synthesized templated form. All the modified structures thus have the same GaPO(4) composition. The lattice energies of all the simulated gallophosphate structures were compared to that of GaPO(4)-quartz as a reference structure. Interestingly, among all predicted calcined structures, various zeolitic topologies were found. The study of the energetics of these zeotypic structures showed a linear dependence of lattice energy upon density. Strikingly, a few simulated structures showed unrealistic structural features, such as important framework distortions, often associated with the occurrence of a hexameric unit in the original as-synthesized structures. Also, those gallophosphates with structural faults were found in the upper part of the energy/density plot. To address the validity of our force field calculations in these special cases, first principles calculations were undertaken on ULM-4, chosen as a typical representative structure. Indeed, the qualitative agreement found between our results and those obtained with the nonlocal density functional theory demonstrates the robustness of our force field. Further minimization also showed that the inclusion of polarizability is crucial for yielding results comparable with those obtained using first principles methods.