The response to small signals with a general time dependence can be obtained from the knowledge of the impulse response. The static zero field mobility is
given by the long time limit of the differential velocity step response. This is exploited to derive a zero field Monte Carlo
algorithm for the mobility tensor from the algorithm presented in the previous section. For a vector-valued physical quantity elements of the
differential step response tensor are related to the differential impulse response tensor in the following manner:
|
(4.52) |
where the differential impulse
and step
response tensors are defined through the
following relations:
|
|
|
|
|
|
|
(4.53) |
In order to obtain the zero field mobility tensor it is necessary to integrate the differential velocity impulse response over a secondary trajectory
for a sufficiently long time. However, the time integration can be stopped after the first velocity randomizing scattering event has occurred, because in
this
case the correlation of the trajectory's initial velocity with the after-scattering velocity is lost. Since in the thermodynamic equilibrium the before
and
after-scattering distributions are equal, the secondary trajectories can be mapped onto the main trajectory. As a result the algorithm schematically
depicted in Fig. 4.6 is obtained.
Figure 4.6:
Zero field Monte Carlo flow chart.
|
- Set , .
- Select initial state arbitrarily.
- Compute a sum of weights:
.
- Select a free-flight time
and add time integral to estimator:
or use the expected value of the time integral:
.
- Perform scattering. If mechanism was isotropic, reset weight: .
- Continue with step 3 until -points have been generated.
- Calculate component of zero field mobility tensor as
.
Especially the diagonal elements can be calculated very efficiently using this algorithm. Consider a system where only isotropic scattering events take
place. Then the product is always positive, independent of the sign of . Therefore, only positive values are added to the estimator, which
leads to a low variance.
S. Smirnov: