Impact-resistance has been one of the target properties of polymer technology from a view point of their mechanical applications. A huge number of investigations have been done for this very important property. One of the key to understand the impact resistance is to investigate the difference in microscopic behavior between brittle and ductile fractures by comparing two extrema of industrial polymers such as PMMA and PC, respectively. However, in spite of an abundance of information about macroscopic behavior of the impact fracture, little has been known about its microscopic picture. This is from difficulty in experiment for the very rapid phenomena.Molecular dynamics (MD) calculation has been one of the powerful tools for investigating microscopic behavior of molecular assemblies at an atomistic resolution. However, so far, poor performance of the supercomputers has prevented us from investigating polymers based on all-atomistic MD calculations since the systems are huge. Then, we have been obliged to depend upon coarse-grained models such as Kremer-Grest model and DPD model. But, they are not good at describing chemical character of polymers which is essential for the investigation of fracture. However, recent very rapid progress of supercomputer technology is going to enable us to investigate polymers based on all-atomistic MD calculations free from phenomenological parameters. The all-atomistic model can describe chemical details of the polymers, which is essential for the investigation of the diversity of polymers. Fracture of the polymers depends much on their chemical details such as main-chain bond breaking and side-chain motions.In particular, it is very important to include bond fracture in the simulation. Since most of the traditional MD calculations have not included bond breaking explicitly, the calculations couldn’t describe brittle fracture of polymers but always resulted in artificial ductile fracture. The bond fracture can be described using potential function such as Morse-type one. The function can be obtained by high quality quantum chemical calculations prior to the MD calculations. In the present study, brittle fracture has been investigated for PMMA including bond fracture explicitly. Systems constructed for the simulation mimicked the real system. There, experimental molecular weight distribution, density, radius of gyration, and entanglement molecular weight were reproduced well in the initial configuration prepared for the impact fracture test calculations. The calculations clearly show characteristics of the brittle fracture at a molecular level such as void formation, chain reactions of main-chain bond breaking, and small Poisson ratio.