function [Physic, Material, Port1Solver, Port2Solver] = configPhysicMaterialFem4()
%CONFIGPHYSICMATERIALFEM4  Fem4 / MeshData2x 基准物理与材料参数
% 对应 COMSOL Fem4.mph：脊形波导，传播方向 z，长度 4 um（MeshData2x 实测 z in [0,4um]）
% 材料与 MainFem2.m 一致：SiO2 + Si + 掺杂 Si(σ=5000)

Physic.lam0 = 1.55e-6;
Physic.c_const = 299792458;
Physic.k0 = 2*pi/Physic.lam0;
Physic.omega = 2*pi*Physic.c_const/Physic.lam0;

Material.epsilonr = [4; 11.9; 11.9; 11.9];
Material.mur = [1; 1; 1; 1];
Material.sigma = [0; 0; 5000; 0];
Material.epsilon0 = 8.854187817e-12;
Material.temp = 1/Physic.k0*120*pi;

Material.epsilonrF = [0    0 4    4    0 0 0 11.9 11.9 0 ...
                      0    0 11.9 11.9 0 0 0 0    0    11.9 ...
                      11.9 0 0    0]';
Material.murF = ones(24, 1);
Material.sigmaF = [0 0 0 0 0 0 0 0 0 0 0 0 5000 5000 ...
                   0 0 0 0 0 0 0 0 0 0]';

Physic.Port1 = [3 8 13 20];
Physic.Port2 = [4 9 14 21];
Physic.PEC = [1 2 5 24];

Port1Solver.targetNeff = 2.6;
Port1Solver.modeNum = 1;
Port1Solver.lam0 = Physic.lam0;
Port1Solver.type = 1;

Port2Solver.targetNeff = 2.6;
Port2Solver.modeNum = 1;
Port2Solver.lam0 = Physic.lam0;
Port2Solver.type = 2;

end
