tutorials for polymer/protein/rna structure simulations

485 views
Skip to first unread message

Olivia Viessmann

unread,
Mar 2, 2022, 10:56:13 AM3/2/22
to hoomd-users
Hi, 
I just started with HOOMD and am interested to simulate polymer/protein/RNA like structures. Is anyone aware of tutorials out there that use HOOMD to simulate these objects? Anything vaguely towards that direction would be greatly appreciated. 

Many thanks
Olivia

Brandon Butler

unread,
Mar 2, 2022, 10:58:33 AM3/2/22
to hoomd...@googlegroups.com
--
You received this message because you are subscribed to the Google Groups "hoomd-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hoomd-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hoomd-users/799a18f6-e05e-45c1-830f-5895ed0a24a8n%40googlegroups.com.
--
Brandon Butler
MolSSI Fellow
PhD Candidate, Chemical Engineering and Scientific Computing | Glotzer Lab, University of Michigan
Email: butl...@umich.edu

Olivia Viessmann

unread,
Mar 2, 2022, 11:25:22 AM3/2/22
to hoomd...@googlegroups.com
Hi,

Thank you for the reply. I am particularly interested in the init.create_random_polymers()
Function that seemed to be available in the hoomd v2 versions. Is there an equivalent in the newer v3 version of hoomd?
I am struggling to figure out how to create a polymer chain from the MD simulation tutorials.

Thanks
Olivia

You received this message because you are subscribed to a topic in the Google Groups "hoomd-users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/hoomd-users/-t0qe4V5PJ0/unsubscribe.
To unsubscribe from this group and all its topics, send an email to hoomd-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hoomd-users/bd2572b3-c739-3b89-2c29-72f496a94490%40umich.edu.

Brandon Butler

unread,
Mar 2, 2022, 11:28:56 AM3/2/22
to hoomd...@googlegroups.com

No there is no such function. Feel free to use Python packages like mbuild to construct polymers. For use see their documentation https://mbuild.mosdef.org/en/stable/index.html.

Olivia Viessmann

unread,
Mar 2, 2022, 11:42:17 AM3/2/22
to hoomd...@googlegroups.com
Hi Brandon,

Thanks for the links. Are you aware of any minimally working example code of a polymer using hoomd v3 that I could start with?

Thanks a ton
Olivia

Jakin Delony

unread,
Mar 2, 2022, 11:45:20 AM3/2/22
to hoomd-users
I had a similar issue in v3 and wrote a small initial configuration generation script for linear homopolymer systems that's attached here.  There are a few caveats, though.  The initial configuration is a naive linear polymer built on a lattice (think pencils in a box) and will require longer relaxation times because the initial state is not very realistic.  Also, in its current form it does not account for complexities such as branching, cross-linking, or copolymerization.  All it does is extend the particle, bond, and angle initialization routine in the tutorial to account for arbitrary chain length (DOP) and number of chains (num_chains).

Hope this helps,
Jakin

gen_lattice_polymer.py

Joshua Anderson

unread,
Mar 2, 2022, 11:50:30 AM3/2/22
to hoomd...@googlegroups.com
Olivia,

The minimal working example was the 2nd link that Brandon sent you: https://hoomd-blue.readthedocs.io/en/latest/howto/molecular.html
------
Joshua A. Anderson, Ph.D.
Research Area Specialist, Chemical Engineering, University of Michigan
> To view this discussion on the web visit https://groups.google.com/d/msgid/hoomd-users/4a204e4c-879f-4e63-b4fb-892b7ed41cb9n%40googlegroups.com.
> <gen_lattice_polymer.py>

Olivia Viessmann

unread,
Mar 2, 2022, 11:56:04 AM3/2/22
to hoomd...@googlegroups.com
Thank you everyone for the input. This is very, very appreciated and super helpful. 
Olivia

Reply all
Reply to author
Forward
0 new messages