The properties of the glass fiber reinforced epoxy composites are greatly influenced by the interaction of the fiber with the epoxy matrix. This interaction is quantified using the interphase traction-separation response, called traction law. Recently, we have developed a molecular dynamics (MD) simulation protocol to predict the interphase thickness and the traction law of the glass fiber-epoxy interphase [Ref. 1]. Interphase properties depend on several factors like silane number density, distribution pattern of the silane molecules on the fiber surface and strain rate, etc. For instances, it can be expected that the interphase with no silane and fully covered silane conditions will give the lower and upper bounds of properties, respectively. In this study, we would like to develop strain rate dependent mixed-mode traction law as a function of silane number density using the developed MD protocol in Ref. 1. Detailed structure and damage modes of the interphase at the atomistic scale will be reported. MD predicted traction laws will be used to bridge length scale in the continuum level micro-mechanics modeling.