Lucia
Casali
a,
Maria
Carta
bc,
Adam A. L.
Michalchuk
ad,
Francesco
Delogu
*bc and
Franziska
Emmerling
*ae
aFederal Institute for Materials Research and Testing, Richard-Willstätter-Straße 11, 12489 Berlin, Germany. E-mail: franziska.emmerling@bam.de
bDepartment of Mechanical, Chemical and Materials Engineering, University of Cagliari, via Marengo 2, 09123 Cagliari, Italy. E-mail: francesco.delogu@unica.it
cCenter for Colloid and Surface Science (CSGI), Cagliari Research unit, Department of Chemistry, University of Florence, via della Lastruccia 3, 50019 – Sesto Fiorentino, FI, Italy
dSchool of Chemistry, University of Birmingham, B15 2TT Edgbaston, Birmingham, UK
eDepartment of Chemistry, Humboldt-Universität zu Berlin, 12489 Berlin, Germany
First published on 26th July 2024
Mechanochemistry is drawing attention from the pharmaceutical industry given its potential for sustainable material synthesis and manufacture. Scaling mechanochemical processes to industrial level remains a challenge due to an incomplete understanding of their underlying mechanisms. We here show how time-resolved in situ powder X-ray diffraction data, coupled with analytical kinetic modelling, provides a powerful approach to gain mechanistic insight into mechanochemical reactions. By using the ibuprofen–nicotinamide co-crystal mechanosynthesis as a benchmark system, we investigate the behaviour of the solids involved and identify the factors that promote the reaction. As mechanochemical mechanisms become increasingly clear, it promises to become a breakthrough in the industrial preparation of advanced pharmaceuticals.
Among the APIs in the WHO (World Health Organisation) model list of essential medicines, R/S ibuprofen (hereafter, ibuprofen) is an excellent example of how co-crystallization can markedly improve API physico-chemical properties. Ibuprofen (Fig. 1a) is an analgesic drug with a market size estimated at around 45000 MT (December 2023), growing at an annual rate of 2%.2 The poor aqueous solubility, bioavailability and thermal stability adversely affect the therapeutic efficacy of ibuprofen. Correspondingly, several strategies have been developed over the decades to overcome these issues,3–5 with co-crystallisation emerging as one of the most effective approaches.6
![]() | ||
Fig. 1 Schematic representation of ibuprofen (left) and nicotinamide (right). The star indicates a chiral carbon centre. |
The cocrystal of ibuprofen with nicotinamide (Fig. 1b) is particularly promising since it presents a 7.5-fold increase in water solubility as compared with pure ibuprofen, along with increased thermal stability.7,8 Given the advantages of the ibuprofen–nicotinamide cocrystal,9 strategies to reliably prepare this compound have been extensively investigated.10 To date, the ibuprofen–nicotinamide cocrystal has been successfully made by solution11,12 and vapor crystallization methods,13 from melts7 and even by mechanochemical methods i.e. extrusion14–16 and ball milling.17,18
Co-crystal syntheses based on mechanochemistry are of particular interest for future industrial applications. The drastically reduced use of solvents and energy, along with reactions that often result in 100% yields of single products, make mechanochemistry a sustainable and eco-friendly method.19 In fact, a recent study demonstrated how the ibuprofen–nicotinamide co-crystal can be prepared mechanochemically without waste,20 and with improved green metrics and reduced environmental impact in comparison with large-scale processes in solution (in batch) or solvent-free in continuous by HME (hot melt extrusion).21
Though mechanochemical methods are very promising, their translation to industry remains hindered by a lack in mechanistic understanding and selectivity. This is exacerbated by the fact that the kinetic and thermodynamic rules of analogous solution-based syntheses tend not to apply to mechanochemical routes. To tackle this challenge, methods for time-resolved in situ (TRIS) monitoring of mechanochemical reactions have been developed,22 paving the way to obtaining otherwise inaccessible information on intermediates or new products, as well as on reaction rates.
TRIS methods for monitoring ball milling reactions via X-ray diffraction (XRD) are particularly powerful to follow mechanochemical co-crystallization,22 and are excellent tools to probe and optimize mechanochemical reactions.23,24 Moreover, the collection of TRIS-XRD data provides access to kinetic profiles, which, when modelled analytically, offer exciting insight into the fundamental behaviors of solids under mechanochemical conditions.25,26 To date, such analyses have shown that mechanochemical kinetics derive from the complex interplay of physical parameters, including changes of the internal and surface energy, the reduction of the crystallite sizes, as well the reduction in the coherence energy of the solids.27–29
Using the synergies of TRIS-XRD and analytical kinetic modelling, we here investigate the mechanochemical synthesis of the ibuprofen–nicotinamide cocrystal and demonstrate that the transformation likely occurs via the mechanically induced melting of ibuprofen.
The first reaction that we investigated by TRIS-PXRD was the mechanochemical co-crystallisation of ibuprofen + nicotinamide performed at 20 Hz with a 4 mm stainless steel ball (Fig. 2, left). From literature,39 we did not expect these mild conditions to yield a complete reaction within 1 h, and indeed under these mild conditions, the first hour of milling was dominated by scattering from the reagents alone. Only early signs of co-crystal formation were seen to appear after 40 min milling.
Aiming to drive the reaction further, we increased the milling frequency to 50 Hz, and doubled the diameter of the milling ball (to 8 mm), Fig. 2, middle. Under these harsher conditions, the co-crystal appeared within the first few minutes of milling, indicating a significant speed-up of the reaction. However, these conditions still did not yield a complete reaction within 1 h; scattering from the reagents remained until the end.
Finally, we increased the ball size to 10 mm, keeping the frequency at 50 Hz, Fig. 2, right. Under these conditions the reaction reached completion within 1 h, providing us a complete dataset on which to perform analytical kinetic analysis (for details see Fig. SI1b, ESI†). To apply our analytical kinetic model to the experimental data, we had to extract quantitative phase information from the TRIS PXRD data.
The most robust approach to extract quantitative phase data for crystalline materials from PXRD data is via Rietveld methods. However, owing to the distortions to Bragg scattering data that are common for TRIS-PXRD datasets,32,40 along with the poor crystallinity of the product, our attempts at performing quantitative Rietveld refinements of our TRIS data were unreliable. As shown in Fig. 3, quantitative phase analysis (QPA) by Rietveld methods on profiles containing reagents was effective, with a good quality of the fitting, and reliable weight percentages of the components. In contrast, with the appearance of products, the Rietveld refinement was no longer suitable.
We subsequently performed QPA via Rietveld methods on ex situ data sets, but in doing so we encountered an unexpected behaviour. Namly, the reaction continued spontaneously once the powder was mechanically activated, even after the mill was stopped. For example, we observed that after ball milling for five minutes, the ex situ PXRD profile measured immediately after milling differed from that measured on the same sample aged for 48 h (see Fig. SI2a, ESI†). This behaviour unfortunately limited the reliability of ex situ data to track the progress of the mechanochemical reaction and required us to identify an alternative strategy for robust QPA of our TRIS-PXRD datasets.
To this end, we applied two alternative approaches to extract phase composition from our TRIS data: (1) we evaluated the variation of the area under a single, isolated Bragg peak of the product phase,41 and (2) we applied the multivariate method MCS-ALS.42 Although these methods are not as robust as QPA via Rietveld methods, reliable kinetic modelling has been performed using them,43–45 especially when non-crystalline or unknown crystalline materials are involved in the mechanochemical reactions. In a recent work, it was shown that MCS-ALS can be used to derive kinetic insights with comparable validity as from Rietveld refinements.42
As shown in Fig. 4, the two analytical methods provide consistent trends in the reaction kinetics, lending significant degree of confidence to the observed trend and the methods themselves (for details see Fig. SI2b, ESI†). In both cases, the relative total amount of reactants is plotted in a semi-logarithmic plot to emphasize the dependence of the fraction of reactants on time. When plotted in this way, the data are approximately linear, revealing exponential kinetics. We note that particle size can influence reaction kinetics in mechanochemical processes, affecting parameters such as rate constants without fundamentally altering the underlying reaction mechanism or compromising the extraction of mechanistic insights from kinetic data.46
![]() | ||
Fig. 4 Comparison between single peak (black dots) and MCS-ALS (red dots) analyses in a semi-logarithmic plot for the data collected the first 30 minutes of the reaction. |
To interpret the kinetic trends and gain deeper insight into the observed transformation, we used a kinetic model as follows. We assumed that the mechanical loading exceeds a certain threshold in a set of small volumes v* that are located at the points of contact between powder particles, or inside them. Once the volumes v* experience such critical loading conditions (CLCs), the reactants undergo effective mixing. The total volume of powder affected by CLCs, v, during an individual impact is simply the sum of the volumes v*. Under the assumption that the powder is effectively stirred and that the involvement of volumes v* in any given impact is stochastic, the statistics of v* affected by CLCs can be described analytically. If we denote with κ the ratio between v and the total volume of powder, V, inside the jar, the volume fraction of powder, χi(m), that has undergone CLCs i times after m impacts can be expressed, to a first approximation, as
χi(m) = [(κm)i/i!]exp(−κm). | (1) |
Although eqn (1) describes how volumes v* undergo CLCs, it does not consider the possible changes induced by CLCs in volumes v*. In this regard, it is reasonable to expect that specific values αi of the degree of chemical conversion in volumes v* after i CLCs can be associated to volume fractions χi(m).
Accordingly, the total degree of chemical conversion, α(m), can be expressed as
![]() | (2) |
Eqn (2) allows us to relate the volume v effectively involved in CLCs to the degree of chemical conversion. Since the total number of impacts, m, undergone by powders can be calculated as the product Nt between the impact frequency, N, and time, t, the fraction of final co-crystal can be expressed as
α(t) = 1 − exp(−κNt) = 1 − exp(−kt) | (3) |
ln[1 − α(t)] = −κNt = −kt, | (4) |
Thus, the slope of the linear plots shown in Fig. 4 is equal to k. Since k represents a measure of the effective amount of powder affected by CLCs per unit time, and then transformed, it can be regarded as the rate constant of the co-crystal formation.
The best-fitted k values are equal to about 0.056 min−1 and 0.059 min−1 for the data obtained from single-peak and multivariate analyses, respectively. The two values are close enough to suggest that we can consider a single average value equal to about 0.057 min−1. By considering that the total mass of reactants in the jar was 120 mg, the mass of powder affected by CLCs per unit time is approximately equal to 6.8 mg min−1.
Since the impact frequency, N, is determined by the details of the jar movement, which are characteristic of each ball mill, and the specific processing conditions, it is more appropriate to refer the kinetics to individual impacts.25 This can be done by estimating N and subsequently extracting κ from k. In this way, we can in principle use the quantity κ, which corresponds to the effective amount of powder affected by CLCs per single impact, to compare the kinetics of co-crystal formation obtained using different ball mills and under different processing conditions.
To this aim, we combined experimental and numerical methods to study the dynamics of the single milling ball inside the moving jar under the selected working conditions.
The ball can be expected to behave like a forced-damped oscillator that undergoes, to a first approximation, harmonic motion along the main jar axis. The amount of powder inside the jar affects the impact elasticity and, in turn, the overall ball dynamics. On a phenomenological basis, the effects of damping can be described by a coefficient of restitution that measures the average degree of elasticity for the impacts between ball and jar, and thus the amount of mechanical energy transferred from the ball to the powder. While details can be found in the ESI,† the most relevant results of numerical simulations performed with a coefficient of restitution of 0.4 are shown in Fig. 5. The sequence of vertical segments has been constructed by plotting the modulus of the velocity vector, vimp, of any detected impact as a function of the time, t, at which the impact was detected. Positive and negative vimp values correspond to impacts on the top and bottom hemispherical caps of the cylindrical jar, respectively. Approximately half of the detected impacts take place at a velocity vimp between ca. 2.0 m s−1 and 5.0 m s−1. Less severe impacts can be considered in the same way as rebounds, probably having minor effects on the compressed powder.
We then performed a DSC analysis on a physical mixture of ibuprofen and nicotinamide (Fig. 6). Along with the information on the thermal behaviour of the compounds, the thermal analysis was instrumental to estimate the effective thermal properties of ibuprofen and nicotinamide. After melting of ibuprofen (78 °C) we observe the melting of the co-crystal at around 90 °C. However, we do not observe melting of nicotinamide (ca. 130 °C). This strongly suggests that nicotinamide reacts entirely with the ibuprofen melt (see Fig. SI4 for more details, ESI†). To better understand the thermal events, we lowered the heating speed from 10 K min−1 to 1 K min−1, and we observed an increase in the co-crystal melting contribution (black circle in Fig. 6), which agrees with our hypothesis.
DSC analysis show a melting enthalpy of ca. 22.3 kJ mol−1 for ibuprofen, in good agreement literature values (25.0 kJ mol−1 to 39.5 kJ mol−1),47,48 and a heat capacity in the range of 200 J K−1 mol−1 to 300 J K−1 mol−1.47 Accordingly, melting 1 mol of ibuprofen (∼206.3 g), requires between 35.6 kJ to 55.4 kJ.47
Our above analysis shows that, per unit time, the volume of powder affected by CLCs is 6.8 mg min−1, and that the effective amount of powder affected by CLCs per single impact, κ, is ∼1.2 μg. The mass of ibuprofen affected by CLCs, and eventually involved in the formation of the co-crystal, is equal to about 0.7 μg min−1. It follows that the amount of energy needed to melt this much ibuprofen is between 0.12 mJ and 0.19 mJ. Notably, these values are well below the amount of energy transferred to powders during a single impact, (ca. 3.0 mJ to 20.0 mJ, see ESI†). The hypothesis of powder melting at impact is supported by high-speed video monitoring, which shows that powder sticks to the jar walls within the first minute of the reaction (Fig. SI5, ESI†).
Mechanical stimulation can induce both bulk amorphization through defect accumulation and surface melting at particle–particle interfaces during milling processes. While bulk amorphization has been the focus of attention in multicomponent reactions, the importance of surface melting, first observed by Bowden in the 1940s,49,50 has been largely overlooked. In this respect, our findings are particularly notable in the field of mechanochemistry. We also note that our experimental results do not indicate formation of amorphous phases (i.e. there is no notable no peak broadening or visible increase in background scattering contributions of the PXRD data). Moreover, our thermal analyses show no indication of glass transitions or recrystallisations of amorphous compositions. Thus, all our experimental data support that the only ‘fluid’ phase being generated in this system is a molten state.
Ultimately, we conclude that processing conditions are compatible with the melting of ibuprofen, which provides a mechanism for co-crystal formation through the dissolution of nicotinamide in molten ibuprofen and the subsequent nucleation of the co-crystal. Though this aspect will be subject to further study, when combined with other observations we hypothesise that ball milling is just a way of promoting a melt-mediated reaction of ibuprofen, consistent with literature data.48
Firstly, we monitored the mechanochemical reactions between ibuprofen and nicotinamide carried out at different grinding conditions to assess the optimal grinding conditions and eventually to detect the formation of transient phases or intermediates. We observed that the use of different milling frequencies and ball sizes only affected the reaction kinetics, so we selected the data set of the most optimised reaction (10 mm steel ball at 50 Hz) for kinetic modelling.
The kinetic analysis on PXRD data, along with a detailed characterization of the ball motion inside the jar, revealed that the formation of the final co-crystal was reasonably facilitated by the presence of a molten phase capable of enhancing the mixing. Further experiments, i.e. thermal analysis (DSC) and monitoring of the mechanochemical reaction with a high-speed camera, supported this suggestion. Therefore, we assumed that ball milling was an effective method for the preparation of the co-crystal as it locally promoted the melting of ibuprofen.
In conclusion, the combination of TRIS-PXRD measurements with kinetic modelling proved to be a powerful tool for the mechanistic understanding of mechanochemical reactions, as well as for their optimisation. We demonstrated that the ibuprofen–nicotinamide system was reasonably formed via the mechanically induced melting of ibuprofen. Although further investigation is required, this study on the benchmark system ibuprofen–nicotinamide can be considered a valuable protocol for understanding and subsequently scaling up mechanochemical reactions. As mechanochemistry becomes clearer, it is believed to represent a breakthrough in the industrial preparation of pharmaceuticals.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cp01457j |
This journal is © the Owner Societies 2024 |