Sang Soo Han

Tuning Metal–Organic Frameworks with Open-Metal Sites and Its Origin for Enhancing CO2 Affinity by Metal Substitution

Reducing anthropogenic carbon emission is a problem that requires immediate attention. Metal-organic frameworks (MOFs) have emerged as a promising new materials platform for carbon capture, of which Mg-MOF-74 offers chemospecific affinity toward CO2 because of the open Mg sites. Here we tune the binding affinity of CO2 for M-MOF-74 by metal substitution (M = Mg, Ca, and the first transition metal elements) and show that Ti- and V-MOF-74 can have an enhanced affinity compared to Mg-MOF-74 by 6-9 kJ/mol. Electronic structure calculations suggest that the origin of the major affinity trend is the local electric field effect of the open metal site that stabilizes CO2, but forward donation from the lone-pair electrons of CO2 to the empty d-levels of transition metals as in a weak coordination bond makes Ti and V have an even higher binding strength than Mg, Ca, and Sc.

ReaxFF Reactive Force Field for the Y-Doped BaZrO3 Proton Conductor with Applications to Diffusion Rates for Multigranular Systems

Proton-conducting perovskites such as Y-doped BaZrO 3 (BYZ) are promising candidates as electrolytes for a proton ceramic fuel cell (PCFC) that might permit much lower temperatures (from 400 to 600 degrees C). However, these materials lead to relatively poor total conductivity ( approximately 10 (-4) S/cm) because of extremely high grain boundary resistance. In order to provide the basis for improving these materials, we developed the ReaxFF reactive force field to enable molecular dynamics (MD) simulations of proton diffusion in the bulk phase and across grain boundaries of BYZ. This allows us to elucidate the atomistic structural details underlying the origin of this poor grain boundary conductivity and how it is related to the orientation of the grains. The parameters in ReaxFF were based entirely on the results of quantum mechanics (QM) calculations for systems related to BYZ. We apply here the ReaxFF to describe the proton diffusion in crystalline BYZ and across grain boundaries in BYZ. The results are in excellent agreement with experiment, validating the use of ReaxFF for studying the transport properties of these membranes. Having atomistic structures for the grain boundaries from simulations that explain the overall effect of the grain boundaries on diffusion opens the door to in silico optimization of these materials. That is, we can now use theory and simulation to examine the effect of alloying on both the interfacial structures and on the overall diffusion. As an example, these calculations suggest that the reduced diffusion of protons across the grain boundary results from the increased average distances between oxygen atoms in the interface, which necessarily leads to larger barriers for proton hopping. Assuming that this is the critical issue in grain boundary diffusion, the performance of BYZ for multigranular systems might be improved using additives that would tend to precipitate to the grain boundary and which would tend to pull the oxygens atoms together. Possibilities might be to use a small amount of larger trivalent ions, such as La or Lu or of tetravalent ions such as Hf or Th. Since ReaxFF can also be used to describe the chemical processes on the anode and cathode and the migration of ions across the electrode-membrane interface, ReaxFF opens the door to the possibility of atomistic first principles predictions on models of a complete fuel cell.

High H2 Uptake in Li-, Na-, K-Metalated Covalent Organic Frameworks and Metal Organic Frameworks at 298 K

The Yaghi laboratory has developed porous covalent organic frameworks (COFs), COF102, COF103, and COF202, and metal-organic frameworks (MOFs), MOF177, MOF180, MOF200, MOF205, and MOF210, with ultrahigh porosity and outstanding H(2) storage properties at 77 K. Using grand canonical Monte Carlo (GCMC) simulations with our recently developed first principles based force field (FF) from accurate quantum mechanics (QM), we calculated the molecular hydrogen (H(2)) uptake at 298 K for these systems, including the uptake for Li-, Na-, and K-metalated systems. We report the total, delivery and excess amount in gravimetric and volumetric units for all these compounds. For the gravimetric delivery amount from 1 to 100 bar, we find that eleven of these compounds reach the 2010 DOE target of 4.5 wt % at 298 K. The best of these compounds are MOF200-Li (6.34) and MOF200-Na (5.94), both reaching the 2015 DOE target of 5.5 wt % at 298 K. Among the undoped systems, we find that MOF200 gives a delivery amount as high as 3.24 wt % while MOF210 gives 2.90 wt % both from 1 to 100 bar and 298 K. However, none of these compounds reach the volumetric 2010 DOE target of 28 g H(2)/L. The best volumetric performance is for COF102-Na (24.9), COF102-Li (23.8), COF103-Na (22.8), and COF103-Li (21.7), all using delivery g H(2)/L units for 1-100 bar. These are the highest volumetric molecular hydrogen uptakes for a porous material under these thermodynamic conditions. Thus, one can obtain outstanding H(2) uptakes with Li, Na, and K doping of simple frameworks constructed from simple, cheap organic linkers. We present suggestions for strategies for synthesis of alkali metal-doped MOFs or COFs.

Fine-Tuning of the Carbon Dioxide Capture Capability of Diamine-Grafted Metal-Organic Framework Adsorbents Through Amine Functionalization.

A combined sonication and microwave irradiation procedure provides the most effective functionalization of ethylenediamine (en) and branched primary diamines of 1-methylethylenediamine (men) and 1,1-dimethylethylenediamine (den) onto the open metal sites of Mg2 (dobpdc) (1). The CO2 capacities of the advanced adsorbents 1-en and 1-men under simulated flue gas conditions are 19 wt % and 17.4 wt %, respectively, which are the highest values reported among amine-functionalized metal-organic frameworks (MOFs) to date. Moreover, 1-den exhibits both a significant working capacity (12.2 wt %) and superb CO2 uptake (11 wt %) at 3 % CO2 . Additionally, this framework showcases the superior recyclability; ultrahigh stability after exposure to O2 , moisture, and SO2 ; and exceptional CO2 adsorption capacity under humid conditions, which are unprecedented among MOFs. We also elucidate that the performance of CO2 adsorption can be controlled by the structure of the diamine ligands grafted such as the number of amine end groups or the presence of side groups, which provides the first systematic and comprehensive demonstration of fine-tuning of CO2 uptake capability using different amines.

Rapid Dye Regeneration Mechanism of Dye-Sensitized Solar Cells

During the light-harvesting process of dye-sensitized solar cells (DSSCs), the hole localized on the dye after the charge separation yields an oxidized dye, D(+). The fast regeneration of D(+) using the redox pair (typically the I(-)/I3(-) couple) is critical for the efficient DSSCs. However, the kinetic processes of dye regeneration remain uncertain, still promoting vigorous debates. Here, we use molecular dynamics simulations to determine that the inner-sphere electron-transfer pathway provides a rapid dye regeneration route of ∼4 ps, where penetration of I(-) next to D(+) enables an immediate electron transfer, forming a kinetic barrier. This explains the recently reported ultrafast dye regeneration rate of a few picoseconds determined experimentally. We expect that our MD based comprehensive understanding of the dye regeneration mechanism will provide a helpful guideline in designing TiO2-dye-electrolyte interfacial systems for better performing DSSCs.

First-Principles Design of Hydrogen Dissociation Catalysts Based on Isoelectronic Metal Solid Solutions

We report an innovative route for designing novel functional alloys based on first-principles calculations, which is an isoelectronic solid solution (ISS) of two metal elements to create new characteristics that are not native to the constituent elements. Neither Rh nor Ag exhibits hydrogen storage properties, whereas the Rh50Ag50 ISS exhibits properties similar to Pd; furthermore, Au cannot dissociate H2, and Ir has a higher energy barrier for the H2 dissociation reaction than Pt, whereas the Ir50Au50 ISS can dissociate H2 in a similar way to Pt. In the periodic table, Pd is located between Rh and Ag, and Pt is located between Ir and Au, leading to similar atomic and electronic structures between the pure metals (Pd and Pt) and the ISS alloys (Rh50Ag50 and Ir50Au50). From a practical perspective, the Ir-Au ISS would be more cost-effective to use than pure Pt, and could exhibit catalytic activity equivalent to Pt. Therefore, the Ir50Au50 ISS alloy can be a potential catalyst candidate for the replacement of Pt.

Atomistics of the lithiation of oxidized silicon (SiOx) nanowires in reactive molecular dynamics simulations

Although silicon oxide (SiOx) nanowires (NWs) are recognized as a promising anode material for lithium-ion batteries (LIBs), a clear understanding of their lithiation mechanism has not been reported yet. We elucidate the lithiation mechanism of SiOx NWs at the atomic scale based on molecular dynamics (MD) simulations employing the ReaxFF reactive force field developed through first-principles calculations. SiOx NWs with crystalline Si (c-Si) core and amorphous SiO2 (a-SiO2) shell structures of ∼1 nm in thickness show smaller volume expansion than pristine Si NWs, as found in previous experiments. Lithiation into SiOx NWs creates two interfaces: c-Si/a-LixSi and a-LixSi/a-LiySiO2. The mobility of the latter, which is located farther toward the outside of the NW, is slower than that of the former, which is one of the reasons why the thin SiO2 layer can suppress the volume expansion of SiOx NWs during lithiation. Another reason can be found from the stress distribution, as the SiOx NWs show stress distribution different from the pristine case. Moreover, the lithiation of SiOx NWs leads to the formation of Li2O and Li4SiO4 compounds in the oxide layer, where several Li atoms (not a majority) in Li4SiO4 can escape from the compound and diffuse into the c-Si, in contrast to the Li2O case. However, Li atoms that pass through the SiO2 layer penetrate into the c-Si preferentially along the 〈110〉 or 〈112〉 direction, similar to the mechanism observed in pristine Si NWs. We expect that our comprehensive understanding of the lithiation mechanism of SiOx NWs will provide helpful guidance for the design of SiOx anodes to obtain better performing LIBs.

Simulation Protocol for Prediction of a Solid-Electrolyte Interphase on the Silicon-based Anodes of a Lithium-Ion Battery: ReaxFF Reactive Force Field

We propose the ReaxFF reactive force field as a simulation protocol for predicting the evolution of solid-electrolyte interphase (SEI) components such as gases (C2H4, CO, CO2, CH4, and C2H6), and inorganic (Li2CO3, Li2O, and LiF) and organic (ROLi and ROCO2Li: R = -CH3 or -C2H5) products that are generated by the chemical reactions between the anodes and liquid electrolytes. ReaxFF was developed from ab initio results, and a molecular dynamics simulation with ReaxFF realized the prediction of SEI formation under real experimental conditions and with a reasonable computational cost. We report the effects on SEI formation of different kinds of Si anodes (pristine Si and SiOx), of the different types and compositions of various carbonate electrolytes, and of the additives. From the results, we expect that ReaxFF will be very useful for the development of novel electrolytes or additives and for further advances in Li-ion battery technology.

Correction and Addition to "Tuning Metal-Organic Frameworks with Open-Metal Sites and Its Origin for Enhancing CO2 Affinity by Metal Substitution".

T geometries referred to in the paper were based on PBE optimization with the binding energies refined using single-point PBE-D2 on the PBE geometries. In the caption of Figure 1 and Table 1, however, it was incorrectly stated that the results were based on the PBE-D2 geometries. The binding energy for V-MOF-74 was also incorrectly cited to be 47.3 kJ/ mol in Table 1, but the correct binding energy is 52.4 kJ/mol. In addition, we have performed the full PBE-D2 optimizations. The full PBE-D2 optimization shifts the previous single-point PBE-D2 binding energies using the PBE geometries up by approximately 1−3 kJ/mol and shortens the M− OCO distances consistently by 0.1−0.2 Å. These new full PBED2 optimizations lead to the CO2 binding energies that are even closer to the experiments by Caskey et al. (J. Am. Chem. Soc. 2008, 130, 10870−10871) compared to the PBE-D2//PBE results cited in our original version. We now added a new table in the Supporting Information (Table S5) that compares the PBE and PBE-D2 optimized results. In summary, the correction and addition described above do not change any of the conclusions presented in the original version but rather strengthen one of our technical points, the validity of PBE-D2 for the MOF-74/CO2 system. ■ ASSOCIATED CONTENT *S Supporting Information Calculation details, calculated structural parameters, electrostatic energy estimations, DOS analyses for all metals before and after the binding of CO2, charge densities for showing orbitals interactions, Bayder charge analysis, forward donation analysis, and comparison of the PBE and PBE-D2 optimized results are described. This material is available free of charge via the Internet at http://pubs.acs.org.

Activity, Selectivity, and Durability of Ruthenium Nanoparticle Catalysts for Ammonia Synthesis by Reactive Molecular Dynamics Simulation: The Size Effect

We report a molecular dynamics (MD) simulation employing the reactive force field (ReaxFF), developed from various first-principles calculations in this study, on ammonia (NH3) synthesis from nitrogen (N2) and hydrogen (H2) gases over Ru nanoparticle (NP) catalysts. Using ReaxFF-MD simulations, we predict not only the activities and selectivities but also the durabilities of the nanocatalysts and discuss the size effect and process conditions (temperature and pressure). Among the NPs (diameter = 3, 4, 5, and 10 nm) considered in this study, the 4 nm NPs show the highest activity, in contrast to our intuition that the smallest NP should provide the highest activity, as it has the highest surface area. In addition, the best selectivity is observed with the 10 nm NPs. The activity and selectivity are mainly determined by the hcp, fcc, and top sites on the Ru NP surface, which depend on the NP size. Moreover, the selectivity can be improved more significantly by increasing the H2 pressure than by increasing the N2 pressure. The durability of the NPs can be determined by the mean stress and the stress concentration, and these two factors have a trade-off relationship with the NP size. In other words, as the NP size increases, its mean stress decreases, whereas the stress concentration simultaneously increases. Because of these two effects, the best durability is found with the 5 nm NPs, which is also in contrast to our intuition that larger NPs should show better durability. We expect that ReaxFF-MD simulations, along with first-principles calculations, could be a useful tool in developing novel catalysts and understanding catalytic reactions.