Monte Carlo simulations of pure liquid substituted benzenes with OPLS potential functions

Intermolecular potential functions have been developed for use in computer simulations of substituted benzenes. Previously reported optimized potentials for liquid simulations (OPLS) for benzene and organic functional groups were merged and tested in Monte Carlo statistical mechanics simulations for the pure liquids of toluene, m‐cresol, anisole, aniline, and benzonitrile at 25°C at 1 atm. The merged potential functions yielded acceptable thermodynamic results for the liquids except in the case of aniline, for which the error in the heat of vaporization was 12%. This was remedied by enhancing the polarity of the model to be more consistent with the observed dipole moment of aniline. Overall, the average errors in computed heats of vaporization and densities were then 2 and 1%, respectively. The structures of the liquids were characterized through energy and radial distribution functions. For m‐cresol and aniline, the molecules participate in averages of 1.6 and 1.4 hydrogen bonds, respectively. Condensed phase effects on the torsional energies for anisole, m‐cresol, and aniline were found to be small; m‐cresol has a slightly enhanced tendency to be nonplanar in the liquid than in the gas phase, while anisole shows the opposite pattern. © 1993 John Wiley & Sons, Inc.

[1]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[2]  D. G. Lister,et al.  The microwave spectrum, structure and dipole moment of aniline , 1974 .

[3]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid alcohols , 1986 .

[4]  J. Evans.  The vibrational assignments and configuration of aniline, aniline-NHD and aniline-ND2 , 1960 .

[5]  J. D. Cox,et al.  Thermochemistry of organic and organometallic compounds , 1970 .

[6]  Douglas A. Smith,et al.  Structural studies of aromatic amines and the dna intercalating compounds m‐amsa and o‐amsa: Comparison of mndo, am1, and pm3 to experimental and ab initio results , 1992 .

[7]  Richard B. Silverman,et al.  The Organic Chemistry of Drug Design and Drug Action , 1992 .

[8]  Peter A. Kollman,et al.  Methoxy group nonplanarity in o-dimethoxybenzenes. Simple predictive models for conformations and rotational barriers in alkoxyaromatics , 1979 .

[9]  D. J. Ruxton,et al.  Thermodynamic properties of organic oxygen compounds. Part 18.—Vapour heat capacities and heats of vaporization of ethyl ketone, ethyl propyl ketone, methyl isopropy1 ketone, and methyl phenyl ether , 1967 .

[10]  William L. Jorgensen,et al.  Monte Carlo simulations of the hydration of substituted benzenes with OPLS potential functions , 1993, J. Comput. Chem..

[11]  William L. Jorgensen,et al.  Aromatic-aromatic interactions: free energy profiles for the benzene dimer in water, chloroform, and liquid benzene , 1990 .

[12]  William L. Jorgensen,et al.  Monte Carlo simulations of liquid acetonitrile with a three-site model , 1988 .

[13]  William L. Jorgensen,et al.  Theoretical studies of medium effects on conformational equilibria , 1983 .

[14]  William L. Jorgensen,et al.  Enhanced view of structure and binding for cyclophane-arene complexes through joint theoretical and experimental study , 1992 .

[15]  William L. Jorgensen,et al.  Relative partition coefficients for organic solutes from fluid simulations , 1990 .

[16]  William L. Jorgensen,et al.  Optimized intermolecular potential functions for liquid hydrocarbons , 1984 .

[17]  F. Diederich.  Supramolecular Cyclophane Chemistry , 1992 .

[18]  Peter A. Kollman,et al.  Theoretical investigations of the rotational barrier in anisole: an ab initio and molecular dynamics study , 1990 .

[19]  William L. Jorgensen,et al.  Temperature and size dependence for Monte Carlo simulations of TIP4P water , 1985 .

[20]  William L. Jorgensen,et al.  Convergence of Monte Carlo simulations of liquid water in the NPT ensemble , 1982 .

[21]  W. L. Jorgensen,et al.  Monte Carlo simulations of liquid alkyl ethers with the OPLS potential functions , 1990 .

[22]  Robert L. Kuczkowski,et al.  Molecular structures of gas‐phase polyatomic molecules determined by spectroscopic methods , 1979 .