Design and synthesis of chromophores with enhanced electro-optic activities in both bulk and plasmonic-organic hybrid devices.

This study demonstrates enhancement of in-device electro-optic activity via a series of theory-inspired organic electro-optic (OEO) chromophores based on strong (diarylamino)phenyl electron donating moieties. These chromophores are tuned to minimize trade-offs between molecular hyperpolarizability and optical loss. Hyper-Rayleigh scattering (HRS) measurements demonstrate that these chromophores, herein described as BAH, show >2-fold improvement in ß versus standard chromophores such as JRD1, and approach that of the recent BTP and BAY chromophore families. Electric field poled bulk devices of neat and binary BAH chromophores exhibited significantly enhanced EO coefficients (r33) and poling efficiencies (r33/Ep) compared with state-of-the-art chromophores such as JRD1. The neat BAH13 devices with charge blocking layers produced very large poling efficiencies of 11.6 ± 0.7 nm2 V-2 and maximum r33 value of 1100 ± 100 pm V-1 at 1310 nm on hafnium dioxide (HfO2). These results were comparable to that of our recently reported BAY1 but with much lower loss (extinction coefficient, k), and greatly exceeding that of other previously reported OEO compounds. 3 : 1 BAH-FD : BAH13 blends showed a poling efficiency of 6.7 ± 0.3 nm2 V-2 and an even greater reduction in k. 1 : 1 BAH-BB : BAH13 showed a higher poling efficiency of 8.4 ± 0.3 nm2 V-2, which is approximately a 2.5-fold enhancement in poling efficiency vs. JRD1. Neat BAH13 was evaluated in plasmonic-organic hybrid (POH) Mach-Zehnder modulators with a phase shifter length of 10 µm and slot widths of 80 and 105 nm. In-device BAH13 achieved a maximum r33 of 208 pm V-1 at 1550 nm, which is ∼1.7 times higher than JRD1 under equivalent conditions.

Electro-Optic Activity in Excess of 1000 pm V-1 Achieved via Theory-Guided Organic Chromophore Design.

High performance organic electro-optic (OEO) materials enable ultrahigh bandwidth, small footprint, and extremely low drive voltage in silicon-organic hybrid and plasmonic-organic hybrid photonic devices. However, practical OEO materials under device-relevant conditions are generally limited to performance of ≈300 pm V-1 (10× the EO response of lithium niobate). By means of theory-guided design, a new series of OEO chromophores is demonstrated, based on strong bis(4-dialkylaminophenyl)phenylamino electron donating groups, capable of EO coefficients (r33 ) in excess of 1000 pm V-1 . Density functional theory modeling and hyper-Rayleigh scattering measurements are performed and confirm the large improvement in hyperpolarizability due to the stronger donor. The EO performance of the exemplar chromophore in the series, BAY1, is evaluated neat and at various concentrations in polymer host and shows a nearly linear increase in r33 and poling efficiency (r33 /Ep , Ep is poling field) with increasing chromophore concentration. 25 wt% BAY1/polymer composite shows a higher poling efficiency (3.9 ± 0.1 nm2 V-2 ) than state-of-the-art neat chromophores. Using a high-ε charge blocking layer with BAY1, a record-high r33 (1100 ± 100 pm V-1 ) and poling efficiency (17.8 ± 0.8 nm2 V-2 ) at 1310 nm are achieved. This is the first reported OEO material with electro-optic response larger than thin-film barium titanate.

The Importance of Assay Imprecision near the Screen Cutoff for Newborn Screening of Lysosomal Storage Diseases.

For newborn screening (NBS) of lysosomal storage diseases, programs measure enzymatic activities in dried blood spots (DBS) and, in most cases, act on samples where the measurement is below a specific cutoff value. The rate of false positives and negatives in any NBS program is of critical importance. The measured values across a population of newborns are governed by many factors, and in this article we focus on assay imprecision. Assay parameters including the Analytical Range and the Z-Factor have been discussed as a way to compare assay performance for NBS of lysosomal storage diseases. Here we show that these parameters are not rigorously connected to the rate of false positives and negatives. Rather, it is the assay imprecision near the screen cutoff that is the most important parameter that determines the rate of false positives and negatives. We develop the theoretical treatment of assay imprecision and how it is linked to screen performance. What emerges is a useful type of parametric plot that allows for rigorous assessment of the effect of assay imprecision on the rate of false positives and false negatives that is independent of the choice of screen cutoff value. Such plots are useful in choosing cutoff values. They also show that a high assay imprecision cannot be overcome by changing the cutoff value or by use of postanalysis, statistical tools. Given the importance of assay imprecision near the cutoff, we propose that quality control DBS are most useful if they span a range of analyte values near the cutoff. Our treatment is also appropriate for comparing the performance of multiple assay platforms that each measure the same quantity (i.e., the enzymatic activity in DBS). The analysis shows that it is always best to use the assay platform that gives the lowest imprecision near the cutoff.

Molecular Engineering of Structurally Diverse Dendrimers with Large Electro-Optic Activities.

To boost electro-optic (EO) performance, a series of multichromophore dendrimers have been developed based on higher hyperpolarizability (CLD-type) chromophore cores that have been used previously (FTC-type dendrimers). The multichromophore dendrimers were molecularly engineered to have either three arms, two arms, or one arm; long or short linkers; and a fluorinated dendron (FD) or tert-butyldiphenylsilyl (TBDPS) shell. The EO performance obtained by FDSD (poling efficiency = 1.60 nm2 V-2), based on succinic diester linkers, was higher than the analogue with longer adipic diester linkers and higher than the analogs with fewer chromophore moieties. Due to the shorter succinic diester linker and improved site isolation, the dendrimer chromophore with TBDPS groups exhibited enhanced glass-transition temperature ( Tg = 108 °C) and comparable poling efficiency (1.62 nm2 V-2) to the FD-containing version. These neat EO dendrimers have a higher index of refraction ( n = 1.75-1.84 at 1310 nm) than guest-host polymeric EO materials ( n ≈ 1.6, 1310 nm) and FTC-type EO dendrimers ( n = 1.73, 1310 nm), which is important, because a key metric for Mach-Zehnder modulators is proportional to n3. In addition, binary chromophore organic glasses (BCOGs) were prepared by doping a secondary EO chromophore at 25 wt % into neat dendrimers. Enhancements of EO performance were found in all BCOG materials compared with neat dendrimers due to the effect of blending. As a result of increased chromophore density, the n values of the BCOGs improved to 1.81-1.92. One BOCG, in particular, displayed the highest poling efficiency (2.35 nm2 V-2) and largest EO coefficient ( r33) value of 275 pm V-1 at 1310 nm, which represents a high n3 r33 figure-of-merit of 1946 pm V-1. The high poling efficiencies and n3 r33 figure-of-merit combined with excellent film forming confirm these neat dendrimers and BCOGs based on them as promising candidates for incorporation into photonic devices.

DANPY (dimethylaminonaphthylpyridinium): an economical and biocompatible fluorophore.

Dyes with nonlinear optical (NLO) properties enable new imaging techniques and photonic systems. We have developed a dye (DANPY-1) for photonics applications in biological substrates such as nucleic acids; however, the design specification also enables it to be used for visualizing biomolecules. It is a prototype dye demonstrating a water-soluble, NLO-active fluorophore with high photostability, a large Stokes shift, and a favorable toxicity profile. A practical and scalable synthetic route to DANPY salts has been optimized featuring: (1) convergent Pd-catalyzed Suzuki coupling with pyridine 4-boronic acid, (2) site-selective pyridyl N-methylation, and (3) direct recovery of crystalline intermediates without chromatography. We characterize the optical properties, biocompatibility, and biological staining behavior of DANPY-1. In addition to stability and solubility across a range of polar media, the DANPY-1 chromophore shows a first hyperpolarizability similar to common NLO dyes such as Disperse Red 1 and DAST, a large two-photon absorption cross section for its size, substantial affinity to nucleic acids in vitro, an ability to stain a variety of cellular components, and strong sensitivity of its fluorescence properties to its dielectric environment.

Detection of mucopolysaccharidosis III-A (Sanfilippo Syndrome-A) in dried blood spots (DBS) by tandem mass spectrometry.

BACKGROUND: With ongoing efforts to develop improved treatments for Sanfilippo Syndrome Type A (MPS-IIIA), a disease caused by the inability to degrade heparan sulfate in lysosomes, we sought to develop an enzymatic activity assay for the relevant enzyme, sulfamidase, that uses dried blood spots (DBS). METHODS: We designed and synthesized a new sulfamidase substrate that can be used to measure sulfamidase activity in DBS using liquid chromatography-tandem mass spectrometry (LC-MS/MS). RESULTS: Sulfamidase activity was readily detected in DBS using the new substrate and LC-MS/MS. Sulfamidase activity showed acceptable linearity proportional to the amount of enzyme and reaction time. Sulfamidase activity in 238 random newborns was well elevated compared to the range of activities measured in DBS from 8 patients previously confirmed to have MPS-IIIA. CONCLUSIONS: This is the first report of an assay capable of detecting sulfamidase in DBS. The new assay could be useful in diagnosis and potentially for newborn screening of MPS-IIIA.

Simple Model for the Benzene Hexafluorobenzene Interaction.

While the experimental intermolecular distance distribution functions of pure benzene and pure hexafluorobenzene are well described by transferable all-atom force fields, the interaction between the two molecules (in a 1:1 mixture) is not well simulated. We demonstrate that the parameters of the transferable force fields are adequate to describe the intermolecular distance distribution if the charges are replaced by a set of charges that are not located at the atoms. The simplest model that well describes the experimental distance distribution, between benzene and hexafluorobenzene, is that of a single ellipsoid for each molecule, representing the van der Waals interactions, and a set of three point charges (on the axis perpendicular to the arene plane) which give the same quadrupole moment as do the all atom charges from the transferable force fields.

Systematic Generation of Anisotropic Coarse-Grained Lennard-Jones Potentials and Their Application to Ordered Soft Matter.

We have developed an approach to coarse-grained (CG) modeling of the van der Waals (vdW) type of interactions among molecules by representing groups of atoms within those molecules in terms of ellipsoids (rather than spheres). Our approach systematically translates an arbitrary underlying all-atom (AA) representation of a molecular system to a multisite ellipsoidal potential within the family of Gay-Berne type potentials. As the method enables arbitrary levels of coarse-graining, or even multiple levels of coarse-graining within a single simulation, we describe the method as a Level of Detail (LoD) model. The LoD model, as integrated into our group's Metropolis Monte Carlo computational package, is also capable of reducing the complexity of the molecular electrostatics by means of a multipole expansion of charges obtained from an AA force field (or directly from electronic structure calculations) of the charges within each ellipsoid. Electronic polarizability may additionally be included. The present CG representation does not include transformation of bonded interactions; ellipsoids are connected at the fully atomistic bond sites by freely rotating links that are constrained to maintain a constant distance. The accuracy of the method is demonstrated for three distinct types of self-assembling or self-organizing molecular systems: (1) the interaction between benzene and perfluorobenzene (dispersion interactions), (2) linear hydrocarbon chains (a system with large conformational flexibility), and (3) the self-organization of ethylene carbonate (a highly polar liquid). Lastly, the method is applied to the interaction of large (∼100 atom) molecules, which are typical of organic nonlinear optical chromophores, to demonstrate the effect of different CG models on molecular assembly.

Dielectric and phase behavior of dipolar spheroids.

The Stockmayer fluid, composed of dipolar spheres, has a well-known isotropic-ferroelectric phase transition at high dipole densities. However, there has been little investigation of the ferroelectric transition in nearly spherical fluids at dipole densities corresponding to those found in many polar solvents and in guest-host organic electro-optic materials. In this work, we examine the transition to ordered phases of low-aspect-ratio spheroids under both unperturbed and poled conditions, characterizing both the static dielectric response and thermodynamic properties of spheroidal systems. Spontaneous ferroelectric ordering was confined to a small region of aspect ratios about unity, indicating that subtle changes in sterics can have substantial influence on the behavior of coarse-grained liquid models. Our results demonstrate the importance of molecular shape in obtaining even qualitatively correct dielectric responses and provide an explanation for the success of the Onsager model as a phenomenological representation for the dielectric behavior of polar organic liquids.

Relation of system dimensionality and order parameters.

Orientational order parameters are useful metrics for characterizing the probability distribution for vector-valued quantities such as the dipole moment or optical axis of molecules in materials such as liquid crystals and organic glasses. These parameters are the moments of the underlying orientational probability distribution. Many molecular systems can be characterized using a single centrosymmetric (even) moment. For dipolar systems, an applied electric or magnetic field can break the symmetry of the system, leading to nonzero acentric (odd) moments. For complex systems, it is difficult to characterize the nature of the bulk structures and to quantitatively understand the relationship between acentric and centrosymmetric moments. We have found that it is useful to relate the moments of the distribution in terms of an apparent dimensionality of the ordering process. Here we show that the idea of noninteger dimensionality, originally introduced by Stillinger, provides a useful method to characterize the relation between centrosymmetric and acentric orientational order parameters. Applying dimensional constraints is equivalent to removing rotational degrees of freedom or constraining rotation within a restricted volume. Simulations based on simple examples­using restoring potentials on arrays of independent dipoles­and on complex many-body Monte Carlo simulations of dipolar spheroids are described. An analysis of the results illustrates the utility of fractional dimensionality to describe ordering in materials.

Optimizing calculations of electronic excitations and relative hyperpolarizabilities of electrooptic chromophores.

CONSPECTUS: Organic glasses containing chromophores with large first hyperpolarizabilities (ß) are promising for compact, high-bandwidth, and energy-efficient electro-optic devices. Systematic optimization of device performance requires development of materials with high acentric order and enhanced hyperpolarizability at operating wavelengths. One essential component of the design process is the accurate calculation of optical transition frequencies and hyperpolarizability. These properties can be computed with a wide range of electronic structure methods implemented within commercial and open-source software packages. A wide variety of methods, especially hybrid density-functional theory (DFT) variants have been used for this purpose. However, in order to provide predictions useful to chromophore designers, a method must be able to consistently predict the relative ordering of standard and novel materials. Moreover, it is important to distinguish between the resonant and nonresonant contribution to the hyperpolarizabiliy and be able to estimate the trade-off between improved ß and unwanted absorbance (optical loss) at the target device's operating wavelength. Therefore, we have surveyed a large variety of common methods for computing the properties of modern high-performance chromophores and compared these results with prior experimental hyper-Rayleigh scattering (HRS) and absorbance data. We focused on hybrid DFT methods, supplemented by more computationally intensive Møller-Plesset (MP2) calculations, to determine the relative accuracy of these methods. Our work compares computed hyperpolarizabilities in chloroform relative to standard chromophore EZ-FTC against HRS data versus the same reference. We categorized DFT methods used by the amount of Hartree-Fock (HF) exchange energy incorporated into each functional. Our results suggest that the relationship between percentage of long-range HF exchange and both ßHRS and λmax is nearly linear, decreasing as the fraction of long-range HF exchange increases. Mild hybrid DFT methods are satisfactory for prediction of λmax. However, mild hybrid methods provided qualitatively incorrect predictions of the relative hyperpolarizabilities of three high-performance chromophores. DFT methods with approximately 50% HF exchange, and especially the Truhlar M062X functional, provide superior predictions of relative ßHRS values but poorer predictions of λmax. The observed trends for these functionals, as well as range-separated hybrids, are similar to MP2, though predicting smaller absolute magnitudes for ßHRS. Frequency dependence for ßHRS can be calculated using time-dependent DFT and HF methods. However, calculation quality is sensitive not only to a method's ability to predict static hyperpolarizability but also to its prediction of optical resonances. Due to the apparent trade-off in accuracy of prediction of these two properties and the need to use static finite-field methods for MP2 and higher-level hyperpolarizability calculations in most codes, we suggest that composite methods could greatly improve the accuracy of calculations of ß and λmax.

Optimum Exchange for Calculation of Excitation Energies and Hyperpolarizabilities of Organic Electro-optic Chromophores.

Organic electro-optic (OEO) materials integrated into silicon-organic hybrid devices afford significant improvements in size, weight, power, and bandwidth performance of integrated electronic/photonic systems critical for current and next generation telecommunication, computer, sensor, transportation, and defense technologies. Improvement in molecular first hyperpolarizability (ß), and in turn electro-optic activity, is crucial to optimizing device performance. Common hybrid density functional theory (DFT) methods, while attractive due to their computational scaling, often perform poorly for optical properties in systems with substantial intramolecular charge-transfer character, such as OEO chromophores. This study evaluates the utility of the long-range corrected (LC) DFT methods for computation of the molecular second-order nonlinear optical response. We compare calculated results for a 14-molecule benchmark set of OEO chromophores with the corresponding experimentally measured ß and one-photon absorption energy, λmax. We analyze the distance dependence of the fraction of exact exchange in LC-DFT methods for accurately computing these properties for OEO chromophores. We also examine systematic tuning of the range-separation parameter to enforce Koopmans'/ionization potential theorem. This tuning method improves prediction of excitation energies but is not reliable for predicting the hyperpolarizabilities of larger chromophores since the tuning parameter value can be too small, leading to instabilities in the computation of ßHRS. Additionally, we find that the size dependence of the optimal tuning parameter for the ionization potential has the opposite size dependence of optimal tuning parameter for best agreement with the experimental λmax, suggesting the tuning for the ionization potential is unreliable for extended conjugated systems.

Direct simulation of magnetic resonance relaxation rates and line shapes from molecular trajectories.

We simulate spin relaxation processes, which may be measured by either continuous wave or pulsed magnetic resonance techniques, using trajectory-based simulation methodologies. The spin-lattice relaxation rates are extracted numerically from the relaxation simulations. The rates obtained from the numerical fitting of the relaxation curves are compared to those obtained by direct simulation from the relaxation Bloch-Wangsness-Abragam-Redfield theory (BWART). We have restricted our study to anisotropic rigid-body rotational processes, and to the chemical shift anisotropy (CSA) and a single spin-spin dipolar (END) coupling mechanisms. Examples using electron paramagnetic resonance (EPR) nitroxide and nuclear magnetic resonance (NMR) deuterium quadrupolar systems are provided. The objective is to compare those rates obtained by numerical simulations with the rates obtained by BWART. There is excellent agreement between the simulated and BWART rates for a Hamiltonian describing a single spin (an electron) interacting with the bath through the chemical shift anisotropy (CSA) mechanism undergoing anisotropic rotational diffusion. In contrast, when the Hamiltonian contains both the chemical shift anisotropy (CSA) and the spin-spin dipolar (END) mechanisms, the decay rate of a single exponential fit of the simulated spin-lattice relaxation rate is up to a factor of 0.2 smaller than that predicted by BWART. When the relaxation curves are fit to a double exponential, the slow and fast rates extracted from the decay curves bound the BWART prediction. An extended BWART theory, in the literature, includes the need for multiple relaxation rates and indicates that the multiexponential decay is due to the combined effects of direct and cross-relaxation mechanisms.

Nano-engineering lattice dimensionality for a soft matter organic functional material.

A high performing electro-optic (EO) chromophore with covalently attached coumarin-based pendant groups exhibits intermolecular correlation of coumarin units through molecular dynamics (MD) simulations. Unique, orthogonal molecular orientations of the chromophore and coumarin units are also evident when investigated optically. Such molecular orientation translates to reduced lattice dimensionality of the bulk C1 soft matter material system, leading to increased acentric order and EO activity. Results are corroborated by nanorheological experimental methods.

Solvents level dipole moments.

The dipole moments of highly polar molecules measured in solution are usually smaller than the molecular dipole moments that are calculated with reaction field methods, whereas vacuum values are routinely calculated in good agreement with available vapor phase data. Whether from Onsager's theory (or variations thereof) or from quantum mechanical methods, the calculated molecular dipoles in solution are found to be larger than those measured. The reason, of course, is that experiments measure the net dipole moment of solute together with the polarized (perturbed) solvent "cloud" surrounding it. Here we show that the reaction field charges that are generated in the quantum mechanical self-consistent reaction field (SCRF) method give a good estimate of the net dipole moment of the solute molecule together with the moment arising from the reaction field charges. This net dipole is a better description of experimental data than the vacuum dipole moment and certainly better than the bare dipole moment of the polarized solute molecule.

Dielectric dependence of the first molecular hyperpolarizability for electro-optic chromophores.

Experimental and computational studies of the solvent dependence of the first molecular hyperpolarizability (ß) for two donor-bridge-acceptor chromophores (CLD-1 and YLD156) are presented. Hyper-Rayleigh scattering (HRS) measurements are performed with 1907 nm excitation in a series of solvents with dielectric constants ranging from ~2 (toluene) to ~36 (acetonitrile). For both chromophores an approximately 2-fold increase in ß is observed by HRS over this range of dielectric constants. Computational studies employing a polarized continuum model to represent the solvent are capable of reproducing this experimental result. The experimental and computational results are compared to the predictions of the widely employed two-state model (TSM) for ß. Surprisingly, for the chromophores studied here the TSM predicts that ß should decrease with increasing dielectric constant over the range investigated. The results presented here demonstrate that the TSM provides neither a quantitative nor qualitative description of the solvent dependence of ß for CLD-1 and YLD156. The enhancement of ß with increased dielectric constant suggests that modification of the dielectric surrounding the chromophore is one path by which the performance of nonlinear optical devices employing these chromophores may be significantly enhanced.

Measuring order in contact-poled organic electrooptic materials with variable-angle polarization-referenced absorption spectroscopy (VAPRAS).

Organic nonlinear electrooptical (ONLO) chromophores must be acentrically ordered for the ONLO material to have electrooptic (EO) activity. The magnitude of the order is characterized by the acentric order parameter, , where ß is the major Euler angle between the main axis of the chromophore and the poling field which imposes the acentric order. The acentric order parameter, which is difficult to measure directly, is related to the centrosymmetric order parameter, defined as = ½(3-1), through the underlying statistical distribution. We have developed a method to determine centrosymmetric order of the ONLO chromophores when the order is low (i.e., < 0.1). We have extended the method (begun by Graf et al. J. Appl. Phys. 1994, 75, 3335.) based on the absorption of light to determine the centrosymmetric order parameter induced by a poling field on a thin film sample of ONLO material. We find that the order parameters, analyzed by two different methods, are similar and also consistent with theoretical estimates from modeling of the system using coarse-grained Monte Carlo statistical mechanical methods.

Reduced dimensionality in organic electro-optic materials: theory and defined order.

Identification of electronic intermolecular electrostatic interactions that can significantly enhance poling-induced order is important to the advancement of the field of organic electro-optics. Here, we demonstrate an example of such improvement achieved through exploitation of the interaction of coumarin pendant groups in chromophore-containing macromolecules. Acentric order enhancement is explained in terms of lattice-symmetry effects, where constraint of orientational degrees of freedom alters the relationship between centrosymmetric and acentric order. We demonstrate both experimentally and theoretically that lattice dimensionality can be defined using the relationship between centrosymmetric order and acentric order. Experimentally: Acentric order is determined by attenuated total reflection measurement of electro-optic activity coupled with hyper-Rayleigh scattering measurement of molecular first hyperpolarizability, and centrosymmetric order is determined by the variable angle polarization referenced absorption spectroscopy method. Theoretically: Order is determined from statistical mechanical models that predict the properties of soft condensed matter.

Dielectric constants of simple liquids: stockmayer and ellipsoidal fluids.

Coarse-grained models of molecular interactions are of interest because they convey the essence of molecular interactions in simple and easy to understand form. However, coarse-grained models fail to adequately predict some material properties, such as the failure of the Stockmayer model to reproduce the dielectric behavior of highly polar liquids. We examine the behavior of the Stockmayer fluid over a range of dipole densities that covers known organic solvents, as well as that of an ellipsoidal Stockmayer-like fluid, using NVT rigid-body Monte Carlo simulations. Both fluids are examined under different electrostatic boundary conditions and ensemble sizes. While the Stockmayer model predicts that liquids of similar dipole density to acetonitrile would be ferroelectric and have a dielectric constant far higher than shown by experiment, the ellipsoidal model provides a better accounting of dielectric behavior. This result bodes well for the use of coarse-grained solvent models for large-scale simulations.

Modeling the optical behavior of complex organic media: from molecules to materials.

For the past three decades, a full understanding of the electro-optic (EO) effect in amorphous organic media has remained elusive. Calculating a bulk material property from fundamental molecular properties, intermolecular electrostatic forces, and field-induced net acentric dipolar order has proven to be very challenging. Moreover, there has been a gap between ab initio quantum-mechanical (QM) predictions of molecular properties and their experimental verification at the level of bulk materials and devices. This report unifies QM-based estimates of molecular properties with the statistical mechanical interpretation of the order in solid phases of electric-field-poled, amorphous, organic dipolar chromophore-containing materials. By combining interdependent statistical and quantum mechanical methods, bulk material EO properties are predicted. Dipolar order in bulk, amorphous phases of EO materials can be understood in terms of simple coarse-grained force field models when the dielectric properties of the media are taken into account. Parameters used in the statistical mechanical modeling are not adjusted from the QM-based values, yet the agreement with the experimentally determined electro-optic coefficient is excellent.