Functional Nanocomposites
concentration c (e.g., 0.034 M for DA-EO
A R T I C L E S
o
5
) much less than the critical
follows. First, a cluster algorithm was developed to identify all
of the surfactants that make up a particular micelle. Next, we
determined the principal axes of the inertia tensor of the cluster.
The axis of the cylindrical micelle is determined as the inertia
tensor’s eigenvector with the smallest eigenvalue. We then
sought to represent each H4T4 surfactant by a vector. Although,
in principle, one could use end-to-end vectors, we decided that
a superior method is to again use the inertia tensor approach,
this time applied to each individual surfactant. In this fashion,
we can represent each DA (i.e., H T ) molecule by a unique
vector that gives the molecule’s orientation and is centered on
the molecule’s center of mass. Finally, a measure of orientational
order of the H4T4 surfactants in the cylindrical micelle was
obtained by forming the inner product of the molecule vector
and the vector describing the direction of cylinder axis. We
collected histograms of this inner product distribution, which
is, apart from a constant, the cosine of the angle between the
molecular axis and the micellar axis.
surfactant micelle concentration cmc, we use evaporative dip-coating,
spin-coating, or casting procedures to prepare thin films on silicon 〈100〉
7
or fused silica substrates. During deposition, preferential evaporation
of THF concentrates the depositing film in water and nonvolatile silica
and surfactant species. The progressively increasing surfactant con-
centration drives self-assembly of diacetylene/silica micelles and their
further organization into ordered liquid crystalline mesophases. Shape
and concentration of the DA surfactants influence the mesophase
obtained (lamellar, hexagonal, or cubic). Ultraviolet (UV) light-initiated
polymerization of the DA units, accompanied by catalyst-promoted
siloxane condensation, topochemically convert the colorless mesophase
into the blue PDA/silica nanocomposite, preserving the highly ordered,
self-assembled architecture (Figure 1).
4
4
Films were prepared by casting, spin-coating at 2000 rpm, or dip-
coating at a rate of 7.6-40 cm/min in an ambient atmosphere (23-25
°C, and relative humidity between 5 and 40%). Polymerization of PDA
to the blue form was done by UV exposure at 254 nm for times ranging
from 30 s to 30 min. Subsequent transformation to the red form was
accomplished by heating at 100 °C for times ranging from 30 s to 2
min or by exposure to solvent. A solvatochromic transition of blue
PDA/silica nanocomposite films to the corresponding red films was
accomplished by immersion of the blue films in the series of solvents,
hexane, 2-propanol, acetone, ethanol, methanol, or dimethylformamide.
Immersion times were 3 s followed by drying in nitrogen at room
temperature.
Molecular Modeling and Quantum Calculations. The
coarse-grained lattice calculations provide valuable insight into
the local packing and orientation of a collection of spontaneously
self-assembled DA molecules. To assess the effects of curvature
on the individual PDA oligomers, we undertook detailed
molecular simulations with detailed force fields to provide
realistic backbone structures whose electronic structures could
then be obtained from semiempirical quantum mechanical
calculations. Present computational limitations do not allow for
a detailed treatment of an entire micelle of (partially) polym-
erized DAs. Instead, we study single PDA oligomers. The
molecular configurations of single PDA oligomers were gener-
ated using molecular mechanics (MM) and dynamics (MD)
performed with the PolyGraf software package (version 3.21)
Modeling Section
Lattice Monte Carlo Simulation. The detailed packing and
orientational information of DA monomer surfactants in cylin-
drical micelles was studied with a 3D three-component lattice
8
model of a water/ethanol/surfactant system. As shown previ-
9,10
ously by Larson, this type of modeling is sufficiently coarse-
grained to allow for the sampling of the spontaneous self-
assembly process and the resulting 3D structures, yet sufficiently
detailed to investigate the orientation of DA within the cylindri-
cal micelle. Our lattice model represents small molecules (i.e.,
ethanol and water) as single sites on a simple cubic lattice with
interaction potentials based on the single fluid phase diagram
2
for isolated DA molecules and with Cerius (version 4.0) for
periodic calculations of molecules confined to silica pores. The
particular force field used for the atomic interactions was
Dreiding II with the dielectric constant set to 2.64. Partial
charges in the PDA oligomer were assigned using the charge
equilibration method. A MD simulation in the NVT ensemble
was performed for 500 ps, during which the temperature (T)
was fixed at 300 K. Structures were written out periodically at
0.1 ps intervals.
11
and employs Lorentz-Berthelot mixing rules for the cross
interactions.11 The surfactant molecules are represented by a
polymer chain, a collection of beads with various patterns of
connectivity. Here we represent the DA molecule by a simple
linear chain consisting of four hydrophilic beads (“heads”, H)
and four hydrophobic beads (“tails”, T), or H4T4. NVT Monte
Semiempirical quantum calculations were performed on
selected saved structures. To reduce the computational effort,
we first “trimmed” the PDA structures. That is, the DA
backbone of the PDA oligomers was terminated at the first
methylene group of each side chain by substituting a hydrogen
atom in its place. On the basis of the trimmed structures obtained
from the MD simulation, we calculated UV-visible spectra of
the PDA oligomers using the transition wavelengths and
3
Carlo (MC) simulations on a lattice of size 32 with periodic
boundary conditions were started from a random configuration.
7
After an initial equilibration of 10 attempted MC moves at high
7
temperature, the temperature was dropped in five stages of 10 -
8
1
0 attempted MC moves, and subsequent self-assembly was
followed until equilibrium was reached after approximately
8
3
× 10 MC moves. The initial composition was located in an
area of the three-component phase diagram of H4T4 that our
previous simulations had identified as the hexagonal phase. The
description of the orientation of H4T4 surfactants in cylindrical
micelles requires an objective analysis and was determined as
oscillator strengths obtained from semiempirical ZINDO/S
12,13
calculations
performed using the HyperChem 5.02 software
package, which provides the orbitals, energies, and transition
energies and dipoles. ZINDO/S is known to give accurate
spectroscopic predictions for these types of molecules when a
sufficient number of HOMOs and LUMOs are included in the
configuration interaction part of the calculation, as is the case
here.
(
7) Lu, Y. F.; Yang, Y.; Sellinger, A.; Lu, M. C.; Huang, J. M.; Fan, H. Y.;
Haddad, R.; Lopez, G.; Burns, A. R.; Sasaki, D. Y.; Shelnutt, J. A.; Brinker,
C. J. Nature 2001, 410, 913.
(8) Rankin, S. E.; Malanoski, A. P.; van Swol, F. Mater. Res. Soc. Proc. 2001,
6
36, D1.2.1.
(
9) Larson, R. G. J. Chem. Phys. 1988, 89, 1642.
(
10) Larson, R. G. J. Phys. II France 1996, 6, 1441.
11) Hansen, J.-P.; McDonald, I. R. Theory of Simple Liquids, 2nd ed.; Academic
Press: London, England, 1986.
(
(12) Ridley, J. E.; Zerner, M. C. Theor. Chim. Acta 1973, 32, 111.
(13) Bacon, A.; Zerner, M. C. Theor. Chim. Acta 1979, 53, 21.
J. AM. CHEM. SOC.
9
VOL. 125, NO. 5, 2003 1271