(color online) (a) Illustration of the NEMD process. Black atoms at the boundaries of the system are fixed in the position. Red and blue areas are denoted as heat bath and heat sink, respectively. (b) Temperature distribution along the heat flux direction in the armchair direction. The red solid line denotes the linear fitting results. Reproduced from Ref. [101] with permission from the Royal Society of Chemistry.