research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoSTRUCTURAL
CHEMISTRY
ISSN: 2053-2296

Origin of phase relative stability and phase transformation in an S-ibu­pro­fen–nicotinamide cocrystal

crossmark logo

aUniversité Lille, CNRS, INRA, ENSCL, UMR 8207 – UMET – Unité Matériaux et Transformations, F-59650 Villeneuve d'Ascq, France, and bSoleil Synchrotron, L'Orme des Merisiers, Saint-Aubin, BP 48, 91192 Gif-sur-Yvette, France
*Correspondence e-mail: [email protected]

Edited by X. Wang, Oak Ridge National Laboratory, USA (Received 30 June 2025; accepted 14 October 2025; online 21 October 2025)

The crystal structure of the metastable form of S-ibu­pro­fen–nicotinamide cocrystals, C13H18O2·C6H6N2O, was solved from powder X-ray diffraction. This form was obtained by melting a molar mixture of S-ibu­pro­fen and nicotinamide at 100 °C, and then cooling. The high-resolution powder X-ray diffraction pattern of this new phase was recorded at room tem­per­a­ture using synchrotron radiation at SOLEIL Synchrotron (France). A hypothetical structure was ob­tained from the Monte-Carlo simulated annealing method and confirmed by Rietveld refinement. The symmetry is monoclinic (space group P21, No. 4) and the unit cell contains four mol­ecules, two of nicotinamide and two of S-ibu­pro­fen. Density functional theory (DFT) energy minimization simulation was performed in order to locate the H atoms. The determination of the crystallographic structure of this metastable form allowed an explanation of the main mechanisms at the origin of the relative stability of the two forms of the S-ibu­pro­fen–nicotinamide cocrystals. This also made it possible to explain the transition mechanism between the two forms with tem­per­a­ture.

1. Introduction

In recent years, the design of functional pharmaceutical mol­ecular materials by the cocrystallization technique has attracted increasing inter­est (Friščić & Jones, 2010View full citation) when other classical approaches based, for example, on salt formation or metastable polymorphs are not possible. Designing cocrystals is mainly based on the concepts used in crystal engineering (Desiraju, 2010View full citation), i.e. on the formation of supra­molecular hy­dro­gen-bonded networks via various hy­dro­gen-bond synthons. The stability of the supra­molecular organization is therefore highly dependant on the large variety of hy­dro­gen-bond strengths. Pharmaceutical cocrystals generally consist of an active pharmaceutical ingredient (API) and a coformer present in the same crystal structure (Friščić & Jones, 2010View full citation; Vishweshwar et al., 2006View full citation; Schultheiss & Newman, 2009View full citation; Brittain, 2013View full citation; Childs et al., 2009View full citation), for example, paracetamol–piperazine (Oswald et al., 2002View full citation), ibu­pro­fen–nicotinamide (Berry et al., 2008View full citation), carbamazepine–saccharin (Fleischman et al., 2003View full citation), etc. These multicom­ponent materials in the crystalline solid state have an obvious inter­est in terms of stability, but also in improving many physicochemical properties of an API, such as aqueous solubility, dissolution, hygroscopicity or bioavailability (Higashi et al., 2017View full citation; Friščić & Jones, 2010View full citation). However, the discovery and preparation of new cocrystals remain empirical, and still based on trial and error (ter Horst et al., 2009View full citation). Cocrystallization can be achieved by many different techniques (Vishweshwar et al., 2006View full citation; Karimi-Jafari et al., 2018View full citation), such as crystallization in solution, grinding, grinding assisted by a solvent, use of supercritical fluids, sono­crystallization, etc., which may lead to different crystalline polymorphs in an uncontrolled manner (Schultheiss & New­man, 2009View full citation; ter Horst et al., 2009View full citation). It is worth noting that using various cocrystallization synthesis methods should pro­mote various hy­dro­gen-bond types, and then various polymorphic or pseudopolymorphic forms of cocrystals, characterized by different degrees of stability.

Ibuprofen (IBP, C13H18O2) is a well-known API used as a nonsteroidal anti-inflammatory drug (NSAID) to reduce fever and to treat pain or inflammation. It can be found in many formulations (powders, capsules, tablets, etc.) and is listed in the World Health Organization's Essential Drugs List (WHO, 2007View full citation). Nevertheless, its low water solubility remains a challenge for oral use (Garzón & Martínez, 2004View full citation). Adding a second mol­ecule, and thus forming a cocrystal, to increase its solubility, is thus of inter­est. The crystallographic structure of S-ibu­pro­fen (S-IBP) is monoclinic (space group P21) and its lattice parameters are a = 12.456 (3), b = 8.0362 (3), c = 13.533 (4) Å and β = 112.86 (3)° [Cambridge Structural Data­base (CSD; Groom et al., 2016View full citation) refcode JEKNOC10 (Freer et al., 1993View full citation)]. The melting point of S-ibu­pro­fen is 53 °C (Dwivedi et al., 1992View full citation). The mol­ecule of S-ibu­pro­fen can be seen in Fig. 1[link].

[Figure 1]
Figure 1
Representation of the ibu­pro­fen (IBP) mol­ecule. The atoms are labelled according to the CIF file for the present structure and Berry et al. (2008View full citation).

Nicotinamide (N, C6H6N2O), one of the com­ponents of vitamin B, is currently used as a food additive. It is also one of the most popular coformers used for designing cocrystals, due to its high aqueous solubility and generally regarded as safe (GRAS) status.

Nicotinamide has a rich polymorphism with nine different polymorphs already known (Fellah et al., 2021View full citation). The crystallographic structure of nicotinamide at atmospheric pressure and ambient tem­per­a­ture is monoclinic (space group P21/c), and its lattice parameters are a = 3.975 (5), b = 15.632 (8), c = 9.422 (4) Å and β = 99.03 (7)° (NICOAM02 and NICOAM03; Miwa et al., 1999View full citation). The mol­ecule of nicotinamide can be seen in Fig. 2[link].

[Scheme 1]
[Figure 2]
Figure 2
Representation of the nicotinamide (N) mol­ecule. The atoms are labelled according to the CIF file for the present structure and Berry et al. (2008View full citation).

Cocrystals of S-ibu­pro­fen and nicotinamide (SOGLAC), and of RS-ibu­pro­fen and nicotinamide (SODDIZ) have been reported by Berry et al. (2008View full citation). The RS-ibu­pro­fen and nicotinamide cocrystal structure has also been confirmed by Alshahateet (2010View full citation). The synthesis of nicotinamide–ibu­pro­fen cocrystals provides the opportunity to combine the therapeutic effect of ibu­pro­fen with the high solubility of nicotinamide. The cocrystals of S-ibuprofen and nicotinamide obtained by Berry et al. (2008View full citation) crystallize in monoclinic sym­metry (space group P21) and the lattice parameters are a = 5.4110 (6), b = 55.883 (6), c = 11.9006 (13) Å and β = 90.004 (2)°.

Recently, we reported a metastable form of this cocrystal, with a different crystallographic structure com­pared to that obtained by Berry et al. (2008View full citation), the metastable form being obtained from melting a mixture of S-IBP and N in equimolar proportions (Guerain, Guinet et al., 2020View full citation). More specifically, cocrystals were synthesized using three different methods: (1) by milling, (2) by recrystallization from the melt and (3) by evaporation from a solution. Methods (1) and (3) lead to the already-known crystallographic structure of Berry et al. (2008View full citation), named form B. Method (2) leads to a new metastable crystallographic form of the cocrystal, named form A. After heating, form A is transformed into form B at 60 °C. Then, form B melts at 85 °C. The present article aims to resolve the structure of the new S-ibu­pro­fen–nicotinamide cocrystal (S-IBP:N, Scheme 1[link]), obtained by melting the molar mixture, and to explain the mechanisms behind the polymorphic transition A→B with tem­per­a­ture and the relative stability between forms A and B. The structure was solved ab initio from powder X-ray diffraction using a direct-space approach (simulated an­neal­ing) and refined by the Rietveld method. The positions of the H atoms were estimated from density functional theory (DFT) energy minimization.

2. Experimental

2.1. Cocrystal synthesis

S-Ibuprofen (mol­ecular weight = 206.28 g mol−1) was pur­chased from Sigma–Aldrich (lot number BCBH0229V, purity 99%). Powder X-ray diffraction and Rietveld analysis have shown that the commercial material is in the monoclinic crystalline form (Freer et al., 1993View full citation), as can be seen in Fig. S1 in the supporting information, and differential scanning calorimetry (DSC, Q1000, TA Instruments, 5 °C min−1) reveals the melting of the material at Tm = 53 °C.

Nicotinamide (mol­ecular weight = 122.12 g mol−1) was purchased from Sigma–Aldrich (lot number BCBV2931, purity >99.5%). Powder X-ray diffraction and Rietveld analysis have shown that the commercial material is in the most stable monoclinic crystalline P21/c form (Miwa et al., 1999View full citation), as can be seen in Fig. S1 in the supporting information, and differential scanning calorimetry (DSC, Q1000, TA Instruments, 5 °C min−1) reveals the melting of the material at Tm = 128 °C.

All the samples were analyzed without further purification.

The new cocrystal was obtained by the following methods.

An equimolar mixture (1:1), i.e. a mass of 61.89 mg of S-ibu­pro­fen and 36.63 mg of nicotinamide, was firstly melted at 100 °C, and secondly cooled:

– either to room tem­per­a­ture for isothermal crystallization

– or below the glass transition tem­per­a­ture, Tg = −20 °C (Guerain, Guinet et al., 2020View full citation), for non-isothermal recrystallization upon heating from 25 °C (at a rate of 0.5 °C min−1).

Each time, the experiment was carried out both under a cover slip and in a capillary. The same results were obtained for both methods. In all cases, a new crystallographic form of S-IBP:N cocrystals was obtained (Guerain, Guinet et al., 2020View full citation).

2.2. Refinement

Crystal data, data collection and structure refinement details are summarized in Table 1[link]. The powder X-ray dif­fraction patterns were measured at the Synchrotron SOLEIL, France, at the high-resolution powder diffraction beamline CRISTAL. The beamline is equipped with a 1D detector, `MYTHEN2 X'. The beamline was set up at an energy of 17 keV and then the energy was calibrated using NIST standard LaB6660a corresponding to a wavelength (λ) of 0.728302 Å. The cocrystal powder was enclosed in a glass capillary (diameter 0.5 mm) and mounted on the diffrac­tom­eter. The capillary was rotated during the experiments to improve the grain-sampling statistics. Data were collected at room temperature (∼20 °C) in the 1.2–40° 2θ range in 3 min.

Table 1
Crystallographic data, profile and structural parameters for S-ibu­pro­fen–nicotinamide cocrystal form A obtained after Rietveld refinement

Crystal data  
Chemical formula C19H24N2O3
Mr 328.4
Cell setting, space group Monoclinic, P21
Temperature (K) 293
a,b,c (Å) 27.9925 (13), 5.5286 (2), 6.0213 (2)
β (°) 93.112 (2
V3) 930.48 (7)
Z 2
F(000) 352
μ (mm−1) 0.079
Specimen shape, size (mm) Cylinder, 0.5
2θ range (°) 1.2–40.284°
   
Data collection  
Beamline CRISTAL (SOLEIL)
Specimen mounting 0.5 mm diameter Lindemann capillary
Data collection mode Transmission
Scan method Continuous scan
Radiation type Synchrotron 17 KeV, λ = 0.728302 Å
Binning size (°2θ) 0.004
   
Refinement  
R factors and goodness-of-fit R = 0.0694, Rwp(nb) = 0.0847, Rexp  = 0.007

The powder X-ray powder diffraction pattern of the new form of the cocrystal was com­pared to the most commonly used pattern databases for therapeutic materials. Coupling of the search/match functionalities of Highscore software with the CSD (Groom et al., 2016View full citation), crystallographic open database (COD; Gražulis et al., 2009View full citation) and the PDF-2 database of the Inter­national Center for Diffraction Data (ICDD) (Kabekkodu et al., 2024View full citation) confirmed that this new form was not already known.

3. Structure determination and Rietveld refinement

First, we com­pared the synchrotron pattern with the reference diffraction patterns obtained in the laboratory (Guerain, Guinet et al., 2020View full citation). These observations showed that traces of pure nicotinamide remained in the sample analysed at the synchrotron. Indeed, as can be seen in Fig. S2 (see supporting information), the most intense peaks of pure nicotinamide can be observed on the cocrystal pattern between 6.8 and 6.9° and at 12.8°. Consequently, these peaks are not considered from the indexation stage.

For the indexation, the profiles of 20 reflections with a 2θ angle lower than 10° were refined individually with the programs WinPlotr (Roisnel & Rodríquez-Carvajal, 2001View full citation) and DASH (David et al., 2006View full citation) in order to obtain their 2θ angular positions. The 2θ values of these reflections were com­puted in the programs DICVOL (Boultif & Louër, 2004View full citation), TOPAS (Markvardsen et al., 2008View full citation; Coelho, 2000View full citation) and McMaille (Le Bail, 2004View full citation), and indexed. The best solution obtained with the three software programs has a monoclinic cell with the following parameters: a = 27.9997, b = 5.5302, c = 6.0231 Å, β = 93.1097° and V = 931.27 Å3. The calculated figures of merit are M20 = 26.7 and F20 = 92.6 (de Wolff, 1968View full citation; Smith & Snyder, 1979View full citation).

The experimental powder X-ray diffraction pattern was used from 1.2 to 24°. The model was refined using the Le Bail method (Le Bail et al., 1988View full citation) within the program DASH. A pseudo-Voight function, a linear combination of a Lorentzian and a Gaussian of the same full width at half maximum (FWHM), was used to fit the Bragg peaks. This FWHM has a θ dependence according to Caglioti's law (Caglioti et al., 1958View full citation). The background was determined with a linear inter­polation between 50 points regularly distributed from 1.2 to 24°. The parameters were refined in the following order: the lattice parameters a, b and c, the zero-shift, the Caglioti profile parameters U, V and W, the mixing parameter η0 of the pseudo-Voight function and its 2θ dependence, and 50 points to define the background. Due to the high quality of the pattern, no parameters for the asymmetry of the Bragg peaks were used.

The DASH probabilistic approach to space-group determination, based on the systematic absences, was used and leads to the space group P21 (No. 4) and a unit cell containing four mol­ecules. In order to obtain a starting structural model, the simulated annealing algorithm of the program DASH was used. The starting configuration of the S-IBP mol­ecule was borrowed from the monoclinic model (Freer et al., 1993View full citation). The starting configuration of the N mol­ecule was borrowed from the monoclinic model (Miwa et al., 1999View full citation). The mol­ecules were introduced randomly in the unit cell. The restraint options used for the calculations did not modify the bond lengths and bond angles. The translation and orientation parameters of the mol­ecule in the cell, as well as the torsion angles, were defined as variables in the calculation. From this structural model, rigid-body Rietveld refinement was performed using DASH. The refinement was performed in three steps: first, the global isotropic tem­per­a­ture factor, second, the translation and orientation parameters of the mol­ecule, and third, the torsion angles. The lattice parameters and the background parameters were set free.

Because H atoms are poorly located in powder X-ray powder diffraction data, such data were not used to independently refine the H-atom positions. They were optimized by density functional theory (DFT). The structural model was minimized using periodic density functional theory with fixed-cell dispersion-corrected density functional theory (DFT-D) (Giannozzi et al., 2009View full citation, 2017View full citation). In this minimization, the positions of the heavy atoms were constrained, while the positions of the H atoms were let free. The Perdew–Burke–Ernzerhof (PBE) functional (Perdew et al., 1996View full citation) was used with projector-augmented wave pseudopotentials and the Grimme D3 correction (Grimme et al., 2010View full citation), as implemented in the pw.x executable of the Quantum Espresso program (Giannozzi et al., 2009View full citation, 2017View full citation).

Then, from the structural model and the H-atom positions, atomic coordinates were introduced in the programs MAUD (Materials Analysis Using Diffraction) (Lutterotti, 2010View full citation) and JANA2020 (Petrícek et al., 2014View full citation), in order to graphically com­pare the calculated and experimental X-ray diffraction diagram (MAUD) and to generate the most accurate and com­plete CIF possible (JANA2020). The lattice parameters are a = 27.9925 (13), b = 5.5286 (2), c = 6.0213 (2) Å, β = 93.112 (2)° and V = 930.48 (7) Å3. The final conventional R factors are R = 0.0694, Rwp(nb) = 0.0847, and Rexp  = 0.007. The experimental and calculated diffraction patterns are shown in Fig. 3[link]. Crystallographic data, profile and structural parameters are given in Table 1[link]. The atomic positions from the final Rietveld refinement and DFT optimization can be seen in Table S1 in the supporting information. A com­parison between the refined and the DFT-optimized structure can also be seen in Fig. S3.

[Figure 3]
Figure 3
Final Rietveld plot of the S-ibu­pro­fen–nicotinamide cocrystal at room tem­per­a­ture. The observed intensities are indicated by dots and solid lines represent the best-fit profile (upper trace) and the difference pattern (lower trace). The vertical bars correspond to the positions of the Bragg peaks.

The R factors are com­parable to those obtained in the case of other cocrystals. In particular, similar factors are obtained with R = 0.1097 and Rwp(nb) = 0.1547 for the RS-ibu­pro­fen–nicotinamide cocrystal, and R = 0.0582 and Rwp(nb) = 0.138 for the S-ibu­pro­fen–nicotinamide form B cocrystal. It is also com­parable to the results obtained for the carbamazepine–tartaric acid cocrystal (Guerain, Derollez et al., 2020View full citation) with the following R factors: R = 0.1207 and Rwp(nb) = 0.1725. The correlation between the calculated and experimental X-ray diagram could be better in the absence of spurious peaks, probably due to the presence of nicotinamide impurities. However, calculations involving the stable form of nicotinamide have been performed but do not significantly improve this correlation. Unfortunately, due to the rich polymorphism of nicotinamide, it is not possible to know unambiguously the polymorphic form of nicotinamide present, and it is also possible that multiple polymorphs are present, which makes the calculation difficult.

4. Discussion

According to our previous experiments (Guerain, Guinet et al., 2020View full citation), the Raman spectroscopy analysis performed on new S-IBP:N cocrystal form A clearly indicates both N—H⋯O mol­ecular associations and O—H⋯N associations. Fig. 4[link] shows the spectrum of X—H stretching vibrations (X = C, N and O) collected at 20 and −100 °C in a previous study (Guerain, Guinet et al., 2020View full citation). It clearly shows three stretching vibrations of bonds involved in inter­molecular associations via hy­dro­gen bonding, easily recognized by the positive tem­per­a­ture dependence of the Raman bands. The N—H stretching band around 3150 cm−1 is distinctive of N—H⋯O mol­ecular associations, while the tem­per­a­ture dependences of O—H stretching bands located at 3390 and 3400 cm−1 reflect O—H⋯N associations.

[Figure 4]
Figure 4
Intra­molecular C—H, N—H and O—H stretching region in form A at 20 and −100 °C (dashed line), and in Form B at 20 °C for S-IBP:N from Guerain, Guinet et al. (2020View full citation). The insert presents a zoom along the y axis of the N—H and O—H stretching region.

The crystallographic structure found here is consistent with these previous results, as can be seen in Fig. 5[link]. Indeed, for the S-IBP:N cocrystal form A obtained here, N—H⋯O associations form between nicotinamide mol­ecules through their primary amide group. A nicotinamide mol­ecule is linked to two others by N—H⋯O hy­dro­gen bonds. The NH2 group of the mol­ecule is linked to the C=O groups of another nicotinamide mol­ecule and its own C=O group is linked to the NH2 group of a third nicotinamide mol­ecule. Moreover, O—H⋯N associations bind the ibu­pro­fen and nicotinamide mol­ecules through bonds between the N atoms of the nicotinamide pyridine ring and the H and O atoms of the ibu­pro­fen carboxyl group.

[Figure 5]
Figure 5
Projection of the unit cell along the [001] direction for S-ibu­pro­fen–nicotinamide cocrystal form A and visualization of the hy­dro­gen-bond network. Colour key: O atoms red, N atoms blue, C atoms black and H atoms white.

The structure of the new S-IBP:N cocrystal form A resolved in this work can also be com­pared to the already known (Berry et al., 2008View full citation) cocrystal obtained by solvent evaporation or milling, namely form B (Guerain, Guinet et al., 2020View full citation).

The lattice parameters of S-IBP:N cocrystal forms A and B are given in Table 2[link]. Both structures exhibit a monoclinic symmetry and crystallize in the same space group P21 (No. 4). Form B has a small lattice parameter (∼5 Å), an inter­mediate one (∼12 Å) and a larger one (∼56 Å), while form A has two small lattice parameter (∼5 and 6 Å) and a larger one (∼28 Å) which is half that of form B. This results in a division by four of the unit-cell volume for form B com­pared to form A reflecting the fact that the cell of form A contains two mol­ecules of ibu­pro­fen and two mol­ecules of nicotinamide, while form B contains eight mol­ecules of each. Moreover, there is a difference in the β angle of the two forms, the cocrystal obtained in this work having a β angle greater than that of form B, which is a little greater than 90°.

Table 2
Lattice unit-cell parameters com­parisons between S-ibu­pro­fen–nicotinamide cocrystal form A and S-ibu­pro­fen–nicotinamide cocrystal form B (Berry et al., 2008View full citation)

Structure a (Å) b (Å) c (Å) β (°) V3) Symmetry Reference
S-IBP:N Form A 27.9925 5.5286 6.0213 93.112 930.48 P21 This work
S-IBP:N Form B 5.4110 55.883 11.9006 90.004 3598.5 P21 Berry et al. (2008View full citation)

These differences are due to a slight difference in the structural arrangement of the mol­ecules (see Figs. 6[link] and 7[link]).

[Figure 6]
Figure 6
Projection of 4 × 2 cells of S-IBP:N cocrystal form A along the [010] direction.
[Figure 7]
Figure 7
Projection of 4 × 1 cells of S-IBP:N cocrystal form B (Berry et al., 2008View full citation) along the [100] direction.

In both cocrystals, mol­ecules of ibu­pro­fen and mol­ecules of nicotinamide are stacked along the c axis without mirror projection or rotation from one mol­ecule to another. In cocrystal S-IBP:N form A, mol­ecules of ibu­pro­fen and mol­ecules of nicotinamide are also stacked along the b axis without mirror projection or rotation from one mol­ecule to another. Thus, in this cocrystal, lattice parameters c and b are similar in value.

In cocrystal S-IBP:N form B, mol­ecules of ibu­pro­fen and mol­ecules of nicotinamide are stacked along the c axis, as mentioned previously, and along the a axis without mirror projection or rotation from one mol­ecule to another.

The difference between the two cocrystals originates from the stacking of the mol­ecules along the larger axis (a for form A and b for form B). For cocrystal form A, along the a axis, there is an alternation of two mol­ecules of ibu­pro­fen, reversed with respect to each other, followed by two mol­ecules of nicotinamide, reversed with respect to each other, and so on. For cocrystal form B, there is also an alternation of two mol­ecules of ibu­pro­fen, reversed relative to each other, followed by two mol­ecules of nicotinamide, reversed relative to each other. However, the next two ibu­pro­fen mol­ecules are not in the same configuration as the two ibu­pro­fen mol­ecules preceding the nicotinamide mol­ecules. Indeed, there is a 180° rotation com­pared to the front one. The same is true for the nicotinamide mol­ecules, which are also rotated by 180° com­pared to the previous nicotinamide mol­ecules. Consequently, the b axis of S-IBP:N form B is twice that of the a axis of S-IBP:N form A. Thus, the structures of forms A and B are similar and differ mainly in the orientations of the benzene rings of ibu­pro­fen.

The small difference between the crystallographic lattice of form A and that of form B is in agreement with the observations made previously during differential scanning calorimetry (DSC) experiments (Guerain, Guinet et al., 2020View full citation). The transformation A→B around 50 °C [revealed by Raman spectroscopy and X-ray diffraction in Guerain, Guinet et al. (2020View full citation)] was difficult to observe in DSC (see Fig. S4 in the supporting information). Thus, it had been hypothesized that this transformation during heating was subtle, with a low energy phase transition and fairly close crystallographic cell for forms A and B. In addition, such small differences between the crystal lattices of the two forms also lead to relatively similar diffraction patterns with a certain number of peaks of forms A and B located at the same positions. This point was also observed in our previous experiments (see Fig. S5 in the supporting information).

Regarding the cocrystal stability, we previously demonstrated that form A, which has been resolved here, is a metastable form, while form B is stable (Guerain, Guinet et al., 2020View full citation). This relative stability and the transition from form A to form B can be explained by the com­parison between the structure obtained here and that obtained by Berry et al. (2008View full citation). For the structure of form A, as can be seen in Fig. 5[link], nicotinamide mol­ecules are linked together by N—H⋯O hy­dro­gen bonds. This is also true for form B, as can be seen in Fig. 8[link]. However, for form A, one mol­ecule of nicotinamide is linked to two different nicotinamide mol­ecules by N—H⋯O hy­dro­gen bonds, whereas in the case of form B, two mol­ecules of nicotinamide are linked together by dimer associations. Moreover, the O2—H2O⋯N1 hy­dro­gen bond exists for both forms, but this bond has a shorter distance between O2 and N1 for form B than for form A (2.652 versus 2.689 Å). Therefore, this bond is more energetic and stable for form B. This result is in agreement with previously published Raman data (Guerain, Guinet et al., 2020View full citation), where we showed that, when the tem­per­a­ture increases, the frequency of Raman bands O—H at 3375 cm−1 in form A is greater than in form B, reflecting weaker mol­ecular inter­actions between nicotinamide and ibu­pro­fen in form A than in form B. In this context, form A is less stable than form B, considering that hy­dro­gen bonds are responsible for the crystalline stability [the Raman spectra analysis performed and published in Guerain, Guinet et al. (2020View full citation) is given in Fig. S6 in the supporting information].

[Figure 8]
Figure 8
Projection of the unit cell along the [100] direction for S-ibu­pro­fen–nicotinamide cocrystal form B (Berry et al., 2008View full citation) and visualization of the hy­dro­gen-bond network.

Finally, form A has a single bond between the mol­ecule of ibu­pro­fen and that of nicotinamide (Fig. 5[link]), while form B has two (Fig. 8[link]). Indeed, the O1⋯H2A—N2 hy­dro­gen bond does not exist in form A (see Figs. 5[link] and 8[link]). Form B therefore has one more type of N—H⋯O hy­dro­gen bond than form A. In addition, the O—H⋯N hy­dro­gen bonds between the nicotinamide mol­ecules and the ibu­pro­fen mol­ecules are more energetic for form B than for form A. This justifies the different degrees of stability of the two forms observed previously (Guerain, Guinet et al., 2020View full citation).

The polymorphic transition A→B could therefore be explained by the fact that in form A, under the effect of mol­ecular agitation, rotation of the –CONH2 group of nicotinamide would lead to the formation of a dimer on the one hand and to the appearance of a new N—H⋯O hy­dro­gen bond between the ibu­pro­fen mol­ecule and that of nicotinamide on the other hand. This new configuration associated with a shorter O—H⋯N hy­dro­gen bond would ensure a stronger inter­action between the nicotinamide and ibu­pro­fen mol­ecules, and would stabilize the structure into form B, the most stable form. However, van der Waals bonds also play an important role in the stability of cocrystals (Cruz-Cabeza et al., 2024View full citation) and it would be inter­esting to carry out calculations involving them to confirm the stated hypothesis.

The energies of forms A and B as determined by the present study and in Berry et al. (2008View full citation) have been com­puted using periodic density functional theory employing the same methods as described in Section 3[link]. As expected, metastable form A possesses a higher energy than stable form B. The energy difference is about 23.6 kcal mol−1, which is quite significant, but well in line with the trend reported in Nyman & Day (2015View full citation). It should be noted that this difference could be a little lower, since the structure of form A has been determined at T = 293 K, while form B was determined at T = 120 K.

Supporting information


Computing details top

(2S)-2-[4-(2-methylpropyl)phenyl]propanoic acid; pyridine-3-carboxamide top
Crystal data top
C13H18O2·C6H6N2OV = 930.48 (7) Å3
Mr = 328.4Z = 2
Monoclinic, P21F(000) = 352
Hall symbol: P 2ybDx = 1.172 Mg m3
a = 27.9925 (13) ÅSynchrotron radiation
b = 5.5286 (2) ÅT = 293 K
c = 6.0213 (2) Åwhite
β = 93.112 (2)°
Data collection top
Synchrotron
diffractometer
2θmin = 2°, 2θmax = 40.284°, 2θstep = 0.004°
Radiation source: synchrotron, synchrotron
Refinement top
Rp = 0.0690 restraints
Rwp = 0.0851 constraint
Rexp = 0.007H-atom parameters constrained
R(F) = 0.293Weighting scheme based on measured s.u.'s
9572 data points(Δ/σ)max = 0.021
Profile function: Pseudo-VoigtBackground function: 8 Legendre polynoms
19 parametersPreferred orientation correction: none
Fractional atomic coordinates and isotropic or equivalent isotropic displacement parameters (Å2) top
xyzUiso*/Ueq
C10.105420.787810.961050.038*
C20.07150.896430.817980.038*
C30.063480.797850.607330.038*
C40.090170.60050.54540.038*
C50.12330.504580.699610.038*
C60.042231.109960.877760.038*
N10.130850.593780.904930.038*
N20.005131.097040.837830.038*
O10.061981.288750.962450.038*
H10.112210.851181.128280.045*
H20.035490.876840.496510.045*
H30.085050.520050.382150.045*
H40.144730.351150.657260.045*
H50.02120.940520.796670.057*
H60.025421.239610.881130.057*
O20.168990.953360.251320.038*
O30.214380.656220.139750.038*
H70.461131.102450.522830.057*
H80.436130.730290.158120.057*
H90.462171.343070.343010.057*
H100.492691.083510.244560.057*
H110.256271.116090.726620.057*
H120.394471.187080.324840.045*
H130.445781.120080.071070.057*
H140.204111.037690.67670.057*
H150.417061.352890.012170.057*
H160.379150.800010.223820.057*
H170.410410.77910.428490.057*
H180.363211.39270.207280.045*
H190.250540.752640.513990.045*
H200.380070.729430.075520.045*
H210.310250.632890.242430.045*
H220.282751.298560.367140.045*
H230.19040.652010.026520.057*
H240.216281.240690.532780.057*
C70.205740.835430.272050.038*
C80.244680.892420.443650.038*
C90.228631.085030.60750.038*
C100.289170.963050.32530.038*
C110.303151.202640.316050.038*
C120.344351.268280.213810.038*
C130.37311.099360.119130.038*
C140.359350.861310.131810.038*
C150.318430.795210.232640.038*
C160.416671.17680.001510.038*
C170.416891.089380.240170.038*
C180.410520.81920.255050.038*
C190.462321.163740.349060.038*
Atomic displacement parameters (Å2) top
U11U22U33U12U13U23
???????
Geometric parameters (Å, º) top
C1—C21.3844H13—C160.9567
C1—N11.3412H14—C90.8629
C1—H11.0736H15—C160.976
C2—C31.388H16—C181.0346
C2—C61.4926H17—C181.0678
C3—C41.3853H18—C120.9973
C3—H21.0926H19—C80.8924
C4—C51.3824H20—C141.0024
C4—H31.0817H21—C150.9291
C5—N11.3375H22—C110.8494
C5—H41.0778H24—C91.0229
C6—N21.3367C7—C81.4947
C6—O11.23C8—C91.5358
N2—H51.0005C8—C101.5195
N2—H61.0141C10—C111.3836
O2—C71.2189C10—C151.3762
O3—H230.931C11—C121.3852
O3—C71.3025C12—C131.3768
H7—C191.0989C13—C141.3749
H8—C181.0254C13—C161.5152
H9—C190.9924C14—C151.375
H10—C191.1216C16—C171.5169
H11—C91.0407C17—C181.5068
H12—C170.9546C17—C191.5192
C2—C1—N1123.09C11—C10—C15116.98
C2—C1—H1121.53H22—C11—C10112.72
N1—C1—H1115.35H22—C11—C12126.16
C1—C2—C3118.04C10—C11—C12120.87
C1—C2—C6124.07H18—C12—C11119.07
C3—C2—C6117.87H18—C12—C13119.18
C2—C3—C4119.52C11—C12—C13121.74
C2—C3—H2118.46C12—C13—C14117.06
C4—C3—H2122.01C12—C13—C16120.72
C3—C4—C5118.26C14—C13—C16122.17
C3—C4—H3121.47H20—C14—C13120.59
C5—C4—H3120.26H20—C14—C15117.86
C4—C5—N1123.12C13—C14—C15121.44
C4—C5—H4120.33H21—C15—C10118.14
N1—C5—H4116.56H21—C15—C14119.94
C2—C6—N2117.82C10—C15—C14121.9
C2—C6—O1119.71H13—C16—H15110.18
N2—C6—O1122.47H13—C16—C13112.11
C1—N1—C5117.95H13—C16—C17106.34
C6—N2—H5121.28H15—C16—C13109.01
C6—N2—H6118.44H15—C16—C17104.77
H5—N2—H6119.05C13—C16—C17114.15
H23—O3—C7108.46H12—C17—C16106.83
O2—C7—O3121.88H12—C17—C18117.2
O2—C7—C8122.55H12—C17—C1999.09
O3—C7—C8115.51C16—C17—C18111.46
H19—C8—C7104.72C16—C17—C19111.68
H19—C8—C9110.39C18—C17—C19109.99
H19—C8—C10108.12H8—C18—H16102.77
C7—C8—C9111.27H8—C18—H17114.98
C7—C8—C10108.4H8—C18—C17111.39
C9—C8—C10113.52H16—C18—H17113.18
H11—C9—H14107.37H16—C18—C17109.66
H11—C9—H24112.67H17—C18—C17104.96
H11—C9—C8108.96H7—C19—H9110.07
H14—C9—H24102.06H7—C19—H10113.29
H14—C9—C8111.6H7—C19—C17110.17
H24—C9—C8113.87H9—C19—H10112.29
C8—C10—C11120.4H9—C19—C17104.44
C8—C10—C15122.56H10—C19—C17106.12
Hydrogen-bond geometry (Å, º) top
D—H···AD—HH···AD···AD—H···A
C1—H1···O2i1.071.812.5933126.42
N2—H5···O1ii1.002.072.6637115.75
N2—H6···C1iii1.012.493.2916135.09
N2—H6···C2iii1.012.443.3031142.51
C18—H16···C141.032.072.8104126.30
C12—H18···C14iv1.002.373.3464166.97
C12—H18···C15iv1.002.193.0067138.27
C14—H20···C181.002.262.8104112.99
C15—H21···C11v0.932.433.3463168.62
C15—H21···C12v0.932.243.0067139.17
C11—H22···C90.852.452.8733111.23
O3—H23···N1vi0.931.812.6890155.66
Symmetry codes: (i) x, y, z+1; (ii) x, y1/2, z+2; (iii) x, y+1/2, z+2; (iv) x, y+1, z; (v) x, y1, z; (vi) x, y, z1.
Crystallographic data, profile and structural parameters for S-ibuprofen–nicotinamide cocrystal form A obtained after Rietveld refinement top
Crystal data
Chemical formulaC19H24N2O3
Mr328.4
Cell setting, space groupMonoclinic, P21
Temperature (K)293
a,b,c (Å)27.9925 (13), 5.5286 (2), 6.0213 (2)
β (°)93.112 (2
V3)930.48 (7)
Z2
F(000)352
µ (mm-1)0.079
Specimen shape, size (mm)Cylinder, 0.5
2θ range (°)1.2–40.284°
Data collection
BeamlineCRISTAL (SOLEIL)
Specimen mounting0.5 mm diameter Lindemann capillary
Data collection modeTransmission
Scan methodContinuous scan
Radiation typeSynchrotron 17 KeV, λ = 0.728302 Å
Binning size (°2θ)0.004
Refinement
R factors and goodness-of-fitR = 0.694, Rwp(nb) = 0.0847, Rexp = 0.007
Lattice cell parameters comparisons between S-ibuprofen–nicotinamide cocrystal form A and S-ibuprofen–nicotinamide cocrystal form B (Berry et al., 2008) top
Structurea (Å)b (Å)c (Å)β (°)V3)SymmetryReference
S-IBP:N Form A27.99255.52866.021393.112930.48P21This work
S-IBP:N Form B5.411055.88311.900690.0043598.5P21Berry et al. (2008)
 

Acknowledgements

This project has received funding from the Inter­reg 2 Seas programme 2014–2020 co-funded by the European Regional Development Fund.

Funding information

Funding for this research was provided by: European Regional Development Fund (subsidy contract No. 2S01-059_IMODE).

References

Return to citationAlshahateet, S. F. (2010). Mol. Cryst. Liq. Cryst. 533, 152–161.  CSD CrossRef Google Scholar
Return to citationBerry, D. J., Seaton, C. C., Clegg, W., Harrington, R. W., Coles, S. J., Horton, P. N., Hursthouse, M. B., Storey, R., Jones, W., Friščić, T. & Blagden, N. (2008). Cryst. Growth Des. 8, 1697–1712.  Web of Science CSD CrossRef CAS Google Scholar
Return to citationBoultif, A. & Louër, D. (2004). J. Appl. Cryst. 37, 724–731.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationBrittain, H. G. (2013). J. Pharm. Sci. 102, 311–317.  Web of Science CrossRef CAS PubMed Google Scholar
Return to citationCaglioti, G., Paoletti, A. & Ricci, F. P. (1958). Nucl. Instrum. 3, 223–228.  CrossRef CAS Web of Science Google Scholar
Return to citationChilds, S. L., Wood, P. A., Rodríguez-Hornedo, N., Reddy, L. S. & Hardcastle, K. I. (2009). Cryst. Growth Des. 9, 1869–1888.  Web of Science CSD CrossRef CAS Google Scholar
Return to citationCoelho, A. A. (2000). J. Appl. Cryst. 33, 899–908.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationCruz-Cabeza, A. J., Spackman, P. R. & Hall, A. V. (2024). Commun. Chem. 7, 284.  PubMed Google Scholar
Return to citationDavid, W. I. F., Shankland, K., van de Streek, J., Pidcock, E., Motherwell, W. D. S. & Cole, J. C. (2006). J. Appl. Cryst. 39, 910–915.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationDesiraju, G. R. (2010). J. Chem. Sci. 122, 667–675.  Web of Science CrossRef CAS Google Scholar
Return to citationDwivedi, S. K., Sattari, S., Jamali, F. & Mitchell, A. G. (1992). Int. J. Pharm. 87, 95–104.  CrossRef Web of Science Google Scholar
Return to citationFellah, N., Zhang, C. J., Chen, C., Hu, C. T., Kahr, B., Ward, M. D. & Shtukenberg, A. G. (2021). Cryst. Growth Des. 21, 4713–4724.  CSD CrossRef Google Scholar
Return to citationFleischman, S. G., Kuduva, S. S., McMahon, J. A., Moulton, B., Bailey Walsh, R. D., Rodríguez-Hornedo, N. & Zaworotko, M. J. (2003). Cryst. Growth Des. 3, 909–919.  Web of Science CSD CrossRef CAS Google Scholar
Return to citationFreer, A. A., Bunyan, J. M., Shankland, N. & Sheen, D. B. (1993). Acta Cryst. C49, 1378–1380.  CSD CrossRef CAS Web of Science IUCr Journals Google Scholar
Return to citationFriščić, T. & Jones, W. (2010). J. Pharm. Pharmacol. 62, 1547–1559.  Web of Science PubMed Google Scholar
Return to citationGarzón, L. C. & Martínez, F. (2004). J. Solution Chem. 33, 1379–1395.  Google Scholar
Return to citationGiannozzi, P., Andreussi, O., Brumme, T., Bunau, O., Buongiorno Nardelli, M., Calandra, M., Car, R., Cavazzoni, C., Ceresoli, D., Cococcioni, M., Colonna, N., Carnimeo, I., Dal Corso, A., de Gironcoli, S., Delugas, P., DiStasio, R. A., Ferretti, A., Floris, A., Fratesi, G., Fugallo, G., Gebauer, R., Gerstmann, U., Giustino, F., Gorni, T., Jia, J., Kawamura, M., Ko, H.-Y., Kokalj, A., Küçükbenli, E., Lazzeri, M., Marsili, M., Marzari, N., Mauri, F., Nguyen, N. L., Nguyen, H.-V., Otero-de-la-Roza, A., Paulatto, L., Poncé, S., Rocca, D., Sabatini, R., Santra, B., Schlipf, M., Seitsonen, A. P., Smogunov, A., Timrov, I., Thonhauser, T., Umari, P., Vast, N., Wu, X. & Baroni, S. (2017). J. Phys. Condens. Matter 29, 465901.  CrossRef PubMed Google Scholar
Return to citationGiannozzi, P., Baroni, S., Bonini, N., Calandra, M., Car, R., Cavazzoni, C., Ceresoli, D., Chiarotti, G. L., Cococcioni, M., Dabo, I., Dal Corso, A., de Gironcoli, S., Fabris, S., Fratesi, G., Gebauer, R., Gerstmann, U., Gougoussis, C., Kokalj, A., Lazzeri, M., Martin-Samos, L., Marzari, N., Mauri, F., Mazzarello, R., Paolini, S., Pasquarello, A., Paulatto, L., Sbraccia, C., Scandolo, S., Sclauzero, G., Seitsonen, A. P., Smogunov, A., Umari, P. & Wentzcovitch, R. M. (2009). J. Phys. Condens. Matter 21, 395502.  Web of Science CrossRef PubMed Google Scholar
Return to citationGražulis, S., Chateigner, D., Downs, R. T., Yokochi, A. F. T., Quirós, M., Lutterotti, L., Manakova, E., Butkus, J., Moeck, P. & Le Bail, A. (2009). J. Appl. Cryst. 42, 726–729.  Web of Science CrossRef IUCr Journals Google Scholar
Return to citationGrimme, S., Antony, J., Ehrlich, S. & Krieg, H. (2010). J. Chem. Phys. 132, 154104.  Web of Science CrossRef PubMed Google Scholar
Return to citationGroom, C. R., Bruno, I. J., Lightfoot, M. P. & Ward, S. C. (2016). Acta Cryst. B72, 171–179.  Web of Science CrossRef IUCr Journals Google Scholar
Return to citationGuerain, M., Derollez, P., Roca-Paixão, L., Dejoie, C., Correia, N. T. & Affouard, F. (2020). Acta Cryst. C76, 225–230.  Web of Science CSD CrossRef IUCr Journals Google Scholar
Return to citationGuerain, M., Guinet, Y., Correia, N. T., Paccou, L., Danède, F. & Hédoux, A. (2020). Int. J. Pharm. 584, 119454.  Web of Science CrossRef PubMed Google Scholar
Return to citationHigashi, K., Ueda, K. & Moribe, K. (2017). Adv. Drug Deliv. Rev. 117, 71–85.  Web of Science CrossRef CAS PubMed Google Scholar
Return to citationKabekkodu, S. N., Dosen, A. & Blanton, T. N. (2024). Powder Diffr. 39, 47–59.  Web of Science CrossRef CAS Google Scholar
Return to citationKarimi-Jafari, M., Padrela, L., Walker, G. M. & Croker, D. M. (2018). Cryst. Growth Des. 18, 6370–6387.  CAS Google Scholar
Return to citationLe Bail, A. (2004). Powder Diffr. 19, 249–254.  Web of Science CrossRef CAS Google Scholar
Return to citationLe Bail, A., Duroy, H. & Fourquet, J. L. (1988). Mater. Res. Bull. 23, 447–452.  CrossRef ICSD CAS Web of Science Google Scholar
Return to citationLutterotti, L. (2010). Nucl. Instrum. Methods Phys. Res. B 268, 334–340.  CrossRef Google Scholar
Return to citationMarkvardsen, A. J., Shankland, K., David, W. I. F., Johnston, J. C., Ibberson, R. M., Tucker, M., Nowell, H. & Griffin, T. (2008). J. Appl. Cryst. 41, 1177–1181.  Web of Science CrossRef CAS IUCr Journals Google Scholar
Return to citationMiwa, Y., Mizuno, T., Tsuchida, K., Taga, T. & Iwata, Y. (1999). Acta Cryst. B55, 78–84.  Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
Return to citationNyman, J. & Day, G. M. (2015). CrystEngComm 17, 5154–5165.  CrossRef Google Scholar
Return to citationOswald, I. D. H., Allan, D. R., McGregor, P. A., Motherwell, W. D. S., Parsons, S. & Pulham, C. R. (2002). Acta Cryst. B58, 1057–1066.  Web of Science CSD CrossRef CAS IUCr Journals Google Scholar
Return to citationPerdew, J. P., Burke, K. & Ernzerhof, M. (1996). Phys. Rev. Lett. 77, 3865–3868.  CrossRef PubMed CAS Web of Science Google Scholar
Return to citationPetrícek, V., Dusek, M. & Palatinus, L. (2014). Z. Kristallogr. Cryst. Mater. 229, 345–352.  Google Scholar
Return to citationRoisnel, T. & Rodríquez-Carvajal, J. (2001). Mater. Sci. Forum 378–381, 118–123.  CrossRef Google Scholar
Return to citationSchultheiss, N. & Newman, A. (2009). Cryst. Growth Des. 9, 2950–2967.  Web of Science CrossRef PubMed CAS Google Scholar
Return to citationSmith, G. S. & Snyder, R. L. (1979). J. Appl. Cryst. 12, 60–65.  CrossRef CAS IUCr Journals Web of Science Google Scholar
Return to citationter Horst, J. H., Deij, M. A. & Cains, P. W. (2009). Cryst. Growth Des. 9, 1531–1537.  Web of Science CrossRef CAS Google Scholar
Return to citationVishweshwar, P., McMahon, J. A., Bis, J. A. & Zaworotko, M. J. (2006). J. Pharm. Sci. 95, 499–516.  Web of Science CrossRef PubMed CAS Google Scholar
Return to citationWHO (2007). Model lists of essential medicines. Geneva: WHO. https://www.who.int/Google Scholar
Return to citationWolff, P. M. de (1968). J. Appl. Cryst. 1, 108–113.  CrossRef IUCr Journals Web of Science Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoSTRUCTURAL
CHEMISTRY
ISSN: 2053-2296
Follow Acta Cryst. C
Sign up for e-alerts
Follow Acta Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds