Open Access Article
Durga Prasad Karothu
*a,
Ibrahim Tahirab,
Sanjit Manohar Majhia,
Ejaz Ahmed
b,
Luca Catalano
c,
Niamh T. Hickey
de,
James Westonf,
Sarah Guerin
de and
Panče Naumov
*abgh
aCenter for Smart Engineering Materials, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates. E-mail: pance.naumov@nyu.edu; dpk3@nyu.edu
bSmart Materials Lab, New York University Abu Dhabi, PO Box 129188, Abu Dhabi, United Arab Emirates
cDynamic Molecular Materials Laboratory, Dipartimento di Scienze della Vita, Università degli Studi di Modena e Reggio Emilia, 41125 Modena, Italy
dDepartment of Chemical Sciences, Bernal Institute, University of Limerick, V94 T9PX, Ireland
eSSPC, The Science Foundation Ireland Research Centre for Pharmaceuticals, University of Limerick, V94 T9PX, Ireland
fCore Technology Platforms, New York University Abu Dhabi, 129188 Abu Dhabi, United Arab Emirates
gResearch Center for Environment and Materials, Macedonian Academy of Sciences and Arts, Bul. Krste Misirkov 2, MK-1000 Skopje, Macedonia
hMolecular Design Institute, Department of Chemistry, New York University, 100 Washington Square East, New York, NY 10003, USA
First published on 17th November 2025
The upper limits of the stiffness of organic crystals, a newly established class of engineering materials, are an underexplored, yet very exciting domain of the material property space, because it could potentially provide ordered materials composed of light atoms with mechanical properties comparable to those of light metals and alloys. At present, however, the discovery of mechanically stiff and hard crystals is entirely serendipitous, since the relationships between the stiffness and crystal structure remain elusive. Here, guided by the very high density of the crystals of mucic (galactaric) acid, we investigated and report that this material is an ultrastiff organic crystal, with a Young's modulus obtained by nanoindentation on the (100)/(
00) pair of faces of E = 50.25 ± 1.55 GPa and a hardness of H = 2.81 ± 0.15 GPa (n = 15). This value of the modulus, which exceeds those of very stiff crystals such as α-glycine, tartaric acid and L-threonine, is corroborated by density functional theory (DFT) calculations, which provide an expectedly even higher value of E = 68.5 GPa. The measured moduli on the other two accessible faces are also above the higher approximate limit (25 GPa) of values that are common for organic crystals, with E = 29.25 ± 0.43 GPa (n = 70) for the (001)/(00
) faces and E = 31.3 ± 0.9 GPa (n = 15) for the (010)/(0
0) faces. Comparison with other organic crystals using material property plots not only underscores the unique combination of high stiffness and high density in this material, rooted in its extensive network of intermolecular hydrogen bonds, but also confirms that it has the highest-to-date measured surface stiffness among the organic crystals. This result is expected to guide the discovery of other extraordinarily stiff organic crystals that could be of interest as mechanically robust, lightweight materials.
Molecular crystals that have high stiffness, as measured by their Young's moduli, are extremely rare. Examples of very stiff organic crystals include five amino acids or peptides with strong hydrogen bonds (α-glycine, γ-glycine, L-alanine, DL-serine, and glycylglycine),54 pointing to the strength of hydrogen bonding as one of the contributing factors. A notable example of an ultrastiff crystal was reported in 2021 with L-threonine, a compound that has an extraordinarily high Young's modulus (40.95 ± 1.03 GPa) and hardness (1.98 ± 0.11 GPa) for an organic crystal.53 The ongoing pursuit in the solid-state research community aimed at rationalizing the crystal stiffness with various structural characteristics such as bonding topology and strength benefits from systematic correlations between the mechanical properties and strength/density of strong intermolecular interactions.15 In this study, we report that crystals of mucic acid (also known as galactaric acid, Fig. 1a), are a rare example of an ultrastiff crystal. The past interest in this material, which is widely used in the chemical, pharmaceutical, cosmetic, and food industries,63 has been focused on its relevance to the synthesis of nylon, pyrone, and polyethylene furanoate.64–66 In 1982, it was reported that its crystals have an unusually high density for an organic crystal of 1.790 g cm−3,67 which exceeds the typical values for carbohydrate crystals (1.4–1.6 g cm−3). Unlike pentaric acids such as D-mannaric acid and D-glucaric acid, mucic acid has very low solubility in water. Given the possible implications of the crystal's mechanical properties, these observations prompted us to characterize this material more extensively, and here we report that it is one of the stiffest organic crystalline materials reported to date.
Analysis of the crystal structure at room temperature (296 K) using single crystal X-ray diffraction (SI Table S1) showed that the crystals are triclinic, space group P
(a = 4.9176(2) Å, b = 5.8208(2) Å, c = 6.8535(2) Å, α = 92.264(1)°, β = 94.208(1)° and γ = 93.549(1)°), with half a molecule in the asymmetric unit. The molecules interact with each other via strong hydrogen bonds involving all available functional groups, as shown in Fig. 1d–g. The hydroxyl groups are engaged in distinct parallelepiped-like hydrogen-bonded patterns with the adjacent molecules (Fig. 1e). The molecules form end-to-end chains by building on the typical carboxylic acid dimer motif across a center of symmetry (Fig. 1d). These molecular chains interact with each other via a network of strong hydrogen bonds between four hydroxyl groups.
The hydrogen bonding distance within the dimers is comparable to that of other carboxylic acids, with an H⋯O distance of 1.845(9) Å and a C⋯O distance of 2.663(1) Å. The structure suggests that the high density of mucic acid is a result of extensive hydrogen bonding and an exceptionally tight packing of hydrogen-bonded molecular chains (Fig. 1). We first utilized nanoindentation to assess the nanomechanical properties of the accessible faces of the crystal (Fig. 2 and SI Fig. S2). This method provided both stiffness (the degree to which an object resists elastic deformation under applied force) and hardness (the material's resistance to pressure or scratching by a sharp object). The crystal faces of mucic acid that were experimentally accessible for indentation were identified by the modelled Bravais–Friedel–Donnay–Harker (BFDH) morphology of the crystal based on the experimental crystal structure (SI Fig. S3). Fig. 2 shows the load-displacement curves at different indentation depths on the (001)/(00
) faces, along with images of the impressions as inspected by atomic force microscopy (AFM). The elastic modulus (E) on the widely accessible (001)/(00
) faces was found to be E = 29.25 ± 0.43 GPa for depths between 500 and 1000 nm with a total of 70 indents (Fig. 2 and SI Fig. S4). On the (010)/(0
0) faces, where the indentation was more difficult, the value was E = 31.3 ± 0.9 GPa for 15 indents across the selected depth range (SI Fig. S2). These values surpass those of many organic molecular crystals, which typically range between 10 and 25 GPa, with only 8% of the reported organic crystals having values above the upper limit.15 They are also significantly higher than the stiffness of materials considered to be stiff, which typically range between 16 and 23 GPa.15,19 When compared to other stiff organic crystals, the Young's moduli on the (001)/(00
) and (010)/(0
0) faces of the crystals of mucic acid are only lower than those of α-glycine (44.00 GPa on the (001) face),54 sucrose (35.96 GPa on the (001) face),15 L-threonine (40.95 GPa on the (001) face),53 and L-tartaric acid (43.34 GPa).19 We note that the elastic moduli reported for β-succinic acid68 and a few peptides69 were not directly comparable to the values measured here due to the different method used to determine the Young's modulus of these compounds, which could systematically overestimate the actual values. The hardness (H) of the (001)/(00
) and (010)/(0
0) faces of mucic acid was determined to be H = 1.74 ± 0.06 GPa and H = 1.85 ± 0.07 GPa, respectively; these values are higher than those of many organic crystals (typically, between 0.1 to 1.0 GPa).15,19 An attempt was then made to establish structure–mechanical property relationships by using energy framework analysis17,70 to rationalize the observed exceptional mechanical properties of the mucic acid crystals. In Fig. 2f and g, the pairwise interaction energies in the crystal structure are depicted as cylinders connecting the individual molecules. The radii of these connecting cylinders are proportional to the relative strengths of the corresponding intermolecular interactions and strong interactions correspond to the thicker and larger cylinders. The analysis revealed that the extended intermolecular hydrogen bonding network results in a tightly packed crystal structure and relatively uniform distribution of energy among these intermolecular interactions, which could be the main factors that contribute to the high stiffness and hardness of the mucic acid crystals. Consistent with several previously reported examples,17,71 the slightly higher Young's modulus observed on the (010)/(0
0) face compared to that of the (001)/(00
) face could be explained by the topology of intermolecular interactions, which point out to the interactions along the [010] crystallographic direction as being the strongest.
Notably, our energy framework analysis revealed that, among the three accessible crystallographic faces—(100), (010), and (001)—the intermolecular interaction topology (as depicted by the energy framework) shows the strongest interaction network along the [100] crystallographic direction, indicating that the most robust intermolecular interactions occur along this direction. This prompted us to determine the nanomechanical properties in that particular direction. Nanoindentation on this specific face, however, quickly proved to be challenging due to the very small available indentation area and the necessity to maintain the crystal in an upright position required to access that face. The sample tilting was minimized by affixing the crystal to the steel stage, which ultimately enabled us to indent that face. Fig. 3a shows the load-displacement curves at a selected indentation depth (1200 nm) on the (100)/(
00) face, and the inset of Fig. 3b shows an image of one of the impressions obtained by using AFM. Several indentations performed on several crystals returned values of the Young's modulus between 38 and 58 GPa. Most of these measurements and indent impressions were inspected after the indentation for any improper impressions and significant tilting, as the crystal had to remain in a straight position along its length (in some cases, very low values were obtained due to sample tilting). The elastic modulus on the (100)/(
00) face was found to be E = 50.25 ± 1.55 GPa for depths between 1000 and 1200 nm with a total of 15 indents (Fig. 3a, b and SI Fig. S5–S7). Fig. S6 and S7 show the indentation measurements performed on multiple crystals at varying penetration depths on (100)/(
00) face. It is noteworthy that, while slightly lower values were occasionally observed, nanoindentation measurements on multiple crystals—performed on the face perpendicular to the long axis—consistently demonstrated the exceptionally high stiffness of mucic acid crystals. These measurements are inherently challenging due to the crystals' anisotropic morphology, limited contact area (only a small fraction of the crystal face provides a sufficiently flat area for indentation, with the remainder exhibiting slope or surface irregularities) and potential for tilting or slippage during indentation. The obtained value is significantly higher than those reported for α-glycine (44.00 GPa on the (001) face),54 sucrose (35.96 GPa on the (001) face),15 L-threonine (40.95 GPa on the (001) face),53 and L-tartaric acid (43.34 GPa),19 and at present it is the highest stiffness reported for an organic crystal by nanoindentation. The hardness (H) of the (100)/(
00) face was determined to be H = 2.81 ± 0.15 GPa.
Due to the experimental challenges and to confirm whether this experimentally observed mechanical stiffness is inherent to the crystal structure, we employed computational approaches to predict the mechanical properties of the (100)/(
00) faces, providing an independent assessment of the intrinsic stiffness of mucic acid crystals. The calculations were performed simultaneously for the remaining two pairs of faces, in order to corroborate the experimentally obtained values. Computational methods such as periodic density functional theory (DFT) calculations have been proven to be able to provide reliable predictions of mechanical properties, even for crystal facets that are experimentally inaccessible.72,73 Fig. 3c and d shows the predicted cjj values for mucic acid, corresponding to the anisotropic elastic stiffness in the direction j. The c22 and c33 values correspond to the directional stiffness on the (010) and (001) faces respectively, showing the same trend as experimental measurements. The calculated moduli, E = 46.7 and 41.7 GPa for the (010) and (001) faces, respectively, are higher than the experimental values of 31.3 and 29.25 GPa. This is to be expected as DFT values are calculated at absolute zero, and the high plane wave cutoff energy (SI Fig. S8) required for convergence results in shorter bond lengths and a denser crystal structure relative to room temperature. The highest predicted stiffness is that on the (100) face, with an exceptionally high value of E = 68.5 GPa. The derived bulk modulus of 27.9 GPa, the Young's modulus of 38.7 GPa, and the shear modulus of 15.2 GPa confirm that mucic acid is an exceptionally stiff crystal. Compared to other molecular crystals of similar size, the longitudinal tensor coefficient values (j = 1–3) are well above the average values,16,54 with the predicted shear stiffness constants (j = 4–6), which vary between 22 and 26 GPa, being in the top 1% of molecular crystals in our recent high-throughput DFT screening.74 The lowest of all the 36 tensor components (SI Table S2), some of which are equal in magnitude due to symmetry, is c24 = c42 = 6 GPa.
Fig. 4 shows a comparison of the Young's modulus, hardness, and density between single crystals of mucic acid and other very stiff organic crystals: sucrose,15 tartaric acid,19 L-ascorbic acid,15 and L-threonine.53 As evident from Fig. 4b, the crystals of mucic acid have higher density compared to other stiff materials. Since the elastic modulus and hardness depend, among other factors, on the crystal density and direction of loading with respect to the strength of intermolecular interactions, this result prompted us to investigate the relationship between these properties and the hydrogen bonds. Even though intermolecular interactions determine the mechanical strength of the crystalline material to a significant extent, the stiffness and hardness cannot be rationalized exclusively by the presence of short intermolecular distances, as shown in Fig. 4d. We infer that, at a qualitative level, the unique molecular arrangement and strong intermolecular interactions in the structure of mucic acid contribute to the compactness of its structural packing and account for its exceptional mechanical properties.
On comparing Fig. 5a and b, it is evident that mucic acid occupies a distinct position in the mechanical property space of organic crystals. Fig. 5a illustrates a broad overview of the relationship between the Young's modulus and density across a wide range of organic materials, where mucic acid stands out due to its high stiffness relative to its high density. This position suggests that this material offers a unique balance between rigidity and weight, making it suitable for applications requiring lightweight, yet structurally robust materials. Fig. 5b narrows the focus to materials with higher Young's modulus values, further highlighting mucic acid's competitive stiffness in comparison to other high-modulus organic crystals. The clustering of mucic acid in Fig. 5b with materials of similar stiffness but varying densities underscores its potential as a strong candidate for engineering applications where high stiffness is desired without a proportional increase in density.
The comparison also emphasizes the exceptional combination of properties found in mucic acid, making it stand out among other well-known stiff organic crystals. The combined analysis of these plots reinforces the utility of this material in applications demanding both structural integrity and efficiency, and turns this and possibly other similar materials into candidates as lightweight organic substitutes for metals or alloys. A comparative summary of the Young's modulus of mucic acid with representative semicrystalline polymers, two-dimensional materials, and densely crosslinked polymers is provided in Table S4 (see the SI) to highlight its distinct mechanical characteristics relative to these material classes. We note that the values reported for mucic acid here approach those of pure aluminium (E ≈ 70 GPa) and are half that of copper (E ≈ 110 GPa). This places mucic acid among materials with significant stiffness, potentially offering advantages in applications requiring high mechanical strength while exhibiting a lower density compared to metals.
), (010)/(0
0) and (100)/(
00) faces of mucic acid crystals display some of the highest reported values for organic molecules. Analysis of the strength and orientation of interacting energies revealed relatively uniform distribution of interaction energies. Notably, the intermolecular interaction topology (energy frameworks) exhibits the strongest supramolecular columns along the [100] crystallographic direction, corroborated by DFT calculations of stiffness along the crystallographic a-axis. These factors collectively contribute to the exceptional mechanical properties of mucic acid single crystals. In the global mechanical property space, this material occupies a distinct position with its high stiffness and relatively high density. Its stiffness, comparable to that of other high-modulus organic crystals, makes it particularly well-suitable for applications requiring materials with both high mechanical strength and long-term durability.
CCDC 2378188 contains the supplementary crystallographic data for this article.75
Supplementary information (SI) is available. See DOI: https://doi.org/10.1039/d5sc05888k.
| This journal is © The Royal Society of Chemistry 2025 |