Modeling benzene with single-site potentials from ab initio calculations: a step toward hybrid models of complex molecules.

Extensive ab initio calculations at the MP2/6-31G* level have been carried out to sample the energy surface for the interactions of the benzene dimers. This database has been used to parameterize two anisotropic single-site models, meant to be used as building blocks in hybrid models of complex, liquid crystal forming molecules. A quadrupolar Gay-Berne (GBQIII) and an S-function (SF) Corner potentials have been obtained in this way. Their ability to reproduce, qualitatively at least, the phase diagram as well as energetic and structural properties of benzene has been tested with Monte Carlo simulations and compared with previous literature potentials, GBQI [S. Gupta et al., Mol. Phys. 65, 961 (1988)] and GBQII [T. R. Walsh, Mol. Phys. 100, 2867 (2002)]. It turned out that GBQI showed no melting transition in the temperature range explored (100-400 K), while GBQII underwent a phase transition from solid to gas, with no liquid phase. Conversely, both models parameterized on our database of ab initio interaction energies (GBQIII and SF) gave rise to a stable liquid phase. Melting has been observed between 100 and 150 K (GBQIII) and in the range 300-350 K (SF), i.e., substantially below and slightly above the experimental value at ambient pressure, 278 K. The description of the crystal structure of benzene at atmospheric pressure is also in better agreement with experimental data if the SF model is used, while positional correlations in the liquid are better described by the GBQIII potential. The S-function potential is also computationally more convenient. These results could be useful in the semirealistic modeling of more complex molecules.

O 1s --> sigma* resonance in O2: inadequacy of only two exchange-split components.

The O 1s-->sigma* transition below the O K-edge in O2 has been investigated by absorption, constant ionic state (CIS) experiments, and extensive configuration interaction calculations. CIS scans of the three lowest-lying final states reached in resonant Auger decay provide a wealth of information on energy range, symmetry, and spin multiplicity of the intermediate states with sigma* character. We conclude that the identification of only two exchange-split components is inadequate because a complex manifold of states with sigma* character exists with no unique energy difference between related states.