the existence of different crystalline forms ( polymorphs , hydrates , and solvates ) represents one of the most challenging phenomena in solid - state chemistry and related sciences , since we are still not able to predict the number of practically relevant forms and the conditions under which these can be grown or exist .
the existence of different solid - state forms of a compound is important as these usually show different physical properties , for example , solubility , density , hardness , melting point , etc .
this is true for pharmaceuticals ( the majority of the active ingredients are used in a crystalline form(2 ) ) , because the solid - state form can profoundly influence the manufacturing process , the long - term stability , and the performance of drug products , and for many other materials used in the chemical industry ( plant protection substances , dyes , explosives , etc . ) .
the present study deals with the solid - state of -resorcylic acid ( 2,4-dihydroxybenzoic acid , ra , figure 1 ) , a small organic molecule exhibiting molecular flexibility and the ability to form different hydrogen bonding motifs .
the compound is used as a starting material for the production of dyestuffs , pharmaceuticals , cosmetic preparations , and fine organic chemicals .
the cambridge structural database ( csd)(5 ) contains entries for three ra solid - state forms , namely two anhydrates ( zzzeeu:(6)p1 , z = 2 and zzzeeu01 to zzzeeu04:(7)p21/n , z = 1 , measured at 90 , 100 , 110 , and 150 k ) , and a hemihydrate ( qivtuk:(8)p1 , z = 1 ) . for the triclinic anhydrate ,
only the lattice dimensions have been reported , and the volume of zzzeeu corresponds to a monohydrate rather than an anhydrous form(9 ) but not to the new monohydrate described in this work . for the monoclinic polymorph ( form ii hereafter ) , the temperature range has been extended very recently down to 20 k.(10 ) furthermore , different hydrate stoichiometries , ranging from 0.5 to 3 mol water per mol of acid can be found in literature reports , but only the crystal structure of the hemihydrate has been determined . a recent study comparing six isomeric dihydroxybenzoic acids failed to crystallize new polymorphs by melt crystallization and sublimation experiments.(12 ) joint experimental and computational studies
have shown that there is no cooperative hydrogen atom disorder in the cooh and o - oh groups in form ii at temperatures up to 150 k. however , ra was neither subjected to a systematic solution crystallization screen nor to a comprehensive solid - state characterization program , and also theoretical predictions of possible crystal structures have not been reported so far .
global ( conf_p1 ) and second lowest conformational minima ( conf_p2 ) of -resorcylic acid ( ra ) .
the intramolecular degrees of freedom ( dihedral angles ) that were optimized within the crystal energy minimizations are indicated with arrows : 1 : c6c1c7o2 , 2 : c3c2o3h , 3 : c5c4o4h and 4 : c1c7o1h
. therefore , our investigation aimed at an efficient screening program , using an experimental and computational(16 ) approach to complement and validate the results and comprehensively characterize all ra solid - state forms at ambient conditions . the experimental screen was based on manual solution crystallizations of the compound in a variety of solvents and crystallization conditions , sublimation and moisture sorption experiments . the thermodynamic and kinetic stability of the solid - state forms were ascertained by hot - stage microscopy , differential scanning calorimetry , thermogravimetic analysis , and solvent - mediated transformation studies .
vibrational spectroscopy ( mid infrared and raman ) and x - ray diffractometry ( powder and single crystal ) were employed to determine the structural features of the phases . however , as neither high - throughput methodologies or other widely applied screening strategies(19 ) guarantee all possible forms will be found , we supported and complemented our manual screen with computational crystal structure prediction ( csp ) . by contrasting the thermodynamically feasible crystal structures with the experimentally observed ones , we discuss the factors that control crystallization and polymorphism of ra .
ra was purchased from fluka ( form ii ) . for the solvent screens , a set of 25 solvents
was chosen ( supporting information , section 1.1 ) , which were all of analytical quality .
crystallization conditions included solvent evaporation , fast and slow cooling crystallization , precipitation with a miscible antisolvent , vapor diffusion , and solvent - mediated transformation . in total ,
more than 150 manual crystallization experiments were performed ( conditions and crystallization outcomes are provided in the supporting information , tables s1s5 ) .
we have named the polymorphs according to the kofler notation using roman numerals in the order of the melting points ( i.e. , the highest melting is named form i ) and flagged the thermodynamically stable form at room temperature with the symbol . form ii was either prepared by slow crystallization from numerous solvents , including n - butanol , n - propanol , i - propanol , acetonitrile , ethyl methyl ketone , ethyl acetate , or by solvent - mediated transformation of any ra form , using water - free solvents that did not form a solvate .
form i could be obtained from solvent crystallization , but it predominantly grew concomitantly with form ii. the easiest way to produce form i was heating any ra form above the transition temperature of the polymorphic transition ii  i ( 150170 c ) .
however , decomposition , although slow compared to the polymorphic transformation , starts at ca .
other methods included sublimation experiments in the same temperature range or the desolvation of the hemihydrate ( hh ) , dimethyl formamide hemisolvate ( sdmf - i ) , dimethyl sulfoxide hemisolvate ( sdmso ) , or dioxane hemisolvate ( sdx ) at temperatures above 60 c .
the two hydrates could be prepared by crystallization from a hot , saturated water solution , with the resulting solid form depending on the cooling rate .
fast crystallization to the final temperature of 0 c ( in ice ) led to the monohydrate ( mh ) , whereas slow cooling ( test tube wrapped in aluminum foil ) produced the hemihydrate ( hh ) .
the dioxane hemisolvate ( sdx ) , dimethyl formamide 0.75-solvate ( sdmf - ii ) , and the dimethylsulfoxide hemisolvate ( sdmso ) were prepared from ii by solvent - mediated transformation experiments in the respective solvent , the acetic acid monosolvate ( saa ) by fast crystallization ( cooling a hot saturated solution in acetic acid to ca . 8 c ) .
finally , the dimethyl formamide hemisolvate ( sdmf - i ) was obtained as an intermediate desolvation product of the sdmf - ii solvate . every crystallization or solvent - assisted grinding experiment with pyridine resulted in the formation of the pyridinium salt.(22 ) for hot - stage thermomicroscopic ( htm ) investigations a reichert thermovar polarization microscope equipped with a kofler hot stage ( reichert , a ) was used .
photographs were taken with a digital camera ( olympus colorview iiiu digital camera , d ) .
dsc was performed with a dsc 7 ( perkin - elmer , norwalk , ct , usa ) using the pyris 2.0 software .
approximately 13  0.0005 mg sample ( um3 ultramicrobalance , mettler , ch ) was weighed into al - pans ( 25 l ) .
dry nitrogen was used as the purge gas ( purge : 20 ml min ) .
the instrument was calibrated for temperature with pure benzophenone ( mp 48.0 c ) and caffeine ( mp 236.2 c ) , and the energy calibration was performed with pure indium ( purity 99.999% , mp 156.6 c , heat of fusion 28.45 j g ) . tga was carried out with a tga7 system ( perkin - elmer , usa ) using the pyris 2.0 software .
two - point calibration of the temperature was performed with ferromagnetic materials ( alumel and ni , curie - point standards , perkin - elmer ) .
heating rates ranging from 10 to 20 k min were applied , and dry nitrogen was used as a purge gas ( sample purge : 20 ml min , balance purge : 40 ml min ) .
the stated error limits of thermochemical data are calculated as confidence intervals at a 95% level .
isothermal ( 25  0.1 c ) moisture sorption isotherms were acquired using a sps-11 moisture sorption analyzer ( projekt messtechnik , d ) .
the samples were gently ground prior to measurement to exclude the influence of particle size and surface area .
sorption and desorption cycles covered the 1090% rh range in 10% steps and the 010% range in 5% steps .
the equilibrium condition for each step was set to a mass constancy of 0.001% over 35 min .
spectra were recorded with a bruker ( bruker optic gmbh , d ) ifs 25 spectrometer connected to a bruker ir microscope i ( 15-cassegrain - objective , spectral range 4000 to 600 cm , resolution 4 cm , 64 scans per spectrum ) .
the samples ( rolled on a znse disk or fused between two znse windows ) were measured in transmission mode .
spectra were recorded with a bruker rfs 100 raman - spectrometer ( bruker analytische messtechnik gmbh , d ) , equipped with a nd : yag laser ( 1064 nm ) as the excitation source and a liquid - nitrogen - cooled , high sensitivity ge - detector . the spectra ( 128 scans per spectrum ) were recorded in aluminum sample holders with a laser power of 200 mw and a resolution of 2 cm .
experiments were performed on an oxford diffraction gemini r ultra ( 4-circle kappa - goniometer , 135 mm ruby ccd detector , mok radiation , monocapillary collimator ) with an oxford cryosystems 700 series cryostream plus low temperature attachment .
the single crystal structures of hh , sdmso , and the pyridinium salt were solved by direct methods using the program package wingx(23 ) ( sir2004(24 ) and shelxl97(25 ) ) .
all hydrogen atoms bonded to carbon atoms were generated by a riding model on idealized geometries with uiso(h ) = 1.2 ueq(c ) .
the polar hydrogens were identified from the difference map and refined isotropically , with the exception of h9 in sdmso , where the position was refined with a constrained oh bond distance . for further details ,
pxrd was used to determine the structure of form i. the sample was loaded in a rotating 1.0 mm borosilicate glass capillary and mounted on a bruker axs d8 powder x - ray diffractometer equipped with primary monochromator ( cuk1 , l = 1.54056  ) and lynxeye position sensitive detector .
data was collected at room temperature using a variable count time scheme ( supporting information , table s6 ) . the diffraction pattern indexed to a monoclinic unit cell ( omitting an impurity peak at 20.8 2 arising from a suspected decomposition product ) and the space group was determined to be p21/a based on a statistical assessment of systematic absences,(26 ) as implemented in the dash structure solution package.(27 ) the data were background subtracted and truncated to 50.5 2 for pawley fitting(28 ) ( pawley  = 15.91 ) .
simulated annealing was used to optimize the form i model against the diffraction data set ( 115 reflections ) in direct space .
the internal coordinate ( z - matrix ) description was derived from the hf/6 - 31g(d , p ) gas phase global conformational minimum ( conf_p1 ) , with oh distances normalized to 0.9  and ch distances to 0.95  .
the structure was solved using 200 simulated annealing runs of 2.5  10 moves per run as implemented in dash , allowing 7 degrees of freedom ( 6 external and 1 internal ) .
4.66 ( profile / pawley  ) and was used as the starting point for a rigid body rietveld refinement(29 ) in topas v4.1.(30 ) the rigid body description was derived from the z - matrix used in the simulated annealing runs and the final refinement included a total of 61 parameters ( 40 profile , 4 cell , 1 scale , 1 uiso , 9 preferred orientation , 3 position , and 3 rotation ) yielding a final rwp = 5.53 ( figure 2 ) .
powder x - ray diffraction pattern and rietveld fit ( rigid body ) of ra form i at 25 c .
we considered all eight planar hf/6 - 31g(d , p ) conformational minima produced by varying the torsion angles 1  3 ( figure 1 ) in the csp searches for anhydrate crystal structures , using a three - stage methodology.(31 ) first , z = 1 and z = 2 crystal structures were generated using crystalpredictor(32 ) in 25 common space groups for organic molecules .
the molecules were held rigid and the lattice energy was evaluated by an exp-6 potential with atomic charges derived using the chelpg scheme(33 ) and minimized .
all crystallographically distinct low energy crystal structures were used as starting points for optimizing the intermolecular lattice energy ( uinter ) , with an improved model for the intermolecular forces .
this was calculated using the fitexp-6 potential parameters and the distributed multipoles(37 ) derived from the pbe0/aug - cc - pvtz charge density using gdma2.(38 ) finally the 28 most stable structures were refined with crystaloptimizer(39 ) to allow small changes in conformation ( torsion angles 1  4 in figure 1 ) by minimizing the lattice energy , elatt = uinter + eintra , where eintra is the conformational energy penalty ( with respect to the global conformational minimum ) paid to improve the intermolecular interactions .
the conformational energy penalty was computed at the pbe0/6 - 31g(d , p ) level .
all isolated - molecule wave function calculations were performed using gaussian03(40 ) and intermolecular lattice energies by dmacrys.(41 ) more details of the conformational analysis , model testing , and search procedure and results are given in supporting information .
the solvent screening program resulted in a new ra polymorph ( form i ) , a new monohydrate ( mh ) , and five new solvates , namely acetic acid ( saa ) , dioxane ( sdx ) , dmso ( sdmso ) , and two dmf solvates ( sdmf - i and sdmf - ii ) .
the previously described anhydrous form ( form ii)(7 ) and the hemihydrate ( hh)(8 ) could be reproduced .
the relationships between the ra solid state forms ( figure 3 ) and their structures are described in experimental section 3.1 , and their relationship to the computed possible structures are described in section 3.2 .
transition pathways of the ra forms at different temperatures and relative humidities ( rh ) ; mh : monohydrate , hh : hemihydrate , i : form i , ii : form ii , saa : acetic acid monosolvate , sdmf - ii ( 0.75 ) and sdmf - i ( hemi ) : dimethyl formamide solvates , sdmso : dimethyl sulfoxide hemisolvate and sdx : dioxane hemisolvate , t : heating , * : removing from the mother liquor ,  mechanical  treatment , for example , grinding
form ii occurs in solvent crystallization / evaporation experiments as long needles , whereas form i occurs as a polycrystalline powder comprising multilayered plates .
htm and dsc experiments revealed that form ii transforms to form i upon heating at temperatures above 150 c .
optically this transformation can be recognized by a darkening of the form ii needles due to cleavage and the formation of small aggregates of form i ( figure 4 ) .
tga experiments and further thermomicroscopic investigations ( silicon oil embedding ) showed that no desolvation process is involved in this process , confirming that the reaction is a single component ( polymorphic ) transition .
the weak endotherm in the dsc curve ( figure 4 ) at an experimental transition temperature ( ttrs , exp ) of 159.8  0.5 c corresponds to the solidsolid transformation ii
 i with a transition enthalpy of trsh = 2.1  0.1 kj mol . from the fact that the transition is endothermic
, it can be concluded that the two polymorphs are enantiotropically related ( heat of transition rule ) .
this implies that the melting point of form ii is lower than that of form i , which melts at 231 c with decomposition ( second endotherm in dsc curve of figure 4 , overlapping with a broader endotherm , indicating the decomposition process ) .
form ii is thus the thermodynamically stable polymorph at room temperature ( rt ) and form i represents the high temperature form , which nevertheless shows a high kinetic stability at rt .
a reversible transformation i  ii was only observed if the sample was mechanically treated ( e.g. , grinding at room temperature ) .
the fact that form ii is thermodynamically stable at 20 c was also confirmed with solvent - mediated transformation experiments .
the energy - temperature diagram of the two polymorphs derived from the thermoanalytical data is provided in supporting information ( figure s2 ) .
hot - stage microscopy photographs of the ra form ii to form i transformation , along with the dsc curve that shows the transformation ( ii  i ) , melting ( mp ) , and decomposition ( dec . ) of form i ( heating rate : 10 k min ) . the hh forms transparent blocklike crystals ( figure 5 ) , which are stable when removed from the aqueous mother liquor . single crystals of hh turn opaque when heated above 70 c , while the original shape of the crystals is more or less maintained ( pseudomorphosis ) . in silicon oil preparations ,
the release of bubbles accompanies the darkening of the crystals and confirms the desolvation process ( see figure s1 , supporting information ) .
the mass loss of 5.55% ( tga , figure 6 ) is consistent with the theoretical water content of a hemihydrate ( 5.52% ) .
the second hydrate ( monohydrate , mh ) crystallizes as needles and was found to be very unstable at ambient conditions .
tga shows ( figure 6 ) that the dehydration of mh occurs below 60 c resulting in form ii. a transformation of mh to hh can only be observed in a water slurry . the measured mass loss ( figure 6 ) of ca .
9.5% corresponds to 0.9 mol water per acid molecule , which is slightly lower than expected for a monohydrate ( theory : 10.47% ) due to some loss of water during the preparation of the unstable mh .
photomicrographs ( dry preparation ) showing the dehydration process of ra hh to form i in a htm experiment ( heating rate ca .
all five solvates form transparent , thin , needle- or plate - like crystals losing their transparency within hours when removed from the mother liquor . upon heating ( figure 6 ) ,
desolvation of saa occurs in one step between about 40 and 60 c resulting in needle - shaped form ii crystals .
sdmf - ii shows two desolvation steps under heating . the first desolvation process ( figure 6 , mass loss ca .
50 and 70 c and corresponds to the transformation of the 0.75-solvate ( sdmf - ii ) to the hemisolvate ( sdmf - i ) , which then desolvates at ca .
90 c to form i. sdmso shows a peritectic melting / decomposition at 90 c . upon further holding the temperature ,
long plates ( form i ) slowly crystallize from the melt ( observed in htm , not visible in shown dsc curve ) .
it was not possible to clearly determine the stoichiometric solvent content from tga experiments , as the desolvation and decomposition process overlaps .
however , the course of the first step of the curve shows that the desolvation process decelerates before the fast decomposition process starts at about 150 c , indicating that the majority of the solvent is released . we can extrapolate the end of the first step to about 80% , suggesting that sdmso is a hemisolvate ( theoretical mass loss : 20.2% )
sdx shows an inhomogeneous melting process consisting of the following events : peritectic melting / decomposition , evaporation of the solvent ( dsc : endotherm at ca .
85 c ) , and the crystallization of form i ( exotherm above 110 c ) .
dsc ( black ) and tga ( red ) curves of ra hydrates and solvates ( pin - holed dsc capsules and a heating rate of 10 k min was used for all thermograms ) .
forms i , ii , and hh were subjected to a gravimetric moisture sorption study at 25 c ( figure 7 ) .
the anhydrous phases absorb water and transform only to hh even at the highest rh value ( 90% ) .
the moisture sorption cycle of form i was completed within 48 h , whereas only ca .
hh dehydrates only under rather dry conditions ( 5% rh and below ) demonstrating its high stability . under these conditions ,
the profile of the moisture sorptiondesorption isotherm ( sharp steps , strong hysteresis between sorption and desorption ) is a very clear indication that hh is a  stoichiometric 
hydrate.(44 ) from the huge hysteresis , we can deduce that the anhydrates and hh can be handled and stored as stable phases provided extremely low or high moisture conditions are avoided .
moisture sorption isotherm of ra performed at 25 c showing a moisture sorption cycle of anhydrous ra to hh and a desorption cycle of hh to a polymorphic mixture of forms i and ii .
the ft - ir spectra ( figure 8) allow an unambiguous identification of all forms , as variations were found in numerous regions , for example , (oh ) , (carh ) , (c = o ) , and (coh ) vibrations .
the most striking difference in the ir spectra of the two polymorphs concerns the (oh ) vibration of the p - oh group . in form
ii a sharp band occurs at 3373 cm , whereas in form i the band is broader and shifted to higher wavenumber ( 3432 cm ) implying weaker and more complex hydrogen bonding .
saa exhibits additional bands in the range of 970880 cm due to the out of plane oho hydrogen deformations(45 ) ( most intensive band : 920 cm ) .
the two dmf solvates showed two distinct vibrations at 1107 and 10941095 cm , whereas all other phases exhibit only one or three bands in this region .
on the basis of the (c = o ) vibrations , the dmf solvates can be discriminated , with sdmf - ii at 1666 cm and sdmf - i at 1652 cm .
the (s = o ) vibrations(46 ) in sdmso occur at 1002 cm ( strong band ) .
finally , the symmetric stretch vibrations of the ether group in sdx can be found at 865 cm .
ft - ir spectra of ra polymorphs ( i , ii ) , hydrates ( hh , mh ) , and solvates ( saa , sdmf - ii , sdmf - i , sdmso , and sdx ) .
ft - raman spectroscopy ( figure 9 ) also allows the identification of the different ra solid - state forms .
the presence of the solvent molecules ( acetic acid , dmf , dmso , and dioxane ) can be clearly seen in the spectral range around 2900 cm ( stretching vibrations of the aliphatic ch3 and ch2 groups , only present in the solvent molecules ) .
moreover , the lattice phonon vibrations below 150 cm differ markedly , indicating that the crystal structures of the ra solid - state forms are different and also that the solvates are not obviously isostructural as frequently observed if compounds form a series of such adducts .
ft - raman spectra of ra polymorphs ( i , ii ) , hydrates ( hh , mh ) , and solvates ( saa , sdmf - ii , sdmf - i , sdmso , and sdx ) .
the clear differences in pxrd patterns from each phase studied are consistent with the very distinct structures for all ra solid - state forms ( supporting information , figure s3 and table s6 ) .
temperature- and humidity - controlled powder x - ray diffraction studies of the dehydration process ( supporting information , figures s4 and s5 ) confirmed the results presented in sections 3.1.2 and 3.1.5 .
form i crystallizes in the monoclinic space group p21/a ( z = 1 ) and form ii in p21/n ( z = 1 ) , both adopting a conformation similar to conf_p1 ( figure 1 ) . the two polymorphs form inversion related r22(8 ) dimers and furthermore share the same 1d arrangement ( supramolecular construct(49 ) ) of inversion related dimeric units ( figure 10a , c )
however , they differ in the way the 1d constructs are linked to neighboring units , the o4ho3 hydrogen bonds , and consequently in the packing . in form ii , the o4ho3 interaction exhibits n glide symmetry ( figure 10b ) , leading to puckered hydrogen - bonded sheets ( figure 10a ) and to a dihedral angle between the 1d constructs ( benzene rings ) of 31. in form i , where 21 symmetry is present a dihedral angle between the benzene rings of 82 is observed .
hirshfeld dnorm surface(50 ) plots ( figure 10e , f ) clearly indicate that the p - oh proton forms only one h - bond in form ii , but that the assumed p - oh conformation in form i corresponds to two hydrogen bonds .
the p - oh proton can form a hydrogen bond to either the o- or p - oh oxygen , without otherwise changing the structure .
the difference in shape and larger volume of the surface around the p - oh in form i in contrast to form ii is consistent with possible disorder of the p - oh proton between the two positions , which due to the weak scattering contribution of hydrogen could not be verified based on the fit to the pxrd data alone but would be consistent with the broad ir band .
it appears that there is an unusual freedom in the position of the p - oh proton in the form i structure and so there will be ambiguities in the proton position .
crystal packing of form ii ( a , b ) and form i ( c , d ) viewed along the a - axis ( a ) ,   direction ( b ) , c - axis ( c ) , and b - axis ( d ) .
black dotted lines indicate the hydrogen bonds . for clarity only selected symmetry symbols are shown .
hirshfeld surfaces for forms ii ( e ) and i ( f ) were generated using the program crystal explorer v. 2.1,(52 ) with ch and oh distances normalized to neutron values , dnorm is mapped on the surfaces over the range 0.48 to 0.78  .
this function highlights contact distances relative to the sum of van der waals radii , with closest contracts shown in red . the circles in ( e , f ) highlight the p - oh group and hirshfeld surface around it ,
hh crystallizes in the triclinic space group p1 with two molecules of ra ( similar to conf_p1 and conf_p2 ) and one water molecule in the asymmetric unit , in agreement with the stoichiometry determined by tga .
the structure comprises columns of r22(8 ) dimers , formed by symmetrically equivalent , inversion related ra molecules ( homodimers ) . each dimeric unit
is directly linked to the adjacent alternate dimers with an o4ho4 hydrogen bond and indirectly via a water molecule bridge , that forms two hydrogen bonds ( o4howater and owaterho4 ) , leading to sheets .
in addition to connecting the acid dimers , water forms the only strong hydrogen bond between the parallel stacked sheets , leading to a 3d network structure ( figure 11 ) . the dehydration of hh involves a structural collapse , as despite the water interactions the o4ho4 hydrogen bonds also break due to the conformational change , which is around 180 in the torsion angle for the p - oh proton in every second molecule .
the units differ from the 1d construct present in the polymorphs in the location of the inversion center ( figure 11a ) .
crystal packing of hh : ( a ) ( 046 ) plane shown , ( b ) viewed along the a - axis showing the parallel sheets .
hydrogen bonds are denoted with dotted lines , and for clarity the water protons were omitted .
two acid molecules and one dmso molecule are present in the asymmetric unit of sdmso , which crystallizes in the orthorhombic space group p212121 .
the two symmetry inequivalent molecules adopt a conformation similar to conf_p1 and form r22(8 ) heterodimers .
strong hydrogen bonds are only present within the alternating heterodimer - dmso bands ( figure 12 ) and not in between the layers or stacking of the bands , which implies that the solvent could escape easily compared with hh and could explain why sdmso is not a long - living solvate .
the solvate and the desolvation product , form i , both contain similar 1d stacks of acids , along a in sdmso ( perpendicular distance : 3.024  , centroid distance : 6.398  ) and c in form i ( perpendicular distance : 3.489  , centroid distance : 5.198  ) .
this implies that on losing the solvent a rearrangement of the acid stacks with respect to adjacent stacks is necessary . on the other hand ,
crystal packing of sdmso : ( a ) corrugated heterodimer  dmso bands , view along the b - axis , ( b ) viewed along a - axis .
the pyridinium cation interacts with ra carboxylate via an ionic nho hydrogen bond ( supporting information , figure s5 ) .
the crystal energy landscape showed that only the two most stable conformations , differing in the position of the p - oh proton , could generate crystal structures within the likely energy range of polymorphism .
the lowest energy crystal structures had the lowest energy conformation ( conf_p1 , figure 1 ) , there were a few with the alternative p - oh conformation ( conf_p2 ) , but the majority of hypothetical structures contained both conformations , conf_p1 and conf_p2 , in a 1:1 ratio producing z = 2 structures .
the lack of any crystal structures varying the conformation of the intramolecular h - bond , involving the o - oh proton and c = o oxygen , is consistent with the experimental cocrystals , solvates , including hydrates and modifications.(7 ) experimentally , only conformations closely related to conf_p1 and conf_p2 , or a mix of the two , were observed , with the largest deviation of up to 10 in the position of the p - oh proton from planarity .
all low - energy structures ( figure 13 ) form one intra- and two intermolecular hydrogen bonds . the majority of computed crystal structures exhibit the experimentally observed r22(8 ) acid dimer ( d ) , but there are also acid catemer ( c ) , coohoh chain ( ch ) , and the combination of d and ch structures . the two experimental anhydrates were found to correspond to the lowest - energy structures within the two lowest energy groups of structures , with form i being the global minimum .
however , there is a low energy catemeric structure that is very competitive in energy with the known forms , and various other high density structures which also appear thermodynamically feasible .
lattice energy landscape for the ra anhydrates ( elatt = uinter + eintra ) after relaxation of the conformation within the crystal structure , classified by the hydrogen - bonding motif .
all structures with the same hydrogen bonding motif ( symbol ) within an ellipse are closely related .
the reproduction of form ii was excellent , given the neglect of thermal effects in our model , with an optimal root - mean square overlay of all non - hydrogen atoms in a 15 molecule coordination cluster ( rmsd15)(56 ) of 0.16  .
in contrast , the reproduction of form i was poor with an rmsd15 of ca . 1  , and it was wrongly predicted to be denser and more stable than form ii. the estimated helmholtz free energies(57 ) derived from the elastic constants(58 ) and k = 0 phonons(59 ) calculated in rigid - body harmonic approximation ( table s11 , supporting information ) for the two polymorphs bring forms i and ii to within 0.2 kj mol in energy .
hence , the static lattice energy model ( figure 13 ) is certainly overestimating the energy differences between the two polymorphs .
since proton disorder is suggested by the ir spectrum and the hirshfeld dnorm plot ( section 3.1.3 and 3.1.4 ) , we computationally generated three alternative ordered versions of form i differing only in the proton conformations which were kept fixed ( supporting information , section 2.3 ) .
this resulted in structures that were slightly less dense , and in similar or better agreement with the experimental structure than the global minimum on the crystal energy landscape .
hence , proton disorder in form i would result in a less dense structure , and some of the destabilization(60 ) may be counteracted by configurational entropy .
the calculations also showed that there was no barrier in this crystal structure to the proton moving by 30 from the lowest energy planar conformation .
unfortunately , we were not able to grow single crystals suitable for structure determination for all solid - state forms . from the ir and raman spectra , it could be assumed that in all phases the ra molecules are dimerized , as the symmetric (c = o ) band appeared in the region 16751625 cm .
furthermore , information about the hydrogen bonding of the p - oh group in the two anhydrates could be derived from the ir spectra .
the shift of (oh ) in form ii to lower wavenumbers indicates a stronger intermolecular oho hydrogen bond compared to form i. the broadening in the form i spectrum was attributed to the undetermined p - oh proton position between adjacent o- and p - oh groups ( see section 3.1.4 ) . on the basis of the experimental as well as theoretical results , we could not rule out the p - oh proton showing static or dynamic disorder .
static proton disorder has been shown to be possible by lattice energy minimization calculations which also show that there is no barrier to changing the proton position within the form i crystal structure , so it could also be dynamically disordered .
further insight into the disorder could be obtained from periodic ab initio calculations,(62 ) but the possibilities of experimental verification are very limited without neutron diffraction data .
the ra solid - state forms differ considerably in their relative stability . at ambient conditions only the two anhydrates and hh are stable .
mh and all five solvates survive only for a short time after harvesting from the mother liquor .
the formation of such unstable solvent adducts is common in organic molecules ; however , these are easily overlooked , especially when the product is dried prior to analysis or when the analysis is not performed immediately after the removal from the mother liquor.(44 ) in contrast to the metastable solvates , hh , in which water plays a structural role , is a very stable phase , as dehydration occurs only at temperatures higher than 60 c or under the driest conditions ( below 10% rh , 25 c ) .
the extreme hysteresis between the sorption and desorption process of hh attests to the high kinetic barrier of the hydration and dehydration processes with respect to vapor pressure .
hh ) exhibit a desolvation temperature above 60 c in the dsc / tga experiments and revert to the metastable form i. the monosolvates ( saa and mh ) desolvate below 60 c to the thermodynamically most stable form ii. therefore , it is obvious that the desolvation temperature is critical for the formation of a specific polymorph .
even though a large number of experiments were performed , we can not guarantee that all possible ra forms were found , as the range of variables in crystallization experiments that could influence the outcome is very large .
the anhydrate crystal energy landscape showed that the two known forms are the most likely polymorphs .
the alternative dimerbased structures ( figure 13 ) are sufficiently similar to the known polymorphs suggesting that they could readily transform to form i or ii in the unlikely event that they were distinct at the nucleation stage .
however , the lattice energy landscape suggests that other polymorphs could exist , which do not have the carboxylic acid dimer motif , but form exclusively cooh catemers or chains .
a csd(5 ) analysis of o - oh benzoic acid derivatives(63 ) showed that the r22(8 ) carboxylic acid dimer motif is the predominant hydrogen bonding motif , as seen in all four structurally characterized ra phases .
the relative stability of the catemeric structures on the crystal energy landscape may be overestimated by the computational model .
we can not exclude the possibility that the formation of the dimer motif is kinetically favored in all the crystallization conditions we have explored , hence ruling out the formation of catemers during our screening experiments .
-resorcylic acid represents another organic molecule exhibiting a complex solid - state behavior , that is , polymorphism , salt , and solvate ( hydrate ) formation .
the experimental search has resulted in seven new solid - state forms ( anhydrate , monohydrate , and five solvates ) , in addition to the previously known , and structurally characterized anhydrate ii and hemihydrate hh.(8 ) we found no evidence for a triclinic polymorph ( zzzeeu(6 ) ) or a hydrate showing a water / acid stoichiometry greater than 1:1.(64 ) before we can answer the question  how exhaustive should a polymorph screen be ?  , perhaps we must add a supplementary question , namely , 
how exhaustive does your knowledge of physical form diversity need to be ?  ; that is , do we need to confirm that the most thermodynamically favorable form at 25 c is already known or do we want to identify all accessible solid forms , perhaps to select a metastable form that has superior properties ? it is worth considering these alongside any limitations in terms of available material , stability , time , etc . when designing a comprehensive screening strategy .
however , recognizing that it is impractical to sample all possible nucleation and growth conditions for a given molecule , we can conclude that as a minimum , a polymorph screen should include sufficient diversity to at least find the thermodynamically most stable polymorph ( non solvated ) and hydrate form(s ) as well as those thermodynamically unstable forms that show sufficient kinetic stability at ambient conditions ( termed metastable forms ) to enable them to be isolated and identified .
furthermore , a polymorph screen should identify all crystallization products , including intermediates in the crystallization process , for example , solvates , amorphous form , etc . and characterize transformation pathways in order to find alternative , or perhaps only , routes to obtain a specific form .
therefore , a robust screening approach requires the combination of a variety of experimental approaches(67 ) including crystallization from solution ( evaporation , cooling , slurry conversions ) , sublimation , crystallization from the melt , thermal and moisture dependent studies and desolvation methods .
the statement by maria kuhnert - brandstaetter,(68 )  probably every substance is potentially polymorphous .
the only question is , whether it is possible to adjust the external conditions in such a way that polymorphism can be realized or not 
, reflects the problem that there are numerous additional strategies , at ambient and nonambient conditions , which might be adequate for the nucleation and growth of further solid - state forms .
it is practically not feasible to cover the whole range of techniques that have been shown to produce new polymorphs for certain systems,(21 ) particularly using all kinds of templates , additives , seeds , or impurities that might generate new forms or stabilize metastable phases . the instability of many solvates of -resorcylic acid , the possibility of considerable variation in the p - oh position in form i , and the computed thermodynamically feasible anhydrates mean that we can not exclude the discovery of other solid - state forms .
therefore , our results for -resorcylic acid emphasize the problems of determining the complete set of solid - state forms when there is evidence of disorder or short - lived intermediate crystallization products whose lifetime is very dependent on conditions .
thus , the consistency between the experimental and computational techniques used in this study adds confidence that the practically most important -resorcylic acid solid - state forms at ambient conditions have been characterized , showing the value of calculating the crystal energy landscape as part of the screening process .
 a combined experimental and computational study was undertaken to establish the solid - state forms of -resorcylic acid ( 2,4-dihydroxybenzoic acid ) . 
 the experimental search resulted in nine crystalline forms : two concomitantly crystallizing polymorphs , five novel solvates ( with acetic acid , dimethyl sulfoxide , 1,4-dioxane , and two with n , n - dimethyl formamide ) , in addition to the known hemihydrate and a new monohydrate . 
 form ii , the thermodynamically stable polymorph at room temperature , was found to be the dominant crystallization product . a new , enantiotropically related polymorph ( form i ) 
 was obtained by desolvation of certain solvates , sublimation experiments , and via a thermally induced solidsolid transformation of form ii above 150 c . to establish their structural features , interconversions , and relative stability , 
 all solid - state forms were characterized with thermal , spectroscopic , x - ray crystallographic methods , and moisture - sorption analysis . 
 the hemihydrate is very stable , while the five solvates and the monohydrate are rather unstable phases that occur as crystallization intermediates . 
 complementary computational work confirmed that the two experimentally observed -resorcylic acid forms i and ii are the most probable polymorphs and supported the experimental evidence for form i being disordered in the p - oh proton position . 
 these consistent outcomes suggest that the most practically important features of -resorcylic acid crystallization under ambient conditions have been established ; however , it appears impractical to guarantee that no additional metastable solid - state form could be found . 