Re: [gmx-users] electron density

2011-05-26 Thread Igor Marques
thanks mate!

  Igor Marques


On Thu, May 26, 2011 at 6:12 PM, Justin A. Lemkul  wrote:

>
>
> Igor Marques wrote:
>
>> dear all,
>>
>> regarding the matter of electron density, can anyone clarify me on the
>> output units?
>> from this discussion i understand the values are presented in
>> electrons/A^3 - is this correct?
>>
>>
> No.  It is electrons/nm^3.  Standard units are defined in Chapter 2 of the
> manual (and printed in the code, of course).
>
> -Justin
>
>  best regards,
>> igor
>>
>>
>> Igor Marques
>>
>> Molecular Modeling Group – Univ. of Aveiro
>> http://molecular-modeling.dq.ua.pt/
>>
>>
>>
>> On Sat, Jan 1, 2011 at 9:01 PM, Jianhui Tian > jianhuit...@gmail.com>> wrote:
>>
>>Problem solved.
>>The "2" in the text is just a typo when I sent the email.
>>So actually there is no partial charge in the electrons.dat file.
>> Thanks, Justin.
>>
>>Jianhui
>>
>>Date: Sat, 01 Jan 2011 13:56:45 -0500
>>
>>From: "Justin A. Lemkul" mailto:jalem...@vt.edu>>
>>Subject: Re: [gmx-users] electron density
>>To: Discussion list for GROMACS users ><mailto:gmx-users@gromacs.org>>
>>Message-ID: <4d1f78ed.5070...@vt.edu <mailto:4d1f78ed.5070...@vt.edu>>
>>
>>
>>Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>
>>
>>
>>Jianhui Tian wrote:
>> > Date: Fri, 31 Dec 2010 19:08:50 -0500
>> > From: "Justin A. Lemkul" ><mailto:jalem...@vt.edu> <mailto:jalem...@vt.edu
>>
>><mailto:jalem...@vt.edu>>>
>> > Subject: Re: [gmx-users] electron density
>> > To: Discussion list for GROMACS users ><mailto:gmx-users@gromacs.org>
>> > <mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>
>>
>> >
>> > Message-ID: <4d1e7092.7060...@vt.edu
>><mailto:4d1e7092.7060...@vt.edu> <mailto:4d1e7092.7060...@vt.edu
>>
>><mailto:4d1e7092.7060...@vt.edu>>>
>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>> >
>> >
>> >
>> > Jianhui Tian wrote:
>> >>  Hi gmx users,
>> >>
>> >>  I tried to calculate both the mass density and electron density
>>for the
>> >>  water layer of a membrane simulation. The mass density is 1000
>>kg/m^3,
>> >>  however the electron density I got is smaller than about 0.33
>>e/A^3. I
>> >
>> >>  How much smaller?
>> > I got about 0.27 e/A^3.
>> >
>> >>  am using a CHARMM force field and the CHARMM version of TIP3P
>>water. Am
>> >>  I doing anything wrong when calculating the electron density?
>>Thanks a
>> > lot.
>> >>
>> >
>> >>  Possibly, but without seeing what you used for input, information
>> > about your
>> >>  simulation, etc it's impossible to say.
>> > In the electrons.dat file, I had
>> > 2
>> > OW = 8.834
>> > HW1 = 0.583
>> > HW2 = 0.583
>> > And then used g_density ... -ei electrons.dat -dens electron ...
>> >
>>
>>There are several problems.  The first line of electrons.dat
>>indicates there are
>>only two unique atom names, but clearly you're supplying three.
>> Second, this
>>file is read in such that at integral number of electrons is
>>detected on each
>>line.  So, in your case, two lines are being read, one that detects
>>8 electrons,
>>and another that detects zero.
>>
>>-Justin
>>
>>--
>>gmx-users mailing listgmx-users@gromacs.org
>><mailto:gmx-users@gromacs.org>
>>http://lists.gromacs.org/mailman/listinfo/gmx-users
>>Please search the archive at
>>http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
>>Please don't post (un)subscribe requests to the list. Use the
>>www interface or send it to gmx-users-requ...@gromacs.org
>><mailto:gmx-users-requ...@gromacs.org>.
>>
>>Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>>
>>
>>
> --
> 
>
> Justin A. Lemkul
>

Re: [gmx-users] electron density

2011-05-26 Thread Justin A. Lemkul



Igor Marques wrote:

dear all,

regarding the matter of electron density, can anyone clarify me on the 
output units?
from this discussion i understand the values are presented in 
electrons/A^3 - is this correct?




No.  It is electrons/nm^3.  Standard units are defined in Chapter 2 of the 
manual (and printed in the code, of course).


-Justin


best regards,
igor


Igor Marques

Molecular Modeling Group – Univ. of Aveiro
http://molecular-modeling.dq.ua.pt/



On Sat, Jan 1, 2011 at 9:01 PM, Jianhui Tian <mailto:jianhuit...@gmail.com>> wrote:


Problem solved.
The "2" in the text is just a typo when I sent the email.
So actually there is no partial charge in the electrons.dat file. 
Thanks, Justin.


Jianhui

Date: Sat, 01 Jan 2011 13:56:45 -0500

From: "Justin A. Lemkul" mailto:jalem...@vt.edu>>
    Subject: Re: [gmx-users] electron density
To: Discussion list for GROMACS users mailto:gmx-users@gromacs.org>>
Message-ID: <4d1f78ed.5070...@vt.edu <mailto:4d1f78ed.5070...@vt.edu>>

Content-Type: text/plain; charset=ISO-8859-1; format=flowed



Jianhui Tian wrote:
 > Date: Fri, 31 Dec 2010 19:08:50 -0500
 > From: "Justin A. Lemkul" mailto:jalem...@vt.edu> <mailto:jalem...@vt.edu
<mailto:jalem...@vt.edu>>>
 > Subject: Re: [gmx-users] electron density
 > To: Discussion list for GROMACS users mailto:gmx-users@gromacs.org>
 > <mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>
 >
 > Message-ID: <4d1e7092.7060...@vt.edu
<mailto:4d1e7092.7060...@vt.edu> <mailto:4d1e7092.7060...@vt.edu
<mailto:4d1e7092.7060...@vt.edu>>>
 > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
 >
 >
 >
 > Jianhui Tian wrote:
 >>  Hi gmx users,
 >>
 >>  I tried to calculate both the mass density and electron density
for the
 >>  water layer of a membrane simulation. The mass density is 1000
kg/m^3,
 >>  however the electron density I got is smaller than about 0.33
e/A^3. I
 >
 >>  How much smaller?
 > I got about 0.27 e/A^3.
 >
 >>  am using a CHARMM force field and the CHARMM version of TIP3P
water. Am
 >>  I doing anything wrong when calculating the electron density?
Thanks a
 > lot.
 >>
 >
 >>  Possibly, but without seeing what you used for input, information
 > about your
 >>  simulation, etc it's impossible to say.
 > In the electrons.dat file, I had
 > 2
 > OW = 8.834
 > HW1 = 0.583
 > HW2 = 0.583
 > And then used g_density ... -ei electrons.dat -dens electron ...
 >

There are several problems.  The first line of electrons.dat
indicates there are
only two unique atom names, but clearly you're supplying three.
 Second, this
file is read in such that at integral number of electrons is
detected on each
line.  So, in your case, two lines are being read, one that detects
8 electrons,
and another that detects zero.

-Justin

--
gmx-users mailing listgmx-users@gromacs.org
<mailto:gmx-users@gromacs.org>
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org
<mailto:gmx-users-requ...@gromacs.org>.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] electron density

2011-05-26 Thread Igor Marques
dear all,

regarding the matter of electron density, can anyone clarify me on the
output units?
from this discussion i understand the values are presented in electrons/A^3
- is this correct?

best regards,
igor


Igor Marques

Molecular Modeling Group – Univ. of Aveiro
http://molecular-modeling.dq.ua.pt/



On Sat, Jan 1, 2011 at 9:01 PM, Jianhui Tian  wrote:

> Problem solved.
> The "2" in the text is just a typo when I sent the email.
> So actually there is no partial charge in the electrons.dat file.
> Thanks, Justin.
>
> Jianhui
>
> Date: Sat, 01 Jan 2011 13:56:45 -0500
>
> From: "Justin A. Lemkul" 
> Subject: Re: [gmx-users] electron density
> To: Discussion list for GROMACS users 
> Message-ID: <4d1f78ed.5070...@vt.edu>
>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>
>
> Jianhui Tian wrote:
> > Date: Fri, 31 Dec 2010 19:08:50 -0500
> > From: "Justin A. Lemkul" mailto:jalem...@vt.edu>>
> > Subject: Re: [gmx-users] electron density
> > To: Discussion list for GROMACS users  > <mailto:gmx-users@gromacs.org>
> >
> > Message-ID: <4d1e7092.7060...@vt.edu <mailto:4d1e7092.7060...@vt.edu>>
> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> >
> >
> >
> > Jianhui Tian wrote:
> >>  Hi gmx users,
> >>
> >>  I tried to calculate both the mass density and electron density for the
> >>  water layer of a membrane simulation. The mass density is 1000 kg/m^3,
> >>  however the electron density I got is smaller than about 0.33 e/A^3. I
> >
> >>  How much smaller?
> > I got about 0.27 e/A^3.
> >
> >>  am using a CHARMM force field and the CHARMM version of TIP3P water. Am
> >>  I doing anything wrong when calculating the electron density? Thanks a
> > lot.
> >>
> >
> >>  Possibly, but without seeing what you used for input, information
> > about your
> >>  simulation, etc it's impossible to say.
> > In the electrons.dat file, I had
> > 2
> > OW = 8.834
> > HW1 = 0.583
> > HW2 = 0.583
> > And then used g_density ... -ei electrons.dat -dens electron ...
> >
>
> There are several problems.  The first line of electrons.dat indicates
> there are
> only two unique atom names, but clearly you're supplying three.  Second,
> this
> file is read in such that at integral number of electrons is detected on
> each
> line.  So, in your case, two lines are being read, one that detects 8
> electrons,
> and another that detects zero.
>
> -Justin
>
> --
> gmx-users mailing listgmx-users@gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-requ...@gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] electron density

2011-01-01 Thread Jianhui Tian
Problem solved.
The "2" in the text is just a typo when I sent the email.
So actually there is no partial charge in the electrons.dat file.
Thanks, Justin.

Jianhui

Date: Sat, 01 Jan 2011 13:56:45 -0500
From: "Justin A. Lemkul" 
Subject: Re: [gmx-users] electron density
To: Discussion list for GROMACS users 
Message-ID: <4d1f78ed.5070...@vt.edu>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed



Jianhui Tian wrote:
> Date: Fri, 31 Dec 2010 19:08:50 -0500
> From: "Justin A. Lemkul" mailto:jalem...@vt.edu>>
> Subject: Re: [gmx-users] electron density
> To: Discussion list for GROMACS users  <mailto:gmx-users@gromacs.org>
>
> Message-ID: <4d1e7092.7060...@vt.edu <mailto:4d1e7092.7060...@vt.edu>>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
>
>
> Jianhui Tian wrote:
>>  Hi gmx users,
>>
>>  I tried to calculate both the mass density and electron density for the
>>  water layer of a membrane simulation. The mass density is 1000 kg/m^3,
>>  however the electron density I got is smaller than about 0.33 e/A^3. I
>
>>  How much smaller?
> I got about 0.27 e/A^3.
>
>>  am using a CHARMM force field and the CHARMM version of TIP3P water. Am
>>  I doing anything wrong when calculating the electron density? Thanks a
> lot.
>>
>
>>  Possibly, but without seeing what you used for input, information
> about your
>>  simulation, etc it's impossible to say.
> In the electrons.dat file, I had
> 2
> OW = 8.834
> HW1 = 0.583
> HW2 = 0.583
> And then used g_density ... -ei electrons.dat -dens electron ...
>

There are several problems.  The first line of electrons.dat indicates there
are
only two unique atom names, but clearly you're supplying three.  Second,
this
file is read in such that at integral number of electrons is detected on
each
line.  So, in your case, two lines are being read, one that detects 8
electrons,
and another that detects zero.

-Justin
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] electron density

2011-01-01 Thread Justin A. Lemkul



Jianhui Tian wrote:

Date: Fri, 31 Dec 2010 19:08:50 -0500
From: "Justin A. Lemkul" mailto:jalem...@vt.edu>>
Subject: Re: [gmx-users] electron density
To: Discussion list for GROMACS users <mailto:gmx-users@gromacs.org>>

Message-ID: <4d1e7092.7060...@vt.edu <mailto:4d1e7092.7060...@vt.edu>>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed



Jianhui Tian wrote:

 Hi gmx users,

 I tried to calculate both the mass density and electron density for the
 water layer of a membrane simulation. The mass density is 1000 kg/m^3,
 however the electron density I got is smaller than about 0.33 e/A^3. I



 How much smaller?

I got about 0.27 e/A^3.


 am using a CHARMM force field and the CHARMM version of TIP3P water. Am
 I doing anything wrong when calculating the electron density? Thanks a 

lot.




 Possibly, but without seeing what you used for input, information 

about your

 simulation, etc it's impossible to say.
In the electrons.dat file, I had 
2 
OW = 8.834

HW1 = 0.583
HW2 = 0.583
And then used g_density ... -ei electrons.dat -dens electron ...



There are several problems.  The first line of electrons.dat indicates there are 
only two unique atom names, but clearly you're supplying three.  Second, this 
file is read in such that at integral number of electrons is detected on each 
line.  So, in your case, two lines are being read, one that detects 8 electrons, 
and another that detects zero.


-Justin



-Justin



--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] electron density

2011-01-01 Thread Jianhui Tian
Date: Fri, 31 Dec 2010 19:08:50 -0500
From: "Justin A. Lemkul" 
Subject: Re: [gmx-users] electron density
To: Discussion list for GROMACS users 
Message-ID: <4d1e7092.7060...@vt.edu>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed



Jianhui Tian wrote:
> Hi gmx users,
>
> I tried to calculate both the mass density and electron density for the
> water layer of a membrane simulation. The mass density is 1000 kg/m^3,
> however the electron density I got is smaller than about 0.33 e/A^3. I

> How much smaller?
I got about 0.27 e/A^3.

> am using a CHARMM force field and the CHARMM version of TIP3P water. Am
> I doing anything wrong when calculating the electron density? Thanks a
lot.
>

> Possibly, but without seeing what you used for input, information about
your
> simulation, etc it's impossible to say.
In the electrons.dat file, I had
2
OW = 8.834
HW1 = 0.583
HW2 = 0.583
And then used g_density ... -ei electrons.dat -dens electron ...


-Justin
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] electron density

2010-12-31 Thread Justin A. Lemkul



Jianhui Tian wrote:

Hi gmx users,

I tried to calculate both the mass density and electron density for the 
water layer of a membrane simulation. The mass density is 1000 kg/m^3, 
however the electron density I got is smaller than about 0.33 e/A^3. I 


How much smaller?

am using a CHARMM force field and the CHARMM version of TIP3P water. Am 
I doing anything wrong when calculating the electron density? Thanks a lot.




Possibly, but without seeing what you used for input, information about your 
simulation, etc it's impossible to say.


-Justin

--


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] electron density

2010-12-31 Thread Jianhui Tian
Hi gmx users,

I tried to calculate both the mass density and electron density for the
water layer of a membrane simulation. The mass density is 1000 kg/m^3,
however the electron density I got is smaller than about 0.33 e/A^3. I am
using a CHARMM force field and the CHARMM version of TIP3P water. Am I doing
anything wrong when calculating the electron density? Thanks a lot.

Jianhui
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Electron Density

2010-01-06 Thread David van der Spoel

afsaneh maleki wrote:

Hi,

I want to obtain electron density for solvent in z direction for bilayer.

I used following options:

 

] g_density -f md.xtc -s md.tpr  -d z  -dens electron  -o 
electrondens.xvg  -ei electrons.dat  -symm -n index.ndx


Atomtapes in md.gro are OH,HW1,HW2

And “electrons.dat”  file contains:

3

  OW= 8

 HW1=1

 HW2=1

 


But I get error:

Fatal error:

Invalid line in datafile at line 1

 


What is my error?


remove = sign



Best wishes,

 Afsaneh Maleki




--
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205. Fax: +4618511755.
sp...@xray.bmc.uu.sesp...@gromacs.org   http://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] Electron Density

2010-01-06 Thread afsaneh maleki
Hi,

I want to obtain electron density for solvent in z direction for bilayer.

I used following options:



] g_density -f md.xtc -s md.tpr  -d z  -dens electron  -o electrondens.xvg  -ei
electrons.dat  -symm -n index.ndx

Atomtapes in md.gro are OH,HW1,HW2

And “electrons.dat”  file contains:

3

  OW= 8

 HW1=1

 HW2=1



But I get error:

Fatal error:

Invalid line in datafile at line 1



What is my error?


Best wishes,
 Afsaneh Maleki
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

[gmx-users] Electron density over area

2008-11-09 Thread anirban polley
Hi,
 I know that electron density can be calculated by g_density. But I want
to calculate electron density over area of the membrane so that I can
analyze the fluctuation of the electron density with the area of the
membrane. Could you tell me how to do it?
Thanking you,
Anirban
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php


Re: [gmx-users] electron density .dat file

2006-11-22 Thread David van der Spoel

Jian Dai wrote:

Hello, all:
I'm not sure how to make the .dat file in calculating the electron 
density, for example,

If I have cholesterol.itp like:
;   nrtype   resnr  residuatomcgnrcharge   ; total charge
 1 CH3   1CHOL C11 0
 2  CB   1CHOL C2 2 0
 3 CH2   1CHOL C33 0
 4 CH2   1CHOL C44 0
 5 CH1   1CHOL C55 0.14 ; charges adopted
 6  OA   1CHOL O6 5 -0.54 ; from Retinol in
 7  HO   1CHOL H  5  0.40 ; ffgmx.rtp
 ...

shall I make electron_CHO.dat like:
29
C1 =6
C2 =6
C3 =6
C4 =6
C5 =6.14
O6 = 7.46
H = 1.4

or like:
29
C1 =6
C2 =6
C3 =6
C4 =6
C5 = 5.76
O6 = 8.54
H =  0.6


I mean, shall I add or subtract the partial charge?
Thank you in advance and happy Thanksgiving.
Jian

What do you want to do with this?

--
David.

David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  75124 Uppsala, Sweden
phone:  46 18 471 4205  fax: 46 18 511 755
[EMAIL PROTECTED]   [EMAIL PROTECTED]   http://folding.bmc.uu.se

___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]

Can't post? Read http://www.gromacs.org/mailing_lists/users.php


[gmx-users] electron density .dat file

2006-11-22 Thread Jian Dai

Hello, all:
I'm not sure how to make the .dat file in calculating the electron density,
for example,
If I have cholesterol.itp like:
;   nrtype   resnr  residuatomcgnrcharge   ; total charge
1 CH3   1CHOL C11 0
2  CB   1CHOL C2 2 0
3 CH2   1CHOL C33 0
4 CH2   1CHOL C44 0
5 CH1   1CHOL C55 0.14 ; charges adopted
6  OA   1CHOL O6 5 -0.54 ; from Retinol in
7  HO   1CHOL H  5  0.40 ; ffgmx.rtp
...

shall I make electron_CHO.dat like:
29
C1 =6
C2 =6
C3 =6
C4 =6
C5 =6.14
O6 = 7.46
H = 1.4

or like:
29
C1 =6
C2 =6
C3 =6
C4 =6
C5 = 5.76
O6 = 8.54
H =  0.6


I mean, shall I add or subtract the partial charge?
Thank you in advance and happy Thanksgiving.
Jian
___
gmx-users mailing listgmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php