defines optical properties of specified materials See also OPT-PROD, Chap. 12}. ForSDUM= WV-LIMIT: defines wavelength range for optical photon transportWHAT(1)> 0.0: minimum wavelength (in cm) for optical photon transport = 0.0: ignored < 0.0: resets to defaultDefault: 2.50E-5 (250 nm)WHAT(2)> 0.0: central wavelength (in cm) for optical photon transport = 0.0: ignored < 0.0: resets to defaultDefault: 5.89E-5 ((589 nm, Na D line)WHAT(3)> 0.0: maximum wavelength (in cm) for optical photon transport = 0.0: ignored < 0.0: resets to defaultDefault: 6.00E-5 (600 nm) WHAT(4-6): assignment to materials, see belowSDUM= WV-LIMIT ForSDUM= OM-LIMIT: defines angular frequency range for optical photon transportWHAT(1)> 0.0: minimum angular frequency for optical photon transport omega = 2pi x frequency (in rad/s) = 0.0: ignored < 0.0: resets to defaultDefault: 3.14E15 rad/s (corresponding to 600 nm)WHAT(2)> 0.0: central angular frequency for optical photon transport omega = 2pi x frequency (in rad/s) = 0.0: ignored < 0.0: resets to defaultDefault: 3.20E15 rad/s (corresponding to 589 nm, Na D line)WHAT(3)> 0.0: maximum angular frequency for optical photon transport omega = 2pi x frequency (in rad/s) = 0.0: ignored < 0.0: resets to defaultDefault: 7.53E15 rad/s (corresponding to 250 nm) WHAT(4-6): assignment to materials, see belowSDUM= OM-LIMIT ForSDUM= RESET: all optical properties are zeroed WHAT(1-3) = no meaning WHAT(4-6): assignment to materials, see belowSDUM= RESET ForSDUM= METAL: flag the material as a metalWHAT(1)= 1^{st}optical property (not used at the moment)WHAT(2)= 2^{nd}optical property (not used at the moment)WHAT(3)= 3^{rd}optical property: (1-r), where r is the reflectivity index at the central wavelength (or at the central angular frequency, depending on which one of the two quantities has been defined). See also Note 2Default= 0.0 WHAT(4-6): assignment to materials, see belowSDUM= METAL ForSDUM= blank:WHAT(1)= 1^{st}optical property: refraction index at the central wavelength (or at the central angular frequency, depending on which one of the two quantities has been defined) < -99: forces to use user routine RFRNDX (see Note 1)Default= 1.0WHAT(2)= 2^{nd}optical property: absorption coefficient (cm^-1) at the central wavelength (or at the central angular frequency, depending on which one of the two quantities has been defined) < -99: forces to use user routine ABSCFF (see Note 1)Default= 0.0WHAT(3)= 3^{rd}optical property: diffusion coefficient (cm^-1) at the central wavelength (or at the central angular frequency, depending on which one of the two quantities has been defined) < -99: forces to use user routine DFFCFF (see Note 1)Default= 0.0 WHAT(4-6): assignment to materials, see belowSDUM= blank ForSDUMcontaining &1 (resp. &2) (resp. &3):WHAT(1)= 4^{th}(resp. 7th) (resp. 10th) optical property of the material (derivatives of the refraction index, see Note 2 below)Default= 0.0WHAT(2)= 5th (resp. 8th) (resp. 11th) optical property of the material (derivatives of the absorption coefficient, see Note 2 below)Default= 0.0WHAT(3)= 6th (resp. 9th) (resp. 12th) optical property of the material (derivatives of the diffusion coefficient, see Note 2 below)Default= 0.0 WHAT(4-6): assignment to materials, see belowSDUM= &1, &2 or &3 For all previousSDUM's:WHAT(4)= lower bound of the indices of materials to which the indicated optical properties refer From materialWHAT(4)...Default= 3WHAT(5)= upper bound of the indices of materials to which the indicated optical properties refer ...To materialWHAT(5)...Default=WHAT(4)WHAT(6)= ... in step ofWHAT(6)Default= 1 ForSDUM= SENSITIV: sets up the optical photon detection sensitivity parameters (See alsoSDUM= WV-SENSI,SDUM= OM-SENSI, Note 3 below and the examples in Chap. 12})WHAT(1)= 0th photon sensitivity parameter < -99: forces to use user routine QUEFFC (see Note 1)WHAT(2)= 1^{st}optical property (not used at the moment)WHAT(3)= 2^{nd}optical property (not used at the moment)WHAT(4)= 3^{rd}optical property: (1-r), where r is the reflectivity index at the central wavelength (or at the central angular frequency, depending on which one of the two quantities has been defined). See also Note 2Default= 0.0WHAT(5)= maximum optical photon sensitivity over the allowed range (must be consistent with the previous values). It can be overestimated.Default= 1.0WHAT(6): not usedSDUM= SENSITIV ForSDUM= WV-SENSI: setup the wavelength of the optical photon sensitivity (See alsoSDUM= SENSITIV,SDUM= OM-SENSI and Note 3 below)WHAT(1)> 0.0: minimum wavelength (in cm) for optical photon sensitivity = 0.0: ignored < 0.0: resets to defaultDefault= 2.5 x 1.0E-5 (250 nm)WHAT(2)> 0.0: central wavelength (in cm) for optical photon sensitivity = 0.0: ignored < 0.0: resets to defaultDefault= 5.89 x 1.0E-5 (589 nm, Na D line)WHAT(3)> 0.0: maximum wavelength (in cm) for optical photon sensitivity = 0.0: ignored < 0.0: resets to defaultDefault= 6.0 x 1.0E-5 (600 nm)WHAT(4)-WHAT(6): not usedSDUM= WV-SENSI ForSDUM= OM-SENSI: setup the angular frequency of the optical photon sensitivity (See alsoSDUM= SENSITIV,SDUM= WV-SENSI and Note 3 below)WHAT(1)> 0.0: minimum angular frequency for optical photon sensitivity omega = 2 pi x frequency in rad/s = 0.0: ignored < 0.0: resets to defaultDefault= 3.14E15 rad/s (corresponding to 600 nm)WHAT(2)> 0.0: central angular frequency for optical photon sensitivity omega = 2 pi x frequency in rad/s = 0.0: ignored < 0.0: resets to defaultDefault= 3.20E15 rad/s (corresponding to 589 nm, Na D line)WHAT(3)> 0.0: maximum angular frequency for optical photon sensitivity omega = 2 pi x frequency in rad/sDefault= 7.53E15 rad/s (corresponding to 250 nm)WHAT(4)-WHAT(6): not usedSDUM= OM-SENSI ForSDUM= SPEC-BDX: flags special boundary crossings for optical photons At the selected boundary crossings special user-defined properties are defined by means of the user routine OPHBDX (see Note 1 below). A maximum of 40 boundaries can be flagged by issuing option OPT-PROP withSDUM= SPEC-BDX as many times as needed.WHAT(1)>= 1.0: special boundary treatment activated for the n-th+1 boundary = 0.0: ignored =< -1.0: special boundary treatment deactivated for the n-th+1 boundaryWHAT(2)= One of the two regions defining the n-th+1 boundaryWHAT(3)= The other region defining the n-th+1 boundaryWHAT(4)>= 1.0: special boundary treatment activated for the n-th+2 boundary = 0.0: ignored =< -1.0: special boundary treatment deactivated for the n-th+2 boundaryWHAT(5)= One of the two regions defining the n-th+2 boundaryWHAT(6)= The other region defining the n-th+2 boundarySDUM= SPEC-BDXDefault(option OPT-PROP not given): no optical photon transportNotes:1) The optional user routines concerning optical photons are: * RFRNDX: to specify a refraction index as a function of wavelength, frequency or energy. Activated by settingWHAT(1)< -99 whenSDUM= blank. * ABSCFF: to specify an absorption coefficient as a function of wavelength, frequency or energy. This is activated by settingWHAT(2)< -99 whenSDUM= blank. * DFFCFF: to specify a diffusion coefficient as a function of wavelength, frequency or energy. Activated by settingWHAT(3)< -99 whenSDUM= blank. * RFLCTV: to specify the reflectivity of a material. This can be activated withSDUM= METAL andWHAT(3)< -99 * OPHBDX: to set optical properties of a boundary surface. The call is activated withSDUM= SPEC-BDX * FRGHNS: to set a possible degree of surface roughness, in order to have both diffusive and specular reflectivity from a given material surface (NOT YET IMPLEMENTED) * QUEFFC: to introduce Quantum efficiency. This is activated by setting the 0-th photon sensitivity parameter to a value < -99 withSDUM= SENSITIV . 2) The 9 material properties input by the user withSDUM= &1, &2 and &3 are derivatives of order 1 to 3 of the three basic quantities which are input withSDUM= blank orSDUM= METAL. The three basic quantities and their derivatives are used to perform three series expansions centreed on the selected central wavelength (ifSDUM= WV-LIMIT has been used for the material concerned) or on the selected central angular frequency (ifSDUM= OM-LIMIT has been used). The variable x used in the expansion is adimensional: x = ( lambda - lambda_central ) / lambda_central or x = ( omega - omega_central ) / omega_central 3) The three basic quantities are: IfSDUM= blank: 1) refraction index n_refr 2) absorption coefficient mu_abs 3) diffusion coefficient mu_diff IfSDUM= METAL: 1) (not implemented yet) 2) (not implemented yet), 3) (1-r) (r = reflectivity index) The 4^{th}, 5th, 6th material properties are: d n_refr/ dx, d mu_abs/dx, d mu_diff/dx the 7th, 8th, 9th material properties are: d2 n_refr/ dx^2, d2 mu_abs/dx^2, d2 mu_diff/dx^2 the 10th, 11th, 12th material properties are: d3 n_refr/ dx^3, d3 mu_abs/dx^3, d3 mu_diff/dx^3 4)SDUM= SENSITIV can be used to set the quantum efficiency as a function of photon energy OVERALL through the problem and it is not material/region dependent. The reason is that it is applied "a priori" at photon generation time. If a quantum efficiency curve is to be introduced AT DETECTION, then a weighting user routine should be used, such as FLUSCW or COMSCW. It is advantageous to reduce computer time avoiding to generate and transport a photon which would have a significant probability to remain undetected. 5) The optical photon sensitivity parameters are: epsilon(0), d epsilon/dx, d2 epsilon/dx2, d3 epsilon/dx3 where x is: x = (lambda - lambda_central)/lambda_central or x = (omega - omega_central)/omega_centralExample 1:* Optical photon transport requested between 3.E15 and 7.E15 rad/s* (4.77E5 and 1.11E6 GHz, or 314 to 628 nm) for materials 6,9,12,15 and 18*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...OPT-PROP 3.E15 6.E15 7.E15 6.0 18.0 3. OM-LIMIT* User routine USRMED called when an optical photon is going to be transported* in materials 6, 12 and 18MAT-PROP 1.0 0.0 0.0 6.0 18.0 6. USERDIREExample 2:* Material 11 has a reflectivity index = 0.32*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...OPT-PROP 0.0 0.0 0.32 11.0 0.0 0. METALExample 3:**...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...* Optical photon transport requested between 300 and 600 nm for water.* (material 9). The optical properties are for the Na D line (589 nm)MATERIAL 1.0 0.0 .0000899 3.0 0.0 1. HYDROGEN MATERIAL 8.0 0.0 0.00143 5.0 0.0 0. OXYGEN MATERIAL 0.0 0.0 1.0 21.0 0.0 0. WATER COMPOUND 2.0 3.0 1.0 5.0 0.0 0. WATER OPT-PROP 3.E-5 5.89E-5 6.E-5 21.0 0.0 0. WV-LIMIT* diffusion coefficient of water from Appl.Opt. 17, 3587 (1978)OPT-PROP 1.33299 0.0013 1.22E-5 21.0 0.0 0.