Adsorption of hydrogen on neutral and charged fullerene: experiment and theory

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.

Adsorption of Helium on Small Cationic PAHs: Influence of Hydrocarbon Structure on the Microsolvation Pattern

The adsorption of up to ∼100 helium atoms on cations of the planar polycyclic aromatic hydrocarbons (PAHs) anthracene, phenanthrene, fluoranthene, and pyrene was studied by combining helium nanodroplet mass spectrometry with classical and quantum computational methods. Recorded time-of-flight mass spectra reveal a unique set of structural features in the ion abundance as a function of the number of attached helium atoms for each of the investigated PAHs. Path-integral molecular dynamics simulations were used with a polarizable potential to determine the underlying adsorption patterns of helium around the studied PAH cations and in good general agreement with the experimental data. The calculated structures of the helium-PAH complexes indicate that the arrangement of adsorbed helium atoms is highly sensitive toward the structure of the solvated PAH cation. Closures of the first solvation shell around the studied PAH cations are suggested to lie between 29 and 37 adsorbed helium atoms depending on the specific PAH cation. Helium atoms are found to preferentially adsorb on these PAHs following the 3 × 3 commensurate pattern common for graphitic surfaces, in contrast to larger carbonaceous molecules like corannulene, coronene, and fullerenes that exhibit a 1 × 1 commensurate phase.

Carbon Nanostructures Doped with Transition Metals for Pollutant Gas Adsorption Systems

The adsorption of molecules usually increases capacity and/or strength with the doping of surfaces with transition metals; furthermore, carbon nanostructures, i.e., graphene, carbon nanotubes, fullerenes, graphdiyne, etc., have a large specific area for gas adsorption. This review focuses on the reports (experimental or theoretical) of systems using these structures decorated with transition metals for mainly pollutant molecules' adsorption. Furthermore, we aim to present the expanding application of nanomaterials on environmental problems, mainly over the last 10 years. We found a wide range of pollutant molecules investigated for adsorption in carbon nanostructures, including greenhouse gases, anticancer drugs, and chemical warfare agents, among many more.

Long-range order in quadrupolar systems on spherical surfaces

The interplay between curvature, confinement and ordering on curved manifolds, with anisotropic interactions between building blocks, takes a central role in many fields of physics. In this paper, we investigate the effects of lattice symmetry and local positional order on orientational ordering in systems of long-range interacting point quadrupoles on a sphere in the zero temperature limit. Locally triangular spherical lattices show long-range ordered quadrupolar configurations only for specific symmetric lattices as strong geometric frustration prevents general global ordering. Conversely, the ground states on Caspar-Klug lattices are more diverse, with many different symmetries depending on the position of quadrupoles within the fundamental domain. We also show that by constraining the quadrupole tilts with respect to the surface normal, which models interactions with the substrate, and by considering general quadrupole tensors, we can manipulate the ground state configuration symmetry.

The Acid-Base Through-the-Cage Interaction as an Example of an Inversion in a Cage Isomerism

We define a new inversion in a cage isomerism (ic): X@C⋯Y₪icY@C⋯X, (₪ is the isomerism relation) as an isomerism in the three-component system of molecules X, Y, and a cage C, in which one of the molecules is located inside and the other outside the cage. The ic isomerism is similar to the endo-exo one, which occurs only if either the interior or exterior of C is empty. By contrast, ic occurs only if neither the interior nor the exterior of C is empty. We also discuss the other closely related types of isomerisms are also discussed. Calculations of the XH⋯NH3@C60 and NH3⋯HX@C60ic isomers were performed at the ωB97XD/Def2TZVP level. The calculated energies demonstrated that the systems with the HX acid outside (X = F, Cl) and the NH3 base inside the cage, XH⋯NH3@C60, are more stable than their ic isomers, NH3⋯HX@C60, by about 4–8 kcal/mol. This is because NH3 is more stabilized inside the cage than HX (a matter of 6.5 kcal/mol). In the studied systems and subsystems, the HX molecules are Lewis acids and the NH3 molecule is always a Lewis base. The C60 molecule with HX inside or outside the cage is also an acid for the NH3 base positioned outside or inside the cage. On the other hand, the C60 cage is truly amphoteric because it is simultaneously an acid and a base.

The Structure of Adamantane Clusters: Atomistic vs. Coarse-Grained Predictions From Global Optimization

Candidate structures for the global minima of adamantane clusters, (C₁₀H₁₆)_N, are presented. Based on a rigid model for individual molecules with atom-atom pairwise interactions that include Lennard-Jones and Coulomb contributions, low-energy structures were obtained up to N = 42 using the basin-hopping method. The results indicate that adamantane clusters initially grow accordingly with an icosahedral packing scheme, followed above N = 14 by a structural transition toward face-centered cubic structures. The special stabilities obtained at N = 13, 19, and 38 are consistent with these two structural families, and agree with recent mass spectrometry measurements on cationic adamantane clusters. Coarse-graining the intermolecular potential by averaging over all possible orientations only partially confirm the all-atom results, the magic numbers at 13 and 38 being preserved. However, the details near the structural transition are not captured well, because despite their high symmetry the adamantane molecules are still rather anisotropic.

Charged Clusters of C60 and Au or Cu: Evidence for Stable Sizes and Specific Dissociation Channels

We have doped helium nanodroplets with C60 and either gold or copper. Positively or negatively charged (C60) mM n± ions (M = Au or Cu) containing up to ≈10 fullerenes and ≈20 metal atoms are formed by electron ionization. The abundance distributions extracted from high-resolution mass spectra reveal several local anomalies. The sizes of the four most stable (C60) mAu n± ions identified in previous calculations for small values of m and n ( m ≤ 2 and n ≤ 2, or m = 1 and n = 3) agree with local maxima in the abundance distributions. Our data suggest the existence of several other relatively stable ions including (C60)2Au3± and (C60)3Au4-. Another feature, namely the absence of bare (C60)2±, confirms the prediction that (C60)2M± dissociates by loss of C60± rather than loss of M. The experimental data also reveal the preference for loss of (charged or neutral) C60 over loss of a metal atom from some larger species such as (C60)3M3+. In contrast to these similarities between Au and Cu, the abundance distributions of (C60)3Au n- and (C60)3Cu n- are markedly different. In this discussion, we emphasize the similarities and differences between anions and cations, and between gold and copper. Also noteworthy is the observation of dianions (C60) mAu n2- for m = 2, 4, and 6.

Peculiar covalent bonding of C60/6H-SiC(0 0 0 1)-(3 × 3) probed by photoelectron spectroscopy

High resolution photoemission with synchrotron radiation was used to study the interface formation of a thin layer of C₆₀ on 6H-SiC(0 0 0 1)-(3 × 3), characterized by protruding Si-tetramers. The results show that C₆₀ is chemisorbed by orbital hybridization between the highest-occupied molecular orbital (HOMO) and the p _z orbital of Si adatom at the apex of the tetramers. The covalent nature of the bonding was inferred from core level as well as valence band spectra. The Si 2p spectra reveal that a large fraction (at least 45%) of the Si adatoms remain unbound despite the reactive character of the associated dangling bonds. This is consistent with a model in which each C₆₀ is attached to the substrate through a single covalent C₆₀-Si bond. A binding energy shift of the core levels associated with sub-surface Si or C atoms indicates a decrease of the SiC band bending caused by a charge transfer from the C₆₀ molecules to the substrate via the formation of donor-like interface states.

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.

Adsorption of molecular hydrogen on coronene with a new potential energy surface

Adsorption of molecular hydrogen on coronene studied with a new potential energy surface. Path integral Monte Carlo and basin-hopping calculations have been performed to investigate energies and structures of the corresponding (H₂)_N-coronene clusters.

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

The coating of various fullerenes by para-hydrogen and ortho-deuterium molecules has been computationally studied as a function of the solvent amount. Rotationally averaged interaction potentials for structureless hydrogen molecules are employed to model their interaction with neutral or charged carbonaceous dopants containing between 20 and 240 atoms, occasionally comparing different fullerenes having the same size but different shapes. The solvation energy and the size of the first solvation shell obtained from path-integral molecular dynamics simulations at 2 K show only minor influence on the dopant charge and on the possible deuteration of the solvent, although the shell size is largest for ortho-D2 coating cationic fullerenes. Nontrivial finite size effects have been found with the shell size varying non-monotonically close to its completion limit. For fullerenes embedded in large hydrogen clusters, the shell size and solvation energy both follow linear scaling with the fullerene size. The shell sizes obtained for C60 (+) and C70 (+) are close to 49 and 51, respectively, and agree with mass spectrometry experiments.

Building Carbon Bridges on and between Fullerenes in Helium Nanodroplets

We report the observation of sequential encounters of fullerenes with C atoms in an extremely cold environment. Experiments were performed with helium droplets at 0.37 K doped with C60 molecules and C atoms derived from a novel, pure source of C atoms. Very high-resolution mass spectra revealed the formation of carbenes of the type C60(C:)n with n up to 6. Bridge-type bonding of the C adatoms to form the known dumbbell C60═C═C60 also was observed. Density functional theory calculations were performed that elucidated the carbene character of the C60(C:)n species and their structures. Mass spectra taken in the presence of water impurities and in separate experiments with added H2 also revealed the formation of the adducts C60C(n)(H2O)n and C60C(n)(H2)n probably by H-OH and H-H bond insertion, respectively, and nonreactivity for the dumbell. So C adatoms that form carbenes C60(C:)n can endow pristine C60 with a higher chemical reactivity.

Experimental evidence for the influence of charge on the adsorption capacity of carbon dioxide on charged fullerenes

The adsorption of CO₂ is sensitive to charge on a capturing model carbonaceous surface, such as C₆₀ fullerenes.

Theoretical study of the interaction between molecular hydrogen and [MC60]+ complexes

In this work we present a density functional theory study of the interaction between a positively charged exohedral metallofullerene and several hydrogen molecules.

A QM:MM model for the interaction of DNA nucleotides with carbon nanotubes

Hybrid materials formed by DNA and carbon nanotubes (CNTs) have shown very interesting properties, but their simulation in solution using quantum mechanical approaches is still a challenge in the computational chemistry community. In this paper, we developed a QM:MM model to study the interactions between charged DNA nucleotides and carbon nanotubes in solution. All four types of DNA nucleotides were taken to interact with two CNTs of similar diameter but different chiralities: (4,4) and (7,0). The nucleotides and CNTs were treated at the QM level, while added water and neutralizing ions were modeled at the MM level. ONIOM simulations were performed at the (M06-2X/6-31G(d):Amber) level for the hybrids, as well as for individually solvated CNTs and nucleotides, which allowed us to evaluate the energy of binding. Our binding energy (BE) values range from 146.60 to 503.43 kJ mol(-1), indicating strong physisorption of nucleotides on CNTs. The relatively large BE, compared with past studies on nucleobase-CNT binding in a vacuum, could be due to the larger size of nucleotides compared with nucleobases, the charges on the nucleotides, and the inclusion of solution which causes the release of water molecules upon hybridization.

Extracting cluster distributions from mass spectra: IsotopeFit

The availability of high resolution mass spectrometry in the study of atomic and molecular clusters opens up challenges for the interpretation of the data. In complex systems each resolved mass peak may contain contributions from multiple species because of the isotope structure of constituent elements and because a multitude of different types of clusters with different compositions are present. A computational procedure which can help to identify a specific cluster from this complex dataset and quantify its relative abundance would be extremely helpful to many who work in this field. Here some new software designed for this purpose, known as IsotopeFit, is described.

On the size and structure of helium snowballs formed around charged atoms and clusters of noble gases

Helium nanodroplets doped with argon, krypton, or xenon are ionized by electrons and analyzed in a mass spectrometer. HenNgx(+) ions containing up to seven noble gas (Ng) atoms and dozens of helium atoms are identified; the high resolution of the mass spectrometer combined with advanced data analysis make it possible to unscramble contributions from isotopologues that have the same nominal mass but different numbers of helium or Ng atoms, such as the magic He20(84)Kr2(+) and the isobaric, nonmagic He41(84)Kr(+). Anomalies in these ion abundances reveal particularly stable ions; several intriguing patterns emerge. Perhaps most astounding are the results for HenAr(+), which show evidence for three distinct, solid-like solvation shells containing 12, 20, and 12 helium atoms. This observation runs counter to the common notion that only the first solvation shell is solid-like but agrees with calculations by Galli et al. for HenNa(+) [J. Phys. Chem. A 2011, 115, 7300] that reveal three shells of icosahedral symmetry. HenArx(+) (2 ≤ x ≤ 7) ions appear to be especially stable if they contain a total of n + x = 19 atoms. A sequence of anomalies in the abundance distribution of HenKrx(+) suggests that rings of six helium atoms are inserted into the solvation shell each time a krypton atom is added to the ionic core, from Kr(+) to Kr3(+). Previously reported strong anomalies at He12Kr2(+) and He12Kr3(+) [Kim , J. H.; et al. J. Chem. Phys. 2006, 124, 214301] are attributed to a contamination. Only minor local anomalies appear in the distributions of HenXex(+) (x ≤ 3). The distributions of HenKr(+) and HenXe(+) show strikingly similar, broad features that are absent from the distribution of HenAr(+); differences are tentatively ascribed to the very different fragmentation dynamics of these ions.

Doubly charged CO2 clusters formed by ionization of doped helium nanodroplets

Helium nanodroplets are doped with carbon dioxide and ionized by electrons. Doubly charged cluster ions are, for the first time, identified based on their characteristic patterns of isotopologues. Thanks to the high mass resolution, large dynamic range, and a novel method to eliminate contributions from singly charged ions from the mass spectra, we are able to observe doubly charged cluster ions that are smaller than the ones reported in the past. The likely mechanism by which doubly charged ions are formed in doped helium droplets is discussed.

4He adsorption on a H(2)-plated C20 molecular surface: the formation of helium buckyballs

We perform path-integral Monte Carlo calculations to study the adsorption of 4He atoms on a H2-plated C20 molecular surface. It is found that 32 H2 molecules form a complete solid layer on C20, where each H2 molecule is located either above one of the 12 pentagon centers or above one of the 20 carbon atoms. The angular density profiles of the first 4He layer on the (H2)32-C20 surface reveal different quantum states as the number of 4He atoms N varies. Especially, the helium layer exhibits an icosidodecahedron structure for N=30, where each 4He atom is located at one of the vertices of 20 corner-sharing triangles. While the 4He density peaks for N=60 constitute a truncated icosahedron with 12 pentagonal and 20 hexagonal faces, the additional atoms beyond N=60 are found to be placed at the hexagon centers of the truncated icosahedron to form a hexakis truncated icosahedron for N=80. The superfluid response of the 4He layer at a temperature of T=0.6 K is found to be completely quenched for N=30 and to be significantly suppressed for N=60 and 80, reflecting the formation of compact buckyball structures.

Ordered phases of ethylene adsorbed on charged fullerenes and their aggregates

In spite of extensive investigations of ethylene adsorbed on graphite, bundles of nanotubes, and crystals of fullerenes, little is known about the existence of commensurate phases; they have escaped detection in almost all previous work. Here we present a combined experimental and theoretical study of ethylene adsorbed on free C₆₀ and its aggregates. The ion yield of [Formula: see text] measured by mass spectrometry reveals a propensity to form a structurally ordered phase on monomers, dimers and trimers of C₆₀ in which all sterically accessible hollow sites over carbon rings are occupied. Presumably the enhancement of the corrugation by the curvature of the fullerene surface favors this phase which is akin to a hypothetical 1 × 1 phase on graphite. Experimental data also reveal the number of molecules in groove sites of the C₆₀ dimer through tetramer. The identity of the sites, adsorption energies and orientations of the adsorbed molecules are determined by molecular dynamics calculations based on quantum chemical potentials, as well as density functional theory. The decrease in orientational order with increasing temperature is also explored in the simulations whereas in the experiment it is impossible to vary the temperature.

Monocarbon cationic cluster yields from N2/CH4 mixtures embedded in He nanodroplets and their calculated binding energies

The formation of monocarbon cluster ions has been investigated by electron ionization mass spectrometry of cold helium nanodroplets doped with nitrogen/methane mixtures. Ion yields for two groups of clusters, CHmN2(+) or CHmN4(+), were determined for mixtures with different molecular ratios of CH4. The possible geometrical structures of these clusters were analyzed using electronic structure computations. Little correlation between the ion yields and the associated binding energies has been observed indicating that in most cases kinetic control is more important than thermodynamic control for forming the clusters.

Decorating (C60) n+, n = 1-3, with CO2 at low temperatures: Sterically enhanced physisorption

Multiple attachment of CO₂ to the monomer, dimer and trimer cations of C₆₀ has been observed in the mass spectra of He nanodroplets sequentially doped with C₆₀ and CO₂ and exposed to electron ionization at 50 eV. Remarkable anomalies were seen in the ion yield for CO₂ coverage for (C₆₀)₂⁺(CO₂)₈ and (C₆₀)₃⁺(CO₂)_1,2. These provide insight into the influence of steric properties on the nature of physisorption. The enhanced stabilities of (C₆₀)₂⁺(CO₂)₈ and (C₆₀)₃⁺(CO₂)_1,2 are attributed to physisorption inside the "groove" of the dimer and the two "dimples" in the trimer cations of C₆₀. Molecular dynamics simulations provide a qualitative assessment of the observed physisorption and a useful visualization of structural aspects.

On the stability of cationic complexes of neon with helium--solving an experimental discrepancy

Helium nanodroplets are doped with neon and ionized by electrons. The size-dependence of the ion abundance of HenNex(+), identified in high-resolution mass spectra, is deduced for complexes containing up to seven neon atoms and dozens of helium atoms. Particularly stable ions are inferred from anomalies in the abundance distributions. Two pronounced anomalies at n = 11 and 13 in the HenNe(+) series confirm drift-tube data reported by Kojima et al. [T. M. Kojima et al., Z. Phys. D, 1992, 22, 645]. The discrepancy with previously published spectra of neon-doped helium droplets, which did not reveal any abundance anomalies [T. Ruchti et al., J. Chem. Phys., 1998, 109, 10679-10687; C. A. Brindle et al., J. Chem. Phys., 2005, 123, 064312], is most likely due to limited mass resolution, which precluded unambiguous analysis of contributions from different ions with identical nominal mass. However, calculated dissociation energies of HenNe(+) reported so far do not correlate with the present data, possibly because of challenges in correctly treating the linear, asymmetric [He-Ne-He](+) ionic core in HenNe(+). Anomalies identified in the distributions of HenNex(+) for x > 1, including prominent ones at He12Ne2(+) and He14Ne2(+), may help to better understand solvation of Ne(+) and Nex(+) in helium.

Adsorption of Polar and Nonpolar Molecules on Isolated Cationic C60 , C70 , and Their Aggregates

Physisorption on graphite, graphene, nanotubes, and other graphitic structures has been the subject of numerous studies, partly driven by interest in the nature of order in two-dimensional systems, their phase transitions, and the use of graphitic scaffolds for reversible storage of hydrogen at high volumetric density and low mass. In contrast, physisorption on individual fullerenes or small aggregates of fullerenes has remained largely unexplored, last but not least, because of technical challenges. A summary of recent progress in identifying specific adsorption sites on positively charged C₆₀ , C₇₀ , and their aggregates is given in this Minireview. Adsorption energies and storage capacities for helium, hydrogen, methane, oxygen, nitrogen, water, and ammonia are determined. Mass spectrometric data reveal the formation of a commensurate phase in which all hollow sites of C₆₀ or C₇₀ are occupied. This phase is identified for all nonpolar molecules, including oxygen, which does not form a commensurate phase on planar graphite. The polar molecules, on the other hand, do not wet fullerenes and they do not form this commensurate phase. A hierarchy of other distinct adsorption sites are identified for nonpolar molecules, namely, groove sites for fullerene dimers and beyond, and dimple sites for fullerene trimers and beyond. Furthermore, evidence is presented for the preferential adsorption of hydrogen and methane in registered sites on fullerene dimers. The interpretation of experimental data that merely count the number of preferred adsorption sites is aided by molecular dynamics simulations, which utilize interaction potentials derived from ab initio calculations to determine adsorption energies.

Methane adsorption on aggregates of fullerenes: site-selective storage capacities and adsorption energies

Methane adsorption on positively charged aggregates of C60 is investigated by both mass spectrometry and computer simulations. Calculated adsorption energies of 118-281 meV are in the optimal range for high-density storage of natural gas. Groove sites, dimple sites, and the first complete adsorption shells are identified experimentally and confirmed by molecular dynamics simulations, using a newly developed force field for methane-methane and fullerene-methane interaction. The effects of corrugation and curvature are discussed and compared with data for adsorption on graphite, graphene, and carbon nanotubes.