How to calculate center of mass of the polymer from atomic coordinates?

Click For Summary

Homework Help Overview

The discussion revolves around calculating the center of mass and radius of gyration of a polymer using atomic coordinates provided in Armstrongs. The original poster seeks assistance with the formula and methodology for this calculation, particularly with a large dataset of over 2000 atomic coordinates.

Discussion Character

  • Exploratory, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants discuss the formula for calculating the center of mass based on atomic positions and masses. There are requests for clarification on how to apply the formula to specific coordinates and how to compute the radius of gyration from the center of mass result.

Discussion Status

Some participants have provided guidance on the calculations, including how to compute the center of mass and the subsequent steps needed to find the radius of gyration. There is an ongoing exploration of the calculations, with participants questioning the accuracy of results and seeking further clarification on the methodology.

Contextual Notes

The original poster mentions using a large set of atomic coordinates and expresses a need for assistance in converting their results into a usable format for further calculations. There are indications of potential errors in the calculations being discussed.

Kimko
Messages
5
Reaction score
0

Homework Statement



Please help me calculate the center of mass of the polymer (radius) from atomic coordinates (more than 2000) given in Armstrongs. Can anyone help with the formula? I need to calculate the radius of gyration of the polymer. Thanks.

Homework Equations





The Attempt at a Solution


 
Physics news on Phys.org
The center of mass is equal to the sum the atoms positions times their masses all divided by the total mass.
 
Can you please show me how to calculate it using the following coordinates?
Coordinates
ATOM 60 -4.062 10.538 11.856 H
ATOM 61 -2.514 15.625 6.177 N
ATOM 62 -2.780 16.719 5.186 C
ATOM 63 -3.809 17.716 5.736 C
ATOM 64 -4.444 17.483 6.747 O

Thank you.
 
(-4.062 10.538 11.856)*mass(H)+(-2.514 15.625 6.177)*mass(N)+ etc etc. Then divide that by mass(H)+mass(N)+mass(C)+mass(C)+mass(O). To multiply or divide a vector by a number just multiply or divide each component. To add two vectors add component by component.
 
I did it the way you sad and my result is -0.46791 20.58064 19.18521 =Ravg.
Please help how I put it into my formula: Rg= [1/N sum/Ri-Ravg/^2]^(1/2)
 
Your result is wrong. For example, you've got z=19.18. That's bigger than all of your other z coordinates. Show me how you computed the z coordinate of the center of mass and I'll tell what you did wrong.
 
I used all 2000 atomic coordinates in excel, not only those 5. I need to know how to convert
my result in one number that I can use in my formula to calculate Rg. Thanks.
 
Well then. So you've got R_avg. Then you just have to sum (R_i-R_avg)^2. If R_i=(xi,yi,zi) and R_avg=(xa,ya,za), (R_i-R_avg)^2=(xi-xa)^2+(yi-ya)^2+(zi-za)^2.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 20 ·
Replies
20
Views
3K
  • · Replies 4 ·
Replies
4
Views
6K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
1
Views
2K
Replies
8
Views
6K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
1
Views
2K