Tips fo pewald3D calculation based on crystal structure

38 views
Skip to first unread message

Paul Hume

unread,
Aug 19, 2020, 8:21:43 PM8/19/20
to votca-ctp
Hi 

I'm a new user of VOTCA, and I'm trying to use the pewald3d calculator for applications in OPV (I have a crystal structure that I'm trying to feed into VOTCA).
Almost everything seems to be working, except I'm not sure how to prepare my system morphology (a triclinic crystal) with a cubic box, as VOTCA requires.
I'm wondering if anyone would be able to share a way of doing this? 
My background is more in quantum chemistry, rather than MD, so I have my fingers crossed that there's an easy way of doing this :)

Thanks very much,

Paul

Denis Andrienko

unread,
Aug 20, 2020, 4:40:42 AM8/20/20
to votc...@googlegroups.com, Paul Hume, Anastasia Markina

Hi Paul,

In principle VOTCA understands triclinic boxes (except the Kinetic Monte Carlo part). Let us skype (crionez) and see if I can help you with specific questions. Anastasia Markina (in cc) did a lot of OPV simulations recently, but she is on a holiday at the moment - working remotely.

Best,
Denis

--
You received this message because you are subscribed to the Google Groups "votca-ctp" group.
To unsubscribe from this group and stop receiving emails from it, send an email to votca-ctp+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/votca-ctp/65a5479e-ce50-4131-ae04-3ccd309bb8f6o%40googlegroups.com.

paulalex...@gmail.com

unread,
Aug 24, 2020, 11:16:43 PM8/24/20
to votca-ctp
Hi Everyone

We resolved this via email, but I'll post the solution in case it's helpful to anyone else.
Most of the calculators in VOTCA are compatible with triclinic boxes (aside from the KMC part).

My problem was that I was using a pdb file that contained unit cell information for my system geometry.
When I ran the command ctp_map -v -t topol.tpr -c test.pdb -s system.xml -f state.sql, I got the following error:

an error occurred:
Non cubical box in pdb file not implemented, yet!

The solution was just to convert my .pdb to a .gro file instead. So I just used the editconf command in Gromacs (specifying my particular periodic boundary conditions).
The calculator then runs fine :)

Thanks again for your help Denis and Anastasia,

Paul
Reply all
Reply to author
Forward
0 new messages