Fix saving lists of arrays with recent versions of numpy
[qpms.git] / notes / kambe_linton_dict.lyx
blob384064368473d0d2fd05c71bc239b8c488eb1021
1 #LyX 2.4 created this file. For more info see https://www.lyx.org/
2 \lyxformat 584
3 \begin_document
4 \begin_header
5 \save_transient_properties true
6 \origin unavailable
7 \textclass article
8 \use_default_options true
9 \maintain_unincluded_children false
10 \language finnish
11 \language_package default
12 \inputencoding utf8
13 \fontencoding auto
14 \font_roman "default" "default"
15 \font_sans "default" "default"
16 \font_typewriter "default" "default"
17 \font_math "auto" "auto"
18 \font_default_family default
19 \use_non_tex_fonts false
20 \font_sc false
21 \font_roman_osf false
22 \font_sans_osf false
23 \font_typewriter_osf false
24 \font_sf_scale 100 100
25 \font_tt_scale 100 100
26 \use_microtype false
27 \use_dash_ligatures true
28 \graphics default
29 \default_output_format default
30 \output_sync 0
31 \bibtex_command default
32 \index_command default
33 \paperfontsize default
34 \use_hyperref false
35 \papersize default
36 \use_geometry false
37 \use_package amsmath 1
38 \use_package amssymb 1
39 \use_package cancel 1
40 \use_package esint 1
41 \use_package mathdots 1
42 \use_package mathtools 1
43 \use_package mhchem 1
44 \use_package stackrel 1
45 \use_package stmaryrd 1
46 \use_package undertilde 1
47 \cite_engine basic
48 \cite_engine_type default
49 \use_bibtopic false
50 \use_indices false
51 \paperorientation portrait
52 \suppress_date false
53 \justification true
54 \use_refstyle 1
55 \use_minted 0
56 \use_lineno 0
57 \index Index
58 \shortcut idx
59 \color #008000
60 \end_index
61 \secnumdepth 3
62 \tocdepth 3
63 \paragraph_separation indent
64 \paragraph_indentation default
65 \is_math_indent 0
66 \math_numbering_side default
67 \quotes_style english
68 \dynamic_quotes 0
69 \papercolumns 1
70 \papersides 1
71 \paperpagestyle default
72 \tablestyle default
73 \tracking_changes false
74 \output_changes false
75 \html_math_output 0
76 \html_css_as_file 0
77 \html_be_strict false
78 \end_header
80 \begin_body
82 \begin_layout Standard
84 \lang english
85 \begin_inset FormulaMacro
86 \newcommand{\vect}[1]{\mathbf{#1}}
87 \end_inset
90 \lang finnish
92 \begin_inset FormulaMacro
93 \newcommand{\Kambe}[1]{#1^{\mathrm{K}}}
94 \end_inset
97 \begin_inset FormulaMacro
98 \newcommand{\Linton}[1]{#1^{\mathrm{L}}}
99 \end_inset
102 \end_layout
104 \begin_layout Standard
105 Here and in Kambe's papers, 
106 \begin_inset Formula $\kappa$
107 \end_inset
109  is the wavenumber (
110 \begin_inset Formula $k$
111 \end_inset
113  in Linton).
114  Here 
115 \begin_inset Formula $\vect K_{p}$
116 \end_inset
118  is a point of the reciprocal lattice (
119 \begin_inset Formula $\vect K_{p}=\Kambe{\vect K_{pt}}=\Linton{\vect{\beta}_{\mu}}$
120 \end_inset
123 \end_layout
125 \begin_layout Section
126 \begin_inset Quotes eld
127 \end_inset
129 Gammas
130 \begin_inset Quotes erd
131 \end_inset
134 \end_layout
136 \begin_layout Standard
137 For 
138 \begin_inset Formula $\kappa$
139 \end_inset
141  positive,
142 \end_layout
144 \begin_layout Standard
145 \begin_inset Formula 
147 \Kambe{\Gamma_{p}}\equiv\begin{cases}
148 \sqrt{\kappa^{2}-\left|\vect K_{p}\right|^{2}} & \kappa^{2}-\left|\vect K_{p}\right|^{2}>0\\
149 i\sqrt{\left|\vect K_{p}\right|^{2}-\kappa^{2}} & \kappa^{2}-\left|\vect K_{p}\right|^{2}<0
150 \end{cases}
153 \end_inset
156 \begin_inset Formula 
158 \Linton{\gamma_{\mu}}\equiv\begin{cases}
159 \sqrt{\left(\frac{\vect K_{p}}{\kappa}\right)^{2}-1} & \kappa-\left|\vect K_{p}\right|\le0\\
160 -i\sqrt{1-\left(\frac{\vect K_{p}}{\kappa}\right)^{2}} & \kappa-\left|\vect K_{p}\right|>0
161 \end{cases}
164 \end_inset
166 hence 
167 \begin_inset Formula 
169 \Kambe{\Gamma_{p}}=-i\kappa\Linton{\gamma_{\mu}},
172 \end_inset
175 \end_layout
177 \begin_layout Standard
178 \begin_inset Formula 
180 \Linton{\gamma_{\mu}}=i\frac{\Kambe{\Gamma_{p}}}{\kappa}.
183 \end_inset
186 \end_layout
188 \begin_layout Section
189 D vs sigma
190 \end_layout
192 \begin_layout Standard
193 In-plane sums [Linton 2009, (4.5)], replacing 
194 \begin_inset Formula $n,m\rightarrow L,M$
195 \end_inset
198 \begin_inset Formula $k\rightarrow\kappa$
199 \end_inset
202 \end_layout
204 \begin_layout Standard
206 \lang english
207 \begin_inset Formula 
208 \begin{eqnarray*}
209 \sigma_{L}^{M(1)} & = & -\frac{i^{L+1}}{2\kappa^{2}\mathscr{A}}\left(-1\right)^{\left(L+M\right)/2}\sqrt{\left(2L+1\right)\left(L-M\right)!\left(L+M\right)!}\times\\
210  &  & \times\sum_{\vect K_{pq}\in\Lambda^{*}}^{'}\sum_{j=0}^{\left[\left(L-\left|M\right|/2\right)\right]}\frac{\left(-1\right)^{j}\left(\beta_{pq}/2\kappa\right)^{L-2j}e^{iM\phi_{\vect{\beta}_{pq}}}\Gamma_{j,pq}}{j!\left(\frac{1}{2}\left(L-M\right)-j\right)!\left(\frac{1}{2}\left(L+M\right)-j\right)!}\left(\frac{\gamma_{pq}}{2}\right)^{2j-1}
211 \end{eqnarray*}
213 \end_inset
215 [Kambe II, (3.17)], replacing 
216 \lang finnish
218 \begin_inset Formula $n\rightarrow j$
219 \end_inset
222 \lang english
224 \lang finnish
226 \begin_inset Formula $A\rightarrow\mathscr{A}$
227 \end_inset
230 \begin_inset Formula $\vect K_{pt}\to\vect K_{p}$
231 \end_inset
234 \begin_inset Formula $\Gamma\left(\frac{1}{2}-j,e^{-i\pi}\Gamma_{p}^{2}\omega/2\right)\to\Gamma_{j,p}$
235 \end_inset
237  and performing little typographic modifications
238 \lang english
240 \begin_inset Formula 
241 \begin{align*}
242 D_{LM} & =-\frac{1}{\mathscr{A}\kappa}i^{\left|M\right|+1}2^{-L}\sqrt{\left(2L+1\right)\left(L+\left|M\right|\right)!\left(L-\left|M\right|\right)!}\times\\
243  & \quad\times\sum_{p}e^{i\vect K_{p}\cdot\vect c_{ijt}}e^{-iM\phi_{K_{p}}}\sum_{j=0}^{\left(L-\left|M\right|\right)/2}\frac{\left(\Gamma_{p}/\kappa\right)^{2j-1}\left(K_{p}/\kappa\right)^{L-2j}\Gamma_{j,p}}{j!\left(\frac{1}{2}\left(L-\left|M\right|\right)-j\right)!\left(\frac{1}{2}\left(L+\left|M\right|\right)-j\right)!}
244 \end{align*}
246 \end_inset
248 Using the relations between 
249 \begin_inset Formula $\Kambe{\Gamma_{p}}=-i\kappa\Linton{\gamma_{\mu}}$
250 \end_inset
252 , we have (also, we replace the 
253 \begin_inset Formula $\mu$
254 \end_inset
256  index with 
257 \begin_inset Formula $p$
258 \end_inset
261 \begin_inset Formula 
262 \begin{align*}
263 D_{LM} & =-\frac{1}{\mathscr{A}\kappa}i^{\left|M\right|+1}2^{-L}\sqrt{\left(2L+1\right)\left(L+\left|M\right|\right)!\left(L-\left|M\right|\right)!}\times\\
264  & \quad\times\sum_{p}e^{i\vect K_{p}\cdot\vect c_{ijt}}e^{-iM\phi_{K_{p}}}\sum_{j=0}^{\left(L-\left|M\right|\right)/2}\frac{\left(-i\gamma_{p}\right)^{2j-1}\left(K_{p}/\kappa\right)^{L-2j}\Gamma_{j,p}}{j!\left(\frac{1}{2}\left(L-\left|M\right|\right)-j\right)!\left(\frac{1}{2}\left(L+\left|M\right|\right)-j\right)!}
265 \end{align*}
267 \end_inset
269 and now, trying to make the exponents look the same as in Linton, 
270 \begin_inset Formula $2^{-1}2^{2j-L}2^{1-2j}=2^{-L}$
271 \end_inset
273  (OK), 
274 \begin_inset Formula $K_{p}^{L-2j}=K_{p}^{L-2j}$
275 \end_inset
277  (OK), 
278 \begin_inset Formula 
279 \begin{align*}
280 D_{LM} & =-\frac{1}{2\kappa\mathscr{A}}i^{\left|M\right|+1}\sqrt{\left(2L+1\right)\left(L+\left|M\right|\right)!\left(L-\left|M\right|\right)!}\times\\
281  & \quad\times\sum_{p}e^{i\vect K_{p}\cdot\vect c_{ij}}e^{-iM\phi_{K_{p}}}\sum_{j=0}^{\left(L-\left|M\right|\right)/2}\frac{\left(-i\right)^{2j-1}\left(K_{p}/2\kappa\right)^{L-2j}\Gamma_{j,p}}{j!\left(\frac{1}{2}\left(L-\left|M\right|\right)-j\right)!\left(\frac{1}{2}\left(L+\left|M\right|\right)-j\right)!}\left(\frac{\gamma_{p}}{2}\right)^{2j-1}
282 \end{align*}
284 \end_inset
286 There are now these differences left:
287 \end_layout
289 \begin_layout Itemize
291 \lang english
292 Additional 
293 \begin_inset Formula $\kappa$
294 \end_inset
296  factor in 
297 \begin_inset Formula $D_{LM}$
298 \end_inset
301 \end_layout
303 \begin_layout Itemize
305 \lang english
306 \begin_inset Formula $i^{L+1}\left(-1\right)^{\left(L+M\right)/2}\left(-1\right)^{j}$
307 \end_inset
309  vs.
311 \begin_inset Formula $i^{\left|M\right|+1}\left(-i\right)^{2j-1}$
312 \end_inset
315 \end_layout
317 \begin_layout Itemize
319 \lang english
320 Opposite phase in the angular part.
321 \end_layout
323 \begin_layout Itemize
325 \lang english
326 Plane wave factor in 
327 \begin_inset Formula $D_{LM}$
328 \end_inset
331 \end_layout
333 \begin_layout Standard
335 \lang english
336 Let's look at the 
337 \begin_inset Formula $i,-1$
338 \end_inset
340  factors (note that 
341 \begin_inset Formula $L+M$
342 \end_inset
344  is odd): 
345 \begin_inset Formula $\left(-i\right)^{2j}=\left(-1\right)^{j},$
346 \end_inset
348  leaving 
349 \begin_inset Formula $i^{L+1}\left(-1\right)^{\left(L+M\right)/2}$
350 \end_inset
352  vs.
354 \begin_inset Formula $i^{\left|M\right|+1}i$
355 \end_inset
358  So there is might be a phase difference due to different conventions, but
359  it does not depend on 
360 \begin_inset Formula $j$
361 \end_inset
363 , so one should be able to transplant the 
364 \begin_inset Formula $z\ne0$
365 \end_inset
367  sum from Kambe without major problems.
368 \end_layout
370 \begin_layout Section
371 Ewald parameter (integration limits)
372 \end_layout
374 \begin_layout Standard
375 \begin_inset Formula 
377 \Linton{\eta}=\sqrt{\frac{1}{2\Kambe{\omega}}}
380 \end_inset
382 (Based on comparison of some function arguments, not checked.)
383 \end_layout
385 \end_body
386 \end_document