Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

I am building a computational model of ellipsoidal cell network formation and I would like to use a particle order parameter to study my model's behavior. I have come across this article Steinhardt, Nelson, Ronchetti, PRB (1983), that was apparently the one that started everything, but the math is quite exotic for me. I have found that the formula:

$$ S = 1/3<cos^2θ - 1> $$

is the order parameter for a 3-D system. This is convenient as S = 0 means no order, while S = 1 means complete alignment. A colleague at the univ has mentioned that for a 2-D particle system the previous equation for S becomes:

$$ S = cos2θ $$

but I don't think this can be right, as in the following example: 2 particles have a 90 degree angle between them, so they are completely unaligned. cos2θ = cos2*90 = -1 which would imply they are completely aligned, which is wrong. Also, the range is [-1:1] instead of 0 < S < 1.

I would really appreciate it, if anyone could point me to the right direction.

share|improve this question
add comment

2 Answers

up vote 0 down vote accepted

I met a paper published in PRL using the orietational order parameter and global bond_orietational order parameter. Although the structures in this paper are about the short chain alkanes, i think the definitions of the two kind of order parameters may provide you some insights into your question. enter image description here

share|improve this answer
    
I have looked into the article you posted and it seems that they are also doing this for 3D molecules. So the formula would hold for 3D. I cannot find a 2D formula anywhere and I am thinking that maybe the same could be used for 2D as well. I need to be careful with the assumptions made for 3D in the original article I posted. I am sorry I cannot upvote you. I need more reputation. :) I am still waiting for other answers. If I don't get any today, I will accept yours. –  Dima1982 May 22 '13 at 14:53
    
It does not matter. We all hope the questions can be desirably solved and we can both learn from the questions and answers. –  Legend_Dyson May 23 '13 at 1:48
add comment

In this article PhysRevE.67.031608 an order parameter is given for a 2D system. Originally, it was introduced by Halperin and Nelson. It is the absolute value of the sixth Fourier component of the bond angle distribution function, which is constant in the isotropic fluid and consists of six equally spaced peaks in the solid phase. It is given by:

$$ \Phi_6 = \Big\langle \frac{1}{N} \sum_{m=1}^{N}\frac{1}{N_b}\sum_{n=1}^{N_b}e^{6i\theta_{mn}} \Big\rangle $$

The angular brackets indicate the configurational average and $\theta_{mn}$ is the angle between some fixed axis (e.g. x or y) and the bond joining the m-th particle with a neighboring n-th particle. $N_b$ is the number of particle-neighbor bonds.

The local bond-orientational order parameter is:

$$ \phi_6 =\Big\langle \frac{1}{N_b}\sum_{n=1}^{N_b}e^{6i\theta_{mn}} \Big\rangle $$

that is the order parameter of a single particle. $|\Phi_6|$ and $|\phi_6|$ produce a number between 0 and 1 that represents the orientational order.

I guess this answers my question.

The only question I have now is if the input angle $\theta_{mn}$ needs to be in radians or degrees. I have tried some mock examples by hand and I couldn't understand my own results! Any help would be appreciated!

EDIT: I have come to this answer by reading about the "hexatic" phase of fluids, nematic liquid crystals and the two dimensional melting. As to how feasible these are to use in my research, remains to be seen.

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.