Template for Electronic Submission to ACS Journals



Computational structure characterization tools for the era of material informatics

Lev Sarkisov1,* and Jihan Kim2

1Institute for Materials and Processes, School of Engineering, The University of Edinburgh, Edinburgh, United Kingdom

2Department of Chemical and Biomolecular Engineering, KAIST, South Korea

*Corresponding author (Lev.Sarkisov@ed.ac.uk)

ABSTRACT: Current advances in synthesis of new porous materials outpace our ability to test them in real adsorption applications. This situation is particularly evident in the area of metal-organics frameworks (MOFs), where hundreds of new MOFs are reported every year and the number of possible MOFs is virtually infinite. How to make sense out of this vast number of existing and possible structures?

In this article, we will review the application of computational structure characterization tools for systematic description and classification of porous materials and their adsorption properties. Using examples from recent research in our groups and others, we will discuss how the information obtained from computational characterization can be used in screening protocols to identify the most promising materials for a specific application before any costly and time consuming experimental effort is committed. We will finally touch upon the need for the tools to systematically organize the information generated in computational studies. These tools combined with the recent impressive advances in synthesis of porous materials may fundamentally change the way we approach material discovery, starting the era of material informatics.

Supplemental Information file:

1. Parameters of the Poreblazer v3.0.2 structural characterization analysis:

2.58, 10.22, 3.314, 298: Helium atom sigma (Å), helium atom epsilon (K), nitrogen atom sigma (Å), temperature (K)

12.8, 1.122, 500: Cutoff distance (Å), accessible surface area coefficient (1.0 for hard sphere

surface, 1.122 for potential minimum surface), number of trials for surface area calculation

0.2: Cubelet size (Å)

20.0, 0.25: Largest anticipated pore diameter (Å), size of the bin for PSD (Å)

21908391: Random number seed

2. Simulation parameters for the results described in Section 4 of the article and related to the breath analysis application

All molecular species are considered as rigid structures. Interaction of a single molecule with the atoms of the structure involves van der Waals and Coulombic components. The van der Waals interactions are described using the Lennard-Jones potential. Lennard-Jones parameters in this study for both MOFs and adsorbates are obtained from the DREIDING (1) and UFF (2) force fields, as described below. All partial charges on atoms are obtained from the B3LYP (3,4) density functional theory calculations using the ChelpG method (5). For MOFs partial charges are taken from the literature (6). All density functional calculations are carried out using Gaussian 09 (7). The adsorbate-adsorbent Coulombic interactions are calculated using a variant of the Wolf pair-wise summation method, proposed by Fennell and Gezelter (8). For both Lennard-Jones and Coulombic interactions the cut-off radius is taken to be 13.0 Å in this work. In the calculation of Henry’s constants of adsorption, the grid size is 0.1 Å and 50 rotations of a rigid molecule per cubelet are tried. The variation of the logarithm of the Henry’s constant with the reciprocal temperature provides a well known route to the isosteric heat of adsorption, or more appropriately, to the differential enthalpy of adsorption:

[pic]

where Na corresponds to the amount adsorbed. In experiments, it is the excess amount adsorbed that is measured, whereas in simulations the absolute amount adsorbed is calculated. In principle, an appropriate conversion from the absolute to excess values must be applied in order to compare the differential enthalpy and other properties obtained in simulations to those in experiments. However, at very low pressures and loadings (the regime of interest here), the difference between the excess and absolute amounts adsorbed is negligible.

2.1 Interaction parameters for atoms of MOFs and adsorbates

|Name |σ, Å |ε, K |Forcefield |

|C |3.480 |41.00 |DREIDING |

|O |3.034 |39.90 |DREIDING |

|H |2.850 |7.64 |DREIDING |

|Zn |2.462 |62.38 |UFF |

|V |2.801 |8.05 |UFF |

|Cu |3.114 |2.51 |UFF |

All adsorbates are treated as rigid.

3. Molecular species and charges from B3LYP/CHELPG

| |x, Å |y, Å |z, Å |q, e |

|Benzaldehyde |  |  |  |

| |14 | | | |

|C |0 |0.57872 |0 |-0.02762 |

|C |-1.03935 |-0.37071 |0 |-0.02427 |

|C |-0.73823 |-1.73207 |0 |-0.11473 |

|C |0.60028 |-2.15469 |0 |-0.05039 |

|C |1.63871 |-1.2145 |0 |-0.06629 |

|C |1.33841 |0.14998 |0 |-0.14916 |

|C |-0.30373 |2.0199 |0 |0.43664 |

|H |0.5783 |2.68734 |0 |-0.01302 |

|O |-1.44911 |2.50103 |0 |-0.47229 |

|H |-2.06491 |-0.01693 |0 |0.09361 |

|H |-1.53779 |-2.46574 |0 |0.09997 |

|H |0.83198 |-3.21523 |0 |0.08818 |

|H |2.6721 |-1.54535 |0 |0.09349 |

|H |2.13671 |0.88785 |0 |0.10588 |

| | | | | |

|Btanone |  |  |  |

| |13 | | | |

|C |-1.39421 |-1.04324 |-0.09635 |-0.30874 |

|C |-0.56773 |0.20983 |0.10534 |0.59862 |

|C |0.82862 |0.04143 |0.68612 |-0.11535 |

|O |-0.99856 |1.33133 |-0.21123 |-0.53137 |

|C |1.87392 |-0.1999 |-0.42955 |-0.04698 |

|H |-0.81569 |-1.82039 |-0.6105 |0.08662 |

|H |-2.29253 |-0.81109 |-0.67129 |0.09238 |

|H |-1.68988 |-1.45594 |0.87775 |0.0803 |

|H |1.08028 |0.95837 |1.22912 |0.04462 |

|H |0.84421 |-0.79624 |1.39447 |0.02097 |

|H |2.87517 |-0.30077 |0.00197 |0.03548 |

|H |1.8856 |0.64099 |-1.13058 |0.03505 |

|H |1.6577 |-1.11418 |-0.9945 |0.00841 |

| | | | | |

|Propanol |  |  |  |  |

| |12 | | | |

|C |1.9113 |0.13512 |0 |-0.27843 |

|C |0.52251 |-0.5218 |0 |0.0784 |

|C |-0.60715 |0.50089 |0 |0.37272 |

|O |-1.86133 |-0.24124 |0 |-0.76465 |

|H |2.05524 |0.76612 |-0.88625 |0.0588 |

|H |2.70283 |-0.62201 |0 |0.07761 |

|H |2.05524 |0.76612 |0.88625 |0.0588 |

|H |0.40413 |-1.16616 |0.87998 |0.01507 |

|H |0.40413 |-1.16616 |-0.87998 |0.01507 |

|H |-0.53489 |1.1464 |-0.88977 |-0.0396 |

|H |-0.53489 |1.1464 |0.88977 |-0.0396 |

|H |-2.62112 |0.37394 |0 |0.44584 |

| | | | | |

|Decane |  |  |  |  |

| |32 | | | |

|C |-5.79404 |-0.28815 |0 |-0.18735 |

|C |-4.49383 |0.53193 |0 |0.13464 |

|C |-3.22725 |-0.34336 |0 |0.01704 |

|C |-1.91978 |0.46928 |0 |-0.05052 |

|C |-0.65366 |-0.4064 |0 |0.06945 |

|C |0.65366 |0.40645 |0 |0.06943 |

|C |1.91977 |-0.46924 |0 |-0.05051 |

|C |3.22726 |0.34338 |0 |0.01705 |

|C |4.49381 |-0.53195 |0 |0.13463 |

|C |5.79406 |0.28808 |0 |-0.18735 |

|H |-6.6775 |0.36081 |0 |0.04319 |

|H |-5.85262 |-0.93361 |0.88556 |0.03949 |

|H |-5.85262 |-0.93361 |-0.88556 |0.03949 |

|H |-4.47572 |1.19245 |-0.87957 |-0.02614 |

|H |-4.47572 |1.19245 |0.87957 |-0.02614 |

|H |-3.24453 |-1.00457 |0.88029 |-0.00694 |

|H |-3.24453 |-1.00457 |-0.88029 |-0.00694 |

|H |-1.90363 |1.13037 |-0.8803 |0.00094 |

|H |-1.90363 |1.13037 |0.8803 |0.00094 |

|H |-0.6697 |-1.06746 |0.88029 |-0.02058 |

|H |-0.6697 |-1.06746 |-0.88029 |-0.02058 |

|H |0.66971 |1.0675 |-0.88029 |-0.02058 |

|H |0.66971 |1.0675 |0.88029 |-0.02058 |

|H |1.90361 |-1.13033 |-0.8803 |0.00094 |

|H |1.90361 |-1.13033 |0.8803 |0.00094 |

|H |3.24456 |1.00458 |0.88029 |-0.00694 |

|H |3.24456 |1.00458 |-0.88029 |-0.00694 |

|H |4.47568 |-1.19246 |0.87957 |-0.02614 |

|H |4.47568 |-1.19246 |-0.87957 |-0.02614 |

|H |6.67749 |-0.36092 |0 |0.04319 |

|H |5.85267 |0.93353 |0.88556 |0.03949 |

|H |5.85267 |0.93353 |-0.88556 |0.03949 |

| | | | | |

|Syrene |  |  |  |  |

| |16 | | | |

|C |0 |0.56233 |0 |0.24701 |

|C |1.34258 |0.13273 |0 |-0.23453 |

|C |1.66795 |-1.22651 |0 |-0.03385 |

|C |0.65337 |-2.18951 |0 |-0.13303 |

|C |-0.68737 |-1.77963 |0 |-0.05652 |

|C |-1.01072 |-0.42219 |0 |-0.18704 |

|C |-0.28376 |2.00866 |0 |-0.18549 |

|C |-1.48574 |2.60773 |0 |-0.28253 |

|H |2.13467 |0.87665 |0 |0.11224 |

|H |2.70953 |-1.53241 |0 |0.08559 |

|H |0.90137 |-3.2461 |0 |0.09729 |

|H |-1.48051 |-2.521 |0 |0.08624 |

|H |-2.05419 |-0.12458 |0 |0.11474 |

|H |0.60467 |2.64016 |0 |0.12045 |

|H |-2.42134 |2.05639 |0 |0.12219 |

|H |-1.57213 |3.68919 |0 |0.12724 |

(1) Mayo, S. L.; Olafson, B. D.; Goddard, W. A. Journal of Physical Chemistry 1990, 94, 8897

(2) Rappe, A. K.; Casewit, C. J.; Colwell, K. S.; Goddard, W. A.; Skiff, W. M. Journal of the American Chemical Society 1992, 114, 10024.

(3) Becke, A. D. Journal of Chemical Physics 1993, 98, 5648.

(4) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Journal of Physical Chemistry 1994, 98, 11623.

(5) Breneman, C. M.; Wiberg, K. B. Journal of Computational Chemistry 1990, 11, 361.

(6) Yazaydin, A. O.; Benin, A. I.; Faheem, S. A.; Jakubczak, P.; Low, J. J.; Willis, R. R.; R. Q. Snurr, Chemistry of Materials 2009, 21, 1425-1430.

(7) Gaussian09.

(8) Fennell, C. J.; Gezelter, J. D.; Journal of Chemical Physics, 2006, 124, 234104.

................
................

In order to avoid copyright disputes, this page is only a partial summary.

Google Online Preview   Download