Simulations with a quantum molecular dynamics (QMD) method ( MD with MO) were demonstrated on the thermal decomposition of lignin monomer at the ground state including excited and positive charged states. Geometry and energy optimized results of the lignin monomer at the singlet and triplet states in single excitation, and at (+2) positive charged state by semi-empirical AM1 MO calculations were used as the initial MD step of QMD calculations. In the QMD calculations, we controlled the total energy of the system using Nose-Hoover thermostats in the total energy range of 0.69-0.95 eV, and the sampling position data with a time step of 0.5 fs were carried out up to 5000 steps at 50 different initial conditions. The calculated neutral, positive and negative charged fragment distributions of the monomer model with 0.82 eV energy control were obtained as 90.6, 3.5, and 5.9% to the total fragments, respectively. The ratios seem to correspond well with to the values observed experimentally in SIMS. Crown Copyright (C) 2008 Published by Elsevier B. V. All rights reserved.