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

Automated solvation shell extraction workflow #8

Draft
wants to merge 45 commits into
base: main
Choose a base branch
from

Conversation

jeevster
Copy link
Collaborator

In this PR, we construct a workflow for extracting diverse solvation shells from MD trajectories of electrolyte systems. This workflow operates on PDB files and outputs a directory containing xyz files organized according to solute type and coordination number.

Note: for low-concentration simulations, this workflow is likely sufficient for extracting the solvation shells we require. However, for high concentration simulations, we will need to do a network/topology analysis to extract more interesting snapshots. This will be completed in a future PR.

electrolytes/run_extraction.sh Outdated Show resolved Hide resolved
electrolytes/solvation_shell_extract.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_extract.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_utils.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_utils.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_utils.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_utils.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_utils.py Outdated Show resolved Hide resolved
electrolytes/solvation_shell_extract.py Show resolved Hide resolved
electrolytes/solvation_shell_extract.py Outdated Show resolved Hide resolved
@jeevster
Copy link
Collaborator Author

@levineds I've switched to using the kabsch_rmsd function defined here.

I also included an example pdb file in the testfiles directory.

@jeevster
Copy link
Collaborator Author

On another note, where should we be tracking pip/conda package dependencies for this code?

@levineds
Copy link
Collaborator

levineds commented May 15, 2024

On another note, where should we be tracking pip/conda package dependencies for this code?

I think that goes in om-data/setup.py (i.e. the setup.py in the parent directory)

lattices = np.array([a, b, c])[None][None]

# Choose solute atom
solu = universe.select_atoms(f"name {solute_atom}")
Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Also, I realized that this doesn't yet support multiple atoms in a single solute. I'll get on this.

@jeevster jeevster marked this pull request as draft May 16, 2024 17:25
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.

3 participants