Open Access Article
Kyle R.
Langlois
ab,
Yangyang
Wang
c,
Fnu
Joshua
ad,
Wesley J.
Newsome
a,
Lei
Zhai
ad,
James K.
Harper
e and
Fernando J.
Uribe-Romo
*ab
aDepartment of Chemistry, University of Central Florida, Orlando, FL 32814, USA. E-mail: fernando@ucf.edu
bREACT: Renewable Energy and Chemical Transformations Cluster, University of Central Florida, Orlando, FL 32814, USA
cCenter for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
dNanoScience Technology Center, University of Central Florida, Orlando, FL 32814, USA
eDepartment of Chemistry, University of Utah, Salt Lake City, UT 84112, USA
First published on 10th June 2025
In this work, we show how the combination of soluble non-polar and polar links allows for the preparation of multivariate metal–organic frameworks (MTV MOFs) that exhibit dipolar solid-solution behavior. We prepared a library of PIZOF-2 MOFs with varied concentrations of p-nitroaniline (PNA)-containing moiety embedded within the MTV links. This library forms a partial solid solution up to 30 mol% with input/output composition ratio of 0.536 ± 0.018. MOFs with x > 5 mol% PNA show concentration-dependent hypsochromic dipole–dipole coupling (H-coupling) in the solid-state UV-visible spectra. The MTV library also exhibits a dielectric relaxation event in the broadband dielectric spectra with activation enthalpies and entropies that vary with temperature and PNA content and follow linear Meyer–Neldel compensation relations.
In this work, we show the incorporation of molecular π-dipoles as multivariate solute links in a PIZOF matrix, producing a library of crystals that show concentration dependent effects (Fig. 1a). The π-dipoles correspond to p-nitroaniline (PNA) incorporated within the central ring of the matrix link (Fig. 1b). PNA is the archetypal example of a dipolar molecule, featuring a donor–π–acceptor (or push/pull) electronic structure that produces strong polarity,14,15 and it has been extensively studied in liquid solutions, polymer blends, and porous solids.16–24 We hypothesized that systematic incorporation of PNA building blocks as solute components in MTV MOF matrices resulted in MOFs that have variable bulk polarity. The multivariate crystals were characterized using solid-state optical and dielectric spectroscopies, observing concentration-dependent dipole–dipole hypsochromic coupling and dielectric relaxation with thermodynamic parameters of activated state that depend on PNA content. We demonstrate how targeting molecular traits in solute links enable the study of their solution-like behaviour in the solid state.
![]() | ||
| Fig. 1 (a) MTV MOF library with varied concentration of dipolar links in the MOF matrix. (b) Synthesis of a dipolar MTV MOF library with p-nitroaniline as solute dipole. | ||
m). Under slightly different synthesis conditions, we crystalized the phase with 100 mol% PEPEP-PNA. It is important to note the absence of the 200 reflections, this absence is characteristic of the interpenetrated nature of the MOF (section S3†).3 Samples observed under scanning electron microscopy (SEM) showed octahedral crystals with particle sizes ranging between 1 to 80 μm across all compositions, also exhibiting an overall small particle size distribution (7.4 to 18.8 μm, see section S6†).
Output compositions were determined via solution 1H NMR spectroscopy from samples digested in D2SO4/DMSO-d6. Integration of the respective peaks for each link provided input/output (I/O) molar ratios shown in Fig. 2b. In an ideal solid-solution the I/O molar ratios should be one, with deviations due to lower miscibility of the solute links in the matrix. We observed a I/O slope of m = 0.536 ± 0.018 with good linear correlation (R2 = 0.9906), for ten compositions, in triplicate, up to 50 mol% input, evidencing the formation of a partial solid-solution up to ca. 30 mol% output. Having a slope less than one indicates that PEPEP-PNA links may not be as miscible in the PEPEP-Me matrix, compared to other MTV libraries.13 Size effects of this matrix-solute immiscibility were noticed when attempting to crystallize a multivariate library using quaterphenyl-dicarboxylate (QPDC) matrix link (PIZOF-1 MOF). In a previous communication,13 we showed that links of different size may form solid-solutions that follow Vegard's law. However, in the case of QPDC and PEPEP-PNA, we observed phase separation, inconsistent solute loadings, and crystal indexing, indicating lower tolerance of PEPEP-PNA in the QPDC MOF matrix.
The molecular structure of the MTV library was assessed using vibrational and solid-state NMR spectroscopies (see ESI section S7†). Solid-state 13C CP MAS NMR spectra (Fig. S43†) displays the presence of the expected chemical signals for carboxylates (ca. 155 ppm), aromatic (150–110 ppm), alkynes (70 ppm), methyl (5 ppm), and some residual THF guest (55, 20 ppm). Similarly, Raman spectra of the multivariates (Fig. S45†) show the presence of symmetric alkyne stretch at ca. 2200 cm−1 and carboxylate at ca. 1600 cm−1. All the prepared MOFs are porous, exhibiting Brunauer–Emmett–Teller (BET) surface areas in the range SBET = 820–1450 m2 g−1. The pore size distribution plots obtained from non-local density functional theory (NLDFT) methods exhibited a change in average pore diameter from 22.5 Å in the 0 mol% MOF matrix to 18 Å in the 100 mol% PEPEP-PNA MOF. This shift can relate to pore crowding by amino and nitro groups. The mixed-link multivariates display both pore sizes, with a small contribution by 18 Å pores and large contribution of 22.5 Å pores, consistent with the trend of PNA content. All the prepared MOFs are thermally stable (under N2) up to around 300–400 °C as determined by thermal gravimetry.
The photo physical properties were probed using solid-state diffuse reflectance UV-visible (DR-UV-vis) spectroscopy. DR-UV-vis measures changes in light reflected by powder samples due to absorption and scattering. A common challenge in the study of powder solids, is to decouple absorption from scattering because it enables the use of solution based theories. Given that we observe similar crystallite size distribution in the SEM (see above), we expect to see minimal or average effect of scattering. According to the theory of light absorption by chromophores in solution, varying the concentration of a solute induces continuous changes in reflectance due to continuous changes in absorption. The changes in absorption are evidenced in the Kubelka–Munk function F(R(λ)) of the MOFs (1):
![]() | (1) |
The solution-state spectrum of the PNA diester link in 1,2-dichloroethane displays at least three absorption bands below 400 nm with λmax = 350 nm, related to π → π* transitions in the dialkyne-PNA electronic structure.26 The λmax signal appears in MTV samples with PEPEP-PNA concentration up to 3.2 mol% (40 mM based on the average unit cell). At higher PEPEP-PNA concentrations, an intense signal at 298 nm appears, displaying significant broadening with increased PNA content and extending up to 700 nm in the 100 mol% sample (1250 mM based on links per unit cell, see section S4†), consistent with broadening due to increased chromophore concentration. Due to the centrosymmetric nature of the MOF matrix (space group Fd
m), it is possible for two dipoles to interact through space and inversion symmetry via hypsochromic dipole–dipole coupling (H-coupling).27–30 In the crystal structure of PIZOF-2, the central ring of the links occupy the 48f Wyckoff position with local symmetry mm2 (C2v). Since the PNA moiety in the PEPEP link has only m symmetry (Cs), positional disorder in the average unit cell induces equal distribution of –NH2 and –NO2 groups in four available positions. This results in two possible geometric configurations of the ground state, each able to produce coupling. With respect to inversion symmetry, two dipoles can interact in a parallel configuration (Fig. 3, right), with Γu representation, and antiparallel with Γg representation.31
In hypsochromic dipole–dipole coupling, a pair of chromophores in either Γg or Γu ground states, interact by producing mixed excited states that arise from the tensor product space of two chromophores, each contributing with a ground and a first excited state as basis vectors. In either geometry, the tensor product space results in two possible mixed states, ∣Ψg〉 and ∣Ψu〉 (Fig. 3, right), of which only the high energy mixed state (∣Ψg〉 in parallel, ∣Ψu〉 in antiparallel) is allowed by symmetry to be accessed from the ground state. This allowed high energy transition results in an intense signal in the F(R) spectrum of the MTV MOFs with x > 5 mol%, with λmax = 298 nm and significant broadening at increased PNA concentration. This signal is shifted by about 32 nm from the main PEPEP-PNA transition a 330 nm. Tauc analysis of the H-band in MOFs with x > 5 mol% produced an optical energy gap of EH = 3.40 ± 0.07 eV over six compositions. From the concentration of PEPEP-PNA in molar units, (mollink dmMOFcrystal−3), the average distance between dipoles correspond to a minimum average distance for detection of H-coupling to be between 28–35 Å (Table S4†). Further evidence of H-coupling was obtained from fluorescence spectroscopy. All MTV samples displayed low solid-state fluorescence (Fig. S26†), and a similar profile with λmax = 539 ± 2 nm. The 0 mol% matrix emmits with a λmax = 519 nm, and the 100 mol% PNA MOF at λmax = 560 nm. H-coupling is characterized by displaying low fluorescence and no fluorescence shift vs. concentration. These observations suggest the MTV library emits from the low-energy mixed state, with no fluorescence shift; in contrast to J-coupling, which shifts the λmax continuously with respect to composition.31 Future studies include transient absorption and lifetime measurements on larger single crystalline samples to further study the proposed electronic structure.
The dielectric properties of the PNA MTV library were measured in the frequency range of 100–106 Hz and temperature range of 163 and 293 K using a Novocontrol broadband dielectric spectrometer. The temperature was controlled by a QUATRO cryogenic system (ESI† section S10). The activated MOF samples were loaded in a custom-made plate–plate capacitor cell. In this work, we focus on analyzing the frequency-dependent complex dielectric permittivity (2):
(f) = ε′(f) − iε′′(f) | (2) |
The frequency-dependent dielectric loss ε′′(f) of all samples (Fig. 4c and S86–S92†), show the appearance of a broad and weak dielectric relaxation signal that peaks at 100 Hz at around 190 K and drifts to ca. 102 Hz at 200 K. Above this temperature, the signal drifts towards 106 Hz at 293 K, displaying a noticeable change in drift rate with respect to the lower temperatures (see white arrows in Fig. 4c). This dielectric behavior is similar to that observed in other simple MOFs,34–37 where the dielectric response is associated with torque around the p-phenylene rings of the links as result of quadrupole-electric field gradient interaction. Fitting the functions with the Havriliak–Negami relation for dielectric relaxation provided the characteristic relaxation time (τ). Due to the bulk powder nature of the prepared MTV MOF library, accurate determination of the induced polarization (ε∞) and dielectric strength (Δε) of the dipoles was not possible.
The dielectric relaxation process follows a thermally activated mechanism as evidenced by the appearance of linear functions in the Eyring plot of relaxation time vs. temperature. Fig. 4d shows the Eyring plot for 22 mol% PNA MOF, exhibiting a linear trend that changes slope at around 200 K, consistent with the observation of signal drift. A change in slope of this kind suggest: a) the appearance of a collective mode of activation at high temperatures, or b) the presence of two relaxation processes that appear at different temperatures. Extraction of the slopes and intercepts for each linear section (LT below 200 K, HT above 200 K) allowed the determination of activation enthalpies, ΔH‡, and entropies, ΔS‡, for each composition (Fig. S95†). The effects of PNA concentration on the dielectric properties of the MTV library can be observed in Fig. 4e (Table S8†). In LT region, the activation enthalpies range between 30–50 kJ mol−1. MTV composition induces an initial drop in activation enthalpy, from 50 kJ mol−1 in 0 mol% matrix to 33 kJ mol−1 in both 5 and 12 mol% samples, showing a slight steady increase in enthalpy to 44 and 45 kJ mol−1 in 27 and 100 mol%, respectively. The HT region ranges between 60–100 kJ mol−1 with similar behavior vs. PNA content, showing a decrease from 80 kJ mol−1 (0% mol MOF) to 60 kJ mol−1 (5 mol% MOF), with also an increase to 80 kJ mol−1 (27 mol%) and 98 kJ mol−1 100 mol%). Further evidence on the activated nature of the relaxation and the entropy effects by multivariate composition can be observed in the Fig. 4f, which displays a linear relation between activation enthalpies and entropies for all compositions in both temperature regimes. This observation is consistent with the Meyer–Neldel compensation rule,38 which states that in thermally activated processes between samples of similar composition, changes in activation energy are proportional to the possible number of pathways to achieve the activated state. This trend has been observed in other small molecules and polymer blends,39,40 and evidences the thermodynamic effects of the MOF solid-solutions.
The average unit cell of the PIZOF-2 MOF (CCDC code OXOLET)3 in Fd
m symmetry, shows positional and dynamic disorder around the phenylene rings in the MOF. Furthermore, we have observed variable rates of phenylene ring rotation depending on their location in the unit cell.41 A crystal model of 100 mol% PNA PIZOF in R3 space group symmetry (Fig. S99†) exhibits steric crowding around the central ring caused by nitro and amino groups, perhaps inhibiting rotational degrees of freedom. A similar situation occurs in the 0 mol% MOF matrix (F4132, Fig. S100†), as the methyl groups also may induce steric crowding reducing molecular motion. Qualitatively, the relative low values of ε′ and ε′′ suggest little permittivity and little losses, indicating that perhaps PNA does not respond to alternating fields due to: (a) dipolar relaxation modes are outside the frequency range, or (b) PNA moiety is not as dynamic in the crystal, and we only observe rotational dynamics of the benzoate phenylene rings and/or methyl groups. Future studies include assessment dielectric spectroscopy studies in larger single-crystalline samples for dielectric strength and dipole moment elucidation.
Footnote |
| † Electronic supplementary information (ESI) available: Synthesis, powder crystallography, materials characterization, optical and dielectric spectra, CCDC 2455083. See DOI: https://doi.org/10.1039/d5ce00299k |
| This journal is © The Royal Society of Chemistry 2025 |