MF研究者総覧

研究者データベース

詳細
researchmap個人ページ
Last Updated :2024/10/24

研究者情報

マスター

アカウント（マスター）

  • 氏名

    菅原　広剛(スガワラ　ヒロタケ), スガワラ　ヒロタケ

所属（マスター）

  • 情報科学研究院　情報エレクトロニクス部門　集積システム分野

所属（マスター）

  • 情報科学研究院　情報エレクトロニクス部門　集積システム分野

独自項目

syllabus

  • 2021, 集積プロセス学特論, Integrated Material Processings, 修士課程, 情報科学研究科, 固体物理，結晶，逆格子，X線回折，格子振動，自由電子ガスモデル，状態密度，フェルミ準位，エネルギーバンド，有効質量，半導体，ドーピング，pn接合
  • 2021, 集積プロセス学特論, Integrated Material Processing, 修士課程, 情報科学院, 固体物理，結晶，逆格子，X線回折，格子振動，自由電子ガスモデル，状態密度，フェルミ準位，エネルギーバンド，有効質量，半導体，ドーピング，pn接合
  • 2021, 集積プロセス学特論, Integrated Material Processings, 博士後期課程, 情報科学研究科, 固体物理，結晶，逆格子，X線回折，格子振動，自由電子ガスモデル，状態密度，フェルミ準位，エネルギーバンド，有効質量，半導体，ドーピング，pn接合
  • 2021, 集積プロセス学特論, Integrated Material Processing, 博士後期課程, 情報科学院, 固体物理，結晶，逆格子，X線回折，格子振動，自由電子ガスモデル，状態密度，フェルミ準位，エネルギーバンド，有効質量，半導体，ドーピング，pn接合
  • 2021, 電気電子工学実験Ⅰ, Electrical and Electronic Engineering Laboratories I, 学士課程, 工学部, バイポーラトランジスタ，増幅回路，MOSFET，CMOS，演算増幅器，論理回路，順序回路，ディジタル回路
  • 2021, 電気電子工学実験Ⅱ, Electrical and Electronic Engineering Laboratories II, 学士課程, 工学部, 偏光，旋光，回折，屈折，結晶構造，X線回折，半導体物性，ホール効果，電気磁気エネルギー，電力
  • 2021, 電気電子工学実験Ⅲ, Electrical and Electronic Engineering Laboratories III, 学士課程, 工学部, 半導体プロセス，MOSFET，集積回路特性評価，SPICE，回路シミュレーション，アセンブリ言語
  • 2021, 数値解析とシミュレーション基礎, Numerical Analysis and Simulations, 学士課程, 工学部, 二進数、誤差、方程式の数値解法、近似法・補間法、数値積分、モンテカルロ法と乱数

researchmap

プロフィール情報

学位

  • 博士(工学)（1997年06月 北海道大学）
  • 修士(工学)（1992年03月 北海道大学）

プロフィール情報

  • プロフィール

    主要研究テーマは放電プラズマのモデリング，電子輸送基礎理論，ならびに関連するシミュレーション技法開発。電磁界下の電子の微視的挙動からプラズマの構造や応答特性を理解することで新たな制御法や用途の開拓を目指す。
    近年の課題はモンテカルロ法他を用いた磁化プラズマの解析。また，ライフワークとして，プラズマ解析の基礎となる電子エネルギー分布や電子輸送係数を電磁界や観測モードなど各種条件に合わせて算出するため，ボルツマン方程式の数値解法であるプロパゲータ法のプログラムを継続的に開発中。
    語学が趣味で海外出張には現地語を学習してから臨む。

  • 姓

    菅原, スガワラ

  • 名

    広剛, ヒロタケ

  • ID各種

    200901037602526746

対象リソース

業績リスト

研究キーワード

  • 電子スオーム   ボルツマン方程式   磁化プラズマ   電子輸送   モンテカルロシミュレーション   プラズマエレクトロニクス   Electric Discharges   Plasma Electronics   

研究分野

  • エネルギー / プラズマ科学
  • ものづくり技術（機械・電気電子・化学工学） / 電気電子材料工学

経歴

  • 2004年 - 2007年 北海道大学大学院情報科学研究科情報エレクトロニクス専攻集積システム講座集積プロセス学研究室 助教授
  • 2004年 - 2007年 Associate Professor
  • 2007年 - 北海道大学大学院情報科学研究科情報エレクトロニクス専攻集積システム講座集積プロセス学研究室 准教授
  • 2007年 - Associate Professor
  • 1998年 - 2004年 北海道大学大学院工学研究科電子情報工学専攻集積材料デバイス工学講座集積電子材料工学分野 助教授
  • 1998年 - 2004年 Associate Professor,Associate Professor, Hokkaido University
  • 1997年 - 1998年 北海道大学大学院工学研究科電子情報工学専攻集積材料デバイス工学講座集積電子材料工学分野 助手
  • 1997年 - 1998年 Research Associate
  • 1995年 - 1997年 北海道大学工学部電子情報工学専攻集積材料デバイス工学講座集積電子材料工学分野 助手
  • 1995年 - 1997年 Research Associate
  • 1992年 - 1995年 北海道大学工学部電気工学科電気応用工学講座 助手
  • 1992年 - 1995年 Research Associate

学歴

  •         - 1992年   北海道大学   工学研究科   電気工学専攻
  •         - 1992年   北海道大学
  •         - 1990年   北海道大学   工学部   電気工学科
  •         - 1990年   北海道大学

委員歴

  • 2018年05月 - 現在   電気学会   編修専門第１部会委員
  • 2018年01月 - 現在   電気学会   プラズマ・パルスパワー技術委員会１号委員
  • 2012年04月 - 2017年12月   社団法人電気学会   プラズマ技術委員会１号委員
  • 2014年10月 - 2017年09月   電気学会   放電技術委員会２号委員
  • 2014年10月 - 2017年09月   電気学会   放電・プラズマ気相シミュレーション技法調査専門委員会委員長
  • 2014年08月 - 2016年07月   日本学術振興会   特別研究員等審査会専門委員
  • 2011年04月 - 2014年03月   社団法人電気学会   原子・分子衝突断面積および放電基礎データ調査専門委員会幹事
  • 2009年 - 2011年   社団法人電気学会   北海道支部協議員   社団法人電気学会
  • 2008年 - 2010年   社団法人電気学会   ２号代議員   社団法人電気学会
  • 2010年   日本工学教育協会   第59回工学・工業教育研究講演会実行委員   日本工学教育協会
  • 2007年 - 2009年   社団法人電気学会   北海道支部総務企画幹事   社団法人電気学会
  • 2006年 - 2009年   社団法人電気学会   放電技術委員会雷放電シミュレーション技術の現状とモデル間の相互比較調査専門委員会委員   社団法人電気学会
  • 2006年 - 2008年   社団法人電気学会   ２号代議員   社団法人電気学会
  • 2008年   社団法人電気学会   平成21年電気学会全国大会実行委員会幹事   社団法人電気学会
  • 2005年 - 2007年   社団法人電気学会   放電技術委員会荷電粒子、励起種、解離種ならびに光子と原子分子ダイナミックス調査専門委員会委員   社団法人電気学会
  • 2003年 - 2005年   社団法人電気学会   放電技術委員会・窒素中における放電とそのプラズマプロセスへの応用技術調査専門委員会委員   社団法人電気学会
  • 2002年 - 2004年   社団法人電気学会   放電技術委員会・窒素中における放電とそのプラズマプロセスへの応用技術調査専門委員会委員   社団法人電気学会
  • 2001年 - 2003年   社団法人電気学会   北海道支部会計幹事   社団法人電気学会
  • 1999年 - 2000年   社団法人電気学会   平成12年電気学会基礎・材料・共通部門大会実行委員会幹事   社団法人電気学会
  • 1997年   社団法人電気学会   論文委員会委員   社団法人電気学会
  • 1994年   社団法人電気学会   平成７年電気学会全国大会実行委員会事務局   社団法人電気学会

受賞

  • 2016年07月 日本学術振興会 平成27年度特別研究員等審査会専門委員（書面担当）表彰
     
    受賞者: 菅原 広剛書面審査における有意義な審査意見
  • 2000年 電気学会 電気学会優秀論文発表賞
     
    平成11年電気学会基礎・材料・共通部門大会 No. 59 (1999)菅原広剛、八幡貴志、酒井洋輔:「一定衝突周波数モデルにおける電磁界下電子群ドリフト速度のベクトル軌跡について」
  • 1997年 電気学会 電気学会優秀論文発表賞
     
    平成８年度電気学会基礎・材料・共通部門総合研究会 ED-96-218 (1996)菅原広剛、田頭博昭、酒井洋輔:「気体中の電子群重心移動速度－速度空間モーメント方程式による導出－」
  • 1995年03月 電気学会 電気学会論文発表賞(賞Ｂ)
     
    平成６年度電気関係学会北海道支部連合大会講演 No. 50 (1994)菅原広剛、酒井洋輔、田頭博昭:「気体の電離係数(CCXXXV)－完全吸収陽極近傍の電子後方拡散の効果－」
  • 1992年04月 電気学会論文発表賞(賞Ａ)
     
    電気学会放電研究会講演 ED-91-69 (1991)菅原広剛、酒井洋輔、田頭博昭:「ボルツマン方程式の直接数値解法（Convective Scheme）によるArガス中電子スオームのシミュレーション」
  • 1990年08月 電気学会 電気学会論文発表賞（賞Ｂ）
     
    受賞者: 菅原 広剛平成２年電気学会全国大会講演 No. 155 (1990)菅原広剛、酒井洋輔、田頭博昭:「気体の電離係数(CLXXXVIII)－ボルツマン方程式の直接数値解法－」

論文

  • Toshiaki Makabe, Hirotake Sugawara
    Plasma Sources Science and Technology 33 9 093001 - 093001 2024年09月01日 [査読有り][招待有り]
     
    Abstract Theoretical study of the electron kinetics (i.e. the velocity distribution and the transport parameter) in gases is generally conducted using the electron Boltzmann equation. The year 2022 marked 150 years since the formulation of the Boltzmann equation. Even in the last several decades, the historical progress has been made synchronously with the development of innovative technologies in gaseous electronics and in combination with the appearance of computers with sufficient speed and memory. Electron kinetic theory based on the Boltzmann equation has mostly been developed as the swarm physics in the hydrodynamic regime in the dc and radio frequency electric fields. In particular, the temporal characteristics are understood in terms of the collisional relaxation times between electron and gas molecule. There are two main theoretical approaches based on the Boltzmann equation for finding the velocity distribution. One is the traditional description of the electron kinetics, starting from the Boltzmann statistics in velocity space under a uniform density or a small density gradient of electrons. The other most recent approach is based on the phase-space tracking of the velocity distribution where the electron transport parameter is given by the moment of the electron density distribution in position space. In the present paper, we will explore the historical development of the electron Boltzmann equation with respect to three key items: collision term, solution method, and intrinsic electron transport in a hydrodynamic regime involved as the key elements in the low-temperature collisional plasma. The important topics listed in a table are briefly noted and discussed.
  • Luca Vialetto, Hirotake Sugawara, Savino Longo
    Plasma 7 1 121 - 145 2024年02月16日 [査読有り][通常論文]
     
    In this review, we detail the commonality of mathematical intuitions that underlie three numerical methods used for the quantitative description of electron swarms propagating in a gas under the effect of externally applied electric and/or magnetic fields. These methods can be linked to the integral transport equation, following a common thread much better known in the theory of neutron transport than in the theory of electron transport. First, we discuss the exact solution of the electron transport problem using Monte Carlo (MC) simulations. In reality we will go even further, showing the interpretative role that the diagrams used in quantum theory and quantum field theory can play in the development of MC. Then, we present two methods, the Monte Carlo Flux and the Propagator method, which have been developed at this moment. The first one is based on a modified MC method, while the second shows the advantage of explicitly applying the mathematical idea of propagator to the transport problem.
  • 小林 司,菅原 広剛,池田 圭
    Japanese Journal of Applied Physics 62 SL SL1020-1 - SL1020-10 2023年05月11日 [査読有り][通常論文]
     
    Abstract In order to accurately evaluate the electron acceleration process in the time evolution calculation of electron velocity distribution function (EVDF) based on the Boltzmann equation, an improved scheme blending upwind and central differences is introduced into the propagator method (PM). While the previous PM based on the upwind scheme needs fine cells to obtain an accurate EVDF at low electric fields, the improved PM is robust against coarse cells, which allows the reduction of the cell resolution. Calculations of the EVDF in Ar under RF electric fields demonstrated that the blending scheme can provide satisfactorily accurate results even with cells about tenfold larger than the upwind case at low reduced electric fields below 1 Td, which leads to much shorter computational time because the reduction of number of cells satisfactorily compensates the complexity of the blending scheme. This technique has been built into a new user-friendly PM software named BOSPROM.
  • Ryota Okazaki, Hirotake Sugawara
    Japanese Journal of Applied Physics 62 SL SL1003-1 - SL1003-8 2023年04月12日 [査読有り][通常論文]
     
    Abstract We investigate the partial resonance as one of the mechanisms of the electron heating to sustain inductively coupled plasmas driven under confronting divergent magnetic fields having an electron confinement function. Dependence of the partial resonance on operating conditions are observed. The electron energy gain (EEG), defined as a measure for evaluation of the partial resonance, was calculated by a Monte Carlo method under various magnetic field strengths and the RF antenna positions. It is confirmed that the partial resonance is observed in the resonant region even in low electric field and strong magnetic field. As an attempt to utilize the partial resonance, we demonstrate enhancement of the EEG and resulting ionization in the resonant region by placing the RF antenna near a region where the magnetic field lines running through the resonant region reach. 
  • 奥山由, 安沢友花, 菅原広剛
    Electrical Engineering in Japan 216 1 2023年02月07日 [査読有り][通常論文]
  • 奥山由, 安沢友花, 菅原広剛
    IEEJ Transactions on Fundamentals and Materials 142 10 397 - 405 2022年10月01日 [査読有り][通常論文]
  • 野原直人, 菅原広剛
    電気学会論文誌Ａ 142 6 288 - 289 2022年06月01日 [査読有り][通常論文]
  • Hirotake Sugawara, Yuya Nakata
    The European Physical Journal D 76 2 32-1 - 32-11 2022年02月 [査読有り][通常論文]
  • Hirotake Sugawara, Hikaru Iwamoto
    Japanese Journal of Applied Physics 60 4 046001 - 046001 2021年04月01日 [査読有り][通常論文]
  • 高橋宏徳,菅原広剛
    Japanese Journal of Applied Physics 59 3 036001-1 - 036001-6 2020年02月17日 [査読有り][通常論文]
  • Yuya Nakata, Takeshi Sato, Hirotake Sugawara
    IEEE Transactions on Plasma Science 49 1 83 - 90 2020年 [査読有り][通常論文]
  • 中島克博, 高橋宏徳, 菅原広剛
    Japanese Journal of Applied Physics 58 11 116001-1 - 116001-7 2019年11月 [査読有り][通常論文]
  • 菅原広剛
    Japanese Journal of Applied Physics 58 10 108002-1 - 108002-4 2019年10月 [査読有り][通常論文]
  • 小澤良輔, 菅原広剛
    電気学会論文誌Ａ 139 5 283 - 284 2019年05月 [査読有り][通常論文]
  • 菅原広剛
    Plasma Science and Technology 21 9 094001-1 - 094001-18 2019年 [査読有り][通常論文]
  • 高橋宏徳, 中島克博, 山本達平, 菅原広剛
    Japanese Journal of Applied Physics 57 12 126101-1 - 126101-7 2018年12月 [査読有り][通常論文]
  • Hirotake Sugawara
    Japanese Journal of Applied Physics 57 3 038001-1 - 038001-3 2018年03月01日 [査読有り][通常論文]
     
    In Monte Carlo simulations of electron swarms, sample electrons were copied periodically so that a sufficient number of samples are obtained in equilibrium after relaxation even under a severe attachment-dominated condition where most electrons vanish during the relaxation. The final sampling results were equivalent to those sampled by a conventional method, and the computational time conventionally wasted for the tracking of vanishing electrons was reduced drastically. The time saved can be utilized for tracking more samples to reduce statistical fluctuation. The efficiency of this technique and its limitation are discussed quantitatively together with details on its implementation.
  • 菅原広剛
    IEEE Transactions on Plasma Science 47 2 1071 - 1082 2018年02月 [査読有り][通常論文]
  • Hirotake Sugawara
    PLASMA SOURCES SCIENCE & TECHNOLOGY 26 4 044002-1 - 044002-16 2017年04月 [査読有り][通常論文]
     
    A propagator method (PM), a numerical technique to solve the Boltzmann equation (BE) for the electron velocity or energy distribution function (EVDF/EEDF) of electron swarms in gases, was customized to obtain the equilibrium solution quickly. The PM calculates the number of electrons in cells defined in velocity space using an operator called the propagator or Green's function. The propagator represents the intercellular transfer of electrons corresponding to the electron velocity change due to the acceleration by the electric field and the collisional events with gas molecules. The relaxation of the EVDF to its drift equilibrium solution proceeds with iterative propagator operations for the EVDF. Merits of the PM are that the series expansion of the EVDF as done in the BE analyses is not required and that time evolution of the electron swarm can be observed if necessary. On the other hand, in case only the equilibrium solution of the EVDF is wanted, the relaxation can be accelerated numerically. A demonstration achieved a shortening of the computational time by about three orders of magnitude. Furthermore, this scheme was applied to calculations of a set of electron transport parameters required in fluid-model simulations, i.e. the effective ionization frequency, the centroid drift velocity and the longitudinal diffusion coefficient, using the zeroth-, first- and second-order moment equations derived from the BE. A detailed description on the PM calculation was presented.
  • 山本達平, 菅原広剛
    電気学会論文誌Ａ 137 6 363 - 364 2017年 [査読有り][通常論文]
     
    Confronting divergent magnetic fields (CDMFs) applied to an inductively coupled plasma have a function to confine electrons within the source region. We discussed judgement criteria for electron escape from the source region to the diffusion region across the separatrix of the CDMFs. Passage of electron gyration center across the separatrix and collisional events in the diffusion region were newly introduced as effective criteria.
  • Hirotake Sugawara, So Ogino
    JAPANESE JOURNAL OF APPLIED PHYSICS 55 7 07LD05-1 - 07LD05-6 2016年07月 [査読有り][通常論文]
     
    The electron motion under confronting divergent magnetic fields (CDMFs) applied to inductively coupled plasmas was simulated using a Monte Carlo method. The CDMFs induced by two coaxial dc coils confined electrons in one side of the separatrix of the CDMFs. However, electrons diffused across the separatrix mainly in two ways. One was the displacement of their gyrocenters due to scattering near the outer part of the separatrix. This process tended to occur for high-energy electrons with correspondingly large gyroradii. The other was passage through the weak magnetic field around the center. This process was allowed for low-energy electrons. The position-dependent selectivity about the electron energy was pointed out as a seed property for applications of the separatrix as a magnetic filter or shutter. The mechanism by which the functions of the separatrix emerge was explained from a viewpoint of electron motion under the CDMFs. (C) 2016 The Japan Society of Applied Physics
  • Yudai Minami, Yusuke Asami, Hirotake Sugawara
    IEEE TRANSACTIONS ON PLASMA SCIENCE 42 10 2550 - 2551 2014年10月 [査読有り][通常論文]
     
    We simulated electron motion in confronting divergent magnetic fields (CDMFs) using a Monte Carlo method. Fundamental plasma structure under a filtering effect of the separatrix of CDMFs was depicted. Hot spots of energy deposition into electrons were observed not only near the RF antenna but also near the equistrength surface of 2B(ECR), where B-ECR is the RF-resonant magnetic field strength.
  • Yusuke Asami, Hirotake Sugawara
    IEEE TRANSACTIONS ON PLASMA SCIENCE 42 10 2540 - 2541 2014年10月 [査読有り][通常論文]
     
    Production and transport of etchant ions in an RF (13.56 MHz) magnetic neutral loop (NL) discharge plasma in the presence and the absence of an ac (1.695 MHz) substrate bias were simulated using a Monte Carlo method. Under the bias, electrons accelerated near the NL were attracted to the substrate and produced more etchant ions near the substrate than in the absence of the bias. Effect of bias leading to enhancement of the ion inflow to the substrate was discussed.
  • Takuhei Yoshida, Yohei Sakurai, Hirotake Sugawara, Akihiro Murayama
    ELECTRICAL ENGINEERING IN JAPAN 185 1 9 - 16 2013年10月 [査読有り][通常論文]
     
    We simulated electron and ion motions in a neutral loop discharge plasma under the control of the foot of the separatrix sweeping over the substrate and the neutral loop moving within a short distance from the RF antenna by a Monte Carlo method. We analyzed the distributions of ion production and ion flux to the substrate. We found that ion production is sensitive to the magnetic field gradient rather than to the electric field strength. Moreover, by superposing the flux distributions weighted by the residence time of the foot of the separatrix on the substrate, we obtained a uniform time-averaged distribution of the ion flux to the substrate in a radius range of r = 4.0 to 14.0 cm with sigma/m = 0.25% (m: mean, sigma: standard deviation). (c) 2013 Wiley Periodicals, Inc. Electr Eng Jpn, 185(1): 9-16, 2013; Published online in Wiley Online Library (wileyonlinelibrary.com). DOI 10.1002/eej.22417
  • Hirotake Sugawara
    JOURNAL OF THE PHYSICAL SOCIETY OF JAPAN 82 6 064501-1 - 064501-6 2013年06月 [査読有り][通常論文]
     
    Electron gyration and scattering in a gas under antiparallel gradient magnetic fields (APGMFs) and rf electric field were simulated using a Monte Carlo method to identify the factors of the electron confinement effect of the APGMFs. The inward shift of the gyrocenter at electron scattering worked as an additional factor to electron meandering, which has been regarded as the primary factor of the electron confinement. This scattering-based inward electron attraction occurs by a simple principle that the gyroradius and residence time of electrons become greater under weaker magnetic fields in the inner region of the APGMFs.
  • Hirotake Sugawara
    JAPANESE JOURNAL OF APPLIED PHYSICS 52 5 056101-1 - 056101-5 2013年05月 [査読有り][通常論文]
     
    Electron conduction in CF4 under a simplified model quadrupole magnetic field (QMF) and rf electric fields at f = 13.56, 27.12, and 40.68 MHz was analyzed using a Monte Carlo method, following the manner of electron swarm analyses. This field configuration simulated current paths of magnetic neutral loop discharge plasmas used for dry etching. The electron conduction therein was focused on as an elemental process of plasma power deposition. The electrons formed a star-shaped distribution with four rays along the separatrices of the QMF. The effective electron conductivity quantified as the amplitude of mean electron velocity was high around the central axis. The time-averaged mean electron velocities in the four regions of the QMF were biased under a rectification effect of the magnetic field. With increasing f, the high-conductivity region expanded and ionization was promoted. An enhancement of plasma power deposition under high driving frequencies was indicated. (C) 2013 The Japan Society of Applied Physics
  • 吉田拓平, 櫻井洋平, 菅原広剛, 村山明宏
    電気学会論文誌 132 4 278 - 283 2012年 [査読有り][通常論文]
     
    We simulated the electron and ion motions in a neutral loop discharge plasma under the control of the foot of separatrix sweeping over a substrate and the neutral loop moving within a short distance from the RF antenna by a Monte Carlo method. We analyzed the distributions of ion production and ion flux to the substrate. We revealed that ion production is sensitive to the gradient of magnetic field rather than the electric field strength. Moreover, by superposing the flux distributions weighted by the passage time of the foot of separatrix on the substrate, we obtained a uniform time-averaged distribution of ion flux to the substrate in a radius range of r = 4.0-14.0 cm with σ/m = 0.25% (m: the average, s: the standard deviation). © 2012 The Institute of Electrical Engineers of Japan.
  • Tsuyoshi Osaga, Hirotake Sugawara, Yohei Sakurai
    PLASMA SOURCES SCIENCE & TECHNOLOGY 20 6 065003-1 - 065003-7 2011年12月 [査読有り][通常論文]
     
    Electron motion in a CF(4) magnetic neutral loop discharge (NLD) plasma was simulated using a Monte Carlo method. The spatial distribution of electrons illustrated the fundamental structure of the NLD plasma, and its dynamics were depicted from the distributions of the mean electron energy, the electron energy gain and the azimuthal electron flux. The peak of mean electron energy appeared at the neutral loop (NL), which confirmed that electrons gain energy near this loop. High mean electron energy was observed not only near the NL but also along the separatrices of the quadrupole magnetic field. Energetic electrons were transported along the separatrices and induced ionization at those locations. However, the electron distribution had valleys along the separatrices, because electrons accelerated near the NL were likely to leave this region where the binding of the magnetic field is weak. The azimuthal electron flux representing the loop plasma current showed that the electron conduction path around the NL, which has conventionally been modelled as a ring conductor, has a particular directionality due to the rectification effect of antiparallel magnetic fields composing the quadrupole magnetic field. The directionality in the upper and lower regions of the quadrupole magnetic field was opposite to that in the inner and outer regions.
  • Yushi Fukuda, Junichi Takayama, Norichika Asaoka, Yoshiyuki Suda, Hirotake Sugawara, Akihiro Murayama
    IEEE TRANSACTIONS ON PLASMA SCIENCE 39 11 3133 - 3139 2011年11月 [査読有り][通常論文]
     
    Remote plasma is an efficient method to synthesize carbon nanomaterials because it can avoid the damage due to ion bombardment and can supply active precursor species for their growth. We introduced dc plasma to the CH(4) thermal chemical vapor deposition (CVD) of carbon nanotubes (CNTs) in order to decompose CH(4) efficiently. The effect of the plasma was investigated through the gas contents observed by a quadrupole mass spectrometer and the quality of CNTs evaluated by Raman spectroscopy. High-quality single-walled CNTs grew on Al(2)O(3)/Fe/Al(2)O(3) substrates, on which CNTs did not grow in conventional thermal CVD. The intensity ratio of the G peak to D peak in Raman spectra (G/D ratio) was 41.8. Since the CNT yield increased with the discharge voltage, the activation of precursors by the plasma seemed to promote the CNT growth significantly. At this time, C(2)H(x) and C(3)H(x), produced in the plasma, also increased. Thus, they seemed to contribute to the CNT growth. The diameter of CNTs, as well as the CNT yield, increased with increasing discharge voltage. This indicated the presence of the optimum precursor supply rate depending on the size of catalyst particles to keep them active.
  • Tsuyoshi Osaga, Yohei Sakurai, Hirotake Sugawara
    IEEE TRANSACTIONS ON PLASMA SCIENCE 39 11 2546 - 2547 2011年11月 [査読有り][通常論文]
     
    We simulated the electron motion in a quadrupole magnetic field that drives a neutral loop discharge (NLD) plasma using a Monte Carlo method. The structure of the electron distribution in the NLD plasma was observed. Electrons underwent reciprocating motion between reflecting boundaries as they gyrated along magnetic field lines. During this motion, they deviated from the neutral loop and concentrated beside the separatrices of the quadrupole magnetic field.
  • Yohei Sakurai, Tsuyoshi Osaga, Hirotake Sugawara
    IEEE TRANSACTIONS ON PLASMA SCIENCE 39 11 2550 - 2551 2011年11月 [査読有り][通常論文]
     
    The production and transport of ions in a neutral loop (NL) discharge plasma under a quadrupole magnetic field are investigated using a Monte Carlo method. The ion production occurred around the separatrices of the quadrupole magnetic field and decreased as the NL shrunk. The ion flux on the substrate had its peak near the foot of a separatrix. It was demonstrated that the time-averaged ion flux on the substrate can be distributed uniformly by controlling the magnetic field so that the foot of the separatrix sweeps over the substrate.
  • Hirotake Sugawara, Tsuyoshi Osaga, Hiromori Yamamoto
    PLASMA SOURCES SCIENCE & TECHNOLOGY 20 5 055002-1 - 055002-7 2011年10月 [査読有り][通常論文]
     
    Electron transport in CF4 under linearly gradient antiparallel magnetic fields was analysed in order to investigate the fundamental properties of magnetic neutral loop discharge plasmas used for material processing. The electron motion was simulated by a Monte Carlo method under a radio-frequency (rf) electric field applied perpendicularly to the directions of both the magnetic field and its gradient. Two typical electron motions, meandering in a weak magnetic field and gyration in a strong magnetic field, were observed with particular directionalities. The meandering electrons drifted forward on average similarly to those under a dc electric field. The gyration induced an electron drift towards the inverse direction. The direction of electron flux was dependent not only on the rf phase but also on the distance from the magnetically neutral midplane between the antiparallel magnetic fields. The electron conduction path formed along the midplane had a structure consisting of forward and inverse lanes. A peculiar result was that the direction of local electron flux was always forward in the strong magnetic field whereas the drift of gyrating electrons was towards the inverse direction. This seemingly paradoxical result can appear in the presence of the density gradient of the electron distribution.
  • Hirotake Sugawara, Tsuyoshi Osaga, Hideo Tsuboi, Kiyoshi Kuwahara, Seiji Ogata
    JAPANESE JOURNAL OF APPLIED PHYSICS 49 8 086001-1 - 086001-5 2010年08月 [査読有り][通常論文]
     
    A Monte Carlo simulation of electron transport in electric and magnetic fields was performed to analyze experimental data of the electron temperature T(e) and electron number density n(e) measured in a magnetic neutral loop discharge (NLD) plasma driven in Ar at 0 13 Pa T(e) and n(e) in the vicinity of the substrate holder were measured with a triple probe, and their radial profiles had peaks at different radial positions The simulation reproduced these peak positions well under the chosen boundary condition that the electron reflectivity of the side wall was lower than that of the reactor ceiling and the substrate holder It was explained that the T(e) peak was formed by high-energy electrons transported from the, neutral loop along a separatrix of the quadrupole magnetic field and that the n(e) peak consisted of electrons undergoing reciprocating motion between the reactor ceiling and the substrate (C) 2010 The Japan Society of Applied Physics
  • 水野巧一朗, 菅原広剛, 村山明宏
    電気学会論文誌Ａ 130 4 319 - 324 2010年 [査読有り][通常論文]
     
    We have composed amorphous fluorocarbon (a-C:F) polymer films, which are low-k insulators, by plasma-enhanced chemical vapor deposition using C8F18 as the feedstock. In our previous attempts to form multi-layered a-C:F films under various combinations of plasma power density P and deposition time T, it was found that a-C:F films deposited at a high P (typically 2.0 W/cm3) in a short T (∼1 min) became wavy when they are deposited on other a-C:F films deposited at a low P (0.2 W/cm3). A heating experiment for the films showed that their thermal tolerance was dependent on P films deposited at 0.2 W/cm3 started to melt at 250-270°C, while those deposited at 2.0 W/cm3 did at 350-370°C. It is considered that species of the precursors for the a-C:F film deposition are determined by P through the degree of C8F18 decomposition. The waving is induced by heating the lower layer during the deposition of the upper layer. The following conditions are necessary for the waving the lower layer is deposited at a low P, the upper layer is deposited at a temperature at which the lower layer expands but does not melt, and the upper layer is not thickened excessively to keep its flexibility. © 2010, The Institute of Electrical Engineers of Japan. All rights reserved.
  • Yoshiyuki Suda, Atsushi Okita, Junichi Takayama, Akinori Oda, Hirotake Sugawara, Yosuke Sakai, Shinichiro Oke, Hirofumi Takikawa
    IEEE TRANSACTIONS ON PLASMA SCIENCE 37 7 1150 - 1155 2009年07月 [査読有り][通常論文]
     
    We have successfully grown carbon nanotubes (CNTs) by plasma-enhanced chemical vapor deposition (PECVD) using alcohol. When 0.01-wt% ferrocene was added to the alcohol, vertically aligned CNTs grew at 650 degrees C. By contrast, a few CNTs and mostly carbon nanoparticles; were obtained by pure alcohol PECVD even though the Fe catalyst was coated on Si substrates. Comparing this PECVD experiment with thermal alcohol CVD showed that only the PECVD method can be used to grow CNTs under the reported experimental conditions. To understand the plasma properties for CNT growth, particularly plasma species contained in a gas phase of alcohol plasma, the plasma was analyzed using optical-emission spectroscopy (OES) and quadrupole mass spectrometry (QMS). From the OES measurement, emission peaks from the excitation states Of C(2), CH, CHO, CH(2)O, CO, H, O(2), C(+), and CO(+) were identified, while the QMS measurement also showed the existence of H(2), O, and CO. These results indicate that, in alcohol plasma, oxidants and reductants exist together and potentially promote/suppress CNT growth depending on the process conditions. The contribution of C(x)H(y) (x >= 1, y >= 3) radicals, which were produced by decomposition reactions in alcohol plasma as a CNT precursor, is discussed.
  • Krishnendu Bhattacharyya, Yoshiyuki Suda, Atsushi Okita, Takeshi Saito, Atsushi Ozeki, Masayuki Maekawa, Junichi Takayama, Hirotake Sugawara, Yosuke Sakai
    JOURNAL OF NANOSCIENCE AND NANOTECHNOLOGY 9 3 1897 - 1903 2009年03月 [査読有り][通常論文]
     
    Based on recent theoretical photoluminescence intensity calculations, the population densities of sixteen semiconducting single-walled carbon nanotubes grown by alcohol catalytic chemical vapor deposition were estimated for two different temperatures. The profiles of population density merely as functions of tube-diameter or chiral angle are found to be widely scattered. However, systematic profiles are detectable when separately split into (2n+m), (n+2m) and (n-m) family arrays. Apart from these well-knit family behaviors, the population densities of the group of nanotubes forming another three possible series [viz. constant-n, constant-m and constant-(n+m)] also show evidence of good correlations. Hence, a two-dimensional chiral-zone selective growth principle is hypothesized.
  • Hirotake Sugawara, Yosuke Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 41 13 135208-1 - 135208-7 2008年07月 [査読有り][通常論文]
     
    Electron motion in CF(4) at 0.67 Pa under crossed electric (E) and magnetic (B) fields was simulated by a Monte Carlo method to investigate fundamental properties of electron transport in neutral loop discharge plasmas for dry etching. As a simplified model of the electron path in the plasma, a magnetically neutral plane was assumed between linearly gradient antiparallel B fields, and a uniform E field was applied along the neutral plane perpendicularly to the B fields. The electron behaviour showed two contrasting modes depending on the direction of the B fields relative to the E field. In the field configuration which confines the electrons near the neutral plane, values of the mean electron energy, the drift velocity and the effective ionization frequency were close to those under the dc E field without B field. On the other hand, in the opposite B field configuration, the electrons hardly drifted along the E field, but instead, they showed a constant lateral diffusion driven by the E x B drift. A reverse-blocking effect of the gradient antiparallel B fields is reported.
  • Youhei Hizume, Atsushi Okita, Akinori Oda, Hirotake Sugawara, Yoshiyuki Suda, Yosuke Sakai
    IEEJ Transactions on Fundamentals and Materials 128 10 5 - 628 2008年 [査読有り][通常論文]
     
    A CH4/H2 radio frequency capacitively coupled plasma used for synthesis of carbon nanotubes (CNTs) was simulated with a 1-dimensional fluid model. The pressure of CH4 was fixed at 1 Torr and that of H2 varied from 0 to 9 Torr. When the H2 mixture ratio was 0.1, the fluxes of H2 + and H3 + to the substrate were respectively about 5 and 10 times those for the case without H2 mixing. In our previous experiment, a slight mixture of H 2 promoted CNT growth. We inferred that this was because H 2 maintained the activity of the catalyst particle from which CNT grew. The present numerical result supports our hypothesis. In addition, C xHy(x, y≥2) accumulated with time. This accumulation is not negligible because it lowers the electron temperature. We estimated the H2 ionization rate in CH4/CxHy is about 2/3 of that in CH4. It is pointed out that we need to consider the accumulation of CxHy to improve the accuracy of the estimation of CNT growth. © 2008 The Institute of Electrical Engineers of Japan.
  • Dual distributions for the metallic and semiconducting single-walled carbon nanotubes observed by Raman spectroscopy
    K. Bhattacharyya, Y. Suda, A. Okita, T. Saito, A. Ozeki, M. Maekawa, J. Takayama, H. Sugawara, Y. Sakai
    Carbon - Sci. Tech. 1 4 - 8 2008年 [査読有り][通常論文]
  • Atsushi Okita, Yoshiyuki Suda, Akinori Oda, Junji Nakamura, Atsushi Ozeki, Krishnendu Bhattacharyya, Hirotake Sugawara, Yosuke Sakai
    CARBON 45 7 1518 - 1526 2007年06月 [査読有り][通常論文]
     
    Carbon nanotubes (CNTs) were synthesized using CH4/H-2 plasmas and plasmas simulated using a one-dimensional fluid model. The thinnest and longest CNTs with the highest number density were obtained using CH4/H-2 = 27/3 sccm at 10 Torr. These conditions allowed CNTs to grow for 90 min without any meaningful loss of catalyst activity. However, an excess H-2 Supply to the CH4/H-2 mixture plasma made the diameter distribution of the CNTs wider and the yield lower. Hydrogen concentration is considered to affect catalyst particle size and activity during the time interval before starting CNT growth (=incubation period). With CH4/H-2 = 27/3 sccm for a growth time of 10 min efficient CNT growth was achieved because the amount of carbon atoms in the CNTs and that calculated from simulation showed good agreement. The effect of hydrogen etching on CNTs was analyzed by scanning electron microscopy and X-ray photoelectron spectroscopy by observing CNTs treated by H-2 plasma after CNT growth. It was confirmed that (a) multi-walled CNTs were not etched by the H-2 plasma, (b) the C 1s XPS spectra of the CNTs showed no chemical shift after the treatment, and (c) C-H bonds were produced in CNTs during their growth. (c) 2007 Elsevier Ltd. All rights reserved.
  • Hirotake Sugawara, Naoki Mori, Yosuke Sakai, Yoshiyuki Suda
    JOURNAL OF COMPUTATIONAL PHYSICS 223 1 298 - 304 2007年04月 [査読有り][通常論文]
     
    Techniques to reduce the computational load for determination of electron-molecule collisions in Monte Carlo simulations of electrical discharges have been presented. By enhancing the detection efficiency of the no-collision case in the decision scheme of the collisional events, we can decrease the frequency of access to time-consuming subroutines to calculate the electron collision cross sections of the gas molecules for obtaining the collision probability. A benchmark test and an estimation to evaluate the present techniques have shown a practical timesaving efficiency. (c) 2006 Elsevier Inc. All rights reserved.
  • Yuki Hayakawa, Yoshiyuki Suda, Tamotsu Hashizume, Hirotake Sugawara, Yosuke Sakai
    JAPANESE JOURNAL OF APPLIED PHYSICS PART 2-LETTERS & EXPRESS LETTERS 46 12-16 L362 - L364 2007年04月 [査読有り][通常論文]
     
    Hydrogen gas sensors based on single-walled carbon nanotube (SWCNT) films were fabricated with Pd comb-like electrodes. The SWCNT films were deposited by Mo/Co catalyst-assisted chemical vapor deposition using alcohol. The film included a high-density network of SWCNTs. The SWCNT film sensor showed a large conductance change and a fast response (20s) upon exposure to 1% H-2 gas in air at room temperature, while the conductance change was negligible for devices using Au electrodes.
  • Atsushi Okita, Atsushi Ozeki, Yoshiyuki Suda, Junji Nakamura, Akinori Oda, Krishnendu Bhattacharyya, Hirotake Sugawara, Yosuke Sakai
    JAPANESE JOURNAL OF APPLIED PHYSICS PART 1-REGULAR PAPERS BRIEF COMMUNICATIONS & REVIEW PAPERS 45 10B 8323 - 8329 2006年10月 [査読有り][通常論文]
     
    We synthesized vertically aligned carbon nanotubes (CNTs) using multilayered catalyst thin films (Fe/Al2O3 and Al2O3/Fe/Al2O3) by RF (13.56 MHz) CH4/H-2/Ar plasma-enhanced chemical vapor deposition. Pretreatment of the catalyst is crucial for CNT growth. In this paper, we analyzed the effect of catalyst reduction on CNT growth. Catalyst thin films on substrates were reduced by H-2 plasma pretreatment at 550 degrees C to form nanometer-sized catalyst particles. The multilayered thin films were analyzed; the chemical composition and oxidation state by X-ray photoelectron spectroscopy (XPS) and the surface morphology by scanning electron microscopy (SEM). The Fe 2p peak of the XPS spectra showed that FexOy in the as-deposited catalyst was effectively reduced to Fe by a pretreatment of duration 4 min. Using this catalyst, we obtained CNTs with an average diameter of 10.7 nm and an average length of 5.3 mu m. However, pretreatment longer than 4 min resulted in shorter CNTs and the Fe peak was shifted from Fe to Fe3O4. These transitions (Fe2O3 -> Fe3O4 -> Fe -> Fe3O4) can be explained by the enthalpy of the oxides. This result indicates the presence of an optimum ratio between Fe and FexOy to maximize the CNT lengths.
  • H Sugawara, Y Sakai
    JAPANESE JOURNAL OF APPLIED PHYSICS PART 1-REGULAR PAPERS BRIEF COMMUNICATIONS & REVIEW PAPERS 45 6A 5189 - 5196 2006年06月 [査読有り][通常論文]
     
    When an electron swarm in gas is composed of component electron swarms individually in drift equilibrium, the higher-order diffusion coefficients (HDCs) of the composite electron swarms are equal to those. of the component electron swarms. We have derived this equality theoretically and have examined it by numerical simulation. The HDCs are the time derivatives of higher-order cumulants, which are quantities characterizing the shape of a spatial electron distribution. This fact seemingly indicates the dependence of the HDCs of the composite electron swarms on the arrangement of the component electron swarms. However, the equality holds irrespective of the relative positions and electron populations of the component electron swarms. We have given a consistent explanation to these facts. We have also discussed electron swarm development from dispersed or multiple electron sources that may appear in practical experiments.
  • MA Bratescu, Y Yoshizaki, Y Suda, Y Sakai, H Sugawara, O Takai
    THIN SOLID FILMS 506 145 - 149 2006年05月 [査読有り][通常論文]
     
    A dielectric barrier discharge (DBD) in Ar/CH4 gas mixtures used for amorphous hydrogenated carbon (a-C:H) film deposition was investigated. It was found that the a-C:H film properties (the ratios of sp(3)/Isp(2) and CH3/CH2, the number density of hydrogen atoms in the film and the energy band gap) correlate with the number densities of Ar(I-S5) atom and CH(chi(2) Pi) radical. A high value of the energy band gap, which corresponds to a high content of sp(3) fraction in the film was obtained, when the CH radical number density in the plasma was relatively small (similar to 10(8) cm(-3)). The DBD breakdown voltage decreased when the a-C:H film was deposited on the dielectric surface of the electrode. (c) 2005 Published by Elsevier B.V.
  • A Okita, Y Suda, A Ozeki, H Sugawara, Y Sakai, A Oda, J Nakamura
    JOURNAL OF APPLIED PHYSICS 99 1 014302-1 - 014302-7 2006年01月 [査読有り][通常論文]
     
    Carbon nanotubes (CNTs) were grown on Si substrates by rf CH4 plasma-enhanced chemical vapor deposition in a pressure range of 1-10 Torr, and then characterized by scanning electron microscopy. At 1 Torr, the CNTs continued growing up to 60 min, while their height at 4 Torr had leveled off at 20 min. CNTs hardly grew at 10 Torr and amorphous carbon was deposited instead. CH4 plasma was simulated using a one-dimensional fluid model to evaluate the production and transport of radicals, ions, and nonradical neutrals. The amount of simulated carbon supplied to the electrode surface via the flux of radicals and ions such as CH3, C2H5, and C2H5+ was consistent with estimations from experimental results.
  • Analysis of CH4 rf plasmas for carbon naotube growth: Comparison between computer modeling and experimental results
    A.Okita, Y.Suda, A.Ozeki, H.Sugawara, Y.Sakai, A.Oda, J.Nakamura, 中村, 潤児
    J.Appl.Phys 99 014302  2006年01月 [査読有り][通常論文]
  • H Sugawara, A Oda, Y Sakai
    IEEE TRANSACTIONS ON PLASMA SCIENCE 33 2 548 - 549 2005年04月 [査読有り][通常論文]
     
    Monte Carlo simulation shows that the electron velocity distribution function (EVDF) in CF4 in crossed electric and magnetic fields (E x B fields) has two peaks at low E/B values. The EVDF peaks are formed around a pair of points in velocity space at which electrons are stationary against the acceleration under the E x B fields and scattering by collisions with CF4 molecules. The stationary points appear when the E/B value is lower than the electron speed at which the total electron collision frequency in CF4 takes its minimum.
  • Investigation of excited species in a carbon ablation plume in nitrogen gas environment
    M. A. Bratescu, Y. Sakai, N. Sakura, D. Yamaoka, Y. Suda, H. Sugawara
    Plasma Processes and Polymers 489 - 498 2005年 [査読有り][通常論文]
  • A. Ozeki, Y. Suda, A. Okita, J. Nakamura, A. Oda, Y. Sakai, H. Sugawara
    Materials Research Society Symposium Proceedings 901 523 - 528 2005年 
    Optical emission spectra of a CH4/H2/Ar gas mixture plasma were observed during carbon nanotube (CNT) growth in RF plasma-enhanced chemical vapor deposition. CNTs with diameters of ∼10-30 nm and length of ∼6 (im were grown on double- and triple-layered films of catalyst/support materials (FexOy/TiO2 and Al2O 3/FexOy/Al2O3) at the total gas pressures of 1-10 Torr with gas flow rates of CH4 - 27 seem, H2 = 3 seem, and Ar = 1 seem. The number density of CNTs increased with the gas pressure, and Al2Oy/Fe xOy/Al2O3 (each thickness of 1 nm) film yielded the thinnest CNTs with a high number density among the present catalysts. The spatial distributions of H atom relative density in the plasma were obtained by actinometry. The H relative density decreased with the pressure, and this suggests the suppression of CH3 radical generation in the plasma. © 2006 Materials Research Society.
  • H Sugawara, T Ishigaki, Y Hirochi, Y Sakai
    JAPANESE JOURNAL OF APPLIED PHYSICS PART 1-REGULAR PAPERS SHORT NOTES & REVIEW PAPERS 43 11A 7705 - 7706 2004年11月 [査読有り][通常論文]
     
    It has been predicted that near-threshold ionization (NTI) in a gaseous dielectric inhibits the development of electron avalanche when the gaseous dielectric has a sufficient capability for low-energy electron attachment. The NTI leaves little energy for the primary and secondary electrons involved in the ionization; thus, both electrons can be captured by dielectric gas molecules without further ionization. A computational estimation indicates that this process can occur in SF6.
  • MA Bratescu, N Sakura, D Yamaoka, Y Sakai, H Sugawara, Y Suda
    APPLIED PHYSICS A-MATERIALS SCIENCE & PROCESSING 79 4-6 1083 - 1088 2004年09月 [査読有り][通常論文]
     
    In this work, we investigated a carbon plasma plume produced by laser ablation of a graphite target in a nitrogen gas environment. The spatial distributions of C and N atoms were measured by time-resolved absorption spectroscopy. The spatial distributions of the relative densities of CN radicals, C-2, and C-3 molecules were measured using time-resolved emission spectroscopy. We determined that nitrogen gas produced an increase in carbon atom and molecule densities in the ablation plume. It was observed that the addition of RF plasma to the plume increased the CN radicals and C atom densities, and decreased the C-2 and C-3 molecule densities. The RF plasma changed the evolution of various fractional species of C, N, CN, C-2, and C-3 in the ablation plume. The chemical reactions with and without RF plasma were explained using the evolution and density of the fractional species of C, N, CN, C-2, and C(3)in the plume.
  • C Biloiu, IA Biloiu, Y Sakai, H Sugawara, A Ohta
    JOURNAL OF VACUUM SCIENCE & TECHNOLOGY A 22 4 1158 - 1165 2004年07月 [査読有り][通常論文]
     
    Amorphous fluorocarbon polymer films (a-C:F) have been grown by plasma enhanced chemical vapor deposition from a new precursor, namely perfluoro-octane (C8F18) vapor. The dielectric and insulating properties of the films have been assessed by means of capacitance-voltage and current-voltage characteristics, breakdown voltage measurements, scanning electron microscopy, and ellipsometric analyses. In the investigated frequency range, 120 Hz-1 MHz, the films have a low dielectric constant (approximate to2.4) and a low dielectric loss (less than or equal to3 x 10(-2)). The electrical conduction mechanism is ohmic in the weak electric field regime, with a bulk resistivity of 4 x 10(15) Omega cm, and a trap modulated space charge limited conduction in the strong electric field regime. For a film thickness between 1 and 12 mum, the dielectric strength varies with film thickness according to E(B)similar toh(-0.53). For a film thickness of 1 mum the film dielectric strength was 2.7 MV/cm while for 12 mum it decreased to 0.9 MV/cm. Measurements of the breakdown voltage in nitrogen gas in the low pressure regime between a-C:F coated cathode and copper anode showed an enhancement of the breakdown voltage by a factor of 2.5 with respect to the case of copper-copper electrodes. (C) 2004 American Vacuum Society.
  • 石垣卓也, 広地祐樹, 菅原広剛, 酒井洋輔
    電気学会論文誌Ａ 124 2 126 - 131 2004年 [査読有り][通常論文]
     
    In a recent trend that reduction of SF6 used for electric power equipment is required to avoid global warming, dilution of SF6 by other gases of less environmental load is considered as a practical solution. Gas mixtures such as SF6/N2 and SF6/c-C4F8 have a desirable feature called synergy effect that the critical values of reduced electric field (E/N)lim are higher than the partial-pressure-weighted values obtained from (E/N)lim of their components. In order to obtain a guide to seek for new insulating gas mixtures with the synergy effect, we have investigated the synergy mechanism from a viewpoint of collaboration between electron energy loss by collisions and electron attachment. The following mechanisms have been indicated as primary factors of the synergy effect (i) a large attachment cross section of SF6 lying in a low electron energy region plays the main role of the electron capture (ii) in order to feed electrons into the low energy region of (i) efficiently, a steep rise of an inelastic collision cross section at its onset energy is required. We have evaluated these factors by a Boltzmann equation analysis, and have synthesized the synergy effect using model gases for demonstration. © 2004, The Institute of Electrical Engineers of Japan. All rights reserved.
  • 吉崎泰直, Maria-Antoaneta Bratescu, 酒井洋輔, 須田善行, 菅原広剛
    電気学会論文誌Ａ 124 2 158 - 163 2004年 [査読有り][通常論文]
     
    A dielectric barrier discharge in Ar and CH4 gas mixtures used for a-C:H film deposition was investigated. The number densities of Ar(1s2: resonant level) and Ar(1s5: metastable levels) atoms were measured using a high sensitive detection method based on plasma modulation and laser absorption spectroscopy. The CH(X2П) radical number density was measured by an actinometry method. The excitation transfer process from the Ar(1s5) atom to the CH(A2Δ) radical was detected by a laser collisionally induced fluorescence method. The a-C:H film properties (the content of sp3 and sp2 hybridization and the content of the CH3 and CH2 hydrocarbon groups) were correlated with the number densities of the Ar(1s5) atom and the CH(X2П) radical. In the a-C:H film, maximum value of the sp3/sp2 was obtained, when the CH(X2П) number density in DBD plasma was minimum. © 2004, The Institute of Electrical Engineers of Japan. All rights reserved.
  • H Sugawara, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 36 16 1994 - 2000 2003年08月 [査読有り][通常論文]
     
    We have simulated electron acceleration under short intense impulse electric fields of 0.1-1 kV cm(-1) strength and similar to1 ns duration in N-2 at 65.5 Pa. The electron acceleration is completed in a practically collisionless condition when the impulse duration is comparable with or shorter than the electron mean free time. On the basis of this feature, we have proposed a way to produce a group of monoenergetic electrons. The electron energy after acceleration is controllable by the impulse strength and duration, and we have demonstrated selective promotion of the excitation producing N-2(A(3)E(u)(+)) metastable excited molecules. The selectivity has been quantified as the efficiency of the electron energy utilized for the target reaction, and efficiencies of 20-30% have been achieved in this simulation. These values are higher than those of electron swarms in DC equilibrium, in which the efficiency is at most 10%. The electron behaviour under impulse electric fields and the dependence of the energy efficiency on the impulse profile have been investigated.
  • SY So, A Oda, H Sugawara, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 35 22 2978 - 2987 2002年11月 [査読有り][通常論文]
     
    CF4 radio frequency (rf) (13.56 MHz) capacitively coupled plasmas were simulated using a two-dimensional fluid model of cylindrical symmetry. We investigated their periodic steady state and transient behaviour after step changes of the amplitude of a power source voltage. The results were compared with those obtained by a one-dimensional model. The two-dimensional model showed that the plasmas in a harmonic steady state have asymmetric profiles due to generation of a dc self-bias through an asymmetric electrode system with grounded wall. The transient response of the plasmas in the two-dimensional model was qualitatively similar to those of the one-dimensional model except in the powered electrode edge region where the radial component of the electric field was significant.
  • MA Bratescu, Y Sakai, D Yamaoka, Y Suda, H Sugawara
    APPLIED SURFACE SCIENCE 197 257 - 262 2002年09月 [査読有り][通常論文]
     
    In a plume produced by excimer ArF laser, the number densities of excited atoms and electrons, and the expansion velocity of ablated carbon atoms in an Ar gas were measured and compared with those obtained in vacuum. A saturation in C(3(1)P(1)(0)) density was observed as the laser fluence increased. This saturation was shown to be closely related to the electron density saturation. The excitation of Ar(4s'[1/2](1)(0)) and C(3(1)P(1)(0)) was interpreted by electron collisions and was dependent on gas pressure and excimer laser fluence. (C) 2002 Elsevier Science B.V All rights reserved.
  • プラズマディスプレイ用混合希ガス中の励起移行効果
    内田諭, 菅原広剛, 酒井洋輔, 渡辺恒雄
    電気学会論文誌Ａ 122-A 8 729 - 735 2002年 [査読有り][通常論文]
  • H Sugawara, K Satoh, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 34 21 3191 - 3196 2001年11月 [査読有り][通常論文]
     
    The electron velocity distribution function (EVDF) of electron swarms in N-2 has a constricted part at an electron energy of around 2-3 eV. This feature has been explained by considering the barrier effect of vibrational excitation collisions that prevents the gain of energy by slow electrons. In order to investigate the barrier effect in detail further, we simulated the electron behaviour around the energy range of the barrier. We found that momentum transfer collisions also contribute to the barrier effect in collaboration with vibrational excitation. In addition, an evaluation of the electron flow in velocity space is thought to indicate a factor by which the presence of the barrier effect in N-2 is revealed in the EVDF; the depression of the electron interchange across the barrier is not only due to the barrier effect itself but also due to the small electron transfer associated with the energy loss by electronic excitation.
  • SY So, A Oda, H Sugawara, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 34 12 1919 - 1927 2001年06月 [査読有り][通常論文]
     
    CF4 radio-frequency (rf) plasmas are often driven in time-modulation modes to obtain higher etching performance. For dynamic and stable control of CF4 plasmas, we investigated the transient behaviour of the plasmas after step up or down of the amplitude of the power source voltage V-s (t). Using a one-dimensional fluid model, we simulated a CF4 capacitively coupled plasma that transfers from a periodic steady state to another at 0.2 Torr in a 20 mm parallel-plate gap. In addition, we compared such plasma responses to the step down to 0 V at the phases at which V-s(t) takes zero and peak values. Fundamental properties of transient CF4 plasmas were discussed.
  • S Uchida, H Sugawara, Y Sakai, T Watanabe, BH Hong
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 34 6 947 - 953 2001年03月 [査読無し][通常論文]
     
    The influence of Ar addition to conventional Plasma Display Panel gases, Xe/Ne and Xe/He mixtures, on the radiation characteristics and the discharge onset voltage is studied using a Boltzmann equation method coupled with rate equations for electrons and excited species. When Ar is added to Xe/Ne mixtures, the efficiency of ultraviolet light radiation from Xe (1S(2,4)) and Xe-2* increases 7-17% for E / p(0) between 3 and 100 V cm(-1) Torr(-1). In Xe/He mixtures, the efficiency increases over 20% for E / p(0) < 10 V cm(-1) Torr(-1) In addition, the visible light radiation of the 600 nm band from excited Ne, which makes the contrast worse, is effectively reduced. On the other hand, the discharge onset voltage increases with Ar addition because of a decrease in the effective secondary electron emission coefficient.
  • CxFy polymer film deposition in rf and dc C7F16 vapor plasmas
    Y. Sakai, M. Akazawa, Y. Sakai, H. Sugawara, M. Tabata, C. P. Lungu, A. M. Lungu
    Transactions of Korean Institute of Electrical and Electronic Material Engineers 2 1 1 - 6 2001年03月 [査読有り][通常論文]
  • CP Lungu, AM Lungu, Y Sakai, H Sugawara, M Tabata, M Akazawa, M Miyamoto
    VACUUM 59 1 210 - 219 2000年10月 [査読有り][通常論文]
     
    Polymer films containing CxFy have been deposited in Fluorinert vapor (C7F16) plasmas. The spatial and temporal evolution of optical emission of CF2 and C-2 radicals was monitored and kinetics of the C7F16 decomposition process was discussed. The films were deposited on stainless steel, glass, molybdenum and silicon wafers at room temperature with the vapor pressures, 40 and 100 Pa. Electrical properties such as the breakdown voltage and dielectric constant as well as the structure of the films were analyzed. The composition and characteristics of the films were determined by Fourier transform infrared spectroscopy and X-ray photoelectron spectroscopy techniques, scanning electron microscopy and electron spin resonance spectroscopy. (C) 2000 Elsevier Science Ltd. All rights reserved.
  • A Oda, H Sugawara, Y Sakai, H Akashi
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 33 12 1507 - 1513 2000年06月 [査読有り][通常論文]
     
    Xe dielectric barrier discharges at different gap lengths under applied pulse voltages with trapezoidal and sinusoidal waveforms were simulated using a self-consistent one-dimensional fluid model. In both waveforms, the light output power depended not only on the amplitude of voltage waveforms but also on the discharge gap length. At the narrower discharge gap, the light output efficiency was improved by increasing the time gradient of the applied voltage when the trapezoidal pulse is applied, and by decreasing the duty ratio in the sinusoidal case. In the present simulation, we adopted a fast numerical method for calculation of electric field introducing an exact expression of the discharge current.
  • H Sugawara, T Yahata, A Oda, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 33 10 1191 - 1196 2000年05月 [査読有り][通常論文]
     
    We calculated the drift velocity vector W Of electron swarms in crossed electric and magnetic fields using a Monte Carlo method. Values of W in F-2, CF4 and constant-collision-frequency (CCF) model gases at various values of the reduced electric and magnetic fields, E/N and B/N (N is the gas molecule number density), were illustrated as vector loci in velocity space. Under a CCF condition, W draws right semicircular loci with varying B and N at a constant E when the number of electrons is conservative. These semicircular loci can be a reference to characterize the W behaviour in real gases. We presented an analytical derivation of W in the CCF condition, extending its expression to that in the presence of ionization and attachment. Features of thr vector loci in F-2 and CF4 were discussed.
  • E Suetomi, M Tanaka, S Kamiya, S Hayashi, S Ikeda, H Sugawara, Y Sakai
    COMPUTER PHYSICS COMMUNICATIONS 125 1-3 60 - 74 2000年03月 [査読有り][通常論文]
     
    Direct current glow discharge plasma sources are being developed as reactors for immobilization of radioactive krypton gas recovered from the off-gas stream of a nuclear fuel reprocessing plant. In order to design these reactors, two-dimensional simulations have been performed using a fluid model coupled with Poisson's equation. Nonlinear governing equations were expressed by a fully implicit scheme and the equations were solved using the Newton-Krylov method. Mr, investigated algorithm performance of two methods from the Krylov projection techniques, namely the Bi-CGSTAB method and the transpose-free quasi-minimal residual (TFQMR) method. The results showed that the Bi-CGSTAB method converged considerably faster than the TFQMR method. Direct current glow discharge plasmas in three types of krypton immobilization apparatus were simulated to study the differences in electrode geometry. Simulation results agreed qualitatively with measurements of electron density. (C) 2000 Elsevier Science B.V. All rights reserved.
  • S Uchida, H Sugawara, Y Sakai, T Watanabe, BH Hong
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 33 1 62 - 71 2000年01月 [査読有り][通常論文]
     
    The electron swarm parameters of Xe/He and Xe/Ne mixtures are analysed with a Boltzmann equation method in which secondary generations of electrons through excited atoms are properly considered. The swarm parameters are examined for a wide range of the electron concentrations from the Townsend discharge region to thr: working condition of plasma display panels (PDPs). Using ionization coefficients obtained in this work and referred values of secondary electron emission coefficients from MgO film, the discharge onset voltage V-b is discussed. The Paschen curves around (p(0)d)(min) at which V-b has a minimum are compared for Xe/He and Xe/Ne mixtures. Under the working conditions for PDP discharges, the radiation rates of ultraviolet light (lambda = 130, 147 and 173 nm) from the Xe resonance levels and Xe excimer are calculated. A suitable discharge condition for designing PDP is discussed with respect to various values of the gas mixture ratio, gas pressure and reduced electric field.
  • CP Lungu, AM Lungu, M Akazawa, Y Sakai, H Sugawara, M Tabata
    JAPANESE JOURNAL OF APPLIED PHYSICS PART 2-LETTERS 38 12B L1544 - L1546 1999年12月 [査読有り][通常論文]
     
    An attempt to fabricate fluorinated carbon (fluoropolymer) films by rf (13.56 MHz) plasma enhanced chemical vapor deposition was made using novel fluorocarbon sourer materials of C7F16, (C3F7)(3)N/(C4F9)(3)N and C8F18/C8F16O. The deposited films were transparent and displayed excellent electrical properties as interlayer dielectrics for LSI using deep-submicron technology, namely, a dielectric constant as low as 2.0 and a dielectric strength higher than 2 MV/cm. The refractive index of these films was 1.38.
  • A Oda, Y Sakai, H Akashi, H Sugawara
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 32 21 2726 - 2736 1999年11月 [査読有り][通常論文]
     
    A self-consistent one-dimensional modelling of a Xe gas discharge between electrodes covered with dielectric barrier is presented for power frequencies from 50 kHz to 1 MHz and gas pressures from 10 to 400 Tan Spatiotemporal profiles of the concentration of electrons, ions, excited atoms and excimers are obtained. Excimers are mainly produced in the sheath regions. The efficiencies of spontaneous emission from excimers and resonance-state atoms increase with an increase in the input powers for gas pressures higher than 50 Torr. A characteristic period during which the barrier wall charge significantly influences the electric field in bulk region was found, and discharge properties in the period are discussed.
  • H Sugawara, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 32 14 1671 - 1680 1999年07月 [査読有り][通常論文]
     
    A simulation technique for the analysis of the transverse evolution of electron swarms in gases was developed based on moment equations derived from the Boltzmann equation. A numerical calculation of the moment equations for an electron swarm was performed using a propagator method and it was demonstrated that the propagator method can be used to calculate the higher-order transverse diffusion coefficients stably. Applying a Hermite expansion technique, the electron distribution in real space and other electron swarm parameters were derived as functions of the transverse position. The calculation result was verified by comparisons with those by a Monte Carlo simulation and other methods. Features of the transverse electron swarm evolution were presented.
  • 内田諭, 菅原広剛, Peter L. G. Ventzek, 酒井洋輔
    電気学会論文誌Ａ 118-A 6 622 - 628 1998年06月 [査読有り][通常論文]
  • H Sugawara, Y Sakai, H Tagashira, K Kitamori
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 31 3 319 - 327 1998年02月 [査読有り][通常論文]
     
    Spatio-temporal development of electron swarms in gases is simulated using a propagator method based on a series of one-dimensional spatial moment equations. When the moments up to a sufficient order are calculated, the spatial distribution function of electrons, p(x), can be constructed by an expansion technique using Hermite polynomials and the weights of the Hermite components are represented in terms of the electron diffusion coefficients. ii is found that the higher order Hermite components tend to zero with time; that is, the normalized form of p(x) tends to a Gaussian distribution. A time constant of the relaxation is obtained as the ratio of the second-and third-order diffusion coefficients, D-3(2)/D-L(3). The validity of an empirical approximation in time-of-flight experiments that treats p(x) as a Gaussian distribution is indicated theoretically. It is also found tl-lat the diffusion coefficient is defined as the derivative of a quantity called the cumulant which quantifies the degree of deviation of a statistical distribution from a Gaussian distribution.
  • J Yang, PLG Ventzek, H Sugawara, Y Sakai, K Kitamori, H Tagashira
    JOURNAL OF APPLIED PHYSICS 82 5 2093 - 2105 1997年09月 [査読有り][通常論文]
     
    Models capturing the periodic steady-state behavior of rf capacitively coupled discharges are now commonplace. New plasma sources have been motivated by selectivity, charge-damage mitigation, and general process control needs in plasma processing of electronic materials. These new sources require models that can accurately capture the transient behavior of the plasma source. Such models are not commonplace because the behavior of transport parameters in transients is still not well understood and because the problem is inherently stiff, i.e., widely disparate time scales are important. In this paper, we present the results of an investigation of the simplest type of transient, known as a step disturbance, in a 2 cm gap parallel-plate argon discharge at 1 Torr. As examples, two classes of step transients are considered: step increases in the peak-to-peak (pp) applied voltage (300 to up to 450 V pp) and step decreases (300 to as low as 150 V pp). The resulting transients are interpreted in terms of time scales representative of electron and ion motion in the sheath, ionization dynamics, and neutral transport processes. The possibility of using these transients as a means of process identification is discussed. (C) 1997 American Institute of Physics.
  • 菅原広剛, 田頭博昭, 酒井洋輔
    電気学会論文誌Ａ 117-A 7 749 - 753 1997年06月 [査読有り][通常論文]
  • H Sugawara, H Tagashira, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 30 3 368 - 373 1997年02月 [査読有り][通常論文]
     
    The centroid drift velocity W-r and the mean electron velocity W-v of electron swarms in gases have different values even under the same electric field when electron-impact ionization and/or electron attachment processes are present. A new evaluation technique for W-r is proposed based on moment equations, in which the quantitative difference between W-r and W-v is described as a function defined in velocity space. W-r values are derived from calculations using a propagator method performed in velocity space without the electron distribution in real space. The effects of ionization and attachment on W-r are quantitatively analysed for F-2 and SF6, and the mechanism producing the difference between W-r and W-v is illustrated.
  • H Sugawara, H Tagashira, Y Sakai
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 29 5 1168 - 1174 1996年05月 [査読有り][通常論文]
     
    An explicit expression for the velocity distribution of ion swarms in warm gas is deduced for a constant-collision-frequency symmetric charge exchange model. The ion velocity distribution, the drift velocity and the diffusion coefficient are calculated under DC and AC electric fields. It is found that the ion velocity distribution normalized with respect to the thermal velocity is a function only of dimensionless parameters fixed by physical quantities given as analysis conditions. These dimensionless parameters are potentially good scaling parameters by which different experimental conditions may be compared.
  • H SUGAWARA, Y SAKAI, H TAGASHIRA
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 28 1 61 - 67 1995年01月 [査読有り][通常論文]
     
    Exponential spatial growth of electron swarms in gases under steady-state Townsend conditions may be observed not only in the downstream region of an electron source but also in the upstream region due to backward diffusion. Relations among swarm parameters in the upstream region were deduced and found to have peculiar characteristics. For example, the sign of the average electron velocity V-d for steady-state Townsend conditions depends on whether the gas is electropositive or electronegative. This property is confirmed quantitatively by calculating the electron energy distribution using a propagator method modified for analysis in the upstream region. As an example of the effect of backward diffusion, the decay in the electron number density in front of the anode in a steady-state Townsend experiment between parallel plane electrodes was studied. The decay is caused by the missing electrons absorbed at the anode. The energy distribution of the missing electrons has a non-equilibrium relaxation region and the equilibrium region lies towards the cathode direction. The number density of the missing electrons showed an exponential spatial dependence in the equilibrium region. The relative gradient of the decrease agreed well with the relative density gradient alpha obtained in the upstream region. It was also found that the electron energy distribution in the decaying region tends to that in the upstream equilibrium region.
  • H SUGAWARA, Y SAKAI, H TAGASHIRA
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 27 1 90 - 94 1994年01月 [査読有り][通常論文]
     
    A new simulation technique for obtaining the equilibrium electron energy distributions of electron swarms in gases under steady-state Townsend conditions is proposed based on a propagator method. The electron energy distribution function could be obtained from calculations in only one slab in real space, assuming exponential spatial growth of an electron swarm and similarity of the electron energy distribution at any position. The electron energy distributions were analysed in argon at E/N = 1414 Td for a case of a positive large value of the effective ionization coefficient alpha and in an argon (90%)/fluorine (10%) mixture at E/N = 56.6, 84.8, 113.1 and 141.1 Td for a case of values of alphaBAR crossing zero. The results agreed excellently with those of a Monte Carlo simulation and a Boltzmann equation analysis. This propagator method was found to deduce stably proper solutions for both positive and negative values of alphaBAR. A relation between the present method for a steady-state Townsend condition and a similar method for a pulsed Townsend condition is discussed.
  • H SUGAWARA, Y SAKAI, H TAGASHIRA
    JOURNAL OF PHYSICS D-APPLIED PHYSICS 25 10 1483 - 1487 1992年10月 [査読有り][通常論文]
     
    The position-dependent swarm parameters and the energy distribution F(epsilon, x) of electrons in argon gas have been analysed using a propagator method in which a proper relation between the trajectories of electrons and the cells in an energy epsilon and position x space were considered. The results showed fine structures of damped fluctuation in the swarm parameters, which had been seen from the results of a Monte Carlo method. When initial electrons with zero energy were ejected from the cathode, F(epsilon, x) at E/N = 283 Td could not reach the equilibrium condition even at the distance at which the equilibrium values of mean energy, drift velocity, etc of electrons seemed to have already been obtained.

MISC

書籍等出版物

  • 電気学会技術報告第1488号「放電・プラズマ気相シミュレーション技法」
    放電・プラズマ気相シミュレーション技法調査専門委員会編, 菅原広剛, 調, 委員会委員長 (担当:分担執筆範囲:1.1節「電子エネルギー分布・速度分布計算法１：ボルツマン方程式解析」1.2節「電子エネルギー分布・速度分布計算法２：モンテカルロ法」1.3節「電子エネルギー分布・速度分布計算法３：プロパゲータ法」)
    社団法人電気学会 2020年06月
  • 電気学会技術報告第1185号「放電基礎パラメータと放電応用技術の最前線」
    菅原 広剛 (担当:分担執筆範囲:2.2.1節「プロパゲータ法による電子スオーム解析」、2.2.2節「モンテカルロ法における電子衝突判定の高速化」)
    社団法人電気学会 2010年04月
  • 電気学会技術報告第1061号「窒素の放電プラズマが見せる多様性と最新応用」
    菅原 広剛 (担当:分担執筆範囲:５章１節「N2およびSF6/N2混合ガス中の電子衝突過程と電子エネルギー分布」)
    社団法人電気学会 2006年08月

講演・口頭発表等

  • H. Takahashi, H. Sugawara
    The 11th Asia-Pacific International Symposium on the Basics and Applications of Plasma Technology 2019年12月 ポスター発表
  • H. Sugawara, T. Sato, Y. Nakata
    The 11th Asia-Pacific International Symposium on the Basics and Applications of Plasma Technology 2019年12月 ポスター発表
  • Investigation of mobility spectra obtained by mobility measurement and simulation in different-purity O2 gases  [通常講演]
    Y. Okuyama, K. Ikeda, H. Sugawara
    The 11th Asia-Pacific International Symposium on the Basics and Applications of Plasma Technology 2019年12月 口頭発表（一般）
  • Response of Average Electron Velocity Vector under AC Electric and DC Magnetic Fields in a Constant-Collision-Frequency Model  [通常講演]
    H. Sugawara
    Proc. 72nd Annual Gaseous Electronics Conference 2019年10月
  • Quantification of electron confinement effect under confronting divergent magnetic fields considering mutual electron transfer across separatrix  [通常講演]
    R. Ozawa, H. Sugawara
    Proc. XXXIV International Conference on Phenomena in Ionized Gases, 10th International Conference on Reactive Plasmas, 37th Symposium on Plasma Processing, and 32nd Symposium on Plasma Science for Materials 2019年07月
  • Electron energy gain mechanisms near chamber wall in inductively coupled magnetized plasmas  [通常講演]
    H. Takahashi, H. Sugawara
    Proc. XXXIV International Conference on Phenomena in Ionized Gases, 10th International Conference on Reactive Plasmas, 37th Symposium on Plasma Processing, and 32nd Symposium on Plasma Science for Materials 2019年07月
  • Evaluation of electron drift velocity vector in gas under RF electric and DC magnetic fields  [通常講演]
    T. Sato, H. Sugawara
    Proc. XXXIV International Conference on Phenomena in Ionized Gases, 10th International Conference on Reactive Plasmas, 37th Symposium on Plasma Processing, and 32nd Symposium on Plasma Science for Materials 2019年07月
  • The comparison between experimental and simulation results of negative ion mobility in O2 in relation to reduced electric field intensity  [通常講演]
    K. Ikeda, Y. Okuyama, H. Sugawara
    Proc. XXXIV International Conference on Phenomena in Ionized Gases, 10th International Conference on Reactive Plasmas, 37th Symposium on Plasma Processing, and 32nd Symposium on Plasma Science for Materials 2019年07月
  • 菅原 広剛
    The 19th Asian Conf. on Electrical Discharge 2018年11月 口頭発表（招待・特別） 咸陽，陝西省，中国
  • 菅原 広剛
    XXIV Europhysics Conf. on Atomic and Molecular Physics of Ionized Gases 2018年07月 ポスター発表 Glasgow, Scotland, UK
  • 菅原 広剛
    20th Gaseous Electronics Meeting 2018年06月 口頭発表（招待・特別） Magnetic Island, Townsville, Australia
  • 菅原 広剛
    10th Asia-Pacific Int. Symp. on the Basics and Applications of Plasma Technol. 2017年12月 口頭発表（招待・特別） 中原大学，桃園，台湾
  • 菅原 広剛
    XXXIII Int. Conf. Phenomena in Ionized Gases 2017年07月 ポスター発表 Estoril/Lisboa, Portugal
  • Fundamental study on filter effect of confronting divergent magnetic fields applied to a low-pressure inductively coupled plasma  [招待講演]
    菅原 広剛
    9th Int. Conf. Reactive Plasmas/68th Ann. Gaseous Electronics Conf./33rd Symp. Plasma Processing 2015年10月 口頭発表（招待・特別）
  • Electron conduction under the rectification effect of antiparallel magnetic fields  [招待講演]
    菅原 広剛
    XVI Int. Workshop on Low-Energy Positron and Positronium Physics & XVII Int. Symp. on Electron-Molecule Collisions and Swarms 2015年07月 口頭発表（招待・特別） Maynooth, Ireland
  • 対向発散磁界中電子運動のモンテカルロシミュレーション  [通常講演]
    菅原 広剛
    第31回電離気体現象国際会議 2013年07月 ポスター発表
  • 反平行傾斜磁界間電子伝導路において内向電子流を生じる確率論的散乱過程  [通常講演]
    菅原 広剛
    第21回電離気体中原子分子物理欧州会議 2012年07月 ポスター発表
  • 不平等磁界下の電子輸送  [招待講演]
    菅原 広剛
    原子分子データ応用フォーラムセミナー・核融合科学研究所原子分子データと原子分子過程モデルの活用・普及研究会 2011年12月 口頭発表（招待・特別）
  • ４極磁界・ラジオ周波数電界下における電子輸送のモンテカルロシミュレーション  [通常講演]
    菅原 広剛
    第30回電離気体現象国際会議 2011年08月 ポスター発表
  • Electron conduction under the rectification effect of antiparallel magnetic fields  [招待講演]
    菅原広剛
    XVI Int. Workshop on Low Energy Positron and Positronium Physics & XVII Int. Symp. Electron-Molecule Collisions and Swarms 2011年07月
  • Electron drift along magnetically neutral channels between gradient magnetic fields  [招待講演]
    菅原 広剛
    15th Gaseous Electronics Meeting 2008年02月 口頭発表（招待・特別） Murramarang Resort, Australia
  • Electron acceleration under short intense impulse electric fields  [招待講演]
    菅原 広剛
    XVth Int. Conf. on Gas Discharges and their Applications 2004年09月 口頭発表（招待・特別） Toulouse, France
  • Control of selectivity in electron-molecule reactions by impulse field electron acceleration  [招待講演]
    菅原 広剛
    Int. Symp. on Electron-Molecule Collisions and Swarms 1999年07月 口頭発表（招待・特別） The Waterfront Town, Tokyo, Japan
  • Electron behavior under characteristic magnetic fields applied to inductively coupled plasmas for control of charged particle transport  [招待講演]
    菅原 広剛
    第18回低エネルギー陽電子およびポジトロニウム物理国際ワークショップ＆第19回電子分子衝突および電子群国際シンポジウム 口頭発表（招待・特別）

担当経験のある科目（授業）

  • 計算機命令実習北海道職業能力開発短期大学校
  • 集積プロセス学特論北海道大学大学院
  • 電気工学実験／電子工学実験／電子情報工学実験北海道大学
  • 電子工学演習北海道大学
  • ユーザ・コンピューティング札幌大学
  • 情報処理概論札幌大学
  • 情報処理札幌大学
  • エネルギー応用工学北海道工業大学
  • プログラム言語実習北海道職業能力開発短期大学校
  • コンピュータ工学北海道職業能力開発短期大学校
  • ソフトウエア工学北海道職業能力開発短期大学校
  • 情報学北海道大学
  • 集積プロセス工学北海道大学
  • 数値解析とシミュレーション基礎北海道大学
  • プラズマエレクトロニクス北海道大学
  • 電気回路北海道大学

所属学協会

  • 応用物理学会   日本工学教育協会   社団法人電気学会   The Institute of Electrical Engineers of Japan   

共同研究・競争的資金等の研究課題

  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2024年04月 -2027年03月 
    代表者 : 菅原 広剛
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2019年04月 -2022年03月 
    代表者 : 菅原 広剛
     
    基盤研究(C)一般、課題番号19K03780
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2016年04月 -2019年03月 
    代表者 : 菅原 広剛
     
    基盤研究(C)一般、課題番号16K05626
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2013年04月 -2016年03月 
    代表者 : 菅原 広剛
     
    基盤研究(C)一般、課題番号25400528
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2010年04月 -2015年03月 
    代表者 : 村山 明宏, 末岡 和久, 菅原 広剛
     
    電子のスピン緩和が抑制されるIII-V族化合物半導体量子ドットを用いて、電子スピン情報をレーザー光の円偏光特性に変換する量子ドットスピンレーザーについて研究した。 量子ドットに対するスピン注入のダイナミクスとスピン緩和、スピン保持に優れた量子ドットの作製等の研究を行い、スピン注入の主要な阻害要因としてパウリブロッキングの影響を明らかにした。そして、ドットの高密度化や二次元電子系とのトンネル効果により、高効率で超高速のスピン注入が可能になることを示した。また、レーザー媒質となる量子ドット活性層と金属強磁性体スピン電極を持つ発光ダイオードを作製し、電子スピンの電流注入による円偏光特性を得た。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2010年 -2012年 
    代表者 : 菅原 広剛
     
    半導体集積回路微細加工用の磁界下環状放電プラズマの計算機解析を行った。環状放電部から伸びる分界線(磁力線の流れが変わる境界)沿いに高エネルギー電子が輸送され、ガスを分解して化学活性種を生成することを解明した。半導体基板上では分界線の到達点に多く化学活性種が入射することを再現し、過去の実験結果を説明した。磁界を変化させ分界線を動かすことで基板上の広い面積を均一に表面処理する技法を提案した。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2007年 -2009年 
    代表者 : 菅原 広剛, 須田 善行
     
    プラズマでフッ化炭素C_8F_<18>分子を分解し電力機器や集積回路の電気絶縁用非晶質フッ化炭素膜を合成した。プラズマ電力が低いと原料分解が進まず分子鎖間の結合が疎な膜ができた。疎膜は従来の高電力で合成した結合が密な膜より耐熱性が低いが低誘電率となり電力損失や電気信号伝達遅延の低減が期待される。また疎膜上に密膜を合成する際の熱による膜の波状化現象を見出した。波状膜下にできる隙間も低誘電率化に有効と期待される。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2007年 -2008年 
    代表者 : 須田 善行, 菅原 広剛, 高山 純一
     
    (1) カーボンナノチューブ(CNT)の原料ガスであるCH_4にCO_2ガスを添加したところ,多層CNTの層数・内径・外径が増加した。また,CO_2添加により単層CNTの直径分布を太い方へ変化できた。(2)CNT成長領域から60cm程度の距離にプラズマを発生させることで,Si基板上に堆積したAl_2O_3/Fe/Al_2O_3の3層構造触媒から欠陥のほとんどない高品質な単層CNTが合成できた。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2005年 -2007年 
    代表者 : 酒井 洋輔, 菅原 広剛, 須田 善行
     
    SF_6ガスは絶縁材料として極めて高い性能を有しているが、地球温暖化係数が大きく、その使用が制限されている。これに代わる一つの絶縁方式として、電源導体部に非晶質フッ化炭素(a-C:F)膜をプラズマ化学気相成長(PCVD)法により被覆し導体部からの二次電子放出を抑制、その結果大気(窒素と酸素)によるガス絶縁に於いてもSF_6ガスと同程度(またはそれ以上)の性能を得ることのできる、新複合絶縁方式を提案した。 結果を要約すると、1)数百ナノメータ程度のa-C:F膜が金属(銅とAl)球表面に物理化学的に安定に堆積できる、2)膜質を表面分析装置(SEM, XPS及びFTIR、等)で評価すると物理化学特性がテトラフルオロエチレン(低誘電率と高絶縁耐力)に近いものが得られる、3)本薄膜を堆積した球電極間の空気中パッシェン特性(pd(気圧と電極間隙の積)領域における絶縁破壊電圧)はSF_6ガスの特性と同程度、4)a-C:F膜表面の特性は、バッファガス(Ar, He)に依存する、等ことを示した。また、簿膜の一部が絶縁破壊したときを想定し、既存の膜上に状にa-C:F)膜をプラズマ堆積した場合、5)本膜の二重、三重層境界を確認するとともに、三重層膜表面では前層のストレス部に対応して現れると思われるに規則的模様の現れることが観測された。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2004年 -2006年 
    代表者 : 菅原 広剛
     
    本研究は、集積回路製造工程で排出されるフルオロカーボンであるCF_4地球温暖化ガスを非平衡プラズマにより効率的に分解する方法として選択的分解反応促進が期待されるプラズマへの衝撃電界重畳を提案し解析するものである。衝撃電界により電子した電子の衝突過程の計算機解析による定量的評価を中心として、以下1〜3の結果を得た。 1.衝撃電界下の電子の挙動解析とエネルギー効率の評価: 衝撃電界で与える電子加速量を変化させ、共にCF_4分解に寄与する中性解離と電離性解離のエネルギー消費割合をモンテカルロ法で計算・比較した。加速量増加に伴い閾値の大きな電離性解離に費やされるエネルギー割合が単調増加し、中性解離の割合は飽和傾向を示した。プラズマ維持と解離促進を同時に図るためには電離性解離の閾値を超える高エネルギー領域にまで電子を加速することが必要で、この点で衝撃電界電子加速が有効と考えられる。 2.CF_4ラジオ周波数プラズマ解析系構築と新高速衝突判定法の提案: 電極を含む系でCF_4プラズマ中の解離生成種の生成頻度と空間分布とを統計的に安定に扱う手法を提案した。また計算負荷となっていたCF_4分子と電子との衝突判定において、無衝突となる場合を効率よく検出し電子衝突断面積計算を省く新高速判定法を提案した。これにより電子追跡に要する計算負荷を著しく軽減させ、追跡電子数増や長時間計算、また混合ガスへの対応を容易にした。 3.CF_4中電子輸送過程: CF_4中の電子加速時初期緩和過程に見られる空間分布の複数ピーク形状と電子輸送係数の緩和過程を解析した。研究者間で想定が異なる電子衝突断面積構成により空間分布、電子輸送係数とも差異が生じた。衝撃電界電子加速の制御性評価の際は特に、加速途中で電子の集団性を失わせる低・中エネルギー部の電子衝突断面積構成への依存性が大きいと考えられる。 主な結果は研究発表欄の論文として公表した。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2002年 -2003年 
    代表者 : 菅原 広剛
     
    地球温暖化防止のためのSF_6絶縁ガス使用量低減に向け、SF_6をN_2他で希釈した際に現われる絶縁特性相乗効果発生機構を混合ガス中の電子エネルギー損失過程と電子捕獲機構の面から解析した。 平成15年度は、電子のエネルギー損失過程および衝突間の電子挙動に関する14年度の解析結果を踏まえ、既に構築済みの計算機プログラムを用いて、電子衝突断面積の組み合わせの効果を主にSF_6/N_2、SF_6/CO_2、SF_6/c-C_4F_8の3つの混合ガスの比較を通じて解析した。解析結果は次のようにまとめられる。 SF_6/N_2混合ガス中では、電子はN_2との衝突によりエネルギーを失って低エネルギーとなり、2〜3eV付辺に存在するN_2の弾性衝突および振動励起衝突の断面積ピークに起因するバリア効果(電子エネルギー上昇を妨げる)によって低エネルギー状態に保たれその間にSF_6により捕獲されることが電子衝突履歴の追跡により確認できた。 SF_6/CO_2でもCO_2のバリア効果が見られるが、CO_2はN_2に比べバリアが低く高エネルギー寄りに位置することから電子を低エネルギー域に閉じ込め付着を促す働きはSF_6/N_2ほど強くない。N_2とCO_2についてはバリア効果の強弱とSF_6との相乗効果の程度は大小関係が一致しており、この種の混合ガスの相乗効果はバリア効果によるものと推測された。 SF_6/c-C_4F_8混合ガスにおいては、c-C_4F_8の電子励起衝突断面積の立ち上がりが急なことが励起の閾値付近での衝突を起こり易くし、励起衝突直後の電子エネルギーを低くする結果をもたらす。SF_6の電子付着は低エネルギーほど起こり易いため、閾値近くの励起衝突が直接電子付着に結びつくためと考えられる。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2001年 -2003年 
    代表者 : 酒井 洋輔, 須田 善行, 赤澤 正道, 菅原 広剛, 中島 昌俊, BRATESCU A. Maria
     
    本研究は地球温暖化ガスの一つとして使用量の削減が迫られているSF_6ガスに替わり窒素や酸素あるいは合成空気と導体部に低誘電率・高絶縁耐力をもつアモルファスC_xF_y(a-C:F)膜を堆積し二次電子放出を抑え、これまで同様あるいはそれ以上の絶縁耐力を得ようとすることである。絶縁膜原料には、C_7F_<16>とC_8F_<18>(常温では液体で、蒸気圧は数〜数十Torr)を用い、RFプラズマCVD法により作成し、これの化学的・電気的特性を評価した。成果は以下のとおりである。 1.Siならびにアルミニウム基板上への堆積速度は、100〜200nm/min以上で、在来原料のCF_4やC_2F_6に比べ数十倍以上の高堆積速度が得られた。 2.堆積a-C:Fは高密度C-CとC-F結合から成り、誘電率(ε_r【approximately equal】2)の低い熱特性にも優れた絶縁膜であった。絶縁耐力は1μm厚で2.7MV/cm。誘電体損はポリテトラクロロエチレンよりは大きかった。 3.本a-C:FをAl電極上に堆積したものを用いてN_2、ArとHeガスのパッセン曲線を測定したが、破壊電圧V_sはAl電極の場合に比べて低pd(気圧x電極間隙)領域で3倍程度上昇した。また、pd<20Torr・cmでは、SF_6のV_sよりも向上することが分かった。 4.この縁耐力向上の原因は、二次電子放出係数が、金属のものに比べ圧倒的に小さく、10^<-5>オーダの小さな値を与えたことによることを明らかにした。 5.プラズマ相で生じた分解種とa-C:F膜の化学組成の関係を調べたが、プラズマ条件が膜質に与える影響が必ずしも明らかになったとは言えず、今後の課題として残った。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 2000年 -2001年 
    代表者 : 菅原 広剛
     
    C_7F_<16>を材料としたプラズマCVD法でアモルファスフルオロカーボン(CF)膜を堆積した。XPSでF/C組成比を、静電容量測定とエリプソメータによる膜厚測定から比誘電率ε_rを、正極性針対平板電極系絶縁破壊試験で絶縁耐力をそれぞれ測定し、熱処理(100〜300℃)後の膜厚減少分から耐熱性を評価した。これらの結果から(1)材料のガス圧依存性(30〜80Pa)、(2)He混合(0〜30Pa)の影響、(3)酸素混合(2〜4Pa)の影響を検討した。 (1)では低圧ほどF/C比が上昇しε_rは2.1程度まで下がったが、同時に堆積率も低下し耐熱性試験時の膜厚減少幅も大きくなった。絶縁耐力は2MV/cm以上であった。材料ガス圧による膜質制御では耐熱性向上と誘電率低減の両立は難しかった。類似の巨大分子材料C_8F_<18>を用いて堆積したCF膜でも比誘電率は2.5以下、絶縁耐力は1MV/cm以上であった。 (2)ではHe増加に伴いε_rが低下し(1)同様の傾向が見られた。加えて、堆積中のプラズマ状態が安定して膜質のばらつきが減り再現性が向上した。 (3)ではCF膜からのC原子引き抜きによるF/C比増加即ちε_r低下を期待したが、酸素活性種のエッチング効果が強く数Paの混合で堆積率が著しく低下し、6Pa以上では膜堆積できなくなった。2〜4Paにおける堆積膜の膜質はばらつきが大きいものの、中には低いε_rが得られた例もあり、多孔質構造のような空疎なCF膜が得られた可能性が示唆される。 続いて、絶縁破壊部再生を模擬し、CF膜で被覆済みの面と基板が露出した面が接する個所へのCF膜再堆積を試みた。膜厚の時間変化から、再堆積開始時は露出面への堆積の方が速やかに進む傾向が見られたが、再堆積が進行した後は、被覆面と露出面に堆積したCF膜の厚さは同程度になった。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1998年 -1999年 
    代表者 : 菅原 広剛
     
    RF(13.56MHz)非平衡放電プラズマに衝撃電界を重畳することによりプラズマ中の電子を短時間のうちに揃って加速し、単一エネルギー的な電子群を生成してフロン12の分解反応の選択性、効率の向上とバルク領域での大体積分解反応の促進を図る技法を流体コードを用いて計算機解析した。 衝撃電界を印加するとプラズマシースの厚さが変化してこれを吸収してしまうため、衝撃電界がバルク領域に到達する時間は限られるが、誘電緩和時間の見積りや数値解析の結果、数ns程度の間は衝撃電界を維持でき、電子エネルギーをフロン12解離のしきい値(中性解離7.2eV、解離性電離12.2eV以上)まで高めることができた。フロン12は高エネルギー領域で中性解離・解離性電離、低エネルギー領域で解離性付着の衝突断面面積を持ち、これらは共に分解反応に寄与する。しかし低エネルギー領域での電子付着が大きいと放電維持が困難になることから、解離性電離を選択的に促進し分解処理と放電維持を同時に図るのが効果的である。 バルク領域へ到達する衝撃電界の大きさは衝撃電界の立ち上がり時間に対して敏感に変化し、衝撃電界の立ち上がりが急峻な程バルクに到達する電界も大きいことが確認された。強度500V/cm、印加時間5nsの衝撃電界について立ち上がり時間を2nsから0.25nsに縮めると、分解処理量は35%程度増加した。この程度の立ち上がり時間は数kV規模の高電圧パルス発生装置で実現しているが、実際の回路では容量と抵抗を小さくしてプラズマ中に大きな変位電流が流れるようにする必要がある。衝撃電界を逆極性で印加しても同程度の電子加速効果が得られたことから考えると、衝撃電界印加は任意の方向で構わないものと思われる。したがって衝撃電界印加のための電極をRF放電駆動用電極とは別に設けて回路定数設定の自由度を高めることも可能と考えられる。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1997年 -1999年 
    代表者 : 酒井 洋輔, 河野 仁志, 五十嵐 龍志, 菅原 広剛, VENTZEK Pete
     
    大気汚染ガスとして問題となっている窒素酸化物と、現在洗浄や電気絶縁に用いられ、その使用後の処理が求められているフルオロカーボン(C_7F_<16>,(C_3F_7)_3N/(C_4F_4)_3N,C_8F_8/C_8F_<16>O)等を直流とRF放電プラズマならびに真空紫外線照射のハイブリッド方式により、分解し大気にもどす、または固化無害化することを目的にその基礎研究を行った。以下に、結果を要約する。 1)1〜3Torrの二酸化窒素ガスは直流放電プラズマにより90%以上窒素と酸素に分解したことを、圧力上昇 (初期封入圧力の約1.5倍に増加)、プラズマ発光スペクトルの変化(二酸化窒素ならびに一酸化窒素の発光スペクトルから窒素スペクトルヘ変化)、ならびに分解生成物の質量分析とガスクロマトグラフ法から確認した。 2)二酸化窒素ガスはキセノンエキシマランプから発する紫外線(波長172nm)を照射した場合も分解すること、また1)の実験に紫外線照射を重畳した場合には、放電プラズマ分解ならびに紫外線分解を別々に行った結果の和よりも窒素と酸素への分解速度と分解率が加速度的に向上する相乗効果が見られた。 3)〜10Torrのフルオロカーボン蒸気をプラズマ分解することにより、電極表面に薄膜の堆積することが確認された。また、プラズマ相においては、分解生成種CF、CF_2、CF_3、C_2F_3、C_2F_5、C_3F_3、C_3F_5、C_3F_6等が質量分析装置で検出され、これらと薄膜の堆積特性の関係が検討された。 4)本研究で得られた薄膜は、意外にも低誘電率で高絶縁耐力があり、半導体集積回路の絶縁層や電気機器の絶縁補助に利用できる可能性を示した。本応用に関しては現在も研究を続けている。 汚染粒子分解用放電プラズマの性能向上に関しては、モデリングと励起種の計測から種々の知見を得た。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1995年 -1996年 
    代表者 : 酒井 洋輔, 五十嵐 龍志, 菅原 広剛, VENTZEK pete, 田頭 博昭
     
    Xeガスエキシマランプの性能向上を目的として、バリア放電プラズマの特性とフィラメント状放電への移行の理解、光出力と効率を決定する放電プラズマ因子、ならびに、これらの限界を検討した。結果を以下にまとめる。 1.Xeガスバリア放電プラズマにおいては、印加電圧周波数が低いほどバリアへの堆積電荷量が増すこと、エキシマは電極近傍で効果的に生成されること、印加電圧波形の立ち上がりが速いもの程エキシマ生成効率が高くなることがわかった。また、エキシマレーザ放電(ArFならびにXeC1)におけるフィラメント状放電が局所的なガス温度上昇から発生することが、一般的に言えることが示唆された。 2.バリア放電プラズマをレーザ分光法による診断を行った場合、信号とプラズマ特性の関係を理解するために行ったモデリングにおいては、通常診断に用いられている光子密度では電子のエネルギー分布関数に有為な攪乱を与えないが、準安定励起原子密度は大きく変わることを明らかにした。 3.DCから数百kHzの電圧を印加した場合のXe/Ne混合ガスプラズマのレーザ分光診断実験においては、Xe準安定励起原子密度はDCでは陰極近傍、ACでは両電極近傍で大きくり、モデリングと同等の結論が得られた。 以上の結果から、ランプとしてはエキシマを効果的に生成する電極シース部を利用する設計が望ましいという結論を得た。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1995年 -1995年 
    代表者 : 田頭 博昭, 中神 芳武, 菅原 広剛, 下妻 光夫, 酒井 洋輔
     
    本研究において、フロンを含有する絶縁材料の使用を排し、「地球にやさしい」新しい高電圧電力機器絶縁方式を確立するための基礎研究を行った。空気と微粉体の混合気体を用いた、フロンを含まない材料による気中ギャップの絶縁破壊電圧を、半球対半球ギャップ、あるいは球対球ギャップについて、実験的に研究した。直流、50Hz交流、標準雷インパルスの3種類の印加電圧波形について、データをとった。粉体材料としては、SiO_2、Si_3N_4、TiO_2、CeO_2、Al_2O_3、CuOを用い、粒径はそれぞれ10μmを平均として分布している。50Hz交流電圧印加時は、微粉体の効果は、ほぼここに記した順序でSiO_2が一番空気との比較において破壊電圧の上昇が大きい。2mmギャップにおいて、上昇率が最も大きく、約70%に達した。一方、Al_2O_3やCuOでは、上昇は極めてわずかであった。直流では一般に交流よりも上昇はわずかで、最も大きかったSiO_2でギャップ長3mmでも上昇率は15%程度であった。インパルス電圧印加の場合はパウダの量を準備する都合でCeO_2を除いた5種類について実験を行ったが、いずれの場合も20〜40%程度の破壊電圧の上昇が見られた。これらとは別に、直流および交流のみについCaCO_3パウダ、及びCa(OH)_2パウダに対しても実験を行った。比較的日常的に我々の周辺にある微粉体に興味がひかれたが、いずれも電圧上昇の効果は小さい上、Ca(OH)_2では破壊電圧は却って低下する結果が得られた。SiO_2パウダで典型的に現れるが、直流電圧印加時にパウダが電極(正電極)に付着する現象が見られる。これはパウダが電子を付着する機能を持つことを示し、破壊電圧の上昇を説明することができる。他の説明としては、パウダの電極への付着による電極の誘電体コーティングによる破壊電圧の上昇がある。現実にはこの両者が破壊電圧の上昇に同時に作用しているものと考えられる。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1994年 -1994年 
    代表者 : 田頭 博昭, 横澤 亜由美, 菅原 広剛, ヴェンシェック ピーター, 酒井 洋輔
     
    本研究は当初Si_3N_4薄膜を形成するためのラジカルを発生するRFプラズマのシミュレーション解析とその薄膜形成過程のシミュレーションに主としてあてられる予定であったが、Si_3N_4ならびにそのプリカーサラジカルに関する研究が不十分のため、急遽方針を転換して、次のような研究を行ない、成果を発表した。 1)アモルファスシリコン膜堆積についてダイヤモンド格子も用いたモデリングを行い、これによってスピン密度の研究を行なった。その結果、スピン密度を与えるもととなるダングリングボンドはボイド中にあるという結論を得たが、これはおそらく相当に堆積速度が大きく十分密な膜となっていないときの結果であろうと思われる。最近のきわめて良質な膜では、ボイドとダングリングボンドは相関が薄いと報告されており、この点の検討が、今後必要であることを明らかにした。 2)電荷転移と弾性衝突がある気体中のイオンの速度分布、スオームパラメータのモンテカルロ法による研究を行ない、電荷転移によるピークをもつ速度分布の特異な形状が弾性衝突により緩和されることを明らかにした。M.Kushnerが与えた周波数をもとにモノシランガス中のシリルイオンの拡散係数を求めた。 3)RF非平衡プラズマの制御を目的として、トランスジェント特性をもとめ、制御を行なうときの基礎的データを得た。また、比例オンオフ制御による制御性をシミュレーション研究し、制御可能性を見いだした。 4)来年度研究予定のTEOS中の電離係数を、北見工大吉田研究室と協力して測定し、初めてその値を与えた。これは、TEOSガス中のプラズマのシミュレーションの基礎を与えるものである。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1993年 -1994年 
    代表者 : 田頭 博昭, ピーター ヴェンツェック, 菅原 広剛, 下妻 光夫, 酒井 洋輔
     
    本研究はミクロンオーダーの直径を持つファインパウダを絶縁空間に浮遊させることにより気体の電気的絶縁破壊電圧を向上し、かつその種類を選ぶことにより、無公害的にこれを行ない、高電圧機器の電気的絶縁をハロゲンフリーな気体絶縁とする可能性を実験的に見いだすことを目的としている。以下に研究成果を報告する。 使用したパウダはSiO_2、TiO_2及びSi_3N_4で、電極間に印加した電圧は直流(DC)、50Hz交流(AC)及び標準雷インパルス(LI)である。LI印加時は、50%フラッシ-オーバ電圧V_<50>を昇降法により求めている。紫外線照射は行なっていない。 まず、SiO_2パウダで半球電極間にLIを印加した時、V_<50>はギャップ長1mmで空気の約1.6倍に上昇した。DC、AC電圧印加時にあったパウダの電極への付着はなかった。 TiO_2パウダではLI印加時で、SiO_2パウダ使用時とほぼ同程度で、ギャップ長1mmで空気での値の約1.5倍となった。しかし、DC及びAC電圧印加時では、絶縁破壊電圧の大幅な上昇はなかった。ACでは、約1.1〜1.2倍の上昇が見られたが、DCでは空気の絶縁破壊電圧とほぼ同じであった。いづれの電圧でも、SiO_2パウダの実験で見られたパウダの電極への付着は観察されなかった。 Si_3N_4パウダを用いて同様に実験を行なった。LI電圧印加時では、そのV_<50>はギャップ長1mmで空気の約1.6倍となった。AC及びDCでは、絶縁破壊電圧の大きな上昇は見られない。また、パウダの電極への付着も観察されていない。 LI電圧印加時では、いずれのパウダでも絶縁破壊電圧は空気より上昇した。しかし、AC及びDC電圧を印加した場合、SiO_2パウダを使用したときが最もパウダの効果が現れている。これは、パウダが帯電し電極に付着するということに関係しているものと思われる。各パウダでの、帯電の状況及び電極に付着した場合の付着量を調べ、絶縁破壊電圧の上昇の要因の検討を行なっている。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1993年 -1993年 
    代表者 : 菅原 広剛
     
    本研究は大気汚染の原因の1つとして考えられている窒素酸化物のコロナ放電による分解をコンピュータシミュレーションにより解析し、分解反応機構の解明ならびに分解反応の高効率化を図ることを目的とする。 コロナ放電のシミュレーションに必要な窒素酸化物分子の衝突断面積の推定は、定常タウンゼント法による電離係数の測定結果に、ボルツマン方程式解析による結果を合わせるよう断面積修正を繰返して行なう計画であったが、実際にボルツマン方程式解析を試みたところ、窒素酸化物のような電子付着性のガスの場合には計算が不安定となる場合があることがわかった。そこで、より安定な解析手法として、研究代表者らが開発しコロナ放電の解析に用いる計画であったプロパゲ-タ法を衝突断面積の推定にも導入することとした。 プロパゲ-タ法を定常タウンゼント法に適用する場合、従来の適用方法では電子の実空間における位置およびエネルギーの両方について空間分割を行なう必要があり、膨大な計算機容量を必要としていた。本研究ではこの点を踏まえ、エネルギーに関する空間分割のみで定常タウンゼント状態の解を求めることができる手法を新たに開発した。新手法の確立に伴って解析に必要な計算機容量が大幅に軽減されたため、従来は大型計算機で行なっていたシミュレーションが本年度購入備品の超高速アクセラレータボードをパーソナルコンピュータに搭載した計算機系で実行可能となった。 本研究で開発した手法は、解析に必要な計算機容量が従来の手法に比較して格段に小さい。また、本研究の対象である窒素酸化物のように、従来の手法では解析困難な場合があるとされている電子付着を含む条件にも対応可能な手法である。その成果は、研究発表覧に示すように英国物理学会誌に論文として掲載された。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1993年 -1993年 
    代表者 : 田頭 博昭, 横澤 亜由美, 菅原 広剛, 酒井 洋輔
     
    1.平成5年度の研究目的は、モノシラン中の非平衡rfプラズマによる水素化アモルファスシリコン堆積時のラジカル発生、輸送、表面反応と膜堆積過程におけるラジカルの挙動を明らかにすることである。 2.平成5年度の研究成果は、以下の(1)〜(3)に要約される。(1)ラジカル発生のシミュレーション;比較的小規模な計算でラジカル発生量を求めることができる連続の式に基づくモデルをつくり、電極間のラジカル発生レート分布を求めた。モンテカルロモデルから得たラジカル発生レートとよい一致が得られた。(2)イオンの拡散係数と移動度の導出;上記のシミュレーション支援のため文献値がないモノシランイオンのドリフト速度および拡散係数、平均エネルギーをモンテカルロ法で求めた。衝突断面積は文献値を用いている。拡散係数は1Torr、E/N=1414Tdで約2500cm^2/sであった。(3)ラジカルによる膜堆積過程のシミュレーション;シリルラジカル等による膜堆積過程を解明するため、堆積サイトがダイヤモンド構造をなしているモデルを用いてモンテカルロ法で、堆積率、水素含有量、スピン密度、膜密度を求めた。従来用いていた立方格子モデルと比較して、膜密度および堆積率に関しては実験的に得られた値(2.1gcm^<-3>)あるいは傾向(基板温度にほとんど依存しない)とほぼ一致した。また水素含有量に関しては、基板温度200〜400°Cにおいて実験値と一致した値が得られたが、基板温度を500°C以上に加熱した場合、水素含有量がほぼ10%一定となり、実験値から得られる傾向より大きな値となった。これは隣接するH終端ボンドとダングリングボンドの結合に必要な活性化エネルギーを大きく見積りすぎたためと思われる。スピン密度に関しては、ラジカルとサイトの反応式の多くがダングリングボンドをつくり出すため、膜中にこれが残る可能性が高く、このため実験値とは大きく異なった値が得られたものと思われる。
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1992年 -1993年 
    代表者 : 酒井 洋輔, 菅原 広剛, 田畑 昌祥, 田頭 博昭
     
    本計画は、微量のNO_xを濃縮し、これを放電ならびに紫外線により効果的に分解除去することを目的とした基礎研究を行うことである。濃縮方法としては、Light-Induced Drift(LID)効果を用いる。分解とLID実験は並行して行っている。本研究の成果を以下に記す。 1)圧力1〜7TorrのNO_2ガスのグロー放電下でNO_2の80〜90%をN_2とO_2に分解できた。分解生成物は、NOを経てN_2とO_2になることが分かった。 2)放電内の紫外線がNO_2の分解に及ぼす効果ならびに紫外線による分解の可能性を検討するために、NO_2の紫外線照射による分解を試みた。放電分解の場合と同様にN_2とO_2に分解されることが分かったが、その過程は放電による分解過程と比べて異なることが示唆された。 3)NO_2ガスの電離係数、二次電離係数は空気のそれらより数十%小さかったが、電子付着係数はO_2のそれより大きな値をとった。NO_2/He混合ガスの電離係数は、100
  • 日本学術振興会：科学研究費助成事業
    研究期間 : 1992年 -1992年 
    代表者 : 田頭 博昭, 菅原 広剛, 下妻 光夫, 酒井 洋輔
     
    本研究では、粒径の異なる4種類のSiO_2ファインパウダ(12μm、36μm50μm、150μm)ならびに異なる材質のパウダとしてSi_3N_4ファインパウダを空気中に浮遊させ、絶縁破壊電圧の測定を行いその諸特性について検訂した。また実験の際、質量の大きいパウダはファンによる浮遊が困難のため、電極間にパウダを散布する装置を作成し、2つの装置を用いて測定を行なった。結果を以下に示す。 (1) 12μmパウダを使用し、ファンを作動させた状態の交流絶縁破壊において、電極間にパウダのブリッジが形成されなくても絶縁破壊電圧は上昇した。この時、両電極に付着したパウダが電極間に流れこむのがみられ、これが絶縁破壊電圧の上昇に寄与しているものと思われる。 36μmパウダで同様の実験を行なったが、12μmパウダで確認されたほどの大きな絶縁破壊電圧の上昇はみられていない。 (2) 散布装置を用いた4種類の粒径の異なるSiO_2パウダに対する測定では、どの粒径のパウダも絶縁破壊電圧の大きな上昇は見られなかったが、粒径の大きいパウダと、小さいパウダではわずかであるが差があると思われる。また、同じ粒径でありながらファンによる結果(12μm、36μm)と差があるのは、電極付近にあるパウダの量の違いが影響しているものと考えられる。 (3) Si_3N_4パウダによる絶縁破壊実験ではパウダが重いため浮遊している量が少なく、絶縁破壊電圧の上昇はあまり見られない。散布の装置に替えて同様に測定を行なっても大きな上昇はみられていない。 (4) 直流での測定において、パウダが陽極側に多量に付着しており、パウダ自身は負に帯電しているものと考えられる。
  • Modelling and Simulation of Weakly Ionized Plasmas

Copyright © MEDIA FUSION Co.,Ltd. All rights reserved.