Solvation of carbonaceous molecules by para-H2 and ortho-D2 clusters. II. Fullerenes

Adsorption of Helium and Hydrogen on Triphenylene and 1,3,5-Triphenylbenzene

The adsorption of helium or hydrogen on cationic triphenylene (TPL, C18H12), a planar polycyclic aromatic hydrocarbon (PAH) molecule, and of helium on cationic 1,3,5-triphenylbenzene (TPB, C24H18), a propeller-shaped PAH, is studied by a combination of high-resolution mass spectrometry and classical and quantum computational methods. Mass spectra indicate that HenTPL+ complexes are particularly stable if n = 2 or 6, in good agreement with the quantum calculations that show that for these sizes, the helium atoms are strongly localized on either side of the central carbon ring for n = 2 and on either side of the three outer rings for n = 6. Theory suggests that He14TPL+ is also particularly stable, with the helium atoms strongly localized on either side of the central and outer rings plus the vacancies between the outer rings. For HenTPB+, the mass spectra hint at enhanced stability for n = 2, 4 and, possibly, 11. Here, the agreement with theory is less satisfactory, probably because TPB+ is a highly fluxional molecule. In the global energy minimum, the phenyl groups are rotated in the same direction, but when the zero-point harmonic correction is included, a structure with one phenyl group being rotated opposite to the other two becomes lower in energy. The energy barrier between the two isomers is very small, and TPB+ could be in a mixture of symmetric and antisymmetric states, or possibly even vibrationally delocalized.

Solvation of coronene oligomers by para-H2 molecules: the effects of size and shape

The stepwise solvation of various cationic coronene oligomers by para-hydrogen (p-H₂) molecules was computationally investigated using a united-atom model for the p-H₂ molecules and the Silvera-Goldman potential, together with a polarizable description for the interaction with the hydrocarbon molecules. A survey of the energy landscape for oligomers containing between 1 and 4 coronene molecules and possible different conformers was carried out using standard global optimization, the hydrocarbon complex being kept as rigid. The most stable structures provided the starting configuration of systematic path-integral molecular dynamics simulations at 2 K. The variations of the geometric and energetic properties of the solvation shell were determined with increasing number of para-hydrogen molecules. The relative stability of the solvation shell is generally found to be more robustly determined by the energy increment (or dissociation energy) than by geometrical indicators, especially when the oligomers have less ordered structures. In agreement with recent mass spectrometry experiments, the size at which the first solvation shell is complete is found to vary approximately linearly with the oligomer size when the coronene molecules stack together, with a slope that is related to the offset between two successive molecules.

Physisorption of H2 on Fullerenes and the Solvation of C60 by Hydrogen Clusters at Finite Temperature: A Theoretical Assessment

The interaction between hydrogen and carbonaceous nanostructures is of fundamental interest in various areas of physical chemistry. In this contribution we have revisited the physisorption of hydrogen molecules and H₂ clusters on fullerenes, following a first-principles approach in which the interaction is quantitatively evaluated for the C₂₀ system using high-level electronic structure methods. Relative to coupled cluster data at the level of single, double, and perturbative triple excitations taken as a benchmark, the results for rotationally averaged physisorbed H₂ show a good performance of MP2 variants and symmetry-adapted perturbation theory, but significant deviations and basis set convergence issues are found for dispersion-corrected density functional theory. These electronic structure data are fitted to produce effective coarse-grained potentials for use in larger systems such as C₆₀-H₂. Using path-integral molecular dynamics, the potentials are also applied to parahydrogen clusters solvated around fullerenes, across the regime where the first solvation shell becomes complete and as a function of increasing temperature. For C₆₀ our findings indicate a sensible dependence of the critical solvation size on the underlying potential. As the temperature is increased, a competition is found between the surface and radial expansions of the solvation shell, with one molecule popping away at intermediate temperatures but getting reinserted at even higher temperatures.