@@ -30,16 +30,26 @@ namespace picongpu
30
30
{
31
31
/* * plane wave (use periodic boundaries!)
32
32
*
33
- * no transversal spacial envelope
33
+ * no transverse spacial envelope
34
34
* based on the electric potential
35
- * Phi = E_0 * exp(0.5 * (t-t_0)^2 / tau^2) * cos(t - t_0 - phi)
36
- * by applying t = x/c, the spatial derivative can be interchanged by the temporal derivative
37
- * resulting in:
38
- * E = E_0 * exp(...) * [sin(...) + t/tau^2 * cos(...)]
39
- * This ensures int_{-infinty}^{+infinty} E(x) = 0 for any phase.
35
+ * Phi = Phi_0 * exp(0.5 * (x-x_0)^2 / sigma^2) * cos(k*(x - x_0) - phi)
36
+ * by applying grad Phi = d/dx Phi = E(x)
37
+ * we get:
38
+ * E = Phi_0 * exp(0.5 * (x-x_0)^2 / sigma^2) * [k*sin(k*(x - x_0) - phi) + x/sigma^2 * cos(k*(x - x_0) - phi)]
40
39
*
41
- * The plateau length needs to be set to a multiple of the wavelength,
42
- * otherwise the integral will not vanish.
40
+ * This approach ensures that int_{-infinity}^{+infinity} E(x) = 0 for any phase
41
+ * if we have no transverse profile as we have with this plane wave train
42
+ *
43
+ * Since PIConGPU requires a temporally defined electric field, we use:
44
+ * t = x/c and (x-x_0)/sigma = (t-t_0)/tau and k*(x-x_0) = omega*(t-t_0) with omega/k = c and tau * c = sigma
45
+ * and get:
46
+ * E = Phi_0*omega/c * exp(0.5 * (t-t_0)^2 / tau^2) * [sin(omega*(t - t_0) - phi) + t/(omega*tau^2) * cos(omega*(t - t_0) - phi)]
47
+ * and define Phi_0*omega/c = E_0
48
+ *
49
+ * Please consider:
50
+ * The above formulae does only apply to a Gaussian envelope. If the plateau length is
51
+ * not zero, the integral over the volume will only vanish if the plateau length is
52
+ * a multiple of the wavelength,
43
53
*/
44
54
namespace laserPlaneWave
45
55
{
@@ -110,7 +120,7 @@ namespace picongpu
110
120
return elong;
111
121
}
112
122
113
- /* * calculates transversal field distribution
123
+ /* * calculates transverse field distribution
114
124
*
115
125
* @param elong
116
126
* @param phase
0 commit comments