Re: [ESPResSo] ESPResSo Digest, Vol 29, Issue 8

[Michael Winokur]

Hello Olaf and all,

I'm sorry but my description wasn't quite clear.  The "dipole" function serves 
only to give a preferred orientation to pairs of bonded oblate particles.  
This orientation would then be used to calculate the canonical dihedral 
interaction between these particles.  (Note: The GB potential has been used in 
modeling the phenyl units in a polystyrene simulation which claimed good 
results.  I note that there was a united atom representation for the CH and 
CH_2 units.)    

My C programming skills are minimal and so simplicity of new modifications is 
key.   Olaf suggested I take a closer look at the FENE  interaction but, after 
thinking it over, I believe there is an even simpler implementation.  The 
standard dihedral interaction is "part p_2 bond bondid p_1 p_3 p_4".  In this 
case p_1 and p_4 are needed for ascertaining the torsion angle across the 
bond.  If one made a conditional on the basis of "part p_2 bond bondid p_2 p_3 
p_3" then whenever the p_1 was set to p_2 or p_4 to p_3 the dihedral 
interaction subroutine could defer to the "dipole orientation" to ascertain 
the needed geometrical constraints.  Is there any reason this should/could not 
be done?

A secondary issue has to do with the subtracted LJ.  My guess is that there is 
no procedure to do the same with the Gay-Berne potential because this is the 
first go at linking GB particles within Espresso.  This all depends on whether 
the subtracted LJ subroutine has been coded to include any GB interactions.

I'm not sure I understand how mixtures of GB particles and spherical LJ 
particles are integrated.   

I also am not sure how the unit vectors u_i and u_j used in the GB interaction 
are established unless they refer back to the dipole orientation.  Am I 
missing something here?

Thanks, again, for any suggestions/feedback.

Michael

On Tuesday 13 October 2009 05:02:03 address@hidden 
wrote:
> Send ESPResSo mailing list submissions to
>       address@hidden
> 
> To subscribe or unsubscribe via the World Wide Web, visit
>       https://fias.uni-frankfurt.de/mailman/listinfo/espresso
> or, via email, send a message with subject or body 'help' to
>       address@hidden
> 
> You can reach the person managing the list at
>       address@hidden
> 
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ESPResSo digest..."
> 
> 
> Today's Topics:
> 
>    1. Re: Coarse Graining with a Gay Berne potential (Olaf Lenz)
> 
> 
> ----------------------------------------------------------------------
> 
> Message: 1
> Date: Mon, 12 Oct 2009 12:21:48 +0200
> From: Olaf Lenz <address@hidden>
> Subject: Re: [ESPResSo] Coarse Graining with a Gay Berne potential
> To: Michael Winokur <address@hidden>
> Cc: address@hidden
> Message-ID: <address@hidden>
> Content-Type: text/plain; charset=UTF-8; format=flowed
> 
> Hi!
> 
> Michael Winokur wrote:
> > I'm now at the point where I'm trying to use the GB potential for
> > representing a phenylene ring in a molecule.  However I do not see how to
> > implement a dihedral interaction between biphenylene pairs.  The simplest
> > method I can imagine is to specify a dipole direction for each GB unit
> > and then use a dipole-dipole interaction for the two particles which maps
> > the more standard dihedral interaction requiring four particle ids.  I do
> > not see anything in the current Espresso interaction list that manages
> > this (although it should not be too difficult to add something like this
> > to the code by someone that is more familiar than I with the source
> > code).  Is there an existing method that I ought use?
> 
> So far, there is no such interaction. To my knowledge, the GB potential
> has not been used in any compounds with bindings, so far, but mostly for
> liquid-crystal-like particles.
> 
> Using the dipole-dipole interaction is not advisable: it is a long-range
> interaction, and therefore much slower than any short-range interaction,
> and it is pretty bad in modeling bindings, would think.
> 
> I guess you will have to implement such a potential yourself. It is not
> very complicated to create a new bond potential yourself. Have a look at
> the implementation of FENE potential and the GB potential itself, and
> you should be able to implement something that works.
> The mathematics of that potential, however, might become tedious.
> 
> Good luck!
> 
>       Olaf
> 
> 
> 
> ------------------------------
> 
> _______________________________________________
> ESPResSo mailing list
> address@hidden
> https://fias.uni-frankfurt.de/mailman/listinfo/espresso
> 
> 
> End of ESPResSo Digest, Vol 29, Issue 8
> ***************************************
> 

-- 
------------------------------------------------------------------------------------------------------
Michael Winokur
Dept. of Physics, University of Wisconsin
1150 University Ave.
Madison, WI 53706
Phone: 608 263 7475 URL: romano.physics.wisc.edu
-----------------------------------------------------------------------------------------------------