A study into the effect of subtle structural details and disorder on the terahertz spectrum of crystalline benzoic acid

The phonon modes of crystalline benzoic acid have been investigated using terahertz time-domain spectroscopy, rigid molecule atom-atom model potential and plane-wave density functional theory lattice dynamics calculations. The simulation results show good agreement with the measured terahertz spectra and an assignment of the terahertz absorption features of benzoic acid is made with the help of both computational methods. Focussing on the strongest interactions in the crystal, we describe each vibration in terms of distortions of the benzoic acid hydrogen bonded dimers that are present in the crystal structure. The terahertz spectrum is also shown to be highly sensitive to the location of the carboxylic acid hydrogen atoms in the cyclic hydrogen-bonded dimers and we have systematically explored the influence of the observed disorder in the hydrogen atom positions on the lattice dynamics.     

Simulated Raman spectra of four tetraphenylbutadiene polymorphs

We report a combined experimental and theoretical investigation on the Raman spectra of the polymorphs α, β, γ, and δ of 1,1,4,4‐tetraphenyl‐1,3‐butadiene (TPB), in the region of the intramolecular modes. The interpretation of the polarized spectra is supported by ab‐initio calculations for the isolated molecules and by lattice dynamics calculations for the crystals. The calculations reproduce the peculiar, and surprisingly large, differences among the spectra of the various polymorphs. The phenyl groups of 1,1,4,4‐tetraphenyl‐1,3‐butadiene may arrange themselves around the butadiene skeleton in 2 stable conformers, which have either inversion (C_i) or 2‐fold (C₂) symmetry and therefore exhibit intramolecular vibrations with quite different Raman selection rules and spectra. The compound forms 4 crystalline polymorphs (α, β, γ, and δ) with different combinations of C_i and C₂ conformers, and correspondingly different intramolecular spectra. The theoretical calculations provide a quantitative analysis of the various spectra.     

Polarized IR spectra of the hydrogen bond in two different oxindole polymorphs with cyclic dimers in their lattices

This article focuses on the problem of remarkably strong changes in the fine structure patterns of the ν(N-H) and ν(N-D) bands attributed to the hydrogen and deuterium bonds accompanying the phase transition, which occurs between two polymorphic forms of oxindole. The lattices of these two different crystals contain hydrogen-bonded cyclic dimers differ in their geometry parameters. The source of these differences in the polymorph spectral properties results from the geometry relations concerning the dimers constituting the lattice structural units. In the case of the "alpha" phase, the hydrogen bond lengths of the dimers differ by 0.18 ?. This leads to the "off-resonance exciton coupling" weakly involving the dimer hydrogen bonds. For the "beta" phase, with practically symmetric dimers in the lattice, the spectra become typical for centrosymmetric hydrogen bond systems due to the full resonance of the proton or deuteron vibrations.     

Quantum chemical study of three polymorphs of the mononuclear spin-transition complex [Fe(DPPA)(NCS)2

The calculations of the high spin (HS) and low spin (LS) states of the [Fe(II)(DPPA)(NCS)(2)] complex have been performed at three experimentally observed geometries corresponding to three synthesized polymorphs with different spin-transition behavior. The structure optimization leads to a single molecular structure, suggesting that the existence of three geometries is not an intrinsic phenomenon but is induced by the crystal lattice. The structural difference between three forms can be reproduced by introducing the Madelung field of the crystal lattice. However, the calculations show that the differences in magnetic behavior of the three polymorphs cannot be attributed only to variations of the energy gap between two spin states.     

Solid-state density functional theory investigation of the terahertz spectra of the structural isomers 1,2-dicyanobenzene and 1,3-dicyanobenzene

The high-resolution waveguide terahertz (THz) time-domain spectra (20-100 cm(-1)) of the two structural isomers 1,2-dicyanobenzene (1,2-DCB) and 1,3-dicyanobenzene (1,3-DCB) have been modeled and assigned using solid-state density functional theory. The THz spectra of these similar molecules are distinctly different in the low-frequency region with the differences being driven by modifications of the crystal packing arrangement between the isomers. Simulations utilizing the hybrid density functionals B3LYP and PBE0 were performed to determine the origins of the observed vibrational features. External lattice vibrations (hindered translations and rotations) are found to dominate these spectra, reinforcing the need for proper solid-state models in the analysis of the THz spectra of organic molecular solids. These calculations were able to account for all of the observed spectral features exhibited by both isomers, even in the case of 1,2-DCB, where the spectrum was found to be the result of two coexisting crystalline polymorphs.     

Simultaneous Quantification of Three Polymorphic Forms of Carbamazepine using Raman Spectroscopy and Multivariate Calibration

Carbamazepine is a pharmaceutical product used to treat epilepsy and bipolar disorder. Some active pharmaceutical ingredients, such as carbamazepine, present polymorphism that may alter the bioavailability. Consequently, the determination of different polymorphic forms has become important for the pharmaceutical industry. In this work, polymorphic forms were synthesized and characterized by differential scanning calorimetry and X-ray diffraction. Raman spectroscopy was used to quantify mixtures of the three common polymorphic forms of carbamazepine. A ternary mixture design was used to create the calibration set of ten samples and six levels of concentration for each polymorph. Partial least squares was performed to build the prediction models. Ten spectra were obtained to obtain representative Raman spectra of the mixtures. The calibration models were built using the average spectra, and an external set of samples was used to evaluate the models. The partial least squares model gave a root mean square error of prediction of 6.2% for carbamazepine I, 6.8% for carbamazepine III, and 11.6% for carbamazepine dihydrate. The results showed that good results were obtained for the solid state characterization of the mixtures of polymorphs using a fast strategy for simultaneous analysis.     

DFT‐GIAO1H NMR chemical shifts prediction for the spectral assignment and conformational analysis of the anticholinergic drugs (−)‐scopolamine and (−)‐hyoscyamine

The relatively large chemical shift differences observed in the ¹H NMR spectra of the anticholinergic drugs (?)‐scopolamine 1 and (?)‐hyoscyamine 2 measured in CDCl₃ are explained using a combination of systematic/molecular mechanics force field (MMFF) conformational searches and gas‐phase density functional theory (DFT) single point calculations, geometry optimizations and chemical shift calculations within the gauge including/invariant atomic orbital (GIAO) approximation. These calculations show that both molecules prefer a compact conformation in which the phenyl ring of the tropic ester is positioned under the tropane bicycle, clearly suggesting that the chemical shift differences are produced by the anisotropic effect of the aromatic ring. As the calculations fairly well predict these experimental differences, diastereotopic NMR signal assignments for the two studied molecules are proposed. In addition, a cursory inspection of the published ¹H and ¹³C NMR spectra of different forms of 1 and 2 in solution reveals that most of them show these diastereotopic chemical shift differences, strongly suggesting a preference for the compact conformation quite independent of the organic or aqueous nature of the solvent. Copyright © 2010 John Wiley & Sons, Ltd.     

Raman spectra of molecular crystals. Carbon disulphide

Raman and far-infrared spectra of polycrytalline CS₂ samples at 79 and 18°K are reported, with assignments of molecular and lattice modes. Observed frequencies are compared with those from lattice dynamics calculations based on atom-atom interactions.     

The structure of water in p-sulfonatocalix[4]arene

Tetrasodium p-sulfonatocalix[4]arene exists as a hydrate with approximately 14 water molecules and has three polymorphic modifications, all of which contain a water molecule in the molecular cavity that is engaged in OH···π interactions. Single-crystal neutron structures are reported for two of these three forms and reveal a "compressed" water molecule with short OH bonds. Partial atomic charges and hardness analysis (PACHA) calculations based on the neutron coordinates give an OH···π interaction energy of 6.9-7.5 kJ mol(-1). The PACHA analysis also reveals the dominance of the charge-assisted hydrogen bonds from the Na(+)-coordinated water molecules. The instability of the crystal towards dehydration can be traced to an uncoordinated lattice water site. The remarkable calixarene-Na(+)-hydrate motif is conserved almost unchanged across all three polymorphs. A single-crystal neutron structure is also reported for pentasodium p-sulfonatocalix[4]arene·12H(2)O, which exhibits an intracavity water molecule that is engaged in both OH···π and OH···O hydrogen bonding. The shorter covalent bond to the hydrogen atom that forms the interaction with the aromatic ring is again apparent.     

13C NMR investigation of solid-state polymorphism in 10-deacetyl baccatin III

To investigate the origins of solid-state NMR shift differences in polymorphs, carbon NMR chemical shift tensors are measured for two forms of solid 10-deacetyl baccatin III: a dimethyl sulfoxide (DMSO) solvate and an unsolvated form. A comparison of ab initio computed tensors that includes and omits the DMSO molecules demonstrates that lattice interactions cannot fully account for the shift differences in the two forms. Instead, conformational differences in the cyclohexenyl, benzoyl, and acetyl moieties are postulated to create the differences observed. X-ray analysis of six baccatin III analogues supports the suggested changes in the cyclohexenyl and benzoyl systems. The close statistical match of the (13)C chemical shifts of both polymorphic forms with those calculated using the X-ray geometry of 10-deacetyl baccatin III supports the contention that the B, C, and D rings are fairly rigid. Therefore, the observed tensor differences appear to arise primarily from conformational variations in ring substituents and the cyclohexenyl ring.     

Infrared and ab initio studies of hydrogen bonding and proton transfer in the complexes formed by pyrazoles

The results of experimental studies and quantum mechanical calculations of vibrational spectra and structure of hydrogen bonded complexes formed by pyrazole (P) and 3,5-dimethylpyrazole (DMP) are presented. IR spectra of pyrazoles in solutions, gas phase, and solid state have been investigated in wide range of concentrations and temperatures. It has been found that in the gas phase both P and DMP reveal the equilibrium between monomers, dimers, and trimers. In solutions the equilibrium between monomers and trimers dominates, no bands, which can be attributed to dimers were detected. DMP retains the trimer structure in solid state, while in the case of pyrazole P, formation of the crystal provides another type of association. Geometrical and spectral characteristics of dimers and trimers, obtained by ab initio calculations, are presented and compared with experimental data.

IR spectra of solutions containing P and DMP with a number of acids (acetic and trifluoroacetic acids, pentachlorophenol, HBr) have been studied in parallel with ab initio calculations. It has been found that pentachlorophenol forms with pyrazoles complexes with one strong hydrogen bond O–HN, while NH pyrazole group remains unbonded. With carboxylic acids DMP forms 1:1 cyclic complexes with two hydrogen bonds. In the case of acetic acid, the complex in CH₂Cl₂ solution reveals molecular structure with OHN and C=OHN bonds, in accordance with results of the calculations. For trifluoroacetic acid, the calculations predict the molecular structure to be energetically more stable in the case of the isolated binary complex (in gas phase), while the experimental spectrum of CH₂Cl₂ solution gives an evidence of the proton transfer with formation of the cyclic ionic pair with two NH⁺O⁻ bonds. The agreement with experimental results can be improved by taking into account the influence of environment in the framework of Onsager or Tomasi models. The shape of proton potential function of the complexes and medium effect on its parameters, resulted from experimental data and calculations, are discussed. It has been found that the number of potential minima and their relative depth depend strongly on the method of calculations and the basic set. Under excess of trifluoroacetic acid, the formation of 2:1 acid–DMP complex has been detected. Spectral characteristics and results of calculations point to the cyclic structure of this complex, which includes homoconjugated bis-trifluoroacetate anion and DMPH⁺ cation. With HBr both studied pyrazoles were found to form ionic complexes including one or two pyrazole molecules per one acid molecule and correspondingly monocation or homoconjugated cation BHB⁺.     

Vibrational spectral analysis of DL-valine DL-valinium and DL-methionine DL-methioninium picrates

A comparative study of infrared and Raman spectra of DL-valine DL-valinium picrate (DL-VVP) and DL-methionine DL-methioninium picrate (DL-MMP) at the room temperature in 4000-50 cm(-1) range helps to determine the effect of hydrogen bonds in these crystals. The existence of the zwitterion and the protonated form in both the crystals has been observed. In DL-MMP, the methionine and the methioninium residues form a dimer through a strong hydrogen bond between them in the crystal. The characteristic bands due to gem-methyl group are observed in the case of DL-VVP. Factor group analysis has been made and the numbers of vibrational modes have been calculated. The tentative assignments of the observed bands are given. The picrate group forms the anion in both the crystals and it is unaffected by the presence of the cations.     

From crystal structure prediction to polymorph prediction: interpreting the crystal energy landscape

Many organic molecules are emerging as having many crystalline forms, including polymorphs and solvates, as more techniques are being used to generate and characterise the organic solid state. The fundamental scientific and industrial interest in controlling crystallisation is inspiring the development of computational methods of predicting which crystal structures are thermodynamically feasible. Sometimes, computing this crystal energy landscape will reveal that a molecule has one way of packing with itself that is sufficiently more favourable than any other so only this crystal structure will be observed. More frequently, there will be many energy minima that are energetically feasible, showing approximately equi-energetic compromises between the various intermolecular interactions allowed by the conformational flexibility. Such cases generally lead to multiple solid forms. At the moment, we usually calculate the lattice energy landscape, an approximation to the real crystal energy landscape at 0 K. Despite its limitations, many studies show that this is a valuable complement to solid form screening, which helps in discovering new structures as well as rationalising the solid forms that are found in experimental searches. The range of factors that can determine which of the thermodynamically feasible crystal structures are observed polymorphs, shows the many further challenges in developing crystal energy landscapes as a tool for control of the organic solid state.     

Crystalline polymorphism and molecular structure of sodium pravastatin

In this work different crystallization processes of sodium pravastatin are explored and a new polymorph is obtained. The analytical results of powder X-ray diffraction (PXRD) and thermal analysis for this new polymorph indicate that it is different from the polymorphs previously reported. This new crystal form shows different physical-chemical properties than the previous forms, such as crystallographic structure, thermal behavior, and melting point, 181.5 degrees C. Besides, all crystallization processes previously reported use an aprotic solvent as antisolvent. However, we propose a new crystallization process for sodium pravastatin that uses only protic solvents, overcoming industrial scaling and environmental problems. Variable-temperature PXRD experiments show a transformation between different crystal forms in the range of 80-120 degrees C. Solid-state 13C NMR, reported in this work for the first time, and Fourier transform infrared (FT-IR) studies of some polymorphs show some differences in intermolecular interactions, especially with carboxylate and hydroxyl groups. Quantum mechanical calculations of the pravastatin molecule are also presented for the first time, obtaining a molecular structure similar to the experimental structure existing within the crystal lattice of the tert-octylamonium salt of pravastatin.     

Calculation of the (29)si NMR chemical shifts of aqueous silicate species

A DFT methodology for calculating (29)Si NMR chemical shifts of silicate species typically present prior to nucleation in zeolite synthesis solutions, incorporating solvent effects through an implicit representation is presented. We demonstrate how our methodology can reproduce the experimentally observed spectra and, by comparison to well characterized peaks in two different experimental studies, demonstrate the transferability and robustness of the methodology. We discuss certain cases in which caution must be exercised when implicit solvent representations are used for calculating silicate cluster geometries: those cases in which intramolecular hydrogen bonding can play a significant role in the geometry. A number of reassignments of previous tentative experimental assignments are proposed, and we also make assignments for the challenging substituted four-ring species. We present all of our computed chemical shift for previously observed species together with a number of other viable silicate clusters to serve as a reference point for future experimental studies.     

Accurate lattice energies of organic molecular crystals from periodic turbomole calculations

Accurate lattice energies of organic crystals are important i.e. for the pharmaceutical industry. Periodic DFT calculations with atom‐centered Gaussian basis functions with the Turbomole program are used to calculate lattice energies for several non‐covalently bound organic molecular crystals. The accuracy and convergence of results with basis set size and k‐space sampling from periodic calculations is evaluated for the two reference molecules benzoic acid and naphthalene. For the X23 benchmark set of small molecular crystals accurate lattice energies are obtained using the PBE‐D3 functional. In particular for hydrogen‐bonded systems, a sufficiently large basis set is required. The calculated lattice energy differences between enantiopure and racemic crystal forms for a prototype set of chiral molecules are in good agreement with experimental results and allow the rationalization and computer‐aided design of chiral separation processes. © 2018 Wiley Periodicals, Inc.     

Neutron spectroscopic study of hydrogen bonding dynamics in L-serine

Inelastic incoherent neutron scattering (IINS) spectra were obtained at 10 K for normal and deuterated L-serine. The geometry of L-serine molecule was optimized for the zwitterion form using ab initio HF, MP2 and DFT (B3LYP) levels with 6-31G* and 6-311 + +G4** basis sets. The theoretical frequencies of normal and d4-L-serine were compared with IINS spectra. Normal coordinate analysis and band assignments based on ab initio calculations and experimental data were presented. IINS frequencies due to the out-of-plane gamma(N-H...O) hydrogen bond motions were observed and identified.     

Hydrogen desorption kinetics from the Si(1-x)Gex(100)-(2x1) surface

We study the influence of germanium atoms upon molecular hydrogen desorption energetics using density functional cluster calculations. A three-dimer cluster is used to model the Si((1-x))Ge(x)(100)-(2x1) surface. The relative stabilities of the various monohydride and clean surface configurations are computed. We also compute the energy barriers for desorption from silicon, germanium, and mixed dimers with various neighboring configurations of silicon and germanium atoms. Our results indicate that there are two desorption channels from mixed dimers, one with an energy barrier close to that for desorption from germanium dimers and one with an energy barrier close to that for desorption from silicon dimers. Coupled with the preferential formation of mixed dimers over silicon or germanium dimers on the surface, our results suggest that the low barrier mixed dimer channel plays an important role in hydrogen desorption from silicon-germanium surfaces. A simple kinetics model is used to show that reasonable thermal desorption spectra result from incorporating this channel into the mechanism for hydrogen desorption. Our results help to resolve the discrepancy between the surface germanium coverage found from thermal desorption spectra analysis, and the results of composition measurements using photoemission experiments. We also find from our cluster calculations that germanium dimers exert little influence upon the hydrogen desorption barriers of neighboring silicon or germanium dimers. However, a relatively larger effect upon the desorption barrier is observed in our calculations when germanium atoms are present in the second layer.     

Phonon dynamics and electron-phonon coupling in pristine picene

The paper reports a complete analysis of the phonon structure of crystalline picene, a recently announced organic semiconductor. Both lattice and intramolecular vibrations are investigated. An exhaustive assignment of lattice phonons is obtained through polarized Raman spectra assisted by lattice dynamics calculations based on a well tested atom-atom potential model. Raman, infrared spectra and density functional (DFT) calculations are used for the characterization of intramolecular modes. Coupling between low-frequency molecular vibrations and lattice phonons is accounted for. Molecule-to-molecule transfer integrals, as well as the Peierls and Holstein (non-local and local) coupling constants, are evaluated through the semiempirical method INDO/S (Intermediate Neglect of Differential Overlap with Spectroscopic parametrization).     

Study of NH stretching vibrations in small ammonia clusters by infrared spectroscopy in He droplets and ab initio calculations

Infrared spectra of the NH stretching vibrations of (NH3)n clusters (n = 2-4) have been obtained using the helium droplet isolation technique and first principles electronic structure anharmonic calculations. The measured spectra exhibit well-resolved bands, which have been assigned to the nu1, nu3, and 2nu4 modes of the ammonia fragments in the clusters. The formation of a hydrogen bond in ammonia dimers leads to an increase of the infrared intensity by about a factor of 4. In the larger clusters the infrared intensity per hydrogen bond is close to that found in dimers and approaches the value in the NH3 crystal. The intensity of the 2nu4 overtone band in the trimer and tetramer increases by a factor of 10 relative to that in the monomer and dimer, and is comparable to the intensity of the nu1 and nu3 fundamental bands in larger clusters. This indicates the onset of the strong anharmonic coupling of the 2nu4 and nu1 modes in larger clusters. The experimental assignments are compared to the ones obtained from first principles electronic structure anharmonic calculations for the dimer and trimer clusters. The anharmonic calculations were performed at the M?ller-Plesset (MP2) level of electronic structure theory and were based on a second-order perturbative evaluation of rovibrational parameters and their effects on the vibrational spectra and average structures. In general, there is excellent (<20 cm(-1)) agreement between the experimentally measured band origins for the N-H stretching frequencies and the calculated anharmonic vibrational frequencies. However, the calculations were found to overestimate the infrared intensities in clusters by about a factor of 4.