May 2023¶
Monday-Tuesday 1st-2nd May¶
I started uploading my data in progress this morning:
- .phd loader
- Up to date fluorophore dataset
These are teaching days though, so don’t expect too much from me.
DAA is in the RCSPC, but the counts are coming in so slowly that it will probably take 10+ hours to accumulate enough data points. I’m currently waiting on a reply form Toby as to whether I can leave it running overnight,
Wednesday 3rd May¶
The fluorometer is booked out for today, but I’m going to see if I can get some absorbance spectra and lifetimes for the ethanol series.
I didn’t get much of this done, but I did a histogram for the DAA in EtOH decay and I got a quick check of all the EtOH species and their excitations/emissions for easy TCSPC work.
I also re-wrote a chunk of my lifetime fitting work, as I realised that I was doing a lot of the maths wrong. it is now fixed at 12:54am on Thursday…
Thursday 4th May¶
To do:¶
- Continue on EtOH series
- Re-implement TRF saving and loading in the notebooks and datasets
I finished off the EtOH series today!
To code on my way in:¶
- Import the rest of the lifetimes
- Add a separate output widget for each tab in spectraImport.ipynb
Friday 5th May¶
Importing the EtOH series, I’ve noticed that I’m going to need to find a way to extrapolate the excitation corrections to beyond 600nm for the excitation spectra of some of the low energy emitters (NR/R800), otherwise I’m going to need to import these without any corrections.
I did manage to get most of the DMF series, but I think I might need to re-do the TCSPC for a few, because the IRF has erroneous peaks:
- DAA
- NDA
- NR
- R800
Saturday 6th May¶
Well today came as a kick in the teeth, as it turns out the IRF for scatter references needs to be collected at the excitation wavelength, not the emission wavelength.
Turns out that the IRF need to be collected at the excitation wavelength, not at the emission wavelength, so all my species need to be re-done.
Sunday 7th May¶
BIIIIIGGGG day spent coding (not tracked because it was in Jupyter )
I compared my lifetime fitting to that coming out of a commercial software package and I noticed that in log space, the data was able to be fitted properly. This makes sense, as the residuals will more heavily be impacted by the larger changes in the linear space, so the tail end won’t get fitted.
In implementing that change, I also noticed that three and four exponents could be fitted when performing it in log space. The fits I’m getting are now gorgeous!
My ฯ values are also making more sense now as well!
This is definitely making me think that some of my cuvettes were dirty though…
Monday 8th May¶
Teaching day.
Tuesday 9th May¶
I have the TCSPC booked out for the afternoon, so I’m going to try and work through some more lifetime data, and try and figure out what the bump in my data is in the tail-end of of the decay.
I’m also going to use the photon power meter from the laser lab to try and approximate some crude excitation lamp corrections beyond 600nm. These won’t be perfect by any means, but they’ll be a good first approximation in order to continue the trend of the data.
They look alright to me!
Wednesday 10th May¶
Slammed out the lifetimes for the DMF series.
Thursday 11th May¶
Today was spent rationalising and cleaning up the TCSPC/confocal table. Was very tiring, and will continue to tomorrow.
Things to bring up with Rosey/Toby tomorrow:¶
- The ND filter approach works fantastically, except that the ND20 is too dense, so the trigger diode has a really low discretisation voltage
- There is a lot of reflection off the coloured filters, can we use that to trigger the diode?
- Only concern is the amount of scatter
To do tomorrow:¶
- White laser head head mount needs to be re-aligned and tightened up
- Add a second post to the white laser head mount to keep it aligned
- Drill a new hole for the white laser?
- Check the scatter on the MCP to see if the white light being filtered later in the box is an issue
Friday 12th May¶
IT IS DONE! (and I am pooped ๐ฎโ๐จ)
Sunday 14th May¶
We seem to be in a bit of a pickle with M3 atm, as our quota just dropped form 5TB to 2TB, and we have a lot of data on there from former group members, so I’ve spent my day trying to move all of their data over to MonARCH, but I’ll need to wait for Philip to chmod the project folder so that we wan delete some old data.
Monday/Tuesday 15th/16th May¶
Teaching days with a big sleep in.
Thursday 18th May¶
Another lab day, I’m going to be working on the new TCSPC rig today (yay) and will be trying to re-do the nhex series.
Looks like Katya wants to write up Laura’s results in a J Phys Chem Letter, so I’ll probably need to:
- Re-analyse the data - possibly my own just to make sure it’s all correct
- Do cam-b3lyp opts and freqs
- DFT SPs need to be re-done at a TZ level
Friday 19th May¶
I realised yesterday as I was in the lab that the dataset is really small on the low polarity end of the spectral series, so I’ve approved it with Katya and we’re adding three new low-polar species:
- 9-fluorenone (ฮต=260, ฯ=18, ฯ=2.7%)
- Fluoranthene (ฮต=10000, ฯ=?, ฯ=30%)
- 9,10-bis(phenylethynyl)anthracene (ฮต=35400, ฯ=3.6, ฯ=100%)
They’ve been ordered and are just waiting for approval. They’re coming form Truganina, so they should arrive very quickly.
Sunday 21st May¶
This week is going to be a bit of a mess and I’m not sure how much time I’ll get in the lab, so I’m heading in today to try and get started on the DMSO series.
DMSO is a pain to work with in winter…
Monday 22nd May¶
Today is a teaching day, but I have also been pushing through the jobs for the new fluorophores.
Tuesday 23rd May¶
Wednesday 24th May¶
MDP session all day today.
And now I have covid again
Thursday 25th May¶
I spent my Thursday/Friday pulling a hyperfocused all nighter getting a standalone lifetime fitting tool written up. I’ve also published it to pypi as my first module.
https://github.com/adreasnow/LifetimeFittingTool
From Today’s Meeting with Katya¶
- Confirmation in late July
- Write up Laura’s results for confirmation
- “here is the fudge factor for this particular method”
- “n-hexane gives a consistent/systematic shift, since it’s not driven by electrostatics”
- Doodlepoll for confirmation dates
- TCSPC benchmarking against reference species for a write-up
- Use this to determine an uncertainty
Chapters¶
- Intro/review
- Experimental dataset dev
- Theoretical breakdown of 2
- Is going to be based on the outcome of 3
- We’ve found the methods we need, we know the models we can modify
- We don’t know how revolutionary the work will be, but we know where it can go
Tuesday 30th May¶
@ 2800V discretiser max ~ 350mV - optimal ~ 125mV
@ 3000V discretiser max >800mV - optimal ~ 400-600 mV
Important Questions:¶
- Can we use the difference between TDDFT Eq and nEq energies to figure out the stokes correction to the 0-0 energy in abs and fluor
- Can we use the shape of the Theoretical PES to determine what the width of the deconvoluting Gaussians should be for fitting?
Experimental work¶
Problematic spectra:¶
- Collecting:
- AZ lifetimes - all solvents. No solvatochromic shift on s2 → s0
- DMSO with azulene seems to have issues, I’m guessing DMSO excites
- Difficulty fitting
- TCSPC
- DAA/nhex
- BSC/etoh
- TCSPC
Experimental to-do:¶
- Re-do TCSPC
- Dry solvents (waiting on molecular sieves)
- EtOH series
- DMF series (sure/seal)
- DMSO series (sure/seal)
- Anisole series (sure/seal)
- THF series (sure/seal)
UV/FL¶
nhex | tol | ans | ether | chcl3 | thf | dcm | c80h | etoh | acn | dmf | dmso | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
nr | ||||||||||||
nda | ||||||||||||
r800 | insol | |||||||||||
c153 | ||||||||||||
bod493 | ||||||||||||
az | ||||||||||||
aaq | ||||||||||||
dapi | insol | |||||||||||
daa | ||||||||||||
bsc | ||||||||||||
fno | ||||||||||||
bpa | ||||||||||||
fat |
TCSPC¶
nhex | tol | ans | ether | chcl3 | thf | dcm | c80h | etoh | acn | dmf | dmso | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
nr | ||||||||||||
nda | ||||||||||||
r800 | insol | |||||||||||
c153 | ||||||||||||
bod493 | ||||||||||||
az | - | - | - | - | - | - | - | - | - | - | - | - |
aaq | ||||||||||||
dapi | insol | |||||||||||
daa | ||||||||||||
bsc | ||||||||||||
fno | ||||||||||||
bpa | ||||||||||||
fat |
Job Matrix¶
CAS (M3)¶
State | Active Spaces | CASSCF/aug-cc-pvdz | QD-NEVPT2/aug-cc-pvdz |
---|---|---|---|
s0 | |||
s1+2 |
DFT - Gas (MonARCH)¶
State | Functional | Opt | Freq | Pol | Ex/Em | ESD |
---|---|---|---|---|---|---|
s0 | ฯB97X-D | |||||
s1+2 | ฯB97X-D | |||||
s0 | CAM-B3LYP | N/A | ||||
s1+2 | CAM-B3LYP | N/A |
DFT - solv (MonARCH)¶
State | Functional | Opt | Freq | Pol | SS | VEM | ESD |
---|---|---|---|---|---|---|---|
s0 | ฯB97X-D | ||||||
s1+2 | ฯB97X-D | ||||||
s0 | CAM-B3LYP | ||||||
s1+2 | CAM-B3LYP |
To do:¶
- Jobs to figure out
- Test QChem to figure out how to do VEM/cLR/LR/IBSF emission
- Code
- Write energy extractors
- pullQChem_Freq
- \(\Delta E_{ZPVE}\)
- \(\Delta E_{elec}\)
- Number of -ve freqs
- Extend and test for CASSCF
- Q-chem ex/em extractors
- ss-ex
- ss-em
- lr-ex
- lr-em
- pullQChem_Freq
- Write energy extractors
GANTT Chart¶
gantt
dateFormat D/M/YY
title Adrea's PhD Plan
%%excludes weekends
section Preparation
Preliminary reading :active, 5/4/22, 10/2/23
Method "Play" :done, 5/4/22, 10/6/22
Software comparisons :done, 10/6/22, 13/7/22
Preliminary benchmarking :done, 12/7/22, 27/10/22
Christmas break :crit, 22/12/22, 3/1/23
Dataset design :active, 12/7/22, 18/12/22
Experimental design :active, 28/11/22, 27/1/23
section IPR/Confirmation
IPR writeup :ipr, 19/12/22, 10/2/23
IPR :milestone, ipr, 23/2/23, 0d
Confirmation writeup : 29/3/23, 10/5/23
Confirmation :crit, milestone, 24/5/23, 0d
section Dataset paper
Data collection :dc, 6/3/23, 12w
Data processing : 6/3/23, 18w
Writeup :w, after dc, 8w
R&R break :crit, after w, 2w
section Computational dataset analysis paper
Computation : 6/3/23, 30w
Analysis : 1/6/23, 20w
Writeup : 1/8/23, 20w
Christmas break :crit, 22/12/23, 3w
section Solvent model extension/development
Computation : 9/1/24, 5/7/24
Writeup : 5/7/24, 6/9/24
R&R break :crit, 6/9/24, 2w
section Progress review
Progress review writeup : 14/3/24, 9/5/24
Progress review :milestone, 23/5/24, 0d
section Fourth paper (possibly benchmarking or deconvolution)
Computation : 27/9/24, 19/12/24
Christmas break :crit, 20/12/24, 2w
Writeup : 6/1/25, 3/3/25
section End of PhD
Final review writeup : 28/3/25, 9/5/25
Final review :milestone, 23/5/25, 0d
Drop everyhting for thesis : 25/4/25, 25/7/25
Thesis submission :crit, milestone, 22/8/25, 0d
section Conferences
APATCC-11 : 01/09/24, 30/11/24
WATOC Kobe : 22/6/25, 6d
Code Object Structure¶
flowchart LR
subgraph Legend
1((Entity))
3[[module]]
4[(database)]
5[/Data Types/]
6[Hidden Code\n12000+ lines]:::backend
7[User Facing Code\n3000+ lines]:::ui
end
classDef backend fill:#bfbdbd;
classDef ui fill:#fa4ba0;
classDef db fill:#00f2ff;
flowchart TD
a -- Output Files --> g
g -- SLURM Scripts -->
a((HPC Clusters))
subgraph Backend Code
s[/Spectrum Object/]:::backend
m[/Fluorophores/]:::backend
n[/Solvents/]:::backend
o[/PCM Models/]:::backend
p[/Software/]:::backend
r[/Methods/]:::backend
q[/Basis Sets/]:::backend
s & m & n & o & p & r & q -.-> l[[Data Types]]:::backend -.- w
g[[HPC Manager]]:::backend -.- w
w[[Resources module]]:::backend
end
c <-- Job Status --> v
c -- QM Properties --> k
subgraph Databases
j[DB Manager]:::ui -.- b & e & k
e[(ฮE/ฮป/ฮฆ/ฯ DB)]:::backend
b[(Spectra DB)]:::backend
k[(Computational Results DB)]:::backend
v[(Computational Progress DB)]:::backend
v <--> j
end
subgraph Experiment
h[Spectra Importer]:::ui
i[Spectra Analyser/Plotter]:::ui
t -- .txt files --> h
d -- .csv files --> h
t((TCSPC))
d((UV-Vis/Fl))
end
c[Job Manager]:::ui
b <-- Fitted spectra --> i
h -.- i -- Spectral Features --> e
b <-- Fitted spectra --> h
w -. Jobs -.- c
w -.- j
k & e --> u[Stats Scripts]:::ui
classDef backend fill:#bfbdbd;
classDef ui fill:#fa4ba0;
classDef db fill:#00f2ff;