Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

2,4-dimethylene-1,3-cyclobutanediyl and 2,4-dimethylenebicyclobutane. Synthesis, spectroscopy, and reactivity of a triplet biradical and its covalent isomer

The preparation and direct observation of triplet 2,4-dimethylene-1,3- cyclobutanediyl (1), the non-Kekule isomer of benzene, is described. The biradical was generated by photolysis of 5,6-dimethylene-2,3- diazabicyclo[2.1.1]hex-2-ene (2) (which was synthesized in several steps from benzvalene) under cryogenic, matrix-isolation conditions. Biradical 1 was characterized by EPR spectroscopy (‌‌‌‌‌│D/hc│ =0.0204 cm^(-1), │E/hc│ =0.0028 cm^(-1)) and found to have a triplet ground state. The Δm_s= 2 transition displays hyperfine splitting attributed to a 7.3-G coupling to the ring methine and a 5.9-G coupling to the exocyclic methylene protons. Several experiments, including application of the magnetophotoselection (mps) technique in the generation of biradical 1, have allowed a determination of the zero-field triplet sublevels as x = -0.0040, y = +0.0136, and z = -0.0096 cm^(-1), where x and y are respectively the long and short in-plane axes and z the out-of-plane axis of 1.

Triplet 1 is yellow-orange and displays highly structured absorption (λ_(max)= 506 nm) and fluorescence (λ_(max) = 510 nm) spectra, with vibronic spacings of 1520 and 620 cm^(-1) for absorption and 1570 and 620 cm^(-1) for emission. The spectra were unequivocally assigned to triplet 1 by the use of a novel technique that takes advantage of the biradical's photolability. The absorption є = 7200 M^(-1) cm^(-1) and f = 0.022, establishing that the transition is spin-allowed. Further use of the mps technique has demonstrated that the transition is x-polarized, and the excited state 1s therefore of B_(1g) symmetry, in accord with theoretical predictions.

Thermolysis or direct photolysis of diazene 2 in fluid solution produces 2,4- dimethylenebicyclo[l.l.0]butane (3), whose ^(l)H NMR spectrum (-80°C, CD_(2)Cl_(2)) consists of singlets at δ 4.22 and 3.18 in a 2:1 ratio. Compound 3 is thermally unstable and dimerizes with second-order kinetics between -80 and -25°C (∆H^(‡) = 6.8 kcal mol^(-1), (∆s^(‡) = -28 eu) by a mechanism involving direct combination of two molecules of 3 in the rate-determining step. This singlet-manifold reaction ultimately produces a mixture of two dimers, 3,8,9- trimethylenetricyclo[5.1.1.0^(2,5)]non-4-ene (75) and trans-3,10-dimethylenetricyclo[6.2.0.0^(2,5)]deca-4,8-diene (76t), with the former predominating. In contrast, triplet-sensitized photolysis of 2, which leads to triplet 1, provides, in addition to 75 and 76t, a substantial amount of trans-5,10- dimethylenetricyclo[6.2.0.0^(3,6)]deca-3,8-diene (77t) and small amounts of two unidentified dimers.

In addition, triplet biradical 1 ring-closes to 3 in rigid media both thermally (77-140 K) and photochemically. In solution 3 forms triplet 1 upon energy transfer from sensitizers having relatively low triplet energies. The implications of the thermal chemistry for the energy surfaces of the system are discussed.

Nonlinear dispersive wave problems

The nonlinear partial differential equations for dispersive waves have special solutions representing uniform wavetrains. An expansion procedure is developed for slowly varying wavetrains, in which full nonlinearity is retained but in which the scale of the nonuniformity introduces a small parameter. The first order results agree with the results that Whitham obtained by averaging methods. The perturbation method provides a detailed description and deeper understanding, as well as a consistent development to higher approximations. This method for treating partial differential equations is analogous to the "multiple time scale" methods for ordinary differential equations in nonlinear vibration theory. It may also be regarded as a generalization of geometrical optics to nonlinear problems.

To apply the expansion method to the classical water wave problem, it is crucial to find an appropriate variational principle. It was found in the present investigation that a Lagrangian function equal to the pressure yields the full set of equations of motion for the problem. After this result is derived, the Lagrangian is compared with the more usual expression formed from kinetic minus potential energy. The water wave problem is then examined by means of the expansion procedure.

Nonlinear optics and wavelength translation via cavity-optomechanics

The field of cavity-optomechanics explores the interaction of light with sound in an ever increasing array of devices. This interaction allows the mechanical system to be both sensed and controlled by the optical system, opening up a wide variety of experiments including the cooling of the mechanical resonator to its quantum mechanical ground state and the squeezing of the optical field upon interaction with the mechanical resonator, to name two.

In this work we explore two very different systems with different types of optomechanical coupling. The first system consists of two microdisk optical resonators stacked on top of each other and separated by a very small slot. The interaction of the disks causes their optical resonance frequencies to be extremely sensitive to the gap between the disks. By careful control of the gap between the disks, the optomechanical coupling can be made to be quadratic to first order which is uncommon in optomechanical systems. With this quadratic coupling the light field is now sensitive to the energy of the mechanical resonator and can directly control the potential energy trapping the mechanical motion. This ability to directly control the spring constant without modifying the energy of the mechanical system, unlike in linear optomechanical coupling, is explored.

Next, the bulk of this thesis deals with a high mechanical frequency optomechanical crystal which is used to coherently convert photons between different frequencies. This is accomplished via the engineered linear optomechanical coupling in these devices. Both classical and quantum systems utilize the interaction of light and matter across a wide range of energies. These systems are often not naturally compatible with one another and require a means of converting photons of dissimilar wavelengths to combine and exploit their different strengths. Here we theoretically propose and experimentally demonstrate coherent wavelength conversion of optical photons using photon-phonon translation in a cavity-optomechanical system. For an engineered silicon optomechanical crystal nanocavity supporting a 4 GHz localized phonon mode, optical signals in a 1.5 MHz bandwidth are coherently converted over a 11.2 THz frequency span between one cavity mode at wavelength 1460 nm and a second cavity mode at 1545 nm with a 93% internal (2% external) peak efficiency. The thermal and quantum limiting noise involved in the conversion process is also analyzed and, in terms of an equivalent photon number signal level, are found to correspond to an internal noise level of only 6 and 4 times 10x^-3 quanta, respectively.

We begin by developing the requisite theoretical background to describe the system. A significant amount of time is then spent describing the fabrication of these silicon nanobeams, with an emphasis on understanding the specifics and motivation. The experimental demonstration of wavelength conversion is then described and analyzed. It is determined that the method of getting photons into the cavity and collected from the cavity is a fundamental limiting factor in the overall efficiency. Finally, a new coupling scheme is designed, fabricated, and tested that provides a means of coupling greater than 90% of photons into and out of the cavity, addressing one of the largest obstacles with the initial wavelength conversion experiment.

Spectroscopic studies of gas-phase molecular clusters

Spectroscopic investigations of hydrogen-bonding and van der Waals' interactions m molecular clusters were studied by the techniques of infrared predissociation and resonance-enhanced multiphoton ionization spectroscopies (REMPI). Ab initio calculations were applied in conjunction for data interpretation.

The infrared predissociation spectroscopy of CN^-•(H_2O)_n (n = 2 - 6) clusters was reported in the region of 2950 - 3850 cm^(-1). The hydrogen bondings for the C-site and N-site binding, and among the water molecules were identified for n = 2 to 4. A spectral transition was observed for n = 5 and 6, implying that the anion was surface-bound onto the water aggregates in larger clusters.

The infrared predissociation spectroscopy of Br^-•(NH_3) and I^-•(NH_3)_n (n =1-3) clusters was reported in the region of 3050-3450 cm^(-1). For the Br^-•(NH_3) complex, a dominating ionic NH stretch appeared at 3175 cm^(-1), and the weaker free NH stretch appeared at 3348 cm^(-1). The observed spectrum was consistent to the structure in which there was one nearly linear hydrogen bond between Br^- and the NH_3 moiety. For the I^- •(NH_3) complex, five distinct IR absorption bands were observed in the spectrum. The spectrum was not consistent with basic frequency patterns of three geometries considered in the ab initio calculations - complex with one, two and three hydrogen bondings between I^- and the NH_3 moiety. Substantial inhomogenous broadening were displayed in the spectra for I^-•(NH_3)_n (n =2-3), suggesting the presence of multiple isomers.

The REMPI spectroscopy of the bound 4p ^2П 1/2 and ^2П 3/2 states, and the dissociative 3d ^2Σ^+ 1/2 state in the Al•Ar complex was reported. The dissociative spectrum at Al^+ channel suggested the coupling of the 4p ^2П 1/2,3/2 states to the repulsive 3d ^2Σ^+ 1/2 state. The spin-electronic coupling was further manifested in the dissociative Al^+ spectrum of the 3d ^2Σ^+ 1/2 state. Using the potential energy curves obtained from ab initio calculations, a bound → continuum Franck-Condon-intensity simulation was performed and compared with the one-photon 3d ^2Σ^+ 1/2 profile. The agreement provided evidence for the petturbation above the Al(3d)Ar dissociation limit, and the repulsive character of the 3d ^2Σ^+ 1/2 state.

The Madden-Julian Oscillation : observation, modeling, and theory

The Madden-Julian Oscillation (MJO) is a pattern of intense rainfall and associated planetary-scale circulations in the tropical atmosphere, with a recurrence interval of 30-90 days. Although the MJO was first discovered 40 years ago, it is still a challenge to simulate the MJO in general circulation models (GCMs), and even with simple models it is difficult to agree on the basic mechanisms. This deficiency is mainly due to our poor understanding of moist convection—deep cumulus clouds and thunderstorms, which occur at scales that are smaller than the resolution elements of the GCMs. Moist convection is the most important mechanism for transporting energy from the ocean to the atmosphere. Success in simulating the MJO will improve our understanding of moist convection and thereby improve weather and climate forecasting.

We address this fundamental subject by analyzing observational datasets, constructing a hierarchy of numerical models, and developing theories. Parameters of the models are taken from observation, and the simulated MJO fits the data without further adjustments. The major findings include: 1) the MJO may be an ensemble of convection events linked together by small-scale high-frequency inertia-gravity waves; 2) the eastward propagation of the MJO is determined by the difference between the eastward and westward phase speeds of the waves; 3) the planetary scale of the MJO is the length over which temperature anomalies can be effectively smoothed by gravity waves; 4) the strength of the MJO increases with the typical strength of convection, which increases in a warming climate; 5) the horizontal scale of the MJO increases with the spatial frequency of convection; and 6) triggered convection, where potential energy accumulates until a threshold is reached, is important in simulating the MJO. Our findings challenge previous paradigms, which consider the MJO as a large-scale mode, and point to ways for improving the climate models.

Test particle motion in a Lorentz gas

This is a two-part thesis concerning the motion of a test particle in a bath. In part one we use an expansion of the operator PLe^it(1-P)LLP to shape the Zwanzig equation into a generalized Fokker-Planck equation which involves a diffusion tensor depending on the test particle's momentum and the time.

In part two the resultant equation is studied in some detail for the case of test particle motion in a weakly coupled Lorentz Gas. The diffusion tensor for this system is considered. Some of its properties are calculated; it is computed explicitly for the case of a Gaussian potential of interaction.

The equation for the test particle distribution function can be put into the form of an inhomogeneous Schroedinger equation. The term corresponding to the potential energy in the Schroedinger equation is considered. Its structure is studied, and some of its simplest features are used to find the Green's function in the limiting situations of low density and long time.

Storm track response to perturbations in climate

This thesis advances our understanding of midlatitude storm tracks and how they respond to perturbations in the climate system. The midlatitude storm tracks are regions of maximal turbulent kinetic energy in the atmosphere. Through them, the bulk of the atmospheric transport of energy, water vapor, and angular momentum occurs in midlatitudes. Therefore, they are important regulators of climate, controlling basic features such as the distribution of surface temperatures, precipitation, and winds in midlatitudes. Storm tracks are robustly projected to shift poleward in global-warming simulations with current climate models. Yet the reasons for this shift have remained unclear. Here we show that this shift occurs even in extremely idealized (but still three-dimensional) simulations of dry atmospheres. We use these simulations to develop an understanding of the processes responsible for the shift and develop a conceptual model that accounts for it.

We demonstrate that changes in the convective static stability in the deep tropics alone can drive remote shifts in the midlatitude storm tracks. Through simulations with a dry idealized general circulation model (GCM), midlatitude storm tracks are shown to be located where the mean available potential energy (MAPE, a measure of the potential energy available to be converted into kinetic energy) is maximal. As the climate varies, even if only driven by tropical static stability changes, the MAPE maximum shifts primarily because of shifts of the maximum of near-surface meridional temperature gradients. The temperature gradients shift in response to changes in the width of the tropical Hadley circulation, whose width is affected by the tropical static stability. Storm tracks generally shift in tandem with shifts of the subtropical terminus of the Hadley circulation.

We develop a one-dimensional diffusive energy-balance model that links changes in the Hadley circulation to midlatitude temperature gradients and so to the storm tracks. It is the first conceptual model to incorporate a dynamical coupling between the tropical Hadley circulation and midlatitude turbulent energy transport. Numerical and analytical solutions of the model elucidate the circumstances of when and how the storm tracks shift in tandem with the terminus of the Hadley circulation. They illustrate how an increase of only the convective static stability in the deep tropics can lead to an expansion of the Hadley circulation and a poleward shift of storm tracks.

The simulations with the idealized GCM and the conceptual energy-balance model demonstrate a clear link between Hadley circulation dynamics and midlatitude storm track position. With the help of the hierarchy of models presented in this thesis, we obtain a closed theory of storm track shifts in dry climates. The relevance of this theory for more realistic moist climates is discussed.

Vibrational pooling and constrained equilibration on surfaces

In this thesis, we provide a statistical theory for the vibrational pooling and fluorescence time dependence observed in infrared laser excitation of CO on an NaCl surface. The pooling is seen in experiment and in computer simulations. In the theory, we assume a rapid equilibration of the quanta in the substrate and minimize the free energy subject to the constraint at any time t of a fixed number of vibrational quanta N(t). At low incident intensity, the distribution is limited to one- quantum exchanges with the solid and so the Debye frequency of the solid plays a key role in limiting the range of this one-quantum domain. The resulting inverted vibrational equilibrium population depends only on fundamental parameters of the oscillator (ω_e and ω_eχ_e) and the surface (ω_D and T). Possible applications and relation to the Treanor gas phase treatment are discussed. Unlike the solid phase system, the gas phase system has no Debye-constraining maximum. We discuss the possible distributions for arbitrary N-conserving diatom-surface pairs, and include application to H:Si(111) as an example.

Computations are presented to describe and analyze the high levels of infrared laser-induced vibrational excitation of a monolayer of absorbed ¹³CO on a NaCl(100) surface. The calculations confirm that, for situations where the Debye frequency limited n domain restriction approximately holds, the vibrational state population deviates from a Boltzmann population linearly in n. Nonetheless, the full kinetic calculation is necessary to capture the result in detail.

We discuss the one-to-one relationship between N and γ and the examine the state space of the new distribution function for varied γ. We derive the Free Energy, F = NγkT − kTln(∑P_n), and effective chemical potential, μn ≈ γkT, for the vibrational pool. We also find the anti correlation of neighbor vibrations leads to an emergent correlation that appears to extend further than nearest neighbor.

Logarithmic Potential Theory on Riemann Surfaces

We develop a logarithmic potential theory on Riemann surfaces which generalizes logarithmic potential theory on the complex plane. We show the existence of an equilibrium measure and examine its structure. This leads to a formula for the structure of the equilibrium measure which is new even in the plane. We then use our results to study quadrature domains, Laplacian growth, and Coulomb gas ensembles on Riemann surfaces. We prove that the complement of the support of the equilibrium measure satisfies a quadrature identity. Furthermore, our setup allows us to naturally realize weak solutions of Laplacian growth (for a general time-dependent source) as an evolution of the support of equilibrium measures. When applied to the Riemann sphere this approach unifies the known methods for generating interior and exterior Laplacian growth. We later narrow our focus to a special class of quadrature domains which we call Algebraic Quadrature Domains. We show that many of the properties of quadrature domains generalize to this setting. In particular, the boundary of an Algebraic Quadrature Domain is the inverse image of a planar algebraic curve under a meromorphic function. This makes the study of the topology of Algebraic Quadrature Domains an interesting problem. We briefly investigate this problem and then narrow our focus to the study of the topology of classical quadrature domains. We extend the results of Lee and Makarov and prove (for n ≥ 3) c ≤ 5n-5, where c and n denote the connectivity and degree of a (classical) quadrature domain. At the same time we obtain a new upper bound on the number of isolated points of the algebraic curve corresponding to the boundary and thus a new upper bound on the number of special points. In the final chapter we study Coulomb gas ensembles on Riemann surfaces.