Re: [AMBER] hexagonal unit cell for membrane simulation with NAMD

From: Brian Radak <>
Date: Fri, 8 Nov 2013 09:32:14 -0500

I'm not an expert on non-orthorhombic boxes, but I did once get AMBER to
work correctly with a rhombic dodecahedron. The solution was to:

1) read a pdb of the water box into leap
2) use "setBox" to set dummy box information
3) edit the box vectors in the inpcrd by hand
4) edit IFBOX to be equal to 3 by hand in the prmtop

I checked all of this by using the "image familiar" command in
ptraj/cpptraj. However, be aware that AMBER will simulate in the triclinic
representation, which will make it quite difficult to visualize what is
actually going on. No guarantees on that working for sure though.


On Thu, Nov 7, 2013 at 7:37 PM, Thomas Evangelidis <>wrote:

> Dear list,
> Does anyone know if I can create a hexagonal membrane patch with lipid11 ff
> to run a simulation in a hexagonal unit cell with NAMD? CHARMM-GUI provides
> an option to create a hexagonal box but I think in the .prmtop file there
> is a pointer named "IFBOX" which determines the box size, and since AMBER
> does not support that kind of shape I suspect there might be a conflict
> with the cell basis vectors in the NAMD configuration file. The same may
> apply with the lattice vectors written in the last line of .inpcrd.
> Other than that, the lipid11 tutorial states that I should choose a
> Rectangular box type in CHARMM-GUI since this geometry is required for
> constant surface tension. Will there be a problem if I choose a Hexagonal
> box type and run the simulation in NAMD provided that I define correctly
> the cell basis vectors?
> Could anyone please clarify these points?
> thanks,
> Thomas
> --
> ======================================================================
> Thomas Evangelidis
> PhD student
> University of Athens
> Faculty of Pharmacy
> Department of Pharmaceutical Chemistry
> Panepistimioupoli-Zografou
> 157 71 Athens
> email:
> website:
> _______________________________________________
> AMBER mailing list

================================ Current Address =======================
 Brian Radak                                             :     BioMaPS
Institute for Quantitative Biology
 PhD candidate - York Research Group       :     Rutgers, The State
University of New Jersey
 University of Minnesota - Twin Cities         :     Center for Integrative
Proteomics Room 308
 Graduate Program in Chemical Physics     :     174 Frelinghuysen Road,
 Department of Chemistry                          :     Piscataway, NJ
08854-8066                                 :
Sorry for the multiple e-mail addresses, just use the institute appropriate
AMBER mailing list
Received on Fri Nov 08 2013 - 07:00:03 PST
Custom Search