Palladium-Catalyzed Decarboxylative and Decarbonylative Transformations in the Synthesis of Fine and Commodity Chemicals

Decarboxylation and decarbonylation are important reactions in synthetic organic chemistry, transforming readily available carboxylic acids and their derivatives into various products through loss of carbon dioxide or carbon monoxide. In the past few decades, palladium-catalyzed decarboxylative and decarbonylative reactions experienced tremendous growth due to the excellent catalytic activity of palladium. Development of new reactions in this category for fine and commodity chemical synthesis continues to draw attention from the chemistry community.

The Stoltz laboratory has established a palladium-catalyzed enantioselective decarboxylative allylic alkylation of β-keto esters for the synthesis of α-quaternary ketones since 2005. Recently, we extended this chemistry to lactams due to the ubiquity and importance of nitrogen-containing heterocycles. A wide variety of α-quaternary and tetrasubstituted α-tertiary lactams were obtained in excellent yields and exceptional enantioselectivities using our palladium-catalyzed decarboxylative allylic alkylation chemistry. Enantioenriched α-quaternary carbonyl compounds are versatile building blocks that can be further elaborated to intercept synthetic intermediates en route to many classical natural products. Thus our chemistry enables catalytic asymmetric formal synthesis of these complex molecules.

In addition to fine chemicals, we became interested in commodity chemical synthesis using renewable feedstocks. In collaboration with the Grubbs group, we developed a palladium-catalyzed decarbonylative dehydration reaction that converts abundant and inexpensive fatty acids into value-added linear alpha olefins. The chemistry proceeds under relatively mild conditions, requires very low catalyst loading, tolerates a variety of functional groups, and is easily performed on a large scale. An additional advantage of this chemistry is that it provides access to expensive odd-numbered alpha olefins.

Finally, combining features of both projects, we applied a small-scale decarbonylative dehydration reaction to the synthesis of α-vinyl carbonyl compounds. Direct α-vinylation is challenging, and asymmetric vinylations are rare. Taking advantage of our decarbonylative dehydration chemistry, we were able to transform enantioenriched δ-oxocarboxylic acids into quaternary α-vinyl carbonyl compounds in good yields with complete retention of stereochemistry. Our explorations culminated in the catalytic enantioselective total synthesis of (–)-aspewentin B, a terpenoid natural product featuring a quaternary α-vinyl ketone. Both decarboxylative and decarbonylative chemistries found application in the late stage of the total synthesis.

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.

Polymerization of functionalized olefins with neutral group ten catalysts

This thesis describes the preparation, characterization, and application of welldefined single-component group ten salicylaldimine complexes for the polymerization of ethylene to high molecular weight materials as well as the copolymerization of ethylene and functionalized olefins. After an initial introduction to the field, Chapter 2 describes the preparation of PPh₃ complexes that contain a series of modified salicylaldimine and naphthaldimine ligands. Such complexes were activated for polymerization by the addition of cocatalysts such as Ni(COD)₂ or B(C₆F₅)₃. As the steric demand of the ligand set increased-the molecular weight, polymerization activity, and lifetime of the catalyst was observed to increase. In fact, complexes containing "bulky" ligands, such as the [^Anthr,HSal] ligand (2.5), were found to be highly-active single component complexes for the polymerization of ethylene. Model hydrido compound were prepared-allowing for a better understanding of both the mechanism of polymerization and one mode of decomposition.

Chapter 3 describes the effect which additives play on neutral Ni^II polymerization catalysts such as 2.5. The addition of excess ethers, esters, ketones, anhydrides, alcohols, and water do not deactivate the catalysts for polymerization. However, the addition of excess acid, thiols, and phosphines was observed to shut-down catalysis. Since excess phosphine was found to inhibit catalysis, "phosphine-free" complexes, such as the acetonittile complex (3.26), were prepared. The acetonitrile complex was found to be the most active neutral polymerization catalyst prepared to date.

Chapter 4 outlines the use of catalyst 2.5 and 3.26 for the preparation of linear functionalized copolymers containing alcohols, esters, anhydrides, and ethers. Copolymers can be prepared with γ-functionalized-α-olefins, functionalized norbornenes, and functionalized tricyclononenes, with up to 30 mol% comonomer incorporation.

Chapter 5 outlines the preparation of a series of Pt^II alkyl/olefin salicylaldimine complexes which serve as models for the active species in the Ni^II-catalyzed polymerization process. Understanding the nature of the M-olefin interaction as a the electronic and steric properties of the salicylaldimine ligand is varied has allowed for a number of predictions about the design of future polymerization systems.

Generalized modal identification of linear and nonlinear dynamic systems

This dissertation is concerned with the problem of determining the dynamic characteristics of complicated engineering systems and structures from the measurements made during dynamic tests or natural excitations. Particular attention is given to the identification and modeling of the behavior of structural dynamic systems in the nonlinear hysteretic response regime. Once a model for the system has been identified, it is intended to use this model to assess the condition of the system and to predict the response to future excitations.

A new identification methodology based upon a generalization of the method of modal identification for multi-degree-of-freedom dynaimcal systems subjected to base motion is developed. The situation considered herein is that in which only the base input and the response of a small number of degrees-of-freedom of the system are measured. In this method, called the generalized modal identification method, the response is separated into "modes" which are analogous to those of a linear system. Both parametric and nonparametric models can be employed to extract the unknown nature, hysteretic or nonhysteretic, of the generalized restoring force for each mode.

In this study, a simple four-term nonparametric model is used first to provide a nonhysteretic estimate of the nonlinear stiffness and energy dissipation behavior. To extract the hysteretic nature of nonlinear systems, a two-parameter distributed element model is then employed. This model exploits the results of the nonparametric identification as an initial estimate for the model parameters. This approach greatly improves the convergence of the subsequent optimization process.

The capability of the new method is verified using simulated response data from a three-degree-of-freedom system. The new method is also applied to the analysis of response data obtained from the U.S.-Japan cooperative pseudo-dynamic test of a full-scale six-story steel-frame structure.

The new system identification method described has been found to be both accurate and computationally efficient. It is believed that it will provide a useful tool for the analysis of structural response data.

Topics in linear and nonlinear dispersive waves

The various singularities and instabilities which arise in the modulation theory of dispersive wavetrains are studied. Primary interest is in the theory of nonlinear waves, but a study of associated questions in linear theory provides background information and is of independent interest.

The full modulation theory is developed in general terms. In the first approximation for slow modulations, the modulation equations are solved. In both the linear and nonlinear theories, singularities and regions of multivalued modulations are predicted. Higher order effects are considered to evaluate this first order theory. An improved approximation is presented which gives the true behavior in the singular regions. For the linear case, the end result can be interpreted as the overlap of elementary wavetrains. In the nonlinear case, it is found that a sufficiently strong nonlinearity prevents this overlap. Transition zones with a predictable structure replace the singular regions.

For linear problems, exact solutions are found by Fourier integrals and other superposition techniques. These show the true behavior when breaking modulations are predicted.

A numerical study is made for the anharmonic lattice to assess the nonlinear theory. This confirms the theoretical predictions of nonlinear group velocities, group splitting, and wavetrain instability, as well as higher order effects in the singular regions.

Numerical methods for ill-posed, linear problems

A means of assessing the effectiveness of methods used in the numerical solution of various linear ill-posed problems is outlined. Two methods: Tikhonov' s method of regularization and the quasireversibility method of Lattès and Lions are appraised from this point of view.

In the former method, Tikhonov provides a useful means for incorporating a constraint into numerical algorithms. The analysis suggests that the approach can be generalized to embody constraints other than those employed by Tikhonov. This is effected and the general "T-method" is the result.

A T-method is used on an extended version of the backwards heat equation with spatially variable coefficients. Numerical computations based upon it are performed.

The statistical method developed by Franklin is shown to have an interpretation as a T-method. This interpretation, although somewhat loose, does explain some empirical convergence properties which are difficult to pin down via a purely statistical argument.

Ring-opening metathesis of bulky norbornene monomers and the radical-mediated hydrophosphonation of olefins

This thesis discusses two major topics: the ring-opening metathesis polymerization (ROMP) of bulky monomers and the radical-mediated hydrophosphonation of olefins. The research into the ROMP of bulky monomers is further divided into three chapters: wedge-shaped monomers, the alternating copolymerization of 1-methyloxanorbornene derivatives with cyclooctene, and the kinetic resolution polymerization of 1-methyloxanorbornene derivatives. The wedge-shaped monomers can be polymerized into diblock copolymers that possess photonic crystal properties. The alternating copolymerization of 1-methyloxanorbornene derivatives with cyclooctene is performed with > 90% alternation via two different routes: typical alternating copolymerization and a sequence editing approach. The kinetic resolution polymerization of these same 1-methyloxanorbornene monomers achieves only modest selectivity (S=4), but there is evidence that the growing polymer chain forms a helix that influences the selectivity of the resolution. The last topic is the radical-mediated hydrophosphonation of olefins. This synthetic method provides access to Wittig reagents that are capable of highly cis-selective olefinations of aldehydes.

Non-linear dispersive waves with a small dissipation

The general theory of Whitham for slowly-varying non-linear wavetrains is extended to the case where some of the defining partial differential equations cannot be put into conservation form. Typical examples are considered in plasma dynamics and water waves in which the lack of a conservation form is due to dissipation; an additional non-conservative element, the presence of an external force, is treated for the plasma dynamics example. Certain numerical solutions of the water waves problem (the Korteweg-de Vries equation with dissipation) are considered and compared with perturbation expansions about the linearized solution; it is found that the first correction term in the perturbation expansion is an excellent qualitative indicator of the deviation of the dissipative decay rate from linearity.

A method for deriving necessary and sufficient conditions for the existence of a general uniform wavetrain solution is presented and illustrated in the plasma dynamics problem. Peaking of the plasma wave is demonstrated, and it is shown that the necessary and sufficient existence conditions are essentially equivalent to the statement that no wave may have an amplitude larger than the peaked wave.

A new type of fully non-linear stability criterion is developed for the plasma uniform wavetrain. It is shown explicitly that this wavetrain is stable in the near-linear limit. The nature of this new type of stability is discussed.

Steady shock solutions are also considered. By a quite general method, it is demonstrated that the plasma equations studied here have no steady shock solutions whatsoever. A special type of steady shock is proposed, in which a uniform wavetrain joins across a jump discontinuity to a constant state. Such shocks may indeed exist for the Korteweg-de Vries equation, but are barred from the plasma problem because entropy would decrease across the shock front.

Finally, a way of including the Landau damping mechanism in the plasma equations is given. It involves putting in a dissipation term of convolution integral form, and parallels a similar approach of Whitham in water wave theory. An important application of this would be towards resolving long-standing difficulties about the "collisionless" shock.

Some approximate solutions of dynamic problems in the linear theory of thin elastic shells

Some aspects of wave propagation in thin elastic shells are considered. The governing equations are derived by a method which makes their relationship to the exact equations of linear elasticity quite clear. Finite wave propagation speeds are ensured by the inclusion of the appropriate physical effects.

The problem of a constant pressure front moving with constant velocity along a semi-infinite circular cylindrical shell is studied. The behavior of the solution immediately under the leading wave is found, as well as the short time solution behind the characteristic wavefronts. The main long time disturbance is found to travel with the velocity of very long longitudinal waves in a bar and an expression for this part of the solution is given.

When a constant moment is applied to the lip of an open spherical shell, there is an interesting effect due to the focusing of the waves. This phenomenon is studied and an expression is derived for the wavefront behavior for the first passage of the leading wave and its first reflection.

For the two problems mentioned, the method used involves reducing the governing partial differential equations to ordinary differential equations by means of a Laplace transform in time. The information sought is then extracted by doing the appropriate asymptotic expansion with the Laplace variable as parameter.

A singularly perturbed linear two-point boundary-value problem

We consider the following singularly perturbed linear two-point boundary-value problem:

Ly(x) ≡ Ω(ε)D_xy(x) - A(x,ε)y(x) = f(x,ε) 0≤x≤1 (1a)

By ≡ L(ε)y(0) + R(ε)y(1) = g(ε) ε → 0^+ (1b)

Here Ω(ε) is a diagonal matrix whose first m diagonal elements are 1 and last m elements are ε. Aside from reasonable continuity conditions placed on A, L, R, f, g, we assume the lower right mxm principle submatrix of A has no eigenvalues whose real part is zero. Under these assumptions a constructive technique is used to derive sufficient conditions for the existence of a unique solution of (1). These sufficient conditions are used to define when (1) is a regular problem. It is then shown that as ε → 0^+ the solution of a regular problem exists and converges on every closed subinterval of (0,1) to a solution of the reduced problem. The reduced problem consists of the differential equation obtained by formally setting ε equal to zero in (1a) and initial conditions obtained from the boundary conditions (1b). Several examples of regular problems are also considered.

A similar technique is used to derive the properties of the solution of a particular difference scheme used to approximate (1). Under restrictions on the boundary conditions (1b) it is shown that for the stepsize much larger than ε the solution of the difference scheme, when applied to a regular problem, accurately represents the solution of the reduced problem.

Furthermore, the existence of a similarity transformation which block diagonalizes a matrix is presented as well as exponential bounds on certain fundamental solution matrices associated with the problem (1).

Diagonal forms, linear algebraic methods and Ramsey-type problems

This thesis focuses mainly on linear algebraic aspects of combinatorics. Let N_t(H) be an incidence matrix with edges versus all subhypergraphs of a complete hypergraph that are isomorphic to H. Richard M. Wilson and the author find the general formula for the Smith normal form or diagonal form of N_t(H) for all simple graphs H and for a very general class of t-uniform hypergraphs H.

As a continuation, the author determines the formula for diagonal forms of integer matrices obtained from other combinatorial structures, including incidence matrices for subgraphs of a complete bipartite graph and inclusion matrices for multisets.

One major application of diagonal forms is in zero-sum Ramsey theory. For instance, Caro's results in zero-sum Ramsey numbers for graphs and Caro and Yuster's results in zero-sum bipartite Ramsey numbers can be reproduced. These results are further generalized to t-uniform hypergraphs. Other applications include signed bipartite graph designs.

Research results on some other problems are also included in this thesis, such as a Ramsey-type problem on equipartitions, Hartman's conjecture on large sets of designs and a matroid theory problem proposed by Welsh.

Interplay of martensitic phase transformation and plastic slip in polycrystals

Inspired by key experimental and analytical results regarding Shape Memory Alloys (SMAs), we propose a modelling framework to explore the interplay between martensitic phase transformations and plastic slip in polycrystalline materials, with an eye towards computational efficiency. The resulting framework uses a convexified potential for the internal energy density to capture the stored energy associated with transformation at the meso-scale, and introduces kinetic potentials to govern the evolution of transformation and plastic slip. The framework is novel in the way it treats plasticity on par with transformation.

We implement the framework in the setting of anti-plane shear, using a staggered implicit/explict update: we first use a Fast-Fourier Transform (FFT) solver based on an Augmented Lagrangian formulation to implicitly solve for the full-field displacements of a simulated polycrystal, then explicitly update the volume fraction of martensite and plastic slip using their respective stick-slip type kinetic laws. We observe that, even in this simple setting with an idealized material comprising four martensitic variants and four slip systems, the model recovers a rich variety of SMA type behaviors. We use this model to gain insight into the isothermal behavior of stress-stabilized martensite, looking at the effects of the relative plastic yield strength, the memory of deformation history under non-proportional loading, and several others.

We extend the framework to the generalized 3-D setting, for which the convexified potential is a lower bound on the actual internal energy, and show that the fully implicit discrete time formulation of the framework is governed by a variational principle for mechanical equilibrium. We further propose an extension of the method to finite deformations via an exponential mapping. We implement the generalized framework using an existing Optimal Transport Mesh-free (OTM) solver. We then model the $\alpha$--$\gamma$ and $\alpha$--$\varepsilon$ transformations in pure iron, with an initial attempt in the latter to account for twinning in the parent phase. We demonstrate the scalability of the framework to large scale computing by simulating Taylor impact experiments, observing nearly linear (ideal) speed-up through 256 MPI tasks. Finally, we present preliminary results of a simulated Split-Hopkinson Pressure Bar (SHPB) experiment using the $\alpha$--$\varepsilon$ model.

Topics in randomized numerical linear algebra

This thesis studies three classes of randomized numerical linear algebra algorithms, namely: (i) randomized matrix sparsification algorithms, (ii) low-rank approximation algorithms that use randomized unitary transformations, and (iii) low-rank approximation algorithms for positive-semidefinite (PSD) matrices.

Randomized matrix sparsification algorithms set randomly chosen entries of the input matrix to zero. When the approximant is substituted for the original matrix in computations, its sparsity allows one to employ faster sparsity-exploiting algorithms. This thesis contributes bounds on the approximation error of nonuniform randomized sparsification schemes, measured in the spectral norm and two NP-hard norms that are of interest in computational graph theory and subset selection applications.

Low-rank approximations based on randomized unitary transformations have several desirable properties: they have low communication costs, are amenable to parallel implementation, and exploit the existence of fast transform algorithms. This thesis investigates the tradeoff between the accuracy and cost of generating such approximations. State-of-the-art spectral and Frobenius-norm error bounds are provided.

The last class of algorithms considered are SPSD "sketching" algorithms. Such sketches can be computed faster than approximations based on projecting onto mixtures of the columns of the matrix. The performance of several such sketching schemes is empirically evaluated using a suite of canonical matrices drawn from machine learning and data analysis applications, and a framework is developed for establishing theoretical error bounds.

In addition to studying these algorithms, this thesis extends the Matrix Laplace Transform framework to derive Chernoff and Bernstein inequalities that apply to all the eigenvalues of certain classes of random matrices. These inequalities are used to investigate the behavior of the singular values of a matrix under random sampling, and to derive convergence rates for each individual eigenvalue of a sample covariance matrix.

Projections in a normed linear space and a generalization of the pseudo-inverse

The concept of a "projection function" in a finite-dimensional real or complex normed linear space H (the function P_M which carries every element into the closest element of a given subspace M) is set forth and examined.

If dim M = dim H - 1, then P_M is linear. If P_N is linear for all k-dimensional subspaces N, where 1 ≤ k < dim M, then P_M is linear.

The projective bound Q, defined to be the supremum of the operator norm of P_M for all subspaces, is in the range 1 ≤ Q < 2, and these limits are the best possible. For norms with Q = 1, P_M is always linear, and a characterization of those norms is given.

If H also has an inner product (defined independently of the norm), so that a dual norm can be defined, then when P_M is linear its adjoint P_M^H is the projection on (kernel P_M)^⊥ by the dual norm. The projective bounds of a norm and its dual are equal.

The notion of a pseudo-inverse F⁺ of a linear transformation F is extended to non-Euclidean norms. The distance from F to the set of linear transformations G of lower rank (in the sense of the operator norm ∥F - G∥) is c/∥F⁺∥, where c = 1 if the range of F fills its space, and 1 ≤ c < Q otherwise. The norms on both domain and range spaces have Q = 1 if and only if (F⁺)⁺ = F for every F. This condition is also sufficient to prove that we have (F⁺)^H = (F^H)⁺, where the latter pseudo-inverse is taken using dual norms.

In all results, the real and complex cases are handled in a completely parallel fashion.

Selectivity, activity, and stability of ruthenium-carbene based olefin metathesis initiators

The olefin metathesis reaction has found many applications in polymer synthesis and more recently in organic synthesis. The use of single component late metal olefin metathesis catalysts has expanded the scope of the reaction to many new applications and has allowed for detailed study of the catalytic species.

The metathesis of terminal olefins of different steric bulk, different geometry as well as electronically different para-substituted styrenes was studied with the ruthenium based metathesis initiators, trans-(PCy₃)₂Cl₂Ru=CHR, of different carbene substituents. Increasing olefin bulk was found to slow the rate of reaction and trans internal olefins were found to be slower to react than cis internal olefins. The kinetic product of a11 reactions was found to be the alkylidene, rather than the methylidene, suggesting the intermediacy of a 2,4-metallacycle. The observed effects were used to explain the mechanism of ring opening cross metathesis and acyclic diene metathesis polymerization. No linear electronic effects were observed.

In studying the different carbene ligands, a series of ester-carbene complexes was synthesized. These complexes were found to be highly active for the metathesis of olefinic substrates, including acrylates and trisubstituted olefins. In addition, the estercarbene moiety is thermodynamically high in energy. As a result, these complexes react to ring-open cyclohexene by metathesis to alleviate the thermodynamic strain of the ester-carbene ligand. However, ester-carbene complexes were found to be thermolytically unstable in solution.

Thermolytic decomposition pathways were studied for several ruthenium-carbene based olefin metathesis catalysts. Substituted carbenes were found to decompose through bimolecular pathways while the unsubstituted carbene (the methylidene) was found to decompose unimolecularly. The stability of several derivatives of the bis-phosphine ruthenium based catalysts was studied for its implications to ring-closing metathesis. The reasons for the activity and stability of the different ruthenium-based catalysts is discussed.

The difference in catalyst activity and initiation is discussed for the bis-phosphine based and mixed N-heterocyclic carbene/phosphine based ruthenium olefin metathesis catalysts. The mixed ligand catalysts initiate far slower than the bis-phosphine catalysts but are far more metathesis active. A scheme is proposed to explain the difference in reactivity between the two types of catalysts.