I. Numerical solution of exact pair equations. II. Numerical solution of first-order pair equations

Part I

Numerical solutions to the S-limit equations for the helium ground state and excited triplet state and the hydride ion ground state are obtained with the second and fourth difference approximations. The results for the ground states are superior to previously reported values. The coupled equations resulting from the partial wave expansion of the exact helium atom wavefunction were solved giving accurate S-, P-, D-, F-, and G-limits. The G-limit is -2.90351 a.u. compared to the exact value of the energy of -2.90372 a.u.

Part II

The pair functions which determine the exact first-order wavefunction for the ground state of the three-electron atom are found with the matrix finite difference method. The second- and third-order energies for the (1s1s)¹S, (1s2s)³S, and (1s2s)¹S states of the two-electron atom are presented along with contour and perspective plots of the pair functions. The total energy for the three-electron atom with a nuclear charge Z is found to be E(Z) = -1.125•Z² +1.022805•Z-0.408138-0.025515•(1/Z)+O(1/Z²)a.u.

A multiple scattering problem

The present work deals with the problem of the interaction of the electromagnetic radiation with a statistical distribution of nonmagnetic dielectric particles immersed in an infinite homogeneous isotropic, non-magnetic medium. The wavelength of the incident radiation can be less, equal or greater than the linear dimension of a particle. The distance between any two particles is several wavelengths. A single particle in the absence of the others is assumed to scatter like a Rayleigh-Gans particle, i.e. interaction between the volume elements (self-interaction) is neglected. The interaction of the particles is taken into account (multiple scattering) and conditions are set up for the case of a lossless medium which guarantee that the multiple scattering contribution is more important than the self-interaction one. These conditions relate the wavelength λ and the linear dimensions of a particle a and of the region occupied by the particles D. It is found that for constant λ/a, D is proportional to λ and that |Δχ|, where Δχ is the difference in the dielectric susceptibilities between particle and medium, has to lie within a certain range.

The total scattering field is obtained as a series the several terms of which represent the corresponding multiple scattering orders. The first term is a single scattering term. The ensemble average of the total scattering intensity is then obtained as a series which does not involve terms due to products between terms of different orders. Thus the waves corresponding to different orders are independent and their Stokes parameters add.

The second and third order intensity terms are explicitly computed. The method used suggests a general approach for computing any order. It is found that in general the first order scattering intensity pattern (or phase function) peaks in the forward direction Θ = 0. The second order tends to smooth out the pattern giving a maximum in the Θ = π/2 direction and minima in the Θ = 0 , Θ = π directions. This ceases to be true if ka (where k = 2π/λ) becomes large (> 20). For large ka the forward direction is further enhanced. Similar features are expected from the higher orders even though the critical value of ka may increase with the order.

The first order polarization of the scattered wave is determined. The ensemble average of the Stokes parameters of the scattered wave is explicitly computed for the second order. A similar method can be applied for any order. It is found that the polarization of the scattered wave depends on the polarization of the incident wave. If the latter is elliptically polarized then the first order scattered wave is elliptically polarized, but in the Θ = π/2 direction is linearly polarized. If the incident wave is circularly polarized the first order scattered wave is elliptically polarized except for the directions Θ = π/2 (linearly polarized) and Θ = 0, π (circularly polarized). The handedness of the Θ = 0 wave is the same as that of the incident whereas the handedness of the Θ = π wave is opposite. If the incident wave is linearly polarized the first order scattered wave is also linearly polarized. The second order makes the total scattered wave to be elliptically polarized for any Θ no matter what the incident wave is. However, the handedness of the total scattered wave is not altered by the second order. Higher orders have similar effects as the second order.

If the medium is lossy the general approach employed for the lossless case is still valid. Only the algebra increases in complexity. It is found that the results of the lossless case are insensitive in the first order of k_imD where k_im = imaginary part of the wave vector k and D a linear characteristic dimension of the region occupied by the particles. Thus moderately extended regions and small losses make (k_imD)² ≪ 1 and the lossy character of the medium does not alter the results of the lossless case. In general the presence of the losses tends to reduce the forward scattering.

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.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Three essays on economics and information shocks

A person living in an industrialized society has almost no choice but to receive information daily with negative implications for himself or others. His attention will often be drawn to the ups and downs of economic indicators or the alleged misdeeds of leaders and organizations. Reacting to new information is central to economics, but economics typically ignores the affective aspect of the response, for example, of stress or anger. These essays present the results of considering how the affective aspect of the response can influence economic outcomes.

The first chapter presents an experiment in which individuals were presented with information about various non-profit organizations and allowed to take actions that rewarded or punished those organizations. When social interaction was introduced into this environment an asymmetry between rewarding and punishing appeared. The net effects of punishment became greater and more variable, whereas the effects of reward were unchanged. The individuals were more strongly influenced by negative social information and used that information to target unpopular organizations. These behaviors contributed to an increase in inequality among the outcomes of the organizations.

The second and third chapters present empirical studies of reactions to negative information about local economic conditions. Economic factors are among the most prevalent stressors, and stress is known to have numerous negative effects on health. These chapters document localized, transient effects of the announcement of information about large-scale job losses. News of mass layoffs and shut downs of large military bases are found to decrease birth weights and gestational ages among babies born in the affected regions. The effect magnitudes are close to those estimated in similar studies of disasters.

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.

π+ photoproduction at angles from 6° to 90° C.M. and energies from 589 to 1269 MeV

The differential cross section for the reaction γp → π⁺n was measured at 32 laboratory photon energies between 589 and 1269 MeV at the Caltech Synchrotron. At each energy, data have been obtained at typically fifteen π⁺ c.m. angles between 6° and 90°. A magnetic spectrometer was used to detect the π⁺ photo-produced in a liquid hydrogen target. Two Cherenkov counters were used to reject the background of positrons and protons. The data clearly show the presence of a pole in the production amplitude due to the one pion exchange. Moravcsik fits to the 32 angular distributions, including data from another experiment, are presented. The extrapolation of these fits to the pole gives a value for the pion-nucleon coupling constant of 14.5 which is consistent with the accepted value. The second and third pion-nucleon resonances are evident as peaks in the total cross section and as changes in the shape of the angular distributions. At the third resonance there is evidence for both a D5/2 and an F5/2 amplitude. The absence of large variations in the 0° and 180° cross sections implies that the second and third resonances are mostly produced from an initial state with helicity ± 3/2.

Topics in vortex motion

Six topics in incompressible, inviscid fluid flow involving vortex motion are presented. The stability of the unsteady flow field due to the vortex filament expanding under the influence of an axial compression is examined in the first chapter as a possible model of the vortex bursting observed in aircraft contrails. The filament with a stagnant core is found to be unstable to axisymmetric disturbances. For initial disturbances with the form of axisymmetric Kelvin waves, the filament with a uniformly rotating core is neutrally stable, but the compression causes the disturbance to undergo a rapid increase in amplitude. The time at which the increase occurs is, however, later than the observed bursting times, indicating the bursting phenomenon is not caused by this type of instability.

In the second and third chapters the stability of a steady vortex filament deformed by two-dimensional strain and shear flows, respectively, is examined. The steady deformations are in the plane of the vortex cross-section. Disturbances which deform the filament centerline into a wave which does not propagate along the filament are shown to be unstable and a method is described to calculate the wave number and corresponding growth rate of the amplified waves for a general distribution of vorticity in the vortex core.

In Chapter Four exact solutions are constructed for two-dimensional potential flow over a wing with a free ideal vortex standing over the wing. The loci of positions of the free vortex are found and the lift is calculated. It is found that the lift on the wing can be significantly increased by the free vortex.

The two-dimensional trajectories of an ideal vortex pair near an orifice are calculated in Chapter Five. Three geometries are examined, and the criteria for the vortices to travel away from the orifice are determined.

Finally, Chapter Six reproduces completely the paper, "Structure of a linear array of hollow vortices of finite cross-section," co-authored with G. R. Baker and P. G. Saffman. Free streamline theory is employed to construct an exact steady solution for a linear array of hollow, or stagnant cored vortices. If each vortex has area A and the separation is L, then there are two possible shapes if A^(1/2)/L is less than 0.38 and none if it is larger. The stability of the shapes to two-dimensional, periodic and symmetric disturbances is considered for hollow vortices. The more deformed of the two possible shapes is found to be unstable, while the less deformed shape is stable.