---
title: 'MSUCOF-4-FeCp: Ferrocene COF for H2 Storage'
url: https://www.emergentmind.com/topics/msucof-4-fecp
type: topic
---

# MSUCOF-4-FeCp: Ferrocene COF for H2 Storage

MSUCOF-4-FeCp is a ferrocene-functionalized covalent-organic framework proposed for room-temperature hydrogen storage within the Multi-binding Sites United in Covalent-Organic Framework (MSUCOF) design paradigm. In the reported first-principles multiscale computational study, it is constructed by incorporating ferrocene (\(\mathrm{FeCp}_2\)) moieties into IRCOF-102 and is predicted to exceed the U.S. Department of Energy ultimate light-duty vehicle hydrogen-storage targets, with total uptake values of \(18.0\) wt % and \(72.6\ \mathrm{g\ H_2\ L^{-1}}\) at \(298\ \mathrm{K}\) and \(700\ \mathrm{bar}\), and deliverable capacities of \(12.2\) wt % and \(52.2\ \mathrm{g\ L^{-1}}\) between \(700\) and \(5\ \mathrm{bar}\) [2602.14927].

## 1. Position within the MSUCOF literature

MSUCOF-4-FeCp occupies a distinct place in the MSUCOF series because it is not part of the original 2023 MSUCOF study. That earlier work investigated only three covalent-organic framework families—MSUCOF-1, MSUCOF-2, and MSUCOF-3—and seven transition-metal chelation variants, namely Co, Cu, Fe, Mn, Ni, Pd, and Pt. It did not report any material labeled “MSUCOF-4,” and it did not include any FeCp or ferrocenyl derivative [2306.10036].

This distinction matters because structural, synthetic, force-field, and adsorption data for MSUCOF-4-FeCp do not exist in the 2023 paper. No FeCp coordination environment, bond lengths, bond angles, fitted isotherm parameters, or uptake numbers for such a material are given there. The 2026 work therefore represents a separate computational development that follows the MSUCOF approach rather than a continuation of the exact materials set previously tabulated [2602.14927].

A common misconception is to treat MSUCOF-4-FeCp as though it were one of the transition-metalated frameworks in the original MSUCOF manuscript. That is incorrect: the earlier study’s GCMC protocol and QM-based force-field philosophy are relevant as methodological antecedents, but the actual ferrocene-functionalized framework is introduced only in the later study [2306.10036].

## 2. Structural design and framework chemistry

The parent framework is IRCOF-102, described as a three-dimensional boroxine-linked COF isoreticular to COF-102 and built from tetrahedral tetrakis(4′-borono-[1,1′-biphenyl]-4-yl)methane nodes connected via boroxine rings. Its topology crystallizes in space group \(I\bar{4}3d\) (No. 220) with cubic cell parameters \(a=b=c=42.955\ \text{\AA}\) and \(\alpha=\beta=\gamma=90^\circ\). The framework contains ca. \(2\ \mathrm{nm}\)-diameter pores arranged in a three-fold-rotational coordination environment around each boroxine linkage, yielding a surface area of \(5\,680\ \mathrm{m^2\ g^{-1}}\), pore volume of \(5.11\ \mathrm{cm^3\ g^{-1}}\), void fraction of \(0.90\), and framework density of \(0.06\ \mathrm{g\ cm^{-3}}\) [2602.14927].

To generate the MSUCOF-4 scaffold, the biphenyl linkers are replaced by indene-derived tetrakis(4-(4-boronoinden-7-yl)phenyl)methane, designated Linker 2. This linker embeds a cyclopentadienyl ring at each node. Condensation of Linker 2 produces MSUCOF-4, an isoreticular variant of IRCOF-102 with one Cp ring per linker but without Fe. Post-synthetic metallation with \(\mathrm{FeCl_2/NaCp}\) then yields MSUCOF-4-FeCp, in which each Cp ring within a tritopic pore region is deprotonated and coordinates one \(\mathrm{Fe^{2+}}\) center together with a second \(\mathrm{Cp^-}\) anion from NaCp [2602.14927].

The resulting ferrocene moieties are positioned at junctions of three struts, forming cooperative binding pockets. Upon FeCp incorporation, the cell expands slightly to \(a=b=c=43.07\ \text{\AA}\) while retaining the same \(I\bar{4}3d\) space group. The study also reports for MSUCOF-4-FeCp a porosity of \(S=4\,780\ \mathrm{m^2\ g^{-1}}\) and \(V_p=2.55\ \mathrm{cm^3\ g^{-1}}\). This suggests a design tradeoff in which some pore volume is sacrificed while introducing a denser distribution of targeted adsorption sites [2602.14927].

## 3. Electronic-structure and simulation methodology

The computational treatment combines fragment DFT, periodic DFT, force-field fitting, and grand-canonical Monte Carlo. For fragment DFT, eclipsed and staggered ferrocene fragments plus single \(\mathrm{H_2}\) molecules were modeled in \(500\ \text{\AA}\) supercells using the M06 functional with the pob-TZVP-rev2 Gaussian basis set in CRYSTAL23. Geometry optimizations used tight SCF convergence of \(10^{-11}\ \mathrm{eV}\), with forces converged to \(1.5\times10^{-3}\ \mathrm{eV\ \AA^{-1}}\); spin polarization and an extra-large integration grid were employed to resolve non-covalent energetics accurately [2602.14927].

For the periodic solid, the full MSUCOF-4-FeCp and IRCOF-102 unit cells were optimized using hybrid HSE06 + D3 dispersion corrections in CRYSTAL23 with the same basis-set quality and convergence criteria. These calculations confirmed the eclipsed Cp conformation in the solid state, the cell expansion associated with metallation, and a band-gap narrowing from \(4.10\ \mathrm{eV}\) in IRCOF-102 to \(3.02\ \mathrm{eV}\) in MSUCOF-4-FeCp [2602.14927].

Hydrogen–framework interactions were represented with a Morse potential,
$$
U_{ij}(r_{ij})=D_0\left[\left(1-e^{-\alpha(r_{ij}-r_0)}\right)^2-1\right],
$$
with parameters for \(\mathrm{H_2}\)–Fe, \(\mathrm{H_2}\)–C, and \(\mathrm{H_2}\)–H in ferrocene fitted in GULP to fragment DFT energy curves. The reported mean absolute errors are \(<0.5\ \mathrm{kJ\ mol^{-1}}\) [2602.14927].

GCMC simulations were carried out in Materials Studio at \(T=298\ \mathrm{K}\) and \(P=1\)–\(700\ \mathrm{bar}\) using the Metropolis algorithm with move ratios translation : rotation : insertion : deletion : regrowth \(=2:1:1:1:0.1\). Each pressure point used \(1\times10^6\) equilibration and \(3\times10^6\) production steps. Bulk \(\mathrm{H_2}\) fugacities were obtained via the van der Waals equation,
$$
\ln(f/P)=\frac{(b-a/(RT))\cdot P}{RT},
$$
with \(a=0.2476\ \mathrm{L^2\ bar\ mol^{-2}}\) and \(b=0.02661\ \mathrm{L\ mol^{-1}}\), giving \(<1\%\) deviation from Peng–Robinson at \(700\ \mathrm{bar}\) [2602.14927].

## 4. Hydrogen binding and adsorption-site hierarchy

The zero-load binding energy is defined in the study as
$$
E_b = E_{\mathrm{COF}+\mathrm{H_2}} - E_{\mathrm{COF}} - E_{\mathrm{H_2}},
$$
including electronic plus zero-point and vibrational enthalpy corrections. For MSUCOF-4-FeCp, the reported primary \(\mathrm{H_2}\) binding energies lie in the \(15\)–\(20\ \mathrm{kJ\ mol^{-1}}\) range, described as an ideal physisorption window for balancing high-pressure uptake with low-pressure release [2602.14927].

The pressure-dependent isosteric heat \(Q_{st}(P)\) remains entirely within \(7\)–\(20\ \mathrm{kJ\ mol^{-1}}\). The study interprets this range as strong enough to enhance adsorption relative to weakly interacting porous frameworks, while still permitting facile desorption under practical delivery conditions. A plausible implication is that the ferrocene installation changes the adsorption problem from one dominated by very high porosity alone to one organized around a controlled spectrum of moderate adsorption energies [2602.14927].

Interaction-energy deconvolution further resolves two main adsorption populations at high pressure: \(\mu_1 \approx -13.2\ \mathrm{kJ\ mol^{-1}}\), assigned to the Fe center, and \(\mu_2 \approx -7.7\ \mathrm{kJ\ mol^{-1}}\), assigned to weaker Cp-ring sites. At \(700\ \mathrm{bar}\), \(83\%\) of \(\mathrm{H_2}\) is reported to occupy the Fe site. This site hierarchy is consistent with the design premise that cooperative pockets centered on ferrocene can create a dominant primary adsorption motif supplemented by weaker secondary filling sites [2602.14927].

## 5. Storage performance and DOE benchmarks

The reported adsorption metrics place MSUCOF-4-FeCp above both the DOE 2025 targets of \(5.5\) wt %, \(40\ \mathrm{g\ L^{-1}}\) and the DOE ultimate targets of \(6.5\) wt %, \(50\ \mathrm{g\ L^{-1}}\) for automotive hydrogen storage [2602.14927].

| Quantity | Value | Context |
|---|---:|---|
| Total gravimetric uptake | \(18.0\) wt % | \(298\ \mathrm{K}, 700\ \mathrm{bar}\) |
| Total volumetric uptake | \(72.6\ \mathrm{g\ H_2\ L^{-1}}\) | \(298\ \mathrm{K}, 700\ \mathrm{bar}\) |
| Working capacity, gravimetric | \(12.2\) wt % | \(700 \rightarrow 5\ \mathrm{bar}\) |
| Working capacity, volumetric | \(52.2\ \mathrm{g\ L^{-1}}\) | \(700 \rightarrow 5\ \mathrm{bar}\) |
| IRCOF-102 total gravimetric uptake | \(15.4\) wt % | identical conditions |
| IRCOF-102 total volumetric uptake | \(32.4\ \mathrm{g\ L^{-1}}\) | identical conditions |

Deliverable, or working, capacity is defined between a charging pressure of \(700\ \mathrm{bar}\) and a delivery pressure of \(5\ \mathrm{bar}\). The net-uptake curves peak at \(P_{\mathrm{opt}} \approx 145\ \mathrm{bar}\), and the study reports broad positive net uptake across \(5\)–\(700\ \mathrm{bar}\). Excess isotherms are stated to show analogous trends when subtracting only pore-volume gas [2602.14927].

The same study compares MSUCOF-4-FeCp with precious-metal-functionalized MSUCOF-1-PtCl\(_2\), for which \(WC_v=41.4\ \mathrm{g\ L^{-1}}\) and \(WC_g=8.7\) wt % are quoted. On that basis, the ferrocene-functionalized framework is presented as outperforming prior MSUCOF variants while using Earth-abundant iron rather than platinum [2602.14927].

## 6. Economic interpretation, limitations, and proposed validation

A central practical claim of the 2026 work is that ferrocene functionalization provides a cost-effective alternative to precious-metal incorporation in COFs. The paper gives approximate raw-material figures of \(\sim\$0.10\ \mathrm{kg^{-1}}\) for iron versus \(>\$40\,000\ \mathrm{kg^{-1}}\) for Pt, and states that iron’s natural abundance exceeds \(8\times10^{11}\ \mathrm{t}\) ore. On this basis, it describes a \(4\times10^5\)-fold economic advantage over platinum [2602.14927].

At the same time, the status of MSUCOF-4-FeCp remains predictive rather than experimentally established in the provided record. The study recommends experimental synthesis of MSUCOF-4, followed by post-synthetic ferrocene installation and gas-adsorption measurements spanning \(5\)–\(700\ \mathrm{bar}\) at \(298\ \mathrm{K}\). It further recommends XPS and Mössbauer spectroscopy to confirm FeCp incorporation, and in situ IR or Raman spectroscopy to probe \(\mathrm{H_2}\) binding [2602.14927].

The proposed future computational directions include molecular dynamics, quantum corrections to \(\mathrm{H_2}\) binding, and extension to alternative metallocenes such as cobaltocene and ruthenocene. This suggests that MSUCOF-4-FeCp functions not only as a specific candidate material but also as a design template: high-porosity COF architectures can be reprogrammed by installing organometallic motifs that create cooperative adsorption pockets with binding enthalpies in the targeted \(15\)–\(20\ \mathrm{kJ\ mol^{-1}}\) range [2602.14927].

Source: https://www.emergentmind.com/topics/msucof-4-fecp