Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Fluorescent material in Wave Optics

Please login with a confirmed email address before reporting spam

Hi all,

For my master thesis research I'm looking into the possibilities of incorporating a luminescent material into an anti-reflective coating. Creating an non-fluorescent AR coating has worked for me with the Wave Optics-Beam envelopes module, however I am quite lost as to how to make one of the AR layers material now flourescent? The internet is also of little help on working with fluorescence in COMSOL, there is little done before it seems. I thought of implementing light-emitting point sources randomly in the layer, however I don't find anything how this can be done in the Wave Optics module. Since point source create singularities in the Wave Optics module. Another idea I had is defining the layer as a sort of 'light source' on its own, but it seems not possible in the Wave Optics module as well.

Anyone with the know-how on how fluorescence can be achieved best?

Many thanks in advance!!


14 Replies Last Post Jul 1, 2020, 2:15 p.m. EDT

Please login with a confirmed email address before reporting spam

Posted: 5 years ago May 19, 2020, 8:13 a.m. EDT
Updated: 5 years ago May 19, 2020, 8:14 a.m. EDT

Hi Karima,

I've never worked with fluorescent materials, but this sounds really interesting! Both of your ideas seems good. Did you had any success applying them on your model?

Another approach that may help you, is to add an imaginary part to the refractive index of the fluorescent material. Normally, this imaginary part accounts for the absorption losses of the material. But, you can switch the sign so it accounts for gain instead of losses. The problem is, this approach will amplify every wavelength impinging on your fluorescent material. But there are tweaks to solve this. So please tell me, which kind of source are you use (or need) to excite the fluorescence?

Best!

-------------------
Felipe Beltran-Mejia
Hi Karima, I've never worked with fluorescent materials, but this sounds really interesting! Both of your ideas seems good. Did you had any success applying them on your model? Another approach that may help you, is to add an imaginary part to the refractive index of the fluorescent material. Normally, this imaginary part accounts for the absorption losses of the material. But, you can switch the sign so it accounts for gain instead of losses. The problem is, this approach will amplify every wavelength impinging on your fluorescent material. But there are tweaks to solve this. So please tell me, which kind of source are you use (or need) to excite the fluorescence? Best!

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 4:47 a.m. EDT

Hii Felipe,

Sorry for this late response! I was having some difficulties with COMSOL and WaveOptics, I did manage to add a complex refractive index to my material. Which did as you predicted, all wavelengths were amplified in this manner. However I also want the material to absorb one part of the spectrum (UV light) and emit the light in another part (red or blue), the so-called Stokes Shift. The conceptual idea is to excite the material with sunlight and have the fluorescent material absorb the UV and re-emit it in the visible region, preferably blue or red. Currently I use a port condition as exciting port and in the pre-defined Wavelength Domain module of the WaveOptics, Frequency domain I use a built in sweep in wavelengths to imitate the solar spectrum, ranging it from 300nm to 700 nm. Hope this paints the picture better of what I'm trying to achieve with COMSOL.

What were your ideas on tweaking this problem with the imaginary k of the refractive index? I'm very curious!!

Hii Felipe, Sorry for this late response! I was having some difficulties with COMSOL and WaveOptics, I did manage to add a complex refractive index to my material. Which did as you predicted, all wavelengths were amplified in this manner. However I also want the material to absorb one part of the spectrum (UV light) and emit the light in another part (red or blue), the so-called Stokes Shift. The conceptual idea is to excite the material with sunlight and have the fluorescent material absorb the UV and re-emit it in the visible region, preferably blue or red. Currently I use a port condition as exciting port and in the pre-defined Wavelength Domain module of the WaveOptics, Frequency domain I use a built in sweep in wavelengths to imitate the solar spectrum, ranging it from 300nm to 700 nm. Hope this paints the picture better of what I'm trying to achieve with COMSOL. What were your ideas on tweaking this problem with the imaginary k of the refractive index? I'm very curious!!

Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 5:02 a.m. EDT

Karima,

you can make the index a function of the wavelength. Would that help?

Cheers Edgar

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Karima, you can make the index a function of the wavelength. Would that help? Cheers Edgar

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 5:32 a.m. EDT

Hi Edgar,

Thank you for your response! I did actually make my refractive index wavelength dependent currently in the non-fluorescent material using experimental data from my research group. The main issue I think for me for the fluorescent material is how to incorporate the fact that the absorption takes place in a different part of the spectrum than the emission.

Hi Edgar, Thank you for your response! I did actually make my refractive index wavelength dependent currently in the non-fluorescent material using experimental data from my research group. The main issue I think for me for the fluorescent material is how to incorporate the fact that the absorption takes place in a different part of the spectrum than the emission.

Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 7:18 a.m. EDT

Could this be done in two subsequent study steps? First step calculates absorption and second step calculates emission using a Stoke's-shifted index derived from the first step. It seems that the problem is one-way coupled only?

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Could this be done in two subsequent study steps? First step calculates absorption and second step calculates emission using a Stoke's-shifted index derived from the first step. It seems that the problem is one-way coupled only?

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 8:51 a.m. EDT
Updated: 4 years ago Jun 23, 2020, 8:52 a.m. EDT

Aha! That sounds like a bright idea Edgar, if I understood it correctly: I would create first a positive imaginary k (maybe as a file under materials-rfi?) which would absorb in the UV region and then create a negative imaginary k which would simulate emission in the red/blue region? However how do I tell each Study Step what k to use?

Aha! That sounds like a bright idea Edgar, if I understood it correctly: I would create first a positive imaginary k (maybe as a file under materials-rfi?) which would absorb in the UV region and then create a negative imaginary k which would simulate emission in the red/blue region? However how do I tell each Study Step what k to use?

Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 23, 2020, 9:35 a.m. EDT

There are different ways to achieve that. You can use a variable k, that may be defined by an interpolation function for absorption and emission in two separate variable nodes for the two ks. Same variable name in both nodes. In the study step you can then tick the 'Modify model configuration for study step' check box and activate/deactivate the respective variable node. You can also set up the physics node twice with suitable k-assignments and run one physics in one node, the other in the second node. There may be more options like using boolean switches, such as: absorptionk1(freq) + emissionk2(freq). You set the parameters absorption and emission to 0 or 1 as needed.

The optimal way may depend on other requirements in your model.

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
There are different ways to achieve that. You can use a variable k, that may be defined by an interpolation function for absorption and emission in two separate variable nodes for the two ks. Same variable name in both nodes. In the study step you can then tick the 'Modify model configuration for study step' check box and activate/deactivate the respective variable node. You can also set up the physics node twice with suitable k-assignments and run one physics in one node, the other in the second node. There may be more options like using boolean switches, such as: absorption*k1(freq) + emission*k2(freq). You set the parameters absorption and emission to 0 or 1 as needed. The optimal way may depend on other requirements in your model.

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 24, 2020, 12:22 p.m. EDT

Hi Edgar,

I tried out one of your ideas today, namely the one where I create two functions for different k's. However I stumbled upon a problem when I created the second study step, since this only allows me to turn on/off the nodes under the Wave Optics node. I think this is because I have defined my luminescent material as a geometrical layer in the section Geometry and specified its refractive index under Materials, both of these nodes (so geometry and materials) I am unable to (de)-activate when I select the 'Modify model configuration for study step' check box in the second study step.

Would you have any other suggestions for me to work around this issue?

Hi Edgar, I tried out one of your ideas today, namely the one where I create two functions for different k's. However I stumbled upon a problem when I created the second study step, since this only allows me to turn on/off the nodes under the Wave Optics node. I think this is because I have defined my luminescent material as a geometrical layer in the section Geometry and specified its refractive index under Materials, both of these nodes (so geometry and materials) I am unable to (de)-activate when I select the 'Modify model configuration for study step' check box in the second study step. Would you have any other suggestions for me to work around this issue?

Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 24, 2020, 2:03 p.m. EDT

Hi Karima,

you can also put the material property directly into the respective physics node. You may need to add a separate physics node for this layer. This will then overide the settings in the material

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Hi Karima, you can also put the material property directly into the respective physics node. You may need to add a separate physics node for this layer. This will then overide the settings in the material

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 25, 2020, 5:41 p.m. EDT

Hi Karima and Edgar,

I'm enjoying a lot reading this thread. Thanks for sharing!

Karima, can you please share the functions you are using with us? And also, can you tell us if these functions come from (i.e. theory or a experimental fit)?

Best

-------------------
Felipe Beltran-Mejia
Hi Karima and Edgar, I'm enjoying a lot reading this thread. Thanks for sharing! Karima, can you please share the functions you are using with us? And also, can you tell us if these functions come from (i.e. theory or a experimental fit)? Best

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 26, 2020, 4:40 p.m. EDT

Hii Felipe and Edgar,

Filipe, welcome back to the thread haha:) I am currently using simplified versions for the absorption/emission functions, my research group has experimental data of the actual refractive index (both n and k). But I first wanted to try with simple Gaussian distributed functions, its parameters (mu and sigma) are not completely made up: the mean is obtained from the experiments I mentioned before and the FWHM as well. The only difference is that the experimental data I will use at a later stage does not follow a perfect Gaussian distribution like I'm simulating now. I used an analytic function to define these k absorbing/emitting, which look like: .

@Edgar creating a second "EM Wave, Frequency Domain" node and only including the fluorescent layer as a domain did not work for me. Since COMSOL does not seems to understand what I mean when I put in mat3.rfi.k_emit as material property under the second physics node, it only understands mat3.rfi.ki_iso. I did define both material functions k_abs and k_emit but chose k_abs as the ki_iso under the material node, since this is needed for the first physics node. So it seems that when I did that, the expression mat3.rfi.k_emit became unaccessable. This made no sense to me and this is when I stopped.

Hii Felipe and Edgar, Filipe, welcome back to the thread haha:) I am currently using simplified versions for the absorption/emission functions, my research group has experimental data of the actual refractive index (both n and k). But I first wanted to try with simple Gaussian distributed functions, its parameters (mu and sigma) are not completely made up: the mean is obtained from the experiments I mentioned before and the FWHM as well. The only difference is that the experimental data I will use at a later stage does not follow a perfect Gaussian distribution like I'm simulating now. I used an analytic function to define these k absorbing/emitting, which look like: k_{abs/emit} = e^{-\frac{1}{2}(\frac{ewfd.lambda0-\mu_{abs/emit}}{\sigma})^2}. @Edgar creating a second "EM Wave, Frequency Domain" node and only including the fluorescent layer as a domain did not work for me. Since COMSOL does not seems to understand what I mean when I put in mat3.rfi.k_emit as material property under the second physics node, it only understands mat3.rfi.ki_iso. I did define both material functions k_abs and k_emit but chose k_abs as the ki_iso under the material node, since this is needed for the first physics node. So it seems that when I did that, the expression mat3.rfi.k_emit became unaccessable. This made no sense to me and this is when I stopped.

Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 28, 2020, 4:42 a.m. EDT

Karima,

not sure what is going wrong there. Things may be easier if you can post your model and the expressions you want to use on the respective domains and study steps.

Cheers Edgar

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Karima, not sure what is going wrong there. Things may be easier if you can post your model and the expressions you want to use on the respective domains and study steps. Cheers Edgar

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jun 29, 2020, 9:32 a.m. EDT
Updated: 4 years ago Jun 29, 2020, 9:33 a.m. EDT

Hi Edgar,

Of course, I attached the file below! Thank you for your help, it means a lot:)

I read up on some literature and I think that maybe using a k_emit is not the accurate way to mimic the actual physical process at hand. Since with this method the impinging light is amplified for a specific wavelength (or range of wavelengths), it will not return any output any if the incoming light does not contain that specified wavelength. I don't want to amplify incoming light, but actually 'create' it. So in theory, when a monochromatic beam with a wavelength in the UV region hits my material slab, it should be absorbed and translated to an emission in the visible spectrum. I don't think this can be achieved with a k_emission as a material property in a second physics node. I thought of maybe having the fluorescent material as a whole still absorb the UV light but then have COMSOL 'create' the emission light with an electric dipole point, do you think this would be possible in COMSOL?

Hi Edgar, Of course, I attached the file below! Thank you for your help, it means a lot:) I read up on some literature and I think that maybe using a k_emit is not the accurate way to mimic the actual physical process at hand. Since with this method the impinging light is amplified for a specific wavelength (or range of wavelengths), it will not return any output any if the incoming light does not contain that specified wavelength. I don't want to amplify incoming light, but actually 'create' it. So in theory, when a monochromatic beam with a wavelength in the UV region hits my material slab, it should be absorbed and translated to an emission in the visible spectrum. I don't think this can be achieved with a k_emission as a material property in a second physics node. I thought of maybe having the fluorescent material as a whole still absorb the UV light but then have COMSOL 'create' the emission light with an electric dipole point, do you think this would be possible in COMSOL?


Edgar J. Kaiser Certified Consultant

Please login with a confirmed email address before reporting spam

Posted: 4 years ago Jul 1, 2020, 2:15 p.m. EDT

Hi Karima,

I modified your model to show how you can run study steps and switch variable and physics nodes in the study step settings. I put the n_layer and k_layer variable for the fluorescent layer definition in two global variable nodes. They are a bit easier to reference then in the wee2 and wee3 physics nodes for the fluorescent layer. The ewfd physics node contains two wee-nodes, one for absorption and one for emission.

I study step1 for absorption you see the respective variable node and wee-physics node for absorption active. Analog for study step2 for emission.

You are right, this not yet solving the complete problem. One approach may be to determine absorption as a function of wavelength in the first step and excite the transmit port in the second step by a function that is a wavelength shifted function some way derived from absorption. This may then as well require to switch between two port definitions in the study sequence. There is certainly more than one option to achieve the goal. I guess you will now need some time to play and try things.

Cheers Edgar

-------------------
Edgar J. Kaiser
emPhys Physical Technology
www.emphys.com
Hi Karima, I modified your model to show how you can run study steps and switch variable and physics nodes in the study step settings. I put the n_layer and k_layer variable for the fluorescent layer definition in two global variable nodes. They are a bit easier to reference then in the wee2 and wee3 physics nodes for the fluorescent layer. The ewfd physics node contains two wee-nodes, one for absorption and one for emission. I study step1 for absorption you see the respective variable node and wee-physics node for absorption active. Analog for study step2 for emission. You are right, this not yet solving the complete problem. One approach may be to determine absorption as a function of wavelength in the first step and excite the transmit port in the second step by a function that is a wavelength shifted function some way derived from absorption. This may then as well require to switch between two port definitions in the study sequence. There is certainly more than one option to achieve the goal. I guess you will now need some time to play and try things. Cheers Edgar

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.