周期性边界和布洛赫边界条件关系及使用情况

Periodic:只使用于周期结构,平面波入射情况下(以免光源被切断而产生衍射边缘效应)。光源正入射,如果光源与坐标轴有一定角度,必须在有角度的平面使用Bloch边界条件。Periodic是Bloch的特殊情况;

Bloch:是Periodic的一般形式;它是一种普遍的边界条件,由于数学上要求它只能针对指定的波长有指定的入射角,其它波长的实际入射角将不同于指定的那个入射角,因此一般情况下,它适合单波长计算。

Symmetric/anti-symmetric:对称/反对称边界条件。要求:结构对称性,光源的偏振也要对称。

  • 3
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
光学布洛赫方程表述了光在周期性介质中的传播,其公式如下: $$i\frac{\partial \psi}{\partial z} = -\frac{\hbar^2}{2n_0}\frac{\partial^2 \psi}{\partial x^2} + \left[V(x) + \frac{\hbar^2 k^2}{2n_0} \right]\psi(x)$$ 其中,$\psi(x)$表示波函数,$z$是传播距离,$x$是坐标,$n_0$是介质的折射率,$k$是波矢,$V(x)$是周期性势场。 Matlab可以通过数值方法求解光学布洛赫方程。例如,使用有限差分法(Finite Difference Method,FDM)可以将微分方程离散化,将空间和时间分为若干个网格点,并通过迭代方法求解方程。下面是一个简单的Matlab程序: ```matlab % Parameters n0 = 1.5; % Refractive index k = 2*pi/633e-9/n0; % Wavevector V0 = 5*k^2*n0; % Potential strength L = 1e-5; % Period length sigma = 0.1e-5; % Pulse width zmax = 1e-3; % Propagation distance Nz = 1000; % Number of z steps Nx = 256; % Number of x steps dx = L/Nx; % Spatial step dt = 1e-17; % Time step Tmax = 1e-14; % Maximum time Nt = round(Tmax/dt); % Number of time steps x = linspace(-L/2, L/2, Nx); % x coordinates % Initialize wavefunction nx0 = round(Nx/2); % Center index psi0 = exp(-x.^2/(2*sigma^2)); % Gaussian wavepacket psi = zeros(Nx,Nz); % Wavefunction % Solve PDE using FDM for j = 1:Nt % Finite difference scheme psi(nx0,1) = exp(-1i*k*x(nx0)*j*dt)*psi0(nx0); for n = 1:Nz-1 psi(1,n+1) = psi(1,n)*(1+1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(1))))) - 1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(2))))*psi(2,n); for k = 2:Nx-1 psi(k,n+1) = psi(k,n)*(1+1i*dt/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(k)))))-1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(k+1))))*psi(k+1,n)-1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(k-1))))*psi(k-1,n); end psi(Nx,n+1) = psi(Nx,n)*(1+1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(Nx))))) - 1i*dt/2/dx^2*(1i*k+n0*sqrt(V0/2/n0-L/2/n0*exp(-1i*2*pi/L*x(Nx-1))))*psi(Nx-1,n); end end % Plot wavefunction figure; plot(x,real(psi(:,end))); xlabel('x (m)'); ylabel('Real(psi)'); ``` 此程序通过FDM解决该方程,并产生波函数的实部的数字表示。程序允许改变各个参数,并可用于研究不同情况下的解决方案。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值