search for


Collisional Ion Energy Transfer Model for Understanding Wrinkles on Poly (dimethylsiloxane) Surfaces after Oxygen Ion Beam Irradiations
Applied Science and Convergence Technology 2018;27:130-134
Published online November 30, 2018;
© 2018 The Korean Vacuum Society.

Seunghun Lee*, Junyoung Yang, Sunghoon Jung, Do-Geun Kim, and Eunyeon Byeon

Advanced Nano-Surface Department, Korea Institute of Materials Science, Gyeongsangnam-do 51510, Republic of Korea
Correspondence to: *Corresponding author: E-mail:
Received September 28, 2018; Revised October 18, 2018; Accepted October 25, 2018.

We introduce a collisional ion energy transfer model to understand wrinkle wavelengths on poly (dimethylsiloxane) surfaces after ion beam irradiations. The collisional ion energy transfer model uses a Monte Carlo simulation and a full molecular method to predict the thickness of a hard skin layer produced by ion beam irradiations on poly (dimethylsiloxane) surfaces. Oxygen ion beam irradiations were conducted to verify the theoretical model. To make wrinkles with wavelength in the range of 0.6–1.0 μm, the model derived appropriate ion energies of 700–1,000 eV. Oxygen ion beams with the average energy of 707, 800, 901, and 1,043 eV were exposed on the polymer surfaces, which showed randomly distributed wrinkle patterns with wavelengths of 0.63 ± 0.25, 0.84 ± 0.30, 0.93 ± 0.32, 1.27 ± 0.42 μm, respectively. The wavelengths observed through the experiments showed consistency with the collisional ion energy transfer model.

Keywords : Ion beam, Polymer, Wrinkle, Collision, Cross-linking
I. Introduction

Nanostructure surfaces on polymers are of considerable interest in nanotechnologies that include tunable optics [1], surface enhanced Raman spectroscopy [2], interlayer adhesion control [3], flexible electronics [4], and biomimetic micro fingerprints for anti-counterfeiting [5]. Wrinkle nanostructures on soft polymers have been intensively studied as wrinkle surface patterning is used in several applications such as stretchable hydrophobic substrates and security applications [5,6]. Wrinkle topologies have been fabricated spontaneously by means of mask-free processes, such as thin film deposition [7], plasma treatment [8], solvent-driven oxidation [9], thermal annealing [10], and ion beam (IB) irradiation [11]. These methods commonly form a hard skin layer that generates local moduli difference causing wrinkles [7]. In particular, the plasma treatment and ion beam irradiation using an ion bombardment is suitable to make the hard skin layer. Because ion-polymer surface interactions change the binding status within limited thickness below several tens nm. In poly(dimethylsiloxane) (PDMS), the IB irradiation results in the formation of hard skin layer within 1–10 nm, which induces wrinkles with wavelength from hundreds nm to several μm because the wrinkle wavelength is a function of the thickness of hard skin layer [12].

Uniform wrinkle formation on PDMS substrates is important for industrial applications. Localized ion irradiation in focused ion beams showed the inhomogeneity of wrinkle patterns due to strain localization [11,1315]. Large area IB irradiation with uniform ion flux and energy is necessary to make homogenous wrinkle patterns. A linear IB irradiation is possible to treat flexible substrates with a web speed of 1–10 m/min when linear ion beam source and roll-to-roll treatment system are utilized [16]. The linear IB source emits linearly uniform gas ion beams such as argon, oxygen, nitrogen with a beam width of 1.5 m and an ion flux uniformity below 5 %. Thus, homogenous wrinkle formation is available in in-line or roll-to-roll processes without inhomogeneity problem due to locally induced strain after IB exposure. And the linear IB process is an attractive candidate to fabricate wrinkle surfaces with high throughput, relatively low cost, and clean process comparing to spin coating, electron beam lithography, optical lithography, and wet chemical processes.

The design of wrinkle wavelength is also important for the various applications. The parameters of IB irradiation such as ion energy and incident angle determine a depth in which the hard skin layer is formed. Recently, ion beam irradiations on PDMS elastomer have been reported to produce wrinkle structures as a function of the IB irradiation parameters [12,17,18]. The wrinkle wavelength has been modulated successfully in range of 0.75–1.95 μm by changing IB energy from 800 eV to 2000 eV [18].

Also the wrinkle wavelength on PDMS have been controlled from 0.93 μm to 2.03 μm by adjusting the incident angle of impinging ions from 15° to 90°. The ion incident angle changes cross-linking density, which affects the Young’s modulus of the PDMS surface [12].

Understanding a wrinkle formation based on ion-polymer surface interactions is imperative to apply the IB technique to the wrinkle surface patterning. Many reports have introduced IB induced wrinkle patterns on PDMS with descriptions based on mechanical and physiochemical analysis, such as modulus measurement and X-ray photoelectron spectroscopy [11,12,17]. Though the post-experiment investigations explained wrinkle formation in detail, a theoretical model which predicts wrinkle formation is required to design a wrinkle wavelength for various purposes. In this work, we demonstrated quantitative interpretations of interactions between impinging ions and PDMS surfaces using a collisional ion energy transfer (CIET) model. The CIET model calculates a depth which is effective to form a the hard skin layer using Monte-Carlo simulations and a full molecular method [19,20]. The wavelength expected by the CIET model were compared with the wrinkle wavelength obtained after linear oxygen IB irradiation on PDMS surfaces.

II. Experimental details

PDMS precursor mixtures were prepared by a combination of Sylgard 184 base and hardener (Sylgard 184, Dow Corning Corporation, USA) in the mass ratio of 10:1, by mixing in a Petri dish for 10 min. The prepared mixtures were placed in vacuum chamber, until no bubbles remained in the bulk of the mixtures. Glass substrates cleaned by a standard method were used for PDMS spin coating. The precursors mixtures were spin coated onto 25 × 25 mm2 glass for 40 s at 500 rpm. This sample preparation procedure is identical to reference to use Young’s modulus and Poisson’s values of bulk PDMS in wavelength calculation [12]. The PDMS coatings were then annealed by heating on hotplate for 120 min at 65 °C. PDMS samples on the glass substrates were exposed to oxygen IB generated by closed drift ion sources [16]. The ion source emits linear IB (width: 400 mm) with the ion current density of 0.1 ± 0.02 mA/cm2 and the ion dose irradiated on PDMS samples were 4.5 ± 0.5 × 1016cm−2. The anode voltage controls the average ion energy: 707, 800, 901, and 1,043 eV. A retarding potential analyzer (RPA) was used to measure the ion energy distribution function (IEDF) of ions irradiated on PDMS samples [21]. The IB irradiation was conducted in vacuum environment of 0.5–1 mTorr. The surface image of the wrinkle patterns was observed by field emission scanning electron microscopy (JSM-6700F, JEOL). The PDMS samples were prepared by 10 nm silver film coating using magnetron sputtering. The acceleration voltage was kept at 5 kV. The SEM images were subjected to 2D FFT analysis to characterize the wrinkle wavenumbers using commercial software (SPIPTM).

III. Results and discussion

Figure 1 demonstrates scanning electron microscope images and wavenumber distributions analyzed by two-dimensional fast Fourier transform power density spectra. As the ion energy increases the PDMS surfaces showed a wrinkle width broadening, which corresponds to the decrease of wrinkle wavenumber: 1.60 ± 0.56, 1.19 ± 0.38, 1.08 ± 0.34, and 0.79 ± 0.24 μm−1 for each ion energy represented in Fig. 1(e). The deviations are the full width at half maximum of the wavenumber distribution. The wrinkle wavelengths that correspond to the wavenumbers are 0.63 ± 0.25, 0.84 ± 0.30, 0.93 ± 0.32, 1.27 ± 0.42 μm, respectively.

IB irradiations make a hard skin layer that determines the wrinkle wavelength. The competition between the bending energy of the hard skin layer and the stretching energy of the soft supporting PDMS determines the wrinkle wavelength [17]. Eq. (1) including the Young’s modulus and Poisson’s ratio describes the competing behaviors:


where λ is the wrinkle wavelength, h is the thickness of the hard skin layer; EHS and EPDMS are the Young’s moduli of the hard skin layer and the PDMS substrate, respectively; and νHS and νPDMS are the Poisson’s ratios of the hard skin layer and the PDMS substrate, respectively. The Young’s modulus of bulk PDMS has been controlled simply by annealing temperature [17]. In this work, the PDMS samples were annealed at 65 °C for the Young’s modulus of 1.5MPa, which simplified Eq. (1) to k=1/147 h using several known values, EHS ≈70GPa, νHS ≈0.2, and νPDMS ≈0.48 [12]. In this procedure, we assumed that the oxygen IB treated PDMS surfaces in Fig. 1 have same EHS and νHS values to predict wrinkle wavelength using the CIET model and Eq. (1).

Collisional behaviors between ions and polymers could explain the thickness of the hard skin layer. Irradiated ions penetrate polymer surfaces to a definite depth accompanying energy transfers by cascade collisions. Incident ions collide with polymer atoms, and make energetic recoils. Ions and recoils lose their energy to make displacements of polymer atoms: nuclear stopping. Displacement occurs when the transferred energy is larger than the threshold energy for displacement. Polymers have 2–10 eV displacement energy, which is smaller than that of semiconductors (~15 eV) and of metal (25–40 eV) [22]. The displacement results in vacancy, replacement, and interstitial productions. If the transferred energy is not sufficient to generate displacement, phonon vibration consumes the recoil energy thermally [22]. The nuclear stopping is a driving force for scission and cross-linking initiating a formation of hard skin layers. Thus, quantitative estimation of the energy transfer as a function of depth could explain the thickness of the hard skin layer after IB irradiations.

The energy transfer as a function of depth was calculated by Monte Carlo based simulation. The stopping and range of ions in matter (SRIM) was used to predict collisional ion energy transfer rates to PDMS components: C, H, O, Si [22]. In the SRIM calculation, the conventional PDMS formulation (C2H6OSi)n was adapted (density: 0.868 g/cm3). The calculation type was “detailed calculation with full damage cascades”. Monoatomic IB irradiation, such as He+ and Ar+ ions, can be precisely calculated, because SRIM applies the binary collision approach. Meanwhile, molecular IBs, such as O2+ and N2+, cannot be directly treated by binary collision. However, the penetration and energy loss of the molecular IBs can be calculated by using a theory of molecular IB implantations: the full molecular method [20]. This assumes that when it enters the target, the molecular ion dissipates immediately as shown in Fig. 1. This is true if the ion energy is much higher than the dissociation energy of molecular ions, which are of the order of 10 eV/atom [20]. For ion energies in several hundreds eV range, this requirement for spontaneous dissociations is certainly met. Thus, an impacting molecular ion can be considered as two monoatomic ions. The incident energy is divided by the mass ratios represented in Eq. (2):


where, Ei is the monoatomic ion’s energy of i species, E is the energy of the molecular IB; Mi and Mj are the mass of each species, respectively; and Nj is the atom number of j species. In the molecular oxygen IB, the CIET model substitutes a molecular oxygen IB (O2+) by two monoatomic oxygen atoms (O) with half of the molecular ion’s energy in PDMS medium. Actually, oxygen plasma generated by low temperature DC discharge consists of mostly molecular oxygen ions (O2+) [23,24]. The few percent of monoatomic oxygen ions having a whole acceleration energy can be ignored in this experiment conditions. Thus, the CIET model assumes that molecular oxygen ions are dissociated to two monoatomic oxygen atoms, which have a half of molecular oxygen ion’s energy.

The oxygen atoms transfer energy to recoils by collisions in Fig. 2. When the transferred energy to recoil is larger than the displacement energy of polymer components, the recoils induce vacancy, replacement, and interstitial to a certain depth, in which scission and cross-linking are occurred. When the transferred energy to recoil is lower than the displacement energy, recoils generate phonon that releases transferred energy thermally. The phonon in amorphous phase cannot deliver heat efficiently and effective heat transfer length is lower than 10 nm. The temperature increases locally in the ion energy transfer depth and the generated heat transfer with a low thermal conductivity of 0.1~0.2W/m-k [25]. Therefore, it is necessary to consider the depth at which the transferred energy effectively change the binding of polymer surface due to displacement and heat generation. SRIM code calculates the collisional energy transfer rate to PDMS recoils (Erecoil) of the monoatomic oxygen atom, which is used to calculate the effective energy transfer depth (heff) in Fig. 1 creating hard skin layer after oxygen IB irradiations.

The Erecoil was calculated by SRIM to predict h in Eq. (1). Ions transfer their incident energy to the PDMS components: C, H, O, Si. Figure 3(a) shows the Erecoil when the energy of atomic oxygen ion is 400 eV at normal injection. Carbon and hydrogen atoms absorb ion energy dominantly at PDMS surface. Note that Methyl group (CH3) dissociation is a main reaction to make the hard skin layer because the hard skin layer has no methyl group. For example, Sylgard 184 base PDMS is composed of dimethylsiloxane oligomers with vinyl-terminated end groups [26]. Pristine PDMS has shown strong infrared (IR) absorptions at 2,850 cm−1 (CH2 symmetric stretch), 2,910 cm−1 (CH3 symmetric stretch), and 2,960 cm.1 (CH3 asymmetric stretch). By oxygen plasma treatment, the methyl (CH3) groups have been transformed into a hydroxyl-terminated surface with negligible residual methyl groups due to hydrogen dissociations [27]. Thus, we assume that the profile of Erecoil to hydrogen atoms relates to the h in Eq. (1).

Figure 3(b) shows the Erecoil of an atomic oxygen to hydrogen atoms in PDMS medium as a function of impinging energy. The maximum energy transfer depth increased from 70 Å to 100 Å as the ion energy is increased from 350 to 500 eV, which corresponds to half of the molecular oxygen ion’s energy, adapting the full molecular method represented in Eq. (2). Considering the hydrogen dissociation energy of a methyl chain in PDMS, monoatomic oxygen transfers enough energy to dissociate a hydrogen atom until a certain depth. We define the certain depth as an effective energy transfer depth (heff) that satisfies the following equation,


where Ediss is the hydrogen dissociation energy, heff is the effective energy transfer depth, Erecoil is the energy transfer rate function in Fig. 3(b), and H is the depth from the PDMS surface. Eq. (3) means that the heff is the thickness when residual energy of penetrating oxygen atom is equal to the energy of methyl dissociation. In other words, heff is the thickness at which a hydrogen dissociation of a methyl chain could be possible by energetic hydrogen recoils. For example, the hydrogen dissociation energy of a methyl chain is 473 kJ/mol, which is same to Ediss = 4.9 eV/atom [28]. The heff of H-CH2 bonding as a function of the monoatomic oxygen ion energy from 350 eV to 500 eV is plotted in Fig. 4. The heff was obtained by integrating the plots of Erecoil(H) in Fig. 3(b) within a range, which satisfies Eq. (3). It shows that monoatomic oxygen ions with energy of 350–500 eV deliver sufficient energy to dissociate H-CH2 bonding to the heff of 45–65 Å. The dissociation depth is linearly increased due to that the Erecoil of hydrogen atom is also linearly increased when monoatomic oxygen ion energy increases. Substituting h in Eq. (1) to the heff in Fig. 4, the wrinkle wavelength expected by heff from the CIET model was 0.65–0.96 μm. The wrinkle wavelength observed in PDMS surfaces after oxygen IB irradiation in Fig. 1 showed consistency with the CIET model within the FWHM of wrinkle wavelength distribution. This shows the possibility that the amount of ion energy transferred to the PDMS recoils can explain the thickness of the hard skin layer using heff based on the assumption that polymer bonds are dissociated when the binding energy is less than the energy delivered to the polymer components.

The prediction of wrinkle wavelength showed deviation in the oxygen ion energy of 500 eV. The PDMS specimens consist of a chain of Methyl (CH3) groups (Ediss = 3.5–4.9 eV/atom), Si-C (Ediss = 4.5 eV/atom), Si-O (Ediss = 8.3 eV/atom) [28]. As the energy of oxygen ion increases, the opportunity that ion beam irradiation initiates the dissociation of bonds that have higher Ediss than that of Methyl groups increases, and it forms different deformed layer. Though the heff in Fig. (4), which is predicted by the correlation between the H-CH2 binding energy (Ediss = 4.9 eV/atom) and the Erecoil to hydrogen atoms, showed reasonable wrinkle wavelength up to 400 eV, another dissociation reaction such as Si-C, Si-O in PDMS should be included for precise wrinkle prediction at the higher ion energy above 400 eV.

IV. Conclusions

Wrinkle wavelength predicted by the proposed collisional ion energy transfer model showed analogous values with the wrinkle wavelength observed after oxygen IB irradiations on PDMS surfaces. In the collisional ion energy transfer model, the effective energy transfer depth (heff) for H-CH2 bond dissociation is only considered to predict the thickness of hard skin layer. This approach showed consistency between the theoretical model and experimental result until a certain ion energy. When the ion energy increases and it causes double or triple dissociation processes, the collisional ion energy transfer model needs to deliberate the complex dissociation and cross-linking processes. In future work, improved CIET model including double or triple bond dissociation processes will be investigated for more precise expectation of the wrinkle wavelength after IB irradiation.


This work was supported by the National Research Foundation of Korea (NRF) funded by the Korean government (MSIP: Ministry of Science, ICT and Future Planning) (No. NRF-2017M2A2A6A01019037, and NRF-2017M1A3A4A03025130).

Fig. 1. Surface images at the ion energy of (a) 707 eV, (b) 800 eV, (c) 901 eV, and (d) 1,043 eV, and (e) wavenumber distribution functions for each PDMS surfaces analyzed by 2D FFT power density spectra.

Fig. 2. Schematic of collisional ion energy transfer model with full molecular method: oxygen molecule ion beam dissipated immediately to two monoatomic oxygen atoms at surface and modify PDMS surface within effective energy transfer depth.

Fig. 3. (a) Energy transfer rate to PDMS recoils (C, H, O, Si) per unit depth by an atomic oxygen ion having an energy of 400 eV. (b) Energy transfer rate to hydrogen atoms per unit depth by an atomic oxygen ion with 350, 400, 450, and 500 eV.

Fig. 4. Effective energy transfer depth to dissociate a hydrogen from a methyl group, the corresponding wrinkle wavelength, and observed wavelength on PDMS surfaces as a function of the monoatomic oxygen energy. The hydrogen dissociation energy of 4.9 eV/atom was used in this plot.

  1. Ohzono, T, and Monobe, H (2010). Langmuir. 26, 6127.
    Pubmed CrossRef
  2. Jeon, TY, Park, S, Kim, D, and Kim, S (2015). Adv Funct Mater. 25, 4681.
  3. Chan, EP, Smith, EJ, Hayward, RC, and Crosby, AJ (2008). Adv Mater. 20, 711.
  4. Khang, DY, Jiang, H, Huang, Y, and Rogers, JA (2006). Science. 311, 208.
    Pubmed CrossRef
  5. Bae, HJ, Bae, S, Park, C, Han, S, Kim, J, Kim, LN, Kim, K, Song, S, Park, W, and Kwon, S (2015). Adv Mater. 27, 2083.
    Pubmed CrossRef
  6. Lee, WK, Jung, WB, Nagel, SR, and Odom, TW (2016). Nano Lett. 16, 3774.
    Pubmed CrossRef
  7. Bowden, N, Brittain, S, Evans, AG, Hutchinson, JW, and Whitesides, GM (1998). Nature. 393, 146.
  8. Huntington, MD, Engel, CJ, Hryn, AJ, and Odom, TW (2013). ACS Appl Mater Interfaces. 5, 6438.
    Pubmed CrossRef
  9. Xie, J, Han, X, Zong, C, Ji, H, and Lu, C (2015). Macromolecules. 48, 663.
  10. Park, JY, Chae, HY, Chung, CH, Sim, SJ, Park, J, Lee, HH, and Yoo, PJ (2010). Soft Matter. 6, 677.
  11. Moon, MW, Lee, SH, Sun, JY, Oh, KH, Vaziri, A, and Hutchinson, JW (2007). Proc Natl Acad Sci. 104, 1130.
    Pubmed KoreaMed CrossRef
  12. Jeong, HC, Park, HG, Jung, YH, Lee, JH, Oh, BY, and Seo, DS (2016). Langmuir. 32, 7138.
    Pubmed CrossRef
  13. Chen, YC, and Crosby, AJ (2013). Soft Matter. 9, 43.
  14. Trujillo, V, Kim, J, and Hayward, RC (2008). Soft Matter. 4, 564.
  15. Moon, MW, Lee, SH, Sun, JY, Oh, KH, Vaziri, A, and Hutchinson, JW (2007). Scr Mater. 57, 747.
  16. Lee, S, Byun, EY, Kim, JK, and Kim, DG (2014). Curr Appl Phys. 14, S180.
  17. Park, HG, Jeong, HC, Jung, YH, and Seo, DS (2015). Sci Rep. 5, 12356.
    Pubmed KoreaMed CrossRef
  18. Jeong, HC, Park, HG, Lee, JH, and Seo, DS (2015). ACS Appl Mater Interfaces. 7, 23216.
    Pubmed CrossRef
  19. Ziegler, JF, and Biersack, JP (1985). Treatise on Heavy-Ion Science: Springer, pp. 93-129
  20. Hössinger, A 2000. . PhD Thesis. Technische Univ. Wien.
  21. Böhm, C, and Perrin, J (1993). Rev Sci Instrum. 64, 31.
  22. Ziegler, JF, Ziegler, MD, and Biersack, JP (2010). Nucl Instruments Methods Phys Res Sect B Beam Interact with Mater Atoms. 268, 1818.
  23. Knewstubb, PF, Dawson, PH, and Tickner, AW (1963). J Chem Phys. 38, 1031.
  24. Vancoppenolle, V, Jouan, PY, Wautelet, M, Dauchot, JP, and Hecq, M (1999). J Vac Sci Technol A Vacuum Surfaces Film. 17, 3317.
  25. Van Krevelen, DW, and Te Nijenhuis, K (2009). Properties of Polymers: Elsevier
  26. Lee, JN, Jiang, X, Ryan, D, and Whitesides, GM (2004). Langmuir. 20, 11684.
    Pubmed CrossRef
  27. Ye, H, Gu, Z, and Gracias, DH (2006). Langmuir. 22, 1863.
    Pubmed CrossRef
  28. Cottrell, TL (1958). The Strengths of Chemical Bonds: Academic Press