24小时热门版块排行榜    

查看: 1976  |  回复: 2

feng初心

新虫 (初入文坛)

[交流] [请教]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折射率的曲线不能够和论文中的一样。

在此向大家求助,我不知道哪里出错了。恳请诸位前辈提供一些建议,谢谢诸位。
[请教]5% MgO-PPLN晶体的折射率计算的问题
我绘制的图片.jpg


[请教]5% MgO-PPLN晶体的折射率计算的问题-1
论文的图片.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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chendequan

铁虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你是对的,论文中图表3出现错误。

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Fig._3.pdf
  • 2016-09-09 18:49:29, 147.49 K
QQ:516477448,真心帮助解决MATLAB相关问题,提供详细资料,Word文档明确具体问题及要求,尽力而为!
2楼2016-09-09 18:49:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feng初心

新虫 (初入文坛)

感谢您的回复。我因为很早注册小木虫,中间间隔了很长一段时间没有用过,所以没有金币,
等我有了金币,我会给您追加金币的,麻烦您的时间了。

我想请教一下,对于论文的更正是如何找到的,如果能够告诉我那就更好了。
3楼2016-09-10 16:14:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 feng初心 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-16 18:00 by i3KBRD19ERkG
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:15 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 16:51 by i3KBRD19ERkG
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 16:27 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 16:00 by i3KBRD19ERkG
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:51 by i3KBRD19ERkG
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-16 15:24 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-16 03:26 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见