A molecular-dynamics simulation study of solvent-induced repulsion between C60 fullerenes in water

Molecular-dynamics simulations of a single C(60) fullerene and pairs of C(60) fullerenes in aqueous solution have been performed for the purpose of obtaining improved understanding of the nature of solvent-induced interactions between C(60) fullerenes in water. Our simulations reveal repulsive solvent-induced interactions between two C(60) fullerenes in aqueous solution in contrast to the associative effects observed for conventional nonpolar solutes. A decomposition of the solvent-induced potential of mean force between fullerenes into entropy and energy (enthalpy) contributions reveals that the water-induced repulsion between fullerenes is energetic in origin, contrasting strongly to entropy-driven association observed for conventional nonpolar solutes. The dominance of energy in the solvent-induced interactions between C(60) fullerenes arises primarily from the high atomic density of the C(60) molecule, resulting in strong C(60)-water van der Waals attraction that is reduced upon association of the fullerenes. The water-induced repulsion is found to decrease with increasing temperature due largely to an increasing contribution from a relatively weak entropy-driven association.     

Supramolecular self-organization in PEO-modified C60 fullerene/water solutions: influence of polymer molecular weight and nanoparticle concentration

Utilizing a first-principles-based coarse-grained implicit solvent model, we have investigated the self-association of C(60) fullerenes that have been symmetrically modified with six grafted poly(ethylene oxide) (PEO) chains in aqueous solution. Despite the highly symmetric nature of the pair interactions between PEO-grafted fullerenes, their supramolecular assemblies are highly anisotropic and resemble the linear clusters formed in Stockmayer fluids. The dipole-like interaction between these symmetrically modified fullerenes results from the shielding of the C(60) fullerenes by PEO, favoring the addition of more PEO-grafted fullerenes to the linear clusters at the relatively unprotected ends. At low nanoparticle concentrations, self-association is dominated by the formation of stable dimers and trimers resulting from fullerene-fullerene contact and favorable PEO-fullerene interactions. With increasing nanoparticle concentration, larger clusters become increasingly probable. The molecular weight of the PEO tethers can be treated as a temperature-like analogue, with a reduction in average cluster size with increasing chain length due to increased steric repulsion, which is qualitatively similar to effects observed in Stockmayer fluids with increasing temperature. The role of PEO in supramolecular self-organization in PEO-modified C(60) fullerene/water solutions is complex, contributing not only to steric stabilization but also to favorable energetic interactions, nanoparticle shielding, and depletion-driven aggregation.     

Selective extraction of higher fullerenes using cyclic dimers of zinc porphyrins

Higher fullerenes (>/=C76) were selectively extracted from a fullerene mixture obtained from a combustion-based industrial production source by cyclic dimers of beta-unsubstituted porphyrin zinc complexes 2C5-2C7 with C5-C7 alkylene spacers as host molecules. Results of single extraction of the fullerene mixture with 2C5-2C7 together with a beta-substituted analogue of 2C6 (1C6) and spectroscopic titration of 2C6 and 1C6 with C60, C70, and C96 indicated that the host selectivity toward higher fullerenes is much dependent on the structure of the porphyrin units and the size of the host cavity. Sequential three-stage extraction of the fullerene mixture with the best-behaved 2C6 resulted in considerable enrichment in very rare fullerenes C102-C110 (<0.1 abs %) up to 82 abs % (C76-C114, 99 abs %) (356 nm) of total fullerenes.     

Stability of silicon-doped C60 dimers

A theoretical investigation on the structure, stability, and thermal behaviors of the smallest polymeric units, the dimers, formed from substitutionally Si-doped fullerenes is presented. A density functional based nonorthogonal tight-binding model has been employed for describing the interatomic interactions. The study focuses on those polymeric structures which involve Si-Si or Si-C interfullerene bonds. The binding energy of the dimers increases with their Si content from about 0.25 eV in C(60)-C(60) to about 4.5 eV in C(58)Si(2)-C(58)Si(2). Moreover, the C(59)SiC(59) dimer, linked through the sharing of the Si atom between the two fullerenes, has been also considered. Upon heating, the dimers eventually fragment into their constituent fullerene units. The fragmentation temperature correlates with the strength of the interfullerene bonds. C(58)Si(2)-C(58)Si(2) exhibits a higher thermal stability (fragmentation temperature of approximately 500 K) than the pure carbon C(60)-C(60) dimer (with a fragmentation temperature of approximately 325 K). Given the higher structural and thermal stabilities of the Si-doped fullerene dimers, the authors propose the use of substitutionally Si-doped fullerenes as the basic units for constructing new fullerene-based polymers.     

Simulating the thermal behavior and fragmentation mechanisms of exohedral and substitutional silicon-doped C60

Structures, thermal behavior, and fragmentation mechanisms of exohedral and substitutional silicon-doped C(60) containing 1-12 Si atoms are investigated by extensive molecular-dynamics simulations. A nonorthogonal tight-binding model is used to mimic the interatomic interactions in the doped fullerenes. Beginning from the minimum-energy structures, the temperature of the doped fullerenes is slowly increased until fragmentation takes place. A correlation can be established between the exohedral and substitutional structures and the corresponding fragmentation mechanisms and fragmentation temperatures. Exohedral C(60)Si(m) fullerenes fragment into two homonuclear pieces, the Si(m) cluster and the C(60) fullerene that remains intact. In contrast, the substitutional C(60-m)Si(m) heterofullerenes undergo structural transformations, including the partial unraveling of the cage, prior to fragmentation. Then, ejection of atoms or small molecules takes place from the distorted structures. The slow heating rate used, combined with long simulation runs, allows us to determine the fragmentation temperature of exohedral and substitutional Si-doped fullerenes as a function of the number of silicon atoms. Substitutional Si-doped fullerenes exhibit much higher fragmentation temperatures (1000-1500 K higher) than the exohedral fullerenes. This can be understood from the different bonding of the Si atoms in both structures.     

Photoabsorption spectra of small fullerenes and Si-heterofullerenes

We study the spectral properties of two kinds of derivatives of the carbon fullerene C(60), small fullerenes and Si-heterofullerenes, by ab initio calculations. The principal method of study is the time-dependent density-functional theory in its full time-propagation form. C(20), C(28), C(32), C(36), and C(50), the most stable small fullerenes in the range of C(20)-C(50), are found to have characteristic features in their optical absorption spectra, originating from the geometry of the molecules in question. The comparison of measured and calculated absorption spectra is found to be a useful tool in differentiating between different, almost isoenergetic ground state structure candidates of small fullerenes. Substitutionally doped fullerenes are of interest due to their enhanced chemical reactivity. It is suggested that the doping degree can be obtained by studying the absorption spectra. For example, it is observed that the spectra gradually change when doping C(60) up to C(48)Si(12) so that absorption in the visible and near infrared regions increases.     

A molecular dynamics simulation study of C60 fullerenes inside a dimyristoylphosphatidylcholine lipid bilayer

We have carried out atomistic molecular dynamics simulations of C60 fullerenes inside a dimyristoylphosphatidylcholine lipid bilayer and an alkane melt. Simulations reveal that the preferred position of a single C60 fullerene is about 6-7 A off of the center plane, allowing the fullerene to take advantage of strong dispersion interactions with denser regions of the bilayer. Further displacement (>8 A) of the fullerene away from the center plane results in a rapid increase in free energy likely due to distortion of the lipid head group layer. The effective interaction between fullerenes (direct interaction plus environment (bilayer)-induced interaction), measured as the potential of mean force (POMF) between two fullerenes as a function of their separation, was found to be significantly less attractive in the lipid bilayer than in an alkane melt of the same molecular weight as the lipid tails. Only part of this difference can be accounted for by the more favorable interaction of the fullerene with the relatively denser bilayer. Additionally, our POMF studies indicate that the bilayer is less able to accommodate the larger aggregated fullerene pair than isolated single fullerenes, again likely due to distortion of the bilayer structure. The implications of these effects on aggregation of fullerenes within lipid bilayer are considered.     

Extraction and high-performance liquid chromatographic analysis of C60, C70, and [6,6]-phenyl C61-butyric acid methyl ester in synthetic and natural waters

Studies have shown that C(60) fullerene can form stable colloidal suspensions in water that result in C(60) aqueous concentrations many orders of magnitude above C(60)'s aqueous solubility; however, quantitative methods for the analysis of C(60) and other fullerenes in environmental media are scarce. Using a 80/20v/v toluene-acetonitrile mobile phase and a 4.6mmx150mm Cosmosil 5mu PYE column, C(60), C(70), and PCBM ([6,6]-phenyl C(61)-butyric acid methyl ester) were fully resolved. Selectivity factors (alpha) for C(60) relative to PCBM and C(70) relative to C(60) were 3.18 and 2.19, respectively. The best analytical wavelengths for the fullerenes were determined to be 330, 333, and 333nm with log molar absorption coefficients (logvarepsilon) of 4.63, 4.82, and 4.60 for PCBM, C(60), C(70), respectively. Extraction and quantitation of all three fullerenes in aqueous suspensions over a range of pH (4-10) and ionic strengths were very good. Whole-method quantification limits for ground and surface suspensions were 2.87, 2.48, and 6.54mug/L for PCBM, C(60), and C(70), respectively.     

Stabilities of multiply charged dimers and clusters of fullerenes

The authors find even-odd variations as functions of r (...+[C60]2(r+)([C60C70](r+)) electron-transfer collisions. This even-odd behavior is in sharp contrast to the smooth one for fullerene monomers and may be related to even-odd effects in dimer ionization energies in agreement with results from an electrostatic model. The kinetic energy releases for dimer dissociations [predominantly yielding intact fullerenes [C60]2(r+)-->C60(r1+)+C60(r2+) in the same (r1=r2) or nearby (r1=r2+/-1) charge states] are found to be low in comparison with the corresponding model results indicating that internal excitations of the separating (intact) fullerenes are important. Experimental appearance sizes for the heavier clusters of fullerenes [C60]n(r+) (n>3 and r=2-5) compare well with predictions from a new nearest-neighbor model assuming that r unit charges in [C60]n(r+) are localized to r C60 molecules such that the Coulomb energy of the system is minimized. The system is then taken to be stable if (i) two (singly) charged C60 are not nearest neighbors and (ii) the r C60(+) molecules have binding energies to their neutral nearest neighbors which are larger than the repulsive energies for the (r-1) C60(+)-C60(+) pairs. Essential ingredients in the nearest-neighbor model are cluster geometries and the present results on dimer stabilities.     

Thermal relaxation mechanism and role of chemical functionalization in fullerene solutions

Using molecular-dynamics simulations we investigate thermal relaxation of C60 and C84 molecules suspended in octane liquid. Pristine fullerenes exhibit relatively slow relaxation due to weak thermal coupling with the liquid. A comparison of the interfacial transport characteristics obtained from relaxation simulations with those obtained from equilibrium simulations and fluctuation-dissipation theorem analysis demonstrates that the relaxation process involves two main steps: (i) energy flow from high- to low-frequency modes within the fullerene, and (ii) energy flow from low-frequency fullerene modes to the liquid. Functionalization of fullerenes with alkene chains leads to significant reduction of the thermal relaxation time. The relaxation time of functionalized fullerenes becomes independent from the functionalizing chain length beyond approximately 10 carbon segments; this can be understood in terms of thermal conductivity along the chain and heat transfer between the chain and the solvent.     

Trifluoromethyl derivatives of insoluble small-HOMO-LUMO-gap hollow higher fullerenes. NMR and DFT structure elucidation of C2-(C74-D3h)(CF3)12, Cs-(C76-Td(2))(CF3)12, C2-(C78-D3h(5))(CF3)12, Cs-(C80-C2v(5))(CF3)12, and C2-(C82-C2(5))(CF3)12

Reaction of a mixture of insoluble higher fullerenes with CF3I at 500 degrees C produced a single abundant isomer of C74(CF3)12, C76(CF3)12, and C80(CF3)12, two abundant isomers of C78(CF3)12 and C82(CF3)12, and an indeterminant number of isomers of C84(CF3)12. Using a combination of 19F NMR spectroscopy, DFT calculations, and the structures and spectra of previously reported fullerene(CF3)n compounds, the most-probable structures of six of the seven isolated compounds were determined to be specific isomers of C2-(C74-D3h)(CF3)12, Cs-(C76-Td(2))(CF3)12), C2-(C78-D3h(5))(CF3)12), Cs-(C80-C2v(5))(CF3)12), C2-(C82-C2(5))(CF3)12), and C2-(C82-C2(3))(CF3)12) containing ribbons and/or loops of edge-sharing para-C6(CF3)2 hexagons. The seventh isolated compound is a C1 isomer of C78(CF3)12 containing two such ribbons. This set of compounds represents only the second reported isolable compound with the hollow C74-D3h cage and the first experimental evidence for the existence of the hollow fullerenes C76-Td(2), C78-D3h(5), C80-C2v(5), and C82-C2(5) in arc-discharge soots.     

Theoretical study of endohedral C36 and its dimers

It is found that atoms of lithium and carbon can be encapsulated in C(36) on the basis of the calculation of their encapsulation energies using density functional theory. Specifically, they can be encapsulated in C(36) better than in C(60) despite the smaller (70%) cavity size of the former. In C@C(36), the encapsulated carbon atom forms covalent bonds with the carbon atoms of the cage, which is in contrast with the case of N@C(60.) Two isomers are expected to be in an equilibrium which involves spin quenching and generation. Li@C(36) and C@C(36) are expected to exist in the form of dimers with nonendohedral fullerenes, i.e., as Li@C(36)-C(36) and C@C(36)-C(36). Three stable isomers were found for the former (A, B, and C). Equilibrium between A and C as well as that between B and C is accompanied by spin transfer between two fullerene units, while that between A and B is not. The two stable isomers in C@C(36)-C(36) form an equilibrium accompanied by spin quenching and generation, allowing the dimer to be potentially useful for molecular devices.     

Preparation of η2-complexes of fullerenes by reduction. Crystal structure and optical properties of {Ni(dppp)·(η2-C70)}·(C6H4Cl2)0.5

The new reduction method for preparation of η(2)-complexes of fullerenes with nickel-1,3-bis(diphenylphosphino)propane has been developed in which Ni(dppp)Cl(2) and C(60)(C(70)) mixtures are reduced with sodium tetraphenylborate. Single crystals of the first η(2)-complex of nickel with fullerene C(70): {Ni(dppp)·(η(2)-C(70))}·(C(6)H(4)Cl(2))(0.5) (1) (C(6)H(4)Cl(2) = o-dichlorobenzene) have been obtained as well as the previously described complex with fullerene C(60): {Ni(dppp)·(η(2)-C(60))}·(Solvent) (2). The crystal structure of 1 has been solved to show the coordination of nickel to the C-C bond of C(70) at the 6-6 ring junction of η(2)-type to form Ni-C(C(70)) bonds of 1.929-1.941(2) ? length, the shortest M-C bonds among those known for η(2)-complexes of fullerenes C(60) and C(70). The length of the C-C bond to which Ni atom is coordinated (1.494(3) ?) is noticeably longer than the average length of these bonds in C(70) (1.381(2) ?). Optical spectra of 1 in the IR- and UV-visible ranges have been analyzed to show the splitting of some C(70) bands due to C(70) symmetry lowering. The complex has a red-brown color in solution and manifests three bands in the visible range at 379, 467 and 680 nm. The solution of 1 is air sensitive since air exposure restores the color and absorption bands of the starting C(70) at 383 and 474 nm.     

Influence of electron doping on the hydrogenation of fullerene C60 : a theoretical investigation

The influence of electron attachment on the stability of the mono- and dihydrogenated buckminsterfullerene C(60) was studied using density functional theory and semiempirical molecular orbital techniques. We have also assessed the reliability of computationally accessible methods that are important for investigating the reactivity of graphenic species and surfaces in general. The B3LYP and M06L functionals with the 6-311+G(d,p) basis set and MNDO/c are found to be the best methods for describing the electron affinities of C(60) and C(60)H(2) . It is shown that simple frontier molecular orbital analyses at both the AM1 and B3LYP/6-31G(d) levels are useful for predicting the most favourable position of protonation of C(60)H(-) , that is, formation of the kinetically controlled product 1,9-dihydro[60]fullerene, which is also the thermodynamically controlled product, in agreement with experimental and previous theoretical studies. We have shown that reduction of exo- and endo-C(60)H makes them more stable in contrast to the reduction of the exo,exo-1,9-C(60)H(2) , reduced forms of which decompose more readily, in agreement with experimental electrochemical studies. However, most other dihydro[60]fullerenes are stabilized by reduction and the regioselectivity of addition is predicted to decrease as the less stable isomers are stabilized more by the addition of electrons than the two most stable ones (1,9 and 1,7).     

The exohedral Diels-Alder reactivity of the titanium carbide endohedral metallofullerene Ti2C2@D(3h)-C78: comparison with D(3h)-C78 and M3N@D(3h)-C78 (M=Sc and Y) reactivity

The chemical functionalization of endohedral (metallo)fullerenes has become a main focus of research in the last few years. It has been found that the reactivity of endohedral (metallo)fullerenes may be quite different from that of the empty fullerenes. Encapsulated species have an enormous influence on the thermodynamics, kinetics, and regiochemistry of the exohedral addition reactions undergone by these species. A detailed understanding of the changes in chemical reactivity due to incarceration of atoms or clusters of atoms is essential to assist the synthesis of new functionalized endohedral fullerenes with specific properties. Herein, we report the study of the Diels-Alder cycloaddition between 1,3-butadiene and all nonequivalent bonds of the Ti(2)C(2)@D(3h)-C(78) metallic carbide endohedral metallofullerene (EMF) at the BP86/TZP//BP86/DZP level of theory. The results obtained are compared with those found by some of us at the same level of theory for the D(3h)-C(78) free cage and the M(3)N@D(3h)-C(78) (M=Sc and Y) metallic nitride EMFs. It is found that the free cage is more reactive than the Ti(2)C(2)@D(3h)-C(78) EMF and this, in turn, has a higher reactivity than M(3)N@D(3h)-C(78). The results indicate that, for Ti(2)C(2)@D(3h)-C(78), the corannulene-type [5,6] bonds c and f, and the type B [6,6] bond 3 are those thermodynamically and kinetically preferred. In contrast, the D(3h)-C(78) free cage has a preference for addition to the [6,6] 1 and 6 bonds and the [5,6] b bond, whereas M(3)N@D(3h)-C(78) favors additions to the [6,6] 6 (M=Sc) and [5,6] d (M=Y) bonds. The reasons for the regioselectivity found in Ti(2)C(2)@D(3h)-C(78) are discussed.     

Preparation and characterization of the fullerene diols 1,2-C(60)(OH)(2), 1,2-C(70)(OH)(2), and 5,6-C(70)(OH)(2)

The simple fullerene diols C(60)(OH)(2) and C(70)(OH)(2) were prepared by addition of RuO(4) followed by acid hydrolysis. The 1,2-C(60)(OH)(2) isomer was formed from C(60), and two isomers (1,2 and 5,6) of C(70)(OH)(2) were formed in the RuO(4) hydroxylation of C(70). These compounds are much more soluble in THF and dioxane than the parent fullerenes. More highly hydroxylated materials are formed as well.     

Synthesis and X-ray or NMR/DFT structure elucidation of twenty-one new trifluoromethyl derivatives of soluble cage isomers of C76, C78, C84, and C90

Adding 1% of the metallic elements cerium, lanthanum, and yttrium to graphite rod electrodes resulted in different amounts of the hollow higher fullerenes (HHFs) C76-D2(1), C78-C2v(2), and C78-C2v(3) in carbon-arc fullerene-containing soots. The reaction of trifluoroiodomethane with these and other soluble HHFs at 520-550 degrees C produced 21 new C76,78,84,90(CF3)n derivatives (n = 6, 8, 10, 12, 14). The reaction with C76-D2(1) produced an abundant isomer of C2-(C76-D2(1))(CF3)10 plus smaller amounts of an isomer of C1-(C76-D2(1))(CF3)6, two isomers of C1-(C76-D2(1))(CF3)8, four isomers of C1-(C76-D2(1))(CF3)10, and one isomer of C2-(C76-D2(1))(CF3)12. The reaction with a mixture of C78-D3(1), C78-C2v(2), and C78-C2v(3) produced the previously reported isomer C1-(C78-C2v(3))(CF3)12 (characterized by X-ray crystallography in this work) and the following new compounds: C2-(C78-C2v(3))(CF3)8; C2-(C78-D3(1))(CF3)10 and C(s)-(C78-C2v(2))(CF3)10 (both characterized by X-ray crystallography in this work); C2-(C78-C2v(2))(CF3)10; and C1-C78(CF3)14 (cage isomer unknown). The reaction of a mixture of soluble higher fullerenes including C84 and C90 produced the new compounds C1-C84(CF3)10 (cage isomer unknown), C1-(C84-C2(11))(CF3)12 (X-ray structure reported recently), D2-(C84-D2(22))(CF3)12, C2-(C84-D2(22))(CF3)12, C1-C84(CF3)14 (cage isomer unknown), C1-(C90-C1(32))(CF3)12, and another isomer of C1-C90(CF3)12 (cage isomer unknown). All compounds were studied by mass spectrometry, (19)F NMR spectroscopy, and DFT calculations. An analysis of the addition patterns of these compounds and three other HHF(X) n compounds with bulky X groups has led to the discovery of the following addition-pattern principle for HHFs: In general, the most pyramidal cage C(sp(2)) atoms in the parent HHF, which form the most electron-rich and therefore the most reactive cage C-C bonds as far as 1,2-additions are concerned, are not the cage C atoms to which bulky substituents are added. Instead, ribbons of edge-sharing p-C6(X)2 hexagons, with X groups on less pyramidal cage C atoms, are formed, and the otherwise "most reactive" fullerene double bonds remain intact.     

Passive transport of C60 fullerenes through a lipid membrane: a molecular dynamics simulation study

To investigate the implications of the unique properties of fullerenes on their interaction with and passive transport into lipid membranes, atomistic molecular dynamics simulations of a C60 fullerene in a fully hydrated di-myristoyl-phoshatidylcholine lipid membrane have been carried out. In these simulations the free energy and the diffusivity of the fullerene were obtained as a function of its position within the membrane. These properties were utilized to calculate the permeability of fullerenes through the lipid membrane. Simulations reveal that the free energy decreases as the fullerene passes from the aqueous phase, through the head group layer and into the hydrophobic core of the membrane. This decrease in free energy is not due to hydrophobic interactions but rather to stronger van der Waals (dispersion) interactions between the fullerene and the membrane compared to those between the fullerene and (bulk) water. It was found that there is no free energy barrier for transport of a fullerene from the aqueous phase into the lipid core of the membrane. In combination with strong partitioning of the fullerenes into the lipidic core of the membrane, this "barrierless" penetration results in an astonishingly large permeability of fullerenes through the lipid membrane, greater than observed for any other known penetrant. When the strength of the dispersion interactions between the fullerene and its surroundings is reduced in the simulations, thereby emulating a nanometer sized hydrophobic particle, a large free energy barrier for penetration of the head group layer emerges, indicating that the large permeability of fullerenes through lipid membranes is a result of their unique interaction with their surrounding medium.     

C60 buckminsterfullerene high yields unraveled

Recently Irle, Morokuma, and collaborators have carried out a series of quantum chemical molecular dynamics simulations of carbon clustering. The results of these computer experiments are that carbon clusters of size greater than 60 atoms are rapidly formed, anneal to giant fullerenes, and then these fullerenes shrink. The simulation could not be carried to long enough times for the shrinking to reach C60, but they propose reasonably that this shrinking process ultimately forms buckminsterfullerene. However, these simulations do not reveal the force driving the shrinking process. Here, this driving force for shrinking is found to be reactions in which C2 is swapped between fullerenes. The key element is that for typical fullerenes the equilibrium constants for such C2 interchanges are near unity, resulting in expansion of the breadth of the fullerene distribution in an annealing process. When fullerenes of 60 or 70 atoms are populated by shrinking, they fall into the local energy minimum of buckminsterfullerene or D5h C70. This simple mechanism accounts for the high yields (>20%) of buckminsterfullerene that can be achieved in pure carbon systems.     

Epitaxial supramolecular assembly of fullerenes formed by using a coronene template on a Au(111) surface in solution

Characteristic properties of the coronene layer formed on Au(111) for the epitaxial growth of various fullerenes are described. The electrochemical behavior of the coronene adlayer prepared by immersing a Au(111) substrate into a benzene solution containing coronene was investigated in 0.1 M HClO4. The as-prepared coronene adlayer on Au(111) revealed a well-defined (4 x 4) structure. Structural changes of the array of coronene molecules induced by potential manipulation were clearly observed by in situ scanning tunneling microscopy (STM). Supramolecularly assembled layers of fullerenes such as C60, C70, C60-C60 dumbbell dimer (C120), C60-C70 cross-dimer (C130), and C60 triangle trimer (C180) were formed on the well-defined coronene adlayer on the Au(111) surface by immersing the coronene-adsorbed Au(111) substrate into benzene solutions containing those molecules. The adlayers thus prepared were characterized by comparison with those which were directly attached to the Au(111) surface. The C60 molecules formed a honeycomb array with an internal structure in each C60 cage on the coronene adlayer, whereas C70 molecules were one-dimensionally arranged with the same orientations. The dimers, C120 and C130 molecules, formed an identical structure with c(11 x 4 radical3)rect symmetry. For the C130 cross-dimer molecule, C60 and C70 cages were clearly recognized at the molecular level. It was difficult to identify the adlayer of the C180 molecule directly attached to Au(111); however, individual C180 molecules could be recognized on the coronene-modified Au(111) surface. Thus, the adlayer structures of those fullerenes were strongly influenced by the underlying coronene adlayer, suggesting that the insertion of a coronene adlayer plays an important role in the formation of supramolecular assemblies of fullerenes.