Message Boards Message Boards


Fluorescence lifetime

Posted 1 year ago
5 Replies
0 Total Likes

Hello, I am new to Mathematica software and I am having problem to start my coding in Mathematica. to be honest, a bit clueless, I want to do simulation for fluorescence lifetime measurement. Here is my equation,

my main equation

Here, ? is the initial fluorescence intensity emitted by the sample when irradiated by laser pulses; ? is the fluorescence lifetime; ? denotes the repetition angular frequency of the short-pulsed laser; T is the period of pulse repetition (T = 2?/?); ? is the phase difference between the excitation pulse sequence and reference signal of the lock-in amplifier. Three phase-difference values of 0, ?/4, and ?/2 were considered.

Next, I want to do

enter image description here

where multiple fluorescence lifetime is used.

Can anyone help me how to do simulation for equation 1 and equation 2? Thank you very much!

5 Replies
Posted 1 year ago

Do these do what you want?

i[tau_,phi_]:=Integrate[alpha E^(-t/tau)Cos[omega t+phi],{t,0,T}];


i[tau1_,tau2_,phi_]:=Integrate[(alpha1 E^(-t/tau1)+alpha2 E^(-t/tau2))Cos[omega t+phi],{t,0,T}]
Posted 1 year ago

thank you so much for replying my question. may i know what is the i[1,0] i[2,Pi/4] for?

Posted 1 year ago

I apologize if I misunderstood or wasn't clear.

In each of those the first line defines a function to do your integration and the second and third lines show example usage.

Asking Mathematica to evaluate this

i[tau_,phi_]:=Integrate[alpha E^(-t/tau)Cos[omega t+phi],{t,0,T}];

then it will perform your integration with those two parameter values and should immediately display this

(alpha*(E^T - Cos[omega*T] + omega*Sin[omega*T]))/(E^T*(1 + omega^2))

Likewise asking Mathematica to evaluate

i[tau1_,tau2_,phi_]:=Integrate[(alpha1 E^(-t/tau1)+alpha2 E^(-t/tau2))Cos[omega t+phi],{t,0,T}];

will perform your integration with those three parameter values and should immediately display

((alpha2*E^T*(-(E^T*omega) + omega*Cos[omega*T] + Sin[omega*T]))/(1 + omega^2) + 
 (alpha1*(-(E^(2*T)*omega) + omega*Cos[omega*T] + 2*Sin[omega*T]))/(4 + omega^2))/E^(2*T)

If I still don't have this correct then please let me know and we will try to resolve this.

Posted 1 year ago

thank you so much for your explanation, I really appreciate it.

if i want to change the value of phi to 0. and then change to to phi = pi/2 and lastly phi = pi/4, i should change the value phi with i[1,0] i[1,pi/2] and i[1,pi/2]?

also, how can i display in on a graph? enter image description here

Posted 1 year ago

A few comments.

If you define and evaluate

i[tau_,phi_]:=Integrate[alpha E^(-t/tau)Cos[2*pi*freq*t*phi),{t,0,T}]

then Mathematica will remember that until you define something it to be something different or until you exit Mathematica. That means you do not need to type that line again.

Defining a constant, like


will be remembered the same way until you change it or exit Mathematica.

The symbol Pi in Mathematica is exactly the value of pi.

If you put a ; at the end of a line it will calculate but not display the result.


i[tau_,phi_]:=Integrate[alpha E^(-t/tau)Cos[2*Pi*freq*t+phi],{t,0,T}];

will calculate the integral and then plot the function of T.



will do the same plot, but will calculate more points to produce a more detailed plot.

Plot will not plot expressions that contain complex numbers, even if the value is very small, and will not plot expressions that have other variables which have not been assigned a value.

Please see if this is sufficient to help you get what you are looking for.

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
or Discard

Group Abstract Group Abstract