This article shows how to evaluate rotational symmetry numbers for different molecular configurations and how to apply them to transition state theory. In general, the symmetry number is given by the ratio of the reactant and transition state rotational symmetry numbers. However, special care is advised in the evaluation of symmetry numbers in the following situations: (i) if the reaction is symmetric, (ii) if reactants and/or transition states are chiral, (iii) if the reaction has multiple conformers for reactants and/or transition states and, (iv) if there is an internal rotation of part of the molecular system. All these four situations are treated systematically and analyzed in detail in the present article. We also include a large number of examples to clarify some complicated situations, and in the last section we discuss an example involving an achiral diasteroisomer.

We have developed an evolutionary algorithm (EA) for the global minimum search of molecular clusters. The EA is able to discover all the putative global minima of water clusters up to (H(2)O)(20) and benzene clusters up to (C(6)H(6))(30). Then, the EA was applied to search for the global minima structures of (C(6)H(6))(n)(+) with n = 2-20, some of which were theoretically studied for the first time. Our results for n = 2-6 are consistent with previous theoretical work that uses a similar interaction potential. Excluding the very symmetric global minimum structure for n = 9, the growth pattern of (C(6)H(6))(n)(+) with n ≥ 7 involves the (C(6)H(6))(2)(+) dimer motif, which is placed off-center in the cluster. Such observation indicates that potentials commonly used in the literature for (C(6)H(6))(n)(+) cannot reproduce the icosahedral-type packing suggested by the available experimental data.

No abstract

Coronene is one of the basic polycyclic aromatic hydrocarbons (PAHs) used to test the reliabilty of a multidimensional potential energy surface (PES) and to assess its influence on the formation dynamics of PAH clusters with defined physical and chemical properties. We report an analytical potential energy surface for modeling the coronene− coronene interaction, whose parameters were fine-tuned on dispersion-corrected DFT estimations performed within the generalized gradient PBE approximation, that is suitable for describing molecular aggregates involving aromatic species. This model was used to build a potential function for coronene clusters (Cor n ) that was then applied in a detailed global geometry optimization study with an evolutionary algorithm. A large variety of low-energy structures were obtained for the Cor n (n = 2−15) clusters ranging from columnar-type to two-stacked in a handshake association motifs. Moreover, it was found that a transition from a single-stack columnar regime to other more complex shapes occurs at n = 6, whereas previous results based on a simpler coarse-grained potential pointed to a transition at n = 8. Geometry reoptimizations were also performed at the DFT level for the most representative low-energy structures of Cor n (n = 3−6), which confirmed the reliability of the present findings.

No abstract