sas3d_example2 - Jupyter Notebook

Menu Python 3 Trusted File Edit View Insert Cell Kernel WidgetsHelp In [1]: import numpy as npimport escape as escfrom escape.utils import show. . .3D structures, potential object and small angle scattering (experimental). Example 2.¶ In this notebook we describe two other entites for sample description: potential and geometry, which are used for creation of structured layered samples. Morphology analysis and metrology of structured samples is of great interest nowadays for electronics industry as well as for scinetific applications. The potential module is under development at the moment and some features are still buggy and do not function properly. Please do not use it for real applications. We continue here with a more complicated example: two trilayer gratings located at a distance 𝐷D from each other. We start with a definition of a structure unit. In [2]: #define widths of the initial boxwx=esc.par("Width X", 50, userlim=[10, 100], units="nm")wy=esc.par("Width Y", 200, userlim=[10, 500], units="nm")wz=esc.par("Width Z", 50, userlim=[10, 100], units="nm")D = esc.par("Distance", 170, userlim=[50, 200], units="nm")T = esc.par("Period", 315, userlim=[50, 400], units="nm")qx=esc.var("Qx")qy=esc.var("Qy")qz=esc.var("Qz")mat1=esc.amorphous("", 1e-5, 0, color="#FFFF00")mat2=esc.amorphous("", 3e-5, 0, color="#00FF00")mat3=esc.amorphous("", 1e-5, 0, color="#FFFF00")V1=esc.box("", qx, qy, qz, mat1, wx, wy, wz, pos=[ (-D/4.0, 0, wz/2.0), (D/4.0, 0, wz/2.0), ])wz2=esc.par("Width Z2", 80, userlim=[10, 160], units="nm")wx2=esc.par("Width X2", 50, userlim=[10, 100], units="nm")V2=esc.box("", qx, qy, qz, mat2, wx2, wy, wz2, pos=[ (-D/4.0, 0, wz+wz2/2.0), (D/4.0, 0, wz+wz2/2.0),])wz3=esc.par("Width Z3", 80, userlim=[10, 160], units="nm")wx3=esc.par("Width X3", 50, userlim=[10, 100], units="nm")V3=esc.box("", qx, qy, qz, mat3, wx3, wy, wz3, pos=[ (-D/4.0, 0, wz+wz2+wz3/2.0), (D/4.0, 0, wz+wz2+wz3/2.0),])V=V1+V2+V3show(V, title="Result shape")Result shapeWidth X[nm]:50.00Width Y[nm]:200.0Width Z[nm]:50.00Distance[nm]:170.0Width X2[nm]:50.00Width Z2[nm]:80.00Width X3[nm]:50.00Width Z3[nm]:80.00 INFO: Refreshed picking FBO to size (337, 168)INFO: Found 0 X3D and nodes...INFO: X3DOM 1.8.2-dev, Build: 7180, Revison: , Date: Tue Jul 21 09:40:22 2020 -0400INFO: Time for setup and init of GL element no. 0: 173 ms.INFO: activate FogBindable null/defaultX3DFogNodeINFO: register FogBindable null/defaultX3DFogNodeINFO: create new Fog for X3DFogNode-stackINFO: activate ViewpointBindable null/Oblique viewINFO: activate first X3DViewpointNode for X3DViewpointNode-stackINFO: activate BackgroundBindable null/INFO: activate first X3DBackgroundNode for X3DBackgroundNode-stackINFO: System ready.INFO: activate NavigationInfoBindable null/INFO: activate first X3DNavigationInfoNode for X3DNavigationInfoNode-stackINFO: activate EnvironmentBindable null/defaultX3DEnvironmentNodeINFO: register EnvironmentBindable null/defaultX3DEnvironmentNodeINFO: create new Environment for X3DEnvironmentNode-stackWARNING: Fallback to inefficient multi-index mode since creaseAngle=0.INFO: register ViewpointBindable null/Oblique viewINFO: register NavigationInfoBindable null/INFO: NavType: undefinedINFO: Switch to examine mode.INFO: register BackgroundBindable null/INFO: webgl2 context found Vendor: Mozilla Google Inc., Renderer: Mozilla ANGLE (Intel(R) UHD Graphics 620 Direct3D11 vs_5_0 ps_5_0), Version: WebGL 2.0, ShadingLangV.: WebGL GLSL ES 3.00, MSAA samples: 4 Extensions: EXT_color_buffer_float, EXT_float_blend, EXT_texture_compression_bptc, EXT_texture_filter_anisotropic, OES_texture_float_linear, OVR_multiview2, WEBGL_compressed_texture_s3tc, WEBGL_compressed_texture_s3tc_srgb, WEBGL_debug_renderer_info, WEBGL_debug_shaders, WEBGL_lose_contextINFO: Initializing X3DCanvas for [x3dom-1595347069442-canvas]INFO: Creating canvas for (X)3D element...INFO: Found 1 X3D and nodes...INFO: X3DOM 1.8.2-dev, Build: 7180, Revison: , Date: Tue Jul 21 09:40:22 2020 -04000 . . .. . .Calculating Small Angle Scattering intensity.¶ For the calculation of small angle scattering intensity map, we should create a grid of coordinates (𝑄𝑥Qx, 𝑄𝑦Qy, 𝑄𝑧Qz) first. In [3]: #Wavelength Cu-Kalphawl = 1.5405 #Size of a gridsz = 512#magnitude of the wave vectork0 = 2.0 * np.pi / 1.5405#angle arrays for the calculation of scttering vector componentsa = np.linspace(-3, 3, sz)b = np.linspace(-3, 3, sz)#scattering vector meshgridav, bv = np.meshgrid(a, b)#caluclating scattering vector component for every point_qy = k0 * np.sin(bv * np.pi / 180.0)_qx = k0 * np.cos(bv * np.pi / 180.0)*np.sin(av * np.pi / 180.0)_qz = k0 * np.cos(bv * np.pi / 180.0)*np.cos(av * np.pi / 180.0)-k0#creating coordinates array of type (x0, y0, z0, x1, y1, z1, ...)c=np.column_stack((_qx.flatten(), _qy.flatten(), _qz.flatten())).flatten(). . .In [4]: #this lattice has a limited size of 10 gratingsT=esc.par("Period X", 100, units="nm")latx = esc.lattice("", qx, T, 10). . .We create Small Angle Cross section kernel as following In [5]: I = esc.sas("", qx, qy, qz, V, latx, multithread=False)show(I, coordinates=c, cblog=True, cbmin=1, rows=sz, plot_type="map")Width X[nm]:50.00Width Y[nm]:200.0Width Z[nm]:50.00Distance[nm]:170.0Width X2[nm]:50.00Width Z2[nm]:80.00Width X3[nm]:50.00Width Z3[nm]:80.00Period X[nm]:100.0Figure . . .In [ ]: . . .