Contents
function [E_ftfxfy_o,E_ftfxfy_e,obj]=fOPA1_1(obj_input,E_ftfxfy_OPA1pump)
_________________create a new instance of class ch23d_________________________
obj = obj_input('tag','OPA1_1');
_________________computational parameters_____________________________________
obj( 'Nz',200, 'propagationLength', 0.012,...
'abortIntensity', inf,...
'showResultEach',25, 'showProgress',1);
_________________material parameters__________________________________________
obj( 'material','LBO','theta',1.5708,'phi',0.186227,'plane','XY','temperature',25);
_________________effects to be included_________________________________________
obj( 'deff1',8.4e-13,'eoo',1,'oeo',0,'ooo',0,...
'deff2',0,'oee',0,'eeo',0,'eee',0,...
'n2',4e-20,'SPMOn',0,'selfFocusingOn',0,...
'walkoffON',1,'noiseOn',0,'noiseFactor',1);
____________________________input pulse definition______________________________
obj({ 'field1','EorI','unchanged','polarization','e',...
'GD',0,'GDD',0,'TOD',0,'Phase',0,...
'shift_x',0,'shift_y',0,...
'alpha_x',0,'alpha_y',0,...
'slant_x',0,'slant_y',0,...
'radiusOfCurvature_x',0,'radiusOfCurvature_y',0});
obj({ 'pulse2','EorI',1e-09,'polarization','o','centerWavelength',8.1e-07,...
'pulseDuration',3.45e-14,'GD',0,'GDD',5e-26,'TOD',0,'Phase',0,'beamShape_t','gauss',...
'beamRadius_x',0.00015,'beamRadius_y',0.00015,...
'shift_x',0,'shift_y',0,...
'alpha_x',0,'alpha_y',0,...
'slant_x',0,'slant_y',0,...
'radiusOfCurvature_x',0,'radiusOfCurvature_y',0,...
'beamShape_x',{'supergauss',2},'beamShape_y',{'supergauss',2}});
____________________________detector properties________________________________
obj.detectors = rmfield(obj.detectors,fieldnames(obj.detectors));
obj.detectors.obs1.plotIntegratedProfiles=1;
obj.detectors.obs1.logScale=0;
obj.detectors.obs1.polarization='o';
obj.detectors.obs1.ftlim=[3.3e+14,4e+14];
obj.detectors.obs1.fxlim=[-32000,31000];
obj.detectors.obs1.fylim=[-inf,inf];
obj.detectors.obs1.spectralData=1;
obj.detectors.obs1.spacetimeData=1;
obj.detectors.obs1.peakpower=1;
obj.detectors.obs2.plotIntegratedProfiles=1;
obj.detectors.obs2.logScale=0;
obj.detectors.obs2.polarization='e';
obj.detectors.obs2.ftlim=[4e+14,7.224e+14];
obj.detectors.obs2.fxlim=[-15000,15000];
obj.detectors.obs2.fylim=[-inf,inf];
obj.detectors.obs2.spectralData=1;
obj.detectors.obs2.spacetimeData=1;
obj.detectors.obs2.peakpower=1;
____________________________run simulation__________________________________
[E_ftfxfy_o,E_ftfxfy_e]=run(obj,{'field1','Eftfxfy',E_ftfxfy_OPA1pump},{'pulse2'});
##### SOURCE BEGIN #####
function [E_ftfxfy_o,E_ftfxfy_e,obj]=fOPA1_1(obj_input,E_ftfxfy_OPA1pump)
% automaticly generated matlab script to simulate a long pulse collinear OPA with chi3D (C) T.Lang (tinolang@chi23d.com)
% obj_input is an instance of chi3D including the same sampling
% (Nt,Nx,Ny,tWindow,xWindow,yWindow)as it was used to generate the input
% pump pulse E_ftfxfy_OPA1pump (see example "SHG with subsequent collinear
% OPA (long pump pulse)")
%% ___________________create a new instance of class ch23d_________________________
obj = obj_input('tag','OPA1_1');
%% ___________________computational parameters_____________________________________
obj( 'Nz',200, 'propagationLength', 0.012,...
'abortIntensity', inf,...
'showResultEach',25, 'showProgress',1);
%% ___________________material parameters__________________________________________
obj( 'material','LBO','theta',1.5708,'phi',0.186227,'plane','XY','temperature',25);
%% ___________________effects to be included_________________________________________
obj( 'deff1',8.4e-13,'eoo',1,'oeo',0,'ooo',0,...
'deff2',0,'oee',0,'eeo',0,'eee',0,...
'n2',4e-20,'SPMOn',0,'selfFocusingOn',0,...
'walkoffON',1,'noiseOn',0,'noiseFactor',1);
%% ______________________________input pulse definition______________________________
obj({ 'field1','EorI','unchanged','polarization','e',...
'GD',0,'GDD',0,'TOD',0,'Phase',0,...
'shift_x',0,'shift_y',0,...
'alpha_x',0,'alpha_y',0,...
'slant_x',0,'slant_y',0,...
'radiusOfCurvature_x',0,'radiusOfCurvature_y',0});
obj({ 'pulse2','EorI',1e-09,'polarization','o','centerWavelength',8.1e-07,...
'pulseDuration',3.45e-14,'GD',0,'GDD',5e-26,'TOD',0,'Phase',0,'beamShape_t','gauss',...
'beamRadius_x',0.00015,'beamRadius_y',0.00015,... %please choose a reasonable values in y
'shift_x',0,'shift_y',0,...
'alpha_x',0,'alpha_y',0,...
'slant_x',0,'slant_y',0,...
'radiusOfCurvature_x',0,'radiusOfCurvature_y',0,...
'beamShape_x',{'supergauss',2},'beamShape_y',{'supergauss',2}});
%% ______________________________detector properties________________________________
obj.detectors = rmfield(obj.detectors,fieldnames(obj.detectors)); %delete default detectors
obj.detectors.obs1.plotIntegratedProfiles=1;
obj.detectors.obs1.logScale=0;
obj.detectors.obs1.polarization='o';
obj.detectors.obs1.ftlim=[3.3e+14,4e+14];
obj.detectors.obs1.fxlim=[-32000,31000];
obj.detectors.obs1.fylim=[-inf,inf];
obj.detectors.obs1.spectralData=1;
obj.detectors.obs1.spacetimeData=1;
obj.detectors.obs1.peakpower=1;
obj.detectors.obs2.plotIntegratedProfiles=1;
obj.detectors.obs2.logScale=0;
obj.detectors.obs2.polarization='e';
obj.detectors.obs2.ftlim=[4e+14,7.224e+14];
obj.detectors.obs2.fxlim=[-15000,15000];
obj.detectors.obs2.fylim=[-inf,inf];
obj.detectors.obs2.spectralData=1;
obj.detectors.obs2.spacetimeData=1;
obj.detectors.obs2.peakpower=1;
%% ______________________________run simulation__________________________________
[E_ftfxfy_o,E_ftfxfy_e]=run(obj,{'field1','Eftfxfy',E_ftfxfy_OPA1pump},{'pulse2'});
% Results=obj.simResults
% fprintf('pulse energy obs2 (J): %.2e \n',Results.obs1.PulseEnergy(end));
##### SOURCE END #####
Tino Lang
Locked
