The time evolution of major pyrolysis products including small-molecule species of a dipropargyl ether bisphenol A based novel boron-containing polymer were examined via ReaxFF-MD simulation to elucidate complicated mechanisms and the dominant reaction pathways on molecular level (Snapshot of a large graphitic precursor structure formed after a simulation run time of 2.Boron-containing polymers have recently attracted worldwide attention due to its ability to affect the thermal and oxidation resistance of the in situ formed material during pyrolytic carbonization. In this work, a novel dipropargyl ether bisphenol A based boron-containing polymer (PB), which can thermally cure via propargyl groups, was synthesized. PB was characterized via Fourier Transform infrared spectroscopy and nuclear magnetic resonance ( 1 H, 11 B and 13 C). Thermogravimetric analysis indicated the outstanding thermo-oxidative stability of PB thermosets with the temperature of 5% weight loss of 362 o C and char yield of 56.7% at 800 o C in air, while those are 416 o C and74.0%,respectively, in nitrogen. The time evolution of major pyrolysis products including pyroysis mechanism of PB thermosets were examined via thermogravimetry-Fourier transform infrared spectra as well as reactive molecular dynamics (ReaxFF-MD) simulations. In addition, the carbonized structure of PB thermosets was analyzed using X-ray photoelectron spectroscopy and X-ray powder diffraction. Formation of the graphitic fragment was simulated with ReaxFF-MD. Combined the ReaxFF-MD simulation with the experimental analysis, the way that boron atoms insert into the graphitic structure was illustrated.