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

Question regarding converting LAMMP data file and dump file to config.Z1 #1

Closed
peter-zhang-chem opened this issue Apr 3, 2024 · 14 comments

Comments

@peter-zhang-chem
Copy link

Dear MK,

I have a LAMMPS data file and a LAMMPs dump file that I'd like to convert to config.Z1 file for analysis. However, inputting the dump file directly in Z1+ program resulted in crashing. I tried to use the script available in the repo to extract the backbone and generate the config.Z1 file; however, I got no output.

Here I attached the test files here. Can you please give me a couple of pointers of how to generate the correct input for the Z1+ program? Thanks!

@peter-zhang-chem
Copy link
Author

Following up on my question, I was able to use Z1+import-lammps.pl to generate a config.Z1 file.

@mkmat
Copy link
Owner

mkmat commented Apr 4, 2024 via email

@peter-zhang-chem
Copy link
Author

Hey!

Thank you for the timely response. I will have to look into the extract-backbone.pl script for to see what might be the issue. I input the same .data and dump file for the extract-backbone.pl and Z1+import-lammps.pl - and extract-backbone.pl was not able to run correctly to output the config.Z1 file. Right now I am trying to analyzing/understand the output I got from the analysis - so I will get back to you at a later time regarding what might have happened with the extract-backbone.pl.

Thanks!

Peter

@mkmat
Copy link
Owner

mkmat commented Apr 4, 2024 via email

@peter-zhang-chem
Copy link
Author

Dear Martin,

Ok, that'd be great!

Here is beginning portion of the data file:

LAMMPS data file. CGCMM style. atom_style full generated by VMD/TopoTools v1.7 on Thu Apr 04 11:43:30 EDT 2024
 9024 atoms
 8960 bonds
 0 angles
 0 dihedrals
 0 impropers
 3 atom types
 1 bond types
 0 angle types
 0 dihedral types
 0 improper types
 -381.353540 1118.646460  xlo xhi
 -325.649424 1174.350576  ylo yhi
 -344.678904 1155.321096  zlo zhi

# Pair Coeffs
#
# 1  A
# 2  C
# 3  G

# Bond Coeffs
#
# 1  

 Masses

 1 39.948002 # A
 2 12.010700 # C
 3 69.723000 # G

 Atoms # full

1 1 2 0.000000 103.589508 164.201080 471.432495 # C CYT
2 2 1 0.000000 102.579193 158.147980 471.438110 # A ADE
3 3 3 0.000000 102.035980 153.521759 474.633545 # G GUA
4 4 2 0.000000 102.872406 151.734192 480.493103 # C CYT
5 5 1 0.000000 102.498383 147.463348 484.110962 # A ADE
6 6 3 0.000000 101.344696 143.881531 488.439941 # G GUA
7 7 2 0.000000 102.650482 139.132507 490.965698 # C CYT
8 8 1 0.000000 106.925201 139.131836 495.002686 # A ADE
9 9 3 0.000000 111.695648 141.219269 498.194336 # G GUA
10 10 2 0.000000 116.481354 140.336823 501.122437 # C CYT
11 11 1 0.000000 121.827545 137.300018 500.431091 # A ADE
12 12 3 0.000000 126.666168 134.139893 502.225952 # G GUA
13 13 2 0.000000 130.583771 130.122253 503.335876 # C CYT
14 14 1 0.000000 129.694427 124.575012 502.921265 # A ADE
15 15 3 0.000000 131.787445 119.330780 502.795776 # G GUA
16 16 2 0.000000 134.016510 114.287109 500.142273 # C CYT
17 17 1 0.000000 139.271881 111.263916 499.137024 # A ADE
18 18 3 0.000000 144.555756 109.422485 500.863159 # G GUA
19 19 2 0.000000 146.606049 106.436554 505.548035 # C CYT
20 20 1 0.000000 144.243561 105.940399 510.782593 # A ADE

and here is the beginning portion of the dump file:

ITEM: TIMESTEP
0
ITEM: NUMBER OF ATOMS
9024
ITEM: BOX BOUNDS pp pp pp
-381.354 1118.65
-325.334 1174.67
-344.769 1155.23
ITEM: ATOMS id type xu yu zu
 1 1 103.59 164.201 471.432

@mkmat
Copy link
Owner

mkmat commented Apr 5, 2024 via email

@peter-zhang-chem
Copy link
Author

Dear Martin,

Hey sorry for the belated response.

  1. Regarding the data file:
    The atoms in data file has the following format:
1 molecule-tag atom-type q x y z nx ny nz  (nx,ny,nz are optional -
  ...                                    see "true flag" input command)
  ...                
  N molecule-tag atom-type q x y z nx ny nz  (N = # of atoms)

Yea, the space before atoms, masses, and bonds is weird - I generated these files from VMD.

  1. Yes, I have a Bonds section in the .data file.

  2. Thank you so much for your time - I have attached the .data file and the dump file with first 2 frames here.
    dump-data-files.zip

@mkmat
Copy link
Owner

mkmat commented Apr 18, 2024 via email

@peter-zhang-chem
Copy link
Author

Dear Martin,

Thank you. I do not think Z1+ returns a lammps data file. Here are the list of files I got from Z1+ output.

  • config.Z1
  • Lpp_values.dat
  • N_values.dat
  • Ree_values.dat
  • Z1+initconfig.dat
  • Z1+SP.dat
  • Z1+summary.dat
  • Z1+summary.html
  • Z_values.dat

I have all the programs you mentioned - currently I am using python to plot values of my interest from the Z1+summary.dat file. It would be awesome to visualize the primitive network used by Z1+ to calculate these values.

And thank you so much, I will wait to see if Clarisse Luap gets back to me.

Peter

@mkmat
Copy link
Owner

mkmat commented Apr 18, 2024 via email

@peter-zhang-chem
Copy link
Author

Thank you so much for your time and this really useful program, Martin! I will check back with you in a few days!

@peter-zhang-chem
Copy link
Author

Dear Martin,

I hope you are well. Just wanted to check in with you regarding the script for creating lammps data file from the Z1+SP.dat file.

Thank you!

Peter

@mkmat
Copy link
Owner

mkmat commented May 3, 2024 via email

@peter-zhang-chem
Copy link
Author

Hey Martin,

Thank you! I was able to visualize the shortest path configuration in vmd using the Z1+SP-to-data.pl.

Peter

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

No branches or pull requests

2 participants