I. Exact solution of some nonlinear evolution equations. II. The similarity solution for the Korteweg-De Vries equation and the related Painleve Transcendent

In Part I, a method for finding solutions of certain diffusive dispersive nonlinear evolution equations is introduced. The method consists of a straightforward iteration procedure, applied to the equation as it stands (in most cases), which can be carried out to all terms, followed by a summation of the resulting infinite series, sometimes directly and other times in terms of traces of inverses of operators in an appropriate space.

We first illustrate our method with Burgers' and Thomas' equations, and show how it quickly leads to the Cole-Hopft transformation, which is known to linearize these equations.

We also apply this method to the Korteweg and de Vries, nonlinear (cubic) Schrödinger, Sine-Gordon, modified KdV and Boussinesq equations. In all these cases the multisoliton solutions are easily obtained and new expressions for some of them follow. More generally we show that the Marcenko integral equations, together with the inverse problem that originates them, follow naturally from our expressions.

Only solutions that are small in some sense (i.e., they tend to zero as the independent variable goes to ∞) are covered by our methods. However, by the study of the effect of writing the initial iterate u_1 = u_(1)(x,t) as a sum u_1 = ^∼/u_1 + ^≈/u_1 when we know the solution which results if u_1 = ^∼/u_1, we are led to expressions that describe the interaction of two arbitrary solutions, only one of which is small. This should not be confused with Backlund transformations and is more in the direction of performing the inverse scattering over an arbitrary “base” solution. Thus we are able to write expressions for the interaction of a cnoidal wave with a multisoliton in the case of the KdV equation; these expressions are somewhat different from the ones obtained by Wahlquist (1976). Similarly, we find multi-dark-pulse solutions and solutions describing the interaction of envelope-solitons with a uniform wave train in the case of the Schrodinger equation.

Other equations tractable by our method are presented. These include the following equations: Self-induced transparency, reduced Maxwell-Bloch, and a two-dimensional nonlinear Schrodinger. Higher order and matrix-valued equations with nonscalar dispersion functions are also presented.

In Part II, the second Painleve transcendent is treated in conjunction with the similarity solutions of the Korteweg-de Vries equat ion and the modified Korteweg-de Vries equation.

Numerical solution of parabolic equations by the box scheme

The box scheme proposed by H. B. Keller is a numerical method for solving parabolic partial differential equations. We give a convergence proof of this scheme for the heat equation, for a linear parabolic system, and for a class of nonlinear parabolic equations. Von Neumann stability is shown to hold for the box scheme combined with the method of fractional steps to solve the two-dimensional heat equation. Computations were performed on Burgers' equation with three different initial conditions, and Richardson extrapolation is shown to be effective.

Reactivity of titanocene methylidene with metal halides, alkene sulfides, and alkene oxides

Titanocene metallacyclobutanes show a wide variety of reactivites with organic and inorganic reagents. Their reactions include methylene transfer to organic carbonyls, formation of enolates, electron transfer from activated alkyl chlorides, olefin metathesis, ring opening polymerization. Recently, preparations of heterobinuclear µ-methylene complexes were reported. In this thesis, mechanistic, synthetic, and structural studies of the heterobinuclear µ-methylene complexes will be described. Also, the reaction of titanocene methylidene trimethylphosphine complex with alkene sulfide and styrene sulfide will be presented.

Heterobinuclear µ-methylene-µ-methyl complexes C_(p2)Ti(µ-CH_2)( µ-CH_3)M(1,5-COD) have been prepared (M = Rh, Ir). X-ray crystallography showed that the methyl group of the complex was bonded to the rhodium and bridges to the titanium through an agostic bond. The ^(1)H,^(13)CNMR, IR spectra along with partial deuteration studies supported the structure in both solution and solid state. Activation of the agostic bond is demonstrated by the equilibration of the µ-CH_3 and µ-CH_2 groups. A nonlinear Arrhenius plot, an unusually large kinetic isotope effect (24(5)), and a large negative activation entropy (-64(3)eu) can be explained by the quantum-mechanical tunneling. Calculated rate constants with Bell-type barrier fitted well with the observed one. This equilibration was best explained by a 4e-4c mechanism (or σ bond metathesis) with the character of quantum-mechanical tunneling.

Heterobinuclear µ-methylene-µ-phenyl complexes were synthesized. Structural study of C_(p2)Ti(µ-CH_(2))(µ-p-Me_(2)NC_(6)H_(4))Rh(l,5-COD) showed that the two metal atoms are bridged by the methylene carbon and the ipso carbon of the p-N,N-dimethylarninophenyl group. The analogous structure of C_(p2))Ti(µ-CH_(2))(µ-o-MeOC_(6)H_(4))Rh(1,5-COD) has been verified by the differential NOE. The aromaticity of the phenyl group observed by ^(1)H NMR, was confirmed by the comparison of the C-C bond lengths in the crystallographic structure. The unusual downfield shifts of the ipso carbon in the ^(13)C NMR are assumed to be an indication of the interaction between the ipso carbon and electron-deficient titanium.

Titanium-platinum heterobinuclear µ-methylene complexes C_(p2)Ti(µ-CH_(2))(µ -X)Pt(Me)(PM_(2)Ph) have been prepared (X= Cl, Me). Structural studies indicate the following:(1) the Ti-CH2 bond possesses residual double bond character, (2) there is a dative Pt→Ti interaction which may be regarded as a π back donation from the platinum atom to the 'Ti=CH_(2)'' group, and (3) the µ-CH_3 group is bound to the titanium atom through a three-center, two-electron agostic bond.

Titanocene (η^(2)-thioformaldehyde)•PMe_3 was prepared from C_(p2)Ti=CH_(2)•PMe_3 and sulfur-containing organic compounds (e.g. alkene sulfide, triphenylphosphine sulfide) including elemental sulfur. Mechanistic studies utilizing trans-styrene sulfide-d_1 suggested the stepwise reaction to explain equimolar mixture of trans- and cis-styrene-d_1 as by-products. The product reacted with methyl iodide to produce cationic titanocene (η_(2)-thiomethoxymethyl) complex. Complexes having less coordinating anion like BF_4 or BPh_4 could be obtained through metathesis. Together with structural analyses, the further reactivities of the complexes have been explored.

The complex C_(p2)TiOCH_(2)CH(Ph)CH_2 was prepared from the compound C_(p2)Ti=CH_(2)-PMe_3 and styrene oxide. The product was characterized with ^(1)H-^(1)H correlated 2-dimensional NMR, selective decoupling of ^(1)H NMR, and differential NOE. Stereospecificity of deuterium in the product was lost when trans-styrene oxide-d_1 was allowed to react. Relative rates of the reaction were measured with varying substituents on the phenyl ring. Better linearity (r = -0.98, p^(+) = -0.79) was observed with σ_(p)^(+)than σ(r = -0.87, p = -1.26). The small magnitude of p^+ value and stereospecificity loss during the formation of product were best explained by the generation of biradicals, but partial generation of charge cannot be excluded. Carbonylation of the product followed by exposure to iodine yields the corresponding β-phenyl γ-lactone.

Structure, chemistry, and energetics of organic and inorganic adsorbates on Ga-rich GaAs and GaP(001) surfaces

The work described in this dissertation includes fundamental investigations into three surface processes, namely inorganic film growth, water-induced oxidation, and organic functionalization/passivation, on the GaP and GaAs(001) surfaces. The techniques used to carry out this work include scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Atomic structure, electronic structure, reaction mechanisms, and energetics related to these surface processes are discussed at atomic or molecular levels.

First, we investigate epitaxial Zn₃P₂ films grown on the Ga-rich GaAs(001)(6×6) surface. The film growth mechanism, electronic properties, and atomic structure of the Zn₃P₂/GaAs(001) system are discussed based on experimental and theoretical observations. We discover that a P-rich amorphous layer covers the crystalline Zn3P2 film during and after growth. We also propose more accurate picture of the GaP interfacial layer between Zn₃P₂ and GaAs, based on the atomic structure, chemical bonding, band diagram, and P-replacement energetics, than was previously anticipated.

Second, DFT calculations are carried out in order to understand water-induced oxidation mechanisms on the Ga-rich GaP(001)(2×4) surface. Structural and energetic information of every step in the gaseous water-induced GaP oxidation reactions are elucidated at the atomic level in great detail. We explore all reasonable ground states involved in most of the possible adsorption and decomposition pathways. We also investigate structures and energies of the transition states in the first hydrogen dissociation of a water molecule on the (2×4) surface.

Finally, adsorption structures and thermal decomposition reactions of 1-propanethiol on the Ga-rich GaP(001)(2×4) surface are investigated using high resolution STM, XPS, and DFT simulations. We elucidate adsorption locations and their associated atomic structures of a single 1-propanethiol molecule on the (2×4) surface as a function of annealing temperature. DFT calculations are carried out to optimize ground state structures and search transition states. XPS is used to investigate variations of the chemical bonding nature and coverage of the adsorbate species.

Intramolecular conflict : conformation and self-assembly of architecturally complex macromolecules in solution

The solution behavior of linear polymer chains is well understood, having been the subject of intense study throughout the previous century. As plastics have become ubiquitous in everyday life, polymer science has grown into a major field of study. The conformation of a polymer in solution depends on the molecular architecture and its interactions with the surroundings. Developments in synthetic techniques have led to the creation of precision-tailored polymeric materials with varied topologies and functionalities. In order to design materials with the desired properties, it is imperative to understand the relationships between polymer architecture and their conformation and behavior. To meet that need, this thesis investigates the conformation and self-assembly of three architecturally complex macromolecular systems with rich and varied behaviors driven by the resolution of intramolecular conflicts. First we describe the development of a robust and facile synthetic approach to reproducible bottlebrush polymers (Chapter 2). The method was used to produce homologous series of bottlebrush polymers with polynorbornene backbones, which revealed the effect of side-chain and backbone length on the overall conformation in both good and theta solvent conditions (Chapter 3). The side-chain conformation was obtained from a series of SANS experiments and determined to be indistinguishable from the behavior of free linear polymer chains. Using deuterium-labeled bottlebrushes, we were able for the first time to directly observe the backbone conformation of a bottlebrush polymer which showed self-avoiding walk behavior. Secondly, a series of SANS experiments was conducted on a homologous series of Side Group Liquid Crystalline Polymers (SGLCPs) in a perdeuterated small molecule liquid crystal (5CB). Monodomain, aligned, dilute samples of SGLCP-b-PS block copolymers were seen to self-assemble into complex micellar structures with mutually orthogonally oriented anisotropies at different length scales (Chapter 4). Finally, we present the results from the first scattering experiments on a set of fuel-soluble, associating telechelic polymers. We observed the formation of supramolecular aggregates in dilute (≤0.5wt%) solutions of telechelic polymers and determined that the choice of solvent has a significant effect on the strength of association and the size of the supramolecules (Chapter 5). A method was developed for the direct estimation of supramolecular aggregation number from SANS data. The insight into structure-property relationships obtained from this work will enable the more targeted development of these molecular architectures for their respective applications.

On the solution of first excursion problems by simulation with applications to probabilistic seismic performance assessment

In a probabilistic assessment of the performance of structures subjected to uncertain environmental loads such as earthquakes, an important problem is to determine the probability that the structural response exceeds some specified limits within a given duration of interest. This problem is known as the first excursion problem, and it has been a challenging problem in the theory of stochastic dynamics and reliability analysis. In spite of the enormous amount of attention the problem has received, there is no procedure available for its general solution, especially for engineering problems of interest where the complexity of the system is large and the failure probability is small.

The application of simulation methods to solving the first excursion problem is investigated in this dissertation, with the objective of assessing the probabilistic performance of structures subjected to uncertain earthquake excitations modeled by stochastic processes. From a simulation perspective, the major difficulty in the first excursion problem comes from the large number of uncertain parameters often encountered in the stochastic description of the excitation. Existing simulation tools are examined, with special regard to their applicability in problems with a large number of uncertain parameters. Two efficient simulation methods are developed to solve the first excursion problem. The first method is developed specifically for linear dynamical systems, and it is found to be extremely efficient compared to existing techniques. The second method is more robust to the type of problem, and it is applicable to general dynamical systems. It is efficient for estimating small failure probabilities because the computational effort grows at a much slower rate with decreasing failure probability than standard Monte Carlo simulation. The simulation methods are applied to assess the probabilistic performance of structures subjected to uncertain earthquake excitation. Failure analysis is also carried out using the samples generated during simulation, which provide insight into the probable scenarios that will occur given that a structure fails.

Studies of ambient organic and inorganic aerosol in Southern California

The negative impacts of ambient aerosol particles, or particulate matter (PM), on human health and climate are well recognized. However, owing to the complexity of aerosol particle formation and chemical evolution, emissions control strategies remain difficult to develop in a cost effective manner. In this work, three studies are presented to address several key issues currently stymieing California's efforts to continue improving its air quality.

Gas-phase organic mass (GPOM) and CO emission factors are used in conjunction with measured enhancements in oxygenated organic aerosol (OOA) relative to CO to quantify the significant lack of closure between expected and observed organic aerosol concentrations attributable to fossil-fuel emissions. Two possible conclusions emerge from the analysis to yield consistency with the ambient organic data: (1) vehicular emissions are not a dominant source of anthropogenic fossil SOA in the Los Angeles Basin, or (2) the ambient SOA mass yields used to determine the SOA formation potential of vehicular emissions are substantially higher than those derived from laboratory chamber studies. Additional laboratory chamber studies confirm that, owing to vapor-phase wall loss, the SOA mass yields currently used in virtually all 3D chemical transport models are biased low by as much as a factor of 4. Furthermore, predictions from the Statistical Oxidation Model suggest that this bias could be as high as a factor of 8 if the influence of the chamber walls could be removed entirely.

Once vapor-phase wall loss has been accounted for in a new suite of laboratory chamber experiments, the SOA parameterizations within atmospheric chemical transport models should also be updated. To address the numerical challenges of implementing the next generation of SOA models in atmospheric chemical transport models, a novel mathematical framework, termed the Moment Method, is designed and presented. Assessment of the Moment Method strengths and weaknesses provide valuable insight that can guide future development of SOA modules for atmospheric CTMs.

Finally, regional inorganic aerosol formation and evolution is investigated via detailed comparison of predictions from the Community Multiscale Air Quality (CMAQ version 4.7.1) model against a suite of airborne and ground-based meteorological measurements, gas- and aerosol-phase inorganic measurements, and black carbon (BC) measurements over Southern California during the CalNex field campaign in May/June 2010. Results suggests that continuing to target sulfur emissions with the hopes of reducing ambient PM concentrations may not the most effective strategy for Southern California. Instead, targeting dairy emissions is likely to be an effective strategy for substantially reducing ammonium nitrate concentrations in the eastern part of the Los Angeles Basin.

Interaction of aqueous inorganic orthophosphate and phosphate rock

An experimental study was made of the interaction of phosphate rock and aqueous inorganic orthophosphate, calcium, and hydroxyl ions. A model of the reaction was developed by observing electron diffraction patterns in conjunction with concentration changes of chemical components. The model was applied in explaining the performance of batch suspensions of powdered phosphate rock and packed columns of granular phosphate rock. In both cases the reaction consisted initially of a rapid nucleation phase that occurred in a time period of minutes. In the batch system the calcium phosphate nuclei then ripened into larger micro-crystals of hydroxyapatite, which eventually became indistinguishable from the original phosphate rock surface. During column operation the high supersaturation ratio that existed after the rapid nucleation phase resulted in a layer of small nuclei that covered a slowly growing hydroxyapatite crystal.

The column steady-state rate constant was found to increase with increasing temperature, pH, and fluoride concentration, and to decrease with increasing concentrations of magnesium sulfate, ammonium chloride, and bicarbonate ion.

An engineering feasibility study indicated that, based on economic considerations, nucleation of apatite on phosphate rock ore has a potential use as a wastewater phosphate removal treatment process.

Oxidation of Volatile Organic Compounds in Aqueous Solution and at the Air-water Interface of Aqueous Microdroplets

Isoprene (ISO),the most abundant non-methane VOC, is the major contributor to secondary organic aerosols (SOA) formation. The mechanisms involved in such transformation, however, are not fully understood. Current mechanisms, which are based on the oxidation of ISO in the gas-phase, underestimate SOA yields. The heightened awareness that ISO is only partially processed in the gas-phase has turned attention to heterogeneous processes as alternative pathways toward SOA.

During my research project, I investigated the photochemical oxidation of isoprene in bulk water. Below, I will report on the λ > 305 nm photolysis of H₂O₂ in dilute ISO solutions. This process yields C₁₀H₁₅OH species as primary products, whose formation both requires and is inhibited by O₂. Several isomers of C₁₀H₁₅OH were resolved by reverse-phase high-performance liquid chromatography and detected as MH⁺ (m/z = 153) and MH⁺-18 (m/z = 135) signals by electrospray ionization mass spectrometry. This finding is consistent with the addition of ·OH to ISO, followed by HO-ISO· reactions with ISO (in competition with O₂) leading to second generation HO(ISO)₂· radicals that terminate as C₁₀H₁₅OH via β-H abstraction by O₂.

It is not generally realized that chemistry on the surface of water cannot be deduced, extrapolated or translated to those in bulk gas and liquid phases. The water density drops a thousand-fold within a few Angstroms through the gas-liquid interfacial region and therefore hydrophobic VOCs such as ISO will likely remain in these relatively 'dry' interfacial water layers rather than proceed into bulk water. In previous experiments from our laboratory, it was found that gas-phase olefins can be protonated on the surface of pH < 4 water. This phenomenon increases the residence time of gases at the interface, an event that makes them increasingly susceptible to interaction with gaseous atmospheric oxidants such as ozone and hydroxyl radicals.

In order to test this hypothesis, I carried out experiments in which ISO(g) collides with the surface of aqueous microdroplets of various compositions. Herein I report that ISO(g) is oxidized into soluble species via Fenton chemistry on the surface of aqueous Fe(II)Cl₂ solutions simultaneously exposed to H₂O₂(g). Monomer and oligomeric species (ISO)1-8H⁺ were detected via online electrospray ionization mass spectrometry (ESI-MS) on the surface of pH ~ 2 water, and were then oxidized into a suite of products whose combined yields exceed ~ 5% of (ISO)1-8H⁺. MS/MS analysis revealed that products mainly consisted of alcohols, ketones, epoxides and acids. Our experiments demonstrated that olefins in ambient air may be oxidized upon impact on the surface of Fe-containing aqueous acidic media, such as those of typical to tropospheric aerosols.

Related experiments involving the reaction of ISO(g) with ·OH radicals from the photolysis of dissolved H₂O₂ were also carried out to test the surface oxidation of ISO(g) by photolyzing H₂O₂(aq) at 266 nm at various pH. The products were analyzed via online electrospray ionization mass spectrometry. Similar to our Fenton experiments, we detected (ISO)1-7H⁺ at pH < 4, and new m/z⁺ = 271 and m/z^- = 76 products at pH > 5.

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.

Solar flare particle propagation--comparison of a new analytic solution with spacecraft measurements

A new analytic solution has been obtained to the complete Fokker-Planck equation for solar flare particle propagation including the effects of convection, energy-change, corotation, and diffusion with ĸ_r = constant and ĸ_Ɵ ∝ r². It is assumed that the particles are injected impulsively at a single point in space, and that a boundary exists beyond which the particles are free to escape. Several solar flare particle events have been observed with the Caltech Solar and Galactic Cosmic Ray Experiment aboard OGO-6. Detailed comparisons of the predictions of the new solution with these observations of 1-70 MeV protons show that the model adequately describes both the rise and decay times, indicating that ĸ_r = constant is a better description of conditions inside 1 AU than is ĸ_r ∝ r. With an outer boundary at 2.7 AU, a solar wind velocity of 400 km/sec, and a radial diffusion coefficient ĸ_r ≈ 2-8 x 10²⁰ cm²/sec, the model gives reasonable fits to the time-profile of 1-10 MeV protons from "classical" flare-associated events. It is not necessary to invoke a scatter-free region near the sun in order to reproduce the fast rise times observed for directly-connected events. The new solution also yields a time-evolution for the vector anisotropy which agrees well with previously reported observations.

In addition, the new solution predicts that, during the decay phase, a typical convex spectral feature initially at energy T_o will move to lower energies at an exponential rate given by T_KINK = T_oexp(-t/Ƭ_KINK). Assuming adiabatic deceleration and a boundary at 2.7 AU, the solution yields Ƭ_KINK ≈ 100h, which is faster than the measured ~200h time constant and slower than the adiabatic rate of ~78h at 1 AU. Two possible explanations are that the boundary is at ~5 AU or that some other energy-change process is operative.