AMBER Archive (2009)

Subject: Re: [AMBER] parameters for O2 molecule

From: Jeffrey (
Date: Thu Oct 29 2009 - 08:01:44 CDT

Dear Jason,
   Thanks very much for your detailed reply.
   I'd like to know which basis set I should use to obtain the bond parameters from ab initio calculation. Is HF/6-31G* ok?
Many thanks.

--- 09年10月22日,周四, Jason Swails <> 写道:

发件人: Jason Swails <>
主题: Re: [AMBER] parameters for O2 molecule
收件人: "AMBER Mailing List" <>
日期: 2009年10月22日,周四,下午10:19

You can probably use the L-J parameters from the carbonyl O.  The bond
length and vibrational frequency (which will directly give you the bond
force constant) should not be hard to find online or in the CRC (and should
be pretty easy to calculate with an ab-initio package anyway).  Also, since
O2 is neutral and a dipole moment is symmetrically forbidden, I'd say that
the partial charge on each atom in a non-polarizable force field is 0.  You
can use parmchk to create an frcmod file, and that'll show you the
parameters you have to fill in (but be careful to give it a name not
recognized by the force field you're using or it'll give partial charges
based on the pre-defined atom type).  The parameter data file is located at
$AMBERHOME/dat/leap/parm so you can get relevant parameters from your
favorite force field.  Good luck!

All the best,

2009/10/22 Jeffrey <>

> Dear AMBER users,
>   We are going to monitor the transport of O2 molecule by MD simulation.
> But I am not sure whether AMBER FF provides parameters of  the oxygen
> molecule? If not, where can I get the parameters?
> Many thanks.
> Jeffery
> ----
> __________________________________________________
> 赶快注册雅虎超大容量免费邮箱?
> _______________________________________________
> AMBER mailing list

Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Graduate Student
AMBER mailing list

___________________________________________________________ 好玩贺卡等你发,邮箱贺卡全新上线! _______________________________________________ AMBER mailing list