Semiconductor Bloch equations in Wannier gauge with well-behaved dephasing

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Thümmler, Martin, Lettau, Thomas, Croy, Alexander, Peschel, Ulf, Gräfe, Stefanie
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909939359809536
author Thümmler, Martin
Lettau, Thomas
Croy, Alexander
Peschel, Ulf
Gräfe, Stefanie
author_facet Thümmler, Martin
Lettau, Thomas
Croy, Alexander
Peschel, Ulf
Gräfe, Stefanie
contents The semiconductor Bloch equations (SBEs) with a dephasing operator for the microscopic polarizations are a well established approach to simulate high-harmonic spectra in solids. We discuss the impact of the dephasing operator on the stability of the numerical integration of the SBEs in the Wannier gauge. It is shown that the standard approach to apply dephasing is ill-defined in the presence of band crossings and leads to artifacts in the carrier distribution. They are caused by rapid changes of the dephasing operator matrix elements in the Wannier gauge, which render the convergence of the simulation in the stationary basis infeasible. In the comoving basis, also called Houston basis, these rapid changes can be resolved, but only at the cost of a largely increased computation time. As a remedy, we propose a modification of the dephasing operator with reduced magnitude in energetically close subspaces. This approach removes the artifacts in the carrier distribution and significantly speeds up the calculations, while affecting the high-harmonic spectrum only marginally. To foster further development, we provide our parallelized source code.
format Preprint
id arxiv_https___arxiv_org_abs_2508_07869
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Semiconductor Bloch equations in Wannier gauge with well-behaved dephasing
Thümmler, Martin
Lettau, Thomas
Croy, Alexander
Peschel, Ulf
Gräfe, Stefanie
Optics
The semiconductor Bloch equations (SBEs) with a dephasing operator for the microscopic polarizations are a well established approach to simulate high-harmonic spectra in solids. We discuss the impact of the dephasing operator on the stability of the numerical integration of the SBEs in the Wannier gauge. It is shown that the standard approach to apply dephasing is ill-defined in the presence of band crossings and leads to artifacts in the carrier distribution. They are caused by rapid changes of the dephasing operator matrix elements in the Wannier gauge, which render the convergence of the simulation in the stationary basis infeasible. In the comoving basis, also called Houston basis, these rapid changes can be resolved, but only at the cost of a largely increased computation time. As a remedy, we propose a modification of the dephasing operator with reduced magnitude in energetically close subspaces. This approach removes the artifacts in the carrier distribution and significantly speeds up the calculations, while affecting the high-harmonic spectrum only marginally. To foster further development, we provide our parallelized source code.
title Semiconductor Bloch equations in Wannier gauge with well-behaved dephasing
topic Optics
url https://arxiv.org/abs/2508.07869