Invited review article: laser vaporization cluster sources

The harpooning mechanism as evidenced in the oxidation reaction of the Al atom

The harpooning model is firstly evidenced by the maximum impact parameter derived from AlO(N_max) products with speed v′ ≈ 0 from the Al + O₂ reaction.

The harpooning mechanism has long been proposed for elementary reaction dynamics involving metals. It is characterized by an initial electron transfer (ET) process from the metal to the oxidant molecule. For the titled reaction Al + O₂, the ET distance can be predicted to be 2.6 Å by simply calculating the energy difference between the ionization energy of the Al atom and the electron affinity of the O₂ molecule. Hereby we experimentally derived the maximum impact parameter b_max of 2.5 ± 0.2 Å for the titled reaction, in consistency with the predicted ET distance. This derivation of b_max was achieved by using the crossed molecular beam experiment at a collision energy of 507 cm^–1 (i.e. 1.45 kcal mol^–1) with a high resolution time-sliced ion velocity imaging detection of the state-selective AlO products based on the (1 + 1) resonance-enhanced multiphoton ionization. The small rotational constant of the AlO(X²Σ⁺) radical (B_e = 0.6413 cm^–1) facilitated the formation of the AlO(v = 0) products in high rotational levels up to the energetically limited state, N_max = 52, with an almost zero velocity mapping. Hence, in this extreme angular momentum disposal case, the collisional orbital angular momentum l was nearly completely channeled into the product rotational angular momentum as a consequence of the conservations of energy and angular momentum, offering a reaction system that breaks the restriction of kinematically favored mass combination in order to obtain information on the impact parameters. The present study yields the first direct derivation of b_max from the maximum rotational level of products under the experimental condition with the recoil energy E′_T ≈ 0. This, in turn, provides solid evidence in supporting the harpooning mechanism.

Stability of small cationic platinum clusters

The relative stability of small cationic platinum clusters is investigated by photofragmentation experiments and density functional theory calculations.

Improved tandem mass spectrometer coupled to a laser vaporization cluster ion source

We describe two improvements to an existing tandem mass spectrometer coupled to a laser vaporization cluster ion source suitable for photodissociation spectroscopy: (i) cooling of the cluster source nozzle and (ii) mass selection prior to the photodissociation region via replacing an octupole ion guide by a quadrupole mass spectrometer. The improved sensitivity and transmission enable the production of larger heteroatomic clusters as well as rare gas solvated clusters. We present two examples demonstrating the new capabilities of the improved setup. In the first application, cooling of the cluster source nozzle produces Si⁺Ar_n and Si₂⁺Ar_n cluster cations with n = 1-25. Magic numbers are extracted from the mass spectrum by applying a transmission function obtained via simulations. In the second example, the vibronic photodissociation spectrum of cold Au₄⁺ cluster ions is recorded with unprecedented detail, resolution, and sensitivity. Such high-resolution optical excitation spectra of metal cluster cations may serve as a benchmark for the performance of Franck-Condon simulations based on quantum chemical calculations for excited states.

Transformation of doped graphite into cluster-encapsulated fullerene cages

An ultimate goal in carbon nanoscience is to decipher formation mechanisms of highly ordered systems. Here, we disclose chemical processes that result in formation of high-symmetry clusterfullerenes, which attract interest for use in applications that span biomedicine to molecular electronics. The conversion of doped graphite into a C₈₀ cage is shown to occur through bottom-up self-assembly reactions. Unlike conventional forms of fullerene, the iconic Buckminsterfullerene cage, I_h-C₆₀, is entirely avoided in the bottom-up formation mechanism to afford synthesis of group 3-based metallic nitride clusterfullerenes. The effects of structural motifs and cluster–cage interactions on formation of compounds in the solvent-extractable C₇₀–C₁₀₀ region are determined by in situ studies of defined clusterfullerenes under typical synthetic conditions. This work establishes the molecular origin and mechanism that underlie formation of unique carbon cage materials, which may be used as a benchmark to guide future nanocarbon explorations.

An understanding of how caged carbon materials self-assemble from doped graphite is a long-standing challenge. Here, the authors show that distinct bottom-up processes lead to the synthesis of high-symmetry clusterfullerenes.

Photoelectron Velocity Map Imaging Spectroscopy of the Beryllium Sulfide Anion, BeS. J Phys Chem A 2017;121:5645-5650. [PMID: 28691819 DOI: 10.1021/acs.jpca.7b04894]

Slow electron velocity map imaging (SEVI) spectroscopy was used to examine the BeS^- anion to neutral ground-state transition, X ²Σ⁺ → X ¹Σ⁺. Rotational constants, vibrational intervals, and the electron binding energy of BeS^- were determined. Partially resolved rotational contours were seen due to the relatively small moment of inertia of beryllium sulfide. Upon analysis of the rotational contours, it was found that changes in the molecular rotational angular momentum, ΔN = -1, -2, -3, and -4, facilitated photodetachment at near-threshold photon energies. The electron affinity of BeS was found to be 2.3346(2) eV. SEVI spectra recorded using photon energies near the threshold for Δv = -1 processes exhibited features that were associated with a dipole-bound state (DBS) of BeS^-. Autodetachment spectroscopy was used to probe this state, and rotationally resolved data were obtained for the DBS ²Σ⁺, v' = 0 - X ²Σ⁺, v″ = 0 transition. Analysis of this structure provided the rotational constants for BeS^- X, v″ = 0, and the electron binding energy of the DBS. Electronic structure calculations, performed at the RCCSD(T) and MRCI levels of theory, gave predictions that were in good agreement with the experimental observations.

A gas aggregation source for the production of heterogeneous molecular clusters

We present the design of a versatile gas aggregation source that allows producing molecular beams of charged clusters containing a controlled amount of chosen impurities. Several examples of clusters production using this source characterized by time of flight mass spectrometry are presented here. We demonstrate the source ability to produce homogeneous clusters, such as pure protonated water and alcohol clusters, as well as inhomogeneous ones such as water clusters containing a few units of uracil, glycine, sulfuric acid, or pyrene.

Photodetachment spectroscopy of the beryllium oxide anion, BeO. J Chem Phys 2017;146:054301. [PMID: 28178838 DOI: 10.1063/1.4974843]

The X²Σ⁺→X¹Σ⁺ anion to neutral ground state photodetachment of BeO^- has been studied by means of photoelectron velocity-map imaging spectroscopy in a newly constructed apparatus. Vibrational intervals, rotational constants, and the electron detachment threshold of BeO^- were determined for the first time. The small moment of inertia of beryllium oxide allowed for the observation of partially resolved rotational contours. Analyses of these contours provided evidence of several detachment channels resulting from changes in molecular rotational angular momenta of ΔN = 0, ±1, ±2, and ±3. The relative intensities of these detachment channels were found to be a function of the electron kinetic energy. Experimental results are compared to the predictions of high level ab initio calculations.

Probing Ionic Complexes of Ethylene and Acetylene with Vacuum-Ultraviolet Radiation

Mixed complexes of acetylene-ethylene are studied using vacuum-ultraviolet (VUV) photoionization mass spectrometry and theoretical calculations. These complexes are produced and ionized at different distances from the exit of a continuous nozzle followed by reflectron time-of-flight mass spectrometry detection. Acetylene, with a higher ionization energy (11.4 eV) than ethylene (10.6 eV), allows for tuning the VUV energy and initializing reactions either from a C2H2(+) or a C2H4(+) cation. Pure acetylene and ethylene expansions are separately carried out to compare, contrast, and hence identify products from the mixed expansion: these are C3H3(+) (m/z = 39), C4H5(+) (m/z = 53), and C5H5(+) (m/z = 65). Intensity distributions of C2H2, C2H4, their dimers and reactions products are plotted as a function of ionization distance. These distributions suggest that association mechanisms play a crucial role in product formation closer to the nozzle. Photoionization efficiency (PIE) curves of the mixed complexes demonstrate rising edges closer to both ethylene and acetylene ionization energies. We use density functional theory (ωB97X-V/aug-cc-pVTZ) to study the structures of the neutral and ionized dimers, calculate their adiabatic and vertical ionization energies, as well as the energetics of different isomers on the potential energy surface (PES). Upon ionization, vibrationally excited clusters can use the extra energy to access different isomers on the PES. At farther ionization distances from the nozzle, where the number densities are lower, unimolecular decay is expected to be the dominant mechanism. We discuss the possible decay pathways from the different isomers on the PES and examine the ones that are energetically accessible.

Recent Advances in the Synthesis, Characterization and Application of Zn+-containing Heterogeneous Catalysts

Monovalent Zn+ (3d104s1) systems possess a special electronic structure that can be exploited in heterogeneous catalysis and photocatalysis, though it remains challenge to synthesize Zn+-containing materials. By careful design, Zn+-related species can be synthesized in zeolite and layered double hydroxide systems, which in turn exhibit excellent catalytic potential in methane, CO and CO2 activation. Furthermore, by utilizing advanced characterization tools, including electron spin resonance, X-ray absorption fine structure and density functional theory calculations, the formation mechanism of the Zn+ species and their structure-performance relationships can be understood. Such advanced characterization tools guide the rational design of high-performance Zn+-containing catalysts for efficient energy conversion.

Understanding ligand effects in gold clusters using mass spectrometry

This review summarizes recent research on the influence of phosphine ligands on the size, stability, and reactivity of gold clusters synthesized in solution. Sub-nanometer clusters exhibit size- and composition-dependent properties that are unique from those of larger nanoparticles. The highly tunable properties of clusters and their high surface-to-volume ratio make them promising candidates for a variety of technological applications. However, because "each-atom-counts" toward defining cluster properties it is critically important to develop robust synthesis methods to efficiently prepare clusters of predetermined size. For decades phosphines have been known to direct the size-selected synthesis of gold clusters. Despite the preparation of numerous species it is still not understood how different functional groups at phosphine centers affect the size and properties of gold clusters. Using electrospray ionization mass spectrometry (ESI-MS) it is possible to characterize the effect of ligand substitution on the distribution of clusters formed in solution at defined reaction conditions. In addition, ligand exchange reactions on preformed clusters may be monitored using ESI-MS. Collision induced dissociation (CID) may also be employed to obtain qualitative insight into the fragmentation of mixed ligand clusters and the relative binding energies of differently substituted phosphines. Quantitative ligand binding energies and cluster stability may be determined employing surface induced dissociation (SID) in a custom-built Fourier transform ion cyclotron resonance mass spectrometer (FT-ICR-MS). Rice-Ramsperger-Kassel-Marcus (RRKM) based modeling of the SID data allows dissociation energies and entropy values to be extracted. The charge reduction and reactivity of atomically precise gold clusters, including partially ligated species generated in the gas-phase by in source CID, on well-defined surfaces may be explored using ion soft landing (SL) in a custom-built instrument combined with in situ time of flight secondary ion mass spectrometry (TOF-SIMS). Jointly, this multipronged experimental approach allows characterization of the full spectrum of relevant phenomena including cluster synthesis, ligand exchange, thermochemistry, surface immobilization, and reactivity. The fundamental insights obtained from this work will facilitate the directed synthesis of gold clusters with predetermined size and properties for specific applications.

Soft- and reactive landing of ions onto surfaces: Concepts and applications

Soft- and reactive landing of mass-selected ions is gaining attention as a promising approach for the precisely-controlled preparation of materials on surfaces that are not amenable to deposition using conventional methods. A broad range of ionization sources and mass filters are available that make ion soft-landing a versatile tool for surface modification using beams of hyperthermal (<100 eV) ions. The ability to select the mass-to-charge ratio of the ion, its kinetic energy and charge state, along with precise control of the size, shape, and position of the ion beam on the deposition target distinguishes ion soft landing from other surface modification techniques. Soft- and reactive landing have been used to prepare interfaces for practical applications as well as precisely-defined model surfaces for fundamental investigations in chemistry, physics, and materials science. For instance, soft- and reactive landing have been applied to study the surface chemistry of ions isolated in the gas-phase, prepare arrays of proteins for high-throughput biological screening, produce novel carbon-based and polymer materials, enrich the secondary structure of peptides and the chirality of organic molecules, immobilize electrochemically-active proteins and organometallics on electrodes, create thin films of complex molecules, and immobilize catalytically active organometallics as well as ligated metal clusters. In addition, soft landing has enabled investigation of the size-dependent behavior of bare metal clusters in the critical subnanometer size regime where chemical and physical properties do not scale predictably with size. The morphology, aggregation, and immobilization of larger bare metal nanoparticles, which are directly relevant to the design of catalysts as well as improved memory and electronic devices, have also been studied using ion soft landing. This review article begins in section 1 with a brief introduction to the existing applications of ion soft- and reactive landing. Section 2 provides an overview of the ionization sources and mass filters that have been used to date for soft landing of mass-selected ions. A discussion of the competing processes that occur during ion deposition as well as the types of ions and surfaces that have been investigated follows in section 3. Section 4 discusses the physical phenomena that occur during and after ion soft landing, including retention and reduction of ionic charge along with factors that impact the efficiency of ion deposition. The influence of soft landing on the secondary structure and biological activity of complex ions is addressed in section 5. Lastly, an overview of the structure and mobility as well as the catalytic, optical, magnetic, and redox properties of bare ionic clusters and nanoparticles deposited onto surfaces is presented in section 6.

Photodissociation of Cerium Oxide Nanocluster Cations

Cerium oxide cluster cations, CexOy(+), are produced via laser vaporization in a pulsed nozzle source and detected with time-of-flight mass spectrometry. The mass spectrum displays a strongly preferred oxide stoichiometry for each cluster with a specific number of metal atoms x, with x ≤ y. Specifically, the most prominent clusters correspond to the formula CeO(CeO2)n(+). The cluster cations are mass selected and photodissociated with a Nd:YAG laser at either 532 or 355 nm. The prominent clusters dissociate to produce smaller species also having a similar CeO(CeO2)n(+) formula, always with apparent leaving groups of (CeO2). The production of CeO(CeO2)n(+) from the dissociation of many cluster sizes establishes the relative stability of these clusters. Furthermore, the consistent loss of neutral CeO2 shows that the smallest neutral clusters adopt the same oxidation state (IV) as the most common form of bulk cerium oxide. Clusters with higher oxygen content than the CeO(CeO2)n(+) masses are present with much lower abundance. These species dissociate by the loss of O2, leaving surviving clusters with the CeO(CeO2)n(+) formula. Density functional theory calculations on these clusters suggest structures composed of stable CeO(CeO2)n(+) cores with excess oxygen bound to the surface as a superoxide unit (O2(-)).

Experimental and theoretical identification of the Fe(vii) oxidation state in FeO4−

Two isomers of iron tetraoxygen anion, dioxoiron peroxide [(η²-O₂)FeO₂]⁻ and tetroxide FeO₄⁻ were characterized by experiment and theoretical calculations, with heptavalent Fe(vii) oxidation state identified in the later.

Binding Hydrated Anions with Hydrophobic Pockets

Using a combination of isothermal titration calorimetry and quantum and molecular dynamics calculations, we demonstrate that relatively soft anions have an affinity for hydrophobic concavity. The results are consistent with the anions remaining partially hydrated upon binding, and suggest a novel strategy for anion recognition.

The size-dependent morphology of Pd nanoclusters formed by gas condensation

Size-selected Pd nanoclusters in the size range from 887 to 10,000 atoms were synthesized in a magnetron sputtering, inert gas condensation cluster beam source equipped with a time of flight mass filter. Their morphologies were investigated using scanning transmission electron microscopy (STEM) and shown to be strongly size-dependent. The larger clusters exhibited elongated structures, which we attribute to the aggregation, through multiple collisions, of smaller clusters during the gas phase condensation process. This was confirmed from the atomically resolved STEM images of the Pd nanoclusters, which showed smaller primary clusters with their own crystalline structures.

Imaging charge transfer in a cation-π system: velocity-map imaging of Ag(+)(benzene) photodissociation

Ag(+)(benzene) complexes are generated in the gas phase by laser vaporization and mass selected in a time-of-flight spectrometer. UV laser excitation at either 355 or 266 nm results in dissociative charge transfer (DCT), leading to neutral silver atom and benzene cation products. Kinetic energy release in translationally hot benzene cations is detected using a new instrument designed for photofragment imaging of mass-selected ions. Velocity-map imaging and slice imaging techniques are employed. In addition to the expected translational energy release, DCT of Ag(+)(benzene) produces a distribution of internally hot benzene cations. Compared with experiments at 355 nm, 266 nm excitation produces only slightly higher translational excitation and a much greater fraction of internally hot benzene ions. The maximum kinetic energy release in the photodissociation sets an upper limit on the Ag(+)(benzene) dissociation energy of 32.8 (+1.4/-1.5) kcal/mol.

Direct experimental observation of weakly-bound character of the attached electron in europium anion

Direct experimental determination of precise electron affinities (EAs) of lanthanides is a longstanding challenge to experimentalists. Considerable debate exists in previous experiment and theory, hindering the complete understanding about the properties of the atomic anions. Herein, we report the first precise photoelectron imaging spectroscopy of europium (Eu), with the aim of eliminating prior contradictions. The measured EA (0.116 ± 0.013 eV) of Eu is in excellent agreement with recently reported theoretical predictions, providing direct spectroscopic evidence that the additional electron is weakly attached. Additionally, a new experimental strategy is proposed that can significantly increase the yield of the lanthanide anions, opening up the best opportunity to complete the periodic table of the atomic anions. The present findings not only serve to resolve previous discrepancy but also will help in improving the depth and accuracy of our understanding about the fundamental properties of the atomic anions.

Design and performance of a high-flux electrospray ionization source for ion soft landing

A high-flux electrospray source enables deposition of micrograms of mass-selected ions for studies in catalysis and materials science.

Bottom-up formation of endohedral mono-metallofullerenes is directed by charge transfer

An understanding of chemical formation mechanisms is essential to achieve effective yields and targeted products. One of the most challenging endeavors is synthesis of molecular nanocarbon. Endohedral metallofullerenes are of particular interest because of their unique properties that offer promise in a variety of applications. Nevertheless, the mechanism of formation from metal-doped graphite has largely eluded experimental study, because harsh synthetic methods are required to obtain them. Here we report bottom-up formation of mono-metallofullerenes under core synthesis conditions. Charge transfer is a principal factor that guides formation, discovered by study of metallofullerene formation with virtually all available elements of the periodic table. These results could enable production strategies that overcome long-standing problems that hinder current and future applications of metallofullerenes.

Ligand-free Ni nanocluster formation at atmospheric pressure via rapid quenching in a microplasma process

The production of metal nanoclusters composed of less than 10(3) atoms is important for applications in energy conversion and medicine, and for fundamental studies of nanomaterial nucleation and growth. Unfortunately, existing synthesis methods do not enable adequate control of cluster formation, particularly at atmospheric pressure wherein formation typically occurs on sub-millisecond timescales. Here, we demonstrate that ligand-free, unagglomerated nickel nanoclusters can be continuously synthesized at atmospheric pressure via the decomposition of bis(cyclopentadienyl)nickel(II) (nickelocene) in a spatially-confined microplasma process that rapidly quenches particle growth and agglomeration. The clusters were measured on line by ion mobility spectrometry (IMS) and further analyzed by atomic force microscopy (AFM). Our results reveal that stable clusters with spherical equivalent mean diameters below 10 Åare produced, and by controlling the nickelocene concentration, the mean diameter can be tuned up to ∼50 Å. Although diameter is often the sole metric used in nanocluster and nanoparticle characterization, to infer the number of atoms in AFM and IMS detected clusters, we compare measured AFM heights and IMS inferred collision cross sections to theoretical predictions based on both bulk matter approximations and density functional theory and Hartree-Fock calculated Ni nanocluster structures (composed of 2-15 atoms for the latter). The calculations suggest that Ni nanoclusters composed of less than 10(2) atoms can be produced repeatably with simple microplasma reactors.

Helium droplet calorimetry of strongly bound species: carbon clusters from C₂ to C₁₂

Helium droplet beam methods are a versatile technique that can be used to assemble a wide variety of atomic and molecular clusters. In recent years, methods have been developed to utilize helium droplets as nano-calorimeters to measure the binding energies of weakly bound complexes assembled within the droplet. In the current investigation we extend the helium droplet calorimetry approach to the study of a very strongly bound system: carbon clusters which are bound by several eV per atom. We utilize laser heating of bulk carbon samples to dope the helium droplets with evaporated carbon species. Depending on the laser target, the vaporization plume is found to consist primarily of C3 alone or C2 and C3. These species are sequentially captured by the droplet and assembled into larger carbon clusters in a stepwise manner. The assembled C(n) clusters are detected via mass spectrometry of the doped droplets and the droplet sizes required to detect the various carbon clusters observed are used to estimate the reaction energies of the associated assembly pathways. The helium droplet data qualitatively reflect the trends in assembly energetics, but at first glance appear to yield energies that differ dramatically from theoretical values. Statistical modeling of the helium droplet calorimetry experiment reconciles the differences quantitatively. Our modeling also generates a calibration curve that relates the assembly/reaction energy and threshold mean droplet size over a range of energies from van der Waals interactions to chemical bonding, enabling helium droplet calorimetry methods to be applied quantitatively to a large number of systems.

In situ SIMS and IR spectroscopy of well-defined surfaces prepared by soft landing of mass-selected ions

Soft landing of mass-selected ions onto surfaces is a powerful approach for the highly-controlled preparation of materials that are inaccessible using conventional synthesis techniques. Coupling soft landing with in situ characterization using secondary ion mass spectrometry (SIMS) and infrared reflection absorption spectroscopy (IRRAS) enables analysis of well-defined surfaces under clean vacuum conditions. The capabilities of three soft-landing instruments constructed in our laboratory are illustrated for the representative system of surface-bound organometallics prepared by soft landing of mass-selected ruthenium tris(bipyridine) dications, [Ru(bpy)3](2+) (bpy = bipyridine), onto carboxylic acid terminated self-assembled monolayer surfaces on gold (COOH-SAMs). In situ time-of-flight (TOF)-SIMS provides insight into the reactivity of the soft-landed ions. In addition, the kinetics of charge reduction, neutralization and desorption occurring on the COOH-SAM both during and after ion soft landing are studied using in situ Fourier transform ion cyclotron resonance (FT-ICR)-SIMS measurements. In situ IRRAS experiments provide insight into how the structure of organic ligands surrounding metal centers is perturbed through immobilization of organometallic ions on COOH-SAM surfaces by soft landing. Collectively, the three instruments provide complementary information about the chemical composition, reactivity and structure of well-defined species supported on surfaces.

Laboratory spectroscopy of astrophysically relevant carbon species

Carbon is one of the most common elements in the solar system, with a fractional abundance of 10(-4) relative to hydrogen. Thus, it is not surprising that over 100 carbon-bearing species have been definitively detected in the interstellar medium via their rotational, infrared, and/or electronic transitions. In order to identify these species, laboratory spectra are needed for comparison to astronomical data. Challenges arise when obtaining laboratory spectra due to the instability of many of these molecules. Over the years, sensitive instrumentation and better techniques for producing these species in situ have been developed to achieve this goal. The use of complementary spectroscopic methods, such as matrix isolation, cavity ringdown, resonance enhanced multiphoton ionization, and ion trapping have led to the identification of several new carbon species at optical and ultraviolet wavelengths. Laboratory spectra have been compared to astronomical data in order to gain further insight into interstellar chemistry. In particular, attempts have been made to identify the carriers of the diffuse interstellar bands, however, with little success. These results are discussed in the following review.

IR spectroscopy of gas phase V(CO2)n+ clusters: solvation-induced electron transfer and activation of CO2

Ion-molecule complexes of vanadium and CO2, i.e., V(CO2)n(+), produced by laser vaporization are mass selected and studied with infrared laser photodissociation spectroscopy. Vibrational bands for the smaller clusters (n < 7) are consistent with CO2 ligands bound to the metal cation via electrostatic interactions and/or attaching as inert species in the second coordination sphere. All IR bands for these complexes are consistent with intact CO2 molecules weakly perturbed by cation binding. However, multiple new IR bands occur only in larger complexes (n ≥ 7), indicating the formation of an intracluster reaction product whose nominal mass is the same as that of V(CO2)n(+) complexes. Computational studies and the comparison of predicted spectra for different possible reaction products allow identification of an oxalate-type C2O4 anion species in the cluster. The activation of CO2 producing this product occurs via a solvation-induced metal→ligand electron transfer reaction.

Metallofullerene and fullerene formation from condensing carbon gas under conditions of stellar outflows and implication to stardust

Carbonaceous presolar grains of supernovae origin have long been isolated and are determined to be the carrier of anomalous (22)Ne in ancient meteorites. That exotic (22)Ne is, in fact, the decay isotope of relatively short-lived (22)Na formed by explosive nucleosynthesis, and therefore, a selective and rapid Na physical trapping mechanism must take place during carbon condensation in supernova ejecta. Elucidation of the processes that trap Na and produce large carbon molecules should yield insight into carbon stardust enrichment and formation. Herein, we demonstrate that Na effectively nucleates formation of Na@C60 and other metallofullerenes during carbon condensation under highly energetic conditions in oxygen- and hydrogen-rich environments. Thus, fundamental carbon chemistry that leads to trapping of Na is revealed, and should be directly applicable to gas-phase chemistry involving stellar environments, such as supernova ejecta. The results indicate that, in addition to empty fullerenes, metallofullerenes should be constituents of stellar/circumstellar and interstellar space. In addition, gas-phase reactions of fullerenes with polycyclic aromatic hydrocarbons are investigated to probe "build-up" and formation of carbon stardust, and provide insight into fullerene astrochemistry.