Effect of Defects on the Mechanical Deformation Mechanisms of Metal

Feb 8, 2018 - Effect of Defects on the Mechanical Deformation Mechanisms of Metal−Organic Framework‑5: A Molecular Dynamics Investigation. Bin Zhe...
4 downloads 9 Views 2MB Size
Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

Article

Effect of Defects on the Mechanical Deformation Mechanisms of Metal-Organic Framework#5: A Molecular Dynamics Investigation Bin Zheng, Fang Fu, Lian Li Wang, Jinlei Wang, Lifei Du, and Hui Ling Du J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b10928 • Publication Date (Web): 08 Feb 2018 Downloaded from http://pubs.acs.org on February 9, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Effect of Defects on the Mechanical Deformation Mechanisms of Metal−Organic Framework‑5: A Molecular Dynamics Investigation Bin Zheng,*,† Fang Fu, † Lian Li Wang,† Jinlei Wang,† Lifei Du,† Huiling Du† †

School of Materials Science and Engineering, Xi’an University of Science and

Technology, Xi’an 710054, PR China

E-mail: [email protected]

TITLE RUNNING HEAD: The role of defects in the mechanical deformation of MOFs was revealed to obtain deformation mechanism and tunable mechanical properties.

ABSTRACT: Defects are likely ubiquitous in metal-organic frameworks (MOFs). Investigation on the effect of defects on the mechanical deformation of MOFs remain at a nascent stage. In this study molecular dynamics simulation was adopted to investigate the deformation mechanisms of defected MOF-5. Our results show that the defects reduce the yield strength via inducing the nucleation of dislocations. The type of dislocations depends on the defect structure, presented as a subject increased scientific interest. Also, the defected MOF-5 was deformed under different loading modes (compression, shear and tension) to understand the deformation mechanism. Revealing the 1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

mechanical deformation of defected MOFs is expected to yield considerable benefits in both controlling and tailoring the mechanical properties of MOFs.

1. INTRODUCTION Metal−organic frameworks (MOFs) are interesting as new types of organic−inorganic hybrid materials, which can incorporate a variety of metal ions/clusters and organic linkers with different chemical or physical properties, constructing a well-designed coordination to form highly ordered spatial networks and diverse porous structures. These kinds of structure and property characteristics emphasizes that MOFs can be promising choices for a variety of applications, including gas adsorption/separation,1 energy storage/conversion,2 biological processes,3, 4 chemic sensors,5 and catalysis.6 The vast majority of studies on MOFs focuses on “ideal crystals” with faultlessly arranged identical groups of atoms in space. Nevertheless, the real MOF structures always contain a significant number of various defects.7-9 One typical example is the zeolitic imidazolate framework-8 (ZIF-8), which demonstrates strong adsorption of water in laboratory test,10 contrary to the theoretical prediction of being hydrophobic.11 This constitutes the strongest evidence, which proved that the hydrophilic defects must be present in ZIF-8, even for the most precise synthesis. The ubiquitous nature of defects in MOFs establishes a critical parameter concerning the overall performance of a material. Thus, besides the composition (metal/metal-oxide nodes and organic linkers) of MOFs, the defects are also capable of tuning and tailoring the chemical and physical properties of the materials for targeted 2

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

applications.7, 8 As an example, the point defects (such as metal/ligand vacancy and dangling linkers) in MOFs can drastically enhance adsorption of gases like H2, CH4, CO2 and CO.12-14 Additionally, the high accessibility of unsaturated metal ions and uncoordinated ligands structures in MOFs constitute a great source of active sites, which can act as reactive centers in many catalytic reactions.9, 15 Special attention has been paid to the stability of defected MOFs.8 Schmidt et al., 16

found that the linkers and metals vacancy in ZIF-8, impose only an insignificant

energy penalty, indicating the negligible thermodynamic effect on the framework stability. However, the challenge concerns the mechanical stability of defected MOFs commonly used in practical application. Much research has examined perfect MOFs structures, which exhibit low mechanical modulus and yield strength,17-19 while the shear-collapse mode dominated the plastic deformation.20 Thus a significant question arises concerning MOFs structure: how do defects in MOFs influence their mechanical properties? The available expertise and theories in deforming metal, ceramic and polymer materials, are not sufficient to be adopted in the mechanical deformation of MOFs. Therefore, it is essential to investigate the effect of defects on the mechanical deformation of MOFs from both theoretical and practical perspective. Experimental observation of the deformation of MOFs can be achieved using an in-situ transmission electron microscopy (TEM) technique;21,

22

however it is

ambitious to directly observe the defects in MOFs. Computational methods can be highly effective in solving these types of science problems, particularly when experimental tools are ineffective. Most of theoretical investigations on the 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

mechanical properties of MOFs have mainly focused on the elastic deformation stage of prefect MOFs.23-25 However, the study of the inelastic deformation of defected MOFs, to our knowledge, is scarce. In this study, molecular dynamics (MD) simulation was adopted to investigate the mechanical deformation of defected MOF-5. The metal/ligand vacancy, owing to its common existence in various MOFs, was chosen as the defect structure. The evolution of internal structure under loading was analyzed to reveal the effect of defects on the nucleation and propagation of dislocations. Different types of deformation modes were examined to clarify the loading dependence of the mechanical deformation mechanism in MOF-5. The understanding of mechanical deformation in defected MOFs is expected to contribute to the design and the usage of MOFs in applications related to the mechanical loading. 2. CALCULATION MODEL AND METHODS The basic unit of MOF-5 consists of tetrahedral Zn4O groups connected by 1,4-benzenedicarboxylate (BDC) ligands. Eight Zn4O groups are placed on the corner position of a simple cubic and the BDC ligands bridge neighbor groups to finally form the unit cell MOF-5 (shaded part in Figure 1). The orientation of all linkers in MOF-5 was along crystal direction. To perform the simulation, a supercell was built by replicating the cubic unit cell 6 × 6 × 6 times along x, y, and z directions. Thus, a dimension of 15.50 × 15.50 × 15.50 nm3 was formed, based on the experimental lattice parameters.26 The perfect structure consisted of 91584 atoms (left part in Figure 1). Periodic boundary conditions (PBC) were applied along the three main 4

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

coordinate directions. The defect was created in the center position of the supercell (shaded part in Figure 1). In the present research, four typical point defects are presented: a) MOF-5 missing 1 linker, b) missing 8 linkers, c) missing 1 Zn4O groups and d) missing 8 Zn4O groups and 12 linkers (right part in Figure 1). The metal/ligand vacancy in MOF-5 has been observed during the experiment.27 Although other types of defect structures such as cracks, grooves, and terraces were reported in the experiment,7 this study only focused on the vacancy defect for two main reasons. First, the point defect is the basic unit constituting any complex defects. Secondly, although experimental observations of the atomic structure for MOFs defects are scarce, the atomic structure of point defects is obvious. The LAMMPS package28, 29 was used to perform the molecular dynamics (MD) simulations. Prior to deformation, an energy minimization based on the conjugate gradient (CG) algorithm was performed at zero temperature to guarantee the atomic positions during geometric optimization. Then, the system was thermally equilibrated at 300 K and zero pressure for 50 ps, using isothermal−isobaric (NPT) ensemble with a time step of 0.1 fs. The used time step is short enough to satisfy the requirement of the ReaxFF model to cope with unusually large interatomic forces. The coupling constants for the Nosé−Hoover thermostat and barostat were set to 0.1 ps and 5 ps respectively. The deformation can be induced by changing the size and/or the shape of the simulation box during a dynamics run. Six adjustable parameters (x, y, z, xy, xz, yz) can be used to change a box. The parameters x, y and z denote the dimensional size, 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

while the xy, xz, and yz parameters are the tilt factors of the simulation box. To perform the compressive (or stretching) deformation, any or all the parameters x, y, and z can be continuously decreased (or increased). The other parameters xy, xz, and yz can be used to induce the shear deformation. In this study, only uniaxial deformation was considered. Then, the box parameters of z and xy was continuously changed to realize the uniaxial compressive (or stretching) and shear deformation respectively. Other parameters were adjusted to reach zero pressure using the NPT ensemble. The linear and angular momentum of this simulation system was set to zero, avoiding the translational and rotational motion of the box during deformation. In every deformation step, the pressure along the deformation direction was collected to construct the loading-strain relation. An engineering strain rate of 3 × 109 s−1 was chosen, based on the reported results.20 The unloading dimensions of the cell during deformation were fully relaxed to zero stress using a Nosé−Hoover barostat. The VMD software30 was applied to visualize the atomic structure. The reliability of any MD simulation results mainly depends on the force field describing the atomic interactions in the materials. The typical semi-empirical force field consists of bonded (bond, angle and torsion) and nonbonded (van der Waals and electrostatic) part. This type of force field has been proved to be reliable in computing the guest adsorption and diffusion in MOF-5, and even the framework lattice dynamics.31 However, due to the bonded nature, the common semi-empirical force field cannot describe the inelastic deformation of MOFs involving bond breaking and forming.32 The reactive force field ReaxFF is also a semi-empirical force field, 6

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

although it owns complicated mathematical expressions with parameters derived from quantum mechanical results.33 In particular, the bond-order dependence, which depends on the local atomic environment, was included in ReaxFF; therefore, the chemical reactions involving bond breaking and forming could be described. Aiming at MOF-5, Adri van Duin et al.34 developed ReaxFF parameters to study the stability of MOF-5 against H2O. In that study, the Zn-O interaction parameters have been confirmed by the prediction of properties of four ZnO crystal structures and their interactions with water.35 Alejandro Strachan et al.

20, 36

used the ReaxFF to

investigate the mechanical deformation and the shockwave energy dissipation of MOF-5 and obtained reliable results. Current research used ReaxFF force field and predicted lattice constant of MOF-5 (Table S1), which was insignificantly overestimated compared to experimental values. Herein, the effect of defects on the lattice constants of MOF-5 was considered negligible (Table S1). 3. RESULTS AND DISCUSSION Figure 2 illustrates the simulated compressive stress versus strain curves for the defected MOF-5. According to figure 2(a), it is noted that the yield strength of the frameworks is slightly reduced with an increasing number of missing linkers. Compared to the absence of linkers, the missing Zn4O groups can further weaken the yield strength of MOF-5 frameworks (Figure 2b). The effect of defects on Young’s modulus for the current MOF-5 samples can be considered negligible. The clear divergence of the position and the strength of the stress peaks (Figure 2b) indicates different plastic deformation behaviors. The evolution of internal atomic structure is 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

required to determine the hidden deformation mechanism. To clearly capture the structure evolution of MOF-5, the unit cell of the original framework was transformed into a simple cubic (Figure 3). The nodes and the linkages in the cubic (located right in Figure 3) denotes the Zn4O groups and the 1,4-benzenedicarboxylate (BDC) ligands, respectively. The linkage breaks if its length changes over 20% of the equilibrium value. In addition, the volume variation of the simple cubic was applied to characterize the degree of local deformation.20 Figure 4 presents the internal structure evolution of defected MOF-5 under compressive loading, employing the simple cubic model (Figure 3). In the initial structure, the strain localization (red area in Figure 4b) can be induced by the defects. The Zn4O group vacancy was found to be more likely to induce strain localization (Figure 4 and S1). External loading can result in the propagation of the strain localization (strain = 2.52% in Figure 4b). The strain localization in MOF-5 with linker vacancy requires larger driving force and scatters throughout the whole system (3.72% and 4.14% strain in Figure 4a), while, the strain localization in MOF-5 with Zn4O vacancy is concentrated in the defect position (Figure 4b). The slip along [100] was observed in the current deformation model (xz plane in Figure 4). One slip consisted of three adjacent atomic layers around the defect position and twin crystal formed as the result of the slip. By switching into the yz plane projection, no trace of a layer slip on the boundary can be observed (strain = 6.84% and 10.14% in Figure 4a and b respectively), excluding the possibility of [110] slip. A highly interesting result is that the three horizontal lines (atomic layers denoted 8

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

as a, b and c in Figure 4a) curved with increasing strain. This point was absent in the MOF-5 sample with missing 8 Zn4O groups and 12 linkers (Figure 4b). Thus, different slipping mechanisms are probable to exist in the two defected MOFs samples. To reveal the slipping mechanism, the three glide layers were tracked during compression (Figure 5). For MOF-5 missing 1 linker (Figure 5a), the atomic layers partially slip along the x-axis. This movement can be traced at the strain of 4.14%. Further loading results in distinct dislocation (denoted by lines 1 and 2 in Figure 5). At strain of 5.70%, the lateral position of the atoms below the dislocation line 2 was almost unchanged, while the glide occurred between the two dislocation lines. The dislocation lines propagated under continuous compression, reaching a strain of 7.74%. The dislocation lines were always parallel to the slipping direction, which is consistent with the character of the screw dislocation. For MOF-5 missing 8 Zn4O and 12 linkers (Figure 5b), the trace of atomic layer slip can be found at a lower strain (2.52 %) compared to the 4.14% strain for MOF-5 missing 1 linker. Further loading can introduce the glide of the whole atomic layer. The dislocation line was perpendicular to the glide direction (Figure 5b), fitting the character of the edge dislocation. The above analysis indicates that the type of dislocation in MOF-5 depends on the defects. The screw dislocation is most probable to appear in MOF-5 samples missing linkers, while the absence of Zn4O groups is likely to induce the nucleation and the propagation of the edge dislocation. The videos (video-1.avi, video-2.avi, 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

video-3.avi, video-4.avi) further support this dependence between dislocations and defects in MOF-5. The appearance of different types of dislocations in the defected MOF-5 can be understood through analysis of the internal strain of dislocations and their interactions with defects (Figure 6). For the screw dislocation, the dislocation planes slide against each other via shear mode (Figure 6a). For the edge dislocation, in the region where there is no extra plane (pink area in Figure 6b), the lattice stretches, while, in the region where there is an extra plane (green area in Figure 6b), the lattice is compressed. Compared to the compressive and tensile lattice strain, the shear lattice strain can be easily induced in MOFs materials.17 Thus, the screw dislocation corresponding to the shear lattice strain prefers to appear in MOF-5. For the defected MOF-5, when only the linkers are absent the defect hardly affects the framework structure and contributes little to the formation of dislocations. For MOF-5 missing Zn4O groups, the node vacancy exists in the lattice and can effectively weaken the strain concentration due to the edge dislocation (Figure 6c). Thus, the edge dislocation becomes easy to nucleate and propagate in MOF-5 missing Zn4O groups (Figure 5 and Figure S1). Under compression, the shear-collapse linkers have been observed in MOF-5.20 Only the linkers in MOF-5, which are parallel to the loading direction, contribute to the deformation of the whole system (Figure 4). The above result was not sufficient to define the deformation in practical systems, where multiple deformation modes are involved. The following study focuses on the effect of the loading mode on the 10

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

deformation mechanism. One deformation mode is the shear, where the loading directly drives the ligand-Zn4O-ligand angle torsion. When the strain reaches approximately ~20% (Figure 7 and S2), the shear strain appears to be homogeneous. This corresponds to significantly higher recoverable strain compared to that in compression (< 3%). Further loading can induce the shear strain localization in the defect position (Figure 7 and video-5.avi). The edge dislocation, instead of the screw dislocation, could be responsible for the motion of the atomic layer. It can be attributed to the shear loading which can directly drive the slipping of the whole atomic layer, fitting the character of the edge dislocation. In addition, only the linkers normal to the loading direction contribute to the whole deformation. Next, the tensile deformation mode was investigated. Figure 8 illustrates the stretching stress-strain relation and the size of the lateral section of MOF-5 that is missing 8 Zn4O groups and 12 linkers under tension, compared to that under compression. The yield stress and strain, were significantly improved, reaching values approximately 8 times and 4 times greater than those for compression, respectively (Figure 8a). Another interesting point is that the cross section of MOF-5 sharply shrinks under tension, especially in the plastic deformation stage (Figure 8b and video-6.avi). Additionally, the lateral size of MOF-5 under compression slightly changes. The evolution of the internal structure can be used to reveal the deformation mechanism under tension (Figure 9). By increasing loading, the framework lines 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

along tensile direction are stretched and broken, corresponding to the yield point (“B” in Figure 8a and 9). The breakage involved the bond linkage between BDC ligands and Zn4O groups and initially occurred at the defect position. Meanwhile, the angle variation between BDC ligands and Zn4O groups corresponded to the shear deformation. Certainly, bond stretching (tensile) was much harder than angle torsion (shear), as explained the large yield strength and strain in tension (Figure 8a). Further tensile loading resulted in enhanced bond breaking, leaving the linkers dangling and creating in the framework hexagonal units (denoted by stars in Figure 9). Consequently, bond breaking, and the vacancy defect gradually expanded (“F” and “G” in Figure 9), and the shear-collapse of the linkers along the lateral section could also be easily induced, because of the reduction of constraints (such as in the hexagonal unit). Thus, the linker shearing contributes to the large contraction of lateral section of MOF-5 during stretching (Figure 8b). 4. CONCLUSIONS MD simulation was employed to reveal the mechanisms responsible for the mechanical deformation of the defected MOF-5. The defects were found to induce the nucleation of dislocations and then to reduce the elastic yield strength in compression. It was demonstrated that the screw dislocation prefers to appear in the MOF-5 framework missing linkers, while the absence of Zn4O group easily induces the nucleation of the edge dislocation. In shear deformation, the metal-ligand-metal angle distortion was easily induced and then the homogenous tilt along the loading direction was obtained, as this contributed to a relative large elastic strain (close to 20%). Still, 12

ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the strain localization mainly appears in the defect position. In tensile deformation, the linkages were stretched, which does not easily occur in shear deformation, and contributes to the improved elastic yield strength of MOF-5. The defects provide the initial position for the linker- Zn4O bond breaking, which in return can induce the expansion of defects. Furthermore, bond breaking can result in the shear-collapse of bonds perpendicular to the stretching direction, as this contributes to a severe shrinkage of the lateral section in MOF-5. The results indicated the importance of defects in the mechanical deformation of MOFs and could help to understand the experimental testing of the stiffness and ductility of MOF materials, while the suggestions obtained from the simulations can enable the design of new MOFs.

ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: Videos showing the defect evolution in MOF-5. The lattice parameters, atomistic snapshots and stress-strain curve of MOF-5.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

Notes 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The authors declare no competing financial interest.

ACKNOWLEDGMENT. This work was supported by the Natural Science Foundation of China under grant 21503165 and 51372197, Shaanxi Province 100 plan, and the Key Innovation Team of Shaanxi Province (2014KCT-04).

REFERENCES

1.

Adil, K.; Belmabkhout, Y.; Pillai, R. S.; Cadiau, A.; Bhatt, P. M.; Assen, A. H.; Maurin, G.;

Eddaoudi, M., Gas/Vapour Separation Using Ultra-Microporous Metal-Organic Frameworks: Insights into the Structure/Separation Relationship. Chem. Soc. Rev. 2017, 46, 3402-3430. 2.

Cao, X.; Tan, C.; Sindoro, M.; Zhang, H., Hybrid Micro-/Nano-Structures Derived from

Metal-Organic Frameworks: Preparation and Applications in Energy Storage and Conversion. Chem. Soc. Rev. 2017, 46, 2660-2677. 3.

Lian, X.; Fang, Y.; Joseph, E.; Wang, Q.; Li, J.; Banerjee, S.; Lollar, C.; Wang, X.; Zhou,

H.-C., Enzyme-MOF (Metal-Organic Framework) Composites. Chem. Soc. Rev. 2017, 46, 3386-3401. 4.

Doonan, C.; Riccò, R.; Liang, K.; Bradshaw, D.; Falcaro, P., Metal–Organic Frameworks at

the Biointerface: Synthetic Strategies and Applications. Acc. Chem. Res. 2017, 50, 1423-1432. 5.

Stassen, I.; Burtch, N.; Talin, A.; Falcaro, P.; Allendorf, M.; Ameloot, R., An Updated

Roadmap for the Integration of Metal-Organic Frameworks with Electronic Devices and Chemical Sensors. Chem. Soc. Rev. 2017, 46, 3185-3241. 6.

Huang, Y.-B.; Liang, J.; Wang, X.-S.; Cao, R., Multifunctional Metal-Organic Framework 14

ACS Paragon Plus Environment

Page 14 of 31

Page 15 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Catalysts: Synergistic Catalysis and Tandem Reactions. Chem. Soc. Rev. 2017, 46, 126-157. 7.

Fang, Z.; Bueken, B.; De Vos, D. E.; Fischer, R. A., Defect-Engineered Metal–Organic

Frameworks. Angew. Chem. Int. Ed. 2015, 54, 7234-7254. 8.

Sholl, D. S.; Lively, R. P., Defects in Metal–Organic Frameworks: Challenge or Opportunity?

J. Phys. Chem. Lett. 2015, 6, 3437-3444. 9.

Chaemchuen, S.; Luo, Z.; Zhou, K.; Mousavi, B.; Phatanasri, S.; Jaroniec, M.; Verpoort, F.,

Defect Formation in Metal–Organic Frameworks Initiated by the Crystal Growth-Rate and Effect on Catalytic Performance. J. Catal. 2017, 354, 84-91. 10. Zhang, K.; Lively, R. P.; Zhang, C.; Koros, W. J.; Chance, R. R., Investigating the Intrinsic Ethanol/Water Separation Capability of ZIF-8: An Adsorption and Diffusion Study. J. Phys. Chem. C 2013, 117, 7214-7225. 11. Ortiz, A. U.; Freitas, A. P.; Boutin, A.; Fuchs, A. H.; Coudert, F.-X., What Makes Zeolitic Imidazolate Frameworks Hydrophobic or Hydrophilic? The Impact of Geometry and Functionalization on Water Adsorption. Phys. Chem. Chem. Phys. 2014, 16, 9940-9949. 12. Wang, Z.; Sezen, H.; Liu, J.; Yang, C.; Roggenbuck, S. E.; Peikert, K.; Fröba, M.; Mavrandonakis, A.; Supronowicz, B.; Heine, T.; Gliemann, H.; Wöll, C., Tunable Coordinative Defects in UHM-3 Surface-Mounted MOFs for Gas Adsorption and Separation: A Combined Experimental and Theoretical Study. Microporous Mesoporous Mater. 2015, 207, 53-60. 13. Wu, H.; Chua, Y. S.; Krungleviciute, V.; Tyagi, M.; Chen, P.; Yildirim, T.; Zhou, W., Unusual and Highly Tunable Missing-Linker Defects in Zirconium Metal–Organic Framework UIO-66 and Their Important Effects on Gas Adsorption. J. Am. Chem. Soc. 2013, 135, 10525-10532. 14. Barin, G.; Krungleviciute, V.; Gutov, O.; Hupp, J. T.; Yildirim, T.; Farha, O. K., Defect 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Creation by Linker Fragmentation in Metal–Organic Frameworks and Its Effects on Gas Uptake Properties. Inorg. Chem. 2014, 53, 6914-6919. 15. Lee, J.; Farha, O. K.; Roberts, J.; Scheidt, K. A.; Nguyen, S. T.; Hupp, J. T., Metal-Organic Framework Materials as Catalysts. Chem. Soc. Rev. 2009, 38, 1450-1459. 16. Zhang, C.; Han, C.; Sholl, D. S.; Schmidt, J. R., Computational Characterization of Defects in Metal–Organic Frameworks: Spontaneous and Water-Induced Point Defects in ZIF-8. J. Phys. Chem. Lett. 2016, 7, 459-464. 17. Tan, J.-C.; Civalleri, B.; Lin, C.-C.; Valenzano, L.; Galvelis, R.; Chen, P.-F.; Bennett, T. D.; Mellot-Draznieks, C.; Zicovich-Wilson, C. M.; Cheetham, A. K., Exceptionally Low Shear Modulus in a Prototypical Imidazole-Based Metal-Organic Framework. Phys. Rev. Lett. 2012, 108, 095502. 18. Hegde, V. I.; Tan, J.-C.; Waghmare, U. V.; Cheetham, A. K., Stacking Faults and Mechanical Behavior Beyond the Elastic Limit of an Imidazole-Based Metal Organic Framework: ZIF-8. J. Phys. Chem. Lett. 2013, 4, 3377-3381. 19. Zheng, B.; Zhu, Y.; Fu, F.; Wang, L. L.; Wang, J.; Du, H., Theoretical Prediction of the Mechanical Properties of Zeolitic Imidazolate Frameworks (ZIFs). RSC Adv. 2017, 7, 41499-41503. 20. Banlusan, K.; Antillon, E.; Strachan, A., Mechanisms of Plastic Deformation of Metal– Organic Framework-5. J. Phys. Chem. C 2015, 119, 25845-25852. 21. Miao, Y.-R.; Su, Z.; Suslick, K. S., Energy Storage During Compression of Metal–Organic Frameworks. J. Am. Chem. Soc. 2017, 139, 4667-4670. 22. Su, Z.; Miao, Y.-R.; Mao, S.-M.; Zhang, G.-H.; Dillon, S.; Miller, J. T.; Suslick, K. S., 16

ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Compression-Induced Deformation of Individual Metal–Organic Framework Microcrystals. J. Am. Chem. Soc. 2015, 137, 1750-1753. 23. Zheng, B.; Wang, L. L.; Hui, J. C.; Du, L.; Du, H.; Zhu, M., Impact of Mechanical Deformation on Guest Diffusion in Zeolitic Imidazolate Frameworks. Dalton Trans. 2016, 45, 4346-4351. 24. Ryder, M. R.; Tan, J.-C., Explaining the Mechanical Mechanisms of Zeolitic Metal-Organic Frameworks: Revealing Auxeticity and Anomalous Elasticity. Dalton Trans. 2016, 45, 4154-4161. 25. Bourg, L. B. d.; Ortiz, A. U.; Boutin, A.; Coudert, F.-X., Thermal and Mechanical Stability of Zeolitic Imidazolate Frameworks Polymorphs. APL Mater. 2014, 2, 124110. 26. Eddaoudi, M.; Kim, J.; Rosi, N.; Vodak, D.; Wachter, J.; O'Keeffe, M.; Yaghi, O. M., Systematic Design of Pore Size and Functionality in Isoreticular MOFs and Their Application in Methane Storage. Science 2002, 295, 469-472. 27. Ravon, U.; Savonnet, M.; Aguado, S.; Domine, M. E.; Janneau, E.; Farrusseng, D., Engineering of Coordination Polymers for Shape Selective Alkylation of Large Aromatics and the Role of Defects. Microporous Mesoporous Mater. 2010, 129, 319-329. 28. http://lammps.sandia.gov. 29. Plimpton, S., Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19. 30. Humphrey, W.; Dalke, A.; Schulten, K., Vmd: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33-38. 31. Amirjalayer, S.; Tafipolsky, M.; Schmid, R., Molecular Dynamics Simulation of Benzene Diffusion in MOF-5: Importance of Lattice Dynamics. Angew. Chem. Int. Ed. 2007, 46, 463-466. 17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

32. Ortiz, A. U.; Boutin, A.; Fuchs, A. H.; Coudert, F.-X., Investigating the Pressure-Induced Amorphization of Zeolitic Imidazolate Framework ZIF-8: Mechanical Instability Due to Shear Mode Softening. J. Phys. Chem. Lett. 2013, 4, 1861-1865. 33. van Duin, A. C. T.; Dasgupta, S.; Lorant, F.; Goddard, W. A., Reaxff:  A Reactive Force Field for Hydrocarbons. J. Phys. Chem. A 2001, 105, 9396-9409. 34. Han, S. S.; Choi, S.-H.; van Duin, A. C. T., Molecular Dynamics Simulations of Stability of Metal-Organic Frameworks against H2O Using the Reaxff Reactive Force Field. Chem. Commun. 2010, 46, 5713-5715. 35. Raymand, D.; van Duin, A. C. T.; Spångberg, D.; Goddard, W. A.; Hermansson, K., Water Adsorption on Stepped ZnO Surfaces from MD Simulation. Surf. Sci. 2010, 604, 741-752. 36. Banlusan, K.; Strachan, A., Shockwave Energy Dissipation in Metal–Organic Framework MOF-5. J. Phys. Chem. C 2016, 120, 12463-12471.

18

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Computational model of defected MOF-5. The ball and stick display style on the right is used to emphasize the defect position.

Figure 2. Compressive stress−strain relations of defected MOF-5.

Figure 3. The process of simplifying the model for visualizing clarification. Blue gray, red and gray denote Zn, O and C atoms, respectively. All H atoms are hidden for clarity.

Figure 4. Atomic snapshots illustrate plastic deformation of defected MOF-5 during compression. The frame lines were explained in Figure 3. Colors indicate the percentage of local volume change compared to the initial value in as-equilibrated 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

systems. For visualizing clarification, the color scale bar on the right varies from 0% to 30%; however, the volume of any group of atoms can change by more than 30%. The compression is along the z axis.

Figure 5. The dislocation motion inside the defected MOF-5 under compression. The layers a, b, and c correspond to the letters in Figure 4 and are colored according to their position along z axis.

Figure 6. Illustration of (a) a screw dislocation, (b) an edge dislocation in perfect lattice and (c) the edge dislocation in lattice with point vacancy.

Figure 7. Structure evolution of defected MOF-5 under shear. The layers a, b, c and d in bottom were colored according to their position along y axis.

Figure 8. Comparison of the stretching and the compressive deformation in defected MOF-5: (a) Stress−strain relations; (b) the variation of lateral size.

Figure 9. Atomic structures of defected MOF-5 under various tensile strains corresponding to the values identified by the letters in Figure 8a. Colors indicate the percentage of local volume change, compared to the initial value in as-equilibrated systems. The linkages between the nearest neighboring nodes were broken if the length was larger than 20% of the equilibrium value. The stars denote the hexagonal 20

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

units in the framework.

TOC Graphic

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1 56x26mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2 43x15mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3 35x12mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4 102x81mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5 103x67mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6 45x13mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7 84x47mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8 48x19mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 9 102x65mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC graphic 73x110mm (300 x 300 DPI)

ACS Paragon Plus Environment