40 |
|
|
41 |
|
|
42 |
|
\title{Nitrile vibrations as reporters of field-induced phase |
43 |
< |
transitions in liquid crystals} |
43 |
> |
transitions in 4-cyano-4'-pentylbiphenyl (5CB)} |
44 |
|
\author{James M. Marr} |
45 |
|
\author{J. Daniel Gezelter} |
46 |
|
\email{gezelter@nd.edu} |
58 |
|
\begin{doublespace} |
59 |
|
|
60 |
|
\begin{abstract} |
61 |
< |
Nitrile Stark shift repsonses to electric fields have been used |
62 |
< |
extensively in biology for the probing of local internal fields of |
63 |
< |
structures like proteins and DNA. Intigration of these probes into |
64 |
< |
different areas of interest are important for studing local structure |
65 |
< |
and fields within confined, nanoscopic |
66 |
< |
systems. 4-Cyano-4'-pentylbiphenyl (5CB) is a liquid crystal with a known |
67 |
< |
macroscopic structure reordering from the isotropic to nematic |
68 |
< |
phase with the application of an external |
69 |
< |
field and as the name suggests has an inherent nitrile group. Through |
70 |
< |
simulations of this molecule where application of |
71 |
< |
large, nanoscale external fields were applied, the nitrile was invenstigated |
72 |
< |
as a local field sensor. It was |
73 |
< |
found that while most computational methods for nitrile spectral |
74 |
< |
calculations rely on a correlation between local electric field and |
75 |
< |
the nitrile bond, 5CB did not have an easily obtained |
76 |
< |
correlation. Instead classical calculation through correlation of the |
77 |
< |
cyanide bond displacement in time use enabled to show a spectral |
78 |
< |
change in the formation of a red |
79 |
< |
shifted peak from the main peak as an external field was applied. This indicates |
80 |
< |
that local structure has a larger impact on the nitrile frequency then |
81 |
< |
does the local electric field. By better understanding how nitrile |
82 |
< |
groups respond to local and external fields it will help |
83 |
< |
nitrile groups branch out beyond their biological |
84 |
< |
applications to uses in electronics and surface sciences. |
61 |
> |
4-cyano-4'-pentylbiphenyl (5CB) is a liquid-crystal-forming compound |
62 |
> |
with a terminal nitrile group aligned with the long axis of the |
63 |
> |
molecule. Simulations of condensed-phase 5CB were carried out both |
64 |
> |
with and without applied electric fields to provide an understanding |
65 |
> |
of the the Stark shift of the terminal nitrile group. A |
66 |
> |
field-induced isotropic-nematic phase transition was observed in the |
67 |
> |
simulations, and the effects of this transition on the distribution |
68 |
> |
of nitrile frequencies were computed. Classical bond displacement |
69 |
> |
correlation functions exhibit a $\sim 40 \mathrm{~cm}^{-1}$ red |
70 |
> |
shift of a portion of the main nitrile peak, and this shift was |
71 |
> |
observed only when the fields were large enough to induce |
72 |
> |
orientational ordering of the bulk phase. Our simulations appear to |
73 |
> |
indicate that phase-induced changes to the local surroundings are a |
74 |
> |
larger contribution to the change in the nitrile spectrum than |
75 |
> |
direct field contributions. |
76 |
|
\end{abstract} |
77 |
|
|
78 |
|
\newpage |
79 |
|
|
80 |
|
\section{Introduction} |
81 |
|
|
82 |
+ |
Nitrile groups can serve as very precise electric field reporters via |
83 |
+ |
their distinctive Raman and IR signatures.\cite{Boxer:2009xw} The |
84 |
+ |
triple bond between the nitrogen and the carbon atom is very sensitive |
85 |
+ |
to local field changes and has been observed to have a direct impact |
86 |
+ |
on the peak position within the spectrum. The Stark shift in the |
87 |
+ |
spectrum can be quantified and mapped into a field value that is |
88 |
+ |
impinging upon the nitrile bond. This has been used extensively in |
89 |
+ |
biological systems like proteins and |
90 |
+ |
enzymes.\cite{Tucker:2004qq,Webb:2008kn} |
91 |
+ |
|
92 |
+ |
The response of nitrile groups to electric fields has now been |
93 |
+ |
investigated for a number of small molecules,\cite{Andrews:2000qv} as |
94 |
+ |
well as in biochemical settings, where nitrile groups can act as |
95 |
+ |
minimally invasive probes of structure and |
96 |
+ |
dynamics.\cite{Lindquist:2009fk,Fafarman:2010dq} The vibrational Stark |
97 |
+ |
effect has also been used to study the effects of electric fields on |
98 |
+ |
nitrile-containing self-assembled monolayers at metallic |
99 |
+ |
interfaces.\cite{Oklejas:2002uq,Schkolnik:2012ty} |
100 |
+ |
|
101 |
+ |
Recently 4-cyano-4'-pentylbiphenyl (5CB), a liquid crystalline |
102 |
+ |
molecule with a terminal nitrile group, has seen renewed interest as |
103 |
+ |
one way to impart order on the surfactant interfaces of |
104 |
+ |
nanodroplets,\cite{Moreno-Razo:2012rz} or to drive surface-ordering |
105 |
+ |
that can be used to promote particular kinds of |
106 |
+ |
self-assembly.\cite{PhysRevLett.111.227801} The nitrile group in 5CB |
107 |
+ |
is a particularly interesting case for studying electric field |
108 |
+ |
effects, as 5CB exhibits an isotropic to nematic phase transition that |
109 |
+ |
can be triggered by the application of an external field near room |
110 |
+ |
temperature.\cite{Gray:1973ca,Hatta:1991ee} This presents the |
111 |
+ |
possiblity that the field-induced changes in the local environment |
112 |
+ |
could have dramatic effects on the vibrations of this particular CN |
113 |
+ |
bond. Although the infrared spectroscopy of 5CB has been |
114 |
+ |
well-investigated, particularly as a measure of the kinetics of the |
115 |
+ |
phase transition,\cite{Leyte:1997zl} the 5CB nitrile group has not yet |
116 |
+ |
seen the detailed theoretical treatment that biologically-relevant |
117 |
+ |
small molecules have |
118 |
+ |
received.\cite{Lindquist:2008bh,Lindquist:2008qf,Oh:2008fk,Choi:2008cr,Waegele:2010ve} |
119 |
+ |
|
120 |
|
The fundamental characteristic of liquid crystal mesophases is that |
121 |
|
they maintain some degree of orientational order while translational |
122 |
|
order is limited or absent. This orientational order produces a |
123 |
|
complex direction-dependent response to external perturbations like |
124 |
< |
electric fields and mechanical distortions. The anisotropy of the |
124 |
> |
electric fields and mechanical distortions. The anisotropy of the |
125 |
|
macroscopic phases originates in the anisotropy of the constituent |
126 |
|
molecules, which typically have highly non-spherical structures with a |
127 |
< |
significant degree of internal rigidity. In nematic phases, rod-like |
127 |
> |
significant degree of internal rigidity. In nematic phases, rod-like |
128 |
|
molecules are orientationally ordered with isotropic distributions of |
129 |
< |
molecular centers of mass, while in smectic phases, the molecules |
130 |
< |
arrange themselves into layers with their long (symmetry) axis normal |
131 |
< |
($S_{A}$) or tilted ($S_{C}$) with respect to the layer planes. |
129 |
> |
molecular centers of mass. For example, 5CB has a solid to nematic |
130 |
> |
phase transition at 18C and a nematic to isotropic transition at |
131 |
> |
35C.\cite{Gray:1973ca} |
132 |
|
|
133 |
< |
The behavior of the $S_{A}$ phase can be partially explained with |
134 |
< |
models mainly based on geometric factors and van der Waals |
135 |
< |
interactions. However, these simple models are insufficient to |
136 |
< |
describe liquid crystal phases which exhibit more complex polymorphic |
137 |
< |
nature. X-ray diffraction studies have shown that the ratio between |
138 |
< |
lamellar spacing ($s$) and molecular length ($l$) can take on a wide |
139 |
< |
range of values.\cite{Gray:1984hc,Leadbetter:1976vf,Hardouin:1980yq} |
140 |
< |
Typical $S_{A}$ phases have $s/l$ ratios on the order of $0.8$, while |
141 |
< |
for some compounds, e.g. the 4-alkyl-4'-cyanobiphenyls, the $s/l$ |
142 |
< |
ratio is on the order of $1.4$. Molecules which form $S_{A}$ phases |
143 |
< |
can exhibit a wide variety of subphases like monolayers ($S_{A1}$), |
144 |
< |
uniform bilayers ($S_{A2}$), partial bilayers ($S_{\tilde A}$) as well |
145 |
< |
as interdigitated bilayers ($S_{A_{d}}$), and often have a terminal |
146 |
< |
cyano or nitro group. In particular lyotropic liquid crystals (those |
147 |
< |
exhibiting liquid crystal phase transition as a function of water |
148 |
< |
concentration) often have polar head groups or zwitterionic charge |
149 |
< |
separated groups that result in strong dipolar |
150 |
< |
interactions.\cite{Collings97} Because of their versatile polymorphic |
151 |
< |
nature, polar liquid crystalline materials have important |
152 |
< |
technological applications in addition to their immense relevance to |
124 |
< |
biological systems.\cite{Collings97} |
133 |
> |
In smectic phases, the molecules arrange themselves into layers with |
134 |
> |
their long (symmetry) axis normal ($S_{A}$) or tilted ($S_{C}$) with |
135 |
> |
respect to the layer planes. The behavior of the $S_{A}$ phase can be |
136 |
> |
partially explained with models mainly based on geometric factors and |
137 |
> |
van der Waals interactions. The Gay-Berne potential, in particular, |
138 |
> |
has been widely used in the liquid crystal community to describe this |
139 |
> |
anisotropic phase |
140 |
> |
behavior.~\cite{Gay:1981yu,Berne72,Kushick:1976xy,Luckhurst90,Cleaver:1996rt} |
141 |
> |
However, these simple models are insufficient to describe liquid |
142 |
> |
crystal phases which exhibit more complex polymorphic nature. |
143 |
> |
Molecules which form $S_{A}$ phases can exhibit a wide variety of |
144 |
> |
subphases like monolayers ($S_{A1}$), uniform bilayers ($S_{A2}$), |
145 |
> |
partial bilayers ($S_{\tilde A}$) as well as interdigitated bilayers |
146 |
> |
($S_{A_{d}}$), and often have a terminal cyano or nitro group. In |
147 |
> |
particular, lyotropic liquid crystals (those exhibiting liquid crystal |
148 |
> |
phase transition as a function of water concentration), often have |
149 |
> |
polar head groups or zwitterionic charge separated groups that result |
150 |
> |
in strong dipolar interactions,\cite{Collings:1997rz} and terminal cyano |
151 |
> |
groups (like the one in 5CB) can induce permanent longitudinal |
152 |
> |
dipoles.\cite{Levelut:1981eu} |
153 |
|
|
154 |
< |
Experimental studies by Levelut {\it et al.}~\cite{Levelut:1981eu} |
155 |
< |
revealed that terminal cyano or nitro groups usually induce permanent |
156 |
< |
longitudinal dipole moments on the molecules. |
154 |
> |
Macroscopic electric fields applied using electrodes on opposing sides |
155 |
> |
of a sample of 5CB have demonstrated the phase change of the molecule |
156 |
> |
as a function of electric field.\cite{Lim:2006xq} These previous |
157 |
> |
studies have shown the nitrile group serves as an excellent indicator |
158 |
> |
of the molecular orientation within the applied field. Lee {\it et |
159 |
> |
al.}~showed a 180 degree change in field direction could be probed |
160 |
> |
with the nitrile peak intensity as it changed along with molecular |
161 |
> |
alignment in the field.\cite{Lee:2006qd,Leyte:1997zl} |
162 |
|
|
163 |
< |
Liquid crystalline materials with dipole moments located at the ends |
164 |
< |
of the molecules have important applications in display technologies |
165 |
< |
in addition to their relevance in biological systems.\cite{LCapp} |
163 |
> |
While these macroscopic fields work well at indicating the bulk |
164 |
> |
response, the atomic scale response has not been studied. With the |
165 |
> |
advent of nano-electrodes and coupling them with atomic force |
166 |
> |
microscopy, control of electric fields applied across nanometer |
167 |
> |
distances is now possible.\cite{citation1} While macroscopic fields |
168 |
> |
are insufficient to cause a Stark effect without dielectric breakdown |
169 |
> |
of the material, small fields across nanometer-sized gaps may be of |
170 |
> |
sufficient strength. For a gap of 5 nm between a lower electrode |
171 |
> |
having a nanoelectrode placed near it via an atomic force microscope, |
172 |
> |
a potential of 1 V applied across the electrodes is equivalent to a |
173 |
> |
field of 2x10\textsuperscript{8} $\frac{V}{M}$. This field is |
174 |
> |
certainly strong enough to cause the isotropic-nematic phase change |
175 |
> |
and as well as Stark tuning of the nitrile bond. This should be |
176 |
> |
readily visible experimentally through Raman or IR spectroscopy. |
177 |
|
|
178 |
< |
Many of the technological applications of the lyotropic mesogens |
179 |
< |
manipulate the orientation and structuring of the liquid crystal |
180 |
< |
through application of local electric fields.\cite{?} |
181 |
< |
Macroscopically, this restructuring is visible in the interactions the |
182 |
< |
bulk phase has with scattered or transmitted light.\cite{?} |
178 |
> |
In the sections that follow, we outline a series of coarse-grained |
179 |
> |
classical molecular dynamics simulations of 5CB that were done in the |
180 |
> |
presence of static electric fields. These simulations were then |
181 |
> |
coupled with both {\it ab intio} calculations of CN-deformations and |
182 |
> |
classical bond-length correlation functions to predict spectral |
183 |
> |
shifts. These predictions made should be easily varifiable with |
184 |
> |
scanning electrochemical microscopy experiments. |
185 |
|
|
186 |
< |
4-Cyano-4'-pentylbiphenyl (5CB), has been a model for field-induced |
187 |
< |
phase changes due to the known electric field response of the liquid |
188 |
< |
crystal\cite{Hatta:1991ee}. It was discovered (along with three |
189 |
< |
similar compounds) in 1973 in an effort to find a LC that had a |
190 |
< |
melting point near room temperature.\cite{Gray:1973ca} It's known to |
191 |
< |
have a crystalline to nematic phase transition at 18 C and a nematic |
192 |
< |
to isotropic transition at 35 C.\cite{Gray:1973ca} |
186 |
> |
\section{Computational Details} |
187 |
> |
The force field used for 5CB was taken from Guo {\it et |
188 |
> |
al.}\cite{Zhang:2011hh} However, for most of the simulations, each |
189 |
> |
of the phenyl rings was treated as a rigid body to allow for larger |
190 |
> |
time steps and very long simulation times. The geometries of the |
191 |
> |
rigid bodies were taken from equilibrium bond distances and angles. |
192 |
> |
Although the phenyl rings were held rigid, bonds, bends, torsions and |
193 |
> |
inversion centers that involved atoms in these substructures (but with |
194 |
> |
connectivity to the rest of the molecule) were still included in the |
195 |
> |
potential and force calculations. |
196 |
|
|
197 |
< |
Nitrile groups can serve as very precise electric field reporters via |
198 |
< |
their distinctive Raman and IR signatures.\cite{Boxer:2009xw} The |
199 |
< |
triple bond between the nitrogen and the carbon atom is very sensitive |
200 |
< |
to local field changes and is observed to have a direct impact on the |
201 |
< |
peak position within the spectrum. The Stark shift in the spectrum |
202 |
< |
can be quantified and mapped into a field value that is impinging upon |
203 |
< |
the nitrile bond. This has been used extensively in biological systems |
204 |
< |
like proteins and enzymes.\cite{Tucker:2004qq,Webb:2008kn} |
197 |
> |
Periodic simulations cells containing 270 molecules in random |
198 |
> |
orientations were constructed and were locked at experimental |
199 |
> |
densities. Electrostatic interactions were computed using damped |
200 |
> |
shifted force (DSF) electrostatics.\cite{Fennell:2006zl} The molecules |
201 |
> |
were equilibrated for 1~ns at a temperature of 300K. Simulations with |
202 |
> |
applied fields were carried out in the microcanonical (NVE) ensemble |
203 |
> |
with an energy corresponding to the average energy from the canonical |
204 |
> |
(NVT) equilibration runs. Typical applied-field runs were more than |
205 |
> |
60ns in length. |
206 |
|
|
207 |
< |
To date, the nitrile electric field response of |
208 |
< |
4-Cyano-4'-n-alkylbiphenyl liquid crystals has not been investigated. |
209 |
< |
While macroscopic electric fields applied across macro electrodes show |
210 |
< |
the phase change of the molecule as a function of electric |
211 |
< |
field,\cite{Lim:2006xq} the effect of a nanoscopic field application |
212 |
< |
has not been probed. These previous studies have shown the nitrile |
213 |
< |
group serves as an excellent indicator of the molecular orientation |
214 |
< |
within the field applied. Lee et. al. showed the 180 degree change in field |
165 |
< |
direction could be probed with the nitrile peak intensity as it |
166 |
< |
decreased and increased with molecule alignment in the |
167 |
< |
field.\cite{Lee:2006qd,Leyte:97} |
207 |
> |
Static electric fields with magnitudes similar to what would be |
208 |
> |
available in an experimental setup were applied to the different |
209 |
> |
simulations. With an assumed electrode seperation of 5 nm and an |
210 |
> |
electrostatic potential that is limited by the voltage required to |
211 |
> |
split water (1.23V), the maximum realistic field that could be applied |
212 |
> |
is $\sim 0.024$ V/\AA. Three field environments were investigated: |
213 |
> |
(1) no field applied, (2) partial field = 0.01 V/\AA\ , and (3) full |
214 |
> |
field = 0.024 V/\AA\ . |
215 |
|
|
216 |
< |
While these macroscopic fields worked well at showing the bulk |
217 |
< |
response, the atomic scale response has not been studied. With the |
218 |
< |
advent of nano-electrodes and coupling them with atomic force |
219 |
< |
microscopy, control of electric fields applied across nanometer |
220 |
< |
distances is now possible\cite{citation1}. This application of |
221 |
< |
nanometer length is interesting in the case of a nitrile group on the |
222 |
< |
molecule. While macroscopic fields are insufficient to cause a Stark |
223 |
< |
effect, small fields across nanometer-sized gaps are of sufficient |
224 |
< |
strength. If one were to assume a gap of 5 nm between a lower |
225 |
< |
electrode having a nanoelectrode placed near it via an atomic force |
226 |
< |
microscope, a field of 1 V applied across the electrodes would |
227 |
< |
translate into a field of 2x10\textsuperscript{8} $\frac{V}{M}$. This |
228 |
< |
field is theoretically strong enough to cause a phase change from |
229 |
< |
isotropic to nematic, as well as Stark tuning of the nitrile |
230 |
< |
bond. This should be readily visible experimentally through Raman or |
231 |
< |
IR spectroscopy. |
216 |
> |
After the systems had come to equilibrium under the applied fields, |
217 |
> |
additional simulations were carried out with a flexible (Morse) |
218 |
> |
nitrile bond, |
219 |
> |
\begin{displaymath} |
220 |
> |
V(r_\ce{CN}) = D_e \left(1 - e^{-\beta (r_\ce{CN}-r_e)}\right)^2 |
221 |
> |
\end{displaymath} |
222 |
> |
where $r_e= 1.157437$ \AA , $D_e = 212.95 \mathrm{~kca~l} / |
223 |
> |
\mathrm{mol}^{-1}$ and $\beta = 2.67566 $\AA~$^{-1}$. These |
224 |
> |
parameters correspond to a vibrational frequency of $2375 |
225 |
> |
\mathrm{~cm}^{-1}$, a bit higher than the experimental frequency. The |
226 |
> |
flexible nitrile moiety required simulation time steps of 1~fs, so the |
227 |
> |
additional flexibility was introducuced only after the rigid systems |
228 |
> |
had come to equilibrium under the applied fields. Whenever time |
229 |
> |
correlation functions were computed from the flexible simulations, |
230 |
> |
statistically-independent configurations were sampled from the last ns |
231 |
> |
of the induced-field runs. These configurations were then |
232 |
> |
equilibrated with the flexible nitrile moiety for 100 ps, and time |
233 |
> |
correlation functions were computed using data sampled from an |
234 |
> |
additional 200 ps of run time carried out in the microcanonical |
235 |
> |
ensemble. |
236 |
|
|
237 |
< |
Herein, we show the computational investigation of these electric field effects through atomistic simulations. These are then coupled with ab intio and classical spectrum calculations to predict changes experiments should be able to replicate. |
237 |
> |
\section{Field-induced Nematic Ordering} |
238 |
|
|
239 |
< |
\section{Computational Details} |
240 |
< |
The force field was mainly taken from Guo et al.\cite{Zhang:2011hh} A |
241 |
< |
deviation from this force field was made to create a rigid body from |
242 |
< |
the phenyl rings. Bond distances within the rigid body were taken from |
243 |
< |
equilibrium bond distances. While the phenyl rings were held rigid, |
244 |
< |
bonds, bends, torsions and inversion centers still included the rings. |
239 |
> |
In order to characterize the orientational ordering of the system, the |
240 |
> |
primary quantity of interest is the nematic (orientational) order |
241 |
> |
parameter. This was determined using the tensor |
242 |
> |
\begin{equation} |
243 |
> |
Q_{\alpha \beta} = \frac{1}{2N} \sum_{i=1}^{N} \left(3 \hat{e}_{i |
244 |
> |
\alpha} \hat{e}_{i \beta} - \delta_{\alpha \beta} \right) |
245 |
> |
\end{equation} |
246 |
> |
where $\alpha, \beta = x, y, z$, and $\hat{e}_i$ is the molecular |
247 |
> |
end-to-end unit vector for molecule $i$. The nematic order parameter |
248 |
> |
$S$ is the largest eigenvalue of $Q_{\alpha \beta}$, and the |
249 |
> |
corresponding eigenvector defines the director axis for the phase. |
250 |
> |
$S$ takes on values close to 1 in highly ordered (smectic A) phases, |
251 |
> |
but falls to zero for isotropic fluids. Note that the nitrogen and |
252 |
> |
the terminal chain atom were used to define the vectors for each |
253 |
> |
molecule, so the typical order parameters are lower than if one |
254 |
> |
defined a vector using only the rigid core of the molecule. In |
255 |
> |
nematic phases, typical values for $S$ are close to 0.5. |
256 |
|
|
257 |
< |
Simulations were with boxes of 270 molecules locked at experimental |
258 |
< |
densities with periodic boundaries. The molecules were thermalized |
259 |
< |
from 0 kelvin to 300 kelvin. To equilibrate, each was first run in NVT |
260 |
< |
for 1 ns. This was followed by NVE for simulations used in the data |
261 |
< |
collection. |
257 |
> |
The field-induced phase transition can be clearly seen over the course |
258 |
> |
of a 60 ns equilibration runs in figure \ref{fig:orderParameter}. All |
259 |
> |
three of the systems started in a random (isotropic) packing, with |
260 |
> |
order parameters near 0.2. Over the course 10 ns, the full field |
261 |
> |
causes an alignment of the molecules (due primarily to the interaction |
262 |
> |
of the nitrile group dipole with the electric field). Once this |
263 |
> |
system started exhibiting nematic ordering, the orientational order |
264 |
> |
parameter became stable for the remaining 50 ns of simulation time. |
265 |
> |
It is possible that the partial-field simulation is meta-stable and |
266 |
> |
given enough time, it would eventually find a nematic-ordered phase, |
267 |
> |
but the partial-field simulation was stable as an isotropic phase for |
268 |
> |
the full duration of a 60 ns simulation. Ellipsoidal renderings of the |
269 |
> |
final configurations of the runs shows that the full-field (0.024 |
270 |
> |
V/\AA\ ) experienced a isotropic-nematic phase transition and has |
271 |
> |
ordered with a director axis that is parallel to the direction of the |
272 |
> |
applied field. |
273 |
|
|
274 |
< |
External electric fields were applied in a simplistic charge-field |
275 |
< |
interaction. Forces were calculated by multiplying the electric field |
276 |
< |
being applied by the charge at each atom. For the potential, the |
277 |
< |
origin of the box was used as a point of reference. This allows for a |
278 |
< |
potential value to be added to each atom as the molecule move in space |
279 |
< |
within the box. Fields values were applied in a manner representing |
280 |
< |
those that would be applable in an experimental set-up. The assumed |
281 |
< |
electrode seperation was 5 nm and the field was input as |
282 |
< |
$\frac{V}{\text{\AA}}$. The three field environments were, 1) no field |
283 |
< |
applied, 2) 0.01 $\frac{V}{\text{\AA}}$ (0.5 V) and 3) 0.024 |
284 |
< |
$\frac{V}{\text{\AA}}$ (1.2 V). Each field was applied in the |
212 |
< |
Z-axis of the simulation box. For the simplicity of this paper, |
213 |
< |
each field will be called zero, partial and full, respectively. |
274 |
> |
\begin{figure}[H] |
275 |
> |
\includegraphics[width=\linewidth]{Figure1} |
276 |
> |
\caption{Evolution of the orientational order parameters for the |
277 |
> |
no-field, partial field, and full field simulations over the |
278 |
> |
course of 60 ns. Each simulation was started from a |
279 |
> |
statistically-independent isotropic configuration. On the right |
280 |
> |
are ellipsoids representing the final configurations at three |
281 |
> |
different field strengths: zero field (bottom), partial field |
282 |
> |
(middle), and full field (top)} |
283 |
> |
\label{fig:orderParameter} |
284 |
> |
\end{figure} |
285 |
|
|
215 |
– |
For quantum calculation of the nitrile bond frequency, Gaussian 09 was |
216 |
– |
used. A single 5CB molecule was selected for the center of the |
217 |
– |
cluster. For effects from molecules located near the chosen nitrile |
218 |
– |
group, parts of molecules nearest to the nitrile group were |
219 |
– |
included. For the body not including the tail, any atom within 6~\AA |
220 |
– |
of the midpoint of the nitrile group was included. For the tail |
221 |
– |
structure, the whole tail was included if a tail atom was within 4~\AA |
222 |
– |
of the midpoint. If the tail did not include any atoms from the ring |
223 |
– |
structure, it was considered a propane molecule and included as |
224 |
– |
such. Once the clusters were generated, input files were created that |
225 |
– |
stretched the nitrile bond along its axis from 0.87 to 1.52~\AA at |
226 |
– |
increments of 0.05~\AA. This generated 13 single point energies to be |
227 |
– |
calculated. The Gaussian files were run with B3LYP/6-311++G(d,p) with |
228 |
– |
no other keywords for the zero field simulation. Simulations with |
229 |
– |
fields applied included the keyword ''Field=Z+5'' to match the |
230 |
– |
external field applied in molecular dynamic runs. Once completed, the central nitrile bond frequency |
231 |
– |
was calculated with a Morse fit. Using this fit and the solved energy |
232 |
– |
levels for a Morse oscillator, the frequency was found. Each set of |
233 |
– |
frequencies were then convolved together with a lorezian lineshape |
234 |
– |
function over each value. The width value used was 1.5. For the zero |
235 |
– |
field spectrum, 67 frequencies were used and for the full field, 59 |
236 |
– |
frequencies were used. |
286 |
|
|
287 |
+ |
\section{Sampling the CN bond frequency} |
288 |
+ |
|
289 |
+ |
The vibrational frequency of the nitrile bond in 5CB depends on |
290 |
+ |
features of the local solvent environment of the individual molecules |
291 |
+ |
as well as the bond's orientation relative to the applied field. The |
292 |
+ |
primary quantity of interest for interpreting the condensed phase |
293 |
+ |
spectrum of this vibration is the distribution of frequencies |
294 |
+ |
exhibited by the 5CB nitrile bond under the different electric fields. |
295 |
+ |
Three distinct methods for mapping classical simulations onto |
296 |
+ |
vibrational spectra were brought to bear on these simulations: |
297 |
+ |
\begin{enumerate} |
298 |
+ |
\item Isolated 5CB molecules and their immediate surroundings were |
299 |
+ |
extracted from the simulations. These nitrile bonds were stretched |
300 |
+ |
and single-point {\em ab initio} calculations were used to obtain |
301 |
+ |
Morse-oscillator fits for the local vibrational motion along that |
302 |
+ |
bond. |
303 |
+ |
\item The potential - frequency maps developed by Cho {\it et |
304 |
+ |
al.}~\cite{Oh:2008fk} for nitrile moieties in water were |
305 |
+ |
investigated. This method involves mapping the electrostatic |
306 |
+ |
potential around the bond to the vibrational frequency, and is |
307 |
+ |
similar in approach to field-frequency maps that were pioneered by |
308 |
+ |
Skinner {\it et al.}\cite{XXXX} |
309 |
+ |
\item Classical bond-length autocorrelation functions were Fourier |
310 |
+ |
transformed to directly obtain the vibrational spectrum from |
311 |
+ |
molecular dynamics simulations. |
312 |
+ |
\end{enumerate} |
313 |
+ |
|
314 |
+ |
\subsection{CN frequencies from isolated clusters} |
315 |
+ |
The size of the periodic condensed phase system prevented direct |
316 |
+ |
computation of the complete library of nitrile bond frequencies using |
317 |
+ |
{\it ab initio} methods. In order to sample the nitrile frequencies |
318 |
+ |
present in the condensed-phase, individual molecules were selected |
319 |
+ |
randomly to serve as the center of a local (gas phase) cluster. To |
320 |
+ |
include steric, electrostatic, and other effects from molecules |
321 |
+ |
located near the targeted nitrile group, portions of other molecules |
322 |
+ |
nearest to the nitrile group were included in the quantum mechanical |
323 |
+ |
calculations. The surrounding solvent molecules were divided into |
324 |
+ |
``body'' (the two phenyl rings and the nitrile bond) and ``tail'' (the |
325 |
+ |
alkyl chain). Any molecule which had a body atom within 6~\AA of the |
326 |
+ |
midpoint of the target nitrile bond had its own molecular body (the |
327 |
+ |
4-cyano-biphenyl moiety) included in the configuration. For the alkyl |
328 |
+ |
tail, the entire tail was included if any tail atom was within 4~\AA |
329 |
+ |
of the target nitrile bond. If tail atoms (but no body atoms) were |
330 |
+ |
included within these distances, only the tail was included as a |
331 |
+ |
capped propane molecule. |
332 |
+ |
|
333 |
+ |
\begin{figure}[H] |
334 |
+ |
\includegraphics[width=\linewidth]{Figure2} |
335 |
+ |
\caption{Cluster calculations were performed on randomly sampled 5CB |
336 |
+ |
molecules (shown in red) from each of the simulations. Surrounding |
337 |
+ |
molecular bodies were included if any body atoms were within 6 |
338 |
+ |
\AA\ of the target nitrile bond, and tails were included if they |
339 |
+ |
were within 4 \AA. Included portions of these molecules are shown |
340 |
+ |
in green. The CN bond on the target molecule was stretched and |
341 |
+ |
compressed, and the resulting single point energies were fit to |
342 |
+ |
Morse oscillators to obtain frequency distributions.} |
343 |
+ |
\label{fig:cluster} |
344 |
+ |
\end{figure} |
345 |
+ |
|
346 |
+ |
Inferred hydrogen atom locations were added to the cluster geometries, |
347 |
+ |
and the nitrile bond was stretched from 0.87 to 1.52~\AA\ at |
348 |
+ |
increments of 0.05~\AA. This generated 13 configurations per gas phase |
349 |
+ |
cluster. Single-point energies were computed using the B3LYP |
350 |
+ |
functional~\cite{Becke:1993kq,Lee:1988qf} and the 6-311++G(d,p) basis |
351 |
+ |
set. For the cluster configurations that had been generated from |
352 |
+ |
molecular dynamics running under applied fields, the density |
353 |
+ |
functional calculations had a field of $5 \times 10^{-4}$ atomic units |
354 |
+ |
($E_h / (e a_0)$) applied in the $+z$ direction in order to match the |
355 |
+ |
molecular dynamics simulations. |
356 |
+ |
|
357 |
+ |
The energies for the stretched / compressed nitrile bond in each of |
358 |
+ |
the clusters were used to fit Morse oscillators, and the frequencies |
359 |
+ |
were obtained from the $0 \rightarrow 1$ transition for the energy |
360 |
+ |
levels for this potential.\cite{Morse:1929xy} To obtain a spectrum, |
361 |
+ |
each of the frequencies was convoluted with a Lorentzian lineshape |
362 |
+ |
with a width of 1.5 $\mathrm{cm}^{-1}$. Available computing resources |
363 |
+ |
limited the sampling to 67 clusters for the zero-field spectrum, and |
364 |
+ |
59 for the full field. Comparisons of the quantum mechanical spectrum |
365 |
+ |
to the classical are shown in figure \ref{fig:spectrum}. |
366 |
+ |
|
367 |
+ |
\subsection{CN frequencies from potential-frequency maps} |
368 |
+ |
One approach which has been used to successfully analyze the spectrum |
369 |
+ |
of nitrile and thiocyanate probes in aqueous environments was |
370 |
+ |
developed by Choi {\it et al.}~\cite{Choi:2008cr,Oh:2008fk} This |
371 |
+ |
method involves finding a multi-parameter fit that maps between the |
372 |
+ |
local electrostatic potential at selected sites surrounding the |
373 |
+ |
nitrile bond and the vibrational frequency of that bond obtained from |
374 |
+ |
more expensive {\it ab initio} methods. This approach is similar in |
375 |
+ |
character to the field-frequency maps developed by Skinner {\it et |
376 |
+ |
al.} for OH stretches in liquid water.\cite{XXXX} |
377 |
+ |
|
378 |
+ |
To use the potential-frequency maps, the local electrostatic |
379 |
+ |
potential, $\phi$, is computed at 20 sites ($a = 1 \rightarrow 20$) |
380 |
+ |
that surround the nitrile bond, |
381 |
+ |
\begin{equation} |
382 |
+ |
\phi_{a} = \frac{1}{4\pi \epsilon_{0}} \sum_{j} |
383 |
+ |
\frac{q_j}{\left|r_{aj}\right|}. |
384 |
+ |
\end{equation} |
385 |
+ |
Here $q_j$ is the partial site on atom $j$, and $r_{aj}$ is the |
386 |
+ |
distance between site $a$ and atom $j$. The original map was |
387 |
+ |
parameterized in liquid water and comprises a set of parameters, |
388 |
+ |
$l_a$, that predict the shift in nitrile peak frequency, |
389 |
+ |
\begin{equation} |
390 |
+ |
\delta\tilde{\nu} =\sum^{20}_{a=1} l_{a}\phi_{a} |
391 |
+ |
\end{equation} |
392 |
+ |
|
393 |
+ |
The simulations of 5CB were carried in the presence of external |
394 |
+ |
electric fields, out without water present, so it is not clear if they |
395 |
+ |
can be applied to this situation without extensive |
396 |
+ |
reparameterization. We do, however, suggest a small modification that |
397 |
+ |
would help |
398 |
+ |
|
399 |
+ |
, so the equations need to be corrected |
400 |
+ |
for the frequency shift caused by the electric field. We attempted to |
401 |
+ |
make small modifications the original $\phi^{water}_{a}$ to |
402 |
+ |
$\phi^{5CB}_{a}$ we get, |
403 |
+ |
\begin{equation} |
404 |
+ |
\phi^{5CB}_{a} = \frac{1}{4\pi \epsilon_{0}} \left( \vec{E}\bullet |
405 |
+ |
\left(\vec{r}_{a}-\vec{r}_{CN}\right) \right) + \phi^{5CB}_{0} |
406 |
+ |
\end{equation} |
407 |
+ |
where $\vec{E}$ is the electric field at each atom, $\left( \vec{r}_{a} - |
408 |
+ |
\vec{r}_{CN} \right)$ is the vector between the nitrile bond and the |
409 |
+ |
cooridinates described by Cho around the bond and $\phi^{5CB}_{0}$ is |
410 |
+ |
the correction factor for the system of parameters. After these |
411 |
+ |
changes, the correction factor was found for multiple values of an |
412 |
+ |
external field being applied. However, the factor was no linear and |
413 |
+ |
was overly large due to the fitting parameters being so small. |
414 |
+ |
|
415 |
+ |
|
416 |
+ |
\subsection{CN frequencies from bond length autocorrelation functions} |
417 |
+ |
|
418 |
|
Classical nitrile bond frequencies were found by replacing the rigid |
419 |
< |
cyanide bond with a flexible Morse oscillator bond |
420 |
< |
($r_0= 1.157437$ \AA , $D_0 = 212.95$ and |
421 |
< |
$\beta = 2.67566$) . Once replaced, the |
422 |
< |
systems were allowed to re-equilibrate in NVT for 100 ps. After |
423 |
< |
re-equilibration, the system was run in NVE for 20 ps with a snapshot |
424 |
< |
spacing of 1 fs. These snapshot were then used in bond correlation |
425 |
< |
calculation to find the decay structure of the bond in time using the |
426 |
< |
average bond displacement in time, |
419 |
> |
cyanide bond with a flexible Morse oscillator bond ($r_0= 1.157437$ |
420 |
> |
\AA , $D_0 = 212.95$ and $\beta = 2.67566$). Once replaced, the |
421 |
> |
systems were allowed to re-equilibrate in the canonical (NVT) ensemble |
422 |
> |
for 100 ps. After re-equilibration, the system was run in the |
423 |
> |
microcanonical (NVE) ensemble for 20 ps. Configurations sampled every |
424 |
> |
fs were then used to compute bond-length autocorrelation functions to |
425 |
> |
find the decay structure of the bond in time using the average bond |
426 |
> |
displacement in time, |
427 |
|
\begin{equation} |
428 |
|
C(t) = \langle \left(r(t) - r_0 \right) \cdot \left(r(0) - r_0 \right) \rangle |
429 |
|
\end{equation} |
437 |
|
with a rigid bond followed by an equilibrium run with the bond |
438 |
|
switched back to a Morse oscillator and a short production run of 20 ps. |
439 |
|
|
260 |
– |
\section{Results} |
440 |
|
|
262 |
– |
In order to characterize the orientational ordering of the system, the |
263 |
– |
primary quantity of interest is the nematic (orientational) order |
264 |
– |
parameter. This is determined using the tensor |
265 |
– |
\begin{equation} |
266 |
– |
Q_{\alpha \beta} = \frac{1}{2N} \sum_{i=1}^{N} \left(3 \hat{e}_{i |
267 |
– |
\alpha} \hat{e}_{i \beta} - \delta_{\alpha \beta} \right) |
268 |
– |
\end{equation} |
269 |
– |
where $\alpha, \beta = x, y, z$, and $\hat{e}_i$ is the molecular |
270 |
– |
end-to-end unit vector for molecule $i$. The nematic order parameter |
271 |
– |
$S$ is the largest eigenvalue of $Q_{\alpha \beta}$, and the |
272 |
– |
corresponding eigenvector defines the director axis for the phase. |
273 |
– |
$S$ takes on values close to 1 in highly ordered phases, but falls to |
274 |
– |
zero for isotropic fluids. In the context of 5CB, this value would be |
275 |
– |
close to zero for its isotropic phase and raise closer to one as it |
276 |
– |
moved to the nematic and crystalline phases. |
277 |
– |
|
278 |
– |
This value indicates phases changes at temperature boundaries but also |
279 |
– |
when a phase change occurs due to external field applications. In |
280 |
– |
Figure 1, this phase change can be clearly seen over the course of 60 |
281 |
– |
ns. Each system starts with an ordering paramter near 0.1 to 0.2, |
282 |
– |
which is an isotropic phase. Over the course 10 ns, the full external field |
283 |
– |
causes a shift in the ordering of the system to 0.5, the nematic phase |
284 |
– |
of the liquid crystal. This change is consistent over the full 50 ns |
285 |
– |
with no drop back into the isotropic phase. This change is clearly |
286 |
– |
field induced and stable over a long period of simulation time. |
287 |
– |
\begin{figure} |
288 |
– |
\includegraphics[trim = 5mm 10mm 3mm 10mm, clip, width=3.25in]{P2} |
289 |
– |
\caption{Ordering of each external field application over the course |
290 |
– |
of 60 ns with a sampling every 100 ps. Each trajectory was started |
291 |
– |
after equilibration with zero field applied.} |
292 |
– |
\label{fig:orderParameter} |
293 |
– |
\end{figure} |
294 |
– |
|
295 |
– |
Interestingly, the field that is needed to switch the phase of 5CB |
296 |
– |
macroscopically is larger than 1 V. However, in this case, only a |
297 |
– |
voltage of 1.2 V was need to induce a phase change. This is impart due |
298 |
– |
to the distance the field is being applied across. At such a small |
299 |
– |
distance, the field is much larger than the macroscopic and thus |
300 |
– |
easily induces a field dependent phase change. |
301 |
– |
|
302 |
– |
In the figure below, this phase change is represented nicely as |
303 |
– |
ellipsoids that are created by the vector formed between the nitrogen |
304 |
– |
of the nitrile group and the tail terminal carbon atom. These |
305 |
– |
illistrate the change from isotropic phase to nematic change. Both the |
306 |
– |
zero field and partial field images look mostly disordered. The |
307 |
– |
partial field does look somewhat orded but without any overall order |
308 |
– |
of the entire system. This is most likely a high point in the ordering |
309 |
– |
for the trajectory. The full field image on the other hand looks much |
310 |
– |
more ordered when compared to the two lower field simulations. |
311 |
– |
\begin{figure} |
312 |
– |
\includegraphics[width=7in]{Elip_3} |
313 |
– |
\caption{Ellipsoid reprsentation of 5CB at three different |
314 |
– |
field strengths, Zero Field (Left), Partial Field (Middle), and Full |
315 |
– |
Field (Right) Each image was created by taking the final |
316 |
– |
snapshot of each 60 ns run} |
317 |
– |
\label{fig:Cigars} |
318 |
– |
\end{figure} |
319 |
– |
|
441 |
|
This change in phase was followed by two courses of further |
442 |
|
analysis. First was the replacement of the static nitrile bond with a |
443 |
|
morse oscillator bond. This was then simulated for a period of time |
465 |
|
\label{fig:twoSpectra} |
466 |
|
\end{figure} |
467 |
|
|
347 |
– |
Before Gaussian silumations were carried out, it was attempt to apply |
348 |
– |
the method developed by Cho et. al. This method involves the fitting |
349 |
– |
of multiple parameters to Gaussian calculated freuencies to find a |
350 |
– |
correlation between the potential around the bond and the |
351 |
– |
frequency. This is very similar to work done by Skinner et. al. with |
352 |
– |
water models like SPC/E. The general method is to find the shift in |
353 |
– |
the peak position through, |
354 |
– |
\begin{equation} |
355 |
– |
\delta\tilde{\nu} =\sum^{n}_{a=1} l_{a}\phi^{water}_{a} |
356 |
– |
\end{equation} |
357 |
– |
where $l_{a}$ are the fitting parameters and $\phi^{water}_{a}$ is the |
358 |
– |
potential from the surrounding water cluster. This $\phi^{water}_{a}$ |
359 |
– |
takes the form, |
360 |
– |
\begin{equation} |
361 |
– |
\phi_{a} = \frac{1}{4\pi \epsilon_{0}} \sum_{m} \sum_{j} |
362 |
– |
\frac{C^{H_{2}O}_{j \left(m \right) }}{r_{aj \left(m\right)}} |
363 |
– |
\end{equation} |
364 |
– |
where $C^{H_{2}O}_{j \left(m \right) }$ indicates the partial charge |
365 |
– |
on the $j$th site of the $m$th water molecule and $r_{aj \left(m\right)}$ |
366 |
– |
is the distance between the site $a$ of the nitrile molecule and the $j$th |
367 |
– |
site of the $m$th water molecule. However, since these simulations |
368 |
– |
are done under the presence of external electric fields and in the |
369 |
– |
absence of water the equations must have a correction factor for the |
370 |
– |
external field change as well as the use of electric field site data |
371 |
– |
instead of charged site points. So by modifing the original |
372 |
– |
$\phi^{water}_{a}$ to $\phi^{5CB}_{a}$ we get, |
373 |
– |
\begin{equation} |
374 |
– |
\phi^{5CB}_{a} = \frac{1}{4\pi \epsilon_{0}} \left( \vec{E}\bullet |
375 |
– |
\left(\vec{r}_{a}-\vec{r}_{CN}\right) \right) + \phi^{5CB}_{0} |
376 |
– |
\end{equation} |
377 |
– |
where $\vec{E}$ is the electric field at each atom, $\left( \vec{r}_{a} - |
378 |
– |
\vec{r}_{CN} \right)$ is the vector between the nitrile bond and the |
379 |
– |
cooridinates described by Cho around the bond and $\phi^{5CB}_{0}$ is |
380 |
– |
the correction factor for the system of parameters. After these |
381 |
– |
changes, the correction factor was found for multiple values of an |
382 |
– |
external field being applied. However, the factor was no linear and |
383 |
– |
was overly large due to the fitting parameters being so small. |
468 |
|
|
469 |
|
Due to this, Gaussian calculations were performed in lieu of this |
470 |
|
method. A set of snapshots for the zero and full field simualtions, |
480 |
|
Since no explicit dependence was observed between the calculated |
481 |
|
frequency and the electric field, it was not a viable route for the |
482 |
|
calculation of a nitrile spectrum. Instead, the frequencies were taken |
483 |
< |
and convolved together. These two spectra are seen below in Figure |
483 |
> |
and convolved together with a lorentzian line shape applied around the |
484 |
> |
frequency value. These spectra are seen below in Figure |
485 |
|
4. While the spectrum without a field is lower in intensity and is |
486 |
< |
almost bimodel in distrobuiton, the external field spectrum is much |
486 |
> |
almost bimodel in distrobution, the external field spectrum is much |
487 |
|
more unimodel. This tighter clustering has the affect of increasing the |
488 |
|
intensity around 2226 cm\textsuperscript{-1} where the peak is |
489 |
|
centered. The external field also has fewer frequencies of higher |
496 |
|
\label{fig:Con} |
497 |
|
\end{figure} |
498 |
|
\section{Discussion} |
499 |
+ |
Interestingly, the field that is needed to switch the phase of 5CB |
500 |
+ |
macroscopically is larger than 1 V. However, in this case, only a |
501 |
+ |
voltage of 1.2 V was need to induce a phase change. This is impart due |
502 |
+ |
to the short distance of 5 nm the field is being applied across. At such a small |
503 |
+ |
distance, the field is much larger than the macroscopic and thus |
504 |
+ |
easily induces a field dependent phase change. However, this field |
505 |
+ |
will not cause a breakdown of the 5CB since electrochemistry studies |
506 |
+ |
have shown that it can be used in the presence of fields as high as |
507 |
+ |
500 V macroscopically. This large of a field near the surface of the |
508 |
+ |
elctrode would cause breakdown of 5CB if it could happen. |
509 |
+ |
|
510 |
|
The absence of any electric field dependency of the freuquency with |
511 |
< |
the Gaussian simulations is strange. A large base of research has been |
512 |
< |
built upon the known tuning of the nitrile bond as local field |
513 |
< |
changes. This differences may be due to the absence of water. Many of |
514 |
< |
the nitrile bond fitting maps are done in the presence of liquid |
419 |
< |
water. Liquid water is known to have a very high internal field which |
511 |
> |
the Gaussian simulations is interesting. A large base of research has been |
512 |
> |
built upon the known tuning of the nitrile bond as the local field |
513 |
> |
changes. This difference may be due to the absence of water or a |
514 |
> |
molecule that induces a large internal field. Liquid water is known to have a very high internal field which |
515 |
|
is much larger than the internal fields of neat 5CB. Even though the |
516 |
< |
application of Gaussian simulations followed by mappying to |
516 |
> |
application of Gaussian simulations followed by mapping it to |
517 |
|
some classical parameter is easy and straight forward, this system |
518 |
|
illistrates how that 'go to' method can break down. |
519 |
|
|
520 |
|
While this makes the application of nitrile Stark effects in |
521 |
< |
simulations of water absent simulations harder, these data show |
521 |
> |
simulations without water harder, these data show |
522 |
|
that it is not a deal breaker. The classically calculated nitrile |
523 |
|
spectrum shows changes in the spectra that will be easily seen through |
524 |
|
experimental routes. It indicates a shifted peak lower in energy |
525 |
< |
should arise. This peak is a few wavenumbers from the larger peak and |
526 |
< |
almost 75 wavenumbers from the center. This seperation between the two |
527 |
< |
peaks means experimental results will have an easily resolved peak. |
525 |
> |
should arise. This peak is a few wavenumbers from the leading edge of |
526 |
> |
the larger peak and almost 75 wavenumbers from the center. This |
527 |
> |
seperation between the two peaks means experimental results will show |
528 |
> |
an easily resolved peak. |
529 |
|
|
530 |
< |
The Gaussian derived spectra do indicate that with an applied field |
530 |
> |
The Gaussian derived spectra do indicate an applied field |
531 |
|
and subsiquent phase change does cause a narrowing of freuency |
532 |
< |
distrobution. |
532 |
> |
distrobution. With narrowing, it would indicate an increased |
533 |
> |
homogeneous distrobution of the local field near the nitrile. |
534 |
|
\section{Conclusions} |
535 |
< |
Field dependent changes in the phase of a system are |
439 |
< |
Jonathan K. Whitmer |
440 |
< |
cho stuff |
535 |
> |
Field dependent changes |
536 |
|
\newpage |
537 |
|
|
538 |
|
\bibliography{5CB} |