光励起下における繊維強化異方性半導体半空間における磁気光熱弾性波伝播への回転効果

サマリー

光励起下で異方性ファイバー強化半導体における回転磁気光熱弾性波の伝播を解析モデルとして提示します。この手法は、高度半導体材料の温度、キャリア密度、変位、応力分布の解析に利用できます。

要約

本研究は、異方性繊維強化半導体半空間の結合した磁気光熱弾性振る舞動を、光励起と磁場を加えた状態を調査します。熱的、弾性的、キャリア密度、電磁的、回転の効果を含む支配方程式は、統一された多物理的枠組みの中で定式化されています。正規モード法は、結合した偏微分方程式を常微分方程式系に変換し、その後一階行列微分系として書き換えるために用いられます。次に、結合問題の解析的解を構築するために固有値ベースの定式化が開発されます。その結果得られる固有値問題、境界条件系、場の量は数値的に評価されます。温度、キャリア密度、変位成分、応力分布に関する数値結果を提示し、回転パラメータと磁場強度が結合した物理場に与える影響を調べます。結果は、異方性半導体媒体内の波動伝搬、減衰特性、場の分布を変化させる回転および電磁相互作用の重要な役割を示しています。シリコンベースの半導体媒体に対して数値計算が行われました。

概要

光熱弾性は、レーザー照射や光熱励起を受ける半導体材料における熱、機械、光学、キャリア密度の相互作用を記述する能力から、活発な研究分野となっています。これらの結合現象は、現代のオプトエレクトロニクスデバイス、半導体技術、レーザー処理システム、マイクロエレクトロニクスの応用において重要な役割を果たしています。Saeedら1は温度依存性を持つ半導体固体における波動伝播^を調査し、Abbasら2は円筒形キャビティと可変熱伝導率を含む半導体における光熱^{相互作用を調査}しました。その後、Ihtisham Ullahらは、半導体媒体中の反射弾性波に対する可変熱伝導率とレーザーパルス^{の影響を解析}しました。さらに、ヤダヴ⁴は二温度多相遅延熱弾性フレームワークの下で拡散半導体における磁気光熱プラズマ波の伝播を研究しました。より最近では、Luteら^5およびLuteら⁶が光熱熱発生、メモリ効果、非局所相互作用を取り入れた先進的な光熱弾性半導体モデルを開発し、AlaofiとEl-Dali⁷は半導体材料向けの多物理的な光熱熱湿弾性の定式化を提案しました。さらに、位相遅延理論や非局所的定式化に基づく複数の一般化された光熱弾性モデルが、半導体^8,9,10,11における結合熱およびキャリア輸送現象の記述を改善することが報告されています。

軽量で高性能な工学構造物の需要増加は、異方性および繊維強化材料に関する広範な研究を促進しています。優れた機械的特性と耐熱性の向上により、繊維強化複合材は航空宇宙、土木、電子工学の分野で広く利用されています。Khanらは温度依存型繊維強化複合材料の統一マイクロポーラー熱弾性^{モデルを開発し、} 微細構造パラメータが波の伝播特性に与える有意な影響を示しました。Pitarresiら¹³ は繊維強化プラスチックの熱弾性応答を調査し、材料の異質性が熱応力分布に与える役割を強調しました。Escalante-Solísら¹⁴ は繊維曲率が強化複合材料の微細機械的挙動に与える影響を調査し、Abo-Dahabら¹⁵ は繊維強化熱弾性媒体における波動反射現象を研究しました。より最近では、PurkaitとKanoria¹⁶ がパルスレーザー励起を受けた回転繊維強化磁気熱弾性媒体を解析し、Akaiら¹⁷ およびQuinlanら¹⁸ は異なる荷重条件下での繊維強化複合構造における熱弾性応答を調査しました。これらの研究は、強化特性と異方性が熱輸送、応力濃度、波の伝搬挙動に大きな影響を与えることを確認しました。それにもかかわらず、利用可能な研究の多くは熱弾性複合材料に集中しており、光熱弾性半導体結合とキャリア輸送機構を同時に考慮していませんでした。

磁場と熱弾性半導体媒体の相互作用も、電磁デバイスや導電性材料における重要性から大きな注目を集めています。^Yadav19 は拡散効果を伴う回転する直交対立磁気熱弾性半空間における波動反射を調査し、Selvamaniらは非局所的分画ナノビームにおける磁気弾性波の伝^播を研究 しました。Abouelregalら²¹ は、記憶依存効果を持つ磁化多孔質構造におけるレーザー熱磁気刺激の結合を研究し、一方Chandelら²² はMoore-Gibson-Thompson熱伝導に基づく非局所的な磁気熱弾性フレームワークを提案しました。半導体の文脈では、Abouelregal²³ は磁場下にある回転半導体半空間の修正分数光熱弾性モデルを開発し、Sur²⁴ は遺伝的特性を示す媒体における磁気・光・熱弾性相互作用を調査しました。より最近では、Saidiら²⁵ およびRashidら²⁶ が一般化回転半導体媒体中の磁気光熱弾性乱乱を解析し、磁気強度が結合波の伝播に与える有意な影響を示しました。同様の観測は、先進的な光熱粘弾性モデルおよびMoore-Gibson-Thompson半導体モデル^27、28、29でも報告されています。

回転効果も一般化熱弾性の重要な側面であり、回転は追加のコリオリ力や向心力を導入し、波の伝播特性を大きく変化させます。Khanらは、レーザー誘起熱荷重を受けた回転熱弾性固体の波動挙動の感度^{解析を行いました} 。Yadav³¹ は拡散効果を持つ回転する直交性媒体における磁気熱弾性波を調査し、Abouelregalら³² は球面空洞を含む回転応力媒体における分数的二相遅延熱弾性応答を調査しました。さらに、Khanらは非局所的な分数三相遅延モデルの下で、回転する多孔質熱弾性半空間における波動伝^播を研究 しました。Abo-Dahabらは半導体熱弾性ナノ構造への回転的影響も調査^し、回転が熱場および力場の分布を大きく変えることを示しました。これらの努力にもかかわらず、回転効果は異方性、ファイバー強化、光励起、キャリア輸送、磁場相互作用と同時に研究されることは稀です。

前述の研究ではかなりの進展が見られましたが、既存の研究の多くは光熱弾性、磁気、回転、または補強効果の選択的な組み合わせに焦点を当てています。光励起、キャリア輸送動力学、異方性、ファイバー強化、磁場相互作用、回転影響を半導体媒体内で同時に組み込む包括的な解析モデルは、依然としてほとんど利用されていません。したがって、これらの物理メカニズム間の相互結合を正確に記述し、波動伝播への影響を明確にする統一された定式化は依然として必要とされています。

これらの観察に基づき、本研究は光励起を受ける回転するファイバー強化異方半導体半空間における磁気光熱弾性波の伝播を調査するための二次元解析モデルを開発します。支配方程式は熱場、機械場、キャリア密度場、電磁場を取り入れた統一された枠組みの中で定式化され、正規モード技術と固有値アプローチを用いて解かれます。注目すべきは、Alaofiら^が最近、 異方性繊維強化シリコンの結合熱光弾性応答を固有値定式化を用いて調査したことです。しかし、磁場効果や回転寄与はモデルでは考慮されていませんでした。本研究の新規性は、磁気相互作用、回転効果、キャリア輸送力学、異方性強化、光励起を単一の磁気・光熱弾性フレームワークに同時に組み込むことで、その定式化を拡張した点にあります。したがって、提案されたモデルは結合波伝播現象をより包括的に記述し、磁強と回転が関心のある物理場に与える複合的な影響を詳細に評価することを可能にします。Aloofiら³⁵が磁場や回転効果を伴わない異方性繊維強化半導体の熱光弾性応答を考慮したのに対し、本定式化は回転、磁場、光励起、キャリア輸送、異方性強化の同時相互作用を統一枠組みの中で組み込んでいます。さらに、固有価値の実装はこのより一般的な結合システムに拡張され、従来のモデルで扱われなかった結合物理効果の解析が可能になりました。さらに、本研究の特徴的な点は固有値アプローチの実装です。多くの関連研究、例えばAlaofiら³⁵は解を得る前に支配系を還元するための消去手続きに依存していますが、本書の定式化は行列ベースの固有値フレームワークを用いて構築・解解されています。この実装は、回転、磁場、光励起、キャリア輸送、ファイバー強化を含むより一般的な結合系に適用され、これまでに発表されたモデルの範囲を拡大しています。

プロトコル

本研究は、解析モデリング、記号計算、数値シミュレーションに完全に基づいており、異方性半導体媒体における結合回転磁気・光・熱弾性相互作用を調査しています。この研究には人間参加者、動物実験、臨床データ、生物標本は含まれていません。したがって、倫理的な承認やインフォームド・コンセントは必要ありませんでした。

回転する繊維強化異方性半導体半空間における磁気光熱弾性問題の数学的定式化

本研究では、光励起と加磁場の下で2次元の回転ファイバー強化異方性半導体半空間を調査します。媒質は半無限領域 x ≥ 0 を占め、ここで x = 0 の境界は外部の光学負荷を受ける露出面を表します。座標系は、x軸が媒質に伸び、y軸が表面に沿って配置されるように選ばれ、構造の面内振る舞いを表しています。加えられた磁場と角速度ベクトルはどちらもz軸に沿って取られます。半導体媒体は均質かつ線形弾性であると仮定され、異方性はx方向に埋め込まれた整列した補強繊維を通じて導入されます。境界での光学吸収は局所的な加熱と過剰電荷キャリアを生み出し、媒質内で熱的、機械的、キャリアの結合相互作用を引き起こします。さらに、磁場は電磁的結合効果をもたらし、一方で回転運動は熱弾性波の伝播や全体的な物理応答に大きな影響を与える慣性効果をもたらします。したがって、媒質の物理状態は、小さな変形を仮定して温度場T(x, y, t)(K)、キャリア密度N(x, y, t)(^m-3)、変位成分u(x, y, t)(m)およびv(x, y, t)(m)で表されます。図1は問題の幾何学的特徴を示しており、座標系、光励起、磁場、回転効果、ファイバーの向きなどが含まれます。本定式は、小変形領域および線形熱弾性の枠組み内で動作する均一異方性繊維強化半導体媒体に適用可能です。このモデルは、繊維の向きが固定され、媒質全体で材料特性が一定であると仮定しています。したがって、非線形材料の挙動、大きな変形、材料損傷、材料特性の空間的変動は本研究では考慮されていません。したがって、提案モデルは応答が線形範囲内に収まる中程度の荷重条件下を想定しています。本研究では、光励起を表面温度と光生成キャリア密度の境界条件を用いてモデル化します。光学吸収、浸透深度、強度分布を含むレーザーと物質の詳細な相互作用過程は明示的に扱われていません。代わりに、その純効果は境界振幅θ₀およびN₀で表され、入射光場によって誘起される熱励起および搬送波励起を特徴づけます。

図1:光励起と外部磁場を受けた回転異方性繊維強化半導体半空間の回路図。 図は、本案で考慮される座標系、光励起、加磁場、光ファイバーの向き、回転効果を含む問題の物理的構成を示しています。 この図の拡大版はこちらをクリックしてご覧ください。

繊維強化異方性熱弾性半導体媒体における応力テンソルの構成関係は、一般化された形で次のように表^{すことができます。}

(1)

ここで、_ij は応力テンソル成分を表し、C σ_ijklは弾性剛性係数、e_klはひずみテンソル、Tは基準温度T ₀に対する温度増分、Nは過剰キャリア密度を表します。テンソルβ_ij、η_ijはそれぞれ熱弾性結合係数とキャリア結合係数に対応します。したがって、繊維強化の明示的影響を含む構成関係は^12,15と書けます。

(2)

この定式化では、λ と _{μ T} はラメ弾性定数であり、μ_L は繊維方向に沿った縦方向のせん断率を表します。パラメータαとβは埋め込み繊維に関連する補強効果を表します。量 _ij δ はクロネッカー・デルタ記号であり、a_i はファイバーの配向を定義する単位ベクトルの成分です。本モデルでは、補強繊維はx方向に沿って整列し、a = (1,0) となります。_ijθ β と _{η ij}N を含む項は、それぞれ熱的結合効果とキャリア結合効果を表します。現在の二次元構成では、支配応力成分は次の形^12,15に還元されます。

.   (3)

。   (4)

.   (5)

ここで、 と はそれぞれ x 方向と y 方向に沿った変位成分を表し、 A_ij は異方性繊維強化媒体の有効弾性係数を表します。熱弾性係数とキャリア結合係数は次のように定義されます

,

,

,

.

上記の関係式では、α_ijは熱膨張係数を表し、ξ_ijは半導体媒体に関連するキャリア膨張係数を表します。繊維強化異方性媒体の有効弾性係数は次のように与えられます。

、 、 、 。

これらの係数は、強化半導体材料の異方性弾性応答を特徴づけ、ファイバーの向きが結合熱弾性挙動にどのように影響するかを記述します。現在の磁気光熱弾性定式化における電磁相互作用の影響を説明するため、変形面x-yに垂直なz方向に均一磁場が加わると仮定します。したがって、磁場ベクトルは^23,25 , の形で考えられ、H₀は一定の磁場強度を表します。本式の定式化は二次元変形に限定されているため、媒質の変位場はとして用いられ、と はそれぞれx方向とy方向の変位成分を表します。

小さな変形とゆっくりと移動する電気伝導性半導体媒質を仮定すると、粒子速度と加磁場との相互作用が誘導電場を生み出します。

移動する導電媒質に関するマクスウェルの電磁関係に基づくと、誘導電場ベクトルは^19,23として表現できます。

.    (6)

ここでμ₀は磁気透磁率、は粒子速度ベクトルを表します。上記の関係に と の式を代入すると、

.   (7)

これにより、次のようになります。

.  (8)

誘導電場の時間微分を取ると、次のようになります

.(9)

半導体媒体の変形によって生成される磁気摂動ベクトルは²³˒²⁵と定義されます

.(10)

上記の式は磁気摂動場のマクスウェル発散条件を自動的に満たします。

.(11)

電流密度を求めるために、まず磁気摂動ベクトルのカールを評価します。行列式の形では、カール演算子は²³と表せます

.(12)

行列式を展開すると次のようになります。

.(13)

電流密度ベクトルはマクスウェルの電磁方程式²³から求められます

.(14)

ここで _{ε 0} は媒質の電気誘電率を表します。

半導体媒体に作用する電磁体力はローレンツ力関係式²³を用いて求められます。

.(15)

ベクトル積 は行列式で次のように評価できます。

.(16)

前の式を式(15)に代入すると、電磁体力ベクトルの成分は次のようになります

.(17)

.(18)

これらの関係は、加えられた磁場が剛性に似た電磁項と、強度に比例した修正慣性項の両方を通じて、支配方程式に追加の結合機構を寄与していることを明確に示している。その結果、磁場は熱弾性波の伝播特性や繊維強化異方性半導体媒体の全体的な動的挙動に大きな影響を与えます。電磁効果に加え、回転の影響も本公式化に組み込まれ、回転座標系から観測された際の媒体の動的応答を記述しています。繊維強化半導体媒体は、一定の角速度ベクトル³³で定数の剛体回転を経験すると仮定されます。ここではz軸周りの一定の角速度を表します。回転軸はx-y平面に垂直であるため、回転座標系で運動方程式を定式化すると、回転系が非慣性であるため追加の慣性加速度が生じます。これらの加速度は主にコリオリ加速度と遠心加速度で構成されています。コリオリ加速度は粒子速度場と関連し、³¹˒³³と表されます 

.(19)

と の式を代入すると、次のようになります。

.(20)

したがって、コリオリ加速度は次のようになります。

.(21)

遠心加速度は変位場自体に直接依存し、^31.35で表されます。

.(22)

まず、ベクトル積 を次のように評価します。

.(23)

次に、得られた結果を遠心加速度関係に代入すると

.(24)

したがって、支配方程式に現れる総回転寄与は次のように表現できます。

.(25)

したがって、 x方向および y方向の回転加速度成分は次のようになります。

, (26)

.(27)

上記の式は、回転運動が遠心加速度による変位依存の慣性効果に加え、コリオリ加速度を介して変位成分間に追加の結合をもたらすことを示しています。したがって、磁場と回転の複合作用により、回転するファイバー強化異方性半導体媒体の動的応答および波伝播特性に大きな変化が生じます。電磁相互作用と回転運動の複合的影響を取り込むために、繊維強化異方性半導体媒体の運動方程式は、電磁体力と回転座標系で生じる追加の慣性加速度の両方を含むように一般化されます。したがって、変形可能な回転連続体に対する一般運動方程式は次のように表せます³¹˒³²

.(28)

ここで、ρは質量密度を表し、F_iは電磁体力成分を表します。さらに、は回転系の角速度ベクトルを表します。エクイズ。(22)と(24)はそれぞれコリオリ加速度と遠心加速度に対応します。x方向とy方向の運動方程式は次のように書けます。

.(29)

.(30)

先に得られた電磁体力成分を上記の式に代入すると、

.(31)

.(32)

次に、繊維強化異方半導体媒体に対応する本構関係を上記の式に代入すると、変位成分、温度場、キャリア密度に関する結合運動方程式が得られます³²

.(33)

.(34)

最後に、式(15)で示された磁気摂動場の式を用いて、これを上記の方程式に代入すると、支配運動方程式は最終的な結合形で次のように書けます。

.(35)

.(36)

これらの方程式は、回転効果と磁場効果が媒質の熱弾性応答に与える連合的影響を示しています。回転項はコリオリ寄与と遠心寄与の両方を考慮し、磁場は追加の電磁結合をもたらし系の動的挙動を変化させます。したがって、支配方程式は回転する磁気光熱弾性繊維強化半導体における波動伝播および多物理的相互作用の統一的な枠組みを確立します。光励起下では、半導体媒体の熱挙動は熱伝導、キャリア輸送、機械的変形の相互作用によって大きく影響を受け、強く結合した熱光弾性過程が形成されます。古典的な熱伝導モデルとは異なり、半導体材料の温度分布は熱拡散だけでなくキャリア再結合や熱弾性結合の影響も受けます。したがって、異方性繊維強化半導体媒体の一般化熱伝導方程式は次のように表現^{できます。}

.(37)

上記の式は、回転する磁気光熱弾性半導体媒体内の熱場が異方性熱伝導、キャリア再結合過程、熱弾性相互作用の複合的な影響によって支配されることを明確に示しています。という用語は光励起下でのキャリア再結合によって生じる熱エネルギーを表し、との結合項は時間依存的な機械的変形が媒体の熱応答に与える影響を示します。その結果、温度場はキャリア密度および弾性場と強く結合し、繊維強化半導体材料における熱弾性波の伝搬特性に重要な役割を果たします。本式では、半導体媒体内のキャリア濃度N(x, y, t)の進化は、キャリア拡散、再結合過程、光励起による熱活性化の複合効果によって制御されます。したがって、非平衡搬送動力学を記述するキャリア輸送方程式は次のように書けます ^4,23

.(38)

ここでD _Eはキャリア拡散係数を表し、はx-y平面上の二次元ラプラシアン作用素を表します。項はキャリア寿命τに関連するキャリア再結合効果に対応します。さらに、κはによって定義される熱キャリア結合パラメータであり、ここでN₀は平衡キャリア濃度を表します。結合項κTは、半導体媒体内の余剰キャリア生成に対する温度場の影響を表します。上記の式は、キャリア濃度が熱活性化されたキャリア生成機構を通じて熱場と強く結合していることを示しています。その結果、キャリアの力学は熱拡散効果と再結合効果の両方に大きく依存し、回転する繊維強化異方性半導体媒体の結合光熱弾性応答に大きな影響を与えます。結合磁気光熱弾性半導体系の支配方程式と数学的定式化は、現在完全に確立されています。シリコン媒質の物理的および材料パラメータは表2にまとめられており、その数値、単位、対応する参照も含まれています。これらのパラメータは、数値計算や無次元化の手法で利用されます。

回転磁気光熱弾性繊維強化半導体モデルの無次元定式化

支配方程式を簡略化し、結合した回転磁気・光熱弾性系のコンパクトな数学的表現を得るために、物理変数を無次元化するための適切な特性尺度が導入されます。この無次元化手順により、物質パラメータの数を減らし、結合方程式の解析的および数値的処理を容易にします。選択された特性量は、半導体媒体の熱弾性、電磁的、回転性、キャリア輸送特性と整合的に選ばれます。^16,21 それに従い、以下の無次元変数が導入されます。

、 、 、 、 1、 、 、 、 、 、 。

ここで、C_Tは特徴的な弾性波速度を表し、t*は結合熱弾性過程に関連する特徴的な熱緩和時間を表します。さらに、パラメータは回転運動が媒体の動的挙動に与える影響を特徴づける無次元回転パラメータを定義します。上記の無次元量を、以前導出された支配方程式に代入すると、結合系は正規化された形に変換されます。この変換により方程式の数学的構造が大幅に簡素化され、磁場、光励起、回転、異方性、キャリア-輸送相互作用の複合効果を調査するための適切な枠組みを提供します。簡便化のため、無次元量に関連する素記法は後の解析で省略されます。したがって、結合回転磁気光熱弾性系の支配方程式は、次の無次元形^16,20で表すことができます。

, (39)

, (40)

, (41)

.(42)

回転異方性繊維強化半導体媒体の対応する無次元応力成分は以下のように得られます。

, (43)

, (44)

.(45)

無次元係数 _i (i = 1,2,...,18) は、結合された半導体媒質の物理的、熱的、電磁的、キャリア的、回転的パラメータの組み合わせを表します。これらの係数は異方性、ファイバー強化、磁場、熱弾性結合、キャリア輸送、回転運動が系全体の挙動に与える影響を特徴づけます。したがって、得られた無次元の支配方程式は、回転する磁気光熱弾性繊維強化半導体媒体における結合波伝播現象や多物理相互作用を解析するためのコンパクトかつ効率的な数学モデルを提供します。無次元パラメータ a_i、γ_i、δ_i は、繊維強化異方性半導体媒体の結合回転磁気・光熱弾性挙動を支配する物理的および材料的性質のコンパクトな組み合わせを表すために導入されました。各係数は結合系内の特定の相互作用メカニズムを反映し、基盤となる物理プロセスの相対的な影響を示す洞察を提供します。明確化のために、無次元パラメータとその定義は表1にまとめられています。

表1:本定式化で用いられている無次元パラメータの定義と物理的解釈。 この表は、支配方程式に現れる無次元パラメータとその物理的意味、そして結合した熱弾性相互作用、電磁的相互作用、キャリア密度相互作用、回転相互作用を記述する役割をまとめています。 この表をダウンロードするには、こちらをクリックしてください。

正規モード技術を用いた解析解法

結合回転磁気光熱弾性系の解析解を導出するために、ノーマルモード技術が用いられます。この方法は、結合された偏微分方程式を常微分方程式の還元系に変換する効果から、一般化熱弾性理論や半導体理論で広く用いられています。この手法は、異方性半導体媒体における波動伝搬、減衰、多物理的相互作用の解析に特に有用です。正規モード解析の後、すべての物理場量は時間と横方向の空間座標に対して調和的に変化すると仮定します。したがって、温度場、キャリア密度、変位成分、応力量は指数関数形式^24,27で表されます。

.(46)

ここでωは 物理場の時間的変動を支配する複素周波数パラメータを表します。 ω の実部は波振幅の時間的減衰(または成長)に関連しており、虚部は伝搬モードの振動振る舞いを表します。これらの解釈は本研究で採用された従来の正規モード解析と整合しており、 a は y方向に沿った空間変動に関連する波数を表します。 と の量は、空間座標 x のみに依存する場の振幅に対応します。上記の正規モード表現を、既に得られた無次元の支配方程式に代入し、得られる式を簡略化することで、元の結合偏微分系は空間座標xに関する常微分方程式の集合に変換 されます。したがって、変換領域内の支配方程式は次の形をとります。

, (47)

, (48)

, (49)

.(50)

さらに、対応する変換された応力成分は次のように得られます。

, (51)

, (52)

.(53)

ここで、 は空間座標に関する微分演算子を表します。得られた変換系は、結合した磁気・光・熱弾性問題の特性方程式を構築し、完全な解析解を導出するための数学的基盤を形成します。変換後の方程式に現れる係数は次のように定義されます。

、 、 、、 、、 、 、 、 、 、、 。

これらの係数には、異方性弾性、磁場相互作用、熱結合、キャリア輸送、回転効果の複合的な寄与が含まれています。したがって、変換されたシステムは、特徴的な根の取得や回転するファイバー強化半導体媒体内での結合波伝搬挙動の調査に適したコンパクトな表現を提供します。

補足ファイル1:行列形式を用いた解析解。このファイルには、詳細な行列の定式化、固有値解法の手順、特性方程式の導出、そして結合磁気・光熱弾性モデルの一般解を得るための中間解析ステップが含まれています。このファイルをダウンロードするには、こちらをクリックしてください。

詳細な行列の定式化、固有値解法の手順、特性方程式の導出は補足ファイル1に記載されています。

境界条件と未知定数の決定

解析的定式化を完成させるために、得られた一般解を表面 x = 0 に課された境界条件に代入しました。この置換により、未知の振幅定数を含む結合代数系が生成。温度、キャリア密度、機械的変位、応力制約に関連する各境界要件は許容される固有モードを用いて表現され、係数を関連付ける線形方程式の集合を導き出しました。便宜上、得られた代数系はコンパクト行列形式に書き換えられ、BC = Dとして表されます。ここでBは境界面で評価された固有ベクトル成分から構成される係数行列、は未知定数のベクトル、Dは熱荷重パラメータθ₀を含む境界条件から生成されるベクトルに対応します。 キャリア励起項_N0と、規定された変位条件。これらの定数を評価した後、完全な解析解を得るために場変数の一般式に代入されました。これらの式はその後、回転するファイバー強化異方性半導体媒体内の熱弾性場、キャリア密度場、変位場の数値計算および図解に用いられました。採用された境界条件は、同時に熱励起とキャリア励起を受ける光学照明された半導体表面を表しています。所定温度条件は入射光場によって生じる熱負荷をモデル化し、キャリア密度条件は光照明によって生成される過剰キャリアを考慮します。さらに、横変位制約は -方向における表面の機械的閉じ込めを表し、ゼロのせん断応力条件は接線的にトラクションフリー境界に対応します。したがって、選択された混合熱、電子、機械的境界条件は、半導体表面における結合された光熱弾性相互作用の物理的に一貫した表現を提供し、溶液の未知定数を決定するために必要な制約を提供します。課される境界条件は以下の通りです:

温度制約:

.(78)

この条件は、周期的な光学加熱によって誘導される調和的に変化する表面温度を表しています。これは媒質内の結合された熱弾性およびキャリア輸送過程を駆動する一次熱励起として機能します。振幅θ0は加えられた熱荷重の強度を示します。

キャリア密度制約:

.(79)

この境界条件は、光照明によって生じる光生成キャリア密度を記述します。これは光子吸収による電子励起と、入射光場と整合した調和変調を反映します。

変位制約:

.(80)

この条件は、境界が横方向に機械的に制約されていることを示します。したがって、表面の -方向に沿った変位は起こりません。

せん断応力制約:

.(81)

この条件はせん断応力に関してトラクションフリー境界に対応します。これは接線方向における機械的に自由な境界と整合的に表面に接力が作用しないことを保証します。 x = 0 での境界条件に加え、無限遠での物理的要件は次のように課されました。 半無限領域内で有界物理解を保証すること。本研究で採用された解析および計算手法の明確な概要を提供するため、解法手法の主なステップは 図2にまとめられています。

図2:本研究で採用された解析解手順のフローチャート。支配方程式の定式化、正規モード解析、固有値解、境界条件の適用、物理場変数の評価を含む。 図は解析解を得るための主要な計算手順をまとめ、その後結合した磁気・光・熱弾性応答を数値的に評価します。 この図の拡大版はこちらをクリックしてご覧ください。

結果

本研究で提示された数値的結果を得るために、シリコン(Si)が半導体および熱弾性工学の応用に広く応用されていることから代表的な半導体材料として選ばれました。計算に用いられる物理的、熱的、弾性的、電磁的、キャリア関連の材料特性は表2にまとめられています 。 これらの物質定数は無次元支配方程式に置き換えられ、後のセクションで議論される数値的結果を生成しました。さらに、数値計算は無次元空間領域 上で空間ステップサイズ を用いて行われました。この選択により、十分な数値分解能が得られ、考慮対象となる物理場の完全な空間的挙動が計算領域内で正確に捉えられることが判明し...

ディスカッション

固有値の許容性、安定性、および有界性解析

得られた固有値の許容性は、すべての場変数が半無限領域 x ≥ 0 内に有界であるという物理的要件によって支配されます。一般的な固有解は、の形の指数モードで表されます。したがって、解の漸近的挙動は固有値の実部に直接依存します。無限遠における有界性条件を満たすため、

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
計算ソフトウェア（記号および数値解析）	Wolfram Research	Wolfram Mathematica（バージョン12.0）が使用されました	固有値の計算、解析解の実装、数値評価に使用
材料パラメータデータセット（シリコン半導体特性）	様々な文献ソース	N/A	計算に使用される物理定数（弾性、熱、キャリア関連）（表1）
パーソナルコンピュータ/ワークステーション	HP	N/A	数値シミュレーションに十分なメモリを持つWindowsOS上で動作する標準的なパーソナルコンピュータで計算が行われました
方程式エディタ	Microsoft  WordとMathType	N/A	原稿中の数学表現の書式設定と表示に使用
参考文献管理ソフトウェア	Elsevier	N/A	参考文献の管理と引用の書式設定（バンクーバースタイル）に使用