Thanks Christoph, I got a python script here in group chat and after modifying it, I can easily make mapp.xml for 500 polymer.
I am running NPT simulation of the polymer with the obtained tabulated potentials. However I have a few queries: 1. For obtaining potentials from IBI, NVT simulations were performed and It reached good convergence (<0.003) but my pressure values were negative. Later I added pressure correction to the settings file. With NVT, values did not change much and although pressure becomes positive, It didn't reach one (1 bar is required pressure). But with NPT and pressure correction, pressure values reached to 1, however the convergence decreased. Should one perform NPT or NVT for potentials with pressure correction. 2. After obtaining potentials, when I run CG NPT simulation. The command is gmx mdrun -s topol.tpr *-tableb table_a0.xvg table_a1.xvg table_a2.xvg table_b0.xvg table_b1.xvg table_b2.xvg* -c confout.gro While we provide bonded tables using -tableb, Is there an option to give non bonded potentials in the command itself. Does this command line will take all non-bonded potentials by itself since we have put energygroup in mdp file. 3. If I want to use GPU for NPT simulations, It gives error that you can't use GPU with user defined tables. You have to give the PME or cut-off. Since we have to use tabulated potentials, does that mean we can't use GPU for such simulations? Thanks in advance On Wed, Feb 24, 2021 at 10:27 PM Christoph Junghans <[email protected]> wrote: > Pragati, > > you could use http://stock.cmm.ki.si/smt/, but for that size of a > system you might need to write a little python script yourself. > > Christoph > > On Wed, Feb 24, 2021 at 3:36 AM Pragati Sharma <[email protected]> > wrote: > > > > Dear all, > > > > I am trying to perform IBI for a 500-mer polymer chain (Each monomer has > 4 carbons and 6 Hs) > > > > I have tabulated potentials from IBI for a 10-mer and 32-mer chain > containing 320 atoms and I will use these potentials to simulate the > 500-mer polymer. > > > > However, a 500-mer polymer contains 5002 atoms. I need conf.gro and > topol.top to run the simulation. And for that I need to map the atomistic > structure files. Making a mapp.xml file for 5002 atoms (Each bead is made > of 3 carbon atoms and 4-5 Hs) and giving individual bonds and angles would > be a tedious task. > > > > > > Is there a way, other than using a mapp file to produce conf.gro for > coarse grained simulations. I have atomistic 500-mer chain simulations and > 10mer-32mer CG simulation files. > > > > Thanks for any idea on this. > > > > -- > > Join us on Slack: https://join.slack.com/t/votca/signup > > --- > > You received this message because you are subscribed to the Google > Groups "votca" group. > > To unsubscribe from this group and stop receiving emails from it, send > an email to [email protected]. > > To view this discussion on the web visit > https://groups.google.com/d/msgid/votca/CAMpF9Q1MAJ1AfUrCCB%2BOv1UAqXNH9z1yAsC%3DC8QuZaEUmqdYig%40mail.gmail.com > . > > > > -- > Christoph Junghans > Web: http://www.compphys.de > > -- > Join us on Slack: https://join.slack.com/t/votca/signup > --- > You received this message because you are subscribed to the Google Groups > "votca" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/votca/CAHG27e5fMSTxy7bWf-My4C1-htb_SyjXBx_N63PXXPwV0L-2fA%40mail.gmail.com > . > -- Join us on Slack: https://join.slack.com/t/votca/signup --- You received this message because you are subscribed to the Google Groups "votca" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/votca/CAMpF9Q16J4EBwuN-WTTqNRC%2B_y%2BCTK-cwiu-nuHrjVM%3DW%3DXMKA%40mail.gmail.com.
