<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Dan.<br>
<br>
Thanks for your hints and info. It is clear now. <br>
<br>
Few hours ago i already guessed how to do it, and i did as you just
explained. In my case i just updated Amber.frc. I also guessed that
there should be 4 quaterions variables, but thank you for
explanation and pointer to the book.<br>
<br>
Today, i already did some animations via blender. First testing
videos are here<br>
<br>
<a class="moz-txt-link-freetext" href="http://alicepc104.jinr.ru/mvala/openmd-test/">http://alicepc104.jinr.ru/mvala/openmd-test/</a><br>
<br>
I am almost done with automatic script which will import xyz file to
blender (there is very nice blender plugin) and render it to movie.
I would like to make script configurable to add some text, camera
movements and so.<br>
<br>
My another questions are :<br>
<br>
1. How can i add other atoms to interact with this Ico rigid body.
And what do i need to add to Martin.frc file to setup interaction.
Is it enough to setup forces? Like it is explained in openmd manual.<br>
2. How can i make my rigid body fixed in positions (no movements)
and when it will interact with other (new atoms), it will not move.<br>
<br>
Thanks again<br>
<br>
Ciao<br>
<br>
Martin<br>
<br>
<br>
<div class="moz-cite-prefix">On 27.11.2014 20:05, Dan Gezelter
wrote:<br>
</div>
<blockquote cite="mid:30A55DFD-27A0-460C-8FBB-C9C1D258B485@nd.edu"
type="cite">
<pre wrap="">Martin,
You do need to define the atom type in a force field file, and the force field can be quite a bit simpler than Amber for testing. The other issue I saw was that the quaternions in this line:
0 pvqj 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 -1.8e-05 1.6e-06 1.9e-05
only had 3 entries. There should be 3 positions, 3 velocities, 4 quaternions, and 3 angular momentum variables. The convention we use for quaternions is defined on page 88 of Allen & Tildesley's "Computer Simulation of Liquids" book. To orient the principal body axis along the z-spatial coordinate, we could choose the following Euler angles:
\theta = 0
\phi = 0
\psi = 0
These correspond to a set of quaternions q = (1, 0, 0, 0).
Anyway, I tested the three attached files, and you should be able to experiment further with these. A couple of changes I made: 1) The _ (underscore) character has special meaning in the selection syntax, so I modified your atom types and rigid body definition. 2) The angular momentum was quite small, so I boosted it a bit so that you can observe the rotation in this simulation.
Try this with the attached files:
openmd Ico.md
Dump2XYZ -i Ico.dump
jmol Ico.xyz then go to "Tools" -> "Animate" -> "Once"
Best regards,
--Dan Gezelter
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">
</pre>
<blockquote type="cite">
<pre wrap="">On Nov 26, 2014, at 11:01 PM, Martin Vala <a class="moz-txt-link-rfc2396E" href="mailto:mvala@saske.sk"><mvala@saske.sk></a> wrote:
Hi all,
i just notice that i didn't define rigid body in my previous configuration. I am attaching my updated config files. I am sorry for wrong configuration files.
Ciao
Martin
On 27.11.2014 11:37, Martin Vala wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi All,
sorry i misstyped my rigid body. It is icosahedron
Ciao
Martin
On 27.11.2014 11:12, Martin Vala wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Dear OpenMD users,
I would like to define my own rigid body with my own atoms (my custom atom is "My"). My rigid body is defined as icohexadral (in RB_My.md file - attached). I am also attaching RigidBody_01.md file.
Problem:
OpenMD FATAL ERROR:
Can not find Matching Atom Type for[My_RB]
I am guessing that "My_RB" atom type is not defined in $FORCE_PARAM_PATH/Amber.frc (in my RB_My.md, i have forceField = "Amber";). So i will have to define my own MyForceField.frc file (or i will update Amber.frc file),where i will define "My" atom and whole force field defined in my this rigid body (probably define atom atom "My" and its interactions with others, e.g. interaction between atom "My" and "My", "My" and others, and so). Is it right?
If yes, then as first task, i would need to configure this rigid body to be fixed in its initial positions. Can you help me to set it up?
Thank you very much
Ciao
Martin
_______________________________________________
Openmd-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Openmd-users@openmd.org">Openmd-users@openmd.org</a>
<a class="moz-txt-link-freetext" href="http://openmd.org/mailman/listinfo.cgi/openmd-users">http://openmd.org/mailman/listinfo.cgi/openmd-users</a>
</pre>
</blockquote>
<pre wrap="">
_______________________________________________
Openmd-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Openmd-users@openmd.org">Openmd-users@openmd.org</a>
<a class="moz-txt-link-freetext" href="http://openmd.org/mailman/listinfo.cgi/openmd-users">http://openmd.org/mailman/listinfo.cgi/openmd-users</a>
</pre>
</blockquote>
<pre wrap="">
<RB_My.md><RigidBody_01.md>_______________________________________________
Openmd-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Openmd-users@openmd.org">Openmd-users@openmd.org</a>
<a class="moz-txt-link-freetext" href="http://openmd.org/mailman/listinfo.cgi/openmd-users">http://openmd.org/mailman/listinfo.cgi/openmd-users</a>
</pre>
</blockquote>
<pre wrap="">
</pre>
</blockquote>
<br>
</body>
</html>