Learning Forum

[Jeremy Brett]

您好，根据文献上得到光学手性的表达式是这样的

因此我写的脚本如下

W=matrix(length(f));

OC=-1/2*mu0*eps0*real(conj(Hx)*Ex+conj(Hy)*Ey+conj(Hz)*Ez);

for (i=1; i<= length(f); i=i+1)

{

W(i)=2*pi*f(i);

OC(:,:,1,i)=OC(:,:,1,i)*W(i);

}

不知道什么原因，这个脚本运行出来的结果就是和文献不对，想问下脚本哪里错了吗

[GuanYo Dong]

你好,

我大概查了下wiki.Electric displacement field - Wikipedia

如果在材料裡面,D用E表示,至少需要permitiffity一起表示
初步看問題可能是您的permitiffity沒有寫到?

 

莎益博,董冠佑

[Jeremy Brett]

您好，我这个前面有乘eps0，OC=-1/2*mu0*eps0*real(conj(Hx)*Ex+conj(Hy)*Ey+conj(Hz)*Ez)，感觉这个系数的问题对于图不怎么影响，我具体的图如下

文献中图是这样的，

其中OC OC1 OC2分别对应文献中沿E方向的optical chirality,OC3 是中间切面，OC4是水平切面

OC OC1 OC2 OC3均能对上

.

OC4是

就是不是全为0，和文章的不符合，不知道是什么原因