Probing C<sub>84</sub>-embedded Si Substrate Using Scanning Probe Microscopy and Molecular Dynamics

Mon-Shu Ho; Chih-Pong Huang; Jyun-Hwei Tsai; Che-Fu Chou; Wen-Jay Lee

doi:10.3791/54235

JoVE Journal > Engineering

Engineering

Probing C₈₄-embedded Si Substrate Using Scanning Probe Microscopy and Molecular Dynamics

Published: September 28, 2016

doi:

10.3791/54235

Mon-Shu Ho¹, Chih-Pong Huang², Jyun-Hwei Tsai³, Che-Fu Chou¹, Wen-Jay Lee³

¹Department of Physics and Institute of Nanoscience,National Chung Hsing University, ²Metallurgy Section, Materials & Electro-Optics Research Division,National Chung-Shan Institute of Science and Technology, ³National Center for High-Performance Computing

Summary

This paper reports the nanomaterial fabrication of a fullerene Si substrate inspected and verified by nanomeasurements and molecular dynamic simulation.

Abstract

This paper reports an array-designed C₈₄-embedded Si substrate fabricated using a controlled self-assembly method in an ultra-high vacuum chamber. The characteristics of the C₈₄-embedded Si surface, such as atomic resolution topography, local electronic density of states, band gap energy, field emission properties, nanomechanical stiffness, and surface magnetism, were examined using a variety of surface analysis techniques under ultra, high vacuum (UHV) conditions as well as in an atmospheric system. Experimental results demonstrate the high uniformity of the C₈₄-embedded Si surface fabricated using a controlled self-assembly nanotechnology mechanism, represents an important development in the application of field emission display (FED), optoelectronic device fabrication, MEMS cutting tools, and in efforts to find a suitable replacement for carbide semiconductors. Molecular dynamics (MD) method with semi-empirical potential can be used to study the nanoindentation of C₈₄-embedded Si substrate. A detailed description for performing MD simulation is presented here. Details for a comprehensive study on mechanical analysis of MD simulation such as indentation force, Young's modulus, surface stiffness, atomic stress, and atomic strain are included. The atomic stress and von-Mises strain distributions of the indentation model can be calculated to monitor deformation mechanism with time evaluation in atomistic level.

Introduction

Fullerene molecules and the composite materials they comprise are distinctive among nanomaterials due to their excellent structural characteristics, electronic conductivity, mechanical strength, and chemical properties^1-4. These materials have proven highly beneficial in a range of fields, such as electronics, computers, fuel cell technology, solar cells, and field emission technology^5,6.

Among these materials, silicon carbide (SiC) nanoparticle composites have received particular attention thanks to their wide band gap, high thermal conductivity and stability, high electrical breakdown ability, and chemical inertness. These benefits are particularly obvious in optoelectronic devices, metal-oxide-semiconductor field-effect transistors (MOSFET), light-emitting diodes (LEDs), and high-power, high-frequency, and high-temperature applications. However, high density defects commonly observed on the surface of conventional silicon carbide can have detrimental effects on the electronic structure, even leading to device failure^7,8. Despite the fact that the application of SiC has been studied since 1960, this particular unresolved problem remains.

The aim of this study was the fabrication of a C₈₄-embedded Si substrate heterojunction and subsequent analysis to obtain a comprehensive understanding of the electronic, optoelectronic, mechanical, magnetic, and field emission properties of the resulting materials. We also addressed the issue of using numerical simulation to predict the characteristics of nanomaterials, through the novel application of molecular dynamics calculations.

Protocol

NOTE: The paper outlines the methods used in the formation of a self-assembled fullerene array on the surface of a semiconducting substrate. Specifically, we present a novel method for the preparation of a fullerene-embedded silicon substrate for use as a field emitter or substrate in microelectromechanical systems (MEMS), and optoelectronic devices in high-temperature, high-power, applications as well as in high-frequency devices^9-13.

1. Fabrication of Hexagonal-closed-packaged (HCP) Overlayer of C₈₄ on Si Substrate

Prepare Clean Si(111) Substrate
1. Subject Si substrate to RCA (Radio Corporation of America) cleaning, involving the application of a solvent followed by heating in an ultra-high vacuum system for the removal of the oxide layer and impurities from the surface of the substrate (see supporting material).
  NOTE: Herein, the term "UHV-ultra high vacuum system" refers to a vacuum below 1 x 10^-8 Pa used in the preparation of a Si(111).
Deposit C₈₄ on Silicon Surface Using Thermal Evaporation in a UHV System
1. Pre-heat a K-cell evaporator with external power supply through heating filaments to 500 °C to promote the outgassing of impurities.
2. Load C₈₄ nanoparticles into a K-cell container. Resistively heat the K-cell to 650 °C. Vaporize C₈₄ nanoparticles as C₈₄ nanoparticles in the container compose vapors. Evaporate C₈₄ nanoparticles in straight lines until nanoparticles strike a Si substrate through a controlled valve at pressure below 5 x 10^-8 Pa.
Embed C₈₄ Molecules within Si Surface via Self-assembly Mechanism
1. Pre-anneal Si(111) substrate in an ultra-high vacuum system at 900 °C to obtain (1×1) structures. Reduce the temperature to 650 °C for 30 min for the deposition of the C₈₄ nanoparticles on the surface of the substrate.
2. Anneal the Si substrate at ~750°C for 12 hr, during which time the powdered-C₈₄ nanoparticles self-assemble into a highly uniform fullerene array on the surface of the Si(111) substrate.
  NOTE: Herein, the term "highly uniform fullerene array" refers to the uniform distribution of fullerene on the substrate, in which most of the nanoparticles are oriented in a compact arrangement perpendicular to the surface of the substrate. This configuration helped to ensure that the vertical height of the fullerene array was essentially identical in all samples.

2. Measurements of Electronic Properties of C₈₄-embedded Si Substrate

Measure Local Electronic Density of States Using UHV-scanning Tunneling Microscopy
1. Measure I-V curves of specific atoms using UHV-SPM
2. Place C₈₄-embedded Si substrate on a SPM sample holder. Introduce the holder into an UHV-STM (scanning tunneling microscope) scanning head system. Sweep applied sample bias from -5 V to 5 V.
3. Click on "I-V" measurement item to measure the tunneling current I at atomic resolution. Choose at least 20 particular locations on the C₈₄-embedded Si substrate for measurements. Calculate the mean value of tunneling current I over 20 particular locations. Derive I as a function of voltage. Plot I-V curves.
4. Calculate the derivative of I(V) with respect to V. Convert the I-V curves to dI/dV as a function of voltage in order to determine the local electronic state of the C₈₄-embedded Si substrate.
Measure Band Gap Energy
1. Obtain I-V curves according to the procedures in 2.1.2 and 2.1.3 from the following: Si(111)-7×7 surface, Si(111)-1×1 surface, single individual C₈₄ nanoparticles on Si, 7-19 C₈₄ clusters on Si, 20-50 C₈₄ clusters on Si, and a monolayer of C₈₄ embedded within Si surface.
2. Calculate the derivative of I(V) with respect to V. Convert the I-V curves to dI/dV curves to measure the HOMO-LUMO energy differences (referred to band gap energy) in each measurement location, as shown in Figure 2a.
Obtain Field Emission (FE) Properties
1. Place C₈₄-embedded Si substrate on a FE sample holder. Insert the holder into FE analysis chamber. Evacuate the chamber to a pressure of approximately 5 x 10^-5 Pa for FE measurement.
  NOTE: The C₈₄-embedded silicon substrate functioned as the cathode and a copper probe with a cross-sectional area of ~0.71 mm² functioned as the anode. The distance between the cathode and anode was approximately 590 µm.
2. Increase applied voltage manually on substrate from 100 V to 1,100 V. Measure the corresponding field emission current as a function of applied voltage using a high-voltage source measurement unit with current amplifier.
3. Calculate the Fowler-Nordheim field emission correlation according to the work function ~5 eV as shown in Figure 2b.
4. Obtain the geometric field enhancement factor (β) as follows: F(field) = β(V/d) with a β value of approximately 4,383.
5. Obtain the electrical breakdown field under vacuum based on the slope of the natural logarithm (J/E²) vs (1/E), which gave us a value of ~4.0 x 10⁶ V/cm for the C₈₄-embedded Si substrate as shown in Figure 2c.
Optoelectronic Properties
1. Transfer testing substrate to an optical emission measurement system. Focus a He-Cd laser source with 325 nm emissions on the substrate that is located in the center of the sample compartment. Set up a spectrometer in a suitable position. Use a spectrometer to acquire the photoluminescence spectrum by collecting and analyzing emitting photons. The optoelectronic result is shown in Figure 2d.

3. Measurements of Surface Magnetism

Obtain MFM (Magnetic Force Microscopy) Topography.
1. Magnetize samples of C₈₄-embedded Si prior to MFM measurements by applying a magnet with a field strength of approximately 2 kOe.
2. Place the magnetized sample on an MFM sample stage. Click on "Obtain MFM topography" item. Observe the microstructure of the fullerene in the magnetic domain embedded within the Si substrate using MFM in lift mode with the application of magnetization perpendicular to the surface of the sample.
3. Use a nano-scale PPP-MFMR cantilever for MFM measurements (Figure 3a). Determine the surface magnetism if MFM topography appears darker(brighter) when the magnetic moment of tip is in the same(opposite) direction of the substrate moment.
SQUID (Superconducting Quantum Interference Device) Measurement
1. Prepare monolayer of C₈₄-embedded Si substrate and C₈₄ clusters on C₈₄ embedded Si substrate.
2. Magnetize samples of C₈₄-embedded Si and C₈₄ clusters on C₈₄ embedded Si substrate prior to SQUID experiments by applying a magnet with a field strength of approximately 2 kOe.
3. Place the sample in an SQUID. Apply a sweeping magnetic field in a range of ~2 kOe. Obtain the magnetization loops plotted versus the external magnetic field in SQUID measurements at room temperature.
  NOTE: The typical M-H curve for a ferromagnetic material can be obtained as shown in Figure 3b.

4. Measurement of Nanomechanical Properties by AFM

NOTE: Atomic force microscopy (AFM) provides a powerful tool for the characterization of material and mechanical properties at the micro- and nano-scales in air as well as in a UHV environment

Measure the Stiffness of C₈₄ Embedded Si Substrate Under Atmospheric Conditions
1. Place the substrate on an AFM sample stage. Drag a sharp tip over the substrates using a scanner. Monitor the displacements of the tip as a measure of tip-sample interaction forces. Record the movements at many tip-sample distances along vertical direction in a certain position by clicking on "force measurement" item.
2. Obtain force measurements using a AFM under atmospheric conditions from an RCA-cleaned Si substrate with 2-3 nm layer of natural oxide as well as from a C₈₄-embedded Si substrate and a Si substrate coated with a thin film of SiC.
3. Using AFM software, plot Force-distance curves under atmospheric conditions.
  NOTE: The AFM cantilever was a Si probe with a tip radius of ~5-20 nm and spring constant of ~40 N/m.
Measure the Stiffness of C₈₄ Embedded Si Substrate in UHV Chamber
1. Obtain force measurements according to the guidance of 4.1.1 using a AFM in a UHV system from an RCA-cleaned Si substrate, a clean Si(111)-7×7 surface, a C₈₄-embedded Si substrate, substrate and a Si substrate coated with a thin film of SiC.
2. Plot Force-distance curves in a UHV system. Note: The AFM cantilever was a Si probe with a tip radius of ~5-20 nm and spring constant of ~40 N/m. Figure 4 presents the force-distance analysis of disordered Si surface, 7 x 7 surface, single self-assembled layer of C₈₄ embedded within Si surface, and Si surface, as determined using UHV-AFM.

5. Measurement of Nanomechanical Properties by MD Simulation

NOTE: In the simulation section, OVITO¹⁶ (open-source visualization software) and, oSSD¹⁷ (Open surface structure database) are used to create the simulation model and results visualization. LAMMPS¹⁴ (an open-source molecular dynamics (MD) simulation package) is employed to perform the nanoindentation simulation and analyze the simulation results¹⁵. All the simulation jobs are performed with parallel computing in the Advanced Large-scale Parallel Supercluster (ALPS) of NCHC.
NOTE: To study the C₈₄ monolayer/Si substrate heterojunction by using MD simulation, one should prepare a simulation model by several steps to obtain a relaxed C₈₄ monolayer embedded into the Si substrate. Note that it is difficult to generate an exactly the same structure from the experimental data, because of the complex of the inter structure between C₈₄ monolayer and Si (111) substrate heterojunction. As a result, we use an artificial way to generate the simulation model with several steps of procedure, which is illustrated in Figure 5. The details are described in the following protocols. We describe how to setup the parameter of MD in LAMMPS, establish a relaxed C₈₄fullerene monolayer embedded into a substrate, perform an indentation procedure, and analyze the simulation results.

Parameter Setting in LAMMPS input File
1. Use boundary command to set the periodic boundary conditions in the x- and the y-directions.
2. Use "fix velocity" command to assign the initial velocity with a Gaussian distribution on each atom of the system, randomly.
3. Use "fix pair_style" command to assign Tersoff¹⁸ and AIREBO¹⁹ potentials to describe the Si-Si and Si-C interaction and the C-C interaction, respectively.
4. Use "fix nvt" and "fix npt" command to adopt the Nosè-Hoover method²⁰ to ensure the system remains at the desired temperature and pressure to generate a canonical and isothermal-isobaric ensemble²⁰, in which system the velocity-Verlet algorithm²⁰ is employed to predict the trajectories of the atoms. Use both "fix nvt" and "run" commands to set a cooling rate of 3 K/psec for annealing process.
5. Use "timestep" command to set a time step of 0.2 fsec as the time integration.
6. Use "fix wall/reflect" command to adopt a reflected wall to confine the degree of freedom (5.3.2).
7. Use "region" and "group" to divide the substrate into different control layers (5.4.3): Newtonian atom layer, a thermal control layer, and a bottom fixed layer, which can be set up by using "fix nve", "fix nvt", and "fix setforce" commands, respectively.
8. Use "region" and "create_atoms" commands to create a spherical probe.
9. Use "fix move" command to embed the C84 monolayer into the substrate (5.4.2) and move the probe during the simulation (5.5.2).
10. Use "run" command to perform MD simulation.
11. Use "compute force" (5.6.1) and "compute stress/atom" (5.6.4) commands to evaluate the atomic stress and indentation force.
  NOTE: In the following, except the structure establishing, all the steps were done by LAMMPS script.
Use oSSD and OVITO to Preparation of Silicon(111) 7 x 7 Surface.
1. Turn on the oSSD software. Click on the "search" button. "Search criteria" panel is presented. Chose Si substrate, elemental type, reconstructed structure, semiconductor elec, diamond lattice, 111 face and 7 x 7 pattern. Click on the "Search" and "Accept" buttons. "Structure list" panel is presented. Click the desired structure (i.e., Si(111) 7 x 7). Click the "File" button. Save the coordination file as .xyz file.
  NOTE: We point out that the structural database extracted from oSSD is not large enough for our indentation simulation. As a result, we rebuild a larger and thicker substrate by the following steps.
2. Turn on the OVITO software. Load the .xyz file into OVITO. Use "Slice" command to capture a supercell of the Si(111) 7 x 7 surface with size of 26.878 x 46.554 Å² in x and y direction. Export the data file. Use "Slice" command to capture a supercell of the bottom Si(111) substrate with size of 26.878 x 46.554 x 9.7 Å³. Use "Show periodic images" command to duplicate the supercell 12 times in z direction. Export the data file.
3. Combine the data files of Si(111) 7 x 7 surface and the Si(111) substrate models by Notepad++ (a free source code editor). Finally, load the combined data into OVITO. Use "Show periodic images" to duplicate a 5 x 3 supercell in x and y directions to enlarge the size of substrate.
4. Use LAMMPS to perform a 20 psec MD simulation time for relaxing the simulation model. In the following, perform a quench process from 1,550 K to room temperature for 500 psec simulation time. Finally, perform a 10 psec simulation time for the final relaxation process.
Preparation of C₈₄Fullerene Monolayer
1. Download the coordination file of the Optimized structure of C₈₄ fullerene from the web²¹ and write a FORTRAN program to replicate 49 C₈₄ fullerenes arranged in a honeycomb structure.
2. Use LAMMPS to setup reflect walls upon and below the C₈₄monolayer to ensure that the molecules stay on a plan. Perform a MD simulation time for 200 psec to relax the simulation model. In the following, perform a quench process from 700 K to room temperature to obtain a glob minimum state for 500 psec simulation time. Finally, perform a 10 psec simulation time for final relaxation process.
Establish the Indentation Model of C₈₄Fullerene Monolayer on Silicon (111) 7 x 7 Surface.
1. Write a FORTRAN code to lay the C₈₄ monolayer on the Si(111) 7 x 7 surface with distance of 3 Å to establish the indentation model.
2. Use LAMMPS to embed the C₈₄ monolayer into the substrate with depth of 2~3 Å. In the following, run a 40 psec simulation time for system relaxation. Finally, anneal the system to room temperature.
3. Divide the silicon substrate into a top Newtonian atom layer, a thermal control layer, and a bottom fixed layer, which are 0.7, 2, and 5.3 nm in thickness, respectively. The C₈₄ monolayers were also modeled as a Newtonian atom.
Indentation Process of MD
1. Use LAMMPS to create a spherical probe with 5nm in diameter upon the C₈₄/Si(111) 7 x 7 surface mode (Figure 5). The probe is set as a rigid body. Specify a constant velocity of 10 m/sec on the probe to move downward toward the specimen in the indentation process.
2. Move the probe downward to the specimen at a constant speed until the specific loading depth (i.e., including the cases of 1.5, 2.5, 4.5, 10, 15, 20, and 30 Å so as to explore the effect of the C₈₄ fullerenes monolayer on the Si substrate, where the size of C₈₄ fullerene is 11 Å) in the loading process. Hold the probe in the substrate in the holding process to allow for the relaxation of atoms. Finally, extract the probe from the substrate at a constant speed in the retraction process.
Calculation and Analysis
1. Calculate the indentation force by summing the vertical force of atoms in the probe according to the following formulas:
  (1)
2. Extracted the reduced modulus and stiffness from the force-distance curve of indentation. Based on Oliver and Pharr's method²², a linear relation can be derived between the Young's modulus and the unloading stiffness. The stiffness (i.e., the slope of the initial portion) of the unloading curve is defined as
  (2)
  where P, h, A, and E_r are the indentation load, elastic displacement of probe, projected area of the indentation, and reduced modulus. β (=1 for circular indenter) is the shape modification factor. The relationship between reduced modulus and Young's modulus can be written as
  (3)
  where E and v are the Young's modulus and Poisson's ratio for the specimen and E_i and v_i are the the Young's modulus and Poisson's ratio for the indenter.
3. Calculate the hardness by definition of H = P_max/A, where P_max and A are the maximum indentation force and projected area of probe.
4. Calculate the virial atomic stress²² on the m plane of the substrate in the n-direction by
  (4)
  where m_i is the mass of atom i; and are the velocity components of atom i in the m– and n-directions, respectively; V_i is the volume assigned around atom i; N_sis the number of particles contained within region S, where S is defined as the region of atomic interaction; Φ(r_ij) is the potential function; r_ij is the distance between atoms i and j, and and are the m– and n-direction components of the vector from atom i to atom j.
5. Use OVITO to show the von-Mises strain of each atom invariant according to the following formulas:
  (5)

Representative Results

A monolayer of C₈₄ molecules on a disordered Si(111) surface was fabricated using a controlled self-assembly process in a UHV chamber. Figure 1 shows a series of topographic images measured by UHV-STM with various degrees of coverage: (a) 0.01 ML, (b) 0.2 ML, (c) 0.7 ML, and (d) 0.9 ML. The electronic and optical properties of the C₈₄ embedded Si substrate were also investigated using a variety of surface analysis techniques, such as STM and PL (Figure 2). The excellent material properties of the resulting samples demonstrate how nanotechnology can be used for the control of matter at the atomic- and nano-scales. The MFM and SQUID results in Figure 3 show the surface magnetism of C₈₄ embedded substrate. Figure 4 presents the UHV-AFM results that refer to the nanomechanics of proposed substrate. Our experimental results demonstrate the potential of the C₈₄ embedded silicon substrate as an alternative to semiconductor carbide in nanoelectronic devices for high-temperature, high-power, high-frequency applications as well as in magnetic and MEMS devices (in Figure 4).

In the simulation section, all the procedures are completed by using LAMMPS to perform the MD simulations. The mechanical properties (indentation force and contact stress) of the fullerene embedded substrate is calculated and shown in Figure 6. The von-Mises strain analysis of atoms at different time step are used to characterize the local deformation. The corresponding snapshots as a function of indentation depth can be seen in the inserts of Figure 6, which were calculated and visualized by OVITO. The results of indentation force as a function of indentation depth are used to calculate the hardness H (Figure 7a), reduced modulus E_r (Figure 7b), and loading stiffness S (Figure 8) of the C₈₄ monolayer. The results can be compare with that determined by experiment and provides a more detail point of view to interpret the variation of the mechanical property.

Figure 1: C₈₄ Embedded Si substrate with different coverage. Series of STM topographic images (40 x 40 nm²) showing C₈₄ molecules adsorbed on Si(111) surface at a negative sample bias of 2 V, as measured by UHV-STM with various degrees of coverage: (a) 0.01 ML, (b) 0.2 ML, (c) 0.7 ML, and (d) 0.9 ML.

Figure 2: Electronic properties measurements on C₈₄ Embedded Si substrate. (a) I-V curves and differential derivative conductance (dI/dV) vs the voltage curve of a single self-assembled layer of C₈₄, as determined by UHV-STM; (b) Field emission current density vs electric field curve; (c) Corresponding F-N plot of surface with embedded C₈₄, as measured using a source-measure unit; (d) Photoluminescence spectrum of single self-assembled layer of C₈₄. Re-print with permission from reference¹². Please click here to view a larger version of this figure.

Figure 3: Surface magnetism on C₈₄ Embedded Si substrate. (a) MFM image of Si Substrate embedded with C₈₄; (b) Magnetization loop plotted against external magnetic field Please click here to view a larger version of this figure.

Figure 4: Nanomechanical investigation on C₈₄ Embedded Si substrate. Force-distance analysis of disordered Si surface, 7 x 7 surface, single self-assembled layer of C₈₄ embedded within Si surface, and Si surface, as determined using UHV-AFM. Re-print with permission from reference¹¹. Please click here to view a larger version of this figure.

Figure 5: Flow chart for establishing simulation model. The dramatic illustrate the setting in MD simulation from a single layer C₈₄ and Si(111) 7 x 7 surface to a C₈₄ monolayer embed into Si(111) 7 x 7 model. The detail procedures can be seen in the section 5 of protocol. Please click here to view a larger version of this figure.

Figure 6: Indentation force and Contact stress analysis. Indentation force (black) and Contact stress (blue) of C₈₄ as a function of indentation depth. Inserts show the corresponding snapshots, where the different color indicates the corresponding von Mises strain (ε_vM) of all atoms. To clear display the strain localization, only the atoms with ε_vM > 0.08 are shown in the snapshot. Please click here to view a larger version of this figure.

Figure 7: Hardness and reduced modulus analysis. (a) Hardness and (b) reduced modulus variation as a function of indentation depth for the C₈₄monolayer on Si surface. Please click here to view a larger version of this figure.

Figure 8: Loading stiffness analysis. Loading stiffness as a function of distance determined by MD simulations compared with that by AFM experiments for C₈₄/Si. Modified from reference¹⁶. Please click here to view a larger version of this figure.

Discussion

In this study, we demonstrate the fabrication of a self-assembled monolayer of C₈₄ on a Si substrate through a novel annealing process (Figure 1). This process can also be used to prepare other kinds of nanoparticle-embedded semiconductor substrates. The C₈₄-embedded Si substrate was characterized at the atomic scale using UHV-STM (Figure 2), field emission spectrometer, photo-luminescence spectroscopy, MFM and SQUID (Figure 3).

The adhesion strength corresponding to nanomechanical properties (i.e., stress) of the C₈₄-embedded Si substrates can be measured using AFM (Figure 4). Our results demonstrate that the hardness of the proposed C₈₄-embedded Si substrate is comparable to that of SiC and Si surfaces, making it applicable as an abrasive material for cutting tools as well as a film in MEMS devices.

In the simulation section, the von Mises strain (ε_vM) analysis is capable to detect the local deformation of atomic structure, which is very difficult to be observed in experiment. However, it is not possible to characterize the phase transformation. Here, we suggest some useful indices such as coordination number and HA index²³ to examine the phase transformation. In the setting of the indentation model, we have to point out that the size of the substrate in plan direction must be at least three times larger than diameter of the probe for eliminating the size effect and the boundary condition limitation, which would affect the dynamics and force flow of atoms.

In addition, due to the time limit of MD simulation, to study the indentation process, the probe should exert the specimen with a very fast speed compared to that in experiment. We note that such a loading speed is too high to come out the long-time atomic diffusion and migration behavior, but it is still suitable to observe and describe the plastic deformation behavior and material properties under mechanical loading²⁴ because the results can be recognized as approximately quasi-static in nature²⁵. An alternative theory, named Parallel Replica dynamics (PRD)²⁶, has been developed to substantially accelerate the simulation time, but it requires heavy computing resources.

The data obtained from the MD simulation study is in agreement with the AFM indentation experiment (Figure 8); additionally, the hardness and reduced modulus of the C₈₄-embedded Si substrate are comparable to other Si substrates. These data suggest that C₈₄-embedded Si substrates can have significant impacts in optoelectronic and dilute magnetic semiconductor (DMS) nanodevices.

Disclosures

The authors have nothing to disclose.

Acknowledgements

The authors would like to thank the Ministry of Science and Technology of Taiwan, for their financial support of this research under Contract Nos. MOST-102-2923-E-492- 001-MY3 (W. J. Lee) and NSC-102- 2112-M-005-003-MY3 (M. S. Ho). Support from the High-performance Computing of Taiwan in providing huge computing resources to facilitate this research is also gratefully acknowledged.

Materials

Silicon wafer			Si(111) Type/Dopant: P/Boron Resistivity: 0.05-0.1 Ohm.cm
Carbon,C₈₄	Legend Star		C84 powder, 98%
Hydrochloric acid	Sigma-Aldrich	84422	RCA,37%
Ammonium	Choneye Pure Chemical		RCA,25%
Hydrogen peroxide	Choneye Pure Chemical		RCA,35%
Nitrogen	Ni Ni Air		high-pressure bottle,95%
Tungsten	Nilaco	461327	wire, diameter 0.3 mm, tip
Sodium hydroxide	UCW	85765	etching Tungsten wire for tip,
Acetone	Marcon Fine Chemicals	99920	suitable for liquid chromatography and UV-spectrophotometry
Methanol	Marcon Fine Chemicals	64837	suitable for liquid chromatography and UV-spectrophotometry
UHV-SPM	JEOL Ltd	JSPM-4500A	Ultrahigh Vacuum Scanning Tunneling Microscope and Ultrahigh Vacuum Atomic Force Microscope
Power supply	Keithley	237	High-Voltage Source-Measure Unit
SQUID	Quantum desigh	MPMS-7	Magnetic field strength: ± 7.0 Tesla, Temperature range: 2 ~ 400 K, Magnetic-dipole range:5 × 10^-7 ~ 300 emu
ALPS	National Center for High-performance Computing, Taiwan		Advanced Large-scale Parallel Supercluster, 177Tflops; 25,600 CPU cores; 73,728 GB RAM; 1074 TB storage

References

Kroto, H. W., Heath, J. R., O’Brien, S. C., Curl, R. F., Smalley, R. E. C60: Buckminsterfullerene. Nature. 318, 162-163 (1985).
Zhu, Z. P., Gu, Y. D. Structure of carbon caps and formation of fullerenes. Carbon. 34, 173-178 (1996).
Margadonna, S., et al. Crystal Structure of the Higher Fullerene C84. Chem. Mater. 10, 1742-1744 (1998).
Diederich, F., et al. The Higher Fullerenes: Isolation and Characterization of C76, C84, C90, C94, and C70O, an Oxide of D5h-C70. Science. 252, 548-551 (1991).
Lv, Z., Deng, Z., Xu, D., Li, X., Jia, Y. Efficient organic light-emitting diodes with C60 buffer layer. Displays. 30, 23-26 (2009).
Tokunaga, K. On the difference in electronic properties between fullerene C60 and C60X2. Chem. Phys. Lett. 476, 253-257 (2009).
Basa, D. K., Smith, F. W. Annealing and crystallization processes in a hydrogenated amorphous Si—C alloy film. Thin Solid Films. 192, 121-133 (1990).
Neudeck, P. G., Powell, J. A. Performance limiting micropipe defects in silicon carbide wafers. IEEE Electron Device Lett. 15, 63-65 (1994).
Huang, C. P., Su, C. C., Ho, M. S. Intramolecular Structures of C60 and C84 Molecules on Si(111)-7×7 Surfaces by Scanning Tunneling Microscopy. Appl. Surf. Sci. 254, 7712-7717 (2008).
Huang, C. P., Su, C. C., Su, W. S., Hsu, C. F., Ho, M. S. Nano-measurements of electronic and mechanical properties of fullerene embedded Si(111) surfaces. Appl. Phys. Lett. 97, 061908 (2010).
Huang, C. P., Hsu, C. F., Ho, M. S. Investigation of Fullerene Embedded Silicon Surfaces with Scanning Probe Microscopy. J. Nanosci. Nanotechnol. 10, 7145-7148 (2010).
Huang, C. P., Su, W. S., Su, C. C., Ho, M. S. Characteristics of Si(111) surface with embedded C84 molecules. RSC Adv. 3 (111), 9234-9239 (2013).
Ho, M. S., Huang, C. P. Method of forming self-assembled and uniform fullerene array on surface of substrate. US Patent. , (2015).
Plimpton, S. J. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comp. Phys. 117, 1-19 (1995).
Su, W. S., et al. Using a functional C84 monolayer to improve the mechanical properties and alter substrate deformation. RSC Adv. 5, 47498-47505 (2015).
Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO- the Open Visualization Tool. Modelling Simul. Mater. Sci. Eng. 18, 015012 (2010).
Tersoff, J. New empirical model for the structural properties of silicon. Phys. Rev. Lett. 56, 632 (1986).
Stuart, S. J., Tutein, A. B., Harrison, J. A. A reactive potential for hydrocarbons with intermolecular interactions. Journal of Chemical Physics. 112, 6472 (2000).
Rapaport, D. C. . The Art of Molecular Dynamics Simulations. , (1997).
Chandra, N., Namilae, S., Shet, C. Local elastic properties of carbon nanotubes in the presence of Stone-Wales defects. Phys. Rev. B. 69, 094101 (2004).
Honeycutt, J. D., Andersen, H. C. Molecular Dynamics Study of Melting and Freezing of Small Lennard-Jones Clusters. J. Phys. Chem. 91, 4950 (1987).
Schulz, M. J., Kelkar, A. D., Sundaresan, M. J. . Nanoengineering of Structural, Functional and Smart Materials. , 736 (2005).
Ju, S. P., Wang, C. T., Chien, C. H., Huang, J. C., Jian, S. R. The Nanoindentation Responses of Nickel Surfaces with Different Crystal Orientations. Molecular Simulation. 33, 905-917 (2007).
Mishin, Y., Suzuki, A., Uberuaga, B. P., Voter, A. F. Stick-slip behavior of grain boundaries studied by accelerated molecular dynamics. Phys. Rev. B. 75, (2007).

Play Video

PDF

DOI

DOWNLOAD MATERIALS LIST

Cite This Article

Ho, M., Huang, C., Tsai, J., Chou, C., Lee, W. Probing C₈₄-embedded Si Substrate Using Scanning Probe Microscopy and Molecular Dynamics. J. Vis. Exp. (115), e54235, doi:10.3791/54235 (2016).

Automatically Generated