FDTD 2D optical trapping

Hi

I am doing a computer program to simulate 2D trapping. Now i have some problems simulating the Gaussian source. Can anyone help me?

Thanks

As all our tweezer collaborators are now on their Christmas break, we have put together some useful links...

You may be interested in joining this group -

SciLife | Optical Tweezers - to help you further your work.

For the latest news -

Blog on Optical Tweezers

- also on Facebook and Google+ along with a number of other tweezer users.

Here's a Java applet demonstrating Tweezer force -

optical tweezers at glasgow

- and here's a Java program program that illustrates how optical tweezers work and has variables you can tweak -

optical-tweezers - Java simulation of optical trapping - Google Project Hosting

What university are you studying at?

Gaussian source can be programmed with following code. Does this apply to your needs?(Taken from 2D FDTD prog by Doug Neubauer)

pi = (acos(-1.0));
cc = 2.99792458e8; //speed of light in free space (meters/second)
muz = 4.0 * pi * 1.0e-7; //permeability of free space
epsz = 1.0 / (cc * cc * muz); //permittivity of free space

freq = 5.0e+9; //center frequency of source excitation (Hz)
lambda = cc / freq; //center wavelength of source excitation
omega = 2.0 * pi * freq; //center frequency in radians

dx = 3.0e-3; //space increment of square lattice (meters)
dt = dx / (2.0 * cc); //time step, seconds, courant limit, Taflove1995 page 177

rtau = 160.0e-12;
tau = rtau / dt;
delay = 3 * tau;
for (i = 0; i < nmax; i++) {
source[i] = 0.0;
} /* iForLoop */

for (n = 0; n < (int )(7.0 * tau); n++) {
temporary = (double )n - delay;
source[n] = sin( omega * (temporary) * dt) * exp(-( (temporary * temporary)/(tau * tau) ) );
} /* forLoop */

What exactly is the problem? The definition of a Gaussian pulse/Inserting the pulse/...?