Keith McLaughlin

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.

Simulations of hydrogen sorption in rht-MOF-1: identifying the binding sites through explicit polarization and quantum rotation calculations

Grand canonical Monte Carlo (GCMC) simulations of hydrogen sorption were performed in rht-MOF-1, a metal–organic framework (MOF) that consists of isophthalate groups joined by copper paddlewheel clusters and Cu3O trimers through tetrazolate moeities. This is a charged rht-MOF that contains extra-framework nitrate counterions within the material. For the simulations performed herein, excellent agreement with experiment was achieved for the simulated hydrogen sorption isotherms and calculated isosteric heat of adsorption, Qst, values only when using a polarizable potential. Thermodynamic agreement is demonstrated via comparing to experimental isotherms and binding sites are revealed by combining simulation and inelastic neutron scattering (INS) data. Simulations involving explicit many-body polarization interactions assisted in the determination of the binding sites in rht-MOF-1 through the distribution of the induced dipoles that led to strong adsorbate interactions. Four distinct hydrogen sorption sites were determined from the polarization distribution: the nitrate ions located in the corners of the truncated tetrahedral cages, the Cu2+ ions of the paddlewheels that project into the truncated tetrahedral and truncated octahedral cages (Cu1 ions), the Cu2+ ions of the Cu3O trimers (Cu3 ions), and the sides of the paddlewheels in the cuboctahedral cage. The simulations revealed that the initial sorption sites for hydrogen in rht-MOF-1 are the nitrate ions; this site corresponds to the high initial Qst value for hydrogen (9.5 kJ mol−1) in the MOF. The radial distribution functions, g(r), about the Cu2+ ions at various loadings revealed that the Cu1 ions are the preferred open-metal sorption sites for hydrogen at low loading, while the Cu3 ions become occupied at higher loadings. The validation of the aforementioned sorption sites in rht-MOF-1 was confirmed by calculating the two-dimensional quantum rotational levels about each site and comparing the levels to the transitions that were observed in the experimental INS spectra for hydrogen in the compound. For each binding site, the rotational transitions from j = 0 to j = 1 were in good agreement to certain transitions that were observed in the INS spectra. From these calculations, the assignment of the peaks in the INS spectra for hydrogen in rht-MOF-1 has been made.

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.

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 Co2+ 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.

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.

Dramatic effect of pore size reduction on the dynamics of hydrogen adsorbed in metal–organic materials

The effects of pore size reduction on the dynamics of hydrogen sorption in metal–organic materials (MOMs) were elucidated by studying SIFSIX-2-Cu and its doubly interpenetrated polymorph SIFSIX-2-Cu-i by means of sorption, inelastic neutron scattering (INS), and computational modeling. SIFSIX-2-Cu-i exhibits much smaller pore sizes, which possess high H2 sorption affinity at low loadings. Experimental H2 sorption measurements revealed that the isosteric heat of adsorption (Qst) for H2 in SIFSIX-2-Cu-i is nearly two times higher than that for SIFSIX-2-Cu (8.6 vs. 4.6 kJ mol−1). The INS spectrum for H2 in SIFSIX-2-Cu-i is rather unique for a porous material, as only one broad peak appears at low energies near 6 meV, which simply increases in intensity with loading until the pores are filled. The value for this rotational transition is lower than that in most neutral metal–organic frameworks (MOFs), including those with open Cu sites (8–9 meV), which is indicative of a higher barrier to rotation and stronger interaction in the channels of SIFSIX-2-Cu-i than the open Cu sites in MOFs. Simulations of H2 sorption in SIFSIX-2-Cu-i revealed two hydrogen sorption sites in the MOM: direct interaction with the equatorial fluorine atom (site 1) and between two equatorial fluorine atoms on opposite walls (site 2). The calculated rotational energy levels and rotational barriers for the two sites in SIFSIX-2-Cu-i are in good agreement with INS data. Furthermore, the rotational barriers and binding energies for site 2 are slightly higher than that for site 1, which is consistent with INS results. The lowest calculated transition for the primary site in SIFSIX-2-Cu is also in good agreement with INS data. In addition, this transition in the non-interpenetrating material is higher than any of the sites in SIFSIX-2-Cu-i, which indicates a significantly weaker interaction with the host as a result of the larger pore size.

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.

Atomistic simulation of frictional anisotropy on quasicrystal approximant surfaces

J. Y. Park et al. [Science 309, 1354 (2005)] have reported eight times greater atomic-scale friction in the periodic than in the quasiperiodic direction on the twofold face of a decagonal Al-Ni-Co quasicrystal. Here we present results of molecular-dynamics simulations intended to elucidate mechanisms behind this giant frictional anisotropy. Simulations of a bare atomic-force-microscope tip on several model substrates and under a variety of conditions failed to reproduce experimental results. On the other hand, including the experimental passivation of the tip with chains of hexadecane thiol, we reproduce qualitatively the experimental anisotropy in friction, finding evidence for entrainment of the organic chains in surface furrows parallel to the periodic direction.