On Thu, Feb 25, 2021 at 11:17 PM Pragati Sharma <[email protected]> wrote: > > 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.
I am not 100% sure IBI in NPT is really well researched. Pressure correction is bit magically especially when you have multiple bead types, in principle one would need to match the partial pressure, but as no MD code has partial pressures implemented, VOTCA cannot do it either. The only advice I can give for the pressure correction is to make sure your trajectories are long enough, the error bars/ fluctuations on the pressure are usually huge. > > 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. Yes, I believe you don't need an exxtra option for the non-bonded tables (at least in gromacs-2019). > > 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? That is unfortunately right, if you have any question on this issue, ask on the gromacs user list ([email protected]) Alternatively you could use LAMMPS or hoomd-blue for your md simulation. Christoph > > 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. -- 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/CAHG27e4mN%3DxeV7F-_qNNLLAnk8tGXppRzWggMUsjSS-viXEEXw%40mail.gmail.com.
