  We simulated modifications to a model of a two-dimensional paramagnetic
semiconductor called the half-BHZ model, also known as the QWZ model, and
simulated a modified full BHZ model, where a time reversal pair is introduced.
Our modifications to the models include adding single and multiple impurities
connected to the lattices or as a connection between the time-reversal pairs.
We employed the Julia programming language to show how to speed up calculations
for time evolutions. By simulating the time evolutions, we could observe the
differences in the effects of these modifications. Our simulations showed the
presence of scattering behavior associated with the infinite QWZ model
topological states. Moreover, we observed scattering and absorption behavior
related to the parameters and placements of impurities and Hamiltonian
imaginary component's symmetry or anti-symmetry. These tools and early results
lay the foundations for developing electronic devices that use the models'
unique scattering and absorption behaviors and explore more complex and
physically accurate modifications to the models.
