////////////////////////////////////////////////////////////////////////////////////////// // // input for lpic // ////////////////////////////////////////////////////////////////////////////////////////// &pulse_front ------------------------------------------------------------------------------------------ Q = 1 # switch ON (Q=1) OFF (Q=0) amplitude = 0 # dimensionless laser field amplitude amplitude2 = 2 # dimensionless laser field amplitude, 2nd harmonic amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic amplituden = 0.0 # dimensionless laser field amplitude, nth harmonic phase2 = 0 # 2nd harmonic's phase with respect to fundamental [degree] phase3 = 0 # 3rd harmonic's phase [degree] phasen = 0 # nth harmonic's phase [degree] delay2 = 0 # delay of 2nd harmonic envelope with respect to main envelope delay3 = 0 # delay of 3rd harmonic envelope with respect to main envelope delayn = 0 # delay of nth harmonic envelope with respect to main envelope seedfrequency = 14.0 # frequency of nth harmonic in fundamental frequencies angle = 45 # in degree polarization = 2 # s=1, p=2, c=3 shape = 3 # linear=1, sin=2, sin^2=3 raise = 15 # pulse raise/fall time in periods duration = 30 # pulse duration in periods delay = 0 # propagation delay in cycles pulse_save = 1 # save pulse shape? yes=1, no=0 pulse_save_step = 0.0004 # time step in periods &pulse_rear ------------------------------------------------------------------------------------------ Q = 0 # switch ON (Q=1) OFF (Q=0) amplitude = 0.0 # dimensionless laser field amplitude amplitude2 = 0.0 # dimensionless laser field amplitude, 2nd harmonic amplitude3 = 0.0 # dimensionless laser field amplitude, 3rd harmonic amplituden = 0.0 # dimensionless laser field amplitude, nth harmonic phase2 = 0 # 2nd harmonic's phase with respect to fundamental phase3 = 0 # 3rd harmonic's phase phasen = 0 # nth harmonic's phase [degree] delay2 = 0 # delay of 2nd harmonic envelope with respect to main envelope delay3 = 0 # delay of 3rd harmonic envelope with respect to main envelope delayn = 0 # delay of nth harmonic envelope with respect to main envelope seedfrequency = 14.0 # frequency of nth harmonic in fundamental frequencies angle = 0 # in degree polarization = 1 # s=1, p=2, c=3 shape = 3 # linear=1, sin=2, sin^2=3 raise = 10 # pulse raise/fall time in periods duration = 20 # pulse duration in periods delay = 0 # propagation delay in cycles pulse_save = 1 # save pulse shape? yes=1, no=0 pulse_save_step = 0.02 # time step in periods &propagate ------------------------------------------------------------------------------------------ prop_start = 0 # start time in periods prop_stop = 50 # stop time in periods &box ------------------------------------------------------------------------------------------ cells_per_wl = 1000 # cells per wavelength, lab frame cells = 8000 # total number of cells cells_left = 3500 # cells vacuum left cells_plasma = 2000 # occupied cells for cells_ramp=0 cells_rampgluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the righter exp-ramp n_ion_over_nc = 225 # maximum density/critical density ux0 = 0.0 # initial momentum in x-direction box_save = 1 # save configuration? yes=1, no=0 &lramp ------------------------------------------------------------------------------------------ length = 0 # linear ramp length (form=0) or righter scale length (form=1 or 2) form = 1 # 0 = linear; 1 = exponential; 2=doubleexponential cutoff = -50 # max ramplength (<0 <=> disabled) gluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the righter exp-ramp length2 = 0 # lefter scale length (form=2) &rramp ------------------------------------------------------------------------------------------ length = 0 # cells in the linear ramp region (form=0) or scalelength form = 0 # 0 = linear; 1 = exponential length2 = 100 # righter scale length (form=2) cutoff = -1 # max ramplength (<0 <=> disabled) gluepos = 250 # glue-point between the two exp-functions (form=2) ... length of the lefter exp-ramp &electrons ------------------------------------------------------------------------------------------ fix = 0 # 0->not fixed 1->fixed ppc = 200 # max. number of MacroParticles per cell vtherm = 0.000 # thermal velocity in units of C &ions ------------------------------------------------------------------------------------------ fix = 1 # 0->not fixed 1->fixed z = 1 # q/e m = 50000 # m/m_e ppc = 200 # max. number of MacroParticles per cell vtherm = 0.00 # thermal velocity in units of C &output ------------------------------------------------------------------------------------------ path = ../data_stufe18_0 # output path &energy Q = 1 # energy plot? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &flux Q = 1 # flux plot? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &reflex Q = 1 # reflectivity plot? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &snapshot Q = 1 # snapshots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &el_phasespace Q = 0 # phasespace plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 0.05 # time step in periods &ion_phasespace Q = 0 # phasespace plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 0.05 # time step in periods &el_velocity Q = 1 # electron velocity distributions? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &ion_velocity Q = 1 # ion velocity distributions? t_start = 0 # start time in periods t_stop = 50 # stop time in periods t_step = 1 # time step in periods &de Q = 1 # electron density plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &di Q = 1 # ion density plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jx Q = 0 # jx plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jy Q = 0 # jy plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &jz Q = 0 # jz plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ex Q = 0 # ex plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ey Q = 0 # ey plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &ez Q = 0 # ez plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &bx Q = 0 # bx plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &by Q = 0 # by plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &bz Q = 0 # bz plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &edens Q = 1 # field energy density plots? t_start = 0 # start time in periods t_stop = 50 # stop time in periods x_start = 0 # left boundary in cells x_stop = 100000 # right boundary in cells &traces Q = 1 # traces? t_start = 0 # start time in periods t_stop = 50 # stop time in periods traces = 111 # of traces at fixed positions x: t0=2, t1=2500, t2=2510, t3=2520, t4=2530, t5=2540, t6=2550, t7=2560, t8=2570, t9=2580, t10=2590, t11=2600, t12=2610, t13=2620, t14=2630, t15=2640, t16=2650, t17=2660, t18=2670, t19=2680, t20=2690, t21=2700, t22=2710, t23=2720, t24=2730, t25=2740, t26=2750, t27=2760, t28=2770, t29=2780, t30=2790, t31=2800, t32=2810, t33=2820, t34=2830, t35=2840, t36=2850, t37=2860, t38=2870, t39=2880, t40=2890, t41=2900, t42=2910, t43=2920, t44=2930, t45=2940, t46=2950, t47=2960, t48=2970, t49=2980, t50=2990, t51=3000, t52=3010, t53=3020, t54=3030, t55=3040, t56=3050, t57=3060, t58=3070, t59=3080, t60=3090, t61=3100, t62=3110, t63=3120, t64=3130, t65=3140, t66=3150, t67=3160, t68=3170, t69=3180, t70=3190, t71=3200, t72=3210, t73=3220, t74=3230, t75=3240, t76=3250, t77=3260, t78=3270, t79=3280, t80=3290, t81=3300, t82=3310, t83=3320, t84=3330, t85=3340, t86=3350, t87=3360, t88=3370, t89=3380, t90=3390, t91=3400, t92=3410,t93=3420, t94=3430, t95=3440, t96=3450, t97=3460; t98=3470, t99=3480, t100=3490, t101=3500, t102=3510, t103=3520, t104=3530, t105=3540, t106=3550, t107=3560, t108=3570, t109=3580, t110=3590 ////////////////////////////////////////////////////////////////////////////////////////// &restart ------------------------------------------------------------------------------------------ Q = 0 # restart from intermediate stage? file = restart # start file Q_save = 1 # save intermediate stages periodically? file_save = restart # save file ¶llel ------------------------------------------------------------------------------------------ N_domains = 1 # number of parallel processes Q_reo = 1 # periodic reorganizations? delta_reo = 1 # laser cycles between reo's ////////////////////////////////////////////////////////////////////////////////////////// thermal velocity v/c vs. thermal energy [eV] M ( vx^2 + vy^2 + vz^2 ) = 3 k T ------------------------------------------------------------------------------------------ 1e-4 0.005 3e-4 0.046 5e-4 0.128 1e-3 0.511 4e-3 8.2 1e-2 51 2e-2 204 3e-2 460 4e-2 818 4.42e-2 1000 1e-1 5110 choose: v_th/c = (lambda_d/dx) * (omega_p/omega) * (2pi dx/lambda) :=1 = sqrt(z_ion * n_ion_over_nc) = 2pi/cells_per_wl ////////////////////////////////////////////////////////////////////////////////////////// source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 p-pol transm. 1 minp-pol transm. 1 maxp-pol transm. 1 source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 source p-pol refl. 1 minp-pol refl. 1 maxp-pol refl. 1 source p-pol refl. 1