| ... | @@ -59,9 +59,13 @@ Static codes allow for defining your own renormalization scheme. You need to pro |
... | @@ -59,9 +59,13 @@ Static codes allow for defining your own renormalization scheme. You need to pro |
|
|
# Regularization scheme and the energy conservation in td calculations
|
|
# Regularization scheme and the energy conservation in td calculations
|
|
|
In publication [arXiv:1606.02225](https://arxiv.org/abs/1606.02225) it was pointed that TDBdG like equations, formally conserve energy only if all quasiparticle states are evolved, see discussion of Eqs.(25)-(26). This situation corresponds to the cubic cutoff. If the space is truncated eg. by introducing a spherical cutoff at some initial time then in general energy maybe not conserved.
|
|
In publication [arXiv:1606.02225](https://arxiv.org/abs/1606.02225) it was pointed that TDBdG like equations, formally conserve energy only if all quasiparticle states are evolved, see discussion of Eqs.(25)-(26). This situation corresponds to the cubic cutoff. If the space is truncated eg. by introducing a spherical cutoff at some initial time then in general energy maybe not conserved.
|
|
|
|
|
|
|
|
In practical applications, we observe that the energy when applied the spherical regularization scheme is conserved only with some accuracy, which is not related to the integrator accuracy. Below we provide an example of (3d calculation), where for the time interval $`te_F<170`$ we apply an external time-dependent potential (we pump energy into the system), and for $`te_F>170`$ the system evolves without any external perturbation.
|
|
In practical applications, we observe that the energy when applied the spherical regularization scheme is conserved only with some accuracy, which is not related to the integrator accuracy. Below we provide an example of (3d calculation), where for the time interval $`te_F<170`$ we apply an external time-dependent potential (we pump energy into the system), and for $`te_F>170`$ the system evolves without any external perturbation.
|
|
|
|
|
|
|
|
|
|
|
|
|
It is clearly visible, that for evolution with the cubic cutoff the energy is conserved up to high accuracy, while for spherical cutoff the quality of the energy conservation is significantly lower.
|
|
It is clearly visible, that for evolution with the cubic cutoff the energy is conserved up to high accuracy, while for spherical cutoff the quality of the energy conservation is significantly lower.
|
|
|
|
|
|
|
|
|
`wlog` files for these runs:
|
|
|
|
* [cubic.wlog](uploads/5aff56cbbd8ee97046cc1c3d10a49867/cubic.wlog)
|
|
|
|
* [spherical.wlog](uploads/e642fe76e9496d2305f4c462ebc61e72/spherical.wlog)
|
|
|
|
|
|
|
In conclusion, we find that typically for trajectories of length $`te_F\approx1000`$ the spherical cutoff provides reasonable accuracy, while for generation of long trajectories $`te_F\gg 1000`$ it is recommended to use the cubic cutoff. |
|
In conclusion, we find that typically for trajectories of length $`te_F\approx1000`$ the spherical cutoff provides reasonable accuracy, while for generation of long trajectories $`te_F\gg 1000`$ it is recommended to use the cubic cutoff. |
|
|
|
\ No newline at end of file |
