Investigating H2 Sorption in a Fluorinated Metal-Organic Framework with Small Pores Through Molecular Simulation and Inelastic Neutron Scattering.

Simulations of H2 sorption were performed in a metal-organic framework (MOF) consisting of Zn(2+) ions coordinated to 1,2,4-triazole and tetrafluoroterephthalate ligands (denoted [Zn(trz)(tftph)] in this work). The simulated H2 sorption isotherms reported in this work are consistent with the experimental data for the state points considered. The experimental H2 isosteric heat of adsorption (Qst) values for this MOF are approximately 8.0 kJ mol(-1) for the considered loading range, which is in the proximity of those determined from simulation. The experimental inelastic neutron scattering (INS) spectra for H2 in [Zn(trz)(tftph)] reveal at least two peaks that occur at low energies, which corresponds to high barriers to rotation for the respective sites. The most favorable sorption site in the MOF was identified from the simulations as sorption in the vicinity of a metal-coordinated H2O molecule, an exposed fluorine atom, and a carboxylate oxygen atom in a confined region in the framework. Secondary sorption was observed between the fluorine atoms of adjacent tetrafluoroterephthalate ligands. The H2 molecule at the primary sorption site in [Zn(trz)(tftph)] exhibits a rotational barrier that exceeds that for most neutral MOFs with open-metal sites according to an empirical phenomenological model, and this was further validated by calculating the rotational potential energy surface for H2 at this site.

Theoretical investigations of CO2 and CH4 sorption in an interpenetrated diamondoid metal-organic material.

Grand canonical Monte Carlo (GCMC) simulations of CO2 and CH4 sorption and separation were performed in dia-7i-1-Co, a metal-organic material (MOM) consisting of a 7-fold interpenetrated net of Co(2+) ions coordinated to 4-(2-(4-pyridyl)ethenyl)benzoate linkers. This MOM shows high affinity toward CH4 at low loading due to the presence of narrow, close fitting, one-dimensional hydrophobic channels-this makes the MOM relevant for applications in low-pressure methane storage. The calculated CO2 and CH4 sorption isotherms and isosteric heat of adsorption, Qst, values in dia-7i-1-Co are in good agreement with the corresponding experimental results for all state points considered. The experimental initial Qst value for CH4 in dia-7i-1-Co is currently the highest of reported MOM materials, and this was further validated by the simulations performed herein. The simulations predict relatively constant Qst values for CO2 and CH4 sorption across all loadings in dia-7i-1-Co, consistent with the one type of binding site identified for the respective sorbate molecules in this MOM. Examination of the three-dimensional histogram showing the sites of CO2 and CH4 sorption in dia-7i-1-Co confirmed this finding. Inspection of the modeled structure revealed that the sorbate molecules form a strong interaction with the organic linkers within the constricted hydrophobic channels. Ideal adsorbed solution theory (IAST) calculations and GCMC binary mixture simulations predict that the selectivity of CO2 over CH4 in dia-7i-1-Co is quite low, which is a direct consequence of the MOM's high affinity toward both CO2 and CH4 as well as the nonspecific mechanism shown here. This study provides theoretical insights into the effects of pore size on CO2 and CH4 sorption in porous MOMs and its effect upon selectivity, including postulating design strategies to distinguish between sorbates of similar size and hydrophobicity.

Efficient calculation of many-body induced electrostatics in molecular systems.

Potential energy functions including many-body polarization are in widespread use in simulations of aqueous and biological systems, metal-organics, molecular clusters, and other systems where electronically induced redistribution of charge among local atomic sites is of importance. The polarization interactions, treated here via the methods of Thole and Applequist, while long-ranged, can be computed for moderate-sized periodic systems with extremely high accuracy by extending Ewald summation to the induced fields as demonstrated by Nymand, Sala, and others. These full Ewald polarization calculations, however, are expensive and often limited to very small systems, particularly in Monte Carlo simulations, which may require energy evaluation over several hundred-thousand configurations. For such situations, it shall be shown that sufficiently accurate computation of the polarization energy can be produced in a fraction of the central processing unit (CPU) time by neglecting the long-range extension to the induced fields while applying the long-range treatments of Ewald or Wolf to the static fields; these methods, denoted Ewald E-Static and Wolf E-Static (WES), respectively, provide an effective means to obtain polarization energies for intermediate and large systems including those with several thousand polarizable sites in a fraction of the CPU time. Furthermore, we shall demonstrate a means to optimize the damping for WES calculations via extrapolation from smaller trial systems.

A molecular H2 potential for heterogeneous simulations including polarization and many-body van der Waals interactions.

A highly accurate aniostropic intermolecular potential for diatomic hydrogen has been developed that is transferable for molecular modeling in heterogeneous systems. The potential surface is designed to be efficacious in modeling mixed sorbates in metal-organic materials that include sorption interactions with charged interfaces and open metal sites. The potential parameters are compatible for mixed simulations but still maintain high accuracy while deriving dispersion parameters from a proven polarizability model. The potential includes essential physical interactions including: short-range repulsions, dispersion, and permanent and induced electrostatics. Many-body polarization is introduced via a point-atomic polarizability model that is also extended to account for many-body van der Waals interactions in a consistent fashion. Permanent electrostatics are incorporated using point partial charges on atomic sites. However, contrary to expectation, the best potentials are obtained by permitting the charges to take on values that do not reproduce the first non-vanishing moment of the electrostatic potential surface, i.e., the quadrupole moment. Potential parameters are fit to match ab initio energies for a representative range of dimer geometries. The resulting potential is shown to be highly effective by comparing to electronic structure calculations for a thermal distribution of trimer geometries, and by reproducing experimental bulk pressure-density isotherms. The surface is shown to be superior to other similarly portable potential choices even in tests on homogeneous systems without strong polarizing fields. The present streamlined approach to developing such potentials allows for a simple adaptation to other molecules amenable to investigation by high-level electronic structure methods.

A Polarizable and Transferable PHAST N2 Potential for Use in Materials Simulation.

A polarizable and transferable intermolecular potential energy function, potentials with high accuracy, speed, and transferability (PHAST), has been developed from first principles for molecular nitrogen to be used in the modeling of heterogeneous processes such as materials sorption and separations. A five-site (van der Waals and point charge) anisotropic model, that includes many-body polarization, is proposed. It is parametrized to reproduce high-level electronic structure calculations (CCSD(T) using Dunning-type basis sets extrapolated to the CBS limit) for a representative set of dimer potential energy curves. Thus it provides a relatively simple yet robust and broadly applicable representation of nitrogen. Two versions are developed, differing by the type of mixing rules applied to unlike Lennard-Jones potential sites. It is shown that the Waldman-Hagler mixing rules are more accurate than Lorentz-Berthelot. The resulting potentials are demonstrated to be effective in modeling neat nitrogen but are designed to also be useful in modeling N2 interactions in a large array of environments such as metal-organic frameworks and zeolites and at interfaces. In such settings, capturing anisotropic forces and interactions with (open and coordinated) metals and charged/polar environments is essential. In developing the potential, it was found that adding a seemingly redundant dimer orientation, slip-parallel (S), improved the transferability of the potential energy surface (PES). Notably, one of the solid phases of nitrogen was not as accurately represented energetically without including S in the representative set. Liquid simulations, however, were unaffected and worked equally well for both potentials. This suggests that accounting for a wide variety of configurations is critical in designing a potential that is intended for use in heterogeneous environments where many orientations, including those not commonly explored in the bulk, are possible. Testing and validation of the potential are achieved via simulations of a thermal distribution of trimer geometries compared to analogous high level electronic structure calculations and molecular simulations of bulk pressure-density isotherms across the vapor, supercritical, and liquid phases. Crystal lattice parameters and energetics of the α-N2 and γ-N2 solid phases are also evaluated and determined to be in good agreement with experiment. Thus the proposed potential is shown to be efficacious for gas, liquid, and solid use, representing both disordered and ordered configurations.

A Polarizable and Transferable PHAST CO2 Potential for Materials Simulation.

Reliable PHAST (Potentials with High Accuracy Speed and Transferability) intermolecular potential energy functions for CO2 have been developed from first principles for use in heterogeneous systems, including one with explicit polarization. The intermolecular potentials have been expressed in a transferable form and parametrized from nearly exact electronic structure calculations. Models with and without explicit many-body polarization effects, known to be important in simulation of interfacial processes, are constructed. The models have been validated on pressure-density isotherms of bulk CO2 and adsorption in three metal-organic framework (MOF) materials. The present models appear to offer advantages over high quality fluid/liquid state potentials in describing CO2 interactions in interfacial environments where sorbates adopt orientations not commonly explored in bulk fluids. Thus, the nonpolar CO2-PHAST and polarizable CO2-PHAST* potentials are recommended for materials/interfacial simulations.

Characterization of Tunable Radical Metal-Carbenes: Key Intermediates in Catalytic Cyclopropanation.

A new class of radical metal-carbene complex has been characterized as having Fischer-like orbital interactions and adjacent π acceptor stabilization. Density Functional Theory (DFT) along with Natural Bond Orbital (NBO) analysis and Charge Decomposition Analysis (CDA) has given insight into the electronics of this catalytic intermediate in an open-shell cobalt-porphyrin, [Co(Por)], system. The complex has a single bond from the metal to the carbene and has radical character with localized spin density on the carbene carbon. In addition, the carbene carbon is found to be nucleophilic and "tunable" through the introduction of different α-carbon substituents. Finally, based on these findings, rational design strategies are proposed which should lead to the enhancement of catalytic activity.