Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Create Boltzmann velocity distributions from Simulations directly #332

Open
wants to merge 12 commits into
base: main
Choose a base branch
from

Conversation

Alexsp32
Copy link
Member

@Alexsp32 Alexsp32 commented Feb 7, 2024

Fixes NQCD/NQCDistributions.jl#9 and requires NQCD/NQCDistributions.jl#12 and NQCD/NQCModels.jl#35 to be merged prior.

Since the VelocityBoltzmann is dependent on the atomic masses, it makes sense to be able to pull both these and the overall dimensions from an existing Simulation. For consistency, a temperature must still be provided separately, but this is quick to do with Simulation.temperature(0).

Also, initialising VelocityBoltzmann from a simulation should check for any atoms not included in mobileatoms, as these should be initialised with zero velocity.

@Alexsp32 Alexsp32 requested review from Nhertl and reinimaurer1 April 9, 2024 10:43
Alexsp32 and others added 5 commits April 14, 2024 10:40
* Add OutputInitial DynamicsOutput
* Stopgap fix for failing unit tests based on memory allocation

* Add allocation test debug messages

* Fixed version logic to be defined for all Julia versions
- Added new DynamicsOutput: OutputDesorptionTrajectory #317
- Added new DynamicsOutput: OutputDesorptionAngle #317
- Created an `Analysis` submodule for common analysis functions

-  Add `Structure` submodule for common utility functions for atomic structure. 

* Add rudimentary tests for NQCDynamics.Structure

* Improved method definitions, possible efficiency gain by using views

* Add unit tests for diatomic analysis functions for desorption

* Add "Analysis" group to unit testing

* Added explanation for Structure module scope
@Alexsp32
Copy link
Member Author

Don't merge until dependent PRs are merged.
Unit tests checked working on my laptop since CI pulls versions without required functionality.

Copy link
Member

@reinimaurer1 reinimaurer1 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

this looks good. I will approve, but can you check that it also works in principle for RingPolymerSimulation objects.

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.

Read parameters for VelocityBoltzmann from NQCDynamics.jl simulation.
2 participants