CN102680503A - Method for determining secondary electron emission characteristics of rough metal surfaces - Google Patents

Method for determining secondary electron emission characteristics of rough metal surfaces Download PDF

Info

Publication number
CN102680503A
CN102680503A CN2012101309564A CN201210130956A CN102680503A CN 102680503 A CN102680503 A CN 102680503A CN 2012101309564 A CN2012101309564 A CN 2012101309564A CN 201210130956 A CN201210130956 A CN 201210130956A CN 102680503 A CN102680503 A CN 102680503A
Authority
CN
China
Prior art keywords
electronics
electron
point
energy
motion
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101309564A
Other languages
Chinese (zh)
Other versions
CN102680503B (en
Inventor
张娜
曹猛
张海波
崔万照
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201210130956.4A priority Critical patent/CN102680503B/en
Publication of CN102680503A publication Critical patent/CN102680503A/en
Application granted granted Critical
Publication of CN102680503B publication Critical patent/CN102680503B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The invention relates to a method for determining the secondary electron emission characteristics of rough metal surfaces. The method is provided for judging whether electrons enter a material, and is utilizing for judging whether the electrons enter again after getting out of the material, consequently, the multiple interactions between the electrons and a rough metal surface are taken into full consideration, and by means of statistical analysis on a large quantity of electrons, the secondary electron emission characteristic of the metal material is obtained. According to the method, since the multiple interactions between the electrons getting away frsom the surface of the material and the rough surface are taken into full consideration, the obtained secondary electron emission characteristic is matched with an experimental result, and moreover, the method can be adopted to analyze the secondary electron emission characteristics of any surface morphologies, and can be used for researching the law between surface morphologies and secondary electron emission characteristics so as to provide theoretical direction for searching surface morphologies capable of inhibiting secondary electron emission.

Description

A kind of method of definite roughened metal surface secondary electron emission characteristic
Technical field
The invention discloses a kind of method of definite roughened metal surface secondary electron emission characteristic, belong to the Physical Electronics field.
Background technology
The interaction of electron beam and material is a very important field in the Physical Electronics research, and the secondary that is wherein excited by incident electron (being also referred to as secondary electron emissions) phenomenon has important use at modern electronic equipments such as scanning electron microscope, Auger electron spectrometer and electron multiplier and device.For example: in electron micrology and surface electron spectrum analysis, adopt the beam bombardment sample; The scattering of in material, passing through electronics produces the various characteristic signals of exosyndrome material character, thus various microscopic properties such as the crystal structure of acquisition material, constituent, electronic structure, surface topography, inherent vice.
Secondary not only relates to the scattering process of electronics at material internal, and electronics also can interact with material surface behind the outgoing material.The physical process of electronics after the material surface outgoing is remarkable to the influence of secondary electron emission characteristic, and this also is the experimental data of many secondary electron emission characteristics and the major reason that notional result can't be coincide.The method of confirming the secondary electron emission characteristic of material at present is divided into two big types: the one, go out analytic formula approximate description secondary electron emission characteristic through experimental fit; The 2nd, through the emission process of DSMC simulation secondary electron, obtain the secondary electron emission characteristic of material.The existing scheme overwhelming majority only is applicable to level and smooth material surface, has promptly only considered the physical process of electronics at material internal.And for real material, the surface is not desirable smooth flat, and it has certain roughness, and promptly the pattern on surface rises and falls.In individual other disposal route, adopt only as the secondary electron emission characteristic of probability method statistics electronics, ignored secondary electron when outgoing with the repeatedly interaction process on surface.
Summary of the invention
The objective of the invention is to overcome the above-mentioned deficiency of prior art; A kind of method of definite roughened metal surface secondary electron emission characteristic is provided; When this method has been considered electronics outgoing material surface and the repeatedly interaction of rough surface, the secondary electron emission characteristic and the experimental result of acquisition are more identical.
Above-mentioned purpose of the present invention mainly is achieved through following technical scheme:
A kind of method of definite roughened metal surface secondary electron emission characteristic comprises the steps:
(1) rough surface morphology to material characterizes, if material surface pattern rule is divided into rectangular node with surface topography, in the rectangular node arbitrarily the fluctuating of some height h adopt analytic expression to represent, and confirm the fluctuating height h that the surface is maximum MaxIf material surface is the irregular surfaces pattern; Adopting microscope to carry out surface topography extracts; Adopt the method for rectangular node point to sample; And set up two Dimension Numerical and describe surface topography, obtain in the rectangular node point of surface topography maximal value h of all height h that rise and fall of the fluctuating of some height h and surface arbitrarily Max
(2) with electronics from A I-1Point moves to A iThe length of the straight-line trajectory of point is designated as t i, with electronics from A iPoint moves to A I+1The length of the straight-line trajectory of point is designated as t I+1, t then I+1Obtain by following recurrence relation, then:
t i + 1 = t i + min ( t x , t y ) n x ≠ 0 , n y ≠ 0 t i + t x n x ≠ 0 , n y = 0 t i + t y n x = 0 , n y ≠ 0 | h i - z i | n x = 0 , n y = 0
t x = ( a - x i ) / n x n x > 0 - x i / n x n x < 0
t y = ( b - y i ) / n y n y > 0 - y i / n y n y < 0
Wherein:
When electronics is in reference position or electronics free path starting point, make t 1=0; I=1,2,3 ...;
(x i, y i) be A iThe local coordinate of spot projection in the corresponding rectangle grid, z iBe A iSpot projection is perpendicular to the coordinate on the z axle of rectangular node;
(x I+1, y I+1) be A I+1The local coordinate of spot projection in the corresponding rectangle grid, z I+1Be A I+1Spot projection is perpendicular to the coordinate on the z axle of rectangular node;
n xAnd n yBe respectively electron motion direction unit vector x and y to projection;
h iBe local coordinate (x i, y i) corresponding surface undulation height;
(3) whether electronics with metal material surface interaction being taken place judges order
f=(z i-h i)(z i+1-h i+1)
If f>=0 and z I+1>h Max, showing that electronics can not interact with material, the energy and the direction of record electronics this moment get into step (8);
If f>0 and z I+1≤h Max, show that electronics does not interact with metal material surface, makes t i=t I+1, return step (2);
If f≤0 shows that electronics is from A iPoint moves to A I+1Can pass material surface during point; If z i-h i>0, show that electronics is from the outside material internal that gets into of material, gets into step (4); If z i-h i<0, show electronics from the material internal outgoing, jump to step (6);
Wherein:
h I+1Be local coordinate (x I+1, y I+1) corresponding surface undulation height;
(4) if electronics from the outside incident material internal of material, according to the incident direction of incident electron
Figure BSA00000709748800031
Surface normal with the incidence point material
Figure BSA00000709748800032
By energy conservation and the momentum conservation that electron motion satisfies, calculate the energy E of electron impact behind the material iAnd direction of motion
Figure BSA00000709748800033
E i=E p+U 0
n i &RightArrow; = E p E p + U 0 n p &RightArrow; + ( E p E p + U 0 cos &theta; - E p cos 2 &theta; + U 0 E p + U 0 ) N i &RightArrow;
Wherein:
E pBe the incident electron energy;
U 0Be built-in potential, U 0=E F+ Φ, E FBe the Fermi level of material, Φ is the work function of material;
Figure BSA00000709748800035
is the incident direction of incident electron;
Figure BSA00000709748800036
is the surface normal of incidence point material;
θ is the angle of surface normal of incident direction and the incidence point material of incident electron,
(5) according to the electron energy E that incides material internal iAnd direction of motion
Figure BSA00000709748800038
Adopt the scattering process of DSMC simulation electronic in material of conventional ideal surfaced, follow the trail of the track of incident electron and secondary electron, write down the electron energy E of each scattering j, the electron motion direction Free path S with electronics j, j=1 wherein, 2 ..., m, to electronics at every section free path S jInterior rectilinear motion is handled according to step 2 and step 3;
(6) when electron motion arrived material surface, recorded electronic was from the position (x of material surface outgoing m, y m, z m), according to the electronic motion direction Normal direction with the eye point material
Figure BSA00000709748800043
Calculate electronics exit probability T:
Figure BSA00000709748800044
Wherein:
Direction when
Figure BSA00000709748800045
arrives material surface for electron motion;
Figure BSA00000709748800046
is the normal direction of eye point material;
Figure BSA00000709748800047
is the angle of electron motion direction and eye point material normal direction,
Figure BSA00000709748800048
E mEnergy during for electronics outgoing material surface;
(7) go up generation random number R ND in interval [0,1],
If T≤RND then makes electronics reflect on the surface, reflection back electronics is still with energy E mContinuation is scattering in material, and direction of motion does
Figure BSA00000709748800049
Figure BSA000007097488000410
Jump to step (5) then;
If T>RND then makes the electronics outgoing, the energy of electronics is E after the outgoing s, direction of motion does
Figure BSA000007097488000411
E s=E m-U 0
Figure BSA000007097488000412
Return step (2)
(8) change the incident electron energy E p, repeating step (2)~(7), the secondary electron direction of motion and the energy of record outgoing are added up the outgoing electron number under the different incident electron energy, obtain the secondary electron yield of said material; The direction and the energy of statistics outgoing electron, the secondary electron power spectrum and the exit direction that obtain said material distribute.
In the method for above-mentioned definite roughened metal surface secondary electron emission characteristic, in the step (1) if material surface is the irregular surfaces pattern, then in the rectangular node point of surface topography arbitrarily the fluctuating of some height h adopt the method acquisition of first-order linear interpolation, that is:
h = h lu + h ru - h lu a m + h ld - h lu b n + h rd + h lu - h ru - h ld ab mn
Wherein:
(m n) is any point coordinate in the rectangular node under the local coordinate;
H is any fluctuating of some height in the rectangular node point;
h LuFluctuating height for the rectangle upper left corner;
h LdFluctuating height for the rectangle lower left corner;
h RuFluctuating height for the rectangle upper right corner;
h RdFluctuating height for the rectangle lower right corner;
A is the length of rectangular node point;
B is the wide of rectangular node point.
The present invention compared with prior art has following beneficial effect:
(1) the present invention has taken into full account the repeatedly interaction of the rough surface of electronics and metal reality, and the secondary electron emission characteristic and the experimental result of acquisition are more identical;
(2) novelty of the present invention provided judge electronics whether with the method for material surface interactions; Utilize this method can judge that electronics is after the material internal outgoing; Incident whether once more; Thereby considered behind the electronics outgoing material surface the repeatedly interaction with rough surface, the secondary electron emission characteristic of the metal material that obtains to coincide more with experimental result;
(3) the present invention has used the interaction of describing electronics and roughened metal surface near the method for electronics actual motion process, has disclosed the microscopic nature of electronics and Rough Metal Surface effect;
(4) the present invention can analyze the secondary electron emission characteristic of arbitrary surfaces pattern, can be used for studying the rule of surface topography and secondary electron emission characteristic, for exploring the surface topography that can suppress secondary theoretical direction is provided.
Description of drawings
Fig. 1 is the local coordinate system of irregular surfaces pattern rectangular node point of the present invention;
Fig. 2 is the synoptic diagram of cutting apart of trajectory of electron motion of the present invention;
Fig. 3 is the energy changing figure of electron impact of the present invention to material internal;
Fig. 4 is the energy changing figure of electronics of the present invention from the material internal outgoing;
Fig. 5 is the structural representation of rectangular channel in the embodiment of the invention;
Fig. 6 is the local coordinate system of the rectangular node point of rectangular channel in the embodiment of the invention;
Fig. 7 is electronic motion stepping in the embodiment of the invention;
Secondary electron yield that Fig. 8 obtains for the inventive method and experimental result, bibliographical information result's contrast synoptic diagram.
Embodiment
Below in conjunction with accompanying drawing and specific embodiment the present invention is described in further detail:
The present invention confirms the flow process of the method for roughened metal surface secondary electron emission characteristic, mainly comprises following step:
Step 1,
The metal material surface pattern is characterized, if regular morphology is divided into rectangular node with surface topography, in the rectangular node arbitrarily the fluctuating of some height h adopt analytic expression to represent, and confirm the fluctuating height h that the surface is maximum Max
If non-regular morphology adopts atomic force microscope or laser microscope to carry out surface topography and characterizes, adopt the method for rectangular node point to sample, and set up two Dimension Numerical and describe surface topography.Each element of array is corresponding to a point on the surface topography, the height h of its value to putting.The method acquisition of the height h of some employing first-order linear interpolation arbitrarily in the rectangular node, as shown in the formula:
h = h lu + h ru - h lu a m + h ld - h lu b n + h rd + h lu - h ru - h ld ab mn - - - ( 1 )
h max=max(h) (2)
h MaxThe maximal value of all height h that rise and fall for the surface;
Be illustrated in figure 1 as the local coordinate system of irregular surfaces pattern rectangular node point of the present invention, true origin is positioned at the rectangle upper left corner; (m n) is any point coordinate in the rectangular node under the local coordinate; H is any fluctuating of some height in the rectangular node point; h LuFluctuating height for the rectangle upper left corner; h LdFluctuating height for the rectangle lower left corner; h RuFluctuating height for the rectangle upper right corner; h RdFluctuating height for the rectangle lower right corner; A is the length of rectangular node point; B is the wide of rectangular node point.
Step 2,
If electronics outside material, is divided into some segments with its straight-line trajectory; If electronics at material internal, is divided into some segments with the straight path between the double scattering of electronics, with electronics from A I-1Point moves to A iThe length of the straight-line trajectory of point is designated as t i, with electronics from A iPoint moves to A I+1The length of the straight-line trajectory of point is designated as t I+1, t then I+1Obtain by following recurrence relation:
t i + 1 = t i + min ( t x , t y ) n x &NotEqual; 0 , n y &NotEqual; 0 t i + t x n x &NotEqual; 0 , n y = 0 t i + t y n x = 0 , n y &NotEqual; 0 | h i - z i | n x = 0 , n y = 0
t x = ( a - x i ) / n x n x > 0 - x i / n x n x < 0 - - - ( 3 )
t y = ( b - y i ) / n y n y > 0 - y i / n y n y < 0
Wherein:
When electronics is in reference position or electronics free path starting point, make t 1=0, i=1,2,3 ...;
(x i, y i) be A iThe local coordinate of spot projection in the corresponding rectangle grid, z iBe A iSpot projection is perpendicular to the coordinate on the z axle of rectangular node; (x I+1, y I+1) be A I+1The local coordinate of spot projection in the corresponding rectangle grid, z I+1Be A I+1Spot projection is perpendicular to the coordinate on the z axle of rectangular node; n xAnd n yBe respectively electron motion direction unit vector x and y to projection; h iBe local coordinate A i(x i, y i), the surface undulation height that adopts step 1 to calculate.
Be illustrated in figure 2 as the synoptic diagram of cutting apart of trajectory of electron motion of the present invention.
Step 3,
Whether electronics with metal material surface interaction is taken place judge order
f=(z i-h i)(z i+1-h i+1) (4)
If f>=0 and z I+1>h Max, showing that electronics can not interact with material, the energy and the direction of record electronics this moment get into step 8;
If f>0 and z I+1≤h Max, show that electronics does not interact with metal material surface, makes t i=t I+1, return step 2;
If f≤0 shows that electronics is from A iPoint moves to A I+1Can pass material surface during point; If z i-h i>0, show that electronics is from the outside material internal that gets into of material, gets into step 4; If z i-h i<0, show that electronics from the material internal outgoing, jumps to step 6.
Wherein:
h I+1Be local coordinate A I+1(x I+1, y I+1), the surface undulation height that adopts step (1) to calculate.
Step 4,
If electronics is from the outside incident material internal of material, according to the incident direction of incident electron
Figure BSA00000709748800081
Surface normal with the incidence point material
Figure BSA00000709748800082
By energy conservation and the momentum conservation that electron motion satisfies, calculate the energy E of electron impact behind the material iAnd direction of motion
Figure BSA00000709748800083
E i=E p+U 0 (5)
n i &RightArrow; = E p E p + U 0 n p &RightArrow; + ( E p E p + U 0 cos &theta; - E p cos 2 &theta; + U 0 E p + U 0 ) N i &RightArrow; - - - ( 6 )
Wherein:
E pBe incident electron energy (with respect to vacuum level);
U 0Be built-in potential, U 0=E F+ Φ, E FBe the Fermi level of material, Φ is the work function of material;
E iBe the energy (with respect to conduction band bottom) of electron impact to material internal;
Figure BSA00000709748800085
is the incident direction of incident electron;
Figure BSA00000709748800086
is the surface normal of incidence point material;
θ is the angle of surface normal of incident direction and the incidence point material of incident electron,
Figure BSA00000709748800087
Be illustrated in figure 3 as the energy changing figure of electron impact of the present invention to material internal.
Step 5,
According to the electron energy E that incides material internal iAnd direction Adopt the scattering process of DSMC simulation electronic in material of conventional ideal surfaced, follow the trail of the track of incident electron and secondary electron, write down the electron energy E of each scattering j, the electron motion direction
Figure BSA00000709748800092
Free path S with electronics j(j=1,2 ..., m), to electronics at every section free path S jInterior rectilinear motion is handled according to step 2 and step 3; E wherein m,
Figure BSA00000709748800093
S mEnergy, direction of motion and free path when being respectively electronics from the material surface outgoing.
Step 6,
When electron motion arrived material surface, recorded electronic was from material surface eye point A mPosition (x m, y m, z m), according to the electronic motion direction
Figure BSA00000709748800094
Normal direction with the eye point material
Figure BSA00000709748800095
Calculate exit probability T:
Figure BSA00000709748800096
Wherein:
Direction when
Figure BSA00000709748800097
arrives material surface for electron motion;
Figure BSA00000709748800098
is the normal direction of eye point material;
Figure BSA00000709748800099
is the angle of electronics exit direction and eye point material normal direction,
Figure BSA000007097488000910
E mEnergy during for electronics outgoing material surface.
Step 7,
Go up generation random number R ND in interval [0,1], if T≤RND then makes electronics reflect on the surface, reflection back electronics is still with energy E mContinuation is scattering in material.And direction of motion is
Figure BSA000007097488000911
Figure BSA000007097488000912
Jump to step 5 then.
If T>RND then makes the electronics outgoing, the energy of electronics is E after the outgoing s, direction of motion does
Figure BSA000007097488000913
E s=E m-U 0 (9)
Jump to step 2 then.
Wherein:
E sBe the energy (with respect to vacuum level) behind the electronics outgoing material;
Figure BSA00000709748800101
is the direction of motion behind the electronics outgoing material.
Be illustrated in figure 4 as the energy changing figure of electronics of the present invention from the material internal outgoing.
Step 8,
Change the incident electron energy E pRepeating step two~step 7; Record can outgoing secondary electron direction of motion and energy; Add up the secondary electron yield that outgoing electron number under the different projectile energies can obtain this metal material, the direction of statistics outgoing electron and secondary electron power spectrum that energy can obtain this metal material and exit direction distribution.
Embodiment
Impinge perpendicularly on the Au material surface with electronics, the Au material surface is the rectangular channel structure, and the parameter of rectangular channel is degree of depth 500nm, the wide 180nm of bottom land, the wide 25.4nm of cell wall, cycle 205.4nm.Fig. 5 has provided the structural representation of rectangular channel.
The rectangular channel structure is a regular morphology, can adopt analytic expression to represent.Divide the rectangular node point according to rectangular node local coordinate shown in Figure 6, the long limit a of rectangular node point is 205.4nm, and broadside b is an endless.The surface undulation of arbitrfary point is in the rectangular node:
h = 0 0 &le; x < 180 nm 500 180 nm &le; x < 205.4 nm
Suppose that electronics initially is positioned at rectangular channel central authorities and vertical incidence surface, A 1Position coordinates be (90, y 1, 500) and (because rectangular channel does not have limit for width, y 1Can be arbitrary value), it is at (the x that is projected as of rectangular node point 1, y 1)=(90, y 1), A 1Spot projection is perpendicular to the coordinate z on the rectangular node z axle 1=500, A 1The fluctuating height h of point 1=0.According to step 2, t 1=0, n xAnd n yBe respectively electron motion direction unit vector x and y to projection because vertical incidence, n x=0 and n y=0, next movement length t of electronics then 2=| h 1-z 1|=500, next movement position A of electronics 2For (90, y 1, 0), A 2Spot projection is perpendicular to the coordinate z on the rectangular node z axle 2=0, A 2The fluctuating height h of point 2=0.Be illustrated in figure 7 as electronic motion stepping in the embodiment of the invention.
According to step 3, f=0, z 1-h 1>0, show that electronics enters into material internal from the material outside.
According to step 4, the Fermi level E of gold copper-base alloy F=5.1eV, work function Φ=5.1eV, built-in potential U 0=10.2eV.The electronics vertical incidence; Angle theta=0 of the surface normal
Figure BSA00000709748800112
of incident direction
Figure BSA00000709748800111
and incidence point material; With incident electron energy 1000eV is example, energy E i and the direction of motion
Figure BSA00000709748800113
of electron impact behind the material
E i=E p+U 0=1010.2(eV)
n i &RightArrow; = E p E p + U 0 n p &RightArrow; + ( E p E p + U 0 cos &theta; - E p cos 2 &theta; + U 0 E p + U 0 ) N i &RightArrow;
= ( 2 E p E p + U 0 - 1 ) n p &RightArrow;
Step 5 adopts the scattering process of DSMC simulation electronic in material of conventional ideal surfaced, writes down the electron energy E of each scattering j, the electron motion direction
Figure BSA00000709748800116
J=1 wherein; 2; ...; M, " the Monte Carlo Modeling for Electron Microscopy and Microanalysis " that concrete grammar can be shown referring to Chinese University of Science and Technology's doctorate paper " scanning electron microscopy learn in the Monte Carlo simulation that produces of secondary electron " and Divid C.Joy.
Suppose through step 5 the energy E of electronics mBe 400eV, the electronic motion direction does
Figure BSA00000709748800117
Figure BSA00000709748800118
Surface normal with the eye point material
Figure BSA00000709748800119
Angle be 80 °,
Figure BSA000007097488001110
X axle and y axle be projected as 0.5 and 0.848, the local coordinate of eye point position be (200, y m, 500) because z m=h m, judge f=0 according to step 3, and electronics has z when material internal moves M-1<h M-1, jump to step 6.
According to step 6, because the probability T of
Figure BSA000007097488001111
electronics outgoing:
Figure BSA000007097488001112
According to step 7, go up generation random number R ND=0.6 in interval [0,1], since T>RND, then electronics outgoing, and the energy of electronics is E after the outgoing s, direction of motion does
Figure BSA000007097488001113
E s=E m-U 0=389.8(eV)
Figure BSA00000709748800121
Figure BSA00000709748800122
The exit direction of electronics and eye point normal direction angle are 68.6 °.
Return step 2, A mThe point coordinate be (200, y m, 500), the exit direction unit vector at x and y to being projected as 0.473 and 0.802.t x=(205.4-200)/and 0.473=11.42, because b is infinitely great, obvious t xMinimum, then t M+1=t m+ t x=11.42, A M+1Position coordinates be (205.4, y m+ 9.16,504.16), A M+1Spot projection is perpendicular to the coordinate z on the rectangular node z axle M+1=504.16, A M+1The fluctuating height h of point M+1=500.
According to step 3, f=0 and z M+1>h M+1, electronics can not interact with material, the energy E of record electronics this moment sAnd direction
Figure BSA00000709748800123
Get into step 8;
The step that provides is to specifications calculated, and makes projectile energy E p=100eV, 200eV, 300eV; ..., 1900eV, 2000eV; 2500eV, 3000eV, repeating step two~step 7; Add up the outgoing electron number under the different projectile energies, can obtain the secondary electron yield of this metal material, the ratio of the secondary electron yield of this secondary electron yield and ideal plane is relative secondary electron yield; Its result is as shown in Figure 8, and Fig. 8 has provided secondary electron yield that the inventive method obtains and experiment, bibliographical information result's contrast synoptic diagram, but by Fig. 8 knowledge capital inventive method because when having considered electronics outgoing material surface and the repeatedly interaction of rough surface; The secondary electron emission characteristic and the experimental result that obtain are more identical, and the inventive method can be confirmed the secondary electron emission characteristic of metal material real surface.
The above; Be merely the best embodiment of the present invention, but protection scope of the present invention is not limited thereto, any technician who is familiar with the present technique field is in the technical scope that the present invention discloses; The variation that can expect easily or replacement all should be encompassed within protection scope of the present invention.
The content of not doing to describe in detail in the instructions of the present invention belongs to this area professional and technical personnel's known technology.

Claims (2)

1. the method for a definite roughened metal surface secondary electron emission characteristic is characterized in that comprising the steps:
(1) rough surface morphology to material characterizes, if material surface pattern rule is divided into rectangular node with surface topography, in the rectangular node arbitrarily the fluctuating of some height h adopt analytic expression to represent, and confirm the fluctuating height h that the surface is maximum MaxIf material surface is the irregular surfaces pattern; Adopting microscope to carry out surface topography extracts; Adopt the method for rectangular node point to sample; And set up two Dimension Numerical and describe surface topography, obtain in the rectangular node point of surface topography maximal value h of all height h that rise and fall of the fluctuating of some height h and surface arbitrarily Max
(2) with electronics from A I-1Point moves to A iThe length of the straight-line trajectory of point is designated as t i, with electronics from A iPoint moves to A I+1The length of the straight-line trajectory of point is designated as t I+1, t then I+1Obtain by following recurrence relation, then:
t i + 1 = t i + min ( t x , t y ) n x &NotEqual; 0 , n y &NotEqual; 0 t i + t x n x &NotEqual; 0 , n y = 0 t i + t y n x = 0 , n y &NotEqual; 0 | h i - z i | n x = 0 , n y = 0
t x = ( a - x i ) / n x n x > 0 - x i / n x n x < 0
t y = ( b - y i ) / n y n y > 0 - y i / n y n y < 0
Wherein:
When electronics is in reference position or electronics free path starting point, make t 1=0; I=1,2,3 ...; (x i, y i) be A iThe local coordinate of spot projection in the corresponding rectangle grid, z iBe A iSpot projection is perpendicular to the coordinate on the z axle of rectangular node;
(x I+1, y I+1) be A I+1The local coordinate of spot projection in the corresponding rectangle grid, z I+1Be A I+1Spot projection is perpendicular to the coordinate on the z axle of rectangular node;
n xAnd n yBe respectively electron motion direction unit vector x and y to projection;
h iBe local coordinate (x i, y i) corresponding surface undulation height;
(3) whether electronics with metal material surface interaction being taken place judges order
f=(z i-h i)(z i+1-h i+1)
If f>=0 and z I+1>h Max, showing that electronics can not interact with material, the energy and the direction of record electronics this moment get into step (8);
If f>0 and z I+1≤h Max, show that electronics does not interact with metal material surface, makes t i=t I+1, return step (2);
If f≤0 shows that electronics is from A iPoint moves to A I+1Can pass material surface during point; If z i-h i>0, show that electronics is from the outside material internal that gets into of material, gets into step (4); If z i-h i<0, show electronics from the material internal outgoing, jump to step (6);
Wherein:
h I+1Be local coordinate (x I+1, y I+1) corresponding surface undulation height;
(4) if electronics from the outside incident material internal of material, according to the incident direction of incident electron
Figure FSA00000709748700021
Surface normal with the incidence point material
Figure FSA00000709748700022
By energy conservation and the momentum conservation that electron motion satisfies, calculate the energy E of electron impact behind the material iAnd direction of motion
Figure FSA00000709748700023
E i=E p+U 0
n i &RightArrow; = E p E p + U 0 n p &RightArrow; + ( E p E p + U 0 cos &theta; - E p cos 2 &theta; + U 0 E p + U 0 ) N i &RightArrow;
Wherein:
E pBe the incident electron energy;
U 0Be built-in potential, U 0=E F+ Φ, E FBe the Fermi level of material, Φ is the work function of material;
is the incident direction of incident electron;
is the surface normal of incidence point material;
θ is the angle of surface normal of incident direction and the incidence point material of incident electron,
Figure FSA00000709748700027
(5) according to the electron energy E that incides material internal iAnd direction of motion
Figure FSA00000709748700031
Adopt the scattering process of DSMC simulation electronic in material of conventional ideal surfaced, follow the trail of the track of incident electron and secondary electron, write down the electron energy E of each scattering j, the electron motion direction
Figure FSA00000709748700032
Free path S with electronics j, j=1 wherein, 2 ..., m, to electronics at every section free path S jInterior rectilinear motion is handled according to step 2 and step 3;
(6) when electron motion arrived material surface, recorded electronic was from the position (x of material surface outgoing m, y m, z m), according to the electronic motion direction
Figure FSA00000709748700033
Normal direction with the eye point material Calculate electronics exit probability T:
Figure FSA00000709748700035
Wherein:
Direction when
Figure FSA00000709748700036
arrives material surface for electron motion;
Figure FSA00000709748700037
is the normal direction of eye point material;
Figure FSA00000709748700038
is the angle of electron motion direction and eye point material normal direction,
Figure FSA00000709748700039
E mEnergy during for electronics outgoing material surface;
(7) go up generation random number R ND in interval [0,1],
If T≤RND then makes electronics reflect on the surface, reflection back electronics is still with energy E mContinuation is scattering in material, and direction of motion does
Figure FSA000007097487000310
Jump to step (5) then;
If T>RND then makes the electronics outgoing, the energy of electronics is E after the outgoing s, direction of motion does
Figure FSA000007097487000312
E s=E m-U 0
Figure FSA000007097487000313
Return step (2)
(8) change the incident electron energy E p, repeating step (2)~(7), the secondary electron direction of motion and the energy of record outgoing are added up the outgoing electron number under the different incident electron energy, obtain the secondary electron yield of said material; The direction and the energy of statistics outgoing electron, the secondary electron power spectrum and the exit direction that obtain said material distribute.
2. the method for a kind of definite roughened metal surface secondary electron emission characteristic according to claim 1; It is characterized in that: in the said step (1) if material surface be the irregular surfaces pattern; The then method acquisition of the fluctuating of some height h employing first-order linear interpolation arbitrarily in the rectangular node point of surface topography, that is:
h = h lu + h ru - h lu a m + h ld - h lu b n + h rd + h lu - h ru - h ld ab mn
Wherein:
(m n) is any point coordinate in the rectangular node under the local coordinate;
H is any fluctuating of some height in the rectangular node point;
h LuFluctuating height for the rectangle upper left corner;
h LdFluctuating height for the rectangle lower left corner;
h RuFluctuating height for the rectangle upper right corner;
h RdFluctuating height for the rectangle lower right corner;
A is the length of rectangular node point;
B is the wide of rectangular node point.
CN201210130956.4A 2012-04-27 2012-04-27 Method for determining secondary electron emission characteristics of rough metal surfaces Active CN102680503B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210130956.4A CN102680503B (en) 2012-04-27 2012-04-27 Method for determining secondary electron emission characteristics of rough metal surfaces

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210130956.4A CN102680503B (en) 2012-04-27 2012-04-27 Method for determining secondary electron emission characteristics of rough metal surfaces

Publications (2)

Publication Number Publication Date
CN102680503A true CN102680503A (en) 2012-09-19
CN102680503B CN102680503B (en) 2014-07-02

Family

ID=46812761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210130956.4A Active CN102680503B (en) 2012-04-27 2012-04-27 Method for determining secondary electron emission characteristics of rough metal surfaces

Country Status (1)

Country Link
CN (1) CN102680503B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103196932A (en) * 2013-02-26 2013-07-10 西安空间无线电技术研究所 Determination method for secondary electron emission coefficient of metal surface of microwave component
CN103366377A (en) * 2013-07-25 2013-10-23 北京航空航天大学 Micro-contact characteristic and image threshold value processing method based on feature characteristics
CN104569014A (en) * 2014-10-27 2015-04-29 西安空间无线电技术研究所 Method and device for testing secondary electron emission coefficient of material under all incidence angles
CN111727489A (en) * 2017-12-15 2020-09-29 德累斯顿莱布尼茨固体材料研究所 Momentum-resolved photoelectron spectrometer and method for momentum-resolved photoelectron spectroscopy

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202066807U (en) * 2011-05-17 2011-12-07 华南理工大学 Electron collector being suitable for intelligently testing secondary electron emission coefficient of solid material

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202066807U (en) * 2011-05-17 2011-12-07 华南理工大学 Electron collector being suitable for intelligently testing secondary electron emission coefficient of solid material

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨春城 等: "离子在金属及介质表面感生的二次电子发射", 《中国空气动力学学会第十届物理气体动力学专利委员会会议论文集》 *
谢爱根 等: "高能原电子能量与金属的有效真二次电子发射系数的关系", 《强激光与粒子束》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103196932A (en) * 2013-02-26 2013-07-10 西安空间无线电技术研究所 Determination method for secondary electron emission coefficient of metal surface of microwave component
CN103196932B (en) * 2013-02-26 2014-11-19 西安空间无线电技术研究所 Determination method for secondary electron emission coefficient of metal surface of microwave component
CN103366377A (en) * 2013-07-25 2013-10-23 北京航空航天大学 Micro-contact characteristic and image threshold value processing method based on feature characteristics
CN103366377B (en) * 2013-07-25 2016-04-27 北京航空航天大学 A kind of micro-contact performance based on shape characteristic and image threshold disposal route
CN104569014A (en) * 2014-10-27 2015-04-29 西安空间无线电技术研究所 Method and device for testing secondary electron emission coefficient of material under all incidence angles
CN104569014B (en) * 2014-10-27 2017-01-25 西安空间无线电技术研究所 Method and device for testing secondary electron emission coefficient of material under all incidence angles
CN111727489A (en) * 2017-12-15 2020-09-29 德累斯顿莱布尼茨固体材料研究所 Momentum-resolved photoelectron spectrometer and method for momentum-resolved photoelectron spectroscopy
CN111727489B (en) * 2017-12-15 2023-11-17 德累斯顿莱布尼茨固体材料研究所 Momentum resolving photoelectron spectrometer and method for momentum resolving photoelectron spectroscopy

Also Published As

Publication number Publication date
CN102680503B (en) 2014-07-02

Similar Documents

Publication Publication Date Title
Belcher Mixing and transport in urban areas
CN102680503B (en) Method for determining secondary electron emission characteristics of rough metal surfaces
Berg et al. The Bolund experiment, part I: flow over a steep, three-dimensional hill
Draper et al. Energy potential of a tidal fence deployed near a coastal headland
CN102955152B (en) Synthetic aperture radar (SAR) signal simulation method for sea waves
Abe et al. Classification of filament formation mechanisms in magnetized molecular clouds
CN106951596B (en) SPH numerical simulation method for interaction of tidal bore and sea pond
Turpin et al. Effect of the topography of an industrial site on dust emissions from open storage yards
CN103196932A (en) Determination method for secondary electron emission coefficient of metal surface of microwave component
Banta et al. Rugged terrain effects on diffusion
Wang et al. POD characterisation of extreme wake patterns of turbulent wind fields past rectangular buildings
CN109723027B (en) Method for designing rain intensity scale Sr in flood discharge atomization model test
Wang et al. On the near-field of a square jet with vortex-generating tabs
Zhang et al. Development and application of a street-level meteorology and pollutant tracking system (S-TRACK)
Grašič et al. Re-evaluation of the Lagrangian particle modelling system on an experimental campaign in complex terrain
Bakhtyar et al. Rapid wave modeling of severe historical extratropical cyclones off the Northeastern United States
CN103106315A (en) Simulation method for calculating value of dust removal efficiency of windproof dust suppression net
Meroney Bluff-body aerodynamics influence on transport and diffusion
Tang et al. Wind-sand movement characteristics and erosion mechanism of a solar photovoltaic array in the middle of the Hobq Desert, Northwestern China
Balczó et al. Wind tunnel and computational fluid dynamics study of wind conditions in an urban square
Manganhar et al. ANALYSIS OF WIND FLOW DIVERTING AND ACCELERATING CHARACTERISTICS OF POROUS AND NON-POROUS CIRCULAR OBSTACLES
Guo-dong et al. Wind-sand movement characteristics and erosion mechanism of a solar photovoltaic array in the middle of the Hobq Desert, Northwestern China
Kellnerova et al. Flow instabilities within an urban intersection
Le Ribault et al. Large eddy simulation of particle transport and deposition over multiple 2D square obstacles in a turbulent boundary layer
Cleynen et al. Performance mapping of ducted freestream hydropower devices

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant