#數值計算

**ChatGPT:**

這裡的「電感片」，沿用我們模型中的**可調表面動能電感薄膜**。它的電感來自電子集體運動的慣性，以 $L_s$ 表示；控制線路另外具有的串聯電感則是 $L_g$。

以下依照「電子運動 → 表面響應 → 集體模態 → 交互能 → 控制功」推導。幾何先取**一片可調薄膜，面對距離為 $a$ 的理想導體**。

1. **電子慣性如何形成表面電感**設載子面密度為 $n$、有效質量為 $m^*$、散射率為 $\gamma$，電子電荷為 $-e$。面內電子速度滿足$$m^*\dot{\mathbf v} +m^*\gamma\mathbf v =-e\mathbf E_\parallel.$$在固定 $n$ 的線性響應下，表面電流密度為$$\mathbf J_s=-en\mathbf v.$$代入得到$$\boxed{ \mathbf E_\parallel =L_s\dot{\mathbf J}_s+R_s\mathbf J_s },$$其中$$\boxed{L_s=\frac{m^*}{ne^2}}, \qquad \boxed{R_s=\gamma L_s}.$$這就是電感片的基本方程。電子需要時間才能被加速，因此電流對電場具有慣性響應。它儲存的電子動能密度為$$u_{\mathrm{kin}} =\frac12nm^*|\mathbf v|^2 =\boxed{\frac12L_s|\mathbf J_s|^2}.$$對長度 $\ell$、寬度 $w$ 的均勻薄膜條，$$I=wJ_s,$$所以$$E_{\mathrm{kin}} =\frac12L_{\mathrm{strip}}I^2, \qquad \boxed{L_{\mathrm{strip}}=L_s\frac{\ell}{w}}.$$$L_s$ 的單位是亨利，通常稱為「每平方的表面電感」。二維電子氣的這種動能電感及其傳輸線響應，已有實驗量測。[CaltechAUTHORS](https://authors.library.caltech.edu/records/15mx4-tv697?utm_source=chatgpt.com)
2. **表面電感如何控制電磁反射**採時間慣例 $e^{-i\omega t}$，上式給出$$Z_s(\omega)=R_s-i\omega L_s,$$$$\boxed{ \sigma_s(\omega) =\frac1{R_s-i\omega L_s} =\frac{ne^2}{m^*(\gamma-i\omega)} }.$$計算平衡交互能時，使用虛頻 $\omega=i\xi$：$$\boxed{ \sigma_s(i\xi)=\frac1{R_s+\xi L_s} }.$$薄膜的 Maxwell 邊界條件是$$\mathbf E_{\parallel,2} =\mathbf E_{\parallel,1},$$$$\hat{\mathbf n}\times (\mathbf H_2-\mathbf H_1) =\mathbf J_s =\sigma_s\mathbf E_\parallel.$$設面內波數為 $k$，並定義$$\kappa_\perp=\sqrt{k^2+\xi^2/c^2}.$$解邊界條件得到真空兩側薄膜的反射係數$$\boxed{ r_{\mathrm{TM}} =\frac{\kappa_\perp\sigma_s} {2\varepsilon_0\xi+\kappa_\perp\sigma_s} },$$$$\boxed{ r_{\mathrm{TE}} =-\frac{\mu_0\xi\sigma_s} {2\kappa_\perp+\mu_0\xi\sigma_s} }.$$這是以面導電率描述薄膜反射的形式；上述公式採 SI 制。[arxiv.org](https://arxiv.org/pdf/1007.1231?utm_source=chatgpt.com)因此，在其他參數固定時：$$n\uparrow \ \Rightarrow\ L_s\downarrow \ \Rightarrow\ \sigma_s(i\xi)\uparrow.$$材料對不同頻率與波數的反射能力，便隨之改變。
3. **電感片如何形成集體諧振子**接著取無耗散、非延遲近似：$$R_s\to0, \qquad \kappa_\perp\simeq k.$$對面電荷的一個波數模態，鏡像法給出的電勢為$$\phi_k =\frac{1-e^{-2ka}}{2\varepsilon_0k}\rho_k.$$定義$$f_k(a)=1-e^{-2ka}.$$這相當於一個模態電容係數$$\boxed{ C_{k,\mathrm{eff}} =\frac{2\varepsilon_0k}{f_k(a)} }.$$另一方面，電荷連續性使電流與電荷變化滿足$$\dot\rho_k+ikJ_k=0.$$因此電子動能提供的模態慣性係數為$$\boxed{\mathcal M_k=\frac{L_s}{k^2}}.$$對一致正規化的實電荷模態 $Q_k$，其哈密頓量就是$$\boxed{ H_k= \frac{P_k^2}{2\mathcal M_k} +\frac{Q_k^2}{2C_{k,\mathrm{eff}}} }.$$所以正常模態頻率為$$\boxed{ \omega_k^2 =\frac1{\mathcal M_kC_{k,\mathrm{eff}}} =\frac{k(1-e^{-2ka})} {2\varepsilon_0L_s} }.$$代入電子密度：$$\boxed{ \omega_k^2(a,n) =\frac{ne^2k}{2\varepsilon_0m^*} (1-e^{-2ka}) }.$$**機制因此很明確：電子慣性提供電感，電荷產生的電場提供恢復作用，另一片導體透過鏡像電荷改變恢復係數。**在長波極限 $ka\ll1$，$$C_{k,\mathrm{eff}}\simeq\frac{\varepsilon_0}{a},$$正好回到平行板電容的單位面積形式。
4. **模態如何產生交互能與壓力**先定義孤立電感片的模態頻率$$\Omega_k^2=\frac{k}{2\varepsilon_0L_s}.$$則$$\omega_k(a)=\Omega_k\sqrt{1-e^{-2ka}}.$$同一結果也能從反射公式得到：$$r_{\mathrm{TM}}(i\xi,k) =\frac{\Omega_k^2}{\xi^2+\Omega_k^2},$$$$1-r_{\mathrm{TM}}e^{-2ka} =\frac{\xi^2+\omega_k^2(a)} {\xi^2+\Omega_k^2}.$$因此，短距離 TM 交互能密度為$$u_C= \frac{\hbar}{4\pi^2} \int_0^\infty k\,dk \int_0^\infty d\xi\, \ln\frac{\xi^2+\omega_k^2(a)} {\xi^2+\Omega_k^2}.$$利用$$\int_0^\infty \ln\frac{\xi^2+x^2}{\xi^2+y^2}\,d\xi =\pi(x-y),$$得到$$\boxed{ u_C= \frac{\hbar}{4\pi} \int_0^\infty k\,dk\, [\omega_k(a)-\Omega_k] }.$$令 $u=ka$，即可整理為$$\boxed{ u_C(a,L_s) =-\frac{\mathcal C\hbar} {\sqrt{2\varepsilon_0L_s}\,a^{5/2}} },$$其中$$\mathcal C= \frac1{4\pi}\int_0^\infty u^{3/2} \left[1-\sqrt{1-e^{-2u}}\right]du>0.$$用載子密度表示：$$\boxed{ u_C(a,n)=-\frac{\eta\sqrt n}{a^{5/2}} }, \qquad \eta=\frac{\mathcal C\hbar e} {\sqrt{2\varepsilon_0m^*}}.$$吸引壓力的大小為$$\boxed{ p_C=\frac{\partial u_C}{\partial a} =\frac{5\eta\sqrt n}{2a^{7/2}} }.$$板受到的力為$$F_a=-Sp_C.$$這段推導限於短距離近似。完整延遲模型需共同計入光子與表面模態；兩者可能有顯著補償。[arxiv.org](https://arxiv.org/pdf/hep-th/0511269?utm_source=chatgpt.com)
5. **切換電感片響應，需要交換多少功**固定板距時，交互能對表面電感的響應是$$\boxed{ \frac{\partial u_C}{\partial L_s} =-\frac{u_C}{2L_s}>0 }.$$所以增加 $L_s$、弱化載子響應時，這部分交互能會升高。以電子密度表示，每增加一個載子造成的交互能變化為$$\boxed{ \mu_C=\frac{\partial u_C}{\partial n} =-\frac{\eta}{2\sqrt n\,a^{5/2}} }.$$但 $u_C$ 只是交互能。孤立薄膜及控制端隨密度改變的能量，必須放在 $\Phi(n)$：$$\boxed{ E(a,n)=S[\Phi(n)+u_C(a,n)] }.$$因此完整微分為$$\boxed{ dE=Sp_C\,da +S[\Phi'(n)+\mu_C]\,dn }.$$這同時給出機械取功與控制功：$$\boxed{P_{\mathrm{gen}}=-Sp_C\dot a},$$$$\boxed{ P_{\mathrm{ctrl}} =S[\Phi'(n)+\mu_C]\dot n }.$$若控制電荷為 $Q_g=eS(n-n_0)$，則$$eV_g=\Phi'(n)+\mu_C, \qquad I_g=eS\dot n.$$所以$$P_{\mathrm{ctrl}}=V_gI_g.$$**電感片的切換功，不能只由 $$L_s=m^*/(ne^2)$$ 判定；還需要材料及載子儲存端的能量函數 $$\Phi$$。**
6. **它在裝置中的作用**在目前模型內，$$n_H>n_L \quad\Longrightarrow\quad L_{s,H}<L_{s,L},$$且同一板距下$$\boxed{ \frac{p_L}{p_H} =\sqrt{\frac{n_L}{n_H}} }.$$因此電感片負責切換強、弱交互作用；機械機組負責把力與位移轉成儲能器能量。當一輪結束時恢復 $n_H$，但保留較小板距 $a_{n+1}<a_n$，在先前的理想控制假設下：$$\boxed{ W_{\mathrm{net},n} =S\eta\sqrt{n_H} \left(a_{n+1}^{-5/2}-a_n^{-5/2}\right)-D_n }.$$電感片因此控制的是**集體響應與功的交換路徑**；淨回收量仍由實際端點交互能下降與耗散決定。

以上的 $n$ 是參與面內響應的有效載子密度。實際薄膜還可能有背景介電與帶間響應，需加入完整導電率。切換採準靜態描述；前面的 Drude 方程是每個固定密度下的線性響應，完整的載子注入過程則須連同控制端與儲存端一起建模。

* * *

**ChatGPT:**

需要。**現在應把幾何、材料模型、控制能量與操作端點固定，才能得到唯一的數值結果。**

下面先建立一個可直接計算的最簡版本：**電感薄膜與理想反射板同時構成控制電容。**這會使先前的控制能 $\Phi(n)$ 具體化為依賴板距的 $\Phi(a,n)$。

1. **固定幾何與近似**

* 兩片平行板，面積 $S$，間距 $a$。
* 一片是理想導體；另一片是零厚度、拋物線能帶的二維載子薄膜。
* 薄膜具有固定正電背景 $+en_0$，電子面密度為 $n$。
* 採零溫、無耗散 Drude 響應、非延遲與大面積近似。
* 機組與控制器共用可雙向交換能量的理想儲能器。

此處的控制電荷是$$Q_g=eS(n-n_0),$$幾何電容直接由間距決定：$$\boxed{C_g(a)=\frac{\varepsilon_0S}{a}}.$$若採獨立控制電極，則須另定其位置、介電層及電磁反射作用，才能建立另一套數值模型。
2. **必要輸入與單位**

| 輸入 | 定義 | SI 單位 |
| --- | --- | --- |
| $S$ | 有效板面積 | $\mathrm{m^2}$ |
| $m_r$ | $m^*/m_e$ | 無因次 |
| $g_d$ | 自旋與谷簡併度 | 無因次 |
| $n_0$ | 中性背景面密度 | $\mathrm{m^{-2}}$ |
| $n_H,n_L$ | 強、弱響應載子面密度 | $\mathrm{m^{-2}}$ |
| $a_i,b,a_f$ | 初始、最小、最終板距 | $\mathrm m$ |
| $T_1,\ldots,T_4$ | 四段操作時間 | $\mathrm s$ |
| $R_g$ | 控制通路串聯電阻 | $\Omega$ |
| $\zeta$ | 機構黏滯阻力係數 | $\mathrm{kg/s}$ |

端點條件為$$0<b\le a_f<a_i, \qquad n_H>n_L>0.$$常數採用$$e=1.602176634\times10^{-19}\ \mathrm C,$$$$\hbar\simeq1.054571817\times10^{-34}\ \mathrm{J\,s},$$$$m_e=9.1093837139\times10^{-31}\ \mathrm{kg},$$$$\varepsilon_0=8.8541878188\times10^{-12}\ \mathrm{F/m}.$$這些數值依據 CODATA 2022。[physics.nist.gov](https://physics.nist.gov/cuu/Constants/Table/allascii.txt?utm_source=chatgpt.com)
3. **將先前的比例常數全部定出**表面電感：$$\boxed{ L_s(n)=\frac{m_rm_e}{ne^2} }.$$若 $n$ 以 $\mathrm{m^{-2}}$ 輸入：$$\boxed{ L_s[\mathrm H] =3.54869119\times10^7\frac{m_r}{n} }.$$模態頻率：$$\boxed{ \omega_k^2(a,n) =\frac{ne^2k}{2\varepsilon_0m_rm_e} (1-e^{-2ka}) },$$其中 $k$ 的單位為 $\mathrm{m^{-1}}$，$\omega_k$ 為 $\mathrm{rad/s}$。先前的模態積分常數，本次積分得到$$\boxed{\mathcal C=0.0098813543181}.$$因此定義$$\eta(m_r) =\frac{\mathcal C\hbar e} {\sqrt{2\varepsilon_0m_rm_e}},$$即$$\boxed{ \eta(m_r) \simeq\frac{4.1568961\times10^{-35}}{\sqrt{m_r}} \ \mathrm{J\,m^{3/2}} }.$$集體模態交互能密度便可直接計算：$$\boxed{ u_C(a,n)=-\eta(m_r)\frac{\sqrt n}{a^{5/2}} }.$$這個冪次式使用短距離表面模態近似；完整延遲模型需共同計入光子與表面模態。[arxiv.org](https://arxiv.org/pdf/hep-th/0511269?utm_source=chatgpt.com)
4. **明定完整能量函數**對簡併度 $g_d$ 的理想二維拋物線能帶，定義$$\boxed{ \alpha=\frac{2\pi\hbar^2}{g_dm_rm_e} }.$$以固定儲存端化學勢 $\mu_0=\alpha n_0$ 為參考，載子能量密度為$$u_{\mathrm{band}} =\frac{\alpha}{2}(n-n_0)^2.$$平均充電的靜電能密度為$$\boxed{ u_{\mathrm{dc}} =\frac{e^2a}{2\varepsilon_0}(n-n_0)^2 }.$$所以這一版的完整計算式是$$\boxed{ \begin{aligned} E(a,n)=S\Bigg[ &\frac{\alpha}{2}(n-n_0)^2\\ &+\frac{e^2a}{2\varepsilon_0}(n-n_0)^2\\ &-\frac{\eta\sqrt n}{a^{5/2}} \Bigg]. \end{aligned} }$$**這是本次具體化最重要的地方：反射板兼作控制電極時，普通靜電儲能也依賴板距，必須和集體交互能一起核算。**
5. **由同一能量式計算力、電壓與切換功**定義吸引壓力的大小$$p_{\mathrm{tot}} =\frac1S\frac{\partial E}{\partial a}.$$則$$\boxed{ p_{\mathrm{tot}} = \underbrace{\frac{e^2(n-n_0)^2}{2\varepsilon_0}}_{p_{\mathrm{dc}}} + \underbrace{\frac{5\eta\sqrt n}{2a^{7/2}}}_{p_C} }.$$板受到的力為$$F_a=-Sp_{\mathrm{tot}}.$$與 $Q_g$ 共軛的控制電壓為$$\boxed{ eV_g(a,n) = \left(\alpha+\frac{e^2a}{\varepsilon_0}\right)(n-n_0) -\frac{\eta}{2\sqrt n\,a^{5/2}} }.$$固定板距下，密度切換功直接是$$\boxed{ W_{\mathrm{sw}}(a;n_1\to n_2) =E(a,n_2)-E(a,n_1) }.$$固定密度下，板距改變所交給機組的功是$$\boxed{ W_m(a_1\to a_2;n) =E(a_1,n)-E(a_2,n) }.$$正值表示機組取功，負值表示機組供功。
6. **完整一輪與損耗**路徑固定為$$(a_i,n_H)\to(b,n_H)\to(b,n_L) \to(a_f,n_L)\to(a_f,n_H).$$理想淨回收為$$W_0=E(a_i,n_H)-E(a_f,n_H),$$展開後：$$\boxed{ \begin{aligned} W_0={}& \underbrace{ \frac{Se^2(n_H-n_0)^2}{2\varepsilon_0}(a_i-a_f) }_{W_{\mathrm{dc}}}\\ &+ \underbrace{ S\eta\sqrt{n_H} (a_f^{-5/2}-a_i^{-5/2}) }_{W_C}. \end{aligned} }$$若每段仍採用$$s(u)=3u^2-2u^3,$$則$$\Delta Q=eS(n_H-n_L),$$$$D_g=\frac65R_g\Delta Q^2 \left(\frac1{T_2}+\frac1{T_4}\right),$$$$D_m=\frac65\zeta \left[ \frac{(a_i-b)^2}{T_1} +\frac{(a_f-b)^2}{T_3} \right].$$最後：$$\boxed{ W_{\mathrm{net}}=W_{\mathrm{dc}}+W_C-D_g-D_m },$$$$\boxed{ P_{\mathrm{avg}} =\frac{W_{\mathrm{net}}}{T_1+T_2+T_3+T_4} }.$$
7. **一組假設參數的計算示例**取$$S=(100\ \mu\mathrm m)^2=10^{-8}\ \mathrm{m^2},$$$$m_r=0.1,\qquad g_d=2,$$$$n_H=10^{16}\ \mathrm{m^{-2}}, \qquad n_L=n_0=2.5\times10^{15}\ \mathrm{m^{-2}},$$$$a_i=100\ \mathrm{nm},\quad b=60\ \mathrm{nm},\quad a_f=80\ \mathrm{nm},$$$$T_j=0.01\ \mathrm s,\quad R_g=10\ \Omega,\quad \zeta=10^{-8}\ \mathrm{kg/s}.$$代入得到：

| 計算量 | 結果 |
| --- | --: |
| $L_s(n_H)$ | $0.355\ \mathrm{nH}$ |
| $L_s(n_L)$ | $1.42\ \mathrm{nH}$ |
| 初始 $p_C$ | $0.104\ \mathrm{Pa}$ |
| 強響應 $p_{\mathrm{dc}}$ | $8.15\times10^4\ \mathrm{Pa}$ |
| 初始控制電壓 | $13.59\ \mathrm V$ |
| 集體交互能貢獻 $W_C$ | $3.10\times10^{-17}\ \mathrm J$ |
| 靜電儲能貢獻 $W_{\mathrm{dc}}$ | $1.63\times10^{-11}\ \mathrm J$ |
| $D_g+D_m$ | $3.49\times10^{-19}\ \mathrm J$ |

這組假設參數顯示：**此兩電極版本的總輸出主要來自原有靜電儲能；集體模態的貢獻須用 $W_C$ 分項辨認。**

以上已是一套可直接代入數值的理想計算模型。若要保留原先「控制端與板距無關」的版本，下一步需要明定獨立控制電極的幾何，再重新求它對靜電能與模態反射的影響；不能僅指定一個控制電容值。

* * *

**ChatGPT:**

下一步可以建立一個**背閘極、電感片、工作反射板分離的三電極模型**。

這個模型能讓控制電容不隨工作板距改變；但背閘極仍會影響電磁模態，因此必須重新計算交互能。

1. **明定三電極幾何**

| 元件 | 位置 | 功能 |
| --- | --- | --- |
| 固定理想導體背閘極 | $z=-h$ | 控制載子密度 |
| 固定二維電感薄膜 | $z=0$ | 提供可調表面響應 |
| 可移動理想反射板 | $z=a$ | 連接機械取功端 |

$h$ 固定，只有 $a$ 隨機組運動。以下先取兩個間隙皆為真空、零溫、局域無耗散 Drude 薄膜，並忽略邊緣與接線寄生效應。**工作間隙的平均靜電場需要另外補償。**令工作板與薄膜的靜電勢相同：$$\phi_{\mathrm w}=\phi_{\mathrm f}.$$此時$$E_{\mathrm{dc,work}}=0, \qquad Q_{\mathrm w}=0.$$載子的平均補償電荷由背閘極承擔：$$Q_g=eS(n-n_0).$$背閘極與薄膜之間的靜電能為$$\boxed{ U_{\mathrm{dc},g} =\frac{Se^2h}{2\varepsilon_0}(n-n_0)^2 }.$$它不隨 $a$ 改變。相應的靜電勢差是$$\boxed{ V_{\mathrm{dc},g} =\frac{eh(n-n_0)}{\varepsilon_0} }.$$若工作間隙仍有殘餘勢差 $\Delta V$，則會留下$$p_{\mathrm{dc,work}} =\frac{\varepsilon_0(\Delta V)^2}{2a^2}.$$所以獨立設置閘極後，仍須維持上述補償條件。
2. **把背閘極納入反射公式**電感片仍使用$$L_s=\frac{m^*}{ne^2}, \qquad \sigma_s(i\xi)=\frac1{\xi L_s}.$$定義$$q=\sqrt{k^2+\xi^2/c^2}, \qquad x=e^{-2qh},$$以及兩個無因次響應量$$\beta_{\mathrm{TM}} =\frac{q\sigma_s}{\varepsilon_0\xi}, \qquad \beta_{\mathrm{TE}} =\frac{\mu_0\xi\sigma_s}{q}.$$對薄膜與背閘極之間的多重反射求和後，得到$$\boxed{ \mathcal R_p = \frac{2x+\beta_p(1-x)} {2+\beta_p(1-x)} }, \qquad p=\mathrm{TE,TM}.$$這裡已將工作反射板的反射符號乘入，所以 $\mathcal R_p$ 是能量公式中使用的非負反射乘積。多層薄膜反射可沿用這類遞迴計算。[arxiv.org](https://arxiv.org/pdf/1007.1231?utm_source=chatgpt.com)這個式子有三個檢查極限：$$h\to\infty: \quad \mathcal R_p\to\frac{\beta_p}{2+\beta_p},$$$$n\to0: \quad \mathcal R_p\to e^{-2qh},$$$$h\to0: \quad \mathcal R_p\to1.$$最後一個極限說明：**背閘極太靠近電感片時，整組結構會接近固定理想反射面，載子調控對總反射的影響會變小。**
3. **改用三電極的交互能**工作板與「電感片＋背閘極」之間的交互能密度為$$\boxed{ \begin{aligned} u(a,h,n) =\frac{\hbar}{4\pi^2} \int_0^\infty d\xi \int_0^\infty k\,dk \sum_p \ln\!\left[ 1-\mathcal R_p e^{-2qa} \right]. \end{aligned} }$$這次保留傳播延遲，數值計算使用此積分，取代先前單片薄膜的 $a^{-5/2}$ 近似。工作板所受吸引壓力為$$\boxed{ \begin{aligned} p(a,h,n) =\frac{\hbar}{2\pi^2} \int_0^\infty d\xi \int_0^\infty k\,dk\,q \sum_p \frac{\mathcal R_p e^{-2qa}} {1-\mathcal R_p e^{-2qa}}. \end{aligned} }$$即使電感片的載子響應消失，背閘極與工作板仍有交互作用：$$\boxed{ u(a,h,0) =-\frac{\pi^2\hbar c}{720(a+h)^3} }.$$這是固定金屬邊界的背景貢獻，需與薄膜可調部分分別列出。
4. **控制能與控制電壓也一起補齊**定義孤立「薄膜＋背閘極」的交互能密度$$\begin{aligned} u_{\mathrm{sg}}(h,n) =\frac{\hbar}{4\pi^2} \int_0^\infty d\xi \int_0^\infty k\,dk \sum_p \ln\!\left[ 1-\rho_p e^{-2qh} \right], \end{aligned}$$其中$$\rho_p=\frac{\beta_p}{2+\beta_p}.$$沿用理想拋物線能帶，$$\alpha=\frac{2\pi\hbar^2}{g_dm^*}.$$將固定部分統一寫成$$\boxed{ \Phi_0(n,h) = \frac12\left( \alpha+\frac{e^2h}{\varepsilon_0} \right)(n-n_0)^2 +u_{\mathrm{sg}}(h,n) }.$$完整能量便是$$\boxed{ E(a,n,h)=S[\Phi_0(n,h)+u(a,h,n)] }.$$以載子儲存端為電化學參考，背閘極控制電壓滿足$$\boxed{ eV_g =\frac{\partial\Phi_0}{\partial n} +\frac{\partial u}{\partial n} }.$$工作板的補償電壓則為$$\boxed{ V_{\mathrm w} =V_g-\frac{eh(n-n_0)}{\varepsilon_0} }.$$在理想補償下 $Q_{\mathrm w}=0$，故工作板補償端的平均電功率為$$V_{\mathrm w}\dot Q_{\mathrm w}=0.$$背閘極端仍交換完整控制功：$$P_{\mathrm{ctrl}} =V_g\,eS\dot n.$$
5. **一輪回收功現在來自哪裡**沿用強響應壓縮、弱響應回推，最後恢復 $n_H$，並保持 $h$ 固定：$$a_i\longrightarrow a_f, \qquad a_f<a_i.$$所有只依賴 $n,h$ 的能量在端點恢復，因此$$\boxed{ W_0 =S[u(a_i,h,n_H)-u(a_f,h,n_H)] }.$$扣除實際損耗：$$\boxed{W_{\mathrm{net}}=W_0-D}.$$**這個版本中，工作間隙的平均靜電力已補償，背閘極的充電儲能也在端點恢復；理想淨功對應的是工作板靠近所降低的場—材料交互能。**
6. **用前例參數計算背閘極距離的影響**取$$S=10^{-8}\ \mathrm{m^2},\quad m^*=0.1m_e,\quad g_d=2,$$$$n_H=10^{16}\ \mathrm{m^{-2}}, \qquad n_L=n_0=2.5\times10^{15}\ \mathrm{m^{-2}},$$$$a_i=100\ \mathrm{nm}, \qquad a_f=80\ \mathrm{nm}.$$以下為指定理想模型的數值積分；力差在 $a=100\ \mathrm{nm}$ 評估：

| $h$ | 背閘極靜電勢差 | $p_H-p_L$ | 未扣損耗的 $W_0$ |
| --- | --: | --: | --: |
| $10\ \mathrm{nm}$ | $1.36\ \mathrm V$ | $0.00860\ \mathrm{Pa}$ | $2.69\times10^{-15}\ \mathrm J$ |
| $100\ \mathrm{nm}$ | $13.6\ \mathrm V$ | $0.0461\ \mathrm{Pa}$ | $2.30\times10^{-16}\ \mathrm J$ |
| $1000\ \mathrm{nm}$ | $136\ \mathrm V$ | $0.0520\ \mathrm{Pa}$ | $3.12\times10^{-17}\ \mathrm J$ |

表中的電壓是靜電勢差；完整端電壓還包含上面的化學勢與交互能修正。

這一步得到的設計關係是：**小 $h$ 降低充電電壓，但固定背閘極背景會削弱薄膜的可調力差；大 $h$ 保留較大的調控幅度，卻提高控制電壓。**近閘極時較大的總回收量，主要來自背閘極—工作板的固定背景交互能。

現在，獨立控制端已對應到明確三電極幾何；新增參數 $h$ 同時決定控制電容與電磁反射，兩者可以在同一模型內數值核算。

* * *

**ChatGPT:**

接著可以把**有限控制電壓與實際能量回收要求**算出來。本例得到兩個明確限制：背閘極距離不能任意增加，控制端的額外損耗也必須低於很小的能量預算。

仍採上一輪三電極模型，工作間隙的平均靜電場保持補償為零。參數取

$$S=10^{-8}\ \mathrm{m^2}, \qquad m^*=0.1m_e,\qquad g_d=2,$$$$n_H=10^{16}\ \mathrm{m^{-2}}, \qquad n_L=n_0=2.5\times10^{15}\ \mathrm{m^{-2}},$$$$a_i=100\ \mathrm{nm},\quad b=60\ \mathrm{nm},\quad a_f=80\ \mathrm{nm}.$$

以下新增一個示範條件：

$$\boxed{V_{\max}=5\ \mathrm V}.$$

這是計算用的控制上限，尚未指定實際材料。

1. **先限制可以使用的背閘極距離**完整控制電壓為$$\boxed{ V_g(a,h,n) = \frac1e\left[ \left(\alpha+\frac{e^2h}{\varepsilon_0}\right)(n-n_0) +\mu_Q(a,h,n) \right] },$$其中$$\alpha=\frac{2\pi\hbar^2}{g_dm^*},$$$$\mu_Q= \frac{\partial}{\partial n} [u_{\mathrm{sg}}(h,n)+u(a,h,n)].$$整條操作路徑須滿足$$\boxed{ \max_{\mathrm{path}}|V_g|\le V_{\max} }.$$忽略本例較小的 $\mu_Q$ 修正時，距離上限近似為$$\boxed{ h_{\max}^{(0)} = \frac{\varepsilon_0}{e(n_H-n_0)} \left[ V_{\max}-\frac{\alpha(n_H-n_0)}e \right] }.$$使用完整反射積分求解後：$$\boxed{h_{\max}\simeq36.71\ \mathrm{nm}}.$$

| $h$ | 路徑最大準靜態電壓 | 初始點的 $p_H-p_L$ |
| --- | --: | --: |
| $10\ \mathrm{nm}$ | $1.375\ \mathrm V$ | $0.00860\ \mathrm{Pa}$ |
| $30\ \mathrm{nm}$ | $4.089\ \mathrm V$ | $0.0256\ \mathrm{Pa}$ |
| $36.71\ \mathrm{nm}$ | $5.000\ \mathrm V$ | $0.0296\ \mathrm{Pa}$ |
| $100\ \mathrm{nm}$ | $13.59\ \mathrm V$ | $0.0461\ \mathrm{Pa}$ |

因此，在這組密度與 $5\,\mathrm V$ 限制下，原定的強、弱響應可在約 $36.7\,\mathrm{nm}$ 以下達到。
2. **選 $$h=30\,\mathrm{nm}$$，完整核算一輪**使用$$E(a,n,h)=S[\Phi_0(n,h)+u(a,h,n)].$$各段功定義為$$W_H=E(a_i,n_H,h)-E(b,n_H,h),$$$$W_R=E(a_f,n_L,h)-E(b,n_L,h),$$$$\begin{aligned} W_{\mathrm{sw}}={}&E(b,n_L,h)-E(b,n_H,h)\\ &+E(a_f,n_H,h)-E(a_f,n_L,h). \end{aligned}$$數值結果為：

| 項目 | 能量 |
| --- | --: |
| 強響應壓縮取得 $W_H$ | $4.0378\times10^{-15}\ \mathrm J$ |
| 弱響應回推需要 $W_R$ | $2.7132\times10^{-15}\ \mathrm J$ |
| 兩次切換淨輸入 $W_{\mathrm{sw}}$ | $2.4423\times10^{-17}\ \mathrm J$ |

因而$$\boxed{ W_0=W_H-W_R-W_{\mathrm{sw}} =1.30018\times10^{-15}\ \mathrm J }.$$它與直接計算端點能差一致：$$W_0=S[u(a_i,h,n_H)-u(a_f,h,n_H)].$$
3. **切換的淨功小，但處理的往返能量較大**定義兩次密度切換實際處理的理想能量總量$$E_{\mathrm{proc}} = |W_{H\to L}|+|W_{L\to H}|.$$本例得到$$\boxed{ E_{\mathrm{proc}} \simeq4.9138\times10^{-11}\ \mathrm J }.$$其中大部分在減少載子時回到儲能器，恢復載子時再取出。因此，理想淨切換功很小，仍要求控制端能高度可逆地處理這筆往返能量。
4. **加入損耗，得到控制端的容許值**沿用$$T_j=0.01\ \mathrm s,\quad R_g=10\ \Omega,\quad \zeta=10^{-8}\ \mathrm{kg/s},$$則$$D_g+D_m\simeq3.49\times10^{-19}\ \mathrm J.$$只計這兩項損耗時：$$\boxed{ W_{\mathrm{net}} \simeq1.29983\times10^{-15}\ \mathrm J }.$$現在用兩個額外參數描述尚未計入的控制消耗：

* $\epsilon$：額外轉換損失占 $E_{\mathrm{proc}}$ 的比例。
* $P_{\mathrm{idle}}$：控制與補償系統的平均待機耗能率。

真正輸出便是$$\boxed{ W_{\mathrm{out}} = W_0-D_g-D_m -\epsilon E_{\mathrm{proc}} -P_{\mathrm{idle}}T_{\mathrm{cyc}} },$$其中$$T_{\mathrm{cyc}}=0.04\ \mathrm s.$$正輸出要求$$\boxed{ \epsilon E_{\mathrm{proc}} +P_{\mathrm{idle}}T_{\mathrm{cyc}} <1.29983\times10^{-15}\ \mathrm J }.$$分別令另一項為零，可得$$\boxed{\epsilon<2.65\times10^{-5}},$$或$$\boxed{ P_{\mathrm{idle}}<3.25\times10^{-14}\ \mathrm W }.$$也就是：本例允許的額外損失比例約低於 $0.00265\%$，平均待機功率預算約為 $32.5\,\mathrm{fW}$。兩者須共同滿足上面的總預算。
5. **另外辨認薄膜調控的貢獻**本例固定金屬邊界的背景端點能差為$$W_{\mathrm{bg}} = \frac{S\pi^2\hbar c}{720} \left[ (a_f+h)^{-3}-(a_i+h)^{-3} \right],$$得到$$W_{\mathrm{bg}}\simeq1.28344\times10^{-15}\ \mathrm J.$$所以總場能回收主要包含這個固定背景。若比較高、低密度在相同起終板距下的能差：$$\boxed{ \Delta W_{HL} = S\!\left[ u(a_i,h,n_H)-u(a_f,h,n_H) -u(a_i,h,n_L)+u(a_f,h,n_L) \right] },$$則$$\boxed{ \Delta W_{HL}\simeq8.37\times10^{-18}\ \mathrm J }.$$這是辨認載子調控效果的比較量，已包含於前面的能量計算中。

**目前已得到一組具體的條件式結果：**在這個理想三電極模型內，$h=30\,\mathrm{nm}$ 可在示範的 $5\,\mathrm V$ 上限內完成原定密度切換，並有有限正回收功；實際能否輸出，關鍵轉為控制端的額外轉換損失與待機消耗能否落在上述預算內。

實際選材時，載子分布與整體層結構的光學反射也須一起求解；已有 MIS–TCO 電性調控 Casimir 力研究採用這類聯合計算，但本例仍是指定理想參數下的結果。[arxiv.org](https://arxiv.org/pdf/2207.02305?utm_source=chatgpt.com)

* * *

**ChatGPT:**

若把三電極的**有效重疊面積 $S$ 放大**，並保持板距 $a$、背閘極距離 $h$、載子面密度 $n_H,n_L$ 與操作時間不變，主要結果是：

**總力與可回收能量增加，但單位面積的壓力、反射響應與理想控制電壓不變。**

這適用於目前使用的大面積平面近似，邊緣效應另行處理。[arxiv.org](https://arxiv.org/pdf/1007.1231?utm_source=chatgpt.com)

令面積放大倍率為

$$r_S=\frac{S}{S_0}.$$

| 物理量 | 公式 | 面積放大後 |
| --- | --- | --- |
| 單位面積交互能 | $u(a,h,n)$ | 不變 |
| 壓力與可調壓差 | $p,\ p_H-p_L$ | 不變 |
| 總力與可調力差 | $F=-Sp,\ \Delta F=-S\Delta p$ | 乘 $r_S$ |
| 理想回收功 | $W_0=S[u(a_i)-u(a_f)]$ | 乘 $r_S$ |
| 控制端處理能量 | $E_{\mathrm{proc}}=S e_{\mathrm{proc}}$ | 乘 $r_S$ |
| 背閘極電容 | $C_g=\varepsilon_0S/h$ | 乘 $r_S$ |
| 切換電荷量 | $\Delta Q=eS(n_H-n_L)$ | 乘 $r_S$ |
| 相同時序的控制電流 | $I_g=eS\dot n$ | 乘 $r_S$ |
| 準靜態控制電壓 | $V_g(a,h,n)$ | 不變 |
| 表面動能電感 | $L_s=m^*/(ne^2)$ | 不變 |

面積增加所帶來的集體效應，可以從模態計數看出：

$$\sum_{\mathbf k} \longrightarrow S\int\frac{d^2k}{(2\pi)^2}.$$

面積越大，參與的橫向模態總數越多；固定波數下的材料響應仍由 $n,m^*,a,h$ 決定。

**控制電路損耗的縮放需要特別計算。**

前面的切換電阻損耗為

$$D_g= \frac65R_g[eS(n_H-n_L)]^2 \left(\frac1{T_2}+\frac1{T_4}\right).$$

因此，若 $R_g$ 與切換時間都固定：

$$\boxed{D_g\propto S^2}, \qquad \boxed{W_0\propto S}.$$

電流增加也會影響驅動器端電壓：

$$V_{\mathrm{drv}} =V_g+R_gI_g+L_g\dot I_g.$$

所以理想 $V_g$ 雖不變，電路的壓降與電流需求仍會增加。

這個平方損耗律取決於驅動配置。若並聯單元使等效電阻滿足

$$R_g(S)\propto\frac1S,$$

則

$$D_g\propto S.$$

若改成保持相同最大電流，則切換時間必須隨面積增加；總能量增加，也須付出較長操作時間。

**面積放大對先前示例的影響**

把邊長從 $100\,\mu\mathrm m$ 增加為 $1\,\mathrm{mm}$，面積增加一百倍。沿用相同密度、板距、$R_g$ 與切換時間：

| 項目 | 原面積 | 面積增加一百倍 |
| --- | --: | --: |
| 理想回收功 $W_0$ | $1.30\times10^{-15}\,\mathrm J$ | $1.30\times10^{-13}\,\mathrm J$ |
| 控制處理能量 $E_{\mathrm{proc}}$ | $4.91\times10^{-11}\,\mathrm J$ | $4.91\times10^{-9}\,\mathrm J$ |
| 切換峰值電流 | $1.80\,\mathrm{nA}$ | $180\,\mathrm{nA}$ |
| 電阻損耗 $D_g$ | $3.47\times10^{-19}\,\mathrm J$ | $3.47\times10^{-15}\,\mathrm J$ |

這裡的機構阻力與其他損耗，還需要依放大後的結構重新給定。

另一方面，理想回收功與控制處理能量都按面積增加，因此

$$\boxed{ \frac{W_0}{E_{\mathrm{proc}}} =\frac{w_0}{e_{\mathrm{proc}}} }$$

不隨面積改變。先前約 $2.65\times10^{-5}$ 的理想比例仍保留；固定待機功率則能由更大的總輸出分攤。

還有兩個需要分清的情況：

* **整片電感不一定隨面積增加。**對長 $\ell$、寬 $w$ 的薄膜條，$$L_{\mathrm{strip}}=L_s\frac{\ell}{w}.$$長寬等比例放大時，$\ell/w$ 不變，因此它的端點動能電感不變。
* **若總電子數固定，前面的面積律會改變。**此時$$n=\frac NS, \qquad L_s=\frac{m^*S}{Ne^2},$$反射響應與壓力也會隨面積改變，必須重新計算$$\boxed{E_C=S\,u(a,h,N/S)}.$$

因此，對目前裝置最直接的面積設計公式是

$$\boxed{ \begin{aligned} W_{\mathrm{out}}(S)={}& S w_0-\epsilon S e_{\mathrm{proc}}\\ &-D_g(S)-D_m(S) -P_{\mathrm{idle}}(S)T_{\mathrm{cyc}}. \end{aligned} }$$

**放大面積能提高總回收量；改善多少淨輸出，則取決於控制電阻、操作時間與機構損耗如何一起縮放。**

* * *

**ChatGPT:**

之前使用的面積是

$$\boxed{S=10^{-8}\ \mathrm{m^2}}.$$

相當於邊長 $100\,\mu\mathrm m$ 的正方形：

$$100\,\mu\mathrm m\times100\,\mu\mathrm m =0.01\ \mathrm{mm^2}.$$

這是電感片、背閘極與工作反射板的有效重疊面積。