File:Microcavity dynamics.gif
Microcavity_dynamics.gif (360 × 223像素,文件大小:1.29 MB,MIME类型:image/gif、循环、215帧)
摘要
描述Microcavity dynamics.gif |
English: Transfer matrix simulation of the dynamic of the electric field when a pulse is shone on a microcavity (in this case a Bragg reflector with a defect in the middle). Most of the pulse is reflected straight away, but the frequencies resonant with the cavity couple to it and the energy is stored in the confined mode. The cavity then relaxes exponentially with a time constant that depends on the Q-factor of the resonance. |
日期 | |
来源 | https://twitter.com/j_bertolotti/status/1075341329817853952 |
作者 | Jacopo Bertolotti |
授权 (二次使用本文件) |
https://twitter.com/j_bertolotti/status/1030470604418428929 |
Mathematica 11.0 code
c = 3 10^8; (*speed of light*) M[n_, k_, d_] := {{Cos[n k d], I c/n Sin[n k d]}, {I n/c Sin[n k d], Cos[n k d]}}; (*transfer matrix*) Mi[n_, k_, d_] := {{Cos[d k n], -((I c Sin[d k n])/n)}, {-((I n Sin[d k n])/c), Cos[d k n]}}; (*Inverse of a transfer matrix*) t[m_, n0_, n2_] := (2 n0/c)/(n2/c m[[1, 1]] - (n0 n2)/c^2 m[[1, 2]] - m[[2, 1]] + n0/c m[[2, 2]]); (*transmission coefficient*) d = 1 10^-6; (*layer thickness in m*) dim = 6; (*number of layers in the Bragg mirror*) s = Join[Table[1., 50], Table[If[EvenQ[j], 1., 2.], {j, 1, dim}], {1, 1}, Table[If[EvenQ[j], 1., 2.], {j, 1, dim}], Table[1., 50]] ;(*Reflective indices of each layer (including some space to show the pulse arrive*) dim = Dimensions[s][[1]]; source = E^(-(1/2) (w - w0)^2 \[Sigma]^2) /. {w0 -> 2.185 10^15, \[Sigma] -> (10 10^-6)/c, a -> 10^12}; nstep = 2000; \[Omega]min = 1.9 10^15; \[Omega]max = 2.8 10^15; sourcel = Table[source, {w, \[Omega]min, \[Omega]max, (\[Omega]max - \[Omega]min)/nstep}]; trasm = Reap[ For[\[Omega] = \[Omega]min, \[Omega] <= \[Omega]max, \[Omega] = \[Omega] + (\[Omega]max - \[Omega]min)/nstep, tm = Apply[Dot, Table[M[s[[j]], \[Omega]/c, d], {j, 1, dim}]]; Sow[N[t[tm, 1, 1]] ]; ];][[2, 1]]; field = trasm*sourcel; (*Field at the last interface*) sexpand = 5; (*increase spatial resolution*) s2 = Flatten@Table[Table[s[[j]], sexpand], {j, 1, dim}]; freq = Table[j, {j, \[Omega]min, \[Omega]max, (\[Omega]max - \[Omega]min)/nstep}]; fn = Transpose[{field, field/c}]; tmp0 = fn; ssm = Reap[For[i = dim*sexpand, i > 0, i--, tmp = Table[((Mi[s2[[i]], freq/c, d/sexpand])[[All, All, j]].tmp0[[j]]), {j, 1, nstep}]; Sow[tmp[[All, 1]]]; tmp0 = tmp; ];][[2, 1]]; fssm = Map[Fourier, ssm]; p1 = Table[ ListPlot[{Re@Reverse@fssm[[All, -j]], Abs@Reverse@fssm[[All, -j]], -Abs@Reverse@fssm[[All, -j]]}, PlotRange -> {-7, 7}, Joined -> True, Axes -> False, PlotStyle -> {Directive[Orange], Directive[Thick, Black], Directive[Thick, Black]}, Epilog -> {Dashed, Black, Thick, Line[{{50*sexpand, -3}, {50*sexpand, 3}}], Line[{{64*sexpand, -3}, {64*sexpand, 3}}], Text[Style["Microcavity", Medium, Bold], 57*sexpand, 5}]} ], {j, -15, 200, 1}]; ListAnimate[Drop[p1, {16}], 10]
许可协议
我,本作品著作权人,特此采用以下许可协议发表本作品:
![]() ![]() |
本作品采用知识共享CC0 1.0 通用公有领域贡献许可协议授权。 |
采用本宣告发表本作品的人,已在法律允许的范围内,通过在全世界放弃其对本作品拥有的著作权法规定的所有权利(包括所有相关权利),将本作品贡献至公有领域。您可以复制、修改、传播和表演本作品,将其用于商业目的,无需要求授权。
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse |
![]() |
管理员或受信用户Ronhjones确认本图片在2018年12月21日可在下列站点找到并符合所选许可证:
https://twitter.com/j_bertolotti/status/1030470604418428929 |
说明
添加一行文字以描述该文件所表现的内容
Simulación de la matriz de transferencia de la dinámica del campo eléctrico cuando se emite un pulso en una microcavidad.
此文件中描述的项目
描绘内容
19 12 2018
文件历史
点击某个日期/时间查看对应时刻的文件。
日期/时间 | 缩略图 | 大小 | 用户 | 备注 | |
---|---|---|---|---|---|
当前 | 2018年12月20日 (四) 14:59 | ![]() | 360 × 223(1.29 MB) | Berto | User created page with UploadWizard |
文件用途
以下页面使用本文件:
全域文件用途
以下其他wiki使用此文件:
- en.wikipedia.org上的用途
- es.wikipedia.org上的用途
- fr.wikipedia.org上的用途
- lv.wikipedia.org上的用途
元数据
此文件含有额外信息,这些信息可能是创建或数字化该文件时使用的数码相机或扫描仪所添加的。如果文件已从其原始状态修改,某些详细信息可能无法完全反映修改后的文件。
GIF文件备注 | Created with the Wolfram Language : www.wolfram.com |
---|