The interface between nanofillers and matrix plays a key role in determining the properties of nanocomposites, but the interfacial characteristics of nanocomposites such as molecular structure and interaction strength are not fully understood yet. In this work, the interfacial features of a typical nanocomposite, namely epoxy resin (EP) filled with boron nitride nanosheet (BNNS) are investigated by utilizing molecular dynamics simulation, and the effect of surface functionalization is analyzed. The radial distribution density (RDD) and interfacial binding energy (IBE) are used to explore the structure and bonding strength of nanocomposites interface. Besides, the interface compatibility and molecular chain mobility (MCM) of BNNS/EP nanocomposites are analyzed by cohesive energy density (CED), free volume fraction (FFV), and radial mean square displacement (RMSD). The results indicate that the interface region of BNNS/EP is composed of three regions including compact region, buffer region, and normal region. The structure at the interfacial region of nanocomposite is more compact, and the chain mobility is significantly lower than that of the EP away from the interface. Moreover, the interfacial interaction strength and compatibility increase with the functional density of BNNS functionalized by CH 3 -(CH 2 ) 4 -O-radicals. These results adequately illustrate interfacial characteristics of nanocomposites from atomic level. Some experiments have been performed to explore the interfacial interaction of nanocomposites [21][22][23][24]. The interfacial properties of nanocomposites such as the chain movement, density, and compatibility can be characterized by Raman spectroscopy, atomic force microscopy, and in situ transmission electron microscopy, which showed that the interface feature has a dominant effect on the properties of composites [25][26][27][28]. The effect of interface characteristics on mechanical and physical properties of CNTs nanocomposites are investigated by measuring its surface microstructures, energy-dispersive spectroscopy, and fracture toughness, which pointed out that improving interfacial properties was a green and promising approach toward preparing high-performance composites [29][30][31]. Q. Li and J. L. He employed a modified Kelvin probe force microscopy method to detect the local polarization property at the matrix/particle interface in ferroelectric nanocomposites, and the results illuminated that the abnormal performance of ferroelectric nanocomposites stems from the interfacial region [32]. To explain these experimental results, several interface theories and models have been proposed to describe the interfacial interaction in nanocomposites [33]. Lewis [34] and Tsagaropoulos [35] considered that the interfacial region could be split into two layers: a tightly bound layer (which does not contribute to the glass transition), and a loosely bound layer (which may exhibit its own glass transition unique from the rest of the polymer), which also indicated that the interfacial interactions...