Modeling Bidirectional Thalamocortical Interactions using NetPyNE