Replicator Dynamics and Evolutionary Stable Strategies in Heterogeneous Games

We generalise and extend the work of Iñarra and Laruelle (2011) by studying two person symmetric evolutionary games with two strategies, a heterogenous population with two possible types of individuals and incomplete information. Comparing such games with their classic homogeneous version vith complete information found in the literature, we show that for the class of anti-coordination games the only evolutionarily stable strategy vanishes. Instead, we find infinite neutrally stable strategies. We also model the evolutionary process using two different replicator dynamics setups, each with a different inheritance rule, and we show that both lead to the same results with respect to stability.

Evolution of Cooperation in the Snowdrift Game with Incomplete Information and Heterogeneous Population

Differently from previous studies of tag-based cooperation, we assume that individuals fail to recognize their own tag. Due to such incomplete information, the action taken against the opponent cannot be based on similarity, although it is still motivated by the tag displayed by the opponent. We present stability conditions for the case when individuals play unconditional cooperation, unconditional defection or conditional cooperation. We then consider the removal of one or two strategies. Results show that conditional cooperators are the most resilient agents against extinction and that the removal of unconditional cooperators may lead to the extinction of unconditional defectors.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

Emergence of Cooperation in Heterogeneous Population: A Discrete-Time Replicator Dynamics Analysis

The emergence of cooperation is analyzed in heterogeneous populations where individuals can be classified in two groups according to their phenotypic appearance. Phenotype recognition is assumed for all individuals: individuals are able to identify the type of every other individual, but fail to recognize their own type, and thus behave under partial information conditions. The interactions between individuals are described by 2 × 2 symmetric games where individuals can either cooperate or defect. The evolution of such populations is studied in the framework of evolutionary game by means of the replicator dynamics. Overlapping generations are considered, so the replicator equations are formulated in discrete-time form. The well-posedness conditions of the system are derived. Depending on the parameters of the game, a restriction may exist for the generation length. The stability analysis of the dynamical system is carried out and a detailed description of the behavior of trajectories starting from the interior of the state-space is given. We find that, provided the conditions of well-posedness are verified, the linear stability of monomorphic states in the discrete-time replicator coincides with the one of the continuous case. Specific from the discrete-time case, a relaxed restriction for the generation length is derived, for which larger time-steps can be used without compromising the well-posedness of the replicator system.

Iso-specific Ziegler-Natta polymerization of α-olefins with a single-component organoyttrium catalyst

The synthesis of a sterically tailored ligand array (M)_2((C_5H_2-2-Si(CH_3)_3-4-C(CH_3)_3)S_2i(CH_3)_2]("M_2Bp") (M = Li, 16; K, 19) is described. Transmetallation of Li_2Bp with YCl_3(THF)_3 affords exclusively the C_2 symmetric product rac-[BpY(µ_2-Cl)_2Li(THF)_2], 20. A X-ray crystal structure of 20 has been determined; triclinic, P1, a= 13.110 (8), b = 17.163 (15), c = 20.623 (14) Å, α = 104.02 (7), β = 99.38 (5), γ = 100.24 (6)° , Z = 4, R = 0.056. Transmetallation of K_2Bp with YCl_3(THF)_3 affords the halide free complex rac-BpYCl, 23. The corresponding rac-BpLaCl, 28, is prepared in an anlogous manner. In all cases the achiral meso isomer is not obtained since only for the racemic isomers are the unfavorable steric interactions between the Si(CH3)_3 groups in the narrow portion of the [Cp-M'-Cp] wedge avoided. Alkylation of 20 or 23 with LiCH(Si(CH_3)_3)_2 affords rac-BpYCH(Si(CH_3)_3)_2, 26 in good yield. Alkylation of 28 with LiCH(Si(CH_3)_3)_2 affords rac-BpLaCH(Si(CH_3)_3)_2 29. Hydrogenation of 26 cleanly affords the bridging hydride species [BpY(µ_2-H)]_2, 27, as the homochiral (R,R) and (S,S) dimeric pairs. 26 is an efficient initiator for the polymerization of ethylene to high molecular weight linear polyethylene. 27 catalyzes the polymerization of propylene (25% v/v in methylcyclohexane) and neat samples of 1-butene, 1-pentene, 1-hexene to moderately high molecular weight polymers: polypropylene (M_n = 4,200, PDI 2.32, T_m 157 °C); poly-1-butene (M_n = 8,500, PDI 3.44, T_m 105 °C); poly-1-pentene (M_n = 20,000, PDI 1.99, T_m 73 °C); poly-1-hexene (M_n = 24,000, PDI 1.75, T_m < 25 °C). ^(13)C NMR spectra at the pentad analysis level indicates that the degree of isotacticity is 99% mmmm for all polymer samples. 27 is the first single component iso-specific α-olefin polymerization catalyst. The presumed origins of the high isospecificity are presented.

The design and synthesis of a true, heterogeneous, asymmetric catalyst

This work describes the design and synthesis of a true, heterogeneous, asymmetric catalyst. The catalyst consists of a thin film that resides on a high-surface- area hydrophilic solid and is composed of a chiral, hydrophilic organometallic complex dissolved in ethylene glycol. Reactions of prochiral organic reactants take place predominantly at the ethylene glycol-bulk organic interface.

The synthesis of this new heterogeneous catalyst is accomplished in a series of designed steps. A novel, water-soluble, tetrasulfonated 2,2'-bis (diphenylphosphino)-1,1'-binaphthyl (BINAP-4S0_3Na) is synthesized by direct sulfonation of 2,2'-bis(diphenylphosphino)-1,1'-binaphthyl (BINAP). The rhodium (I) complex of BINAP-4SO_3Na is prepared and is shown to be the first homogeneous catalyst to perform asymmetric reductions of prochiral 2-acetamidoacrylic acids in neat water with enantioselectivities as high as those obtained in non-aqueous solvents. The ruthenium (II) complex, [Ru(BINAP-4SO_3Na)(benzene)Cl]Cl is also synthesized and exhibits a broader substrate specificity as well as higher enantioselectivities for the homogeneous asymmetric reduction of prochiral 2-acylamino acid precursors in water. Aquation of the ruthenium-chloro bond in water is found to be detrimental to the enantioselectivity with some substrates. Replacement of water by ethylene glycol results in the same high e.e's as those found in neat methanol. The ruthenium complex is impregnated onto a controlled pore-size glass CPG-240 by the incipient wetness technique. Anhydrous ethylene glycol is used as the immobilizing agent in this heterogeneous catalyst, and a non-polar 1:1 mixture of chloroform and cyclohexane is employed as the organic phase.

Asymmetric reduction of 2-(6'-methoxy-2'-naphthyl)acrylic acid to the non-steroidal anti-inflammatory agent, naproxen, is accomplished with this heterogeneous catalyst at a third of the rate observed in homogeneous solution with an e.e. of 96% at a reaction temperature of 3°C and 1,400 psig of hydrogen. No leaching of the ruthenium complex into the bulk organic phase is found at a detection limit of 32 ppb. Recycling of the catalyst is possible without any loss in enantioselectivity. Long-term stability of this new heterogeneous catalyst is proven by a self-assembly test. That is, under the reaction conditions, the individual components of the present catalytic system self-assemble into the supported-catalyst configuration.

The strategies outlined here for the design and synthesis of this new heterogeneous catalyst are general, and can hopefully be applied to the development of other heterogeneous, asymmetric catalysts.

Bis(thiophenolate)pyridine pincer ligands and trivalent zirconocene complexes relevant to early transition metal polymerization catalysts

Two major topics are covered: the first chapter is focused on the development of post-metallocene complexes for propylene polymerization. The second and third chapters investigate the consequences of diisobutylaluminum hydride (HAlⁱBu₂) additives in zirconocene based polymerization systems.

The synthesis, structure, and solution behavior of early metal complexes with a new tridentate LX₂ type ligand, bis(thiophenolate)pyridine ((SNS) = (2-C₆H₄S)₂-2,6-C₅H₃N) are investigated. SNS complexes of Ti, Zr, and Ta having dialkylamido coligands were synthesized and structurally characterized. The zirconium complex, (SNS)Zr(NMe₂)₂, displays C₂ symmetry in the solid state. Solid-state structures of tantalum complexes (SNS)Ta(NMe₂)₃ and (SNS)TaCl(NEt₂)₂ also display pronounced C₂ twisting of the SNS ligand. 1D and 2D NMR experiments show that (SNS)Ta(NMe₂)₃ is fluxional with rotation about the Ta N(amide) bonds occurring on the NMR timescale. The fluxional behavior of (SNS)TaCl(NEt₂)₂ in solution was also studied by variable temperature ¹H NMR. Observation of separate signals for the diastereotopic protons of the methylene unit of the diethylamide indicates that the complex remains locked on the NMR timescale in one diastereomeric conformation at temperatures below -50 °C.

Reduction of Zr(IV) metallocenium cations with sodium amalgam (NaHg) produces EPR signals assignable to Zr(III) metallocene complexes. Thus, chloro-bridged heterobinuclear ansa-zirconocenium cation [((SBI))Zr(μ-Cl)₂AlMe₂]⁺B(C₆F₅)_4¯ (SBI = rac-dimethylsilylbis(1-indenyl)), gives rise to an EPR signal assignable to the complex (SBI)Zr^III(μ-Cl)₂AlMe₂, while (SBI)Zr^III-Me and (SBI)Zr^III(-H)2AlⁱBu₂ are formed by reduction of [(SBI)Zr(μ-Me)₂AlMe₂]⁺B(C₆F₅)_4¯ and [(SBI)Zr(μ-H)₃(AlⁱBu₂)₂]⁺B(C₆F₅)4¯, respectively. These products are also formed, along with (SBI)Zr^III-ⁱBu and [(SBI)Zr^III]⁺ AlR4¯ when (SBI)ZrMe₂ reacts with HAlⁱBu₂, eliminating isobutane en route to the Zr(III) complex. Studies concerning the interconversion reactions between these and other (SBI)Zr(III) complexes and reaction mechanisms involved in their formation are also reported.

The addition of HAlⁱBu₂ to precatalyst [(SBI)Zr(µ-H)₃(AlⁱBu₂)₂]⁺ significantly slows the polymerization of propylene and changes the kinetics of polymerization from 1st to 2nd order with respect to propylene. This is likely due to competitive inhibition by HAlⁱBu₂. When the same reaction is investigated using [(ⁿBuCp)₂Zr(μ-H)₃(AlⁱBu₂)₂]⁺, hydroalumination between propylene and HAlⁱBu₂ is observed instead of propylene polymerization.

Bimetallic olefin polymerization catalysis : mechanisms and applications of proximal effects

This dissertation covers progress with bimetallic polymerization catalysts. The complexes we have designed were aimed at expanding the capabilities of homogeneous polymerization catalysts by taking advantage of multimetallic effects. Such effects were examined in group 4 and group 10 bimetallic complexes; proximity and steric repulsion were determined to be major factors in the effects observed.

Chapters 2 and 3 introduce the rigid p-terphenyl dinucleating framework utilized in most of this thesis. The permethylation of the central arene allows for the separation of syn and anti atropisomers of the terphenyl compounds. Kinetic studies were carried out to examine the isomerization of the dinucleating bis(salicylaldimine) ligand precursors. Metallation of the syn and anti bis(salicylaldimine)s using Ni(Me)2(tmeda) and excess pyridine afforded dinickel bisphenoxyiminato complexes with a methyl and a pyridyl ligand on each nickel. The syn and anti atropisomers of the dinickel complexes were structurally characterized and utilized in ethylene and ethylene/α-olefin polymerizations. Monometallic analogues were also synthesized and tested for polymerization activity. Ethylene polymerizations were performed in the presence of primary, secondary, and tertiary amines – additives that generally deactivate nickel polymerization catalysts. Inhibition of this deactivation was observed with the syn atropisomer of the bimetallic species, but not with the anti or monometallic analogues. A mechanism was proposed wherein steric repulsion of the substituents on proximal nickel centers disfavors simultaneous ligation of base to both of the metal centers. The bimetallic effect has been explored with respect to size and binding ability of the added base.

Chapter 4 presents the optimization of the bisphenoxyimine ligand synthesis and synthesis of syn and anti m-terphenyl analogues. Metallation with NiClMe(PMe₃)₂ yielded phosphine-ligated dinickel complexes, which have been structurally characterized. Ethylene/1-hexene copolymerizations in the presence of amines using Ni(COD)₂ as a phosphine scavenger showed significantly improved activity relative to the pyridine-ligated analogues. Incorporation of amino olefins in copolymerizations with ethylene was accomplished, and a mechanism was proposed based on proximal effects. Copolymerization trials with a variety of amino olefins and ethylene/1-hexene/amino olefin terpolymerizations were completed.

Early transition metal complexes based on the rigid p-terphenyl framework were designed with a variety of donor sets (Chapter 5 and Appendix B). Chapter 5 details the use of syn dizirconium di[amine bis(phenolate)] complexes for isoselective 1-hexene and propylene homopolymerizations. Ligand variation and monometallic complexes were studied to determine the origin of tacticity control. A mechanistic proposal was presented based on the symmetry at zirconium and the steric effects of the proximal metal center. Appendix B covers additional studies of bimetallic early transition metal complexes based on the p-terphenyl. Dititanium, dizirconium, and asymmetric complexes with bisphenoxyiminato ligands and derivatives thereof were targeted. Progress toward the synthesis of these complexes is described along with preliminary polymerization data. 1-hexene/diene copolymerizations and attempted polymerizations in the presence of ethers and esters with the syn dizirconium di[amine bis(phenolate)] complexes demonstrate the potential for further applications of this system in catalysis.

Appendix A includes work toward palladium catalysts for insertion polymerization of polar monomers. These complexes were based on dioxime and diimine frameworks with the intent of binding Lewis acidic metals at the oxime oxygens, at pendant phenolic donors, or at pendant aminediol moieties. The synthesis and structural characterization of a number of palladium and Lewis acid complexes is presented. Due to the instability of the desired species, efforts toward isolation of the desired complexes proved unsuccessful, though preliminary ethylene/methyl acrylate copolymerizations using in situ activation of the palladium species were attempted.

Zirconocenes as models for homogeneous Ziegler-Natta olefin polymerization catalysts

Using density functional theory, we studied the fundamental steps of olefin polymerization for zwitterionic and cationic Group IV ansa-zirconocenes and a neutral ansa- yttrocene. Complexes [H₂E(C₅H₄)₂ZrMe]ⁿ (n = 0: E = BH₂ (1), BF₂ (2), AlH₂(3); n = +: E = CH₂(4), SiH₂(5)) and H₂Si(C₅H₄)₂YMe were used as computational models. The largest differences among these three classes of compounds were the strength of olefin binding and the stability of the β-agostic alkyl intermediate towards β-hydrogen elimination. We investigated the effect of solvent on the reaction energetics for land 5. We found that in benzene the energetics became very similar except that a higher olefin insertion barrier was calculated for 1. The calculated anion affinity of [CH₃BF₃]^- was weaker towards 1 than 5. The calculated olefin binding depended primarily on the charge of the ansa linker, and the olefin insertion barrier was found to decrease steadily in the following order: [H₂C(C₅H₄)₂ZrMe]⁺ > [F₂B(C₅H₄)₂ZrMe] ≈ [H₂B(C₅H₄)₂ZrMe] > [H₂Si(C₅H₄)₂ZrMe]⁺ > [H₂Al(C₅H₄)₂ZrMe].

We prepared ansa-zirconocene dicarbonyl complexes Me₂ECp₂Zr(CO)₂ (E = Si, C), and t-butyl substituted complexes (t-BuCp)₂Zr(CO)₂, Me₂E(t-BuCp)₂Zr(CO)₂ (E = Si, C), (Me₂Si)₂(t-BuCp)₂Zr(CO)₂ as well as analogous zirconocene complexes. Both the reduction potentials and carbonyl stretching frequencies follow the same order: Me₂SiCp₂ZrCl₂> Me₂CCp₂ZrCl₂> Cp₂ZrCl₂> (Me₂Si)₂Cp₂ZrCl₂. This ordering is a result of both the donating abilities of the cyclopentadienyl substituents and the orientation of the cyclopentadiene rings. Additionally, we prepared a series of analogous cationic zirconocene complexes [LZrOCMe₃][MeB(C₆F₅)₃] (L = CP₂, Me₂SiCp₂, Me₂CCP₂, (Me₂Si)₂Cp₂) and studied the kinetics of anion dissociation. We found that the enthalpy of anion dissociation increased from 10.3 kcal•mol^-1 to 17.6 kcal•mol^-1 as exposure of the zirconium center increased.

We also prepared series of zirconocene complexes bearing 2,2-dimethyl-2-sila-4-pentenyl substituents (and methyl-substituted olefin variants). Methide abstraction with B(C₆F₅) results in reversible coordination of the tethered olefin to the cationic zirconium center. The kinetics of olefin dissociation have been examined using NMR methods, and the effects of ligand variation for unlinked, singly [SiMe₂]-linked and doubly [SiMe₂]-linked bis(cyclopentadienyl) arrangements has been compared (ΔG‡ for olefin dissociation varies from 12.8 to 15.6 kcal•mol^-1). Methide abstraction from 1,2-(SiMe₂)₂(η⁵-C₅H₃)₂Zr(CH₃)-(CH₂CMe₂CH₂CH = CH₂) results in rapid β-allyl elimination with loss of isobutene yielding the allyl cation [{1,2-(SiMe₂)₂(η⁵-C₅H₃)₂Zr(η³-CH₂CH=CH₂)]⁺.