--- trunk/ripplePaper/ripple.tex 2005/03/28 20:54:31 2142 +++ trunk/ripplePaper/ripple.tex 2005/05/31 20:43:32 2254 @@ -1,12 +1,12 @@ -\documentclass[11pt]{article} -\usepackage{amsmath} -\usepackage{amssymb} +\documentclass[12pt]{article} \usepackage{endfloat} +\usepackage{amsmath} \usepackage{epsf} -\usepackage{berkeley} +\usepackage{times} +\usepackage{setspace} +\usepackage{tabularx} \usepackage{graphicx} \usepackage[ref]{overcite} -\usepackage{tabularx} \pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm \evensidemargin 0.0cm @@ -18,34 +18,43 @@ \begin{document} -\title{Ripple Phase of the Lipid Bilayers: A Monte Carlo Simulation} +\title{Symmetry breaking and the Ripple phase} \author{Xiuquan Sun and J. Daniel Gezelter\footnote{Corresponding author. Email: gezelter@nd.edu} \\ Department of Chemistry and Biochemistry \\ University of Notre Dame \\ Notre Dame, Indiana 46556} +\date{\today} + \maketitle \begin{abstract} -The molecular explanation for the origin and properties of the ripple -phase is addressed in this paper. A model which contains the surface -tension and dipole-dipole interactions is used to describe the -potential for a monolayer of simple point dipoles. The simulations are -carried out using Monte Carlo method. It is shown asymmetry of the -translational freedom of the dipoles breaks the symmetry of the -hexagonal lattice and allow antiferroelectric ordering of the -dipoles. The existence of the ripples only depends on the dipolar -property of the system. The structure of the ripples is affected by -surface tension. Only close to the hexagonal lattice, can the ripple -phase be reached. Surface has the lowest transition temperature on -hexagonal lattice elucidates the reason of the existence of the ripple -phase in organism. A mechanism for the phase transition of the lipid -bilayer is proposed. +The ripple phase in phosphatidylcholine (PC) bilayers has never been +explained. Our group has developed some simple (XYZ) spin-lattice +models that allow spins to vary their elevation as well as their +orientation. The extra degree of freedom allows hexagonal lattices of +spins to find states that break out of the normally frustrated +randomized states and are stabilized by long-range anti-ferroelectric +ordering. To break out of the frustrated states, the spins must form +``rippled'' phases that make the lattices effectively non-hexagonal. Our +XYZ models contain surface tension and dipole-dipole interactions to +describe the interaction potential for monolayers and bilayers of +model lipid molecules. The existence of the ripples depends primarily +on the strength and lattice spacing of the dipoles, while the shape +(wavelength and amplitude) of the ripples is only weakly sensitive to +the applied surface tension. Additionally, the wave vector for the +ripple is always perpendicular to the director axis for the +dipoles. Non-hexagonal lattices of dipoles are not inherently +frustrated, and are therefore less likely to form ripple phases +because they can easily form low-energy anti-ferroelectric states. +Indeed, we see that the dipolar order-disorder transition is +substantially lower for hexagonal lattices and the ordered phase for +this lattice is clearly rippled. \end{abstract} \section{Introduction} \label{Int} -\indent + Fully hydrated lipids will aggregate spontaneously to form bilayers which exhibit a variety of phases according to temperature and composition. Among these phases, a periodic rippled @@ -79,7 +88,9 @@ The model used in our simulations is shown in Fig. \re \section{Model and calculation method} \label{Mod} -The model used in our simulations is shown in Fig. \ref{fmod1} and Fig. \ref{fmod2}. +The model used in our simulations is shown in Fig. \ref{fmod1} and +Fig. \ref{fmod2}. + \begin{figure} \centering \includegraphics[width=\linewidth]{picture/lattice.eps} @@ -89,7 +100,9 @@ potentials.} potentials.} \label{fmod1} \end{figure} + \begin{figure} +\centering \includegraphics[width=\linewidth]{picture/xyz.eps} \caption{The 6 coordinates describing the state of a 2-dipole system in our extended X-Y-Z model. $z_i$ is the height of dipole $i$ from the initial x-y plane, $\theta_i$ is the angle that the dipole is away @@ -97,6 +110,7 @@ The lipids are represented by the simple point-dipole. the dipole on x-y plane with the x axis.} \label{fmod2} \end{figure} + The lipids are represented by the simple point-dipole. During the simulations, dipoles are locked (in the x-y plane) to lattice points of hexagonal (or distorted) lattice. Each dipole can move freely out @@ -220,12 +234,15 @@ We use the last configuration of this simulation as th %\end{figure} From the results of the simulation at $T = 300$ for hexagonal lattice, the system is in an antiferroelectric state. Every $3$ or $4$ arrows of the dipoles form a strip whose direction is opposite to its neighbors. $P_2$ is about $0.7$. The ripple is formed clearly. The simulation results shows the ripple has equal opportunity to be formed along different directions, this is due to the isotropic property of the hexagonal lattice. We use the last configuration of this simulation as the initial condition to increase the system to $T = 400$ every $10\ kelvin$, at the same time decrease the temperature to $T = 100$ every $10\ kelvin$. The trend that the order parameter varies with temperature is plotted in Fig. \ref{t-op}. + \begin{figure} -\begin{center} +\centering \includegraphics[width=\linewidth]{picture/hexorderpara.eps} -\caption{ The orderparameter $P_2$ vs temperature T at hexagonal lattice.\label{t-op}} -\end{center} +\caption{ The orderparameter $P_2$ vs temperature T at hexagonal +lattice.} +\label{t-op} \end{figure} + The $P_2 \approx 0.9$ for $T = 100$ implies that the system is in a highly ordered state. As the temperature increases, the order parameter is decreasing gradually before $T = 300$, from $T = 310$ the @@ -248,12 +265,16 @@ Fig. \ref{kr-a-hf}. the further insight into this problem. At first, we fixed the value of $s = 7$, and vary $k_r$, the results are shown in Fig. \ref{kr-a-hf}. + \begin{figure} -\begin{center} +\centering \includegraphics[width=\linewidth]{picture/kr_amplitude.eps} -\caption{ The amplitude $A$ of the ripples at $T = 300$ vs $k_r$ for hexagonal lattice. The height fluctuation $h_f$ vs $k_r$ is shown inset for the same situation.\label{kr-a-hf}} -\end{center} +\caption{ The amplitude $A$ of the ripples at $T = 300$ vs $k_r$ for +hexagonal lattice. The height fluctuation $h_f$ vs $k_r$ is shown +inset for the same situation.} +\label{kr-a-hf} \end{figure} + When $k_r < 0.1$, due to the small surface tension part, the dipoles can go far away from their neighbors, lots of noise make the ripples undiscernable. So, we start from $k_r = 0.1$. However, even when $k_r @@ -299,12 +320,16 @@ system is shown in Fig. \ref{phase}. lattice, the distorted hexagonal lattices prefer a particular director axis due to their anisotropic property. The phase diagram for this system is shown in Fig. \ref{phase}. + \begin{figure} -\begin{center} +\centering \includegraphics[width=\linewidth]{picture/phase.eps} -\caption{ The phase diagram with temperature $T$ and lattice variable $\gamma$. The enlarged view near the hexagonal lattice is shown inset.\label{phase}} -\end{center} +\caption{ The phase diagram with temperature $T$ and lattice variable +$\gamma$. The enlarged view near the hexagonal lattice is shown +inset.} +\label{phase} \end{figure} + $T_c$ is the transition temperature. The hexagonal lattice has the lowest $T_c$, and $T_c$ goes up with lattice being more distorted. There is only two phases in our diagram. When we do @@ -349,12 +374,16 @@ fairly stable. The results are shown in Fig. \ref{samp $\gamma = 1.875$, $k_r = 0.1$, $T = 260$. Under this situation, the system reaches the equilibrium very quickly, and the ripples are fairly stable. The results are shown in Fig. \ref{samplitude}. + \begin{figure} -\begin{center} +\centering \includegraphics[width=\linewidth]{picture/samplitude.eps} -\caption{ The amplitude of ripples $A$ at $T = 260$ vs strength of dipole $s$ for $\gamma = 1.875$. The orderparameter $P_2$ vs $s$ is shown inset at the same situation.\label{samplitude}} -\end{center} +\caption{ The amplitude of ripples $A$ at $T = 260$ vs strength of +dipole $s$ for $\gamma = 1.875$. The orderparameter $P_2$ vs $s$ is +shown inset at the same situation.} +\label{samplitude} \end{figure} + For small $s$, there is no long range ordering in the system, so, we start from $s = 7$, and we use the rippled state as the initial configuration for all the simulations to reduce the noise. There is no