You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Because I am hoping to get hired in the computer science field when I'm done here, I think its best if my progress is written in the language of the people. Thus, I have a private repo of my work on Github. Some parts of the project need to be hidden from the repo because they are Dr. Nascimento's work (*Bohr*) and I want to respect his work without accidentally publishing something private. That being said, when the project is done, I will release the repo as public to be seen by the masses.
3
+
Because I am hoping to get hired in the computer science field when I'm done here, I think its best if my progress is written in the language of the people. Thus, I have a private repo of my work on GitHub. Some parts of the project need to be hidden from the repo because they are Dr. Nascimento's work (*Bohr*) and I want to respect his work without accidentally publishing something private. That being said, when the project is done, I will release the repo as public to be seen by the masses.
4
4
5
5
### Workflow
6
-
As I have combined Bohr and MEEP into one project, I have create this chart to keep track of the overall workflow of the code.
6
+
As I have combined Bohr and Meep into one project, I have created this chart to keep track of the overall workflow of the code.
Although it is possible to use TDDFT to calculate the polarizability of a molecule under an external field, and it is possible to find dielectric response from polariability, it does not make much sense to simulate something as small as a single molecule as a material with a dielectric response. Thus, simulating a molecule as a point source will do.
11
+
Although it is possible to use TDDFT to calculate the polarizability of a molecule under an external field, and it is possible to find dielectric response from polarizability, it does not make much sense to simulate something as small as a single molecule as a material with a dielectric response. Thus, simulating a molecule as a point source will do.
12
12
13
13
#### Capturing Absorption Spectrum from Molecule
14
14
The TDDFT code assumes that all energy absorbed by the molecule will be re-emitted by it in the form of an induced dipole, we are not able to validate a portion of the simulation by comparing it to experimental absorption spectra.
@@ -17,9 +17,10 @@ The TDDFT code assumes that all energy absorbed by the molecule will be re-emitt
17
17
Since we need to calculate a derivative of the molecule's wavefunction as it relates to the electric field, we collect three timeframes of electric field. Three is the minimum requirement to perform the Runge-Kutta method of derivation.
18
18
19
19
#### Time Scale Between Meep and Bohr
20
-
You are able to calculate the timestep of Meep by `dT = sim.Courant/resolution` but this is in units of Meep time. One time unit in meep is 3.33E-15 seconds. One time unit in bohr is 2.4188843265857E-17 seconds
20
+
You are able to calculate the timestep of Meep by `dT = sim.Courant/resolution`, but this is in units of Meep time. One time unit in Meep is 3.33E-15 seconds. One time unit in Bohr is 2.4188843265857E-17 seconds
21
21
so `dTbohr = dT * (333.3333333333333/2.4188843265857)`
22
22
23
-
#### Induced Dipole to Electric Field
24
-
Since we found it very difficult to calculate the electric field intensity of an electric dipole at its center, we simply calculate the electric field a few picometers from the center of the molecule.
25
-
> **NOTE:** this is the plan, but has not been implemented yet.
23
+
#### Nuances of Simulating a Molecule as a Point Source
24
+
Though the original incident electric field in Meep will be polarized in one coordinate direction (default: `z`), the objects simulated are three-dimensional, so their field response will illicit electric fields in all three directions. All that being said, the electric field collection at the position of the molecule needs to survey all three direction components and pass each into Bohr individually. ([#1](https://github.com/kombatEldridge/plasmonMolecule/issues/1))
25
+
26
+
According to Physics, the electric field due to a point dipole will irradiate in a spherical fashion, but the strength along that sphere will not be equal. We need to create an overall dipole moment from Bohr's calculations and use the equation the issue cited below to find the field strength at every point on the sphere. ([#2](https://github.com/kombatEldridge/plasmonMolecule/issues/2))
0 commit comments