Skip to content

Commit 53b47a2

Browse files
committed
correct sign of potential definition as suggested by @BeyondEspresso
1 parent cd2c917 commit 53b47a2

File tree

1 file changed

+4
-4
lines changed

1 file changed

+4
-4
lines changed

src/picongpu/include/fields/laserProfiles/laserPlaneWave.hpp

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -31,19 +31,19 @@ namespace picongpu
3131
* no transverse spacial envelope
3232
* based on the electric potential
3333
* Phi = Phi_0 * exp(0.5 * (x-x_0)^2 / sigma^2) * cos(k*(x - x_0) - phi)
34-
* by applying grad Phi = d/dx Phi = E(x)
34+
* by applying -grad Phi = -d/dx Phi = E(x)
3535
* we get:
36-
* 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)]
36+
* 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)]
3737
*
3838
* This approach ensures that int_{-infinity}^{+infinity} E(x) = 0 for any phase
3939
* if we have no transverse profile as we have with this plane wave train
4040
*
4141
* Since PIConGPU requires a temporally defined electric field, we use:
4242
* 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
4343
* and get:
44-
* 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)]
44+
* 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)]
4545
* and define:
46-
* E_0 = Phi_0*omega/c
46+
* E_0 = -Phi_0*omega/c
4747
* integrationCorrectionFactor = t/(omega*tau^2)
4848
*
4949
* Please consider:

0 commit comments

Comments
 (0)