1 |
\documentclass[journal = jpccck, manuscript = article]{achemso} |
2 |
\setkeys{acs}{usetitle = true} |
3 |
|
4 |
\usepackage{caption} |
5 |
\usepackage{float} |
6 |
\usepackage{geometry} |
7 |
\usepackage{natbib} |
8 |
\usepackage{setspace} |
9 |
\usepackage{xkeyval} |
10 |
\usepackage{amsmath} |
11 |
\usepackage{amssymb} |
12 |
\usepackage{times} |
13 |
\usepackage{mathptm} |
14 |
\usepackage{setspace} |
15 |
%\usepackage{endfloat} |
16 |
\usepackage{tabularx} |
17 |
\usepackage{longtable} |
18 |
\usepackage{graphicx} |
19 |
\usepackage{multirow} |
20 |
\usepackage{multicol} |
21 |
\usepackage{achemso} |
22 |
\usepackage{subcaption} |
23 |
\usepackage[colorinlistoftodos]{todonotes} |
24 |
\usepackage[version=3]{mhchem} % this is a great package for formatting chemical reactions |
25 |
% \usepackage[square, comma, sort&compress]{natbib} |
26 |
\usepackage{url} |
27 |
\pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm |
28 |
\evensidemargin 0.0cm \topmargin -21pt \headsep 10pt \textheight |
29 |
9.0in \textwidth 6.5in \brokenpenalty=10000 |
30 |
|
31 |
% double space list of tables and figures |
32 |
%\AtBeginDelayedFloats{\renewcomand{\baselinestretch}{1.66}} |
33 |
\setlength{\abovecaptionskip}{20 pt} |
34 |
\setlength{\belowcaptionskip}{30 pt} |
35 |
|
36 |
% \bibpunct{}{}{,}{s}{}{;} |
37 |
|
38 |
%\citestyle{nature} |
39 |
% \bibliographystyle{achemso} |
40 |
|
41 |
|
42 |
\title{Nitrile vibrations as reporters of field-induced phase |
43 |
transitions in liquid crystals} |
44 |
\author{James M. Marr} |
45 |
\author{J. Daniel Gezelter} |
46 |
\email{gezelter@nd.edu} |
47 |
\affiliation[University of Notre Dame]{251 Nieuwland Science Hall\\ |
48 |
Department of Chemistry and Biochemistry\\ |
49 |
University of Notre Dame\\ |
50 |
Notre Dame, Indiana 46556} |
51 |
|
52 |
\date{\today} |
53 |
|
54 |
\begin{document} |
55 |
|
56 |
\maketitle |
57 |
|
58 |
\begin{doublespace} |
59 |
|
60 |
\begin{abstract} |
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 various contributions to the Stark shift of the terminal |
66 |
nitrile group. A field-induced isotropic-nematic phase transition |
67 |
was observed in the simulations, and the effects of this transition |
68 |
on the distribution of nitrile frequencies were computed. Classical |
69 |
bond displacement correlation functions exhibited a ($\sim 40 |
70 |
\mathrm{cm}^{-1}$ red shift of a fraction of the main nitrile peak, |
71 |
and this shift was observed only when the fields were large enough |
72 |
to induce orientational ordering of the bulk phase. Our simulations |
73 |
appear to indicate that phase-induced changes to the local |
74 |
surroundings are a larger contribution to the change in the nitrile |
75 |
spectrum than the direct field contribution. |
76 |
\end{abstract} |
77 |
|
78 |
\newpage |
79 |
|
80 |
\section{Introduction} |
81 |
The Stark shift of nitrile groups in response to applied electric |
82 |
fields have been used extensively in biology for probing the internal |
83 |
fields of structures like proteins and DNA. Integration of these |
84 |
probes into different materials is also important for studying local |
85 |
structure in confined fluids. This work centers on the vibrational |
86 |
response of the terminal nitrile group in 4-Cyano-4'-pentylbiphenyl |
87 |
(5CB), a liquid crystalline molecule with an isotropic to nematic |
88 |
phase transition that can be triggered by the application of an |
89 |
external field. |
90 |
|
91 |
The fundamental characteristic of liquid crystal mesophases is that |
92 |
they maintain some degree of orientational order while translational |
93 |
order is limited or absent. This orientational order produces a |
94 |
complex direction-dependent response to external perturbations like |
95 |
electric fields and mechanical distortions. The anisotropy of the |
96 |
macroscopic phases originates in the anisotropy of the constituent |
97 |
molecules, which typically have highly non-spherical structures with a |
98 |
significant degree of internal rigidity. In nematic phases, rod-like |
99 |
molecules are orientationally ordered with isotropic distributions of |
100 |
molecular centers of mass, while in smectic phases, the molecules |
101 |
arrange themselves into layers with their long (symmetry) axis normal |
102 |
($S_{A}$) or tilted ($S_{C}$) with respect to the layer planes. |
103 |
|
104 |
The behavior of the $S_{A}$ phase can be partially explained with |
105 |
models mainly based on geometric factors and van der Waals |
106 |
interactions. However, these simple models are insufficient to |
107 |
describe liquid crystal phases which exhibit more complex polymorphic |
108 |
nature. X-ray diffraction studies have shown that the ratio between |
109 |
lamellar spacing ($s$) and molecular length ($l$) can take on a wide |
110 |
range of values.\cite{Gray:1984hc,Leadbetter:1976vf,Hardouin:1980yq} |
111 |
Typical $S_{A}$ phases have $s/l$ ratios on the order of $0.8$, while |
112 |
for some compounds, e.g. the 4-alkyl-4'-cyanobiphenyls, the $s/l$ |
113 |
ratio is on the order of $1.4$. Molecules which form $S_{A}$ phases |
114 |
can exhibit a wide variety of subphases like monolayers ($S_{A1}$), |
115 |
uniform bilayers ($S_{A2}$), partial bilayers ($S_{\tilde A}$) as well |
116 |
as interdigitated bilayers ($S_{A_{d}}$), and often have a terminal |
117 |
cyano or nitro group. In particular, lyotropic liquid crystals (those |
118 |
exhibiting liquid crystal phase transition as a function of water |
119 |
concentration), often have polar head groups or zwitterionic charge |
120 |
separated groups that result in strong dipolar |
121 |
interactions.\cite{Collings97} Because of their versatile polymorphic |
122 |
nature, polar liquid crystalline materials have important |
123 |
technological applications in addition to their immense relevance to |
124 |
biological systems.\cite{Collings97} |
125 |
|
126 |
Experimental studies by Levelut {\it et al.}~\cite{Levelut:1981eu} |
127 |
revealed that terminal cyano or nitro groups usually induce permanent |
128 |
longitudinal dipole moments on the molecules. Liquid crystalline |
129 |
materials with dipole moments located at the ends of the molecules |
130 |
have important applications in display technologies in addition to |
131 |
their relevance in biological systems.\cite{LCapp} |
132 |
|
133 |
Many of the technological applications of the lyotropic mesogens |
134 |
manipulate the orientation and structuring of the liquid crystal |
135 |
through application of external electric fields.\cite{?} |
136 |
Macroscopically, this restructuring is visible in the interactions the |
137 |
bulk phase has with scattered or transmitted light.\cite{?} |
138 |
|
139 |
4-Cyano-4'-pentylbiphenyl (5CB), has been a model for field-induced |
140 |
phase changes due to the well-studied electric field |
141 |
response,\cite{Hatta:1991ee} and the fact that it has a set of phase |
142 |
transitions near room temperature.\cite{Gray:1973ca} The have a solid |
143 |
to nematic phase transition at 18 C and a nematic to isotropic |
144 |
transition at 35 C.\cite{Gray:1973ca} Recently there has been renewed |
145 |
interest in 5CB in nanodroplets to understand defect sites and |
146 |
nanoparticle structuring.\cite{PhysRevLett.111.227801} |
147 |
|
148 |
Nitrile groups can serve as very precise electric field reporters via |
149 |
their distinctive Raman and IR signatures.\cite{Boxer:2009xw} The |
150 |
triple bond between the nitrogen and the carbon atom is very sensitive |
151 |
to local field changes and is observed to have a direct impact on the |
152 |
peak position within the spectrum. The Stark shift in the spectrum |
153 |
can be quantified and mapped into a field value that is impinging upon |
154 |
the nitrile bond. This has been used extensively in biological systems |
155 |
like proteins and enzymes.\cite{Tucker:2004qq,Webb:2008kn} |
156 |
|
157 |
To date, the nitrile electric field response of |
158 |
4-Cyano-4'-n-alkylbiphenyl liquid crystals has not been investigated. |
159 |
While macroscopic electric fields applied across macro electrodes show |
160 |
the phase change of the molecule as a function of electric |
161 |
field,\cite{Lim:2006xq} the effect of a nanoscopic field application |
162 |
has not been probed. These previous studies have shown the nitrile |
163 |
group serves as an excellent indicator of the molecular orientation |
164 |
within the field applied. Lee {\it et al.}~showed the 180 degree |
165 |
change in field direction could be probed with the nitrile peak |
166 |
intensity as it decreased and increased with molecule alignment in the |
167 |
field.\cite{Lee:2006qd,Leyte:97} |
168 |
|
169 |
While these macroscopic fields worked well at showing the bulk |
170 |
response, the atomic scale response has not been studied. With the |
171 |
advent of nano-electrodes and coupling them with atomic force |
172 |
microscopy, control of electric fields applied across nanometer |
173 |
distances is now possible\cite{citation1}. This application of |
174 |
nanometer length is interesting in the case of a nitrile group on the |
175 |
molecule. While macroscopic fields are insufficient to cause a Stark |
176 |
effect, small fields across nanometer-sized gaps are of sufficient |
177 |
strength. If one were to assume a gap of 5 nm between a lower |
178 |
electrode having a nanoelectrode placed near it via an atomic force |
179 |
microscope, a field of 1 V applied across the electrodes would |
180 |
translate into a field of 2x10\textsuperscript{8} $\frac{V}{M}$. This |
181 |
field is theoretically strong enough to cause a phase change from |
182 |
isotropic to nematic, as well as Stark tuning of the nitrile |
183 |
bond. This should be readily visible experimentally through Raman or |
184 |
IR spectroscopy. |
185 |
|
186 |
In the rest of this paper, we outline a series of classical molecular |
187 |
dynamics simulations of 5CB that were done in the presence of static |
188 |
electric fields. These simulations were then coupled with both {\it ab |
189 |
intio} calculations of CN-deformations and classical correlation |
190 |
functions to predict spectral shifts. These predictions should be |
191 |
easily varifiable with scanning electrochemical microscopy |
192 |
experiments. |
193 |
|
194 |
\section{Computational Details} |
195 |
The force field used for 5CB was taken from Guo {\it et |
196 |
al.}\cite{Zhang:2011hh} However, for most of the simulations, each |
197 |
of the phenyl rings was treated as a rigid body to allow for larger |
198 |
time steps and very long simulation times. The geometries of the |
199 |
rigid bodies were taken from equilibrium bond distances and angles. |
200 |
Although the phenyl rings were held rigid, bonds, bends, torsions and |
201 |
inversion centers included in these bodies (but with connectivity to |
202 |
the rest of the molecule) were still included in the potential and |
203 |
force calculations. |
204 |
|
205 |
Periodic simulations cells contained 270 molecules and were locked at |
206 |
experimental densities. Electrostatic interactions were computed |
207 |
using damped shifted force (DSF) electrostatics.\cite{Fennell:2006zl} |
208 |
The molecules were equilibrated for 1~ns at a temperature of 300K. |
209 |
Simulations with applied fields were carried out in the microcanonical |
210 |
(NVE) ensemble with an energy corresponding to the average energy from |
211 |
the canonical (NVT) equilibration runs. Typical applied-field runs |
212 |
were more than 60ns in length. |
213 |
|
214 |
Static electric fields with magnitudes similar to what would be |
215 |
available in an experimental setup were applied to the different |
216 |
simulations. With an assumed electrode seperation of 5 nm and an |
217 |
electrostatic potential that is limited by the voltage required to |
218 |
split water (1.23V), the maximum realistic field that could be applied |
219 |
is $\sim 0.024 V / \AA$. Three field environments were investigated: |
220 |
(1) no field applied, (2) $0.01 V / \AA$ (0.5 V), and (3) $0.024 V / |
221 |
\AA$ (1.2 V). Each field was applied along the $z$-axis of the |
222 |
simulation cell. For simplicity, these field strengths will be |
223 |
referred to as zero, partial, and full field. |
224 |
|
225 |
After the systems had come to equilibrium under the applied fields, |
226 |
additional simulations were carried out with a flexible (harmonic) |
227 |
nitrile bond with an equilibrium bond distance of XXX \AA and a force |
228 |
constant of XXX kcal / mol $\AA^2$, corresponding to a vibrational |
229 |
frequency of YYYY $\mathrm{cm}^{-1}$. The flexible nitrile moiety |
230 |
required simualtion time steps of 1fs, so the additional flexibility |
231 |
was introducuced only after the rigid systems had come to equilibrium |
232 |
under the applied fields. Whenever time correlation functions were |
233 |
computed from the flexible simulations, statistically-independent |
234 |
configurations were sampled from the last ns of the induced-field |
235 |
runs. These configurations were then equilibrated with the flexible |
236 |
nitrile moiety for 100 ps, and time correlation functions were |
237 |
computed using data sampled from an additional 200 ps of run time |
238 |
carried out in the microcanonical ensemble. |
239 |
|
240 |
\section{Field-induced Nematic Ordering} |
241 |
|
242 |
In order to characterize the orientational ordering of the system, the |
243 |
primary quantity of interest is the nematic (orientational) order |
244 |
parameter. This was determined using the tensor |
245 |
\begin{equation} |
246 |
Q_{\alpha \beta} = \frac{1}{2N} \sum_{i=1}^{N} \left(3 \hat{e}_{i |
247 |
\alpha} \hat{e}_{i \beta} - \delta_{\alpha \beta} \right) |
248 |
\end{equation} |
249 |
where $\alpha, \beta = x, y, z$, and $\hat{e}_i$ is the molecular |
250 |
end-to-end unit vector for molecule $i$. The nematic order parameter |
251 |
$S$ is the largest eigenvalue of $Q_{\alpha \beta}$, and the |
252 |
corresponding eigenvector defines the director axis for the phase. |
253 |
$S$ takes on values close to 1 in highly ordered (smectic A) phases, |
254 |
but falls to zero for isotropic fluids. In nematic phases, typical |
255 |
values are close to 0.5. |
256 |
|
257 |
In Figure \ref{fig:orderParameter}, the field-induced phase change can |
258 |
be clearly seen over the course of a 60 ns equilibration run. 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 landed in the nematic-ordered state, it became stable for the |
264 |
remaining 50 ns of simulation time. It is possible that the |
265 |
partial-field simulation is meta-stable and given enough time, it |
266 |
would eventually find a nematic-ordered phase, but the partial-field |
267 |
simulation was stable as an isotropic phase for the full duration of a |
268 |
60 ns simulation. |
269 |
\begin{figure} |
270 |
\includegraphics[trim = 5mm 10mm 3mm 10mm, clip, width=3.25in]{P2} |
271 |
\caption{Ordering of each external field application over the course |
272 |
of 60 ns with a sampling every 100 ps. Each trajectory was started |
273 |
after equilibration with zero field applied.} |
274 |
\label{fig:orderParameter} |
275 |
\end{figure} |
276 |
|
277 |
In figure \ref{fig:Cigars}, the field-induced isotropic-nematic |
278 |
transition is represented using ellipsoids aligned along the long-axis |
279 |
of each molecule. The vector between the nitrogen of the nitrile |
280 |
group and the terminal tail atom is used to orient each |
281 |
ellipsoid. Both the zero field and partial field simulations appear |
282 |
isotropic, while the full field simulations has clearly been |
283 |
orientationally ordered |
284 |
\begin{figure} |
285 |
\includegraphics[width=7in]{Elip_3} |
286 |
\caption{Ellipsoid reprsentation of 5CB at three different |
287 |
field strengths, Zero Field (Left), Partial Field (Middle), and Full |
288 |
Field (Right) Each image was created by taking the final |
289 |
snapshot of each 60 ns run} |
290 |
\label{fig:Cigars} |
291 |
\end{figure} |
292 |
|
293 |
\section{Analysis} |
294 |
|
295 |
For quantum calculation of the nitrile bond frequency, Gaussian 09 was |
296 |
used. A single 5CB molecule was selected for the center of the |
297 |
cluster. For effects from molecules located near the chosen nitrile |
298 |
group, parts of molecules nearest to the nitrile group were |
299 |
included. For the body not including the tail, any atom within 6~\AA |
300 |
of the midpoint of the nitrile group was included. For the tail |
301 |
structure, the whole tail was included if a tail atom was within 4~\AA |
302 |
of the midpoint. If the tail did not include any atoms from the ring |
303 |
structure, it was considered a propane molecule and included as |
304 |
such. Once the clusters were generated, input files were created that |
305 |
stretched the nitrile bond along its axis from 0.87 to 1.52~\AA at |
306 |
increments of 0.05~\AA. This generated 13 single point energies to be |
307 |
calculated. The Gaussian files were run with B3LYP/6-311++G(d,p) with |
308 |
no other keywords for the zero field simulation. Simulations with |
309 |
fields applied included the keyword ''Field=Z+5'' to match the |
310 |
external field applied in molecular dynamic runs. Once completed, the central nitrile bond frequency |
311 |
was calculated with a Morse fit. Using this fit and the solved energy |
312 |
levels for a Morse oscillator, the frequency was found. Each set of |
313 |
frequencies were then convolved together with a lorezian lineshape |
314 |
function over each value. The width value used was 1.5. For the zero |
315 |
field spectrum, 67 frequencies were used and for the full field, 59 |
316 |
frequencies were used. |
317 |
|
318 |
Classical nitrile bond frequencies were found by replacing the rigid |
319 |
cyanide bond with a flexible Morse oscillator bond |
320 |
($r_0= 1.157437$ \AA , $D_0 = 212.95$ and |
321 |
$\beta = 2.67566$) . Once replaced, the |
322 |
systems were allowed to re-equilibrate in NVT for 100 ps. After |
323 |
re-equilibration, the system was run in NVE for 20 ps with a snapshot |
324 |
spacing of 1 fs. These snapshot were then used in bond correlation |
325 |
calculation to find the decay structure of the bond in time using the |
326 |
average bond displacement in time, |
327 |
\begin{equation} |
328 |
C(t) = \langle \left(r(t) - r_0 \right) \cdot \left(r(0) - r_0 \right) \rangle |
329 |
\end{equation} |
330 |
% |
331 |
where $r_0$ is the equilibrium bond distance and $r(t)$ is the |
332 |
instantaneous bond displacement at time $t$. Once calculated, |
333 |
smoothing was applied by adding an exponential decay on top of the |
334 |
decay with a $\tau$ of 6000. Further smoothing |
335 |
was applied by padding 20,000 zeros on each side of the symmetric |
336 |
data. This was done five times by allowing the systems to run 1 ns |
337 |
with a rigid bond followed by an equilibrium run with the bond |
338 |
switched back to a Morse oscillator and a short production run of 20 ps. |
339 |
|
340 |
\section{Results} |
341 |
|
342 |
|
343 |
|
344 |
|
345 |
This change in phase was followed by two courses of further |
346 |
analysis. First was the replacement of the static nitrile bond with a |
347 |
morse oscillator bond. This was then simulated for a period of time |
348 |
and a classical spetrum was calculated. Second, ab intio calcualtions |
349 |
were performed to investigate if the phase change caused any change |
350 |
spectrum through quantum effects. |
351 |
|
352 |
The classical nitrile spectrum can be seen in Figure 2. Most noticably |
353 |
is the position of the two peaks. Obviously the experimental peak |
354 |
position is near 2226 cm\textsuperscript{-1}. However, in this case |
355 |
the peak position is shifted to the blue at a position of 2375 |
356 |
cm\textsuperscript{-1}. This shift is due solely to the choice of |
357 |
oscillator strength in the Morse oscillator parameters. While this |
358 |
shift makes the two spectra differ, it does not affect the ability to |
359 |
qualitatively compare peak changes to possible experimental changes. |
360 |
With this important fact out of the way, differences between the two |
361 |
states are subtle but are very much present. The first and |
362 |
most notable is the apperance for a strong band near 2300 |
363 |
cm\textsuperscript{-1}. |
364 |
\begin{figure} |
365 |
\includegraphics[width=3.25in]{2Spectra} |
366 |
\caption{The classically calculated nitrile bond spetrum for no |
367 |
external field application (black) and full external field |
368 |
application (red)} |
369 |
\label{fig:twoSpectra} |
370 |
\end{figure} |
371 |
|
372 |
Before Gaussian silumations were carried out, it was attempt to apply |
373 |
the method developed by Cho {\it et al.}~\cite{Oh:2008fk} This method involves the fitting |
374 |
of multiple parameters to Gaussian calculated freuencies to find a |
375 |
correlation between the potential around the bond and the |
376 |
frequency. This is very similar to work done by Skinner {\it et al.}~with |
377 |
water models like SPC/E. The general method is to find the shift in |
378 |
the peak position through, |
379 |
\begin{equation} |
380 |
\delta\tilde{\nu} =\sum^{n}_{a=1} l_{a}\phi^{water}_{a} |
381 |
\end{equation} |
382 |
where $l_{a}$ are the fitting parameters and $\phi^{water}_{a}$ is the |
383 |
potential from the surrounding water cluster. This $\phi^{water}_{a}$ |
384 |
takes the form, |
385 |
\begin{equation} |
386 |
\phi_{a} = \frac{1}{4\pi \epsilon_{0}} \sum_{m} \sum_{j} |
387 |
\frac{C^{H_{2}O}_{j \left(m \right) }}{r_{aj \left(m\right)}} |
388 |
\end{equation} |
389 |
where $C^{H_{2}O}_{j \left(m \right) }$ indicates the partial charge |
390 |
on the $j$th site of the $m$th water molecule and $r_{aj \left(m\right)}$ |
391 |
is the distance between the site $a$ of the nitrile molecule and the $j$th |
392 |
site of the $m$th water molecule. However, since these simulations |
393 |
are done under the presence of external fields and in the |
394 |
absence of water, the equations need a correction factor for the shift |
395 |
caused by the external field. The equation is also reworked to use |
396 |
electric field site data instead of partial charges from surrounding |
397 |
atoms. So by modifing the original |
398 |
$\phi^{water}_{a}$ to $\phi^{5CB}_{a}$ we get, |
399 |
\begin{equation} |
400 |
\phi^{5CB}_{a} = \frac{1}{4\pi \epsilon_{0}} \left( \vec{E}\bullet |
401 |
\left(\vec{r}_{a}-\vec{r}_{CN}\right) \right) + \phi^{5CB}_{0} |
402 |
\end{equation} |
403 |
where $\vec{E}$ is the electric field at each atom, $\left( \vec{r}_{a} - |
404 |
\vec{r}_{CN} \right)$ is the vector between the nitrile bond and the |
405 |
cooridinates described by Cho around the bond and $\phi^{5CB}_{0}$ is |
406 |
the correction factor for the system of parameters. After these |
407 |
changes, the correction factor was found for multiple values of an |
408 |
external field being applied. However, the factor was no linear and |
409 |
was overly large due to the fitting parameters being so small. |
410 |
|
411 |
Due to this, Gaussian calculations were performed in lieu of this |
412 |
method. A set of snapshots for the zero and full field simualtions, |
413 |
they were first investigated for any dependence on the local, with |
414 |
external field included, electric field. This was to see if a linear |
415 |
or non-linear relationship between the two could be utilized for |
416 |
generating spectra. This was done in part because of previous studies |
417 |
showing the frequency dependence of nitrile bonds to the electric |
418 |
fields generated locally between solvating water. It was seen that |
419 |
little to no dependence could be directly shown. This data is not |
420 |
shown. |
421 |
|
422 |
Since no explicit dependence was observed between the calculated |
423 |
frequency and the electric field, it was not a viable route for the |
424 |
calculation of a nitrile spectrum. Instead, the frequencies were taken |
425 |
and convolved together with a lorentzian line shape applied around the |
426 |
frequency value. These spectra are seen below in Figure |
427 |
4. While the spectrum without a field is lower in intensity and is |
428 |
almost bimodel in distrobution, the external field spectrum is much |
429 |
more unimodel. This tighter clustering has the affect of increasing the |
430 |
intensity around 2226 cm\textsuperscript{-1} where the peak is |
431 |
centered. The external field also has fewer frequencies of higher |
432 |
energy in the spectrum. Unlike the the zero field, where some frequencies |
433 |
reach as high as 2280 cm\textsuperscript{-1}. |
434 |
\begin{figure} |
435 |
\includegraphics[width=3.25in]{Convolved} |
436 |
\caption{Lorentzian convolved Gaussian frequencies of the zero field |
437 |
system (black) and the full field system (red)} |
438 |
\label{fig:Con} |
439 |
\end{figure} |
440 |
\section{Discussion} |
441 |
Interestingly, the field that is needed to switch the phase of 5CB |
442 |
macroscopically is larger than 1 V. However, in this case, only a |
443 |
voltage of 1.2 V was need to induce a phase change. This is impart due |
444 |
to the short distance of 5 nm the field is being applied across. At such a small |
445 |
distance, the field is much larger than the macroscopic and thus |
446 |
easily induces a field dependent phase change. However, this field |
447 |
will not cause a breakdown of the 5CB since electrochemistry studies |
448 |
have shown that it can be used in the presence of fields as high as |
449 |
500 V macroscopically. This large of a field near the surface of the |
450 |
elctrode would cause breakdown of 5CB if it could happen. |
451 |
|
452 |
The absence of any electric field dependency of the freuquency with |
453 |
the Gaussian simulations is interesting. A large base of research has been |
454 |
built upon the known tuning of the nitrile bond as the local field |
455 |
changes. This difference may be due to the absence of water or a |
456 |
molecule that induces a large internal field. Liquid water is known to have a very high internal field which |
457 |
is much larger than the internal fields of neat 5CB. Even though the |
458 |
application of Gaussian simulations followed by mapping it to |
459 |
some classical parameter is easy and straight forward, this system |
460 |
illistrates how that 'go to' method can break down. |
461 |
|
462 |
While this makes the application of nitrile Stark effects in |
463 |
simulations without water harder, these data show |
464 |
that it is not a deal breaker. The classically calculated nitrile |
465 |
spectrum shows changes in the spectra that will be easily seen through |
466 |
experimental routes. It indicates a shifted peak lower in energy |
467 |
should arise. This peak is a few wavenumbers from the leading edge of |
468 |
the larger peak and almost 75 wavenumbers from the center. This |
469 |
seperation between the two peaks means experimental results will show |
470 |
an easily resolved peak. |
471 |
|
472 |
The Gaussian derived spectra do indicate an applied field |
473 |
and subsiquent phase change does cause a narrowing of freuency |
474 |
distrobution. With narrowing, it would indicate an increased |
475 |
homogeneous distrobution of the local field near the nitrile. |
476 |
\section{Conclusions} |
477 |
Field dependent changes |
478 |
\newpage |
479 |
|
480 |
\bibliography{5CB} |
481 |
|
482 |
\end{doublespace} |
483 |
\end{document} |