Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.
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.
How-to : Transmission Spectra vs Wavelength
Posted 2013年11月13日 GMT-5 20:55 RF & Microwave Engineering, Wave Optics Version 4.2, Version 4.2a, Version 4.3, Version 4.3a, Version 4.3b, Version 4.4 9 Replies
Please login with a confirmed email address before reporting spam
I had attached the concept of my simulation as picture where you can understand my idea.
For obtaining Wavelength vs Transmission spectra I need to include all Refractive index (n and K values of GaAs). I had imported Refractive index values via interpolation. I had converted n and k values to Relative permitivity say GaAs_real and GaAs_imag.
Previously in Parameters I had defined Lambda = 1000 and frequency = c_const/Lambda[nm]. I had used nm in Lambda since imported (ex. GaAs_real) wavelength values as say 800 ..... without nm.
In GaAs Material properties I had defined Relative permittivity as = GaAs_real(Lambda)-j*GaAs_imag(Lambda). I think this is correct definition.
So In Parametric sweep now I can do Lambda dependent.
My Goal is to get PL spectra ( Transmission ) vs wavelength. I had placed Point in GaAs and used that as Dipole source. My target is to achieve Transmission spectra coming out from GaAs vs wavelength (say 800-1200nm). I m really confused here how-to-do. For obtaining transmission spectra do I need to use port condition or what ?
Attachments:
Please login with a confirmed email address before reporting spam
Certainly, you could put a boundary probe on the top air boundary to record the transmitted field, but there are several questions to determine. It is far enough away from the GaAs/Au structure to form a uniform field appropriate for retrieval? How do you determine the baseline "incident" power from a point dipole inside a small structure (i.e. in the near field)?
I've set up a short example that is probably not very accurate, but at least demonstrates the challenges of trying to do what you have requested. I set up a boundary probe on both the air and the GaAs/Air interfaces, as well as a point domain probe placed near the dipole. The results of the normalized and averaged power transmittances are shown in the global 1D plots, and show the discrepancies between using the GaAs/air boundary and the point probe.
Keep in mind I am working in "freq" terms, but either you can switch it around for "lamb", or you can just edit my "Parameteric 1" component to sweep what wavelengths you want.
Attachments: GaAs.mph
Reply | Reply with Quote | Edit
Vasanthan Devaraj
November 15, 2013 5:09am EST in response to Bryan Adomanis
Re: Adding Refractive Index
Hi Bryan Adomanis
First of all thanks for your model and try. Its really helpful since I m just beginner in COMSOL. I had attached the geometry of my structure . Previously i used Scattering BC, but now changed to PML as per your ex.
For obtaining transmission(hereafter T or PL) spectra, i think we can adjust the height of transmission monitor to get the best from GaAs top surface.
Coming to your 2nd point - //How do you determine the baseline "incident" power from a point dipole inside a small structure (i.e. in the near field)?// I couldn't understand this ? Why we have to determine the baseline incident power ?
I understood that you used Boundary probes for obtaining T. In Relative permittivity it should be square of refractive index right ?
And in Dipole, you had chosen Dipole moment . How about choosing Magnitude and direction and is it possible for getting T??
Please login with a confirmed email address before reporting spam
So I mention the incident field/power because that is what is necessary to compute transmittance. When you say "transmission monitor", I am not sure what you are indicating. Certainly the boundary probe records the full field, but it does not compute a transmittance. In order to get transmittance, you must determine a ratio of your recorded field at some region of space over the source field. Normally this is done via a set of ports, where one is both an active "source" boundary AND a "monitor" for the reflected field, while the other is only the monitor for the transmitted field. The source is generated as a plane wave with some known field strength, orientation, power and phase at the boundary, and so the data recorded at the transmission port references these inputs to compute the spectral power transmittance.
I do not know of a way that this is possible for a point source in the near field. Physically, it doesn't make sense to record any field data near the dipole for comparison to the field at the air boundary to get some transmission. First of all, you are measuring the full field, and have no knowledge of the original source field, other than the analytical solution for a field from a point dipole. Secondly, the field is a near field, and thus cannot give you any real comparison to a far-field plane wave needed for a transmittance determination.
Maybe I am still not understanding your problem sufficiently. Could you explain what is it you are trying to acquire, and what recorded data you expect to need to get it?
Please login with a confirmed email address before reporting spam
Coming to this topic, I will explain my focus. Actually I did FDTD simulation with this structure. Why i m preferring dipole source means, I m dealing with Quantum dots (QD). In FDTD, we use dipole source as in meaning of QD. So trying apply the same in COMSOL. I had calculated the transmission spectra vs wavelength. I had attached the copy of result for my structure using FDTD (black line spectra for our current target). It is similar with my experiment results. I think you know my structure now. I think now it will be more clear for you.
So for each wavelength n and k values differ and thats y we are importing n and k values via interpolation. So in COMSOL we have to prefer QD by means of Dipole source right? So with this we have calculate Transmission spectra vs wavelength.
Attachments:
Please login with a confirmed email address before reporting spam
In FDTD, you generally have good information about the source pulse, as well as the information as to how it was transformed to a continuous wave approximation, so recording transmission data requires less consideration in setup. In Lumerical FDTD, for example, you record transmission simply by setting up a transmission "monitor", which takes the complex field and power data and automatically normalizes it to the source field/power to give you the complex transmission coefficient(t) and the spectral transmittance (T). You do not need to explicitly reference the source information for the t or T calculations.
In COMSOL, you must explicitly reference the source data. This is why you normally set up TWO ports--one to define the plane-wave source and record complex reflection data, another to record complex transmission data. Without this source port--i.e. with only the dipole source--I do not know how you plan to identify the input field/power and phase to perform the t and T calculations. Certainly you could extend the air boundary or do a far-field projection to get a plane wave solution for your "output", but how do you handle the near-field "input" to get your normalization factor? Where do you place your monitor?
Please login with a confirmed email address before reporting spam
I clearly got your point this time regarding COMSOL in dealing with dipole sources. Since I m a beginner in COMSOL, I just thought placing monitor kind of study will be similar to get data as like Lumerical FDTD. May be i think its better I will do this wavelength vs T calculations in FDTD. I m really thankful for your responses, explanations and your solutions. I learned in this progress how to use COMSOL.
Actually i think finding electric field with Point source will be easy using COMSOL. I think for this we dont need port conditions and all. I m now interested in this which will be useful for my studies.
I want to simulate the Vertical and Horizontal polarization of dipole source. May be better I can post this as new topic of discussion which will be easier for us to discuss. I would be glad and thankful if you could help me.
Anyways I can give a description here too.
Its the same structure with electric dipole source. For this we need to import n and k values for both GaAs and Gold right ? I did it and imported via interpolation. Since I m dealing in Lambda I made arguments in nm. The data source file for GaAs and Au is attached below. I got that from refractive index.info site. I think for each wavelength, the values are not listed.So this would be fine for simulation?? And I had defined the relative permittivity in GaAs and Au. Can you check if thats right ?? Can we use Scat B.C or PML for this ?
For Vertical and Horizontal Polarization, I chose Dipole : Magnitude and direction. So in this x = 0 and y=1 for V-Pol and vice versa for H-Pol right ? Correct me if I m wrong.
Thank you.
Attachments:
Please login with a confirmed email address before reporting spam
Second, it looks like the material properties are properly set: the indexes (and thus the permittivities) are functions of "lambda" and this is the variable you are parametrizing. You can check if you have done it properly by selecting the interpolation and hitting "plot", and look at the associated curve on RefractiveIndex.Info site. You ask " I think for each wavelength, the values are not listed.So this would be fine for simulation??" This is precisely what the interpolation has done for you--converted discrete data points (from published data) into a continuous curve that can be sampled for ANY wavelength in that range!
Now, back to the dipole. Yes, you can control the polarization in the manner in which you state. You could even place multiple dipoles of various magnitudes and directions on top of each other to simulate coherent, unpolarized emission (I think there is a laborious way to do incoherent emission). However, the question remains as to what kind of data you would like to retrieve. You mention "finding electric field with Point source will be easy using COMSOL"--you are correct, it is! Just place the boundary/point probe where you want to sample (average, integrate, etc.) the field generated (and scattered) from the point source. But so what? Is the full-field output all the information you need? You will not get a transmittance that way! You won't even isolate the scattered field that way. So what is it you WANT to get out of the simulation?
Please login with a confirmed email address before reporting spam
Thanks for your reply.
If you dont mind, I had posted the reply in my relevant topic since this is getting diverted here.
www.comsol.com/community/forums/general/thread/40765/#p110687
Please login with a confirmed email address before reporting spam
I have set up n and k interpolations for the wavelenght range of interest similar to the description above.
When I plot the reflectance I am getting a highly oscillatory result. Does anyone know why this might be happening? When I dont have any structures and keep the surface flat, I get the right reflectance values from literature.
More details on my problem are on this thread www.comsol.com/community/forums/general/thread/80041/#p184421
I have attached images of the results and the structure of interest.
Any help on this would be great. Thank you.
Attachments:
Please login with a confirmed email address before reporting spam
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.