scholarly journals The Use of Quantum-Chemical Semiempirical Methods to Calculate the Lattice Energies of Organic Molecular Crystals. Part III: The Lattice Energy of Borazine (B3N3H6) and its Packing in the Solid State*

2004 ◽  
Vol 59 (9) ◽  
pp. 609-614 ◽  
Author(s):  
Gerhard Raabe

A previously presented quantum-chemical scheme has been used to calculate the lattice energies of borazine (B3N3H6), the low pressure polymorph of benzene (C6H6), and of borazine in the lowpressure benzene lattice utilizing some frequently used semiempirical methods (CNDO/2, INDO, MINDO/3, MNDO, AM1, PM3, MSINDO). With all methods the lattice energy of the title compound was found to be less favourable than that of isoelectronic benzene, which offers an explanation of the significantly lower melting point of B3N3H6. Calculation of the lattice energy of borazine in the crystal lattice of the low-pressure modification of benzene revealed that the interactions between the molecules in this environment are not so stabilizing as those in its own lattice. This is predominantly due to a less favourable contribution of the dispersion energy. The semiempirical results have qualitatively been confirmed by quantum-chemical calculations on small molecular clusters at the MP2/6-31+G*//HF/6-31+G* level of ab initio theory. In these calculations we assumed pairwise additivity of the intermolecular interactions and calculated the energy of interaction between a reference molecule and all those neighbours to which the shortest intermolecular distance does not exceed 3Å.

2000 ◽  
Vol 55 (6-7) ◽  
pp. 609-615 ◽  
Author(s):  
Gerhard Raabe

Abstract A method to calculate the lattice energies of organic molecular crystals is described. It is based on the semiempirical quantum chemical MINDO/3 approximation but might also be used within the framework of any other quantum chemical method. The lattice energy is approximated by the sum of dispersion-, induction-, exchange repulsion-, and electrostatic energy. Different, however, from other schemes employed in this field, like for example the atom-atom-potential method, the variables in the expression for the lattice energy have not been fitted to reproduce experimental values and, therefore, the single contributions retain their original physical meaning. Moreover, the method offers the advantage that it may be directly applied to all compounds that can be treated within the framework of the underlying quantum chemical method. Thus, time consuming readjustment of the entire parameter set upon extension of the group of target molecules by another class of compounds becomes obsolete. As an example, the lattice energies of the three polymorphs of glycine are calculated.


2002 ◽  
Vol 57 (12) ◽  
pp. 961-966 ◽  
Author(s):  
Gerhard Raabe

The lattice energies (ΔE lat) of α- and β-oxalic acid ((COOH)2) have been calculated using a recently introduced semiempirical quantum-chemical procedure. Within the framework of this method the lattice energy (ΔE lat) is evaluated as the sum of the semiempirically calculated intermolecular dispersion (ΔEdis), induction (ΔEind), repulsion (ΔErep), and electrostatic energy (ΔEels). The lattice energies of the two polymorphs of oxalic acid obtained in this way correlate not only with the results of other calculations but also with the experimentally determined heats of sublimation in that the α-modification, which has a somewhat higher heat of sublimation, is slightly more stable than the β-polymorph. However, additional quantum-chemical calculations at the non-empirical ab initio level (e. g. ZPE+MP2(FC)/6-311++G**//MP2(FC)/6-311++G**) revealed that the absolute values of the lattice energies and the heats of sublimation are not directly comparable to each other because the structures of the (COOH)2 molecules in the crystal lattices of both polymorphs differ significantly from that of the most stable form of the free molecule in the gasphase. At about 4.3 kcal/mol the calculated energy difference between the structure of the molecule in the solid state and the energetically most favourable conformation of the free (COOH)2 molecule in the gasphase is much smaller than that in the recently described case of α-glycine (28±2 kcal/mol). However, even such a small difference might be the source of serious problems if the heats of sublimation are employed to fit parameters to be used in the optimization of crystal packings.


1999 ◽  
Vol 54 (10-11) ◽  
pp. 611-616 ◽  
Author(s):  
Gerhard Raabe

A new value for the lattice energy of a-glycine was determined by combination of the experimentally measured heat of sublimation taken from literature and the quantum-chemically calculated energy difference Etot,gp - Etot,cry, where Etot,gp is the total energy of the most stable form of the compound in the gas phase (carboxylic acid) and Etot,cry the total energy of the molecule as it occurs in its crystalline form (betaine). At the highest levels of ab initio theory employed in this study this energy difference is -(28±2) kcal/mol, indicating that older work overestimated this difference significantly. The reason for the overestimation of this energy difference was determined by means of additional ab initio calculations. The lattice energy of -(67±2)kcal/mol obtained using the new value for Etot,gp - Etot,cry is significantly more positive than an older value of -103 kcal/mol frequently cited in the literature.


2017 ◽  
Author(s):  
Robson de Farias

In the present work, the reliability of the volume-based thermodynamics (VBT) methods in the calculation of lattice energies is investigated by applying the “traditional” Kapustinskii equation [8], as well as Glasser-Jenkins [3] and Kaya [5] equations to calculate the lattice energies for Na, K and Rb pyruvates [9-11] as well as for the coordination compound [Bi(C<sub>7</sub>H<sub>5</sub>O<sub>3</sub>)<sub>3</sub>C<sub>12</sub>H<sub>8</sub>N<sub>2</sub>] [17] (in which C<sub>12</sub>H<sub>8</sub>N<sub>2</sub> = 1,10 phenathroline and C<sub>7</sub>H<sub>5</sub>O<sub>3</sub><sup>-</sup>= <i>o</i>-hyddroxybenzoic acid anion). As comparison, the lattice energies are also calculated using formation enthalpy values for sodium pyrivate and [Bi(C<sub>7</sub>H<sub>5</sub>O<sub>3</sub>)<sub>3</sub>C<sub>12</sub>H<sub>8</sub>N<sub>2</sub>]. For the pyruvates, is verified that none of the considered approach, Kapustinskii, Glasser, Kaya or density, provides values that agrees in an acceptable % difference, with the lattice energy values calculated from the formation enthalpy values. However, it must be pointed out that Kaya approach, with deals with a chemical hardness approach is the better one for such kind of inorganic-organic salts. Based on data obtained for [Bi(C<sub>7</sub>H<sub>5</sub>O<sub>3</sub>)<sub>3</sub>C<sub>12</sub>H<sub>8</sub>N<sub>2</sub>] is concluded that the only one VBT method that provides reliable lattice energies for compounds with bulky uncharged ligands is that one based on density values (derived by Glasser-Jenkins).


Author(s):  
Edward O. Pyzer-Knapp ◽  
Hugh P. G. Thompson ◽  
Graeme M. Day

We present a re-parameterization of a popular intermolecular force field for describing intermolecular interactions in the organic solid state. Specifically we optimize the performance of the exp-6 force field when used in conjunction with atomic multipole electrostatics. We also parameterize force fields that are optimized for use with multipoles derived from polarized molecular electron densities, to account for induction effects in molecular crystals. Parameterization is performed against a set of 186 experimentally determined, low-temperature crystal structures and 53 measured sublimation enthalpies of hydrogen-bonding organic molecules. The resulting force fields are tested on a validation set of 129 crystal structures and show improved reproduction of the structures and lattice energies of a range of organic molecular crystals compared with the original force field with atomic partial charge electrostatics. Unit-cell dimensions of the validation set are typically reproduced to within 3% with the re-parameterized force fields. Lattice energies, which were all included during parameterization, are systematically underestimated when compared with measured sublimation enthalpies, with mean absolute errors of between 7.4 and 9.0%.


2014 ◽  
Vol 92 (11) ◽  
pp. 1111-1117
Author(s):  
Xueli Zhang ◽  
Xuedong Gong

Nitrogen-rich 1,2,4-triazole (1) and 1,2,3-triazole (2) react as bases with the oxygen-rich acids HNO3 (a), HN(NO2)2 (b), and HClO4 (c) to produce energetic salts (1a, 1b, and 1c and 2a, 2b, and 2c, respectively) potentially applicable to composite explosives and propellants. In this study, these salts were studied with the dispersion-corrected density functional theory. For the isomers such as 1a and 2a, the more negative ΔrGm of the formation reaction leads to a higher thermally stable salt. The ability to form intramolecular hydrogen bonds predicted with the quantum theory of atoms in molecules has the order of 2 > 1. Different hydrogen bonds result in different second-order perturbation energies, redshifts in IR, and electron density differences. The charge transfer, binding energy, dispersion energy, lattice energy, and energy gap between frontier orbits in the salts of 1 are larger than those of 2, which is helpful for stabilizing the former, and 1 is more obviously stabilized than 2 by formation of salts. Different conformations of 1 and 2 hardly affect the frontier orbital distributions. Base 1 is a more preferred base than 2 to form salts.


The thermochemical data of part I, the heats of formation and solution of the alkali-metal (group 1 a ) azides, are used in conjunction with other data to derive values for the lattice energies of alkali-metal azides, the heat of formation of the azide radical, for the electron affinity and hydration heat of the azide ion. Calculations by previous workers of these magnitudes, which are not of course susceptible to direct measurement, have generally been erroneous. The lattice energies of the alkali azides (kcal mole -1 ) are: LiN 3 , 194; NaN 3 , 175; KN 3 , 157; RbN 3 , 152; CsN 3 , 146. For potassium, rubidium and caesium azides a term-by-term theoretical calculation of the lattice energy which allows for the non-spherical character of the azide ion supports these figures, which are based on experimental data of part I. The standard enthalpy of formation of the azide radical ∆ H 0 f (N 3G ) is estimated to be 116 kcal mole -1 . The electron affinity of the azide radical E (N 3G ) is 81 kcal mole -1 . These figures permit the evaluation of other lattice energies and the following values (kcal mole -1 ) have been obtained: NH 4 N 3 , 175; CuN 3 , 227; AgN 3 , 205; TlN 3 , 163·5; CaN 6 , 517; SrN 6 , 494; BaN 6 , 469 and PbN 6 , 516. From the enthalpy of formation of the azide radical the bond dissociation energies D ( X — N 3 ) in some covalent azides may be derived. D (H — N 3 ) is 96 kcal mole -1 and D (C— N 3 ) in aliphatic azides is about 83 kcal mole -1 .


Sign in / Sign up

Export Citation Format

Share Document