Analysis of Specific Problems in Laser Scanning Optical System Design
 Author: Joo Won Don
 Organization: Joo Won Don
 Publish: Current Optics and Photonics Volume 15, Issue1, p22~29, 25 March 2011

ABSTRACT
We analyze aberrations in an optical laser printer system in order to know how to determine an allowable nonuniformity of the movement of a light spot , how to determine allowed variation of spot sizes, and how to minimize the influence of these deviations on technological errors In this paper, the correction and the tolerance of distortion are analyzed by using the concept of zonal and global distortions. The tolerance of field curvature is also obtained from Gaussian beam properties. In order to reduce the change of the entrance pupil position and to make a more compact laser printer system the minimum size of the rotator is exactly derived from the geometry with the introduction of the shift angle of the input beam.

KEYWORD
Zonal distortion , Binary spot size , Polygon mirror , (080.2740) Geometric optical design , (220.3620) Lens system design , (080.4298)Nonimaging optics

I. INTRODUCTION
Recently much research and many inventions for the development of laser scanning optical systems have been carried out in order to improve their resolution in keeping with the increase in the use of computers in a modern society. The greater part of the results are published in patents[111]and in books,[1214] but these works need more general and theoretical analyses.
The goal of this paper is to theoretically and synthetically analyze the distortion, spot size and field curvature required for a laser printer system in order to get a smaller and more uniform spot size. Also, we show theoretically how to minimize of the radius of the rotator by adjustment of the entrance pupil position in order to make a more compact laser printer.
A general type of laser printer consists of a semiconductor laser as a source, a collimator lens, a rotating mirror (in general a polygonal mirror), and a socalled "
fθ ” lens. The laser beam, assumed to be Gaussian, passes first through a collimator, and the parallel beam that is formed by the collimator makes a spot on a recording plane after passing through the polygonal mirror and the “fθ ” lens. All elements are important, but the "fθ ” lens is the most complicated and main part of the laser scanning optical system.When we design the “
fθ ” lens, two important things are considered: (1) the movement of the light spots cast on a recording plane by a polygon with a constant angular velocity has to be proportional to a scan angle, i.e.y' =f?θ , wherey' is the image height,f is the focal length of the “fθ ” lens andθ is the scan angle. In order to solve this problem we will define the “fθ ” distortion and investigate its allowable and optimum values. (2) The sizes of these spots must be uniform on the recording plane. The spot size is primarily determined by field curvature in practice, so we will investigate the tolerance of field curvature with the analysis of binary spot size and the Gaussian beam property. These topics will addressed in detail, and we will analyze not only the properties of these aberrations but also their allowable values.Note that other aberrations such as spherical aberration and coma, etc., can be disregarded due to the small aperture of the beam passing through the “
fθ ” lens. The effect of these aberrations will be discussed in the section 4.When the polygonal mirror rotates, as is known, the entrance pupil position travels along the scan angle
θ in the laser printer system[12, 15]. (In general we put the entrance pupil position on the facet of the rotator.) The smaller the size of the rotator, the smaller the change of entrance pupil position and the more compact the laser printer system. In order to minimize the size of the rotator we will introduce the shifting angle of the input beam along the facet of the rotator. Thus, the minimum size of the rotator and the changing value of the entrance pupil position will be derived without approximation.In order to reduce the influence of errors introduced in the manufacture of the rotator on the image shift in the perpendicular section of the recording plane, i.e. in the sagittal section, we can choose an optical system in which the sagittal section is different from the meridional. However, in this paper we will limit our analysis to the meridional section.
II. ANALYSIS OF DISTORTION
An ordinary distortionfree lens has the image height which follows the rule
y _{0}'=f?tanθ , wherey _{0}' is the ideal image height. Therefore, as is known, a common fractional distortionGD _{t} is usually calculated from the following equation[15]:where
y ' is the actual image height formed by a chief ray.In the case of a perfect “
fθ ” lens, i.e.y '=f?θ , from Eq.(1) we have the following common fractional global distortion:where
ξ =θ/θ _{max} (0≤ξ ≤1) is a fractional field coordinate. The plot of this distortion versusξ is shown on Fig. 1 (in the case ofθ _{max}=45°).According to the general Eq.(1) the distortion of the “
fθ ” lens must have a large negative value (large barrel distortion). Therefore, during the design of the “fθ ” lens we have to introduce a large controlled distortion value along the curve shown in Fig. 1, but it is very inconvenient to do so.The alternative way is to use a lens design program which is able to calculate distortion by using
y _{0}'=f?θ rather thany _{0}'=f?tanθ . In this case the fractional distortion called the global “fθ ” distortion has the following form:When we use the global “
fθ ” distortion Eq.(3) instead of Eq.(1), the distortion of the “fθ ” lens system does not need to be corrected with a large negative value and simply must have the smallest value available. The global “fθ ”distortion will be called global distortion in this paper.Let us investigate the distortion of the “
fθ ” lens in detail. Many authors use the socalled calibrated focal length, instead of the paraxial focal length, to reduce maximum global distortion in the “fθ ” lens[12]. But in order to obtain the optimal form of the correction and the tolerance of distortion in the laser printer system, we must note that the distortion has to be considered as nonuniform velocity of the scanning spot, i.e. as nonuniformity of the sizes of the groups consisting of several spots on the recording plane. This indicates that we have to consider not the global distortion according to Eq.(3) but local or zonal distortion(ZD ), which is determined as change not iny' but in the differential ofy' as follows:where d
y' and dy0' are the differential heights ofy' andy0' corresponding to differential scan angledθ , respectively.Only the
ZD is important for the design and evaluation of the laser scanning system and the proper specification must be done for its value. The question is what value ofZD must be considered. The simple analysis shows that we must use not the maximum value ofZD but its peakvalley value (the difference between maximum and minimum value). If we useZD during the design of the “fθ ” lens, we only have to minimize the peakvalley value ofZD .However, the problem is thatZD is not convenient for optical designers, who commonly useGD according to Eq.(3).Also, all commercial computer programs for optical system design useGD . In fact, when we useGD , it is not enough to minimize its value; it is very important to have the optimum form of distortion curve. To solve this problem we will determine the relation betweenGD andZD . From Eqs.(3) and (4) it is easy to find the relation of these distortions as follows:In general, as it is enough to consider only the 3rd and 5th order distortion for the “
fθ ” lens, we can writeGD as follows:where
a andb are the coefficients of the 3rd and 5th order distortions,u =ξ ^{2}. From Eqs.(5) and (6)ZD is expressed as follows:With the factor
=a/b we can express Eqs.(6) and (7) as follows:Κ The factor
shows the relation between the 3rd and 5th order distortions, hence it determines the form of the distortion curve. Our task is to find the optimal value ofΚ , i.e. the optimal form of global distortion correction. To complete this task, we introduce the coefficient C=Κ ZD _{pv}/GD _{pv} in the ratio of the peakvalley value ofZD toGD . It is clear that the smaller the value of coefficientC , the better the form of the distortion correction. Table 1. shows the peakvalley values ofZD andGD along , and Fig. 2 shows the plot C versusΚ , all from Eqs.(8a) and (8b).Κ In Fig. 2 we can divide
into two regions with boundary atΚ =5/3. The coefficientΚ C shows a comparatively gradual slope and a small value in the range ≤5/3. A sharp slope and large value in the range,Κ >5/3, characterizes coefficientΚ C . Thus, the former is a good region for the correction ofZD and the latter is inferior. In particular, in the best case the coefficientC has a minimum value1.8 in the range 2< <5/3, but in the worst case C has a maximum valueΚ 9.8 at =1. For example, ifΚ GD _{pv} is 0.2% and =1, thenΚ C =9.8 andZD _{pv} is obtained1.96 %(ZD _{pv}=C?GD _{pv}). That is even worse than the case ofGD _{pv}=0.5 % and =5/3 ( (Κ ZD _{pv}=0.9 %). If =±infinity, we have only 3rd order distortion and in this case from Table 1, we seeΚ C =3.During the design of an “
fθ ” lens we have to know the valueC or for the correction ofΚ ZD . Its value can be obtained from the shape ofGD because it depends on the valueu _{m}, and this valueu _{m} is decided by the factor (Κ u _{m}= /Κ 2 ). The valueu _{m} along is given in Table 1 and Fig. 3 shows the shapes ofΚ GD andZD along various .Κ In order to obtain a satisfactory
ZD of the “fθ ” lens we suggest the shape ofGD must have the forms (a), (b)or (c), in the left parts of Fig. 3, and must avoid the forms (d), (e) and (f), in the right parts of Fig. 3.Ultimately, we have to consider not only the value of
GD , but also its shape for the correction ofZD during the design of the “fθ ” lens.Note that from our analysis the calibrated focal length that is used to reduce the maximum value of
GD [12] cannot reduce the peakvalley value ofZD .III. DEFINITION OF SPOT SIZE ON A RECORDING PLANE
From the theory of Gaussian beams[16] the normalized irradiance of a Gaussian beam is expressed by the following form in any plane:
where
w_{x} andw_{y} are the beam radii at which the irradiance isexp (2 ) of the maximum value inx andy coordinates respectively, i.e. the beam has86.5 % of total beam power in the beam radius (w _{x} ^{2}+w _{y} ^{2})^{1/2} from properties of Gaussian function. Let us focus on the meridional section only, i.e.x =0 .The analysis of spot size on the image plane depends on the kind of scanning system. If any system is used to scan an image for computer image processing, we have to analyze the aerial spot image and evaluate the spot size at the level of beam power
86.5 %. For more detailed analysis in this situation we have to use MTF of scanning system. Nevertheless, our consideration deals with the laser printer and laser processing scanning systems. Such systems, laser printer and laser cutting device and so on, are characterized as the nonlinear process of spot image detection. If the exposure is higher than a reference levelI _{e} (I >I _{e}), it is recognized as a white image (for the printer) or cut material(for the cutting device) but ifI <I_{e} , then it is black or fails to cut This image transformation may be called binary according to the following equation:where
I (x ,y ) is the aerial image andI' (x ,y ) is the detected(binary) image. (See Fig. 4).For such a process the use of MTF is not proper and the binary spot size has to be determined in a different manner from the aerial image.
The binary spot size for onaxis beam in meridional section can be expressed as follows on the recording plane:
where
δ _{0}' is the binary spot size of onaxis beam on the recording plane in the meridional section,v_{0}' is the angular beam divergence in the image space and can be considered as the exit numerical apertureA_{0}' of the “fθ ” lens,w_{0}' is the beam radius at which the irradiance isexp (2 ) of the maximum value on the recording plane,J_{y} is Lagrange invariant in the meridional section andk is a coefficient to depend onI_{e} . Here we assume the beam waist onaxis is located on the recording plane. For example, in the case ofI_{e} =I_{0}exp (1/2 ),thenk is0.5 , whereI_{0} is the maximum irradiance of the onaxis beam.In general, since every laser does not have perfect coherence, we must consider Lagrange invariance in deciding the spot size of the Gaussian beam. Therefore, Eq.(10) is a more general equation than the known equation
δ _{0}' =aF λ[13], wherea is the aperture shape factor,F is the F number and λ is the wavelength.The value
k which decides the binary spot size may be arbitrarily chosen, but we must select the optimum valuek . We suggestk =0.5 , and in this caseI _{e}=I_{0}exp (1/2 ) is put at the point where the Gaussian curve has a maximum slope. It indicates that the nonuniformity of the binary spot sizes is very small on the recording plane under an instability of the reference irradianceI_{e} . Also, the effective binary spot size is two times less than in the case ofI_{e} =I_{0}exp (2 ) ork =1 .In the case of offaxis, the beam radius
w' on the recording plane is larger thanw_{0}' in spite of the fact that total beam energy is conserved, so the maximum irradiance is smaller than it is in the case of the onaxis (see Fig. 4). On the other hand, since the reference irradianceI_{e} is fixed, the binary image size of offaxis differs from that of the onaxis. In order to know how the binary image size of offaxis changes, let us introduce two parametersη andγ :where
η andγ are the fractional changes of aerial spot radii and binary spot radii respectively, which are determined by an arbitrary reference irradianceI_{e} .w_{0}' andw' are the Gaussian beam radii of the onaxis and the offaxis respectively, that contain86.5% of beam power (aerial spot radii). Also w_{b0}' andw'_{b} are the binary spot radii of the onaxis and the offaxis beam with a defined power respectively, and w_{b0}' =k_{w0}' from Eq.(10). In fact, from the parametersη andγ we can see the nonuniformity of aerial and binary spot size during the movement of the scanning spot on the recording plane.From Eqs.(11) and (12) and with the relation between the reference level
I and the Gaussian curve of offaxis, i.e. exp(2e k ^{2})/w'_{0} =exp{2(w'_{b}/w' )^{2}}/w' , we obtainγ as follows:For a small
η , i.e.η <<1, for the case ofk =0.5 Eq.(13) is approximately expressed as follows:From the negative sign of Eq.(14) we see the binary spot size of the offaxis beam is always smaller than that of the onaxis (since
w_{b}' =w_{0}' (1 η ^{2})) and the change of the offaxis spot size is proportional to the square of a small valueη . For example, ifη =20% , thenγ =4% according to both Eq.(13) and Eq.(14). Ifη =40% , thenγ =20% from Eq.(13) andγ =16% according to simple Eq.(14). From the above example we can use the simple Eq.(14) as an approximation for the case ofη <40% . For a greaterη we can reduce the value γ by the choice of the proper valuek from Eq.(13). It is interesting to note from Eq.(14) the binary process of detection favors spot size uniformity better than other processes, and with the binary detection we have greater tolerance for aberration.In the result, when we choose
k =0.5 the change of the offaxis binary spot size has the small value under the change of the offaxis Gaussian beam envelope, and the change of the onaxis binary spot size Δδ_{0}' has a minimum value under an irradiance instability ΔI_{e} .In addition, the length
l of scanning line on the recording plane is expressed as follows:Also the exit numerical aperture
A_{0}' of the “fθ ” lens is connected with the focal length f of the “fθ ” lens and the incident parallel beam diameterD :From Eqs.(10), (15) and (16) we have the following relation:
The parameters of the right hand side of Eq.(17) (the length of scanning line, Lagrange invariant of a laser beam and spot size on a recording plane) are already given. The parameters of the left hand side (the beam diameter and the maximum scan angle) can be arbitrarily chosen under the condition of a constant preservation of their product. Lagrange invariant depends on the quality of the laser beam, and for the perfect beam we have the minimum λ/π. For example, in the case of
1200 dpi .δ_{0}' =21.2 ㎛ ,l =220 mm ,J_{y} =λ/π (λ=7800A ),k =0.5 andθ_{max} =45 °, from Eq.(17) we can easily ascertain a beam diameter ofD =3.28 mm ., but if we chooseJ_{y} =λ, thenD =10.3 mm . We can see the quality of laser beam (degree of coherence) is very important in a laser printer system because the smaller the beam diameter, the better the resolution of a printed image.IV. ANALYSIS OF MERIDIONAL FIELD CURVATURE AND PRELIMINARY “
fθ ” LENS DESIGNNow we consider the requirements for the correction of a meridional field curvature within an allowable nonuniformity of spot sizes on the recording plane. From the fact that the laser beam has Gaussian beam properties, the beam radius is decided by the following form both in meridional and sagittal section on the recording plane:[16]
where
w' is the beam radius on the recording plane,w_{0}' is the beam radius in the waist,z' is the distance from the recording plane to the beam waist, andz'_{0} =w_{0}' /v_{0}' =w_{0}'^{2} /J is the Rayleigh range. In this section we assume that the beam waist radius of the offaxis is equal to that of the onaxis since we want to see the effect of the variation of spot size due only to the meridional field curvature. For analysis in the meridional section, we usez'_{m} andJ_{y} instead of z' andJ in Eq.(18) respectively, wherez'_{m} is the meridional field curvature andJ_{y} is the meridional Lagrange invariant.The fractional nonuniformity of the aerial spot sizes due to the meridional field curvature is defined by the following expression from Eqs.(11) and (18):
From Eq.(14) and
δ_{0}' =w_{0}' (in the case ofk =0.5 ) we can determine a tolerance limit of a meridional field curvature:where
γ is an allowable fractional nonuniformity of binary spot sizes. For example, in the case ofγ =10 %,δ_{0}' =21.2㎛ (1200dpi ) andJ_{y} =λ/π (λ=7800A ), thenz'_{m} <1.55 mm but in the case ofJ_{y} =λ we havez'_{m} <0.49 mm . Eq.(20) as well as Eq.(14) are comparatively accurate in the case of smallη . However, ifη is larger than 0.4, we cannot get a simple expression and the problem may be solved by numerical analysis. We limit our analysis only to small values ofη .In fact, we have to use
w'/cosθ' instead ofw' in Eq.(19) since the direction of the offaxis outgoing beam is not consistent with the normal direction of the image plane, whereθ' is the angle between a chief ray and an optical axis in the image space. In this casez'_{m} (the tolerance limit of meridional field curvature) is expressed as follows:Let us consider an example of the design of a simple“
fθ ” lens according to the previous analysis of distortion and meridional field curvaturez'_{m} (in section 2 and 4):Table 2. shows the parameters of “
fθ ” lens:r is the radius of curvature,d is the corresponding axial distance,n is the name of the glass as given in the Shott catalog andS_{p} is the entrance pupil position, i.e. the distance between a center of the rotator facet (in the caseθ =0) and a first surface of the “fθ ” lens. We designed this system with two spherical lenses so that the nonuniformity of the spot movements(ZD ) and the meridional field curvature do not exceed1 %and0.3 mm respectively (see Fig. 5).In addition, let us investigate the influence of coma in the “
fθ ” lens on the spot position. In the case of the above preliminary “fθ ” lens the maximum value of coma is about0.0325 mm for the beam diameter of10 mm at the end of the field, and the shift of the effective spot center due to this effect is less than0.0325 mm . In fact, the shift of an effective spot center is0.0074 mm according to wave front criterion (minimum RMS wave front deformation) and0.0111 mm according to geometrical criterion (minimum RMS transverse aberration). When we compare these values withGD (0.459 mm ) at the end of a field, we can neglect the effect of coma. In general, the effect of coma can be disregarded in the system with the small aperture.V. MINIMIZATION OF THE RADIUS OF ROTATOR AND THE CHANGE OF ENTRANCE PUPIL POSITION
In order to minimize the change of an entrance pupil position and to make a compact system, it is very important to minimize the radius of the rotator. But, it is impossible to have an extremely small size of rotator because the radius depends on the input beam diameter
D .We want to find the minimum radius of rotator under the given conditions: a maximum scan angle
θ_{max} , a vertical angle of rotatorψ =π/n (wheren is the number of rotator facets), a shifting angleβ of the incident chief ray from the center of the rotator facet along the facet in the case ofθ =0 , the incident beam widthD , and a variation angleα of optical axis direction in the case ofθ =0 . From the geometry in Fig. 6 we can find the minimum radii of rotatorr _{+} andr _{} at which the incident beam reflects from the mirror of rotator without losing it, wherer _{+},r _{} are the minimum radii of rotator when the polygonal mirror rotates counterclockwise and clockwise respectively:From Eq.(22) we can see the smaller
β , the largerr _{+} and the smallerr _{}, but the size of rotator is determined by the greater value ofr _{+} andr _{} (see Fig. 7).In the case
β =0 ,r _{+} is always larger thanr _{} and the radius of rotator has to be chosen asr _{+}. From Fig. 7 we can see that the radius of rotator has a minimum value whenr _{+} is equal tor _{} and in this casetanβ can be expressed by the following form:From Eqs.(22) and (23) we can get a minimum radius
r :Of course, the smaller the
α , ther , i.e. we can reduce the radius of the rotator by the proper choice of parameters.In a design process, it is assumed that the direction of the incident beam changes centering around one fixed point on the rotator facet, therefore the entrance pupil position is usually located at this point. In general, the entrance pupil’s position is fixed at an intersection point between a chief ray and an optical axis in an optical system. However, in the “
fθ ” lens the entrance pupil position cannot be constant and identical to all field coordinates because of the geometry of beams deflected by the rotating polygonal mirror[12, 15]. So, we must consider the fact that the entrance pupil position of the “fθ ” lens is changed by the rotating polygonal mirror in the real laser scanning system (see Fig. 8).From Fig. 8 the change of the entrance pupil position Δ
S_{p} can be exactly expressed by the following form:where
r is the perpendicular distance from the center of the rotator to its facet, andθ is the scan angle. Here the scan angleθ is positive if the polygonal mirror rotates counterclockwise, and negative if it rotates clockwise.From the introduction of the angle
β we can very effectively reduce not only the radius of the rotator but also ΔS_{p} . For example, whenD =10 mm ,θ_{max} =45 °,ψ =60 ° andα =45 °,if we introduceβ of Eq.(23), thenβ is3.43 °. In this case, the minimum radiusr is35.91 mm and ΔS_{p} changes from1.18 mm to3.82 mm , but ifβ =0 , thenr =59.65 mm and ΔS_{p} changes from1.17 mm to5.46 mm .Note that Eqs.(22), (23), (24) and (25) are derived for the general case and without any approximation.
VI. CONCLUSIONS
We have shown that the distortion in the “
fθ ” lens should be evaluated with the peakvalley value ofZD instead of the maximum value ofGD . We derived how to correctZD of the “fθ ” lens with the relation between the 3rd and 5th order distortion and the ratio of the peakvalley value ofZD toGD . Also, we focused on the importance of the form of the global distortion curve and suggested several good forms for the correction ofZD in the “fθ ” lens.The spot size that is detected on the recording plane of a laser scanning system was analyzed not by aerial image but by binary image, and we have pointed out that the reference irradiance for binary image has to be selected as
I_{e} =I_{0}exp (1/2 ), i.e.k =0.5 when the change of the aerial spot sizeη is comparatively small.In order to decide a tolerance limit of meridional field curvature we utilized the analysis on the change of the binary spot size and the Gaussian beam properties. From these tools we obtained the general forms of Eqs.(20) and (21).
We know that the smaller the radius of the rotator in a laser printer system, the smaller the Δ
S_{p} and the more compact the system size. Therefore, it is natural that the radius of the rotator must be as small as possible. From Eqs.(23) and (24) we may always find the minimum radius of the rotator.

[FIG. 1.] Common fractional global distortion curve of thedistortionfree “fθ ” lens.

[TABLE 1.] TABLE 1. Peakvalley values of zonal and global distortion and coefficient C along Κ where um is the value at which the global distortion has an extreme value where we assume b=1

[FIG. 2.] Ratio of peakvalley value of zonal distortion to global versus Κ.

[FIG. 3.] FIG. 3. Global (dotted line) and zonal (solid line) distortioncurves: (a) Κ=5/3; (b) Κ=2; (c) Κ=10/3 (d) Κ=1; (e) Κ=0; (f) Κ=1.

[FIG. 4.] Beam spot size by binary image: (a) onaxis and (b)offaxis.

[TABLE 2.] Constructive data of “fθ” lens

[FIG. 5.] Global distortion and field curvature of the “fθ” lens.

[FIG. 6.] Scheme of rotator.

[FIG. 7.] Radii of rotator according to the shifting angle β where θmax=45° ψ=60°and α=90°.

[Fig. 8.] Change of the entrance pupil position ΔSp in rotatingthe polygonal mirror.