
Web Release Date: December 8,
On the Formation of Single-Walled Carbon Nanotubes in Pulsed-Laser-Assisted Chemical Vapor Deposition
IFW Dresden, P.O. Box 270116, D-01171 Dresden, Germany, and Centre of Knowledge Based Technologies and Nanomaterials, Institute of Chemical and Environment Engineering, Szczecin University of Technology, 70-322 Szczecin, Poland
Received August 10, 2007
Revised Manuscript Received September 28, 2007
Abstract:
We report on a combined approach for single-wall carbon nanotube (SWCNT) synthesis merging pulsed-laser ablation with floating catalyst chemical vapor deposition. This hybrid route allows the synthesis of SWCNT with a defined mean diameter and diameter distribution. Further, a careful investigation of the resulting yield and obtained SWCNT diameter behavior reveals that controlling both the reaction temperature and the carbon availability (via pressure, flow rate, or type of carbon feedstock) is crucial. The complex influence of reaction temperature and carbon availability on the resulting yield and diameter behavior are well-explained within our previously described catalyst volume to surface area formation mechanism of SWCNT.
Single-walled carbon nanotubes (SWCNT)1 are of great interest to science and industry because of their electronic and vibronic properties,2 which allow the study of one-dimensional physics as well as applications in numerous devices. Among other things, the high aspect ratio and the tuneability of the electronic structure is of advantage. Synthesized SWCNT of high quality, with a narrow diameter distribution and a well-defined mean diameter are desired for several proposed applications as well as fundamental studies. Today, there exist three primary routes for their synthesis. In arc-discharge3, an arc is ignited between two electrodes consisting of carbon mixed with some metal catalyst. Although the yield is relatively good, the as-produced material includes many catalyst particles as well as amorphous carbon. In laser ablation (LA)4, a laser beam is used to evaporate a carbon/metal target. LA allows for narrow diameter distributions with appropriate tuning while maintaining a high yield, but the technique is expensive and not suited for mass production. Chemical vapor deposition (CVD)5 is the most promising method for a continuous high yield nanotube synthesis. There has been recent progress in CVD synthesis using tailored catalyst particles in order to grow carbon nanotubes in a specific diameter range6, 7. Alexandrescu et al.8 did some work toward combining catalytic CVD and lasers in using a floating catalyst method with a continuous wave CO2 laser beam for heating. Another form of laser-assisted chemical vapor deposition (LACVD) is the pyrolysis of metal−organic compounds9. These routes, however, do not use pulsed evaporation. Nishide et al.10 introduced a pulsed laser evaporation method in which a metal target is evaporated in a hydrocarbon atmosphere, namely, cobalt in ethanol. We further develop this promising approach, which blends the well-known diameter control from pulsed-laser evaporation systems with floating catalyst CVD11. In this study, we present a new approach to achieving enhanced control over the diameter distribution combined with an improved yield in the investigated temperature region. The control over an advanced set of process parameters as well as the combination of two well-established synthesis routes, namely, CVD and LA, allows accurate and stable control of the mean diameter. Furthermore, a study of the mean diameter changes with temperature, pressure, and flow rate provides valuable experimental input for elaborating a consistent microscopic description of the formation mechanism of SWCNT.
The targets for LA are mounted in the center of a 40 mm wide quartz reactor tube with an optional 16 mm wide constriction around the center. The reactor tube passes through the hot zone of a horizontal oven with possible temperatures up to 1100 °C. A Q-switched high-power Nd:YAG laser (2.5 GW per pulse, pulse width 8 ns)12 is used at its fundamental wavelength of 1064 nm and focused to a spot of ca. 0.5 cm2. The laser beam is in line with the applied gas stream in the reaction zone. Type, pressure, and flow rate of the hydrocarbons (serving as the carbon feedstock), the range of investigated synthesis temperatures, and the different targets are provided in Table 1. The targets were mounted to a molybdenum holder, which in turn is attached to a water-cooled copper coldfinger. The soot is collected on this coldfinger. A Bruker IFS100 Fourier transform−Raman spectrometer (excitation wavelength λ = 1064 nm) and a Renishaw In Via Raman microscope spectrometer (λ = 514 and 785 nm) were used for sample characterization. The samples were prepared by dispersing the products in acetone and then thoroughly sonicating them to ensure homogeneity. Subsequently, the solution was dropped onto heated plates of Al, resulting in a thin film of homogeneous material. For TEM imaging, the solution was dropped on a KBr single crystal and the resulting film was floated off in distilled water and collected on a standard Cu TEM grid. A FEI Tecnai F30 300 kV transmission electron microscope (TEM) was used to assess both the yield and quality of the samples on a local scale. The samples were also characterized by optical absorption spectroscopy (OAS) on a Bruker IFS113V/88 FT-spectrometer.
Optical techniques such as OAS and resonant Raman spectroscopy are popular tools to characterize SWCNT because they allow a broad range of information on yield, diameter, diameter distribution and quality to be obtained. A typical SWCNT Raman spectrum contains the tangential modes viz. G-mode, which is reminiscent of in-plane modes in graphite, the disorder induced D-band, and the radial breathing modes (RBM). The diameters of the SWCNT are related to the RBM modes via
14, where ν is the Raman shift and d is the diameter of the SWCNT. To determine the yield of the SWCNT, in general, we can use the normalized intensity of the G-mode13. Additionally, the ratio of the G-mode intensity to the D-mode intensity (G:D) provides a means to assess the quality of the SWCNT. OAS also provides information on SWCNT yield and on their diameter distribution via analysis of the characteristic absorption peaks15, 16. In these studies we use the first absorption peak, E11s, (which is corrected for its excitonic nature)16 to determine the SWCNT yield, the mean diameter and the diameter distribution16. The yield corresponds to the area of the background-strapped E11s peak, whereas information on the mean diameter and diameter distribution can be obtained from the fit parameters of a Gaussian distribution15. The bottom left panel in Figure 1 shows that yield determination from OAS and Raman spectroscopy agree well. Thus the use of the G-mode intensity to determine the yield is particularly useful when the SWCNT yield is insufficient to be accurately determined via OAS. However, because of the resonance enhancement of the G-line, caution should be exercised in using the G-mode on its own for yield determination. However, as we show, in our case, it is useful to show trends. These trends are supported by TEM and OAS data.
We investigated a broad variety of different carbon feedstocks; methane, methanol, ethanol, propane, cyclohexane, toluene, and benzylamine with iron as a catalyst17. Propane, cyclohexane, and ethanol18, 19 were also tested with a nickel−cobalt−molybdenum (50:40:10 wt %) composite catalyst. Figure 2 shows the yield as a function of temperature for selected combinations of catalyst and carbon feedstock. For synthesis of these samples, the 16 mm wide reactor tube was used and the flow rate was kept constant at 1 lpm. The partial pressure of methane was set to 50 mbar for all samples shown on the left-hand side (A−C). In A, a composite target consisting of iron and 25 wt % carbon was used instead of pure iron. In C, pure methane was used instead of argon with 5% methane. The optimum temperature is always 900 °C. The extra carbon in A strongly increases the available carbon in the reaction process. The yield is highest in A relative to B and C. The upper two panels on the right-hand column (D, E) compare the yield of samples synthesized using a pure iron target with different carbon feedstocks. The samples in D and E were grown from cyclohexane and propane, respectively. The pressure of cyclohexane was chosen to be 40 mbar and the pressure of propane was set to 80 mbar to provide the same total amount of carbon per unit time. Although the maximum yield is quite comparable in these two cases, the optimum temperature is lower for both cyclohexane and propane as compared to methane. Finally, in part F of Figure 2, we show the yield vs temperature profile of an alloyed metal target (NiCoMo 50:40:10 wt %) under the same conditions as the pure iron target in E. The temperature profile is significantly smoother in the case of the compound target. The TEM images in Figure 2 belong to samples grown at the respective optimized temperatures (marked with an asterisk on panels D−F). The typical overall morphology of all high yield SWCNT material consists of bundles of nanotubes as well as a fair fraction of carbon coated catalyst material. However, for typical low-yield samples like in B and C (not shown), more individual nanotubes are observed than bundles. Figure 3 shows the yield and quality versus temperature profiles for propane and ethanol as carbon sources and NiCoMo (50:40:10 wt %) as the catalyst. Solid symbols are used for the wider 40 mm reaction tube, and open symbols denote the use of the 16 mm wide constriction. The optimum temperature for maximum yield is higher for ethanol than for propane. The wider reaction tube clearly shows improved characteristics regarding the overall yield as well as the quality of the nanotubes. Resonance Raman spectroscopy, as depicted for propane and selected temperatures in Figure 4, reveals a clear increase in the SWCNT mean diameter with increasing temperature (decrease in the RBM frequencies). The two columns on the left-hand side (A, 800 °C; B, 1000 °C) show the RBM at different laser excitation energies. The two columns on the right-hand side show the RBM and high-frequency region (G-line, D-line) for various temperatures at infrared excitation (1.16 eV). The shift in mean diameter can easily be tracked via the relative ratio of the RBM responses at 163 and 265 cm−1 (I163:I265) in the spectra at 1.16 eV. This specific excitation energy is chosen to compare the SWCNT diameters because its resonances cover the widest spread of diameters. OAS studies as shown in Figure 1 for SWCNT synthesized from ethanol at different temperatures show the same increase in the mean diameter with increasing reaction temperature. In addition, OAS reveals a concomitant increase in the width of the diameter distribution (Figure 1, left side). The well resolved fine structures in the E11s peak (Figure 1 right side) may be a signature of the predominant diameters16 or stem from a chirality dependence of the matrix elements. Neither can be excluded, because a specific (n, m) assignment requires a detailed analysis via Raman spectroscopy using multiple-frequency excitation lasers or fluorescence spectroscopy. Figure 6 shows the mean diameter (left panels) and yield (right panels) as a function of either pressure (upper panels) or flow rate (lower panels) for SWCNT synthesized from propane at 800 °C. An increasing pressure decreases the mean diameter and diminishes the yield at the same time. An increasing flow rate increases the mean diameter up to a saturation value, and the yield runs through a clear maximum just before the mean diameter starts to saturate.
| | Figure 5. Proposed nucleation mechanism after ref 21. Right: Overlap of the nucleation window with the precipitating cluster distribution determines the SWCNT diameter distribution. The change in the upper and lower limits arises from changes in the carbon availability and therefore different V:A constraints. Left: With a shift in the nucleation window with increasing carbon availability, the size of the nucleating catalyst clusters is reduced. Clusters with larger size (V:A ratio) are encapsulated by carbon, whereas small clusters cannot precipitate enough carbon to form a nucleating cap.21 |
The various experimental observations on the yield and diameter distribution reported above are brought together in a detailed microscopic model as discussed in the following section.
Going back to Figure 2, we want to explain the differences in yield using the various target-carbon feedstock combinations. These experiments were conducted with the 16 mm diameter tube. It is evident that Ar:5% CH4 does not provide enough carbon to allow high-yield carbon nanotube nucleation and growth and that this can be compensated for by adding carbon to the target. The effect of using pure methane at 50 mbar versus diluted methane (partial pressure 50 mbar, total pressure 1 bar) is not so easy to understand, as the use of 50 mbar pure methane does not change the carbon to metal ratio in the reaction zone at first sight. However, it does change the dynamics of the evaporation plume because a much lower background gas pressure leads to a larger expansion of the plume20, slower cooling and better intermixing with the hydrocarbon gas. This changes both the nucleation point in time and the respective numbers and sizes of catalyst particles, which nucleate SWCNT. It was recently suggested that not all catalyst particles can nucleate SWCNT21. In brief, particles that are too small cannot take up enough carbon to form a stable SWCNT cap because of the low volume to surface area ratio (V:A ratio). Particles that are too big take up too much carbon and become encapsulated by carbon before they can nucleate SWCNT. This imposes a V:A ratio constraint on the catalyst particles available for nucleation, thus defining a nucleation window. SWCNT nucleation can occur only when the nucleation window and the precipitating catalyst particle distribution overlap21. The compounds cyclohexane and propane have relatively higher molar carbon content. In addition, they also have C−C bonds, which are energetically easier to break as compared to the C−H bonds in methane. A radical induced decomposition is also likely, as the carbon to hydrogen ratio is much higher in these compounds compared to methane22. The lower optimum temperatures for these compounds can be explained as a reduction of available carbon at high temperatures because of self-pyrolysis and self-polymerization and condensation to light and heavy oils23, 24. Propane proved to be a little bit more stable than cyclohexane. The increased yield−temperature curve width when using a NiCoMo catalyst (Figure 2F) can be explained by a wider spread of particle compositions21 and thus slightly different eutectic temperatures. For each of these temperatures, a part of the precipitating catalyst distribution lies in the nucleation window, and these catalyst particles are able to nucleate SWCNT21. Figure 3 shows the yield and quality results for two different carbon feedstocks, propane and ethanol. It is easy to see that the yield is much higher when using the 40 mm diameter reactor tube, in contrast to what one might expect25. This can be explained by the larger volume available for plume expansion and, more importantly, by the slower gas velocity in the reaction zone. The carbon required for nanotube growth and nucleation is not directly available, as found in LA, but has to be obtained from the decomposition of hydrocarbons, which requires some time. The optimum temperature does not change with the width of the reactor, pointing to a strict dependence on thermal aspects of catalyst particle formation21. The optimum temperature using ethanol (with NiCoMo) as a carbon feedstock is restored to 900 °C, comparing well with the work of Nishide et al.10 for a pure Co target. This is because self-polymerization and self-pyrolysis do not play a significant role for ethanol below 1000 °C, an explanation being OH radical formation26. OH radicals are also known to etch amorphous carbon species27, resulting in a better G:D ratio for the ethanol-synthesized SWCNT (Figure 3). OAS and Raman studies show that the width of the diameter distribution and the mean diameter increase with temperature (Figure 1). Their evolution with temperature corresponds well with similar LA studies at these temperatures;12, 21, 28 however, the width of the diameter distribution is expected to saturate at roughly twice the width of a comparative LA sample. The shift in mean diameter can be explained by a shift in the nucleation window against the precipitation window and therefore a different overlap21. At higher temperatures, the cooling is slower because of a smaller temperature gradient and the resulting SWCNT diameters will increase because the catalysts particles have more time to coalesce/ripen. If the overlap is large, the width of the distribution is large21. This model is also applicable to the changes in the mean diameter upon a change of pressure and flow rate. These changes are not as pronounced as compared to the changes with temperature. Figure 6 shows these results, where a decrease of the mean diameter (small I163:I265) with increasing pressure and decreasing flow is evident. On the right side of Figure 6, the yield of the respective samples is shown. Increased pressure means a higher carbon:metal ratio because of a higher thermal conductivity and higher precursor density at the point of nucleation. The higher carbon availability shifts the nucleation window to smaller catalyst particles, as bigger ones are encapsulated because of an excess of carbon (Figure 5, right panel). The yield is found to decrease with increasing carbon content, which is contrary to what one might expect. This is because the overlap of the nucleation window and precipitation window is shifted toward the tail of the precipitating catalyst cluster distribution (Figure 5, left panel). Concomitant with these results, a decrease in SWCNT diameters with increasing pressure has also been reported for ferrocene-based floating catalyst CVD29. The flow-rate dependence of diameter and yield follows a similar argumentation. Low gas velocities mean that the residence time of the catalyst particles in the active zone is prolonged, leading to a higher carbon:metal ratio because there is more time available for precursor decomposition, resulting in CNT with smaller diameters. In turn, the mean diameter increases with increasing gas flow. This mechanism is known from LA21. However, at very high gas velocities, it is expected that the gas is not fully heated (lower gas temperature). This effect counteracts the increase in diameter and leads to a stabilization of diameters in the high-flow region investigated. The decrease in yield at too high flow rates can again be explained by the fact that the overlap of the two respective windows diminishes.
Hence the pulsed-laser-assisted CVD system inherits the dynamics of catalyst particle formation of LA and extends them by an independently tunable carbon feedstock. This hybrid synthesis approach leads to SWCNT yields smaller than state-of-the-art LA but allows synthesis at lower temperatures. Further, for our system, no growth promoters like sulfur31 or water vapor32 were used. The relative purity of samples synthesized with our system as determined with the method developed by Itkis et al.30 (NiCoMo 50:40:10 wt % in 50 mbar ethanol) is about 12%, which is approximately 25% of those obtained from LA with a total of 10 wt % catalyst. Further, this hybrid route allows lower synthesis temperatures than in LA and affords controlled mean diameters and narrow diameter distributions, which other CVD routes do not offer. A study of the effects of the laser power and power density as well as the use of growth promoters known from CVD like sulfur31 or water32 would provide further insight into the nucleation and growth processes and clarify the role of these remaining optimization parameters. If we assume this to be known, a scale-up is in principle possible in the same manner as LA may be scaled up. This route has a further advantage, in that the use of pure metal targets have a much longer lifetime as compared to mixed carbon:metal targets as used in conventional LA.
We have demonstrated that it is possible to synthesize SWCNT by laser evaporation in a hydrocarbon gas atmosphere. The hybrid route inherits the positive aspects of LA and CVD systems. SWCNT with narrow diameter distributions, as found in LA, are obtained. The use of a hydrocarbon feedstock allows in principle for continuous synthesis through a constant supply of fresh catalyst particles. The obtained data show a correlation between the diameter of the nanotubes and the carbon availability, enabling the mean diameter of the SWCNT to be fine-tuned via the background pressure and flow rate of the precursor. The results are fully explained by our previously suggested volume to surface area constraint growth model21, which we now extend to incorporate the effects of carbon availability on the SWCNT diameter. Our findings provide new insight into understanding the nucleation processes governing the synthesis of carbon nanotubes. Furthermore, the use of this hybrid synthesis route provides a means to bridge the gap between CVD and laser ablation, enabling improved yield and diameter control as compared to conventional supported catalyst CVD in this temperature region.
The work was partly conducted under and financed by the “Pakt für Forschung und Innovation”. C.K. acknowledges the IMPRS Dynamical Processes in Atoms, Molecules, and Solids for financial support. M.H.R. acknowledges DFG PI440/3-/4 for financial support. The authors are grateful to R. Hübel, R. Schönfelder, and S. Leger for technical assistance.
* Corresponding author. E-mail: markus.loeffler@ifw-dresden.de.
1. Iijima, S.; Ichihashi, T. Nature 1993, 363, 603.
2. Saito, R.; Dresselhaus, G.; Dresselhaus, M. S. Physical Properties of Carbon Nanotubes; Imperial College Press: London, 1998
3. Journet, C.; Maser, W. K.; Bernier, P.; Loiseau de la Chapelle, M.; Lefrant, S.; Deniard, P.; Lerr, R.; Fischer, J. E. Nature 1997, 388, 756–758.
4. Maser, W. K.; Munoz, E.; Benito, A. M.; Martinez, M. T.; de la Fuente, G. F.; Maniette, Y.; Anglaret, E.; Sauvajol, J. L. Chem. Phys. Lett. 1998, 292, 587–593.
5. Colomer, J. F.; Bister, G.; Willems, I.; Kónya, Z.; Fonseca, A.; Van Tendeloo, G.; Nagy, J. B. Chem. Commun. 1999, 14, 1343–1344.
6. Kim, S. H.; Zachariah, M. R. Mater. Lett. 2007, 61, 2079–2083.
7. Schäffel, F.; Kramberger, C.; Rümmeli, M. H.; Grimm, D.; Mohn, E.; Gemming, T.; Pichler, T.; Rellinghaus, B.; Büchner, B.; Schultz, L. Chem. Mater. 2007, 19, 5006–5009.
8. Alexandrescu, R.; Crunteanu, A.; Morjan, R. E.; Morjan, I.; Rohmund, F.; Falk, L.K.L.; Ledoux, G.; Husiken, F. Infrared Phys. Technol. 2003, 44, 43–50.
9. Rümmeli, M. H.; Kramberger, C.; Löffler, M.; Kalbác, M.; Hübers, H-W.; Grüneis, A.; Barreiro, A.; Grimm, D.; Ayala, P.; Gemming, T.; Schäffel, F.; Dunsch, L.; Büchner, B.; Pichler, T. Nanotechnology 2006, 17, 5469–5473.
10. Nishide, D.; Kataura, H.; Suzuki, S.; Okubo, S.; Achiba, Y. Chem. Phys. Lett. 2004, 392, 309–313.
11. Cheng, H. M.; Li, F.; Su, G.; Pan, H. Y.; He, L. L.; Sun, X.; Dresselhaus, M. S. App. Phys. Lett. 1998, 72 (25), 3282.
12. Rümmeli, M. H.; Löffler, M.; Kramberger, C.; Ferenc, S.; Ferenc, F.; Jost, O.; Schönfelder, R.; Grüneis, A.; Gemming, T.; Pompe, W.; Büchner, B.; Pichler, T. J. Phys. Chem. C 2007, 111, 4094–8.
13. Williams, K. A.; Tachibana, M.; Allen, J. L.; Grigorian, L.; Cheng, S.-C.; Fang, S. L.; Sumanasekera, G. U.; Loper, A. L.; Williams, J. H.; Eklund, P. C. Chem. Phys. Lett. 1999, 310, 31–7.
14. Kramberger, C.; Pfeiffer, R.; Kuzmany, H.; Zólyomi, V.; Kürti, J. Phys. Rev. B 2003, 68, 235404.
15. Jost, O.; Gorbunov, A. A.; Pompe, W. Appl. Phys. Lett. 1999, 55 (15), 2217–2219.
16. Liu, X.; Pichler, T.; Knupfer, M.; Golden, M. S.; Fink, J.; Kataura, H.; Achiba, Y Phys. Rev. B 2002, 66, 045411.
17. Grüneis, A.; Kramberger, C.; Grimm, D.; Gemming, T.; Rümmeli, M. H.; Barreiro, A.; Ayala, P.; Pichler, T.; Schaman, C.; Kuzmany, H.; Schumann, J.; Büchner, B. Chem. Phys. Lett. 2006, 425, 301–305.
18. Mizuno, K.; Hata, K.; Saito, T.; Ohshima, S.; Yumura, M.; Iijima, S. J. Phys. Chem. B 2005, 109, 2632–2637.
19. Inoue, S.; Kikuchi, Y. Chem. Phys. Lett. 2005, 410, 209–12.
20. Harilal, S. S.; Bindhu, C. V.; Tillack, M. S.; Najmabadi, F.; Gaeris, A. C. J. Appl. Phys. 2003, 93 (5), 2381–2388.
21. Rümmeli, M. H.; Kramberger, C.; Löffler, M.; Jost, O.; Bystrzejewski, M.; Grüneis, A.; Gemming, T.; Pompe, W.; Büchner, B.; Pichler, T. J. Phys.Chem. B 2007, 111, 8234–8241.
22. Reilly, P. T. A.; Whitten, W. B. Carbon 2006, 44, 1653–60.
23. Descamps, C.; Vignoles, L.; Féron, O.; Langlais, F.; Lavenac, J. J. Electrochem. Soc. 2001, 148 (10), C695–C708.
24. Rodgers, R. P.; Reilly, P. T. A.; Whitten, W. B.; Ramsey, J. M. Carbon 2003, 41, 1469–75.
25. Scott, C. D.; Arepalli, S.; Nikolaev, P.; Smalley, R. E. Appl. Phys. A: Mater. Sci. Process. 2001, 72, 573–80.
26. Park, J.; Zhu, R. S.; Lin, M. C. J. Chem. Phys. 2002, 117 (7), 3224–3231.
27. Maruyama, S.; Kojima, R.; Miyauchi, Y.; Chiashi, S.; Kohno, M. Chem. Phys. Lett. 2002, 360, 229–234.
28. Bandow, S.; Asaka, S.; Saito, Y.; Rao, A.; Grigorian, L.; Richter, E.; Eklund, P. C. Phys. Rev. Lett. 1998, 80, 3779.
29. Barreiro, A.; Hampel, S.; Rümmeli, M. H.; Kramberger, C.; Grüneis, A.; Biedermann, K.; Leonhardt, A.; Gemming, T.; Büchner, B.; Bachtold, A.; Pichler, T. J. Phys. Chem. B 2006, 110, 20973–20977.
30. Itkis, M. E.; Perea, D. E.; Jung, R.; Sandip, N.; Haddon, R. C. J. Am. Chem. Soc. 2005, 127, 3439–3448.
31. Barreiro, A.; Kramberger, C.; Rümmeli, M. H.; Grüneis, A.; Grimm, D.; Hampel, S.; Gemming, T.; Büchner, B.; Pichler, T. Carbon 2007, 45, 55–61.
32. Hata, K.; Futaba, D. N.; Mizuno, K.; Namai, T.; Yumura, M.; Iijima, S. Science 2004, 306, 1362–1364.
| Table 1. Overview of the Used Catalysts and Carbon Feedstock As Well As the Respective System Parameters Explored | ||||||||||||||||||||||||||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|