dipole_matrix Interface

public interface dipole_matrix

Form either a real or a complex dipole matrix.

The real-typed version is equivalent to $\mathbf q=0$.

$$\boldsymbol{\mathcal A}:=[\mathbf a_1\mathbf a_2\mathbf a_3],\qquad\boldsymbol{\mathcal B}:=[\mathbf b_1\mathbf b_2\mathbf b_3] \\ \boldsymbol{\mathcal B}=2\pi(\boldsymbol{\mathcal A}^{-1})^\mathrm T,\qquad \partial\boldsymbol{\mathcal B}=-\big((\partial\boldsymbol{\mathcal A})\boldsymbol{\mathcal A}^{-1}\big)^\mathrm T\boldsymbol{\mathcal B} \\ \mathbf R_\mathbf n=\boldsymbol{\mathcal A}\mathbf n,\qquad\partial\mathbf R_\mathbf n=(\partial\boldsymbol{\mathcal A})\mathbf n, \\ \mathbf G_\mathbf m=\boldsymbol{\mathcal B}\mathbf m,\qquad \partial\mathbf G_\mathbf m=-\big((\partial\boldsymbol{\mathcal A})\boldsymbol{\mathcal A}^{-1}\big)^\mathrm T\mathbf G_\mathbf m, \\ \frac{\partial G_{\mathbf ma}}{\partial A_{bc}}=-\mathcal A^{-1}_{ca}G_{\mathbf mb}$$

$$\begin{gathered} \mathbf T_{ij}(\mathbf q)=\sum_{\mathbf n}\mathbf T(\mathbf R_{\mathbf nij})\mathrm e^{-\mathrm i\mathbf q\cdot\mathbf R_{\mathbf nij}},\quad\mathbf R_{\mathbf nij}=\mathbf R_j+\mathbf R_\mathbf n-\mathbf R_i \\ \frac{\mathrm d\mathbf R_{\mathbf nij}}{\mathrm d\mathbf R_k}=(\delta_{jk}-\delta_{ik})\mathbf I \\ \mathbf{T}_{ij}(\mathbf{q})\approx\mathbf{T}^\text{Ew}_{ij}(\mathbf{q}) =\sum_\mathbf n^{|\mathbf R_{\mathbf nij}|<R_\text c}\mathbf T^\text{erfc}(\mathbf R_{\mathbf nij};\gamma)\mathrm e^{-\mathrm i\mathbf q\cdot\mathbf R_{\mathbf nij}} +\frac{4\pi}{V_\text{uc}}\sum_{\mathbf m}^{0<|\mathbf k_\mathbf m|<k_\text c}\mathbf{\hat k}_\mathbf m\otimes\mathbf{\hat k}_\mathbf m\,\mathrm e^{-\frac{k_\mathbf m^2}{4\gamma^2}-\mathrm i\mathbf G_\mathbf m\cdot\mathbf R_{ij}} \\ -\frac{4\gamma^3}{3\sqrt\pi}\delta_{ij}\mathbf I +\delta(\mathbf q)\frac{4 \pi}{3 V_\text{uc}}\mathbf I,\qquad \mathbf k_\mathbf m=\mathbf G_\mathbf m+\mathbf q \end{gathered}$$

$$\partial\left(\frac{4\pi}{V_\text{uc}}\right)=-\frac{4\pi}{V_\text{uc}}\frac{\partial V_\text{uc}}{V_\text{uc}},\qquad\frac{\partial V_\text{uc}}{\partial\boldsymbol{\mathcal A}}=V_\text{uc}(\boldsymbol{\mathcal A}^{-1})^\mathrm T \\ \partial(k^2)=2\mathbf k\cdot\partial\mathbf k \\ \mathbf{\hat k}\otimes\partial\mathbf{\hat k}=\frac{\mathbf k\otimes\partial\mathbf k}{k^2}-\frac{\mathbf k\otimes\mathbf k}{2k^4}\partial(k^2)$$

$$\gamma:=\frac{2.5}{\sqrt[3]{V_\text{uc}}},\quad R_\text c:=\frac6\gamma,\quad k_\text c:=10\gamma$$

Contents

Module Procedures

dipole_matrix_real dipole_matrix_complex

Module Procedures