On Mon, Jun 12, 2023, 03:37 Kankana Bhattacharjee <
[email protected]> wrote:

> And, sir while using the csg_call command for single propane, then what
> settings.xml file has to be used to convert into table.xvg format ? Because
> for single propane folder, there is no settings.xml file.
>
The settings file is just a simple file that contains some more information
about the interaction, this file should work:
https://github.com/votca/votca/blob/master/csg-tutorials/propane/ibi/settings.xml

When you run csg_call it will tell you which items are read.


>
>
> [image: Mailtrack]
> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>  Sender
> notified by
> Mailtrack
> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>  06/12/23,
> 03:03:46 PM
>
> On Mon, Jun 12, 2023 at 3:00 PM Kankana Bhattacharjee <
> [email protected]> wrote:
>
>> Ok. Thanks a lot sir. But, what would be the maximum value for specifying
>> grid in case of bonds ?
>>
>>
>>
>> [image: Mailtrack]
>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>  Sender
>> notified by
>> Mailtrack
>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>  06/12/23,
>> 02:59:32 PM
>>
>> On Mon, Jun 12, 2023 at 8:09 AM Christoph Junghans <[email protected]>
>> wrote:
>>
>>>
>>>
>>> On Fri, Jun 9, 2023 at 5:29 AM Kankana Bhattacharjee <
>>> [email protected]> wrote:
>>>
>>>> Sir, I am not understanding the Resampling and extrapolation steps for
>>>> getting bonded.xvg files. For resampling options, in the command line:
>>>>   csg_resample --in table.pot --out table_resample.pot \ --grid
>>>> min:step:max
>>>>
>>>> How to choose min, step and max values here ? Kindly help me with
>>>> this.
>>>>
>>> There are some conditions that gromacs wants:
>>> - bonds type need to start at 0 to max (in nm)
>>> - angle type need to from 0 to 180 degrees
>>> - dihedrals type need to go from -18o t0 180 degrees
>>> (VOTCA internally uses radians, so you will need to convert that)
>>>
>>> Most of that is implemented in:
>>>
>>> https://github.com/votca/votca/blob/master/csg/share/scripts/inverse/potential_to_gromacs.sh
>>>
>>> You can actually try to use that with:
>>> $ csg_call --ia-type bond --ia-name bond1 --settings.xml table.xml
>>> convert_potential gromacs table.pot table.xvg
>>>
>>> Also see: https://www.votca.org/csg/preparing.html#exporting-the-table
>>>
>>> Christoph
>>>
>>>>
>>>>
>>>>
>>>>
>>>> [image: Mailtrack]
>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>  Sender
>>>> notified by
>>>> Mailtrack
>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>  06/09/23,
>>>> 04:58:12 PM
>>>>
>>>> On Fri, Jun 9, 2023 at 4:05 PM Kankana Bhattacharjee <
>>>> [email protected]> wrote:
>>>>
>>>>> Dear Sir,
>>>>>
>>>>> When converting the Boltzmann inverted potential (angle.pot.ib) into
>>>>> table.xvg using the command line : *csg_call --ia-type angle
>>>>> --ia-name angle1 --options table_a1.xml \ convert_potential gromacs
>>>>> angle_rad.pot.ib table_a1.xvg*
>>>>> then, getting the error like this:
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> *table_linearop.pl <http://table_linearop.pl>: angle_rad.pot.ib to
>>>>> angle_rad.pot.scale.3tqbn with x' = 57.2957795*x + 0readin_table: Not
>>>>> enought columns in line 1 in file angle_rad.pot.ib, if you don't have 
>>>>> flags
>>>>> in your table add --sloppy-tables option to
>>>>> csg_callCallstack:/usr/bin/csg_call - linenumber 199    do_external -
>>>>> linenumber 176 in /usr/share/votca/scripts/inverse/functions_common.sh
>>>>>   /usr/share/votca/scripts/inverse/potential_to_gromacs.sh - linenumber
>>>>> 124            do_external - linenumber 22 (see 'csg_call --cat function
>>>>> do_external')                die - linenumber 2 (see 'csg_call --cat
>>>>> function
>>>>> die')######################################################################################################################
>>>>>
>>>>>                                       ## ERROR:
>>>>>
>>>>>    ## do_external: subscript
>>>>>                                              ##
>>>>> /usr/share/votca/scripts/inverse/table_linearop.pl
>>>>> <http://table_linearop.pl> --on-x angle_rad.pot.ib
>>>>> angle_rad.pot.scale.3tqbn 57.2957795 0 ## (from tags table linearop) 
>>>>> failed
>>>>>
>>>>>     ## Details can be found above
>>>>>                                              ##
>>>>>
>>>>>
>>>>> ######################################################################################################################Terminated*
>>>>>
>>>>> But, for bonded distributions table.xvg files generated. For angle
>>>>> potential, not getting the table and not able to figure out the error.
>>>>> Kindly help me with this.
>>>>>
>>>>> Thanks & Regards
>>>>> Kankana Bhattacharjee
>>>>>
>>>>>
>>>>>
>>>>> [image: Mailtrack]
>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>  Sender
>>>>> notified by
>>>>> Mailtrack
>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>  06/09/23,
>>>>> 04:02:39 PM
>>>>>
>>>>> On Fri, Jun 9, 2023 at 3:10 PM Kankana Bhattacharjee <
>>>>> [email protected]> wrote:
>>>>>
>>>>>> Dear Sir,
>>>>>>
>>>>>> For bonded distribution, if I want to convert the table into xvg
>>>>>> format, then for the "Resampling step" the command line is, I want to 
>>>>>> know
>>>>>> what would be the "min:step:max" options for bonded distributions ? I am
>>>>>> not understanding this part.
>>>>>>
>>>>>> csg_resample --in table.pot --out table_resample.pot \             
>>>>>> --grid min:step:max
>>>>>>
>>>>>> I want to know what would be the *"min:step:max"* options for bonded 
>>>>>> distributions ? I am not understanding this part.
>>>>>>
>>>>>>
>>>>>> Thanks & Regards
>>>>>>
>>>>>> Kankana Bhattacharjee
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> [image: Mailtrack]
>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>  Sender
>>>>>> notified by
>>>>>> Mailtrack
>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>  06/09/23,
>>>>>> 03:08:17 PM
>>>>>>
>>>>>> On Fri, Jun 9, 2023 at 1:19 PM Kankana Bhattacharjee <
>>>>>> [email protected]> wrote:
>>>>>>
>>>>>>> Sir, for getting the xvg table, at first I directly created
>>>>>>> table.xml file for bond length distribution and then used csg_call 
>>>>>>> command
>>>>>>> to convert the table into xvg format. But, I excluded the "Resampling" 
>>>>>>> and
>>>>>>> "Extrapolation" options. Is it the correct path ?
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> [image: Mailtrack]
>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>  Sender
>>>>>>> notified by
>>>>>>> Mailtrack
>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>  06/09/23,
>>>>>>> 01:14:13 PM
>>>>>>>
>>>>>>> On Thu, Jun 8, 2023 at 11:07 PM Christoph Junghans <
>>>>>>> [email protected]> wrote:
>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, Jun 8, 2023 at 07:54 Kankana Bhattacharjee <
>>>>>>>> [email protected]> wrote:
>>>>>>>>
>>>>>>>>> Sir, to convert the bonded potentials into xvg format, at first I
>>>>>>>>> have made table_b1.xml for bond1 type and then, used the command line:
>>>>>>>>>  csg_call --ia-type bond --ia-name bond1 --options table_b1.xml \
>>>>>>>>> convert_potential gromacs angle.pot.ib table_b1.xvg. But, the output 
>>>>>>>>> file:
>>>>>>>>> table_b1.xvg looks somewhat different from the file given in the ibi
>>>>>>>>> folder. Kindly provide some suggestions. Sir, my table_b1.xvg is 
>>>>>>>>> looking
>>>>>>>>> somewhat different from that of table_b1.xvg provided in the ibi 
>>>>>>>>> folder.
>>>>>>>>> For converting into table.xvg, should "Resampling", and 
>>>>>>>>> "Extrapolation"
>>>>>>>>> options also have to be followed in case of bonded distributions ?
>>>>>>>>>
>>>>>>>> It depends on the system, but usually yes!
>>>>>>>>
>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Thu, Jun 8, 2023 at 6:02 PM Christoph Junghans <
>>>>>>>>> [email protected]> wrote:
>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Thu, Jun 8, 2023, 06:29 Kankana Bhattacharjee <
>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>
>>>>>>>>>>> Sorry sir, bonded section is also there. To obtain a bond table,
>>>>>>>>>>> the used command line is:
>>>>>>>>>>>
>>>>>>>>>>> csg_call --ia-type bond --ia-name XXX --options table.xml \
>>>>>>>>>>> convert_potential gromacs table_extrapolate.pot table.xvg
>>>>>>>>>>>
>>>>>>>>>>> I would like to know, in --ia-name flag, what I should I put
>>>>>>>>>>> there, would it foolow the name which is mentioned in the 
>>>>>>>>>>> propane.xml file
>>>>>>>>>>> in bonded section (bond1, bond2) ?
>>>>>>>>>>>
>>>>>>>>>> Yes, it is just the same names you defined in table.xml.
>>>>>>>>>>
>>>>>>>>>> Christoph
>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>  Sender
>>>>>>>>>>> notified by
>>>>>>>>>>> Mailtrack
>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>  06/08/23,
>>>>>>>>>>> 03:57:07 PM
>>>>>>>>>>>
>>>>>>>>>>> On Thu, Jun 8, 2023 at 3:46 PM Kankana Bhattacharjee <
>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> But, sir for this process, a table.xml file is given for
>>>>>>>>>>>> non-bonded intercations. But, how to do the same job for bonded
>>>>>>>>>>>> distributions ?
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>  Sender
>>>>>>>>>>>> notified by
>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>  06/08/23,
>>>>>>>>>>>> 03:44:48 PM
>>>>>>>>>>>>
>>>>>>>>>>>> On Thu, Jun 8, 2023 at 3:42 PM Christoph Junghans <
>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Jun 8, 2023, 06:01 Kankana Bhattacharjee <
>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Ok. Thank you sir. I understood that, table number concept
>>>>>>>>>>>>>> i.e, table_a1.xvg, table_b1.xvg. But, how can one generate those 
>>>>>>>>>>>>>> files to
>>>>>>>>>>>>>> supply in bonded.itp file ? I can see those table (bond1, 
>>>>>>>>>>>>>> angle1).xvg files
>>>>>>>>>>>>>> in the ibi directory. But, how those files get generated ?
>>>>>>>>>>>>>>
>>>>>>>>>>>>> https://www.votca.org/csg/preparing.html#exporting-the-table
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>> ---------- Forwarded message ---------
>>>>>>>>>>>>>> From: Christoph Junghans <[email protected]>
>>>>>>>>>>>>>> Date: Thu, Jun 8, 2023 at 3:27 PM
>>>>>>>>>>>>>> Subject: Re: [votca] Regarding-Getting-itp-files:
>>>>>>>>>>>>>> To: Kankana Bhattacharjee <[email protected]>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> On Thu, Jun 8, 2023, 01:43 Kankana Bhattacharjee <
>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I would like to know, how in the case of CG propane.itp
>>>>>>>>>>>>>>> file, I have found that, "func type 8" is used which is 
>>>>>>>>>>>>>>> tabulated
>>>>>>>>>>>>>>> functional form. And, in the case of "func type 8", there are 
>>>>>>>>>>>>>>> two
>>>>>>>>>>>>>>> parameters "table number" and "spring constant" values. And, 
>>>>>>>>>>>>>>> how two
>>>>>>>>>>>>>>> parameters (table number, spring constant) are obtained (which 
>>>>>>>>>>>>>>> is equal to
>>>>>>>>>>>>>>> 1 in topol.top file) ? Like, as I am doing it manually, and 
>>>>>>>>>>>>>>> bonded values
>>>>>>>>>>>>>>> obtained from BI method, then how can I mention spring 
>>>>>>>>>>>>>>> constant, and table
>>>>>>>>>>>>>>> number values in bonded.itp file ?
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>> The spring constant is yet another number the table gets
>>>>>>>>>>>>>> multipled with, but for BI you don't need another factor that 
>>>>>>>>>>>>>> hence we set
>>>>>>>>>>>>>> it to 1.0!
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> The table number is arbitrary, but needs to consistent with
>>>>>>>>>>>>>> the files name, table_b<#>.xvg and table_a<#>.xvg. (b for bond, 
>>>>>>>>>>>>>> a for
>>>>>>>>>>>>>> angle, see the Gromacs manual for details).
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thanks and regards
>>>>>>>>>>>>>>> Kankana Bhattacharjee
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>  06/08/23,
>>>>>>>>>>>>>>> 11:08:09 AM
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 7:13 PM Kankana Bhattacharjee <
>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Ok sir
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023, 19:00 Christoph Junghans <
>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 7:11 AM Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Ok sir.
>>>>>>>>>>>>>>>>>> If i use CG MD only using gromacs without using votca, is
>>>>>>>>>>>>>>>>>> it possible to do with gromacs 2020 or later version? If i 
>>>>>>>>>>>>>>>>>> don't use
>>>>>>>>>>>>>>>>>> iterative methods
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Well even the non-iterative methods like Boltzmann
>>>>>>>>>>>>>>>>> inversion and Force Matching use tables.
>>>>>>>>>>>>>>>>> You would have to fit them to a Lennard-Jones form to use
>>>>>>>>>>>>>>>>> Gromacs 2020.
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023, 18:16 Christoph Junghans <
>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 6:35 AM Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Ok sir. What should I do to fix it ?
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Build VOTCA against gromacs-2019. VOTCA even has an
>>>>>>>>>>>>>>>>>>> option to build its own gromacs (BUILD_OWN_GROMACS).
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Is it true for all iterative methods ?
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Yes, but of course you can always use lammps to do the
>>>>>>>>>>>>>>>>>>> CG MD.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>  06/05/23,
>>>>>>>>>>>>>>>>>>>> 06:05:14 PM
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 6:04 PM Christoph Junghans <
>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 6:17 AM Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> This is the Fatal error obtained from inverse.log
>>>>>>>>>>>>>>>>>>>>>> file:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> *Fatal error:The group cutoff scheme has been removed
>>>>>>>>>>>>>>>>>>>>>> since GROMACS 2020. Please use theVerlet cutoff scheme.*
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> That is right to use tabulated interactions you will
>>>>>>>>>>>>>>>>>>>>> to use gromacs 2019 (or lammps).
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> VOTCA also prints a warning about that issue then you
>>>>>>>>>>>>>>>>>>>>> build it.
>>>>>>>>>>>>>>>>>>>>> Status of tables in newer gromacs versions here:
>>>>>>>>>>>>>>>>>>>>> https://gitlab.com/gromacs/gromacs/-/issues/1347, but
>>>>>>>>>>>>>>>>>>>>> unfortunately not much has moved.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Thanks & Regards
>>>>>>>>>>>>>>>>>>>>>> Kankana Bhattacharjee
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>  06/05/23,
>>>>>>>>>>>>>>>>>>>>>> 05:46:29 PM
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> On Mon, Jun 5, 2023 at 5:32 PM Christoph Junghans <
>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> On Sun, Jun 4, 2023, 23:15 Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> RDF for A-A, B-B, A-B computed from atomistic
>>>>>>>>>>>>>>>>>>>>>>>> simulation using csg_stat tool. Now, I am doing ibi 
>>>>>>>>>>>>>>>>>>>>>>>> method. But, error is
>>>>>>>>>>>>>>>>>>>>>>>> coming during step 1 like this:
>>>>>>>>>>>>>>>>>>>>>>>> Appending to existing logfile inverse.log
>>>>>>>>>>>>>>>>>>>>>>>> We are doing Method: ibi
>>>>>>>>>>>>>>>>>>>>>>>> step 0 is already done - skipping
>>>>>>>>>>>>>>>>>>>>>>>> Doing iteration 1 (dir step_001)
>>>>>>>>>>>>>>>>>>>>>>>> Simulation with gromacs
>>>>>>>>>>>>>>>>>>>>>>>> Automatically added 'cutoff-scheme = Group' to
>>>>>>>>>>>>>>>>>>>>>>>> grompp.mdp, tabulated interactions only work with 
>>>>>>>>>>>>>>>>>>>>>>>> Group cutoff-scheme!
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ##################################################################################################################
>>>>>>>>>>>>>>>>>>>>>>>> #
>>>>>>>>>>>>>>>>>>>>>>>>                                                        
>>>>>>>>>>>>>>>>>>>>>>>>       #
>>>>>>>>>>>>>>>>>>>>>>>> # ERROR:
>>>>>>>>>>>>>>>>>>>>>>>>                                                        
>>>>>>>>>>>>>>>>>>>>>>>>        #
>>>>>>>>>>>>>>>>>>>>>>>> # critical: 'gmx grompp -n index.ndx -f grompp.mdp
>>>>>>>>>>>>>>>>>>>>>>>> -p topol.top -o topol.tpr -c conf.gro' failed          
>>>>>>>>>>>>>>>>>>>>>>>>        #
>>>>>>>>>>>>>>>>>>>>>>>> # For details see the logfile
>>>>>>>>>>>>>>>>>>>>>>>> /mnt/c/Users/KANKANA/Downloads/votca/csg-tutorials/propane/ibi-kankana/inverse.log
>>>>>>>>>>>>>>>>>>>>>>>> #
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Please look up and post the detailed error message
>>>>>>>>>>>>>>>>>>>>>>> from inverse.log.
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> #
>>>>>>>>>>>>>>>>>>>>>>>>                                                        
>>>>>>>>>>>>>>>>>>>>>>>>       #
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> ##################################################################################################################
>>>>>>>>>>>>>>>>>>>>>>>> Terminated
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Kindly help me with this.
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> Thanks & Regards
>>>>>>>>>>>>>>>>>>>>>>>> kankana Bhattacharjee
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>  06/05/23,
>>>>>>>>>>>>>>>>>>>>>>>> 10:44:51 AM
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023 at 8:42 PM Christoph Junghans <
>>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023 at 9:00 AM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <[email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> Thank you sir. But IBI would be for multiple
>>>>>>>>>>>>>>>>>>>>>>>>>> propane molecules. Isn't it ? It won't be for single 
>>>>>>>>>>>>>>>>>>>>>>>>>> propane i think
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>> Yes, IBI (and IMC and RE) are for multiple
>>>>>>>>>>>>>>>>>>>>>>>>> molecules.
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023, 19:47 Christoph Junghans <
>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023 at 3:16 AM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <[email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> Sir, i performed atomistic simulation of
>>>>>>>>>>>>>>>>>>>>>>>>>>>> single  propane in vaccum. But, if we use multiple 
>>>>>>>>>>>>>>>>>>>>>>>>>>>> propane simulation in
>>>>>>>>>>>>>>>>>>>>>>>>>>>> vaccum, then csg_stat can give RDF which can be 
>>>>>>>>>>>>>>>>>>>>>>>>>>>> compared with RDF of
>>>>>>>>>>>>>>>>>>>>>>>>>>>> ibi/imc methods ?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> The single molecule in vacuum is a common
>>>>>>>>>>>>>>>>>>>>>>>>>>> separation strategy to determine the intertra 
>>>>>>>>>>>>>>>>>>>>>>>>>>> molecular interactions. In
>>>>>>>>>>>>>>>>>>>>>>>>>>> short get the bonded interactions from BI for the 
>>>>>>>>>>>>>>>>>>>>>>>>>>> single chain run and then
>>>>>>>>>>>>>>>>>>>>>>>>>>> run IBI for the non-bonded interactions.
>>>>>>>>>>>>>>>>>>>>>>>>>>> see the Tschoep paper from 1998 (
>>>>>>>>>>>>>>>>>>>>>>>>>>> https://www.votca.org/csg/bibliography.html#tschoep-1998
>>>>>>>>>>>>>>>>>>>>>>>>>>> )
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> I want to know how to choose which method (ibi,
>>>>>>>>>>>>>>>>>>>>>>>>>>>> imc) would be accurate to give accurate structural 
>>>>>>>>>>>>>>>>>>>>>>>>>>>> properties?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>> That depends on the system. IBI, IMC & Re are
>>>>>>>>>>>>>>>>>>>>>>>>>>> all structure-based methods and try to reproduce 
>>>>>>>>>>>>>>>>>>>>>>>>>>> the structure.
>>>>>>>>>>>>>>>>>>>>>>>>>>> Usually IBI is more robust, but needs more
>>>>>>>>>>>>>>>>>>>>>>>>>>> iterations. But IMC converges faster, but needs 
>>>>>>>>>>>>>>>>>>>>>>>>>>> longer iterations. Marvin
>>>>>>>>>>>>>>>>>>>>>>>>>>> (who is on the mailing list as well) wrote an 
>>>>>>>>>>>>>>>>>>>>>>>>>>> interesting paper about that
>>>>>>>>>>>>>>>>>>>>>>>>>>> recently:
>>>>>>>>>>>>>>>>>>>>>>>>>>> https://doi.org/10.1021/acs.jctc.2c00665
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023, 09:48 Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <[email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Sir, actually it was given by professor so
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> that, for any arbitrary molecule also I can make 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> itp file to run CG-MD
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> simulation. It is just for my learning purpose.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Suppose, in the tutorial for CG systems
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> topol.top file is there. Like for ibi, imc method 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> etc. topol.top is there.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> But, for any unknown organic CG molecule how can 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> one generate itp file to
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> run CGMD simulation ?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  06/03/23,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 09:45:05 AM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023 at 9:42 AM Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Junghans <[email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 21:00 Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Actually sir, it was asked me to prepare
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> propane.itp (for CG).
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Sorry, can you explain a bit more why you
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> need an itp file? All the files (except for the 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> potentials) are in the
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> tutorial already.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  06/03/23,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 08:29:51 AM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Sat, Jun 3, 2023 at 4:15 AM Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Junghans <[email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 12:34 PM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I have done boltzmann inversion of
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> propane. now, for getting nonbonded 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> parameters what protocol should I
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> follow ? Because, my aim is to generate 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> GROMACS compatible itp file for CG
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> propane so that, can perform CG-MD using 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Votca.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Why can't you use the topology files from
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> the propane tutorial?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> https://github.com/votca/votca/tree/master/csg-tutorials/propane/ibi
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I am asking because most of the CG methods
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> VOTCA implement work on tabulated potentials 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> only.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thanks & Regards
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Kankana Bhattacharjee
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  06/03/23,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 12:00:35 AM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 9:13 PM Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Junghans <[email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Forgot to CC the mailing list, so here is
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> the answer for others.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ---------- Forwarded message ---------
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> From: Christoph Junghans <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Date: Fri, Jun 2, 2023 at 8:22 AM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Subject: Re: [votca]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Regarding-Getting-itp-files:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> To: Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 8:06 AM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Sir, after boltzmann inversion,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> bond.dist.ib, angle.dist.ib and 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> bond.pot.ib, angle.pot.ib (after using tab
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> command)files generated. So, where will I 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> get sigma, epsilon values
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I am not a 100% sure what you mean.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> When you use nbfunc=1 and comb-rule=1,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> the sigma and epsilon column in the 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> atomtypes block are actually the C6 and
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> C12 parameters (i.e. the prefactor of 1/r**6 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> and 1/r**12).
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> But if you then use tables (vdwtype =
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> user in the mdp file), 1/r**6 and 1/r**12 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> gets replaced with the function
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> from files you provided (e.g. obtained by 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> boltzmann inversion) And hence in
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> VOTCA we just set C6=C12=1, so that the 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> table is used unmodified.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Hope that helps,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  06/02/23,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 07:34:36 PM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 7:27 PM Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Junghans <[email protected]> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 7:44 AM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Sir, I would like to know how sigma,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> epsilon and default section was decided 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> fpr CG-MD  simulation ?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> no sigma and epsilon, it is all
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> tabulated interactions.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> see
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> https://github.com/votca/votca/blob/master/csg-tutorials/hexane/ibi_nonbonded/grompp.mdp#L73
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> and those tables you get from boltzmann
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> inversion (or force matching).
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ---------- Forwarded message ---------
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> From: Christoph Junghans <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Date: Fri, Jun 2, 2023 at 7:10 PM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Subject: Re: [votca]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Regarding-Getting-itp-files:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> To: Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 3:40 AM Kankana
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I would like to know that, in votca
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> beads are defined using "A, B" etc. But, 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> how to understand what type of
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> bead it is i.e, polar, aploar, charged 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> etc. Using A, B beads can anyone run
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> GROMACS MD simulation ?
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> In the topol.top of the coarse
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> grained simulation you can set the 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> charge, e.g. see:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> https://github.com/votca/votca/blob/master/csg-tutorials/hexane/ibi_nonbonded/topol.top#L5_L8
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Christoph
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Thanks & Regards
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Kankana Bhattacharjee
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  Sender
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> notified by
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Mailtrack
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>  06/02/23,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 03:05:36 PM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> On Fri, Jun 2, 2023 at 11:06 AM
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Kankana Bhattacharjee <
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [email protected]>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Dear Sir,
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> I have included the individual bond
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> section in propane.xml file. Now, 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> distribution is generating.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> But, for getting non-bonded
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> parameters, which procedure should I 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> have to follow ? I can use either
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Force Matching or Iterative methods to 
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> get non-bonded parameters.
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> [image: Mailtrack]
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=signaturevirality11&;>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

-- 
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/CAHG27e7WohuaNePWwAWeYNP8vZygEYjvpj70hXcGC3xR-kh67w%40mail.gmail.com.

Reply via email to