24小时热门版块排行榜

返回列表

feng初心

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 2
在线: 5小时
虫号: 3514626
注册: 2014-11-02

[交流] [请教]5% MgO-PPLN晶体的折射率计算的问题已有1人参与

各位科研前辈，大家晚上好。

我是一名硕士生，我在阅读以下论文的时候，遇到了困难，在此向大家求助关于
5% MgO-doped congruent LiNbO3的折射率计算的问题。

论文是springer出版社出版的，名字如下：
G. Gayer et. al. Temperature and wavelength dependent refractive index
equations for MgO-doped congruent and stoichiometric LiNbO3;
Appl. Phys. B . 2008

链接地址如下：
link.springer.com/article/10.1007/s00340-008-2998-2

主要目的：
我想要复现论文的Figure 3绘制的折射率随着波长变化的曲线。

采用的方法：
我根据论文Table 1中描述的参数，以及论文公式(2),(3)函数，编写了MATLAB代码。

我选择的波长范围是0.5um-4um，温度按照论文描述的室内温度选择（即T=25摄氏度）
但是绘制出来的5% MgO PPLN的ne折射率的曲线不能够和论文中的一样。

在此向大家求助，我不知道哪里出错了。恳请诸位前辈提供一些建议，谢谢诸位。

我绘制的图片.jpg

论文的图片.png

附录1：
计算代码
clc;
clear;
close all;
N=1000;    % 1000 Points
lambda = linspace(0.5e-6,4e-6,N); % unit: m
T = 20;    % degree

ne = zeros(1,N);

for i = 1:N
ne(i) = F_n_e(lambda(i),T);
end
figure();
plot(lambda,ne);
grid on;
xlim([0.5e-6,4e-6]);
ylim([2,2.35]);

附录2：
编写的计算折射率的代码
function ne_result = F_n_e(lambda_m, T )
%Sellmeier Equation's coefficient

%Sellmeier Equation  for ne in Congruent LN  %reference: Temperature and wavelength dependent refractive index equations for MgO-doped congruent and stoichiometric LiNbO3;
%Appl. Phys. B . 2008. G. Gayer et. al.
%5% MgO Doped congruent LiNbO3(CLN), HCP Corporation(Which is purchased from Crystal Technology Inc.)

% Input:
% lambda_m: Wavelength lambda in meter, unit: meter
% T: Temperature, unit: degree
% Output:
% ne_result: the refractive index

a1 = 5.756;
a2 = 0.0983;
a3 = 0.2020;
a4 = 189.32;
a5 = 12.52;
a6 = 1.32*10^(-2);

b1 = 2.860*10^(-6);
b2 = 4.700*10^(-8);
b3 = 6.113*10^(-8);
b4 = 1.516*10^(-4);

fT =(T-24.5)*(T+570.82);

lambda = lambda_m * 10^6; ;%将m表示的波长，转为 um表示的波长
ne_2  =  a1 + b1*fT + ( a2 + b2*fT) /(lambda^2-(a3 + b3*fT)^2) + ( a4 + b4 * fT )/ (lambda^2 -a5^2) - a6 * lambda^2 ;
ne_result = sqrt(ne_2);
end

回复此楼