Noncovalent interactions in paired DNA nucleobases investigated by terahertz spectroscopy and solid-state density functional theory

Cocrystallized adenine and thymine derivatives, along with the pure monomeric crystals, were investigated by terahertz spectroscopy and solid-state density functional theory (DFT). The methylated nucleobase derivatives crystallize in planar hydrogen-bonded adenine-thymine pairs similar to the manner found in DNA. The spectra obtained for 1-methylthymine, 9-methyladenine, and the 1:1 cocrystal in the range of 10-100 cm(-1) clearly demonstrate that absorptions in this spectral range originate from the uniquely ordered assembly and the intermolecular interactions found in each individual crystal system. The quality of spectral reproduction for the DFT simulations of each system was clearly improved by the inclusion of an empirical correction term for London-type dispersion forces to the calculations. Notably, it was found that these weak dispersion forces in the adenine-thymine cocrystal were necessary to produce a properly converged crystal structure and meaningful simulation of the terahertz vibrational spectrum.     

Noncovalent interactions between modified cytosine and guanine DNA base pair mimics investigated by terahertz spectroscopy and solid-state density functional theory

Modified cytosine and guanine nucleobases cocrystallize in a hydrogen bonding configuration similar to that observed in native DNA. The noncovalent interactions binding these base pairs in the crystalline solid were investigated using terahertz (THz) spectroscopy and solid-state density functional theory (DFT). While stronger hydrogen bonding interactions are responsible for the general molecular orientations in the crystalline state, it is the weaker dipole-dipole and dispersion forces that determine the overall packing arrangement. The inclusion of dispersion interactions in the DFT calculations was found to be necessary to accurately simulate the unit cell structure and THz vibrational spectrum. Using properly modeled intermolecular potentials, the lattice vibrational motions of the cytosine and guanine derivatives were calculated. The vibrational characters of the modes exhibited by the DNA base pair mimic in the THz region were primarily rotational motions and are indicative of the energies and the nature of vibrations that would likely be observed between similar base pairs in DNA molecules.     

Importance of accurate spectral simulations for the analysis of terahertz spectra: citric acid anhydrate and monohydrate

The terahertz (THz) spectra of crystalline solids are typically uniquely sensitive to the molecular packing configurations, allowing for the detection of polymorphs and hydrates by THz spectroscopic techniques. It is possible, however, that coincident absorptions may be observed between related crystal forms, in which case careful assessment of the lattice vibrations of each system must be performed. Presented here is a THz spectroscopic investigation of citric acid in its anhydrous and monohydrate phases. Remarkably similar features were observed in the THz spectra of both systems, requiring the accurate calculation of the low-frequency vibrational modes by solid-state density functional theory to determine the origins of these spectral features. The results of the simulations demonstrate the necessity of reliable and rigorous methods for THz vibrational modes to ensure the proper evaluation of the THz spectra of molecular solids.     

Molecular analysis of water clusters: Calculation of the cluster structures and vibrational spectrum using density functional theory

Density functional theory (DFT) is applied to obtain absorption spectra at THz frequencies for molecular clusters of H₂O. The vibrational modes of the clusters are calculated. Coupling among molecular vibrational modes explains their spectral features associated with THz excitation. THz excitation is associated with vibrational frequencies which are here calculated within the DFT approximation of electronic states. This is done for both isolated molecules and collections of molecules in a cluster. The principal result of the paper is that a crystal-like cluster of 38 water molecules together with a continuum solvent background is sufficient to replicate well the experimental vibrational frequencies.     

Theoretical analysis of the terahertz spectrum of the high explosive PETN.

The experimental solid-state terahertz (THz) spectrum (3 to 120 cm(-1)) of the high explosive pentaerythritol tetranitrate (PETN, C(5)H(6)N(4)O(12)) has been modeled using solid-state density functional theory (DFT) calculations. Solid-state DFT, employing the BP density functional, is in best qualitative agreement with the features in the previously reported THz spectrum. The crystal environment of PETN includes numerous intermolecular hydrogen-bonding interactions that contribute to large (up to 80 cm(-1)) calculated shifts in molecular normal-mode positions in the solid state. Comparison of the isolated-molecule and solid-state normal-mode calculations for a series of density functionals reveals the extent to which the inclusion of crystal-packing interactions and the relative motions between molecules are required for correctly reproducing the vibrational structure of solid-state THz spectra. The THz structure below 120 cm(-1) is a combination of both intermolecular (relative rotations and translations) and intramolecular (torsions, large amplitude motions) vibrational motions. Vibrational-mode analyses indicate that the first major feature (67.2 cm(-1)) in the PETN THz spectrum contains all of the optical rotational and translational cell modes and no internal (molecular) vibrational modes.     

Molecular vibration mode assignment of nematic liquid crystal 5CB on Terahertz spectra

Using DFT/B3LYP/6-311++G** method, the molecular structure and absorption spectra in terahertz (THz) range of liquid crystal 5CB are investigated. In a frequency range 0–15 THz, an assignment of the vibrational modes corresponding to absorption frequencies is performed using potential energy distribution for the first time. It is found that the cyano group radical (–CN) do actively take part in the strongest THz absorption of 1.743, 3.942, 5.169 and 14.769 THz in different vibration modes. The results suggest that the strong polar group should be avoided in designing liquid crystal molecule and mixtures in order to reduce the absorption intensity in THz range.     

Energy ranking of molecular crystals using density functional theory calculations and an empirical van der waals correction

By combination of high level density functional theory (DFT) calculations with an empirical van der Waals correction, a hybrid method has been designed and parametrized that provides unprecedented accuracy for the structure optimization and the energy ranking of molecular crystals. All DFT calculations are carried out using the VASP program. The van der Waals correction is expressed as the sum over atom-atom pair potentials with each pair potential for two atoms A and B being the product of an asymptotic C(6,A,B)/r(6) term and a damping function d(A,B)(r). Empirical parameters are provided for the elements H, C, N, O, F, Cl, and S. Following Wu and Yang, the C(6) coefficients have been determined by least-squares fitting to molecular C(6) coefficients derived by Meath and co-workers from dipole oscillator strength distributions. The damping functions d(A,B)(r) guarantee the crossover from the asymptotic C(6,A,B)/r(6) behavior at large interatomic distances to a constant interaction energy at short distances. The careful parametrization of the damping functions is of crucial importance to obtain the correct balance between the DFT part of the lattice energy and the contribution from the empirical van der Waals correction. The damping functions have been adjusted to yield the best possible agreement between the unit cells of a set of experimental low temperature crystal structures and their counterparts obtained by lattice energy optimization using the hybrid method. On average, the experimental and the calculated unit cell lengths deviate by 1%. To assess the performance of the hybrid method with respect to the lattice energy ranking of molecular crystals, various crystal packings of ethane, ethylene, acetylene, methanol, acetic acid, and urea have been generated with Accelrys' Polymorph Predictor in a first step and optimized with the hybrid method in a second step. In five out of six cases, the experimentally observed low-temperature crystal structure corresponds to the most stable calculated structure.     

Semi-Empirical Model to Retrieve Finite Temperature Terahertz Absorption Spectra using Morse Potential

Terahertz (THz) absorption is a fingerprint property of materials, due to the underlying low-frequency vibration/phonon modes being strongly dependent on the chemical constitutions and microscopic structures. The low excitation energies (0.414-41.4 meV) are related to two intrinsic properties of THz vibrations: the potential energy surfaces (PESs) are shallow, and the vibrationally excited states are usually populated via thermal fluctuations. The shallow PESs make the vibrations usually anharmonic, leading to redshifted vibrational excited state absorption; combined with considerable vibrational excited states population, characteristic THz signals are usually redshifted and congested with varying degrees at different temperatures. Combining existing experimental THz spectra at low temperatures, first principles vibration analysis, and the Morse potential, we developed a semi-empirical model to evaluate the anharmonicity of the low-frequency modes. The model was benchmarked with purine molecular crystal to generate THz spectra at different temperatures, the results were consistent with experiments. The good agreement suggests this model would facilitate the application of THz spectroscopy in molecular crystal characterization.     

Angle-dependent terahertz time-domain spectroscopy of amino acid single crystals

The measurement of absorption spectra using angle-dependent terahertz (THz) time-domain spectroscopy for amino acid single crystals of l-cysteine and l-histidine is reported for the first time. Linearly polarized THz radiation enables us to observe angle-dependent far-infrared absorption spectra of amino acid single crystals and determine the direction of the oscillating dipole of the molecules in the 20-100 cm(-1) range. By comparing the THz spectra of a single crystal and powder, we found that there was a clear hydrogen-bond peak in the crystal spectrum as a result of the larger hydrogen-bond network. The low-temperature THz spectra of amino acid microcrystals showed more intermolecular vibrational modes than those measured at room temperature. An ab initio frequency calculation of a single amino acid molecule was used to predict the intramolecular vibrational modes. The validity of the calculation models was confirmed by comparing the results with experimentally obtained data in the Raman spectral region.     

Application of London-type dispersion corrections to the solid-state density functional theory simulation of the terahertz spectra of crystalline pharmaceuticals

The effects of applying an empirical dispersion correction to solid-state density functional theory methods were evaluated in the simulation of the crystal structure and low-frequency (10 to 90 cm(-1)) terahertz spectrum of the non-steroidal anti-inflammatory drug, naproxen. The naproxen molecular crystal is bound largely by weak London force interactions, as well as by more prominent interactions such as hydrogen bonding, and thus serves as a good model for the assessment of the pair-wise dispersion correction term in systems influenced by intermolecular interactions of various strengths. Modifications to the dispersion parameters were tested in both fully optimized unit cell dimensions and those determined by X-ray crystallography, with subsequent simulations of the THz spectrum being performed. Use of the unmodified PBE density functional leads to an unrealistic expansion of the unit cell volume and the poor representation of the THz spectrum. Inclusion of a modified dispersion correction enabled a high-quality simulation of the THz spectrum and crystal structure of naproxen to be achieved without the need for artificially constraining the unit cell dimensions.     

Spin crossover transition of Fe(phen)2(NCS)2: periodic dispersion-corrected density-functional study

Periodic dispersion corrected DFT calculations have been performed to study the spin-crossover transition of Fe(phen)(2)(NCS)(2) in the molecular and in the crystalline state. We show that London dispersion interactions play a crucial role in the cohesion of the crystals. Based on calculations of vibrational eigenstates of the isolated molecule and of the crystalline phase in both the low- and high-spin states, the transition entropies and enthalpies have been calculated. We demonstrate that, due to the stabilization of the low-spin state by intermolecular dispersion forces, the transition enthalpy at the transition temperature is larger for the crystalline phase in comparison with an isolated molecule. The effective coordination number of the nitrogen atoms of the ligands around the iron atom has been identified as the order parameter driving the quasi-reversible low-spin to high-spin transition in the crystal. Finally, using constrained geometry relaxations at fixed values of the coordination number, we computed the energy barrier of the LS to HS transition and found it to be in a reasonable agreement with the experimental value.     

甲烷晶体的晶格能和弹性性质: 不同方法及泛函的评估

通过对甲烷晶体进行结构、晶格能和弹性特性的研究, 评估了不包含和包含色散能量修正的密度泛函理论的性能. 我们分别利用不包含色散能量修正的密度泛函理论(DFT) (包含不同的标准泛函和杂化泛函)和包含色散能量修正的密度泛函理论(DFT-D)计算了甲烷晶体特性, 并与实验作对比. 尽管DFT-D 与传统密度泛函理论及杂化密度泛函理论相比, 修正了甲烷晶体中的范德华(vdW)相互作用, 但是一些修正方案过分修正了这种相互作用. 因此, 人们在使用DFT-D方法时务必谨慎.     

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.     

Investigation of (1R,2S)‐(−)‐Ephedrine by Cryogenic Terahertz Spectroscopy and Solid‐State Density Functional Theory

The terahertz (THz) spectrum of the pharmaceutical (1R,2S)‐(?)‐ephedrine from 8.0 to 100.0 cm^?1 is investigated at liquid‐nitrogen (78.4 K) temperature. The spectrum exhibits several distinct features in this range that are characteristic of the crystal form of the compound. A complete structural analysis and vibrational assignment of the experimental spectrum is performed using solid‐state density functional theory (DFT) and cryogenic single‐crystal X‐ray diffraction. Theoretical modeling of the compound includes an array of density functionals and basis sets with the final assignment of the THz spectrum performed at a PW91/6‐311G(d,p) level of theory, which provides excellent solid‐state simulation agreement with experiment. The solid‐state analysis indicates that the seven experimental spectral features observed at low temperature consist of 13 IR‐active vibrational modes. Of these modes, nine are external crystal vibrations and provide approximately 57 % of the predicted spectral intensity. This study demonstrates that the THz spectra of complex pharmaceuticals may be well reproduced by solid‐state DFT calculations and that inclusion of the crystalline environment is necessary for realistic and accurate simulations.     

Vibrational spectra of copper sulfate hydrates investigated with low-temperature Raman spectroscopy and terahertz time domain spectroscopy

In this paper, the vibrational spectra of copper sulfate hydrates (CuSO(4)·xH(2)O, x = 5, 3, 1, 0) have been investigated with low-temperature Raman spectroscopy and terahertz time domain spectroscopy (THz-TDS). It is found that the four groups of Raman bands between 90 and 4000 cm(-1) can be assigned to lattice vibration as well as intramolecular vibrations of a copper complex, sulfate group, and water molecules. The variation of vibrational spectra during the dehydrated process are discussed in detail considering the transformation of the crystal structure, especially the bands between 3000 and 3500 cm(-1), which are attributed to the ν(1) and ν(3) modes of water molecules. In addition, as a complement of Raman spectra, the THz spectra at 0.1-3 THz indicate the absorption due to the low-frequency lattice vibration and hydrogen bond.     

PEFF: A program for the development of empirical force fields

PEFF is a new computer program designed to assist in the development of empirical force fields used in molecular mechanics calculations. Its main features are: constrained and unconstrained energy minimization available with four different techniques, rigid group refinement, crystal lattice summations, calculation of normal modes, thermodynamic functions and crystallographic temperature factors, vibrational corrections of calculated crystal structures, and a multidimensional driver to scan the energy hypersurface. Used in force field optimization mode, the program employs a least-squares method to fit the force field parameters to a set of experimental data.