Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Reactions of platinum(II) complexes with dioxygen: progress toward alkane functionalization

Whereas stoichiometric activation of C-H bonds by complexes of transition metals is becoming increasingly common, selective functionalization of alkanes remains a formidable challenge in organometallic chemistry. The recent advances in catalytic alkane functionalization by transition-metal complexes are summarized in Chapter I.

The studies of the displacement of pentafluoropyridine in [(tmeda)Pt(CH_3)(NC_5F_5)][BAr^f_4] (1) with γ- tetrafluoropicoline, a very poor nucleophile, are reported in Chapter II. The ligand substitution occurs by a dissociative interchange mechanism. This result implies that dissociative loss of pentafluoropyridine is the rate-limiting step in the C-H activation reactions of 1.

Oxidation of dimethylplatinum(II) complexes (N-N)Pt(CH_3)_2 (N-N = tmeda(1), α-diimines) by dioxygen is described in Chapter III. Mechanistic studies suggest a two-step mechanism. First, a hydroperoxoplatinum(IV) complex is formed in a reaction between (N-N)Pt(CH_3)_2 and dioxygen. Next, the hydroperoxy complex reacts with a second equivalent of (N-N)Pt(CH_3)_2 to afford the final product, (N-N)Pt(OH)(OCH_3)(CH_3)_2. The hydroperoxy intermediate, (tmeda)Pt(OOH)(OCH_3)(CH_3)_2 (2), was isolated and characterized. The reactivity of 2 with several dime thylplatinum(II) complexes is reported.

The studies described in Chapter IV are directed toward the development of a platinum(II)-catalyzed oxidative alkane dehydrogenation. Stoichiometric conversion of alkanes (cyclohexane, ethane) to olefins (cyclohexene, ethylene) is achieved by C-H activation with [(N-N)Pt(CH_3)(CF_3CH_2OH)]BF_4 (1, N-N is N,N'-bis(3,5-di-t- butylphenyl)-1,4-diazabutadiene) which results in the formation of olefin hydride complexes. The first step in the C-H activation reaction is formation of a platinum(II) alkyl which undergoes β-hydrogen elimination to afford the olefin hydride complex. The cationic ethylplatinum(II) intermediate can be generated in situ by treating diethylplatinum(II) compounds with acids. Treatment of (phen)PtEt_2 with [H(OEt_2)_2]Bar^f_4 at low temperatures resulted in the formation of a mixture of [(phen)PtEt(OEt_2)]Bar^f_4 (8) and [(phen)Pt(C_2H_4)H] Bar^f_4 (7). The cationic olefin complexes are unreactive toward dioxygen or hydrogen peroxide. Since the success of the overall catalytic cycle depends on our ability to oxidize the olefin hydride complexes, a series of neutral olefin complexes of platinum(II) with monoanionic ligands (derivatives of pyrrole-2-carboxyaldehyde N-aryl imines) was prepared. Unfortunately, these are also stable to oxidation.

Ionization induced scattering of femtosecond intense laser pulses in cluster plasmas

The 45 degrees scattering of a femtosecond (60 fs) intense laser pulse with a 20 nm FWHM (the full width at half maximum) spectrum centered at 790 nm has been studied experimentally while focused in argon clusters at intensity similar to 10(16) W/cm(2). Scattering spectra under different backing pressures and laser-plasma interaction lengths were obtained, which showed spectral blueshifting, beam refraction and complex modulation. These ionization-induced effects reveal the modulation of laser pulses propagating in plasmas and the existing obstacle in laser cluster interaction at high laser intensity and high electron density.

Chelation-enforced metal-arene interactions : insights into substrate binding and catalysis by late transition metal complexes

Understanding and catalyzing chemical reactions requiring multiple electron transfers is an endeavor relevant to many outstanding challenges in the field of chemistry. To study multi-electron reactions, a terphenyl diphosphine framework was designed to support one or more metals in multiple redox states via stabilizing interactions with the central arene of the terphenyl backbone. A variety of unusual compounds and reactions and their relevance toward prominent research efforts in chemistry are the subject of this dissertation.

Chapter 2 introduces the para-terphenyl diphosphine framework and its coordination chemistry with group 10 transition metal centers. Both mononuclear and dinuclear compounds are characterized. In many cases, the metal center(s) are stabilized by the terphenyl central arene. These metal–arene interactions are characterized both statically, in the solid state, and fluxionally, in solution. As a proof-of-principle, a dinickel framework is shown to span multiple redox states, showing that multielectron chemistry can be supported by the coordinatively flexible terphenyl diphosphine.

Chapter 3 presents reactivity of the terphenyl diphosphine when bound to a metal center. Because of the dearomatizing effect of the metal center, the central arene of the ligand is susceptible to reactions that do not normally affect arenes. In particular, Ni-to-arene H-transfer and arene dihydrogenation reactions are presented. Additionally, evidence for reversibility of the Ni-to-arene H-transfer is discussed.

Chapter 4 expands beyond the chelated metal-arene interactions of the previous chapters. A dipalladium(I) terphenyl diphosphine framework is used to bind a variety of exogenous organic ligands including arenes, dienes, heteroarenes, thioethers, and anionic ligands. The compounds are structurally characterized, and many ligands exhibit unprecedented bindng modes across two metal centers. The relative binding affinities are evaluated spectroscopically, and equilibrium binding constants for the examined ligands are determined to span over 13 orders of magnitude. As an application of this framework, mild hydrogenation conditions of bound thiophene are presented.

Chapter 5 studies nickel-mediated C–O bond cleavage of aryl alkyl ethers, a transformation with emerging applications in fields such as lignin biofuels and organic methodology. Other group members have shown the mechanism of C–O bond cleavage of an aryl methyl ether incorporated into a meta-terphenyl diphosphine framework to proceed through β-H elimination of an alkoxide. First, the electronic selectivity of the model system is examined computationally and compared with catalytic systems. The lessons learned from the model system are then applied to isotopic labeling studies for catalytic aryl alkyl ether cleavage under dihydrogen. Results from selective deuteration experiments and mass spectrometry draw a clear analogy between the mechanisms of the model and catalytic systems that does not require dihydrogen for C–O bond cleavage, although dihydrogen is proposed to play a role in catalyst activation and catalytic turnover.

Appendix A presents initial efforts toward heterodinuclear complexes as models for CO dehydrogenase and Fischer Tropsch chemistry. A catechol-incorporating terphenyl diphosphine is reported, and metal complexes thereof are discussed.

Appendix B highlights some structurally characterized terphenyl diphosphine complexes that either do not thematically belong in the research chapters or proved to be difficult to reproduce. These compounds show unusual coordination modes of the terphenyl diphosphine from which other researchers may glean insights.

Fundamental Studies of Carbon Oxygen Bond Activation in Nickel Diphosphine Ether Complexes. And, Metallomacrocycles as Ligands: Synthesis and Characterization of Aluminum-Bridged Bisglyoximato Complexes of Iron and Cobalt

In order to develop better catalysts for the cleavage of aryl-X bonds fundamental studies of the mechanism and individual steps of the mechanism have been investigated in detail. As the described studies are difficult at best in catalytic systems, model systems are frequently used. To study aryl-oxygen bond activation, a terphenyl diphosphine scaffold containing an ether moiety in the central arene was designed. The first three chapters of this dissertation focus on the studies of the nickel complexes supported by this diphosphine backbone and the research efforts in regards to aryl-oxygen bond activation.

Chapter 2 outlines the synthesis of a variety of diphosphine terphenyl ether ligand scaffolds. The metallation of these scaffolds with nickel is described. The reactivity of these nickel(0) systems is also outlined. The systems were found to typically undergo a reductive cleavage of the aryl oxygen bond. The mechanism was found to be a subsequent oxidative addition, β-H elimination, reductive elimination and (or) decarbonylation.

Chapter 3 presents kinetic studies of the aryl oxygen bond in the systems outlined in Chapter 2. Using a series of nickel(0) diphosphine terphenyl ether complexes the kinetics of aryl oxygen bond activation was studied. The activation parameters of oxidative addition for the model systems were determined. Little variation was observed in the rate and activation parameters of oxidative addition with varying electronics in the model system. The cause of the lack of variation is due to the ground state and oxidative addition transition state being affected similarly. Attempts were made to extend this study to catalytic systems.

Chapter 4 investigates aryl oxygen bond activation in the presence of additives. It was found that the addition of certain metal alkyls to the nickel(0) model system lead to an increase in the rate of aryl oxygen bond activation. The addition of excess Grignard reagent led to an order of magnitude increase in the rate of aryl oxygen bond activation. Similarly the addition of AlMe3 led to a three order of magnitude rate increase. Addition of AlMe₃ at -80 °C led to the formation of an intermediate which was identified by NOESY correlations as a system in which the AlMe3 is coordinated to the ether moiety of the backbone. The rates and activation parameters of aryl oxygen bond activation in the presence of AlMe3 were investigated.

The last two chapters involve the study of metalla-macrocycles as ligands. Chapter 5 details the synthesis of a variety of glyoxime backbones and diphenol precursors and their metallation with aluminum. The coordination chemistry of iron on the aluminum scaffolds was investigated. Varying the electronics of the aluminum macrocycle was found to affect the observed electrochemistry of the iron center.

Chapter 6 extends the studies of chapter 5 to cobalt complexes. The synthesis of cobalt dialuminum glyoxime metal complexes is described. The electrochemistry of the cobalt complexes was investigated. The electrochemistry was compared to the observed electrochemistry of a zinc analog to identify the redox activity of the ligand. In the presence of acid the cobalt complexes were found to electrochemically reduce protons to dihydrogen. The electronics of the ancillary aluminum ligands were found to affect the potential of proton reduction in the cobalt complexes. These potentials were compared to other diglyoximate complexes.

Models of the Oxygen-Evolving Complex of Photosystem II

In the five chapters that follow, I delineate my efforts over the last five years to synthesize structurally and chemically relevant models of the Oxygen Evolving Complex (OEC) of Photosystem II. The OEC is nature’s only water oxidation catalyst, in that it forms the dioxygen in our atmosphere necessary for oxygenic life. Therefore understanding its structure and function is of deep fundamental interest and could provide design elements for artificial photosynthesis and manmade water oxidation catalysts. Synthetic endeavors towards OEC mimics have been an active area of research since the mid 1970s and have mutually evolved alongside biochemical and spectroscopic studies, affording ever-refined proposals for the structure of the OEC and the mechanism of water oxidation. This research has culminated in the most recent proposal: a low symmetry Mn₄CaO₅ cluster with a distorted Mn₃CaO₄ cubane bridged to a fourth, dangling Mn. To give context for how my graduate work fits into this rich history of OEC research, Chapter 1 provides a historical timeline of proposals for OEC structure, emphasizing the role that synthetic Mn and MnCa clusters have played, and ending with our Mn₃CaO₄ heterometallic cubane complexes.

In Chapter 2, the triarylbenzene ligand framework used throughout my work is introduced, and trinuclear clusters of Mn, Co, and Ni are discussed. The ligand scaffold consistently coordinates three metals in close proximity while leaving coordination sites open for further modification through ancillary ligand binding. The ligands coordinated could be varied, with a range of carboxylates and some less coordinating anions studied. These complexes’ structures, magnetic behavior, and redox properties are discussed.

Chapter 3 explores the redox chemistry of the trimanganese system more thoroughly in the presence of a fourth Mn equivalent, finding a range of oxidation states and oxide incorporation dependent on oxidant, solvent, and Mn salt. Oxidation states from Mn^II₄ to Mn^IIIMn^IV₃ were observed, with 1-4 O^2– ligands incorporated, modeling the photoactivation of the OEC. These complexes were studied by X-ray diffraction, EPR, XAS, magnetometry, and CV.

As Ca²⁺ is a necessary component of the OEC, Chapter 4 discusses synthetic strategies for making highly structurally accurate models of the OEC containing both Mn and Ca in the Mn₃CaO₄ cubane + dangling Mn geometry. Structural and electrochemical characterization of the first Mn₃CaO₄ heterometallic cubane complex— and comparison to an all-Mn Mn₄O₄ analog—suggests a role for Ca²⁺ in the OEC. Modification of the Mn₃CaO₄ system by ligand substitution affords low symmetry Mn₃CaO₄ complexes that are the most accurate models of the OEC to date.

Finally, in Chapter 5 the reactivity of the Mn₃CaO₄ cubane complexes toward O- atom transfer is discussed. The metal M strongly affects the reactivity. The mechanisms of O-atom transfer and water incorporation from and into Mn₄O₄ and Mn₄O₃ clusters, respectively, are studied through computation and ¹⁸O-labeling studies. The μ3-oxos of the Mn₄O₄ system prove fluxional, lending support for proposals of O^2– fluxionality within the OEC.

Speculation-Aware Resource Allocation for Cluster Schedulers

Real-time demand response is essential for handling the uncertainties of renewable generation. Traditionally, demand response has been focused on large industrial and commercial loads, however it is expected that a large number of small residential loads such as air conditioners, dish washers, and electric vehicles will also participate in the coming years. The electricity consumption of these smaller loads, which we call deferrable loads, can be shifted over time, and thus be used (in aggregate) to compensate for the random fluctuations in renewable generation.

In this thesis, we propose a real-time distributed deferrable load control algorithm to reduce the variance of aggregate load (load minus renewable generation) by shifting the power consumption of deferrable loads to periods with high renewable generation. The algorithm is model predictive in nature, i.e., at every time step, the algorithm minimizes the expected variance to go with updated predictions. We prove that suboptimality of this model predictive algorithm vanishes as time horizon expands in the average case analysis. Further, we prove strong concentration results on the distribution of the load variance obtained by model predictive deferrable load control. These concentration results highlight that the typical performance of model predictive deferrable load control is tightly concentrated around the average-case performance. Finally, we evaluate the algorithm via trace-based simulations.

Towards a standard methodology to evaluate internal cluster validity indices

The evaluation and comparison of internal cluster validity indices is a critical problem in the clustering area. The methodology used in most of the evaluations assumes that the clustering algorithms work correctly. We propose an alternative methodology that does not make this often false assumption. We compared 7 internal cluster validity indices with both methodologies and concluded that the results obtained with the proposed methodology are more representative of the actual capabilities of the compared indices.

I. Electron spin relaxation studies of manganese(II) complexes in acetonitrile. II. Nuclear spin relaxation studies of bromine in tetrabromomanganese(II) in acetonitrile

Magnetic resonance techniques have given us a powerful means for investigating dynamical processes in gases, liquids and solids. Dynamical effects manifest themselves in both resonance line shifts and linewidths, and, accordingly, require detailed analyses to extract desired information. The success of a magnetic resonance experiment depends critically on relaxation mechanisms to maintain thermal equilibrium between spin states. Consequently, there must be an interaction between the excited spin states and their immediate molecular environment which promote changes in spin orientation while excess magnetic energy is coupled into other degrees of freedom by non-radiative processes. This is well known as spin-lattice relaxation. Certain dynamical processes cause fluctuations in the spin state energy levels leading to spin-spin relaxation and, here again, the environment at the molecular level plays a significant role in the magnitude of interaction. Relatively few electron spin relaxation studies of solutions have been conducted and the present work is addressed toward the extension of our knowledge in this area and the retrieval of dynamical information from line shape analyses on a time scale comparable to diffusion controlled phenomena.

Specifically, the electron spin relaxation of three Mn⁺²3d⁵ complexes, Mn(CH₃CN)₆⁺², MnCl₄^-2 in acetonitrile has been studied in considerable detail. The effective spin Hamiltonian constants were carefully evaluated under a wide range of experimental conditions. Resonance widths of these Mn⁺² complexes were studied in the presence of various excess ligand ions and as a function of concentration, viscosity, temperature and frequency (X-band, ~9.5 Ԍ Hz and K-band, ~35 Ԍ Hz).

A number of interesting conclusions were drawn from these studies. For the Et₄NCl-₄^-2 system several relaxation mechanisms leading to resonance broadening were observed. One source appears to arise through spin-orbit interactions caused by modulation of the ligand field resulting from transient distortions of the complex imparted by solvent fluctuations in the immediate surroundings of the paramagnetic ion. An additional spin relaxation was assigned to the formation of ion pairs [Et₄N⁺…MnCl₄^-2] and it was possible to estimate the dissociation constant for this specie in acetonitrile.

The Bu₄NBr-MnBr₄^-2 study was considerably more interesting. As in the former case, solvent fluctuations and ion-pairing of the paramagnetic complex [Bu₄N⁺…MnBr₄^-2] provide significant relaxation for the electronic spin system. Most interesting, without doubt, is the onset of a new relaxation mechanism leading to resonance broadening which is best interpreted as chemical exchange. Thus, assuming that resonance widths were simply governed by electron spin state lifetimes, we were able to extract dynamical information from an interaction in which the initial and final states are the same

MnBr₄^-2 + Br^- = MnBr₄^-2 + Br^-.

The bimolecular rate constants were obtained at six different temperatures and their magnitudes suggested that the exchange is probably diffusion controlled with essentially a zero energy of activation. The most important source of spin relaxation in this system stems directly from dipolar interactions between the manganese 3d⁵ electrons. Moreover, the dipolar broadening is strongly frequency dependent indicating a deviation between the transverse and longitudinal relaxation times. We are led to the conclusion that the 3d⁵ spin states of ion-paired MnBr₄^-2 are significantly correlated so that dynamical processes are also entering the picture. It was possible to estimate the correlation time, ^Td, characterizing this dynamical process.

In Part II we study nuclear magnetic relaxation of bromine ions in the MnBr4-2-Bu4NBr-acetonitrile system. Essentially we monitor the 79Br and 81Br linewidths in response to the [MnBr4-2]/[Br-] ratio with the express purpose of supporting our contention that exchange is occurring between "free" bromine ions in the solvent and bromine in the first coordination sphere of the paramagnetic anion. The complexity of the system elicited a two-part study: (1) the linewidth behavior of Bu4NBr in anhydrous CH3CN in the absence of MnBr4-2 and (2) in the presence of MnBr4-2. It was concluded in study (1) that dynamical association, Bu4NBr k1= Bu4N+ + Br-, was modulating field-gradient interactions at frequencies high enough to provide an estimation of the unimolecular rate constant, k1. A comparison of the two isotopic bromine linewidth-mole fraction results led to the conclusion that quadrupole interactions provided the dominant relaxation mechanism. In study (2) the "residual" bromine linewidths for both 79Br and 81Br are clearly controlled by quadrupole interactions which appear to be modulated by very rapid dynamical processes other than molecular reorientation. We conclude that the "residual" linewidth has its origin in chemical exchange and that bromine nuclei exchange rapidly between a "free" solvated ion and the paramagnetic complex, MnBr4-2.