Ab Initio Investigation of Structure and Cohesive Energy of Crystalline Urea

Exploring Co-Amorphous Formulations Of Nevirapine: Insights From Computational, Thermal, And Solubility Analyses

This study aimed to assess the formation of nevirapine (NVP) co-amorphs systems (CAM) with different co-formers (lamivudine-3TC, citric acid-CAc, and urea) through combined screening techniques as computational and thermal studies, solubility studies; in addition to develop and characterize suitable NVP-CAM. NVP-CAM were obtained using the quench-cooling method, and characterized by differential scanning calorimetry (DSC), X-ray diffractometry (XRD), Fourier Transform Infrared Spectroscopy (FTIR), and polarized light microscopy (PLM), in addition to in vitro dissolution in pH 6.8. The screening results indicated intermolecular interactions occurring between NVP and 3TC; NVP and CAc, where shifts in the melting temperature of NVP were verified. The presence of CAc impacted the NVP equilibrium solubility, due to hydrogen bonds. DSC thermograms evidenced the reduction and shifting of the endothermic peaks of NVP in the presence of its co-formers, suggesting partial miscibility of the compounds. Amorphization was proven by XRD and PLM assays. In vitro dissolution study exhibited a significant increase in solubility and dissolution efficiency of NVP-CAM compared to free NVP. Combined use of screening studies was useful for the development of stable and amorphous NVP-CAM, with increased NVP solubility, making CAM promising systems for combined antiretroviral therapy.

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Hybrid density functional theory (DFT) remains intractable for large periodic systems due to the demanding computational cost of exact exchange. We apply the tensor hypercontraction (THC) (or interpolative separable density fitting) approximation to periodic hybrid DFT calculations with Gaussian-type orbitals using the Gaussian plane wave approach. This is done to lower the computational scaling with respect to the number of basis functions (N) and k-points (Nk) at a fixed system size. Additionally, we propose an algorithm to fit only occupied orbital products via THC (i.e., a set of points, NISDF) to further reduce computation time and memory usage. This algorithm has linear scaling cost with k-points, no explicit dependence of NISDF on basis set size, and overall cubic scaling with unit cell size. Significant speedups and reduced memory usage may be obtained for moderately sized k-point meshes, with additional gains for large k-point meshes. Adequate accuracy can be obtained using THC-oo-K for self-consistent calculations. We perform illustrative hybrid density function theory calculations on the benzene crystal in the basis set and thermodynamic limits to highlight the utility of this algorithm.

Controlling the aqueous growth of urea crystals with different growth inhibitors: a molecular-scale study

Molecular scale understanding of the mechanism of solution-mediated nucleation and the growth of crystalline materials in the presence of growth inhibitors together with the process parameters continues to attract the interest of the scientific community though much headway has been made in recent years. Growth inhibitors can be added to solution of a crystallizing parent molecule to alter the rate of growth of different crystal faces, size and shape of the crystalline materials. In this work, we investigated the effects of a number of shape-controlling inhibitors, such as acetone, biuret and biurea, on the growth kinetics of the various faces of aqueous-grown urea crystals as a means to predictably control the crystal growth morphology. We combined the adsorption energy landscape of various auxiliaries with the kinetics of the molecular growth processes to develop an analytical model to compute the rate of growth as a function of supersaturation and the additive concentration. The model relates the kinetic and thermodynamic aspects of the adsorption of the solute, solvent and additive to provide a quantitative description of the crystal growth. Ab initio periodic dispersion-corrected density functional theory using the hybrid exchange-correlation functional was employed to determine the interfacial structure of the adsorption of various auxiliaries at crystalline surfaces. The calculated adsorption energies of different auxiliaries were employed to examine the role played by these auxiliaries during the aqueous crystallization of urea crystals containing small amounts of additives. Our results showed that the growth of (110), (111) and (1̄1̄1̄) faces were nearly unaltered by the addition of moderate amounts of acetone as it has lower adsorption energies with the surfaces of these faces. Nevertheless, the presence of acetone in the solution reasonably impeded the growth of the (001) face. The addition of biuret or biurea in the solution led to a higher adsorption energy at (001) and (111) faces. Consequently, the low concentration of these additives severely obstructed the growth of (001) and (111) faces as most of the adsorption sites were occupied by these additives. On the other hand, these additives were weakly adsorbed at the (110) face and, hence, the growth of the (110) face largely remained unaltered. Moreover, unlike biuret, biurea considerably inhibited the growth of the (1̄1̄1̄) face. Our results are in agreement with the experimental and computational results reported in the literature.

Crystal structure of propionitrile (CH3CH2CN) determined using synchrotron powder X-ray diffraction

The structure and thermal expansion of the astronomical molecule propionitrile have been determined from 100 to 150 K using synchrotron powder X-ray diffraction. This temperature range correlates with the conditions of Titan's lower stratosphere, and near surface, where propionitrile is thought to accumulate and condense into pure and mixed-nitrile phases. Propionitrile was determined to crystallize in space group, Pnma (No. 62), with unit cell a = 7.56183 (16) Å, b = 6.59134 (14) Å, c = 7.23629 (14), volume = 360.675 (13) ų at 100 K. The thermal expansion was found to be highly anisotropic with an eightfold increase in expansion between the c and b axes. These data will prove crucial in the computational modelling of propionitrile-ice systems in outer Solar System environments, allowing us to simulate and assign vibrational peaks in the infrared spectra for future use in planetary astronomy.

Cohesive properties of the crystalline phases of twenty proteinogenic α-aminoacids from first-principles calculations

Cohesive properties (lattice and cohesive energy of the crystal and corresponding sublimation enthalpy) of the complete set of twenty enantiopure anhydrous proteinogenic amino acids are investigated using first-principles calculations. In contrast to neutral amino acid molecules in the vapor phase, all amino acids form crystals in their zwitterionic form. Therefore, reliable ab initio calculations of the proton transfer energy are an indispensable step of such calculations. Simplifying procedures, designed to rationalize the computational cost of the quasi-harmonic approximation, which proves too demanding if performed fully at the given quantum level of theory, are presented and tested. For this purpose, atomic multipoles (up to the quadrupoles) for the amoeba force field are parametrized for all amino acid zwitterions. While the calculated lattice energies of the amino acids range from 235-458 kJ mol-1 in absolute value, the proton transfer energies typically amount to 100-220 kJ mol-1, which translates to sublimation enthalpies ranging from 117-202 kJ mol-1, appreciably exceeding the sublimation enthalpy values common for nonionic molecular crystals. Critically assessed experimental data on sublimation enthalpies are used as a benchmark for comparison of the data calculated in this work. Cohesive properties of most amino acids calculated in this work, combining the PBE-D3(BJ)/PAW and CCSD(T)-F12/aug-cc-pVDZ levels of theory used for predictions of the lattice energies and of the proton transfer energies, respectively, exhibit a reasonable agreement with the experiment. At the same time, this work contains the first published data on cohesive properties for several enantiopure amino acids.

Recent developments in the Inorganic Crystal Structure Database: theoretical crystal structure data and related features

The Inorganic Crystal Structure Database (ICSD) is the world's largest database of fully evaluated and published crystal structure data, mostly obtained from experimental results. However, the purely experimental approach is no longer the only route to discover new compounds and structures. In the past few decades, numerous computational methods for simulating and predicting structures of inorganic solids have emerged, creating large numbers of theoretical crystal data. In order to take account of these new developments the scope of the ICSD was extended in 2017 to include theoretical structures which are published in peer-reviewed journals. Each theoretical structure has been carefully evaluated, and the resulting CIF has been extended and standardized. Furthermore, a first classification of theoretical data in the ICSD is presented, including additional categories used for comparison of experimental and theoretical information.

Supercell calculations of the geometry and lattice energy of α-glycine crystal

Evidence about the presence of glycine in the interstellar medium (ISM) has been motivating studies aiming the understanding of the chemical behavior of this amino acid in such environment. Since glycine is expected to be predominantly found in the ISM in solid phase, this work focuses on the search for a theoretical methodology for obtaining a molecular cluster for α-glycine that provides a good description of the geometry of the unit cell and lattice energy. Calculations have been performed using the B3LYP-D3, PBE0-D3, and WB97X-D3 functionals, with def2-SVP, def2-TZVP, def2-TZVPP, and def2-QZVPP basis sets for two models: (a) the unit cell, containing 4 glycine units, and (b) the 2 × 1 × 2 expanded cell, with 16 glycine units. Corrections for the basis set superposition error have also been applied. No significant changes in geometries and lattice energy predictions from the different functionals and basis sets have been observed for each model. Nevertheless, results obtained for the larger molecular cluster are in better agreement with the experimental data. The best lattice energy prediction, obtained for the 2 × 1 × 2 supercell at the B3LYP-gCP-D3/def2-TZVPP level, is - 15.35 kcal mol^-1, with a root mean square deviation of the predicted Cartesian coordinates of the inner molecules (with respect to the experimental α-glycine unit cell geometry) of 0.966 Å. This methodology is finally recommended for future studies of similar molecular cluster, and the predicted geometry is proposed for further studies aiming to describe glycine surface reactions in the ISM.

Impact of dihydrogen bonding on lattice energies and sublimation enthalpies of crystalline [H2GaNH2]3, [H2BNH2]3 and [H2GeCH2]3

Calculated lattice energies and sublimation enthalpies provide quantitative measures of the importance of intermolecular dihydrogen bonds.

Towards reliable ab initio sublimation pressures for organic molecular crystals - are we there yet?

Knowledge of molecular crystal sublimation equilibrium data is vital in many industrial processes, but this data can be difficult to measure experimentally for low-volatility species. Theoretical prediction of sublimation pressures could provide a useful supplement to experiment, but the exponential temperature dependence of sublimation (or any saturated vapor) pressure curve makes this challenging. An uncertainty of only a few percent in the sublimation enthalpy or entropy can propagate to an error in the sublimation pressure exceeding several orders of magnitude for a given temperature interval. Despite this fundamental difficulty, this paper performs some of the first ab initio predictions of sublimation pressure curves. Four simple molecular crystals (ethane, methanol, benzene, and imidazole) have been selected for a case study showing the currently achievable accuracy of quantum chemistry calculations. Fragment-based ab initio techniques and the quasi-harmonic approximation are used for calculations of cohesive and phonon properties of the crystals, while the vapor phase is treated by the ideal gas model. Ab initio sublimation pressure curves for model compounds are compared against their experimental counterparts. The computational uncertainties are estimated, weak points of the computational methodology are identified, and further improvements are proposed.

A variable-temperature X-ray diffraction and theoretical study of conformational polymorphism in a complex organic molecule (DTC)

Two conformational crystal polymorphs of 3-diethylamino-4-(4-methoxyphenyl)-1,1-dioxo-4H-1λ⁶,2-thiazete-4-carbonitrile (DTC) have been analyzed in the 100 K-room temperature range by single crystal X-ray diffraction and high quality DFT calculations. DTC has strongly polar as well as aliphatic substituents but no hydrogen bonding groups, and thus qualifies as a test molecule for the relative importance of electrostatic vs. dispersion–repulsion terms. The two polymorphs have the same P2₁/n space group and differ by a flipping of the –OCH₃ group, the two conformations being almost equi-energetic and separated by a low barrier. The system is monotropic in the observed temperature range with nearly identical thermal expansion coefficients and energy–temperature slopes, one phase consistently predicted to be more stable in agreement with the relative ease of appearance. Energy decompositions show that the electrostatic term is dominant and stabilizes with decreasing temperature. Dispersion and repulsion show the expected behavior, the former becoming more stabilizing at lower temperature in contrast with increasing repulsion at higher density. Absolute values and trends are very similar in the two phases, explaining the small total energy difference. Geometrical analyses of intermolecular contacts using fingerprint plots, as well as the study of molecular dipole moments as a function of T in the framework of the Quantum Theory of Atoms in Molecules, reveal more details of phase stability.

A small conformational change in the asymmetric unit has a significant effect on how non-covalent interactions determine (i) the crystal packing and (ii) the effect of T on the relative balance of electrostatics and dispersion–repulsions.

Specific features of supramolecular organisation and hydrogen bonding in proline cocrystals: a case study of fenamates and diclofenac

New zwitterionic cocrystals of fenamate drugs and diclofenac with the naturally occurring amino acid l-proline have been obtained and thoroughly characterised by a variety of experimental and theoretical techniques.

The prediction of far-infrared spectra for planetary nitrile ices using periodic density functional theory with comparison to thin film experiments

The application of high accuracy periodic DFT calculations to replicate laboratory thin film spectra allowing the assignment of vibrational far-infrared modes in nitrile-bearing planetary ices.

Role of tailor-made additives in controlling vapour growth asymmetry along the polar axis of α-resorcinol crystals: a molecular-scale study

The role of different tailor-made additives and conformers of resorcinol have been investigated to elucidate the mechanism of unidirectional vapour growth of α-resorcinol crystal. Our results show the growth anisotropy is linked with the surface stabilization of the slow growing face.

State-of-the-Art Calculations of Sublimation Enthalpies for Selected Molecular Crystals and Their Computational Uncertainty

A computational methodology for calculation of sublimation enthalpies of molecular crystals from first principles is developed and validated by comparison to critically evaluated literature experimental data. Temperature-dependent sublimation enthalpies for a set of selected 22 molecular crystals in their low-temperature phases are calculated. The computational methodology consists of several building blocks based on high-level electronic structure methods of quantum chemistry and statistical thermodynamics. Ab initio methods up to the coupled clusters with iterative treatment of single and double excitations and perturbative triples correction with an estimated complete basis set description [CCSD(T)/CBS] are used to calculate the cohesive energies of crystalline phases within a fragment-based additive scheme. Density functional theory (DFT) calculations with periodic boundary conditions (PBC) coupled with the quasi-harmonic approximation are used to evaluate the thermal contributions to the enthalpy of the solid phase. The properties of the vapor phase are calculated within the ideal-gas model using the rigid-rotor harmonic-oscillator model with correction for internal rotation using a one-dimensional hindered rotor approximation and a proper treatment of the molecular rotational degrees of freedom in the vicinity of 0 K. All individual terms contributing to the sublimation enthalpy as a function of temperature are discussed and their uncertainties estimated by comparison to critically evaluated experimental data.

Thermal properties of molecular crystals through dispersion-corrected quasi-harmonic ab initio calculations: the case of urea

An ab initio quantum-mechanical theoretical framework is presented to compute the thermal properties of molecular crystals. The present strategy combines dispersion-corrected density-functional-theory (DFT-D), harmonic phonon dispersion, quasi-harmonic approximation to the lattice dynamics for thermal expansion and thermodynamic functions, and quasi-static approximation for anisotropic thermo-elasticity. The proposed scheme is shown to reliably describe thermal properties of the urea molecular crystal by a thorough comparison with experimental data.

Simulating growth morphology of urea crystals from vapour and aqueous solution

The prediction of solvent-induced urea crystal growth shape from the internal and interfacial atomic structure was investigated. A computational model was used to calculate crystal growth morphologies from vapour and aqueous solution as a function of supersaturation; predicted growth shapes agree well with experimental observations.

The applicability of the dimeric heterosynthon concept to molecules with equivalent binding sites. A DFT study of crystalline urea–H2O2

The limited applicability of the dimeric heterosynthon concept to a two-component urea–H₂O₂crystal is reported.

Cl···Cl interactions in molecular crystals: insights from the theoretical charge density analysis

The structure, IR harmonic frequencies and intensities of normal vibrations of 20 molecular crystals with the X-Cl···Cl-X contacts of different types, where X = C, Cl, and F and the Cl···Cl distance varying from ~3.0 to ~4.0 Å, are computed using the solid-state DFT method. The obtained crystalline wave functions have been further used to define and describe quantitatively the Cl···Cl interactions via the electron-density features at the Cl···Cl bond critical points. We found that the electron-density at the bond critical point is almost independent of the particular type of the contact or hybridization of the ipso carbon atom. The energy of Cl···Cl interactions, E(int), is evaluated from the linking E(int) and local electronic kinetic energy density at the Cl···Cl bond critical points. E(int) varies from 2 to 12 kJ/mol. The applicability of the geometrical criterion for the detection of the Cl···Cl interactions in crystals with two or more intermolecular Cl···Cl contacts for the unique chlorine atom is not straightforward. The detection of these interactions in such crystals may be done by the quantum-topological analysis of the periodic electron density.

A post-Hartree-Fock study of pressure-induced phase transitions in solid nitrogen: the case of the α, γ, and ε low-pressure phases

We have performed fully periodic ab initio post-Hartree-Fock calculations on three low-temperature molecular phases of solid nitrogen, namely α, γ, and ε. The aim is to characterize the pressure-induced transitions among these phases, which are prototypical molecular crystals, through a periodic local-MP2 method, implemented in the CRYSCOR program. Cohesive energies are computed using extended correlation consistent molecular basis sets, up to quintuple-ζ quality. The MP2 description of equilibrium volumes, pressure-volume curves, and transition pressures is found to be in good agreement with the experiments.

Conformational polymorphism in a Schiff-base macrocyclic organic ligand: an experimental and theoretical study

Polymorphism in the highly flexible organic Schiff-base macrocycle ligand 3,6,9,17,20,23-hexa-azapentacyclo(23.3.1.1(11,15).0(2,6).0(16,20))triaconta-1(29),9,11,13,15(30),23,25,27-octaene (DIEN, C(24)H(30)N(6)) has been studied by single-crystal X-ray diffraction and both solid-state and gas-phase density functional theory (DFT) calculations. In the literature, only solvated structures of the title compound are known. Two new polymorphs and a new solvated form of DIEN, all obtained from the same solvent with different crystallization conditions, are presented for the first time. They all have P\bar 1 symmetry, with the macrocycle positioned on inversion centres. The two unsolvated polymorphic forms differ in the number of molecules in the asymmetric unit Z', density and cohesive energy. Theoretical results confirm that the most stable form is (II°), with Z' = 1.5. Two distinct molecular conformations have been found, named `endo' or `exo' according to the orientation of the imine N atoms, which can be directed towards the interior or the exterior of the macrocycle. The endo arrangement is ubiquitous in the solid state and is shared by two independent molecules which constitute an invariant supramolecular synthon in all the known crystal forms of DIEN. It is also the most stable arrangement in the gas phase. The exo form, on the other hand, appears only in phase (II°), which contains both the conformers. Similarities and differences among the occurring packing motifs, as well as solvent effects, are discussed with the aid of Hirshfeld surface fingerprint plots and correlated to the results of the energy analysis. A possible interconversion path in the gas phase between the endo and the exo conformers has been found by DFT calculations; it consists of a two-step mechanism with activation energies of the order of 30-40 kJ mol(-1). These findings have been related to the empirical evidence that the most stable phase (II°) is also the last appearing one, in accordance with Ostwald's rule.

Density functional study of the proton transfer effect on vibrations of strong (short) intermolecular O-H...N/O-...H-N+ hydrogen bonds in aprotic solvents

The structure and spectroscopic properties of the 1:1 complexes of substituted pyridines with benzoic acid and phenol derivatives in aprotic solvents are studied using B3LYP functional combined with the polarizable continuum model approximation. Two extreme structures are investigated: the state without (HB) and with proton transfer (PT). In the presence of an external electric field the O...N distance is contracted and the PT state does appear. The PT state of both the pyridine-benzoic and the pyridine-phenol complexes displays the only IR-active band in the 2800-1800 frequency region, which is located around 2000 cm(-1). However, the nature of the band is different for these two complexes. In the pyridine-benzoic acid complex it is practically a pure stretching vibration of the HN(+) group, while in the pyridine-phenol complex it is the mixed vibration of the bridging proton. A specific feature of the PT state in the pyridine-phenol complex is an IR-intensive band near 600 cm(-1), associated with the asymmetric stretching vibrations of the O(-)...HN(+) fragment. Its intensity is reciprocally proportional to the O...N distance. The appearance of this band provides an efficient criterion to differentiate between the HB and PT states of the 1:1 complexes of phenols with pyridines in aprotic solvents.

DFT Study of Proton Dynamics in the Potassium Hydrogen Maleate Crystal: the Infrared Versus the Inelastic Neutron Scattering Spectra

The structure, harmonic frequencies and infrared (IR) intensities of the fundamental transitions in potassium hydrogen maleate (KHM) crystal have been computed using the BLYP/6-31G** method with periodic boundary conditions. The inelastic neutron scattering (INS) intensities of the fundamental transitions were evaluated and compared with the experimental INS spectra of hydrogenated sample and isotopic analogues of the KHM crystal; the agreement is reasonable good. Spectroscopic manifestations of the proton dynamics in the strong intramolecular H-bond were clarified by comparison of the computed frequencies of the IR and INS-active fundamental transitions. Due to several reasons, the number of the bands and their relative intensity were found to differ between the INS and IR spectra of the KHM crystal. Reduced masses for the asymmetric O…H…O stretching vibrations of the hydrogen maleate ion appear to be larger 1 a.m.u. because of the strong coupling between the bridging proton motion and the various intramolecular vibrations.