S. H. Hilal

Calculating physical properties of organic compounds for environmental modeling from molecular structure

Mathematical models for predicting the transport and fate of pollutants in the environment require reactivity parameter values – that is the value of the physical and chemical constants that govern reactivity. Although empirical structure–activity relationships have been developed that allow estimation of some constants, such relationships are generally valid only within limited families of chemicals. The computer program, SPARC, uses computational algorithms based on fundamental chemical structure theory to estimate a large number of chemical reactivity parameters and physical properties for a wide range of organic molecules strictly from molecular structure. Resonance models were developed and calibrated using measured light absorption spectra, whereas electrostatic interaction models were developed using measured ionization pKas in water. Solvation models (i.e., dispersion, induction, H-bonding, etc.) have been developed using various measured physical properties data. At the present time, SPARC’s physical property models can predict vapor pressure and heat of vaporization (as a function of temperature), boiling point (as a function of pressure), diffusion coefficient (as a function of pressure and temperature), activity coefficient, solubility, partition coefficient and chromatographic retention time as a function of solvent and temperature. This prediction capability crosses chemical family boundaries to cover a broad range of organic compounds.

Demonstration of a consensus approach for the calculation of physicochemical properties required for environmental fate assessments

Eight software applications are compared for their performance in estimating the octanol-water partition coefficient (Kow), melting point, vapor pressure and water solubility for a dataset of polychlorinated biphenyls, polybrominated diphenyl ethers, polychlorinated dibenzodioxins, and polycyclic aromatic hydrocarbons. The predicted property values are compared against a curated dataset of measured property values compiled from the scientific literature with careful consideration given to the analytical methods used for property measurements of these hydrophobic chemicals. The variability in the predicted values from different calculators generally increases for higher values of Kow and melting point and for lower values of water solubility and vapor pressure. For each property, no individual calculator outperforms the others for all four of the chemical classes included in the analysis. Because calculator performance varies based on chemical class and property value, the geometric mean and the median of the calculated values from multiple calculators that use different estimation algorithms are recommended as more reliable estimates of the property value than the value from any single calculator.

Estimating the melting point, entropy of fusion, and enthalpy of fusion of organic compounds via SPARC.

Abstract The entropy of fusion, enthalpy of fusion, and melting point of organic compounds can be estimated through three models developed using the SPARC (SPARC Performs Automated Reasoning in Chemistry) platform. The entropy of fusion is modelled through a combination of interaction terms and physical descriptors. The enthalpy of fusion is modelled as a function of the entropy of fusion, boiling point, and flexibility of the molecule. The melting point model is the enthalpy of fusion divided by the entropy of fusion. These models were developed in part to improve SPARC’s vapour pressure and solubility models. These models have been tested on 904 unique compounds. The entropy model has a RMS of 12.5 J mol−1 K−1. The enthalpy model has a RMS of 4.87 kJ mol−1. The melting point model has a RMS of 54.4°C.

Development of Monopole Interaction Models for Ionic Compounds. Part I: Estimation of Aqueous Henry’s Law Constants for Ions and Gas Phase pKa Values for Acidic Compounds

The SPARC (SPARC Performs Automated Reasoning in Chemistry) physicochemical mechanistic models for neutral compounds have been extended to estimate Henry’s Law Constant (HLC) for charged species by incorporating ionic electrostatic interaction models. Combinations of absolute aqueous pKa values, relative pKa values in the gas phase, and aqueous HLC for neutral compounds have been used to develop monopole interaction models that quantify the energy differences upon moving an ionic solute molecule from the gas phase to the liquid phase. Inter‐molecular interaction energies were factored into mechanistic contributions of monopoles with polarizability, dipole, H‐bonding, and resonance. The monopole ionic models were validated by a wide range of measured gas phase pKa data for 450 acidic compounds. The RMS deviation error and R2 for the OH, SH, CO2H, CH3 and NR2 acidic reaction centers (C) were 16.9 kcal/mol and 0.87, respectively. The calculated HLCs of ions were compared to the HLCs of 142 ions calculated by quantum mechanics. Effects of inter‐molecular interaction of the monopoles with polarizability, dipole, H‐bonding, and resonance on acidity of the solutes in the gas phase are discussed.