Optimization of pyrolysis tube cracking efficiency

Abstract

A method for optimizing a cracking efficiency with which a pyrolysis tube of a deposition apparatus cracks a precursor material into reactive species is disclosed, including measuring an input pressure at an entrance to the pyrolysis tube, outside of the pyrolysis tube; measuring an output pressure at an exit from the pyrolysis tube, outside of the pyrolysis tube; measuring a pyrolysis temperature within the pyrolysis tube; calculating a cracking efficiency based on the input pressure, the output pressure and the pyrolysis temperature; and determining an adjustment to be made to at least one of the input pressure, the output pressure and the pyrolysis temperature to increase the cracking efficiency.

Claims

1. A method for optimizing a cracking efficiency with which a pyrolysis tube of a deposition apparatus cracks a precursor material into reactive species, comprising: measuring an input pressure of the pyrolysis tube; measuring an output pressure of the pyrolysis tube; measuring a pyrolysis temperature within the pyrolysis tube; calculating a cracking efficiency based on the input pressure, the output pressure, and the pyrolysis temperature; and adjusting at least one of the input pressure, the output pressure, or the pyrolysis temperature to maintain a particular cracking efficiency, wherein adjusting comprises increasing or decreasing one of the input pressure of the pyrolysis tube, the output pressure of the pyrolysis tube, or the pyrolysis temperature within the pyrolysis tube to maintain the particular cracking efficiency.

2. The method of claim 1, further comprising: measuring a vaporization temperature at the entrance to the pyrolysis tube, outside of the pyrolysis tube; calculating the cracking efficiency based on the input pressure, the output pressure, the vaporization temperature, or the pyrolysis temperature.

3. The method of claim 1, further comprising: determining a pressure differential between the input pressure and the output pressure.

4. The method of claim 3, further comprising: calculating a residence time that the precursor material remains within the pyrolysis tube with the pressure differential.

5. The method of claim 4, wherein calculating the residence time comprises employing the following formula:
Δt=t.sub.0(P.sub.1−P.sub.2).sup.−β, where Δt is the residence time, (P.sub.1−P.sub.2) is pressure differential between the input pressure and the output pre, to is an experimentally determined pre-factor for the residence time and corresponds specifically to the pyrolysis tube and β is an experimentally determined exponent that corresponds specifically to the pyrolysis tube.

6. The method of claim 5, wherein calculating the cracking efficiency comprises employing the following formula:
η=1−e.sup.−λΔt, where η is cracking efficiency, e is Euler's number and λ is a reaction rate at which molecules of the precursor material are cracked into reactive species.

7. The method of claim 1, further comprising identifying an adjustment to be made that comprises identifying a residence time that is at least as great as a predetermined minimum threshold residence time.

8. The method of claim 7, further comprising identifying an adjustment to be made that comprises identifying a residence time that is at least 0.1 second.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIGS. 1A and 1B are representations of the processes used to produce a polymer in accordance with the prior art.

(2) FIG. 2 illustrates material deposition equipment in accordance with an exemplary embodiment of the disclosed subject matter.

(3) FIG. 3 is graph illustrating cracking efficiency as a function of pyrolyzer temperature in accordance with an exemplary embodiment of the disclosed subject matter.

(4) FIG. 4 is a graph illustrating cracking efficiency as a function of pyrolyzer temperature in accordance with an exemplary embodiment of the disclosed subject matter.

(5) FIG. 5 is a graph illustrating cracking efficiency and conversion efficiency as a function of pyrolyzer temperature in accordance with an exemplary embodiment of the disclosed subject matter.

(6) FIG. 6 is a graph illustrating conversion efficiency as a function of pyrolyzer temperature in accordance with an exemplary embodiment of the disclosed subject matter.

(7) FIG. 7 is a graph illustrating the dependence of cracking efficiency on temperature at different residence times in accordance with an exemplary embodiment of the disclosed subject matter.

(8) FIG. 8 illustrates examples of variation of cracking efficiency during actual depositions of polyp-xylylene) in accordance with an exemplary embodiment of the disclosed subject matter.

(9) FIG. 9 illustrates a representation of pressures and residence times from experimental data in accordance with an exemplary embodiment of the disclosed subject matter.

DETAILED DESCRIPTION

(10) Reference will now be made in detail to exemplary embodiments of the disclosed subject matter, examples of which are illustrated in the accompanying drawings. The method and corresponding steps of the disclosed subject matter will be described in conjunction with the detailed description of the system.

(11) During deposition of the polymer, the preference is for the deposition of a strong and clear coating. However, dynamics within the pyrolyzer may affect the characteristics of the polymer.

(12) The performance of the pyrolyzer may be characterized in terms of its cracking efficiency, η, defined as the ratio of the amount of dimer that is pyrolyzed to the total dimer load for a given deposition process,

(13) $\begin{array}{cc}\eta =\frac{M_i-M_f}{M_i}& \left(\mathrm{Eq}.1\right)\end{array}$
where M.sub.i and M.sub.f are the amounts of dimer introduced into the pyrolyzer and exiting the pyrolyzer, respectively. In an ideal case, cracking efficiency η should be 100% while a polymer is being deposited onto a surface of a substrate.

(14) In practical operation, there are difficulties associated with operation of the pyrolyzer to achieve and maintain 100% cracking efficiency. Two undesired scenarios may arise, namely “under-cracking” and “over-cracking.”

(15) When under-cracking occurs, a mixture of un-cracked dimer and monomer exits the pyrolyzer, with the un-cracked dimer quickly condensing onto surfaces between the pyrolyzer and the deposition chamber and onto cooler surfaces in the deposition chamber. The heavy dimer deposits that are formed on these surfaces can adversely affect performance of the material deposition equipment, and undesirably increase the frequency with which the material deposition equipment must be cleaned and maintained.

(16) Under-cracking also can result in low dimer utilization, meaning that the polymer coating may be undesirably thin or that the formation of a polymer coating of desired thickness may be undesirably slow. Additionally, when the dimer is not sufficiently cracked by the pyrolyzer, dimer may be present within, and contaminate, the polymer coating, which may undesirably affect the quality of the polymer coating.

(17) When over-cracking occurs, the monomer molecules are further dissociated into degraded monomers with lower molecular weights, as well unwanted products such as chlorine gas (Cl.sub.2). Chlorine gas and moisture can readily form hydrochloric acid (HCl), which can corrode the surfaces of the material deposition equipment that the hydrochloric acid contacts. In addition, degraded monomers may deposit onto surfaces of the material deposition equipment, which can adversely affect its performance and undesirably increase the frequency with which the equipment must be cleaned and maintained.

(18) An embodiment of material deposition equipment 10 in accordance with an embodiment of this disclosure is depicted by FIG. 2. The material deposition equipment 10 is useful in connection with embodiments of the methods described herein. The material deposition equipment includes a vaporizer 12, a pyrolyzer 14 and a deposition chamber 16. The vaporizer 12 receives a precursor material, such as a parylene dimer, and vaporizes or sublimates the precursor material. The pyrolyzer 14, which includes a heater 18, cracks the precursor material into active species, such as monomers. The deposition chamber 16 holds one or more substrates, and receives active species from the pyrolyzer 14 to enable the active species to deposit onto surfaces of the one or more substrates and to form a polymer on surfaces of the one or more substrates.

(19) In FIG. 2, P.sub.1 represents the input pressure at an entrance to the pyrolysis tube 14, outside of the pyrolysis tube, or the “inlet pressure,” and P.sub.2 represents the output pressure at an exit from the pyrolysis tube 14, outside of the pyrolysis tube, or “exit pressure.” T.sub.v represents the vaporization temperature at the entrance to the pyrolysis tube, outside of the pyrolysis tube, or the “inlet temperature.” T.sub.p represents the temperature within the pyrolyzer tube 14, or the “pyrolyzer temperature.” Temperature may be measured in any suitable manner. Without limitation, the pyrolyzer temperature, T.sub.p, may be measured by thermocouples that are in close proximity to the wall of the pyrolyzer.

(20) The cracking efficiency η of the pyrolyzer is represented by the following equation:
η=1−e.sup.−λΔt   (Eq. 2)
where Δt is the residence time of dimer molecules in the pyrolyzer 14 and λ is the reaction rate. In terms of the process parameters identified in FIG. 2 (i.e., P.sub.1, P.sub.2, T.sub.v, and T.sub.p), Δt and λ are represented by the following equations:

(21) $\begin{array}{cc}\Delta t={t_0\left(P_1-P_2\right)}^{-\beta }& \left(\mathrm{Eq}.3a\right)\\ \lambda =4\mathrm{\rho \pi }{d}^2\frac{P_1}{{\mathrm{kT}}_v}\sqrt{\frac{8k{T}_p}{\pi \mu }}{e}^{-\frac{E_a}{{\mathrm{kT}}_p}}& \left(\mathrm{Eq}.3b\right)\end{array}$
where: t.sub.0 a pre-factor for residence time, is dependent on specific material deposition equipment (see exemplary embodiment at APPENDIX B); β is an experimentally determined exponent, which is also tool dependent (e.g., β≈0.9 on the material deposition equipment used to provide the data disclosed herein) (see exemplary embodiment at APPENDIX B); ρ a steric factor (e.g., ρ˜1); d an equivalent diameter of the dimer molecule (e.g., d≈1 nm); k Boltzmann constant; μ reduced mass of dimer and monomer (e.g., μ=92.3 a.u. for parylene C dimer); E.sub.a activation energy (e.g., E.sub.a≈1.25 eV for parylene C dimer; E.sub.a may be about 0.98 eV for other precursor materials; E.sub.a may be 0.1 to 2.0 for various other precursor materials.)

(22) As indicated by Eq. 3b and Eq. 2, above, the cracking efficiency, η, is dependent upon the pyrolyzer temperature, T.sub.p. This relationship is further illustrated by the graph of FIG. 3, which shows cracking efficiency data obtained using WB5000 material deposition equipment available from HZO, Inc., of Draper, Utah, with parylene C at a residence time, Δt, of 0.12 s. As illustrated by the dotted lines in FIG. 3, when the pyrolyzer of that material deposition equipment is heated to a pyrolyzer temperature, T.sub.p, of 590° C., the pyrolyzer and, thus, of the material deposition equipment cracks parylene C with an efficiency of about 90%.

(23) The cracking efficiency may change with a change in residence time, Δt. As an example, when residence time, Δt, is increased to 0.4 s., the cracking efficiency η may increase at a variety of different pyrolyzer temperatures, T.sub.p. This result is illustrated by the graph of FIG. 4, in which the curve represents a calculated cracking efficiency η (i.e., “percent conversion”) and the symbols are experimental data points from Fortin, J. B., et al., J. Vac. Sci. Technol., A 18, 2459 (2000), wherein the cracking efficiency was determined with different materials as the liner in their pyro tube. The cracking efficiency η may be almost 100% at a pyrolyzer temperature, T.sub.p, of 550° C. In some embodiments, the cracking efficiency may be set at a value of less than 100% (e.g., 98%, 95%, 92%, 90%, etc.). Even at a pyrolyzer temperature, T.sub.p, of 450° C., nearly 20% of dimer is converted into monomers.

(24) Notably, Eq. 2 does not provide any indication regarding the extent over-cracking. In fact, when cracking efficiency η approaches 100%, significant over-cracking is expected to occur. The effective conversion efficiency χ can be defined as

(25) $\begin{array}{cc}\chi =\frac{\mathrm{total}\mathrm{amount}\mathrm{of}\mathrm{non}\text{-}\mathrm{degraded}\mathrm{monomer}}{2\times \mathrm{total}\mathrm{amount}\mathrm{of}\mathrm{dimer}}& \left(\mathrm{Eq}.4\right)\end{array}$
APPENDIX A (eq. A18) provides an approximation of χ as
χ=(1−e.sup.−λΔt)e.sup.−λ′Δt   (Eq. 5)
where λ is defined in Eq. 3b, and

(26) $\begin{array}{cc}{\lambda }^{\prime }=4\mathrm{\rho \pi }{d}^2\frac{P_1}{{\mathrm{kT}}_v}\sqrt{\frac{8{\mathrm{kT}}_p}{\mathrm{\pi \mu }}}{e}^{-\frac{E_b}{{\mathrm{kT}}_p}}& \left(\mathrm{Eq}.6\right)\end{array}$

(27) FIGS. 5 and 6 show the effect of over-cracking on the effective conversion efficiency:

(28) FIG. 5 shows the difference between cracking efficiency η and effective conversion efficiency χ. Line 50 represents cracking efficiency η with no over-cracking and line 52 represents cracking efficiency η with over-cracking. (Cracking efficiency η denoted on Y-axis on left side of FIG. 5) Line 54 represents the difference between no over-cracking and with over-cracking, which is the conversion efficiency χ. (Conversion efficiency χ denoted on Y-axis on right side of FIG. 5.)

(29) FIG. 6 shows the effect of residence time when over-cracking is considered. It can be seen that: (a) over-cracking occurs at a wide range of pyrolyzer temperatures; (b) the longer the residence time Dt, the lower the onset temperature of over-cracking; and (c) the maximum effective conversion efficiency is usually lower than 100%.

(30) FIG. 7 shows the dependence of cracking efficiency on temperature at different residence times Δt due to different pressure drops P.sub.1−P.sub.2 across the pyrolyzer. It can be seen that, at a given pyrolyzer temperature T.sub.p when the residence time Δt decreases, the cracking efficiency also decreases. For example, when the pyrolyzer temperature T.sub.p is 600° C., at a residence time Δt of 1 s, the cracking efficiency is 100%; at a residence time Δt of 0.1 s, the cracking efficiency is about 90%; and at a residence time, Δt, of 0.01 s, =0.1 s, the cracking efficiency is only about 25%.

(31) FIG. 8 illustrates examples of variation of cracking efficiency during actual depositions of polyp-xylylene) using WB5000 material deposition equipment available from HZO, Inc., of Draper, Utah. The cracking efficiency were determined over the course of three different deposition courses, a first at a pyrolyzer temperature, T.sub.p, of 550° C., a second at a pyrolyzer temperature, T.sub.p, of 585° C. and a third at a pyrolyzer temperature, T.sub.p, of 625° C. Due to uncertainty over the rate constant for over-cracking, the calculations represented by the curves of FIG. 8 are calculated without correction for over-cracking. If over-cracking were taken into account, the effective conversion efficiency would be lower than those illustrated by FIG. 8.

(32) When polyp-xylylene) is deposited, even at a fixed pyrolysis temperature, T.sub.p, as depicted by FIG. 8, under-cracking and over-cracking will occur at different points or stages during the deposition process. Accordingly, in some embodiments, a model according to this disclosure may be used to repeatedly or even continuously optimize cracking efficiency throughout the course of a deposition process, and one or both of the pyrolyzer temperature and the residence time, Δt, (i.e., the inlet pressure, P.sub.1, and/or the exit pressure, P.sub.2) may be adjusted to maintain or substantially maintain a particular (e.g., a maximum, etc.) cracking efficiency.

(33) A semi-empirical model for optimizing cracking efficiency may be based on collision theory. This model can be applied to process optimization by real time monitoring and adjustment of process parameters.

APPENDIX A

Derivation of Cracking Efficiency

(34) The process of pyrolysis can be understood within the general framework of the collision theory of gas phase chemical reactions. The vapor in a pyrolyzer is usually a mixture of dimer and monomer molecules. For C-dimer, denoting the dimer and monomer as D and M respectively, the possible reactions are
D+D.fwdarw.4 M   (A1)
M+D.fwdarw.M+2M   (A2)
M+M.fwdarw.D   (A3)
M+M.fwdarw.M+M*+Cl   (A4)

(35) Reaction A4 refers to over-cracking, in which the C—Cl bond in a monomer is broken. M* denotes degraded monomer due to over-cracking.

(36) For parylene C dimer, at a pyrolysis temperature, T.sub.p, of more than 500° C., the re-dimerization (reaction A3) can be neglected. Also, since the bond strength of C—Cl bond is much stronger than the strained C—C bond, it can be assumed that no over-cracking occurs within this framework.

(37) Under these assumptions, the rate equation becomes,

(38) $\begin{array}{cc}\frac{d\left[M\right]}{\mathrm{dt}}=4{k_{\mathrm{DD}}\left[D\right]}^2+2k_{\mathrm{MD}}\left[D\right]\left[M\right]& \left(\mathrm{eq}.\mathrm{A1}\right)\end{array}$
where [D] and [M] are the molar concentration of the dimer and monomer molecules. They are related as

(39) $\begin{array}{cc}\left[D\right]+\frac{\left[M\right]}{2}={\left[D\right]}_0& \left(\mathrm{eq}.\mathrm{A2}\right)\end{array}$
where [D].sub.0 is the initial dimer concentration at pyrolyzer entry. Combining Eq. A1 and Eq. A2, the rate equation becomes

(40) $\begin{array}{cc}\frac{d\left[M\right]}{\mathrm{dt}}=4\left\{\left(k_{\mathrm{DD}}-k_{\mathrm{MD}}\right)\left[D\right]+{k_{\mathrm{MD}}\left[D\right]}_0\right\}.Math.\left[D\right]& \left(\mathrm{eq}.\mathrm{A3}\right)\end{array}$

(41) The two rate constants, k.sub.MD and k.sub.DD, are

(42) 0$\begin{array}{cc}{k}_{\mathrm{MD}}=\mathrm{\rho \pi }{d}_{\mathrm{MD}}^2\sqrt{\frac{8{\mathrm{kT}}_p}{{\mathrm{\pi \mu }}_{\mathrm{MD}}}}{e}^{-\frac{E_a}{{\mathrm{kT}}_p}}.Math.N_A& \left(\mathrm{eq}.4A\right)\end{array}$

(43) $\begin{array}{cc}{k}_{\mathrm{DD}}=\mathrm{\rho \pi }{d}_{\mathrm{DD}}^2\sqrt{\frac{8{\mathrm{kT}}_p}{{\mathrm{\pi \mu }}_{\mathrm{DD}}}}{e}^{-\frac{E_a^{\prime }}{{\mathrm{kT}}_p}}.Math.N_A& \left(\mathrm{eq}.\mathrm{A5}\right)\end{array}$
where: ρ is a steric factor (e.g., ρ˜1); d.sub.MD is a radius of the cross section for a dimer-monomer collision; d.sub.DD is a radius of the cross-section for a dimer-dimer collision; T.sub.p is the pyrolyzer temperature; μ.sub.MD is a reduced mass of dimer-monomer pair; μ.sub.DD is a reduced mass of dimer-dimer pair; E.sub.a is an activation energy for a dimer-monomer collision; E′.sub.a is an activation energy for a dimer-dimer collision; k is the Boltzmann constant; and N.sub.A is the Avagadro constant.

(44) For parylene C dimer, k.sub.MD and k.sub.DD are rather similar, therefore the first term in Eq. A3 is much smaller than the second term, i.e.,
(k.sub.DD−k.sub.MD)[D]«k.sub.MD[D].sub.0   (eq. A6)
Therefore, the rate equation reduces to

(45) $\begin{array}{cc}\frac{d\left[M\right]}{\mathrm{dt}}=4{k_{\mathrm{MD}}\left[D\right]}_0.Math.\left[D\right]& \left(\mathrm{eq}.\mathrm{A7}\right)\end{array}$

(46) Note that

(47) $\begin{array}{cc}\frac{d\left[M\right]}{\mathrm{dt}}=-2\frac{d\left[D\right]}{\mathrm{dt}}& \left(\mathrm{eq}.\mathrm{A8}\right)\end{array}$
We have,

(48) $\begin{array}{cc}\frac{d\left[D\right]}{\mathrm{dt}}=-2{k_{\mathrm{MD}}\left[D\right]}_0.Math.\left[D\right]& \left(\mathrm{eq}.\mathrm{A9}\right)\end{array}$
[D].sub.0 can be calculated from process parameters as,

(49) $\begin{array}{cc}{\left[D\right]}_0=\frac{P_1}{{\mathrm{RT}}_v},& \left(\mathrm{eq}.\mathrm{A10}\right)\end{array}$
where R is the ideal gas constant, R=8.31 J.Math.K.sup.−1.Math.mol.sup.−1. T.sub.v is the vapor temperature at the pyrolyzer entry. Let λ=2 k.sub.MD[D].sub.0, the rate equation becomes

(50) $\begin{array}{cc}\frac{d\left[D\right]}{\mathrm{dt}}=-\lambda \left[D\right]& \left(\mathrm{eq}.\mathrm{A11}\right)\end{array}$

(51) By solving this equation and integrating from t.sub.0 to (t.sub.0+Δt), where Δt is the residence time of the vapor molecules, we have
[D]=[D].sub.0e.sup.−λΔt   (eq. A12)

(52) The cracking efficiency, η, can be calculated as

(53) $\begin{array}{cc}\eta =\frac{{\left[D\right]}_0-\left[D\right]}{{\left[D\right]}_0}=1-e^{-\mathrm{\lambda \Delta }t}& \left(\mathrm{eq}.\mathrm{A13}\right)\end{array}$

(54) When deriving the cracking efficiency (Eq. A13), the loss of monomer to over-cracking is not considered. This is justified because the C—Cl bond is much stronger than the strained C—C bond, and requires much higher activation energy in the rate constant. However, when cracking efficiency is high, the vapor near the pyrolyzer exit consists of dominantly monomer molecules, and over-cracking can be significant. Here, an estimate of over-cracking efficiency, χ, is given.

(55) Following similar procedures, the rate equation for over-cracking is

(56) $\begin{array}{cc}\frac{d\left[M\right]}{\mathrm{dt}}=4{k_{\mathrm{DD}}\left[D\right]}^2+2k_{\mathrm{MD}}\left[D\right]\left[M\right]-2{k_{\mathrm{MM}}\left[D\right]}_0\left[M\right]& \left(\mathrm{eq}.\mathrm{A14}\right)\\ k_{\mathrm{MM}}=\mathrm{\rho \pi }{d}_{\mathrm{MM}}^2\sqrt{\frac{8{\mathrm{kT}}_p}{{\mathrm{\pi \mu }}_{\mathrm{MM}}}}{e}^{-\frac{E_b}{{\mathrm{kT}}_p}}.Math.N_A& \left(\mathrm{eq}.A15\right)\end{array}$

(57) Here, only the collisions between monomers are considered to contribute to over-cracking, since collisions between dimer and monomer will most likely produce monomers. Also, it is assumed that both “good” monomer and degraded monomer contribute to over-cracking. E.sub.b in Eq. A15 is the activation energy of over-cracking, and E.sub.b>E.sub.a.

(58) Specifically, for parlyene C dimer:
k.sub.MM«k.sub.MD   (eq. A16)
Under this condition, an approximate form of [M] can be written as
[M]=2[D].sub.0(1−e.sup.−λΔt)e.sup.−λ′t   (eq. A17)
and the effective conversion efficiency is

(59) $\begin{array}{cc}\chi =\frac{\left[M\right]}{{2\left[D\right]}_0}=\left(1-e^{-\lambda t}\right)e^{-{\lambda }^{\prime }t}& \left(\mathrm{eq}.\mathrm{A18}\right)\end{array}$
Where
λ′=2k.sub.MM[D].sub.0   (eq. A19)

APPENDIX B

Determination of the Values of t.SUB.0 .and β in Eq. (3a) on R&D WB5000

(60) The residence time at different dimer loads was determined, and the data is fitted according to Eq. 3a, above. The data analyzed included the residence time Δt of a plurality of different precursor material load amounts in the pyrolysis tube. A logarithm of the residence time corresponding to each precursor material load amount of the plurality of precursor material load amounts was then calculated. A pressure drop (P.sub.1−P.sub.2) within the pyrolysis tube corresponding to each precursor material load amount of the plurality of precursor material load amounts was determined, and a logarithm of the pressure drop within the pyrolysis tube corresponding to each precursor material load amount was calculated.

(61) Subsequently, data points for each precursor material load amount were plotted Each data point was a plot of the logarithm of the residence time Δt corresponding to the precursor material load amount against the logarithm of the pressure drop (P1−P2) within the pyrolysis tube corresponding to the precursor material load amount. A best fit line was determined through a plurality of data points corresponding to the plurality of precursor material load amounts.

(62) An equation for the best fit line was set equal to
Δt=t.sub.0(P.sub.1−P.sub.2).sup.−β  (Eq. 3a)
where t0 is the pyrolysis tube dependent pre factor for residence time and β is the pyrolysis tube dependent exponential constant.

(63) Table 1 shows the experimental data. FIG. 9 shows the fit to the data according to Eq. 3a. The values of t.sub.0 and β are found to be:
t.sub.0=41.7 s
β=0.9

(64) TABLE-US-00001 TABLE 1 Residence time as a function of the pressure drop, P.sub.1-P.sub.2. Res. log Time log P1-P2 1/(P1-P2) (P1-P2) (s) (Res time) Fit 751 0.001332 2.875639937 0.104533 −0.980745302 −0.96808 759 0.001318 2.880241776 0.108339 −0.965213308 −0.97222 799 0.001252 2.902546779 0.099307 −1.003018425 −0.99229 833.5 0.0012 2.920905604 0.083313 −1.079288996 −1.00882 617 0.001621 2.790285164 0.118616 −0.925856391 −0.89126 871 0.001148 2.940018155 0.089701 −1.04720038 −1.02602 874.5 0.001144 2.941759814 0.080277 −1.095408567 −1.02758 973.2 0.001028 2.9882021 0.081653 −1.088028426 −1.06938 981.7 0.001019 2.991978791 0.081087 −1.091047949 −1.07278 876.5 0.001141 2.94275192 0.089249 −1.049398169 −1.02848 472.78 0.002115 2.674659096 0.150622 −0.822110256 −0.78719 679.5 0.001472 2.832189461 0.132871 −0.876569915 −0.92897 642.3 0.001557 2.807737922 0.13505 −0.869505354 −0.90696 620.64 0.001611 2.792839762 0.131052 −0.88255775 −0.89356 652.4 0.001533 2.814513952 0.119542 −0.922477871 −0.91306 784.8 0.001274 2.894758994 0.122159 −0.91307577 −0.98528 703.5 0.001421 2.847264102 0.121932 −0.913882493 −0.94254 529 0.00189 2.723455672 0.147983 −0.829787 −0.83111 522 0.001916 2.717670503 0.152662 −0.816267722 −0.8259 547 0.001828 2.737987326 0.135026 −0.869583368 −0.84419 499 0.002004 2.698100546 0.146723 −0.833500528 −0.80829 586 0.001706 2.767897616 0.111628 −0.95222858 −0.87111 636 0.001572 2.803457116 0.111754 −0.951738657 −0.90311 655 0.001527 2.8162413 0.123517 −0.908272963 −0.91462

(65) While the disclosed subject matter is described herein in terms of certain preferred embodiments, those skilled in the art will recognize that various modifications and improvements may be made to the disclosed subject matter without departing from the scope thereof.

(66) In addition to the specific embodiments claimed below, the disclosed subject matter is also directed to other embodiments having any other possible combination of the dependent features claimed below and those disclosed above. As such, the particular features presented in the dependent claims and disclosed above can be combined with each other in other manners within the scope of the disclosed subject matter such that the disclosed subject matter should be recognized as also specifically directed to other embodiments having any other possible combinations. Thus, the foregoing description of specific embodiments of the disclosed subject matter has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosed subject matter to those embodiments disclosed.