Electron beam lithography is the most promising technology to realize the next generation LSI, but its low throughput is the most crucial issue. The simplest solution of the issue is to increase the beam current. However, as the current increases, the Coulomb interaction among electrons becomes severe, and a blur at the exposure pattem increases. In order to find the most preferable condition, it is necessary to quantify the Coulomb effect. In order to formulate the Coulomb effects, several approaches have been proposed. (l):By slicing the beam optics along the axis, the amount of the electron deviation in a segment is given by an equation. [l] (2):Considering only nearest-neighbor electron and a statistical equation is derived to express electron deviation after passing through the optics.[2] (3):Instead of deriving analytical equations, numerical treatment using Monte Carlo simulation can be flexibly applied.[3] However, since it calculates the Coulomb force between all combinations of existent electrons in the optics, as the beam current increases, the number of combination becomes large and the computation time is quite large. It is known that it is effective to reduce the calculation time by adopting a tree-code programming in high current beam, such as more than 20pA. [4] It is understood that the Coulomb interaction between electrons in the optical system can be categorized into two effects: One is the global Coulomb effect, and the other is the stochastic Coulomb effect. Electron beam in the optics produces a potential distribution as a whole, and every individual electron is influenced by the distribution as the global Coulomb effect. On the other hand, electrons in the beam may occasionally meet the other electron in a very short distance, and they will be recoiled each other, and large energy and angular deviations are found, as the stochastic Coulomb effect. In the present study we propose a new approach to consider the Coulomb effect by separating the global and the stochastic mechanisms, and then, they are combined in a Monte Carlo electron trajectory simulation. As the global Coulomb effect, the potential distribution is obtained in the beam, and the electron deflection is calculated due to the distribution. As the stochastic Coulomb effect, it is assumed that electrons interact with only their nearest-neighbor electron in the beam. All electron trajectories in the whole optical system are traced three-dimensionally considering both the potential distribution in the beam and the recoiled electron caused by the nearest neighbor electron. Since the Coulomb interaction is calculated for only one electron instead of all other electrons around the electron, the simulation time is remarkably reduced.In this model calculation, a round pattem beam with the radius of 200pm is emitted at 200" above the wafer with the energy of 2keV at 1 pA, and it is simply projected as a parallel beamThe validity of the present simulation is demonstrated in Fig.1. 112