statistical mechanical modeling of chromophore interactions at high concentration
Skip this Video
Download Presentation
Statistical Mechanical Modeling of Chromophore Interactions at High Concentration

Loading in 2 Seconds...

play fullscreen
1 / 70

Statistical Mechanical Modeling of Chromophore Interactions at High Concentration - PowerPoint PPT Presentation

  • Uploaded on

Alvin L Kwiram Symposium Electrical, Optical and Magnetic Properties of Organic and Hybrid Materials. Statistical Mechanical Modeling of Chromophore Interactions at High Concentration. Bruce H. Robinson Robert D. Nielsen Harrison Rommel. Outline.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about ' Statistical Mechanical Modeling of Chromophore Interactions at High Concentration' - stacie

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
statistical mechanical modeling of chromophore interactions at high concentration

Alvin L Kwiram Symposium

Electrical, Optical and Magnetic Properties of Organic and Hybrid Materials

Statistical Mechanical Modeling of Chromophore Interactions at High Concentration

Bruce H. Robinson

Robert D. Nielsen

Harrison Rommel

  • Fundamental Concepts: Matter organization is a many body problem
  • Why is order important for NLO materials
  • Many different patterns of arrangement are possible in Materials
  • Treating molecules as classical objects.
  • The principles of Statistical Mechanics to determine organization
  • Monte Carlo Methods to determine overall order in frustrated systems
  • Analytic Methods to compute average order
  • Application to EO materials and the EO coefficient
organization of molecules in condensed phases
Organization of Molecules in Condensed phases
  • Put the NLO molecules (aka chromophores) into a host.
  • Hosts: neat, solvent, polymer, block-copolymer
  • Condensed phases means liquid or solid states
  • Most chromophores in a polymer slurry are glassy and have a solid to glass transition temperature (Tg).
  • Usually: No real lattice or lattice order (not crystals)
why is order important for nlo materials the electro optic effect

Why is order important for NLO materials The Electro-Optic Effect

To generate a voltage across the material when shining light we need to have the molecules somewhat aligned or asymmetrically arranged.

We now proceed to explain

how and why this is so.

how do we get asymmetric arrangements
How do we get asymmetric arrangements:
  • Must start with asymmetric molecules
  • use lager scale organization (dendrimers host-polymers)
  • Grow crystals
  • Generate ferro-electric phases
  • Generate ordered phases that give higher order than “normal" para-electric or anti-ferroelectric order
  • Generate some order using a poling or ordering Electric (E) field (place sample in capacitor):
why do we get an e o effect for organic nlo materials
Why do we get an E.O. effect for Organic NLO materials
  • The electric field of the light pushes and pulls electrons in the NLO chromophore.
  • The movement of electrons sets up a field (as a voltage) around the material
  • The distribution of charges can be detected by a capacitor type arrangement.
  • The response of the material to light (or a DC voltage) is equivalent to a change in the index of refraction of the material
electro optic eo coupling
Electro Optic (EO) Coupling

Explaining the Phenomenon:

  • Molecules (arrows in the box) are somewhat aligned on z pointing up.
  • r is the EO coupling coefficient – coupling electric field effects to the index of refraction changes.
  • The applied E field (along z) changes the index of refraction on either z or x depending on the nature of the coupling coefficient r.
  • Coupling along z causes a change in the index of refraction that then causes a phase delay in the light traveling along the x direction (which has its E field in the z direction).
measurement and functionality
Measurement and Functionality
  • Phase of the light in the two arms is changed by the driving E field.
  • When the phases differ by 180, or pi, no light out on z.
  • Low value of driving E field or
  • Large dynamic range
  • Low insertion loss (<1 db)
large angle fast response spatial light modulator slm
Large Angle, Fast Response Spatial Light Modulator (SLM)

Experimental Results

Schematic Diagram

12 mrad/


Literature Citations

• Steier, et al., “Polymeric waveguide prism based electro-optic beam deflector,” Opt. Eng., 40, 1217-22 (2001)

• Steier, et al., “Beam deflection with electro-optic polymer waveguide prism array,” Proc. SPIE, 3950, 108-116 (2000)

• Steier, et al., “Polymeric waveguide beam deflector for electro-optic switching,” Proc. SPIE, 4279, 37-44 (2001)

phased array radar with photonic phase shifter
Phased Array Radar with Photonic Phase Shifter

Steier, et al., IEEE mW & Guided Wave Lett. 9, 357 (1999)

  • Engineering goal is to optimize macroscopic E.O. effect

1) Bandwidth limitation

2) Push-Pull modulation

  • What do we optimize on molecular scale

to yield best effect?

microscopic basis of e o effect for organic nlo materials
Microscopic basis of E.O. effect for organic NLO materials
  • Large molecular Hyperpolarizability
  • Classical perturbation picture
  • Quantum perturbation picture
how a single molecule responds to light
How a single molecule responds to light
  • E field of the light interacts with the electron distribution in the molecule to change charge distributions.
  • Molecules have no net charge, but a very large dipole moment. For two charges (plus and minus q) separated by a distance r: the dipole moment is
microscopic basis of e o effect
Microscopic basis of E.O. effect

The E field causes polarization and hyperpolarization

For the most part, consider effect on z, in molecule

molecular hyperpolarizability b
Molecular Hyperpolarizability (b)

Quantum Mechanical description

In the molecule, along z.

microscopic to macroscopic
Microscopic to Macroscopic

Net acentric ordering to give macroscopic Hyperpolarizability

The EO coefficient is proportional to molecular hyperpolarizability and the net order, which, we will see, is related to the intrinsic dipole moment.

summary organic nlo materials
Summary: Organic NLO materials
  • Couple applied electric fields to optical transmission
  • Solid state, Guest/Host, Chromophore/polymer, matrix
  • High molecular electric dipole moment ( >5 D)
  • Large molecular Hyperpolarizability ( )
  • Organic Chromophores aligned by “Poling” with electric field above glass transition temperature of Guest Polymer followed by cooling
  • Experimental observation: At high Chromophore concentrations poling becomes inefficient
how does a poling field give order
How does a poling field give order?
  • The dipole of the molecule interacts with the E field of the poling field
  • How a single molecule interacts with E:


Probability of being at a particular orientation q, is given by Boltzmann’s formula

modeling an nlo chromophore
Modeling an NLO Chromophore

We need to model the structures with simpler, classical rules. When considering many molecular interactions Q. M. too time consuming.

Chromophore: DMC3-97

  • Model Chromophores as Ellipsoids
  • Each ellipsoid has a dipole moment at its center along the major axis.
  • Energy due of E-field induced polarization will be neglected.
  • Leonard Jones type ellipsoidal repulsions will be included.
bolzmann s rule
Bolzmann’s Rule
  • G is the Gibbs-Distribution or the Generating Function.
  • Probability of being at a particular orientation q:
  • Need the integral at the bottom to be sure probabilities sum to 1.

The strength of the interaction is relative to temperature, as a single parameter, x

how much order can we get
How much order can we get
  • Problem: At low T a molecule would align but it could not move to align
  • At T>Tg molecules can rotate but will not perfectly align.
  • Example:
  • One electron displaced by 1 Angstrom from a proton has a dipole of 4.8 Debye
  • For a 10 Debye dipole: F=0.7
  • So you get some, but not a lot of order
boltzmann probability distribution
Boltzmann Probability Distribution
  • Probability as a function of the cosine of the angle shows most probable orientation is to be aligned but but it is not a lot.
  • For F=0.7

and 2.0

Area under the two curves are the same.



computing average quantities
Computing Average Quantities
  • We want to know how much, on average, is aligned.
  • With the generating function, G, one can compute any quantity.
  • Example, the mean cos angle raised to the nth power.
averages how much order can you get
Averages: How much order can you get?
  • Order as a function of dipole strength:
  • Maximum possible average would be 1
  • In practice won’t get above 0.6

Low Strength limit: F/3


the influence of number density
The influence of Number Density
  • Dipoles form anti-ferroelectric state at high density.
  • Local Dipole-Dipole interaction fights poling-field induced order
  • Local steric forces restrict order
  • We investigated Number Density effect by:

Analytic mean field theory

MC simulations of dipole orientations

analytic theory of piekara
Analytic Theory of Piekara

Piekara, A. (1939) Proc. R. Soc. London A 172, 360-381

  • The general Statistical Mechanical Problem is a complicated (N-body) interaction.
  • Probability distribution of each chromophore, G(1), in a field of many other chromophores, and poling field.
analytic theory of piekara cont d
Analytic Theory of Piekara (cont’d)
  • Treat poling field E as interacting with dipoles:
  • Have an effective interaction with all other chromophores, within a “crystallite”
  • Average over all “crystallites” uniformly.
comparison with low chromophore density langevintheory
Comparison with Low Chromophore Density LangevinTheory
  • Langevin Theory has the partition function:
  • Piekara has crystallite partition function:
piekara s w
Piekara’s W
  • F and W directly compete
  • Problem: W is not known, nor is it self-consistent
  • Hypothesis: W is related to dipole-dipole energy of interaction:

Piekara’s limit

maximizing e o effect1
Maximizing E.O. Effect

There is always a maximum loading, which is very simply related to the dipole moment strength, and nearly independent of the poling field at typical fields

ratio of e o coeffcients
Ratio of E.O. Coeffcients

Propagation along controlling field

  • Increasing dipole moment alone does not give larger E.O. effect
  • Steric forces indicate spherical shape is optimal
  • is not constant, but is field dependent
  • Aggregation of Chromophores not necessary to explain E.O. roll-off
  • Can the choice of W be justified?
monte carlo simulations

Monte Carlo Simulations

  • A straightforward way to compute averages
  • A Justification for Piekara’s W
monte carlo calculations
Monte Carlo Calculations
  • Use Monte Carlo methods to determine the effect of dipolar interactions between chromophores.
  • Place dipoles on a grid (simple cubic lattice and body centered cubic lattice)
  • M by M by M array with r as nearest neighbor distance.
  • Periodic boundary condition

A 5 by 5 two dimensional array

Randomly oriented dipoles

how monte carlo works
How Monte Carlo Works
  • Choose a dipole
  • Rotate dipole by: a rotation axis and angle, selected randomly
  • Compare the energy before and after rotation.

If the energy is higher, compare Boltzmann Probability with a [0,1] random number, and keep if larger.

If the energy is lower, keep the move

summary of mc calculation results
Summary of MC Calculation Results
  • Planes Statistically Oppose One Another for cubic lattice, and represent long range order.
  • The dominant angular distribution depends on dipole moment to the Fourth Power
  • The dominant angular distribution depends on chromophore density to the Second Power.
  • MC predictions not sensitive to Cubic vs. BCC
  • MC statistics depend upon lattice distortion and are a first step to understanding steric effects
  • Scaling of W is needed to compare with analytic theory
  • Strongly anti-ferro-electric interactions resist poling field in condensed matter. This is the source of roll-off of EO response with chromophore number density.
  • Aggregation model not necessary
  • Piekara’s crystallite interaction, W, can be interpreted as a field and competes with poling in a simple geometric way.
  • As the dipole moment increases, the maximum loading decreases.
  • For a given dipole moment there is an optimum number density that is easily calculated from the analytic theory.
  • Spherical shaped chromophores are the optimum geometry (better than either prolate or oblate ellipsoids)
  • MC shows the expansion of N body interactions has a cosine term as the dominant distribution. You will never get this from a small number of chromophores.
  • MC results reproduce W dependence on ,N,T.
  • Include steric effect in MC calculations
  • Generalize from point Dipoles to charges
  • Explore extended multi-chromophore structures
microscopic basis of e o effect for organic nlo materials1
Microscopic basis of E.O. effect for organic NLO materials
  • Large molecular Hyperpolarizability
  • Net acentric ordering to give macroscopic Hyperpolarizability