Understanding the influence of polymorphism on phonon spectra: lattice dynamics calculations and terahertz spectroscopy of carbamazepine

Rigid molecule atomistic lattice dynamics calculations have been performed to predict the phonon spectra of the four polymorphs of carbamazepine, and these calculations predict that there should be differences in the spectra of all four forms. Terahertz spectra have been measured for forms I and III, and there are clearly different features between polymorphs' spectra, that are accentuated at low temperature. While carbamazepine adopts the same hydrogen bonded dimers in all of its known polymorphs, the calculations show that differences in packing arrangements of the dimers lead to changes in the frequency ranges for each type of hydrogen bond vibration, giving a physical explanation to the observed differences between the spectra. Although the agreement between calculated and observed spectra does not allow a definitive characterization of the spectra, it is possible to make tentative assignments of many of the observed features in the terahertz region for the simpler form III; we can only make some tentative assignments of specific modes in the more complex spectrum of form I. While harmonic rigid molecule lattice dynamics shows promise for understanding the differences in spectra between polymorphs of organic molecules, discrepancies between observed and calculated spectra suggest areas of improvement in the computational methods for more accurate modeling of the dynamics in molecular organic crystals.     

Extended tunnelling states in the benzoic acid crystal: infrared and Raman spectra of the OH and OD stretching modes

We compare Raman and infrared spectra of the nuOH/OD modes in benzoic acid crystal powders at 7 K. The extremely sharp Raman bands contrast to the broad infrared profiles and suggest adiabatic separation of hydrogen (deuterium) dynamics from the crystal lattice. There is no evidence of any proton-proton coupling term. The assignment scheme is consistent with a quasisymmetric double-minimum potential, largely temperature independent. Tunnel splitting is a major band shaping mechanism, in addition to anharmonic coupling with lattice modes. The proton/deuteron dynamics are rationalized with nonlocal pseudoparticles and extended states. We propose a symmetry-related damping mechanism to account for the broad infrared profiles, as opposed to the sharp Raman bands. We assign spectral features to distinct interconversion mechanisms based on either pseudoparticle transfer or adiabatic pairwise transfer. We establish close contacts with theoretical models based on first-principles calculations.     

Theoretical scrutinization of nine benzoic acid dimers: Stability and energy decomposition analysis

Aromatic carboxylic acids are able to form diverse dimers and multimers due to their hydrogen bond donor and acceptor cites, as well as the aromatic rings. In this work, we examine nine benzoic acid dimers stabilized by hydrogen bonding and stacking interactions. Interacting quantum atoms methodology revealed that dominant attractive interactions in all of them, including hydrogen bonded systems, are due to exchange-correlation. Coulomb interactions are significant only in the most stable dimer with a double hydrogen bond, although the corresponding energy term is almost two times lower compared to the nonclassical one. Since interacting quantum atoms approach treats monomers binding by considering electronic energy only, in order to examine dissociation kinetics we performed density functional theory-based molecular dynamics simulations of selected stacked dimers: in 40% of the studied systems at 300 K thermal energy was sufficient to overpower barrier for dissociation within 1 ps, which resulted in the separation of the monomers, whereas 20% of them remained in the stacked position even after 5 ps. These results highlight the importance of noncovalent interactions, particularly weak stacking interactions, on the structure and dynamics of carboxylic acids and their derivatives.     

Substituent Effects on the Low-Frequency Vibrational Modes of Benzoic Acid and Related Compounds

Well-resolved absorption spectra of benzoic acid and its derivatives with one hydrogen atom replaced by a substituent group CH3, OH, NH2 or NO2 were reported in the frequency region between 6 and 67 cm^-1 at room temperature with terahertz time-domain spectroscopy （THz-TDS）. These substances can be distinguished easily based on the terahertz absorption spectra. The measurements suggested that even minor changes in the molecular configuration and chemical composition lead to distinct differences in THz spectrum. Density functional theory （DFT） method was used to assist the analysis and assignment of the individual THz absorption spectra of benzoic acid and its methyl derivatives. Observed THz responses of samples can be assigned to the collective vibrations associated with intermolecular hydrogen bonds.     

Isotope effects associated with tunneling and double proton transfer in the hydrogen bonds of benzoic acid

The isotope effects associated with double proton transfer in the hydrogen bonds of benzoic acid (BA) dimers have been measured using field-cycling (1)H NMR relaxometry and quasielastic neutron scattering. By studying mixed isotope (hydrogen and deuterium) samples, the dynamics of three isotopologues, BA-HH, BA-HD, and BA-DD, have been investigated. Low temperature measurements provide accurate measurements of the incoherent tunneling rate, k(0). This parameter scales accurately with the mass number, m, according to the formula k(0)=(E/m)e(-Fm) providing conclusive evidence that the proton transfer process is a strongly correlated motion of two hydrons. Furthermore, we conclude that the tunneling pathway is the same for the three isotopologue species. Measurements at higher temperatures illuminate the through barrier processes that are mediated via intermediate or excited vibrational states. In parallel with the investigation of proton transfer dynamics, the theoretical and experimental aspects of studying spin-lattice relaxation in single crystals of mixed isotope samples are investigated in depth. Heteronuclear dipolar interactions between (1)H and (2)H isotopes contribute significantly to the overall proton spin-lattice relaxation and it is shown that these must be modeled correctly to obtain accurate values for the proton transfer rates. Since the sample used in the NMR measurements was a single crystal, full account of the orientation dependence of the spin-lattice relaxation with respect to the applied B field was incorporated into the data analysis.     

Vibrational spectra of anhydrous and monohydrated caffeine and theophylline molecules and crystals

Density functional theory and classical molecular dynamics simulations are used to investigate the vibrational spectra of caffeine and theophylline anhydrous and monohydrate molecules and those of their crystalline anhydrous and monohydrated states, with emphasis in the terahertz region of the spectra. To better understand the influence of water in the monohydrate crystal spectra, we analyze the vibrational spectra of water monomer, dimer, tetramer, and pentamer, and also those of liquid water at two different temperatures. In small water clusters, we observe the progressive addition of translational and librational modes to the terahertz region of the spectra. The water spectra predicted by rigid and flexible water models is examined with classical molecular dynamics, and the respective peaks, especially in the terahertz region, are compared with those found in the small clusters. Similar analysis done for caffeine and theophylline monohydrate molecules using density functional theory clearly shows the presence of water modes in the librational states and in the water stretching region. Molecular dynamics of caffeine and theophylline anhydrous and monohydrate crystals reveal the influence of vibrations from the molecule-molecule (caffeine or theophylline) crystal stacks and those from the water-molecule interactions found in the monohydrate molecules and new modes from molecule-molecule, water-molecule, and water-water hydrogen bonding interactions arising from collective effects in the crystal structure. Findings illustrate challenges of terahertz technology for the detection of specific substances in condensed phases.     

Theoretical modeling of the O-H stretching IR bands of hydrogen-bonded dimers of benzoic acid in S(0) and S(1) electronic states

Theoretical model for vibrational interactions in the hydrogen-bonded dimer of benzoic acid is presented. The model takes into account anharmonic-type couplings between the high-frequency O-H and the low-frequency O[cdots, three dots, centered]O stretching vibrations in two hydrogen bonds, resonance interactions (Davydov coupling) between two hydrogen bonds in the dimer, and Fermi resonance between the O-H stretching fundamental and the first overtone of the O-H in-plane bending vibrations. The vibrational Hamiltonians and selection rules for the C(2h) geometry in the S(0) state and for the C(s) in-plane bent geometry in the S(1) state of the dimer are derived. The model is used for theoretical simulation of the O-H stretching IR absorption bands of benzoic acid dimers in the gas phase in the electronic ground and first excited singlet states. Ab initio CIS and CIS(D)6-311++G(d,p) calculations have been performed to determine geometry, frequencies, and excited state energies of benzoic acid dimer in the S(1) state.     

The ground-state tunneling splitting of various carboxylic acid dimers

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.     

Liquid crystalline behaviour of benzoic acid derivatives containing alkoxyazobenzene

Liquid crystal trimers based on the hydrogen bonding dimerization of 4‐{n‐[4‐(4‐m‐alkoxy‐phenylazo)phenoxy]alkoxy}benzoic acid (BAm‐n) have been synthesized and characterized. Temperature‐dependent FTIR spectroscopic studies showed that the carboxylic acid groups in BAm‐n are associated to form H‐bonded cyclic dimers both in their crystalline and liquid crystalline phases. The trimers exhibited enantiotropic liquid crystalline behaviour except for BA1‐3 which showed monotropic behaviour, and the mesophases changed from nematic to smectic phase, with the increase of length of the spacer and the terminal substituents. Pronounced odd–even effects in the melting temperatures, clearing temperatures and nematic–isotropic enthalpy changes were observed.     

Ground and excited states of the monomer and dimer of certain carboxylic acids

The ground-state properties of the monomer and the dimer of formic acid, acetic acid, and benzoic acid have been investigated using Hartree-Fock (HF) and density functional theory (DFT) methods using the 6-311++G(d,p) basis set. Some of the low-lying excited states have been studied using the time-dependent density functional theory (TDDFT) with LDA and B3LYP functionals and also employing complete-active-space-self-consistent-field (CASSCF) and multireference configuration interaction (MRCI) methodologies. DFT calculations predict the ground-state geometries in quantitative agreement with the available experimental results. The computed binding energies for the three carboxylic acid dimers are also in accord with the known thermodynamic data. The TDDFT predicted wavelengths corresponding to the lowest energy n-pi* transition in formic acid (214 nm) and acetic acid (214 nm) and the pi-pi* transition in benzoic acid (255 nm) are comparable to the experimentally observed absorption maxima. In addition, TDDFT calculations predict qualitatively correctly the blue shift (4-5 nm) in the excitation energy for the pi-pi* transition in going from the monomer to the dimer of formic acid and acetic acid and the red shift (approximately 19 nm) in pi-pi* transition in going from benzoic acid monomer to dimer. This also indicates that the electronic interaction arising from the hydrogen bonds between the monomers is marginal in all three carboxylic acids investigated.     

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.     

Intermolecular potential and lattice dynamics of orthorhombic acetylene

The crystal structure and lattice dynamics of orthorhombic acetylene have been calculated with an intermolecular potential consisting of atom-atom and multipole-multipole interactions and including a hydrogen bond. A new assignment of the Raman lattice vibrations is discussed and utilized in the refinement of the potential parameters. The non-transferability of the potential to the cubic phase is attributed to the breaking of the hydrogen bond at the phase transition, and to the large anharmonicity expected for the high-temperature phase.     

Li+3 ion beam irradiation induced changes in the thermodynamic and electrical parameters of 4-n-(nonyloxy) benzoic acid

Li⁺³ ion beam irradiation studies on 4-n-(nonyloxy) benzoic acid (NOBA) have been carried out. Thermodynamic measurements demonstrate that all the phase transition temperatures are depressed for low value of irradiation fluence as compared to those for the pure NOBA and thereafter increase linearly with increase in irradiation fluence. It has also been observed that the temperature range of N and SmC phases in the heating and cooling cycles is improved due to irradiation. Dielectric investigations show that both the transverse as well as the longitudinal components of the dielectric permittivity are increased as compared with those of the pure sample due to irradiation. The observed change in the thermodynamic and dielectric parameters is attributed to the increased concentration of NOBA monomers over the NOBA dimers due to fragmentation of the hydrogen bonding in some of the NOBA dimers.     

Car-parrinello molecular dynamics simulation of liquid formic acid

First-principles molecular dynamics has been used to investigate the structural, vibrational, and energetic properties of formic acid, formic acid-formate anion dimers, and liquid formic acid in a periodically repeated box with 32 formic acid molecules. We found that in liquid formic acid the hydrogen-bonded clusters mainly consist of linear branching chains. From our simulation, we got good agreement with the available structural and dynamical data. We also studied the proton transfer in the cis-formic acid-formate anion dimer, and we showed that this proton transfer does not have any potential barrier. The hydrogen bonding statistics as well as the mean lifetime of the hydrogen bonds are analyzed.     

An NMR study of trimethyl acetic acid in a nematic phase

Proton NMR spectra of trimethyl acetic acid, dissolved in a benzoic acid type nematic liquid crystal, have been measured and analysed. Making use of the hydrogen bonding capabilities of this sort of nematic phase appears to be an excellent method to obtain very high solute orientations.     

X‐(Diiodoarsanyl)benzoic acid (X = 3 and 4)

The title compounds, 4‐(diiodoarsanyl)benzoic acid, (I), and 3‐(diiodoarsanyl)benzoic acid, (II), both [As(C₇H₅O₂)I₂], which possess a –COOH coordinating group, form molecular crystal structures composed of hydrogen‐bonded dimers, the packing differences of which are caused by the relative position of the diiodoarsanyl groups. The para isomer, with Z′ = 1, crystallizes in a layered structure with shortened contacts of the As atoms to only the arene rings of adjacent molecules. In contrast, the meta isomer, with Z′ = 3, forms separate rectangular blocks of three ribbons, each composed of dimeric molecular units positioned almost directly above each other and with the As atoms possessing only two As...I contacts to the I atoms of neighbouring molecules.     

腺嘌呤与富马酸共晶体的太赫兹光谱分析

利用太赫兹时域光谱(THz-TDS)技术在室温下对腺嘌呤、富马酸及两者的共晶体进行测量, 实验结果显示腺嘌呤与富马酸共晶体在0.92、1.24、1.52 THz处有明显的吸收峰, 与腺嘌呤和富马酸不同, 表明共晶体物相结构不同于原料. 根据腺嘌呤分子氢键供体与受体的结构特点, 使用密度泛函理论(DFT)对腺嘌呤与富马酸三种可能的共晶体结构进行模拟. 结果显示其中一种可能的共晶体结构在0.89、1.16、1.41 THz处存在特征吸收峰, 与实验结果较好吻合. 由此判断腺嘌呤与富马酸共晶体氢键形成位置为腺嘌呤的氨基与富马酸其中一个羧酸的碳氧双键形成氢键, 而此羧酸的羟基与腺嘌呤六元环上的邻位氮原子形成第二处氢键. 本文还结合理论模拟的结果对腺嘌呤与富马酸共晶体的特征吸收峰对应的相关振动模式进行了归属.     

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.     

Influence of lattice dynamics on charge transport in the dianthra[2,3-b:2',3'-f]-thieno[3,2-b]thiophene organic crystals from a theoretical study

The influence of lattice dynamics on carrier mobility has received much attention in organic crystalline semiconductors, because the molecular components are held together by weak interactions and the transfer integrals between neighboring molecular orbitals are extremely sensitive to small nuclear displacements. Recently, it has been shown that the dynamic disorder has little effect on hole mobility in the ab plane of pentacene, but a reasonable explanation is absent for such a puzzle. To better understand the effect of lattice vibrations on carrier transport, a further study is required for other organic materials. In this work, a mixed molecular dynamic and quantum-chemical methodology is used to assess the effect of nuclear dynamics on hole mobility in the dianthra[2,3-b:2',3'-f]-thieno[3,2-b]thiophene (DATT) crystals which exhibit high air stability with the hole mobility as large as that in rubrene-based devices. It is found that the lattice vibrations lead to an increasing encumbrance for hole transport in the ab plane of the DATT crystals as the temperature increases. By comparing the crystal structures of DATT and pentacene, the reduced hole mobility in DATT is attributed to the unsymmetric arrays of nearest-neighboring molecular dimers in the ab plane, because the electronic coupling exhibits unbalanced thermal fluctuations for the nearest-neighboring dimers which then induces a stronger oscillation for carriers along the directions with asymmetric packing. To further relate the dynamic disorder with hole transport, the variations of anisotropic mobilities are also analyzed. As a result, the negligible effect of lattice dynamics on the hole mobility in pentacene is explained by the centrosymmetric molecular packing of the nearest-neighboring molecular pairs in the ab plane.