| 1 |
< |
\documentclass[11pt]{article} |
| 2 |
< |
\usepackage{amsmath} |
| 3 |
< |
\usepackage{amssymb} |
| 1 |
> |
\documentclass[12pt]{article} |
| 2 |
|
\usepackage{endfloat} |
| 3 |
+ |
\usepackage{amsmath} |
| 4 |
|
\usepackage{epsf} |
| 5 |
< |
\usepackage{berkeley} |
| 5 |
> |
\usepackage{times} |
| 6 |
> |
\usepackage{setspace} |
| 7 |
> |
\usepackage{tabularx} |
| 8 |
|
\usepackage{graphicx} |
| 9 |
|
\usepackage[ref]{overcite} |
| 9 |
– |
\usepackage{tabularx} |
| 10 |
|
\pagestyle{plain} |
| 11 |
|
\pagenumbering{arabic} |
| 12 |
|
\oddsidemargin 0.0cm \evensidemargin 0.0cm |
| 24 |
|
University of Notre Dame \\ |
| 25 |
|
Notre Dame, Indiana 46556} |
| 26 |
|
|
| 27 |
+ |
\date{\today} |
| 28 |
+ |
|
| 29 |
|
\maketitle |
| 30 |
|
|
| 31 |
|
\begin{abstract} |
| 47 |
|
|
| 48 |
|
\section{Introduction} |
| 49 |
|
\label{Int} |
| 50 |
< |
\indent |
| 50 |
> |
|
| 51 |
|
Fully hydrated lipids will aggregate spontaneously to form bilayers |
| 52 |
|
which exhibit a variety of phases according to temperature and |
| 53 |
|
composition. Among these phases, a periodic rippled |
| 81 |
|
\section{Model and calculation method} |
| 82 |
|
\label{Mod} |
| 83 |
|
|
| 84 |
< |
The model used in our simulations is shown in Fig. \ref{fmod1} and Fig. \ref{fmod2}. |
| 84 |
> |
The model used in our simulations is shown in Fig. \ref{fmod1} and |
| 85 |
> |
Fig. \ref{fmod2}. |
| 86 |
> |
|
| 87 |
|
\begin{figure} |
| 88 |
|
\centering |
| 89 |
|
\includegraphics[width=\linewidth]{picture/lattice.eps} |
| 93 |
|
potentials.} |
| 94 |
|
\label{fmod1} |
| 95 |
|
\end{figure} |
| 96 |
+ |
|
| 97 |
|
\begin{figure} |
| 98 |
+ |
\centering |
| 99 |
|
\includegraphics[width=\linewidth]{picture/xyz.eps} |
| 100 |
|
\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 |
| 101 |
|
the initial x-y plane, $\theta_i$ is the angle that the dipole is away |
| 103 |
|
the dipole on x-y plane with the x axis.} |
| 104 |
|
\label{fmod2} |
| 105 |
|
\end{figure} |
| 106 |
+ |
|
| 107 |
|
The lipids are represented by the simple point-dipole. During the |
| 108 |
|
simulations, dipoles are locked (in the x-y plane) to lattice points |
| 109 |
|
of hexagonal (or distorted) lattice. Each dipole can move freely out |
| 227 |
|
%\end{figure} |
| 228 |
|
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. |
| 229 |
|
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}. |
| 230 |
+ |
|
| 231 |
|
\begin{figure} |
| 232 |
< |
\begin{center} |
| 232 |
> |
\centering |
| 233 |
|
\includegraphics[width=\linewidth]{picture/hexorderpara.eps} |
| 234 |
< |
\caption{ The orderparameter $P_2$ vs temperature T at hexagonal lattice.\label{t-op}} |
| 235 |
< |
\end{center} |
| 234 |
> |
\caption{ The orderparameter $P_2$ vs temperature T at hexagonal |
| 235 |
> |
lattice.} |
| 236 |
> |
\label{t-op} |
| 237 |
|
\end{figure} |
| 238 |
+ |
|
| 239 |
|
The $P_2 \approx 0.9$ for $T = 100$ implies that the system is in a |
| 240 |
|
highly ordered state. As the temperature increases, the order |
| 241 |
|
parameter is decreasing gradually before $T = 300$, from $T = 310$ the |
| 258 |
|
the further insight into this problem. At first, we fixed the value |
| 259 |
|
of $s = 7$, and vary $k_r$, the results are shown in |
| 260 |
|
Fig. \ref{kr-a-hf}. |
| 261 |
+ |
|
| 262 |
|
\begin{figure} |
| 263 |
< |
\begin{center} |
| 264 |
< |
\includegraphics[width=\linewidth]{picture/kr_amplitude.eps} |
| 265 |
< |
\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}} |
| 266 |
< |
\end{center} |
| 263 |
> |
\centering |
| 264 |
> |
\includegraphics[width=\linewidth]{picture/kr_amplitude.eps} |
| 265 |
> |
\caption{ The amplitude $A$ of the ripples at $T = 300$ vs $k_r$ for |
| 266 |
> |
hexagonal lattice. The height fluctuation $h_f$ vs $k_r$ is shown |
| 267 |
> |
inset for the same situation.} |
| 268 |
> |
\label{kr-a-hf} |
| 269 |
|
\end{figure} |
| 270 |
+ |
|
| 271 |
|
When $k_r < 0.1$, due to the small surface tension part, the dipoles |
| 272 |
|
can go far away from their neighbors, lots of noise make the ripples |
| 273 |
|
undiscernable. So, we start from $k_r = 0.1$. However, even when $k_r |
| 313 |
|
lattice, the distorted hexagonal lattices prefer a particular director |
| 314 |
|
axis due to their anisotropic property. The phase diagram for this |
| 315 |
|
system is shown in Fig. \ref{phase}. |
| 316 |
+ |
|
| 317 |
|
\begin{figure} |
| 318 |
< |
\begin{center} |
| 318 |
> |
\centering |
| 319 |
|
\includegraphics[width=\linewidth]{picture/phase.eps} |
| 320 |
< |
\caption{ The phase diagram with temperature $T$ and lattice variable $\gamma$. The enlarged view near the hexagonal lattice is shown inset.\label{phase}} |
| 321 |
< |
\end{center} |
| 320 |
> |
\caption{ The phase diagram with temperature $T$ and lattice variable |
| 321 |
> |
$\gamma$. The enlarged view near the hexagonal lattice is shown |
| 322 |
> |
inset.} |
| 323 |
> |
\label{phase} |
| 324 |
|
\end{figure} |
| 325 |
+ |
|
| 326 |
|
$T_c$ is the transition temperature. The hexagonal lattice has the |
| 327 |
|
lowest $T_c$, and $T_c$ goes up with lattice being more |
| 328 |
|
distorted. There is only two phases in our diagram. When we do |
| 367 |
|
$\gamma = 1.875$, $k_r = 0.1$, $T = 260$. Under this situation, the |
| 368 |
|
system reaches the equilibrium very quickly, and the ripples are |
| 369 |
|
fairly stable. The results are shown in Fig. \ref{samplitude}. |
| 370 |
+ |
|
| 371 |
|
\begin{figure} |
| 372 |
< |
\begin{center} |
| 372 |
> |
\centering |
| 373 |
|
\includegraphics[width=\linewidth]{picture/samplitude.eps} |
| 374 |
< |
\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}} |
| 375 |
< |
\end{center} |
| 374 |
> |
\caption{ The amplitude of ripples $A$ at $T = 260$ vs strength of |
| 375 |
> |
dipole $s$ for $\gamma = 1.875$. The orderparameter $P_2$ vs $s$ is |
| 376 |
> |
shown inset at the same situation.} |
| 377 |
> |
\label{samplitude} |
| 378 |
|
\end{figure} |
| 379 |
+ |
|
| 380 |
|
For small $s$, there is no long range ordering in the system, so, we |
| 381 |
|
start from $s = 7$, and we use the rippled state as the initial |
| 382 |
|
configuration for all the simulations to reduce the noise. There is no |
