Complete rework of how SRH is calculated

[jmmshn]

The WSWQ is only associated with the initial state so that should just be included as an optional attr.
Then you can call a function on the initial and defect HarmonicDefect objects to compute the SRH coefficient.

Now you can parse an additional folder that contains WSWQ.i files which will allow electron-phonon matrix elements to be calculated for the initial state.

[codecov]

Codecov Report

Base: 90.57% // Head: 88.97% // Decreases project coverage by -1.59% ⚠️

Coverage data is based on head (115a8a6) compared to base (3f0782a).
Patch coverage: 56.00% of modified lines in pull request are covered.

❗ Current head 115a8a6 differs from pull request most recent head 32dd868. Consider uploading reports for the commit 32dd868 to get more accurate results

Additional details and impacted files

@@            Coverage Diff             @@
##             main      #73      +/-   ##
==========================================
- Coverage   90.57%   88.97%   -1.60%     
==========================================
  Files          10       10              
  Lines        1454     1470      +16     
  Branches      233      234       +1     
==========================================
- Hits         1317     1308       -9     
- Misses         78      104      +26     
+ Partials       59       58       -1     

Impacted Files	Coverage Δ
pymatgen/analysis/defects/ccd.py	73.01% <56.00%> (-10.71%)	⬇️
pymatgen/analysis/defects/utils.py	89.95% <0.00%> (-0.84%)	⬇️
pymatgen/analysis/defects/thermo.py	92.36% <0.00%> (-0.04%)	⬇️

Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here.

☔ View full report at Codecov.
📢 Do you have feedback about the report comment? Let us know in this issue.