CN117744451B - 基于有限元和分子动力学的sers基底仿真优化方法 - Google Patents
基于有限元和分子动力学的sers基底仿真优化方法 Download PDFInfo
- Publication number
- CN117744451B CN117744451B CN202410185682.1A CN202410185682A CN117744451B CN 117744451 B CN117744451 B CN 117744451B CN 202410185682 A CN202410185682 A CN 202410185682A CN 117744451 B CN117744451 B CN 117744451B
- Authority
- CN
- China
- Prior art keywords
- model
- sers substrate
- adsorption
- substrate
- rhodamine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000758 substrate Substances 0.000 title claims abstract description 189
- 238000004416 surface enhanced Raman spectroscopy Methods 0.000 title claims abstract description 145
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 238000000329 molecular dynamics simulation Methods 0.000 title claims abstract description 32
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000001179 sorption measurement Methods 0.000 claims abstract description 98
- VYXSBFYARXAAKO-WTKGSRSZSA-N chembl402140 Chemical compound Cl.C1=2C=C(C)C(NCC)=CC=2OC2=C\C(=N/CC)C(C)=CC2=C1C1=CC=CC=C1C(=O)OCC VYXSBFYARXAAKO-WTKGSRSZSA-N 0.000 claims abstract description 65
- 230000005684 electric field Effects 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000012805 post-processing Methods 0.000 claims abstract description 16
- 238000010219 correlation analysis Methods 0.000 claims abstract description 12
- 239000000463 material Substances 0.000 claims description 67
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 45
- 239000010931 gold Substances 0.000 claims description 36
- 229910052737 gold Inorganic materials 0.000 claims description 36
- 230000003993 interaction Effects 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 20
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 14
- 230000000007 visual effect Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 239000013078 crystal Substances 0.000 claims description 6
- VYXSBFYARXAAKO-UHFFFAOYSA-N ethyl 2-[3-(ethylamino)-6-ethylimino-2,7-dimethylxanthen-9-yl]benzoate;hydron;chloride Chemical compound [Cl-].C1=2C=C(C)C(NCC)=CC=2OC2=CC(=[NH+]CC)C(C)=CC2=C1C1=CC=CC=C1C(=O)OCC VYXSBFYARXAAKO-UHFFFAOYSA-N 0.000 claims description 6
- 239000002086 nanomaterial Substances 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 6
- 238000002493 microarray Methods 0.000 claims description 4
- 238000013433 optimization analysis Methods 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000002939 conjugate gradient method Methods 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 230000010363 phase shift Effects 0.000 claims description 3
- 238000005381 potential energy Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 12
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000001069 Raman spectroscopy Methods 0.000 abstract description 7
- 238000010276 construction Methods 0.000 abstract description 2
- 238000001514 detection method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- 241000682719 Adina Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000001237 Raman spectrum Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 238000007431 microscopic evaluation Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于有限元仿真优化技术领域,具体涉及基于有限元和分子动力学的SERS基底仿真优化方法,步骤包括通过有限元仿真软件构建包括上部介质区域和下部SERS基底区域的整体物理模型;通过模型开发器的接口和模块设置相关参数;在模型开发器中划分网格,进行求解计算和后处理得到电场增强因子;构建分子动力学模型和吸附模型并求解得到吸附信息;进行相关性分析获得参数优化模型。本发明充分考虑了实际SERS基底结构的有限元物理模型构建与电磁分析,并且充分考虑了罗丹明6G分子构型和基底结构,进而建立了分子吸附的分子动力学模型,能够较为准确的解析SERS基底的拉曼增强机制,并为基底结构优化提供了量化数据。
Description
技术领域
本发明属于有限元仿真优化技术领域,具体涉及基于有限元和分子动力学的SERS基底仿真优化方法。
背景技术
表面增强拉曼散射(Surface Enhanced Raman Scattering, SERS)技术是一种新兴的表面光谱分析技术。SERS技术在分析过程中无需对样品进行复杂的预处理,通过拉曼光谱的指纹特性可以对物质进行快速检测与识别。同时,SERS基底表面的“热点”区域能够将待测物质的信号强度放大百万倍以上,有效解决待测物质信号弱和待测组分含量低的问题。
近年来,SERS技术已成为环境检测、食品安全、物证检验等痕量检测领域极具应用潜力和研究价值的方法之一。此外,它在医学和生物技术领域具有广阔的前景。然而对于SERS技术,高性能活性SERS基底是实现高精度痕量检测的关键,因此如何构建金属基底结构和优化成为SERS领域的研究焦点。为了实现具有高灵敏性、均匀性、稳定性和再现性的SERS活性基底结构优化,需要开发一种基于有限元分析和分子动力学模拟的SERS基底结构优化方法,从宏观电磁增强和微观分析吸附两个角度对活性基底的拉曼增强机理进行表征,进而提出SERS基底优化策略。
为了揭示活性基底形貌特征对SERS增强程度的影响,需要对活性基底的拉曼增强机制进行模拟。目前对于SERS基底的拉曼增强机理研究多基于理想化的、简化的二维或三维模型进行单一的物理或者化学的定性分析模拟,不能够为优化SERS基底结构提供参数优化的定量数据。因此对于制备具有高性能的SERS基底,从实际的基底结构出发构造物理模型,进而通过宏观有限元方法和微观分子动力学方法分别研究其电磁增强和分子吸附机制,最终提出SERS基底结构优化方案是目前亟待解决的问题。
对于SERS基底结构优化,其中涉及到的电磁增强机制宏观仿真以及分子吸附机制微观仿真。
电磁增强机制宏观仿真可以采用有限元分析方法,有限元分析(FEA,FiniteElement Analysis)是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。利用简单而又相互作用的元素(即单元),就可以用有限数量的未知量去逼近无限未知量的真实系统。有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段,被广泛的应用于求解电磁场、热传导、流体力学等连续性问题。有限元分析软件目前最流行的有ANSYS、ADINA、ABAQUS、MSC、Hyperworks、Fluent、Comsol等。
分子吸附机制微观仿真可采用分子动力学模拟方法,分子动力学是一门结合物理、数学和化学的综合技术。该方法主要是依靠牛顿力学来模拟分子体系的运动,以在由分子体系的不同状态构成的系统中抽取样本,从而计算体系的构型积分,并以构型积分的结果为基础进一步计算体系的热力学量和其他宏观性质。
发明内容
根据以上现有技术中的不足,本发明的目的在于提供一种基于有限元和分子动力学的SERS基底仿真优化方法,充分考虑了实际SERS基底结构的有限元物理模型构建与电磁分析,并且充分考虑了罗丹明6G分子构型和基底结构,进而建立了分子吸附的分子动力学模型,能够较为准确的解析SERS基底的拉曼增强机制,并为基底结构优化提供了量化数据。
为达到以上目的,本发明提供了基于有限元和分子动力学的SERS基底仿真优化方法,包括以下步骤:
S1、利用原子力显微镜,通过形貌扫描获得SERS基底的三维结构特征数据,包括SERS基底的微米结构形状、微米阵列结构间距和纳米结构尺寸;
S2、通过有限元仿真软件,基于获得的SERS基底的三维结构特征数据,应用可视化点云方法构建三维表面,进而应用有限元仿真软件的模型开发器中的实体转换构造SERS基底三维实体物理模型,随后构建包括上部介质区域和下部SERS基底区域的整体物理模型;
S3、通过模型开发器的材料接口,设置整体物理模型中上部介质和下部SERS基底的材料,并通过模型开发器的虚构域接口设置完美匹配层;
S4、在模型开发器的波动光学物理场电磁波-频域模块中,设置仿真参数,仿真参数包括入射光波长和功率;
S5、在模型开发器中划分网格,并完成求解器设置,进行求解计算,得到SERS基底电场分布云图,并进行求解结果后处理,得到电场增强因子;
S6、通过Materials studio软件,依据S2中的SERS基底三维实体物理模型,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型;
S7、通过Materials studio软件,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中,对吸附模型进行边界条件设置;
S8、选取吸附模型的势函数,并进行Lammps软件的求解设置,求解计算后,进行结果后处理,得到吸附模型的吸附信息;
S9、结合SERS基底的三维结构特征数据,对S5中得到的电场增强因子和S8中得到的吸附信息进行相关性分析,并对SERS基底进行结构和尺寸优化分析,得到SERS基底特征尺寸参数的最优解,进而获得参数优化模型。
最终获取的参数优化模型即为SERS基底最优的优化方案,依据其中的特征尺寸参数即可获得最优的基底结构设计方案。
所述的S1中,利用原子力显微镜扫描SERS基底xy平面的90μm×90μm区域(选取区域的尺寸也可以根据需求进行调整,下同),获得SERS基底表面三维结构的坐标数据。
所述的S2中,通过COMSOL有限元仿真软件,基于SERS基底表面三维结构的坐标数据,应用可视化点云方法构建三维表面,随后利用COMSOL的模型开发器,右击“几何”按钮,选择“转换实体”构造SERS基底三维实体物理模型,其几何总长、宽、高分别为90μm、90μm、5μm;
利用模型开发器,右击“几何”按钮,建立90μm×90μm×20μm的三维区域,即为包括上部介质区域和下部SERS基底区域的整体物理模型,其中下部SERS基底区域即为SERS基底三维实体物理模型占据的区域。
也可采用ANSYS、ADINA、ABAQUS、MSC、Hyperworks、Fluent等有限元仿真软件。
所述的S3中,对于上部介质区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为空气;对于下部SERS基底区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为金;
在模型开发器的虚构域接口中,设置整体物理模型各个面的最外层为完美匹配层。
所述的S4中,通过模型开发器的波动光学物理场电磁波-频域模块,分别设置入射光波长和功率,并使用近轴高斯光束近似,光束沿Z轴负向传播,其背景场方程如下:
(1);
式中,Eb为入射光背景场;Ebg0为高斯光束电场振幅;w0为光束半径;w(z)为光束在传播轴不同位置处的光斑半径;k为波数;R(z)为光束的相位曲率半径;η(z)为Gouy相位移;j为虚数单位;x,y,z即为空间直角坐标;
其中,w(z),R(z)和η(z)的表达式为:
(2);
式中,p0为Z轴上的焦平面;z0为瑞利长度,其获取公式为:
(3);
散射电场的求解取决于麦克斯韦波动方程,该方程表示为:
(4);
式中,为梯度算子;Esca为散射总电场;μr为介质的相对渗透率;k0为自由空间中的波数;εr为介质的介电常数,ε0为真空中介电常数,σ为电导率,ω为散射光波角频率。
所述的S5中,在模型开发器中划分网格、进行求解计算以及求解结果后处理的步骤为:
S51、在模型开发器中划分网格:在构建的上部介质区域和下部SERS基底区域内划分非均匀网格,划分过程中,对下部SERS基底区域表面进行网格细化,同时对该区域进行区域边界粗化;计算网格时,在“模型开发器”的“网格”中选择“自适应网格”;
S52、完成求解器设置并进行求解计算:在模型开发器的“研究1”节点下,右击“参数化扫描”在“波长”列表下设置波长为300nm到1200nm,并设置步长,单击“计算”,展开“结果”列表,得到SERS基底电场分布云图;
S53、求解后处理:对SERS基底电场分布云图进行处理得到电场增强因子,采用的公式为:
(5);
式中,I为电场增强因子;E0为初始电场。
所述的S6中,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型的步骤为:
S61、构建SERS基底分子动力学模型:在Materials studio软件的晶体库中导入“金”晶体模型,使用“创建”模块中“对称”列表下的“扣除”功能进行扩胞,依据S2中的SERS基底三维实体物理模型的尺寸(90μm、90μm、5μm),将其缩小3300倍,并结合金晶胞大小,得到体积为274.246Å×274.246Å×125Å的金材料基底;
分子动力学模型是利用了S1中形貌扫描获得的SERS基底表面三维结构的相似结构,即形貌结构一样,但应用的尺寸从微米(μm)级别缩小到了埃(Å)级别,因为分子动力学方法很难根据实际尺寸来模拟模型,实际尺寸模拟会极大地消耗计算机资源并且受到其复杂性的限制。因此,在本发明中,罗丹明6G分子在具有三维结构基底上的吸附模型为尺寸按比例缩小到埃级的模型。
使用“编辑”模块中“原子选择”列表下的“径向距离”功能选择半径为36Å的半球,以距离x和y方向为48Å的原子为初始中心原子进行雕刻,其余中心原子距离初始中心原子的距离为36Å,雕刻过程与初始位置雕刻过程一致,得到具有实际三维结构相似的金材料基底;
S62、构建罗丹明6G分子在具有三维结构基底上的吸附模型:通过Materialsstudio软件中的“3D原子”模块,结合罗丹明6G分子的分子式,构建罗丹明6G分子的分子模型,并将构建完成的50个罗丹明6G分子随机分布在具有三维结构的金材料基底上方7Å的空间内,获得总范围为274.246Å×274.246Å×25Å的空间,得到罗丹明6G分子在具有三维结构基底上的吸附模型,点击“力场设置”模块将吸附模型的力场设置为“CVFF力场”,电荷设置为“中性”。
所述的S7中,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中对吸附模型进行边界条件设置的步骤为:
S71、转化Lammps软件识别文件:通过Materials studio软件,将吸附模型文件输出为CAR格式进行保存,通过Lammps软件中的“msi2lmp”程序,将CAR格式的吸附模型文件转化为Lammps软件可识别的data格式文件;
S72、对吸附模型进行边界条件设置:将data格式的吸附模型文件导入Lammps软件中,对吸附模型在x、y方向上设置周期性边界条件,z方向上设置固定边界条件并设置反射墙,对三维结构的金材料基底进行固定,并对罗丹明6G分子采用共轭梯度法进行能量最小化处理。能量最小化处理可以使罗丹明6G分子体系的能量降到一个较低的水平,便于后续计算。
所述的S8中,选取吸附模型的势函数,并进行Lammps软件的求解设置和求解计算,以及求解结果后处理的步骤为:
S81、势函数的选取:用混合势来描述不同种原子间的相互作用,选取“Morse”势函数来模拟金原子之间的相互作用,选取“CVFF力场”模拟罗丹明6G分子之间的相互作用及金原子与罗丹明6G分子之间的相互作用,采用“Lennard-Jones”势函数混合规则计算金原子和罗丹明6G分子的CVFF势的能量和距离参数;
S82、求解设置:求解时,采用Nosé-Hoover热浴法进行温度控制,采用PPPM方法对粒子之间的长程静电作用进行处理,采用速度verlet算法求解粒子牛顿运动方程;
S83、求解计算:使用“lammps compute”命令计算整组罗丹明6G分子与三维结构的金材料基底的相互作用能,使用“compute msd”命令计算整组罗丹明6G分子的均方位移,其中相互作用能的计算公式如式(6)所示,均方位移的计算公式如式(7)所示:
(6);
式中,ER6G+substrate为罗丹明6G分子吸附在三维结构的金材料基底上的总能量;ER6G和Esubstrate分别为罗丹明6G分子的能量和三维结构的金材料基底的能量;
(7);
式中,ri、ri,ref分别为吸附模型中原子i的瞬时位置向量及其参考位置向量,N为吸附模型中的原子总数;
计算相互作用能和均方位移后,采用“NVT”系综进行分子动力学模拟,并得到吸附模型的相互作用能、均方位移和轨迹文件;
S84、求解后处理:根据得到的相互作用能、均方位移和轨迹文件,获取罗丹明6G分子在三维结构的金材料基底上的吸附现象,得到罗丹明6G分子在三维结构的金材料基底上的吸附信息,使用ovito软件对lammps输出的轨迹文件进行可视化处理,得到罗丹明6G分子在不同时间的模拟结果的图片。
所述的S9中,进行相关性分析的方法为,通过改变SERS基底的三维结构,并在每次改变时,均重复S1-S8步骤,获取对应的SERS基底的三维结构特征数据、电磁增强因子和吸附信息,将SERS基底的微米结构形状、微米阵列结构间距、纳米结构尺寸3个结构特征的参数变化,与电磁增强因子和吸附信息进行相关性分析,并得出结构参数对电磁增强和罗丹明6G分子吸附的影响程度。
本发明使用的有限元仿真软件和分子动力学软件可以通过电子设备执行,电子设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,通过处理器执行软件实现上述的仿真、模拟和计算。
本发明所具有的有益效果是:
本发明应用了有限元分析软件和分子动力学软件,从SERS基底电磁增强机制和罗丹明6G分子吸附SERS基底机制两个方面对SERS基底结构进行仿真优化分析,实现了SERS基底电磁增强“热点”可视化,以及罗丹明6G分子吸附SERS基底过程可视化,并通过相关性分析确定了对SERS基底结构性能影响明显的特征参数,从而可以更加精确的提供优化设计的方向,进而提出SERS基底结构优化策略。
本发明充分考虑了SERS基底形貌对拉曼增强的影响,分别从宏观电磁增强和微观分子吸附两个角度进行仿真,能够较为准确计算出各种结构对表面增强拉曼的影响。
本发明可以为分析SERS基底电磁增强“热点”分布、增强因子、吸附强度等提供解决方案,从理论上提供了SERS基底结构设计方向和优化设计的合理方法,为提高SERS检测精度、拓宽SERS应用提供了理论依据。
附图说明
图1是本发明的流程原理图;
图2是本发明实施例中包括上部介质区域和下部SERS基底区域的整体物理模型的示意图,图2中的(a)为整体物理模型的立体示意图,图2中的(b)为整体物理模型的侧面示意图;
图3是本发明实施例中SERS基底电场分布结果图;
图4是本发明实施例中罗丹明6G分子在具有三维结构基底上的吸附模型的示意图;图4中的(a)为吸附模型的立体示意图,图4中的(b)为吸附模型的侧面示意图;
图5是本发明实施例中罗丹明6G分子在不同时刻的吸附情况;
图6是本发明实施例中SERS基底参数化特征的相关性分析结果图。
具体实施方式
下面结合附图对本发明的实施例做进一步描述:
如图1所示,基于有限元和分子动力学的SERS基底仿真优化方法包括以下步骤:
S1、利用原子力显微镜,通过形貌扫描获得SERS基底的三维结构特征数据,包括SERS基底的微米结构形状、微米阵列结构间距和纳米结构尺寸;
S2、通过有限元仿真软件,基于获得的SERS基底的三维结构特征数据,应用可视化点云方法构建三维表面,进而应用有限元仿真软件的模型开发器中的实体转换构造SERS基底三维实体物理模型,随后构建包括上部介质区域和下部SERS基底区域的整体物理模型;
S3、通过模型开发器的材料接口,设置整体物理模型中上部介质和下部SERS基底的材料,并通过模型开发器的虚构域接口设置完美匹配层;
S4、在模型开发器的波动光学物理场电磁波-频域模块中,设置仿真参数,仿真参数包括入射光波长和功率;
S5、在模型开发器中划分网格,并完成求解器设置,进行求解计算,得到SERS基底电场分布云图,并进行求解结果后处理,得到电场增强因子;
S6、通过Materials studio软件,依据S2中的SERS基底三维实体物理模型,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型;
S7、通过Materials studio软件,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中,对吸附模型进行边界条件设置;
S8、选取吸附模型的势函数,并进行Lammps软件的求解设置,求解计算后,进行结果后处理,得到吸附模型的吸附信息;
S9、结合SERS基底的三维结构特征数据,对S5中得到的电场增强因子和S8中得到的吸附信息进行相关性分析,并对SERS基底进行结构和尺寸优化分析,得到SERS基底特征尺寸参数的最优解,进而获得参数优化模型。
S1中,利用原子力显微镜扫描SERS基底xy平面的90μm×90μm区域,获得SERS基底表面三维结构的坐标数据。
S2中,通过COMSOL有限元仿真软件,基于SERS基底表面三维结构的坐标数据,应用可视化点云方法构建三维表面,随后利用COMSOL的模型开发器,右击“几何”按钮,选择“转换实体”构造SERS基底三维实体物理模型,其几何总长、宽、高分别为90μm、90μm、5μm;
利用模型开发器,右击“几何”按钮,建立90μm×90μm×20μm的三维区域,即为包括上部介质区域和下部SERS基底区域的整体物理模型,如图2中的(a)所示,其中下部SERS基底区域即为SERS基底三维实体物理模型占据的区域。
S3中,对于上部介质区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为空气;对于下部SERS基底区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为金;
在模型开发器的虚构域接口中,设置整体物理模型各个面的最外层为完美匹配层,完美匹配层的厚度为1500nm,如图2中的(b)所示,为了便于计算,将整体物理模型沿x、y方向上的四周边界均设置为周期性边界。
S4中,通过模型开发器的波动光学物理场电磁波-频域模块,分别设置入射光波长和功率为633nm和0.6MW,并使用近轴高斯光束近似,光束沿Z轴负向传播,其背景场方程如下:
(1);
式中,Eb为入射光背景场;Ebg0为高斯光束电场振幅;w0为光束半径;w(z)为光束在传播轴不同位置处的光斑半径;k为波数;R(z)为光束的相位曲率半径;η(z)为Gouy相位移;j为虚数单位;x,y,z即为空间直角坐标;
其中,w(z),R(z)和η(z)的表达式为:
(2);
式中,p0为Z轴上的焦平面;z0为瑞利长度,其获取公式为:
(3);
散射电场的求解取决于麦克斯韦波动方程,该方程表示为:
(4);
式中,为梯度算子;Esca为散射总电场;μr为介质的相对渗透率;k0为自由空间中的波数;εr为介质的介电常数,ε0为真空中介电常数,σ为电导率,ω为散射光波角频率。
S5中,在模型开发器中划分网格、进行求解计算以及求解结果后处理的步骤为:
S51、在模型开发器中划分网格:在构建的上部介质区域和下部SERS基底区域内划分非均匀网格,划分过程中,对下部SERS基底区域表面进行网格细化,同时对该区域进行区域边界粗化;计算网格时,在“模型开发器”的“网格”中选择“自适应网格”;
S52、完成求解器设置并进行求解计算:在模型开发器的“研究1”节点下,右击“参数化扫描”在“波长”列表下设置波长为300nm到1200nm,并设置步长为20nm,单击“计算”,展开“结果”列表,得到SERS基底电场分布云图,如图3所示;图3中,深色代表电场强度高,即电场被增强的部位,也即“热点部位”;
S53、求解后处理:对SERS基底电场分布云图进行处理得到电场增强因子,采用的公式为:
(5);
式中,I为电场增强因子;E0为初始电场。
S6中,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型的步骤为:
S61、构建SERS基底分子动力学模型:在Materials studio软件的晶体库中导入“金”晶体模型,使用“创建”模块中“对称”列表下的“扣除”功能进行扩胞,依据S2中的SERS基底三维实体物理模型的尺寸,将其缩小3300倍,并结合金晶胞大小,得到体积为274.246Å×274.246Å×125Å的金材料基底;
使用“编辑”模块中“原子选择”列表下的“径向距离”功能选择半径为36Å的半球,以距离x和y方向为48Å的原子为初始中心原子进行雕刻,其余中心原子距离初始中心原子的距离为36Å,雕刻过程与初始位置雕刻过程一致,得到具有实际三维结构相似的金材料基底;
S62、构建罗丹明6G分子在具有三维结构基底上的吸附模型:通过Materialsstudio软件中的“3D原子”模块,结合罗丹明6G分子的分子式,构建罗丹明6G分子的分子模型,并将构建完成的50个罗丹明6G分子随机分布在具有三维结构的金材料基底上方7Å的空间内,获得总范围为274.246Å×274.246Å×25Å的空间,得到罗丹明6G分子在具有三维结构基底上的吸附模型,如图4中的(a)和图4中的(b)所示;点击“力场设置”模块将吸附模型的力场设置为“CVFF力场”,电荷设置为“中性”。
S7中,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中对吸附模型进行边界条件设置的步骤为:
S71、转化Lammps软件识别文件:通过Materials studio软件,将吸附模型文件输出为CAR格式进行保存,通过Lammps软件中的“msi2lmp”程序,将CAR格式的吸附模型文件转化为Lammps软件可识别的data格式文件;
S72、对吸附模型进行边界条件设置:将data格式的吸附模型文件导入Lammps软件中,对吸附模型在x、y方向上设置周期性边界条件,z方向上设置固定边界条件并设置反射墙,对三维结构的金材料基底进行固定,并对罗丹明6G分子采用共轭梯度法进行10000步的能量最小化处理。
S8中,选取吸附模型的势函数,并进行Lammps软件的求解设置和求解计算,以及求解结果后处理的步骤为:
S81、势函数的选取:用混合势来描述不同种原子间的相互作用,选取“Morse”势函数来模拟金原子之间的相互作用,选取“CVFF力场”模拟罗丹明6G分子之间的相互作用及金原子与罗丹明6G分子之间的相互作用,采用“Lennard-Jones”势函数混合规则计算金原子和罗丹明6G分子的CVFF势的能量和距离参数;其中,“Morse”势函数和“Lennard-Jones”势函数的截止距离分别为8Å和10Å;
S82、求解设置:求解时,采用Nosé-Hoover热浴法进行温度控制,模拟温度为300K,采用PPPM方法对粒子之间的长程静电作用进行处理,采用速度verlet算法求解粒子牛顿运动方程;
S83、求解计算:使用“lammps compute”命令计算整组罗丹明6G分子与三维结构的金材料基底的相互作用能,使用“compute msd”命令计算整组罗丹明6G分子的均方位移,其中相互作用能的计算公式如式(6)所示,均方位移的计算公式如式(7)所示:
(6);
式中,ER6G+substrate为罗丹明6G分子吸附在三维结构的金材料基底上的总能量;ER6G和Esubstrate分别为罗丹明6G分子的能量和三维结构的金材料基底的能量;
(7);
式中,ri、ri,ref分别为吸附模型中原子i的瞬时位置向量及其参考位置向量,N为吸附模型中的原子总数;
计算相互作用能和均方位移后,采用“NVT”系综进行500ps分子动力学模拟,并得到吸附模型的相互作用能、均方位移和轨迹文件;
S84、求解后处理:根据得到的相互作用能、均方位移和轨迹文件,获取罗丹明6G分子在三维结构的金材料基底上的吸附现象,得到罗丹明6G分子在三维结构的金材料基底上的吸附信息,使用ovito软件对lammps输出的轨迹文件进行可视化处理,得到罗丹明6G分子在不同时间的模拟结果的图片,如图5所示。
S9中,进行相关性分析的方法为,通过改变SERS基底的三维结构,并在每次改变时,均重复S1-S8步骤,获取对应的SERS基底的三维结构特征数据、电磁增强因子和吸附信息,将SERS基底的微米结构形状、微米阵列结构间距、纳米结构尺寸3个结构特征的参数变化,与电磁增强因子和吸附信息进行相关性分析,并得出结构参数对电磁增强和罗丹明6G分子吸附的影响程度,如图6所示,影响程度即为图6中的敏感程度。
Claims (4)
1.一种基于有限元和分子动力学的SERS基底仿真优化方法,其特征在于包括以下步骤:
S1、利用原子力显微镜,通过形貌扫描获得SERS基底的三维结构特征数据,包括SERS基底的微米结构形状、微米阵列结构间距和纳米结构尺寸;
S2、通过有限元仿真软件,基于获得的SERS基底的三维结构特征数据,应用可视化点云方法构建三维表面,进而应用有限元仿真软件的模型开发器中的实体转换构造SERS基底三维实体物理模型,随后构建包括上部介质区域和下部SERS基底区域的整体物理模型;
S3、通过模型开发器的材料接口,设置整体物理模型中上部介质和下部SERS基底的材料,并通过模型开发器的虚构域接口设置完美匹配层;
S4、在模型开发器的波动光学物理场电磁波-频域模块中,设置仿真参数,仿真参数包括入射光波长和功率;
S5、在模型开发器中划分网格,并完成求解器设置,进行求解计算,得到SERS基底电场分布云图,并进行求解结果后处理,得到电场增强因子;
S6、通过Materials studio软件,依据S2中的SERS基底三维实体物理模型,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型;
S7、通过Materials studio软件,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中,对吸附模型进行边界条件设置;
S8、选取吸附模型的势函数,并进行Lammps软件的求解设置,求解计算后,进行结果后处理,得到吸附模型的吸附信息;
S9、结合SERS基底的三维结构特征数据,对S5中得到的电场增强因子和S8中得到的吸附信息进行相关性分析,并对SERS基底进行结构和尺寸优化分析,得到SERS基底特征尺寸参数的最优解,进而获得参数优化模型;
所述的S2中,通过COMSOL有限元仿真软件,基于SERS基底表面三维结构的坐标数据,应用可视化点云方法构建三维表面,随后利用COMSOL的模型开发器,右击“几何”按钮,选择“转换实体”构造SERS基底三维实体物理模型,其几何总长、宽、高分别为90μm、90μm、5μm;
利用模型开发器,右击“几何”按钮,建立90μm×90μm×20μm的三维区域,即为包括上部介质区域和下部SERS基底区域的整体物理模型,其中下部SERS基底区域即为SERS基底三维实体物理模型占据的区域;
所述的S4中,通过模型开发器的波动光学物理场电磁波-频域模块,分别设置入射光波长和功率,并使用近轴高斯光束近似,光束沿Z轴负向传播,其背景场方程如下:
(1);
式中,Eb为入射光背景场;Ebg0为高斯光束电场振幅;w0为光束半径;w(z)为光束在传播轴不同位置处的光斑半径;k为波数;R(z)为光束的相位曲率半径;η(z)为Gouy相位移;j为虚数单位;x,y,z即为空间直角坐标;
其中,w(z),R(z)和η(z)的表达式为:
(2);
式中,p0为Z轴上的焦平面;z0为瑞利长度,其获取公式为:
(3);
散射电场的求解取决于麦克斯韦波动方程,该方程表示为:
(4);
式中,为梯度算子;Esca为散射总电场;μr为介质的相对渗透率;k0为自由空间中的波数;εr为介质的介电常数,ε0为真空中介电常数,σ为电导率,ω为散射光波角频率;
所述的S5中,在模型开发器中划分网格、进行求解计算以及求解结果后处理的步骤为:
S51、在模型开发器中划分网格:在构建的上部介质区域和下部SERS基底区域内划分非均匀网格,划分过程中,对下部SERS基底区域表面进行网格细化,同时对该区域进行区域边界粗化;计算网格时,在“模型开发器”的“网格”中选择“自适应网格”;
S52、完成求解器设置并进行求解计算:在模型开发器的“研究1”节点下,右击“参数化扫描”在“波长”列表下设置波长为300nm到1200nm,并设置步长,单击“计算”,展开“结果”列表,得到SERS基底电场分布云图;
S53、求解后处理:对SERS基底电场分布云图进行处理得到电场增强因子,采用的公式为:
(5);
式中,I为电场增强因子;E0为初始电场;
所述的S6中,构建SERS基底分子动力学模型以及罗丹明6G分子在具有三维结构基底上的吸附模型的步骤为:
S61、构建SERS基底分子动力学模型:在Materials studio软件的晶体库中导入“金”晶体模型,使用“创建”模块中“对称”列表下的“扣除”功能进行扩胞,依据S2中的SERS基底三维实体物理模型的尺寸,将其缩小3300倍,并结合金晶胞大小,得到体积为274.246Å×274.246Å×125Å的金材料基底;
使用“编辑”模块中“原子选择”列表下的“径向距离”功能选择半径为36Å的半球,以距离x和y方向为48Å的原子为初始中心原子进行雕刻,其余中心原子距离初始中心原子的距离为36Å,雕刻过程与初始位置雕刻过程一致,得到具有实际三维结构相似的金材料基底;
S62、构建罗丹明6G分子在具有三维结构基底上的吸附模型:通过Materials studio软件中的“3D原子”模块,结合罗丹明6G分子的分子式,构建罗丹明6G分子的分子模型,并将构建完成的50个罗丹明6G分子随机分布在具有三维结构的金材料基底上方7Å的空间内,获得总范围为274.246Å×274.246Å×25Å的空间,得到罗丹明6G分子在具有三维结构基底上的吸附模型,点击“力场设置”模块将吸附模型的力场设置为“CVFF力场”,电荷设置为“中性”;
所述的S7中,将建立的吸附模型文件转化为Lammps软件识别文件,并将其导入Lammps软件中对吸附模型进行边界条件设置的步骤为:
S71、转化Lammps软件识别文件:通过Materials studio软件,将吸附模型文件输出为CAR格式进行保存,通过Lammps软件中的“msi2lmp”程序,将CAR格式的吸附模型文件转化为Lammps软件可识别的data格式文件;
S72、对吸附模型进行边界条件设置:将data格式的吸附模型文件导入Lammps软件中,对吸附模型在x、y方向上设置周期性边界条件,z方向上设置固定边界条件并设置反射墙,对三维结构的金材料基底进行固定,并对罗丹明6G分子采用共轭梯度法进行能量最小化处理;
所述的S8中,选取吸附模型的势函数,并进行Lammps软件的求解设置和求解计算,以及求解结果后处理的步骤为:
S81、势函数的选取:用混合势来描述不同种原子间的相互作用,选取“Morse”势函数来模拟金原子之间的相互作用,选取“CVFF力场”模拟罗丹明6G分子之间的相互作用及金原子与罗丹明6G分子之间的相互作用,采用“Lennard-Jones”势函数混合规则计算金原子和罗丹明6G分子的CVFF势的能量和距离参数;
S82、求解设置:求解时,采用Nosé-Hoover热浴法进行温度控制,采用PPPM方法对粒子之间的长程静电作用进行处理,采用速度verlet算法求解粒子牛顿运动方程;
S83、求解计算:使用“lammps compute”命令计算整组罗丹明6G分子与三维结构的金材料基底的相互作用能,使用“compute msd”命令计算整组罗丹明6G分子的均方位移,其中相互作用能的计算公式如式(6)所示,均方位移的计算公式如式(7)所示:
(6);
式中,ER6G+substrate为罗丹明6G分子吸附在三维结构的金材料基底上的总能量;ER6G和Esubstrate分别为罗丹明6G分子的能量和三维结构的金材料基底的能量;
(7);
式中,ri、ri,ref分别为吸附模型中原子i的瞬时位置向量及其参考位置向量,N为吸附模型中的原子总数;
计算相互作用能和均方位移后,采用“NVT”系综进行分子动力学模拟,并得到吸附模型的相互作用能、均方位移和轨迹文件;
S84、求解后处理:根据得到的相互作用能、均方位移和轨迹文件,获取罗丹明6G分子在三维结构的金材料基底上的吸附现象,得到罗丹明6G分子在三维结构的金材料基底上的吸附信息,使用ovito软件对lammps输出的轨迹文件进行可视化处理,得到罗丹明6G分子在不同时间的模拟结果的图片。
2.根据权利要求1所述的基于有限元和分子动力学的SERS基底仿真优化方法,其特征在于:所述的S1中,利用原子力显微镜扫描SERS基底xy平面的90μm×90μm区域,获得SERS基底表面三维结构的坐标数据。
3.根据权利要求1所述的基于有限元和分子动力学的SERS基底仿真优化方法,其特征在于:所述的S3中,对于上部介质区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为空气;对于下部SERS基底区域,在模型开发器中右击“材料”按钮,选择“从库中添加材料”,材料设置为金;
在模型开发器的虚构域接口中,设置整体物理模型各个面的最外层为完美匹配层。
4.根据权利要求1所述的基于有限元和分子动力学的SERS基底仿真优化方法,其特征在于:所述的S9中,进行相关性分析的方法为,通过改变SERS基底的三维结构,并在每次改变时,均重复S1-S8步骤,获取对应的SERS基底的三维结构特征数据、电磁增强因子和吸附信息,将SERS基底的微米结构形状、微米阵列结构间距、纳米结构尺寸3个结构特征的参数变化,与电磁增强因子和吸附信息进行相关性分析,并得出结构参数对电磁增强和罗丹明6G分子吸附的影响程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410185682.1A CN117744451B (zh) | 2024-02-20 | 2024-02-20 | 基于有限元和分子动力学的sers基底仿真优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410185682.1A CN117744451B (zh) | 2024-02-20 | 2024-02-20 | 基于有限元和分子动力学的sers基底仿真优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117744451A CN117744451A (zh) | 2024-03-22 |
CN117744451B true CN117744451B (zh) | 2024-05-14 |
Family
ID=90257783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410185682.1A Active CN117744451B (zh) | 2024-02-20 | 2024-02-20 | 基于有限元和分子动力学的sers基底仿真优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117744451B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002071013A1 (en) * | 2001-03-01 | 2002-09-12 | New Mexico State University Technology Transfer Corporation | Optical devices and methods employing nanoparticles, microcavities, and semicontinuous metal films |
CN101487832A (zh) * | 2009-02-10 | 2009-07-22 | 哈尔滨工业大学 | 一种表征环氧树脂复合材料界面机理的方法 |
CN102026719A (zh) * | 2008-03-11 | 2011-04-20 | 免疫之光有限责任公司 | 从外源进行内部能量活化的等离子体激元辅助系统和方法 |
CN102680453A (zh) * | 2011-11-21 | 2012-09-19 | 南开大学 | 一种涂覆增益介质的拉曼光谱高电磁增强基底及制备 |
CN103862032A (zh) * | 2014-02-26 | 2014-06-18 | 国家纳米科学中心 | 四方超晶格的核壳贵金属纳米棒及其自组装方法 |
CN106645083A (zh) * | 2016-11-18 | 2017-05-10 | 吉林大学 | 一种可变激发角的集成化等离子体增强拉曼光谱检测装置 |
KR101752838B1 (ko) * | 2016-03-03 | 2017-06-30 | 울산과학기술원 | 화학증강 라만산란용 기재, 이의 제조방법, 및 이를 이용한 화학증강 라만산란 방법 |
CN108460227A (zh) * | 2018-03-21 | 2018-08-28 | 长沙大衡智能科技有限公司 | 一种柔轮结构设计与优化方法 |
JP2020173185A (ja) * | 2019-04-11 | 2020-10-22 | Toyo Tire株式会社 | 構造物のfem解析方法、システム及びプログラム |
CN113096749A (zh) * | 2021-06-10 | 2021-07-09 | 武汉大学深圳研究院 | n型共掺杂金刚石半导体材料制备的多尺度耦合仿真方法 |
CN115859548A (zh) * | 2022-01-04 | 2023-03-28 | 天津大学 | 一种用于对防水密封胶渗透特性进行评估的方法及系统 |
CN117288735A (zh) * | 2022-06-16 | 2023-12-26 | 华北理工大学 | 一种核壳式sers基底的优化制备 |
CN117350137A (zh) * | 2023-12-05 | 2024-01-05 | 山东理工大学 | 一种放电等离子体通道瞬态特性的有限元仿真方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10710005B2 (en) * | 2017-07-26 | 2020-07-14 | Agency For Science, Technology And Research | Adsorbent material |
-
2024
- 2024-02-20 CN CN202410185682.1A patent/CN117744451B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002071013A1 (en) * | 2001-03-01 | 2002-09-12 | New Mexico State University Technology Transfer Corporation | Optical devices and methods employing nanoparticles, microcavities, and semicontinuous metal films |
CN102026719A (zh) * | 2008-03-11 | 2011-04-20 | 免疫之光有限责任公司 | 从外源进行内部能量活化的等离子体激元辅助系统和方法 |
CN101487832A (zh) * | 2009-02-10 | 2009-07-22 | 哈尔滨工业大学 | 一种表征环氧树脂复合材料界面机理的方法 |
CN102680453A (zh) * | 2011-11-21 | 2012-09-19 | 南开大学 | 一种涂覆增益介质的拉曼光谱高电磁增强基底及制备 |
CN103862032A (zh) * | 2014-02-26 | 2014-06-18 | 国家纳米科学中心 | 四方超晶格的核壳贵金属纳米棒及其自组装方法 |
KR101752838B1 (ko) * | 2016-03-03 | 2017-06-30 | 울산과학기술원 | 화학증강 라만산란용 기재, 이의 제조방법, 및 이를 이용한 화학증강 라만산란 방법 |
CN106645083A (zh) * | 2016-11-18 | 2017-05-10 | 吉林大学 | 一种可变激发角的集成化等离子体增强拉曼光谱检测装置 |
CN108460227A (zh) * | 2018-03-21 | 2018-08-28 | 长沙大衡智能科技有限公司 | 一种柔轮结构设计与优化方法 |
JP2020173185A (ja) * | 2019-04-11 | 2020-10-22 | Toyo Tire株式会社 | 構造物のfem解析方法、システム及びプログラム |
CN113096749A (zh) * | 2021-06-10 | 2021-07-09 | 武汉大学深圳研究院 | n型共掺杂金刚石半导体材料制备的多尺度耦合仿真方法 |
CN115859548A (zh) * | 2022-01-04 | 2023-03-28 | 天津大学 | 一种用于对防水密封胶渗透特性进行评估的方法及系统 |
CN117288735A (zh) * | 2022-06-16 | 2023-12-26 | 华北理工大学 | 一种核壳式sers基底的优化制备 |
CN117350137A (zh) * | 2023-12-05 | 2024-01-05 | 山东理工大学 | 一种放电等离子体通道瞬态特性的有限元仿真方法 |
Non-Patent Citations (3)
Title |
---|
Analog bandwidth optimization in photonic sampling system based on cavity-less optical pulse source;Zhengkai Li,;IEEE;20211231;全文 * |
氧化石墨烯吸附亚甲基蓝的分子动力学模拟;马莹;王恬;张恒;;高等学校化学学报;20191210(12);全文 * |
马莹 ; 王恬 ; 张恒 ; .氧化石墨烯吸附亚甲基蓝的分子动力学模拟.高等学校化学学报.2019,(12),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN117744451A (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sparre et al. | The physics of multiphase gas flows: fragmentation of a radiatively cooling gas cloud in a hot wind | |
Jandhyala et al. | A combined steepest descent-fast multipole algorithm for the fast analysis of three-dimensional scattering by rough surfaces | |
Godeferd et al. | Structure and dynamics of rotating turbulence: a review of recent experimental and numerical results | |
Vikram et al. | A novel wideband FMM for fast integral equation solution of multiscale problems in electromagnetics | |
Chourou et al. | HipGISAXS: a high-performance computing code for simulating grazing-incidence X-ray scattering data | |
CN103353675B (zh) | 基于微结构的移频超分辨显微成像方法和装置 | |
Wang et al. | Constructing shapes and mixing structures of black carbon particles with applications to optical calculations | |
CN104567682B (zh) | 液态环境下微粒三维位置纳米级分辨力测量方法 | |
Krucker-Velasquez et al. | Underscreening and hidden ion structures in large scale simulations of concentrated electrolytes | |
Litvinenko et al. | Computation of electromagnetic fields scattered from objects with uncertain shapes using multilevel Monte Carlo method | |
Kaminsky et al. | An efficient batch K-fold cross-validation voronoi adaptive sampling technique for global surrogate modeling | |
Banerjee et al. | Modelling nearest neighbour distributions of biased tracers using hybrid effective field theory | |
CN117744451B (zh) | 基于有限元和分子动力学的sers基底仿真优化方法 | |
CN111899329A (zh) | 全各向异性媒质的仿真方法、装置、计算机设备及存储介质 | |
Bista et al. | Using GPUs for realtime prediction of optical forces on microsphere ensembles | |
Hammack et al. | Experiments on deep-water waves with two-dimensional surface patterns | |
Fang et al. | Mechanics of bio-sediment transport | |
CN113625379B (zh) | 一种菲涅尔波带片的设计方法、制作方法和设计装置 | |
Fiala et al. | Using a multiscale toroidal element to model a hydrogen atom | |
Patro et al. | Speeding up particle trajectory simulations under moving force fields using graphic processing units | |
Huang et al. | Parallel Performance and Optimization of the Lattice Boltzmann Method Software Palabos Using CUDA | |
Bilbao et al. | Numerical analysis of electromagnetic fields | |
Fushimi et al. | A digital twin approach for experimental acoustic hologram optimization | |
Polak et al. | Fourth-order accurate numerical modeling of the multi-fluid plasma equations with adaptive mesh refinement | |
Hore | Predicting the Optical and Electrical Properties of Polymer Nanocomposites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |