Cprime-params

Cprime-params [<dk>=1E-4] [<degeneracyThreshold>=1E-6] [<vThreshold>=1E-4] [<realSpaceTruncated>=yes]

Control dC/dk calculation for L and Q output. Here, <dk> in a0^-1 controls the finite difference used for dC/dk, while <degeneracyThreshold> specifies the energy range within unitary rotations are accounted for in comparing wavefunctions between k. Within degenerate subspaces of energy, rotations are first revolved by the velocity operator and then obtained by best match between wavefunctions for sub-subspaces that have the same velocity within <vThreshold> (in Eh-a0 atomic units). If <realSpaceTruncated> is yes (default), then truncated directions are computed by a real space multiplication by r, instead of dC/dk.

This is used for the calculation of orbital angular momenta L, output in the same binary format as the momenta and selected by 'dump End L'.

Additionally, 'dump End Q' selects output of electric quadrupole matrix elements, defined as traceless symmetric tensor of (rj pk + pj rk). Specifically, the output is a nStates x 5 x nBands x nBands complex binary file, where the 5 components in order are xy, yz, zx, xxr, yyr (where xxr = xx - rr/3). Note that zzr = -(xxr + yyr) because the trace of the tensor is removed, and this redundant component is excluded.

**Requires:** (None)

**Forbids:** (None)

**Allow multiple:** no

**Default:** (None)

Back to: Input file documentation or Index of commands