and Shobatake K., “ Multistage Gas–Surface Interaction Model for the Direct Simulation Monte Carlo Method,” Physics of Fluids, Vol. 11, No. 11, Nov. 1999, pp. 3540–3552. doi: JPCAFH 1089-5639 Crossref Google Scholar and Yakovlev A., “ Molecular Dynamics Simulations of Laser-Induced Incandescence of Soot Using an Extended ReaxFF Reactive Force Field,” Journal of Physical Chemistry A, Vol. 114, No. 48, Dec. 2010, pp. 12561–12572. and Liu F., “ Molecular Dynamics Simulations of Translational Thermal Accommodation Coefficients for Time-Resolved LII,” Applied Physics B, Vol. 94, No. 1, Jan. 2009, pp. 39–49. doi: IJHMAK 0017-9310 Crossref Google Scholar Duan K., “ Thermal Accommodation Coefficients Between Polyatomic Gas Molecules and Soot in Laser-Induced Incandescence Experiments,” International Journal of Heat and Mass Transfer, Vol. 52, Nos. 21–22, Oct. 2009, pp. 5081–5089. J., Computer Simulation of Liquids, Oxford Science Publications, Great Clarendon St., Oxford, 1989, pp. 1–385. doi: JPCHAX 0022-3654 Crossref Google Scholar O., “ Thermal Accommodation Coefficients,” Journal of Physical Chemistry, Vol. 84, No. 12, June 1980, pp. 1431–1445. Baule B., “ Theoretische Behandlung der Erscheinungen in Verdunnten Gasen,” Annalen der Physik, Vol. 349, No. 9, 1914, pp. 145–176. Y., “ Formula for Thermal Accommodation Coefficients,” Journal of Chemical Physics, Vol. 46, No. 6, 1967, pp. 2376–2386. Y., Dynamics of Gas–Surface Scattering, Academic Press, New York, 1976, pp. 1–327. doi: JFLSA7 0022-1120 Crossref Google Scholar G., “ Direct Simulation Monte Carlo Calculations of Rarefied Flows with Incomplete Surface Accommodation,” Journal of Fluid Mechanics, Vol. 239, Jan. 1992, pp. 449–459. doi: AIAJAH 0001-1452 Link Google Scholar and Lampis M., “ New Scattering Kernel for Gas–Surface Interaction,” AIAA Journal, Vol. 35, No. 6, June 1997, pp. 1000–1011. doi: CESCAC 0009-2509 Crossref Google Scholar R., “ Knudsen’s Cosine Law and Random Billiards,” Chemical Engineering Science, Vol. 59, No. 7, March 2004, pp. 1541–1556. R., “ The Kinetic Theory of Gases: Some Modern Aspects,” Journal of Physical Chemistry, Vol. 39, No. 2, 1935, p. 307. doi: JPPOEL 0748-4658 Link Google Scholar and Onofri M., “ Chemical Erosion of Carbon-Phenolic Rocket Nozzles with Finite-Rate Surface Chemistry,” Journal of Propulsion and Power, Vol. 29, No. 5, Sept. 2013, pp. 1220–1230. and Rakich J., “ Catalytic Surface Effects Experiment on Space Shuttle,” AIAA 16th Thermophysics Conference, AIAA, New York, June 1981, pp. 1–16. and Gordeev A., “ Catalysis Effects on Quartz Surface in High Enthalpy Sub-Sonic Oxygen and Carbon-Dioxide Flows,” Proceedings of the 3rd European Symposium of Aerodynamics for Space Applications, Vol. 3, ESA Special Publ., Noordwijk, The Netherlands, Dec. 1998, pp. 537–544. and Muylaert J., “ Comparative Study of Surface Catalycity Under Subsonic Air Test Conditions,” Proceedings of the 4th European Symposium of Aerodynamics for Space Applications, Vol. 4, ESA Special Publ., Capua, Italy, Oct. 2001, pp. 481–488. A., Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford Science Publications, Great Clarendon St., Oxford, 1994, pp. 1–458. A direct-velocity-sampling gas–surface interaction model was implemented in direct-simulation Monte Carlo, which was found to predict flow and surface properties comparable to the Maxwell gas–surface interaction model when the molecular dynamics data used for velocity sampling were obtained from a rough (diffuse) surface. The collision statistics for the three categories were found to depend on the incidence angle and the speed, and a strong correlation among incidence speed, angle, and surface topology in the probability distribution of the energy accommodation coefficients was observed. Based on the postcollision behavior of the nitrogen molecule, the molecular dynamics trajectories were classified into three categories, namely, single, and multiple with and without escape. The trajectory molecular dynamics simulations were performed on these two surfaces to study the difference between atomistically smooth and rough cases. In the present work, molecular dynamics is used to study molecular nitrogen impinging on multilayered graphene and fused-quartz surfaces at different incidence speeds and angles to obtain energy accommodation coefficients for use in direct-simulation Monte Carlo.