Skip to content

relaxonVelocities branch: correcting lifetimes & viscosity iiii contributions written to hdf5#273

Open
elliotperviz wants to merge 1 commit into
phoebe-team:relaxonVelocitiesfrom
elliotperviz:relaxonVelocities
Open

relaxonVelocities branch: correcting lifetimes & viscosity iiii contributions written to hdf5#273
elliotperviz wants to merge 1 commit into
phoebe-team:relaxonVelocitiesfrom
elliotperviz:relaxonVelocities

Conversation

@elliotperviz

@elliotperviz elliotperviz commented Jun 23, 2026

Copy link
Copy Markdown

Two small changes:

  • i) apply 2pi division to relaxon lifetimes written to hdf5 files during output (src/observable/relaxons.cpp)
  • ii) move iiiiContrib out of nested for loops to prevent duplicated relaxon contributions to each viscosity component (src/observable/transport_coefficients.cpp)

-> means that kappa or viscosity contributions reconstructed using relaxon transport properties (written during output to hdf5) match within numerical tolerance the kappa / viscosity contributions (and totals) written to json.

Note that in i), the division by 2pi is only performed if isPhonon() evaluates to true. Right now I only have phonon BTE inputs/outputs at hand, so I did not check in my calculations if the same issue exists for the electronBTE. Probably the twoPi division should be applied in both cases.

…nd move viscosity iiii contribution out of i,j for loops to be called once per relaxon mode
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant