Part 1: Longitudinal static and dynamic effects in semiconductor lasers. Part II: Spectral characteristics of passively mode-locked quantum well lasers

In the first part of this thesis a study of the effect of the longitudinal distribution of optical intensity and electron density on the static and dynamic behavior of semiconductor lasers is performed. A static model for above threshold operation of a single mode laser, consisting of multiple active and passive sections, is developed by calculating the longitudinal optical intensity distribution and electron density distribution in a self-consistent manner. Feedback from an index and gain Bragg grating is included, as well as feedback from discrete reflections at interfaces and facets. Longitudinal spatial holeburning is analyzed by including the dependence of the gain and the refractive index on the electron density. The mechanisms of spatial holeburning in quarter wave shifted DFB lasers are analyzed. A new laser structure with a uniform optical intensity distribution is introduced and an implementation is simulated, resulting in a large reduction of the longitudinal spatial holeburning effect.

A dynamic small-signal model is then developed by including the optical intensity and electron density distribution, as well as the dependence of the grating coupling coefficients on the electron density. Expressions are derived for the intensity and frequency noise spectrum, the spontaneous emission rate into the lasing mode, the linewidth enhancement factor, and the AM and FM modulation response. Different chirp components are identified in the FM response, and a new adiabatic chirp component is discovered. This new adiabatic chirp component is caused by the nonuniform longitudinal distributions, and is found to dominate at low frequencies. Distributed feedback lasers with partial gain coupling are analyzed, and it is shown how the dependence of the grating coupling coefficients on the electron density can result in an enhancement of the differential gain with an associated enhancement in modulation bandwidth and a reduction in chirp.

In the second part, spectral characteristics of passively mode-locked two-section multiple quantum well laser coupled to an external cavity are studied. Broad-band wavelength tuning using an external grating is demonstrated for the first time in passively mode-locked semiconductor lasers. A record tuning range of 26 nm is measured, with pulse widths of typically a few picosecond and time-bandwidth products of more than 10 times the transform limit. It is then demonstrated that these large time-bandwidth products are due to a strong linear upchirp, by performing pulse compression by a factor of 15 to a record pulse widths as low 320 fs.

A model for pulse propagation through a saturable medium with self-phase-modulation, due to the a-parameter, is developed for quantum well material, including the frequency dependence of the gain medium. This model is used to simulate two-section devices coupled to an external cavity. When no self-phase-modulation is present, it is found that the pulses are asymmetric with a sharper rising edge, that the pulse tails have an exponential behavior, and that the transform limit is 0.3. Inclusion of self-phase-modulation results in a linear upchirp imprinted on the pulse after each round-trip. This linear upchirp is due to a combination of self-phase-modulation in a gain section and absorption of the leading edge of the pulse in the saturable absorber.

Dirac spectra, summation formulae, and the spectral action

Noncommutative geometry is a source of particle physics models with matter Lagrangians coupled to gravity. One may associate to any noncommutative space (A, H, D) its spectral action, which is defined in terms of the Dirac spectrum of its Dirac operator D. When viewing a spin manifold as a noncommutative space, D is the usual Dirac operator. In this paper, we give nonperturbative computations of the spectral action for quotients of SU(2), Bieberbach manifolds, and SU(3) equipped with a variety of geometries. Along the way we will compute several Dirac spectra and refer to applications of this computation.

Characterization of exoplanet atmospheres : spectral retrieval and chemistry

The study of exoplanets is rapidly evolving into an important and exciting field of its own. My investigations over the past half-decade have focused on understanding just a small sliver of what they are trying to tell us. That small sliver is their atmospheres. Atmospheres are the buffer between the bulk planet and the vacuum of space. The atmosphere is an important component of a planet as it is the most readily observable and contains the most information about the physical processes that can occur in a planet. I have focused on two aspects of exoplanetary atmospheres. First, I aimed to understand the chemical mechanisms that control the atmospheric abundances. Second, I focused on interpreting exoplanet atmospheric spectra and what they tell us about the temperatures and compositions through inverse modeling. Finally, I interpreted the retrieved temperature and abundances from inverse modeling in the context of chemical disequilibrium in the planetary atmospheres.

A variational framework for spectral discretization of the density matrix in Kohn Sham density functional theory

Kohn-Sham density functional theory (KSDFT) is currently the main work-horse of quantum mechanical calculations in physics, chemistry, and materials science. From a mechanical engineering perspective, we are interested in studying the role of defects in the mechanical properties in materials. In real materials, defects are typically found at very small concentrations e.g., vacancies occur at parts per million, dislocation density in metals ranges from $10^{10} m^{-2}$ to $10^{15} m^{-2}$, and grain sizes vary from nanometers to micrometers in polycrystalline materials, etc. In order to model materials at realistic defect concentrations using DFT, we would need to work with system sizes beyond millions of atoms. Due to the cubic-scaling computational cost with respect to the number of atoms in conventional DFT implementations, such system sizes are unreachable. Since the early 1990s, there has been a huge interest in developing DFT implementations that have linear-scaling computational cost. A promising approach to achieving linear-scaling cost is to approximate the density matrix in KSDFT. The focus of this thesis is to provide a firm mathematical framework to study the convergence of these approximations. We reformulate the Kohn-Sham density functional theory as a nested variational problem in the density matrix, the electrostatic potential, and a field dual to the electron density. The corresponding functional is linear in the density matrix and thus amenable to spectral representation. Based on this reformulation, we introduce a new approximation scheme, called spectral binning, which does not require smoothing of the occupancy function and thus applies at arbitrarily low temperatures. We proof convergence of the approximate solutions with respect to spectral binning and with respect to an additional spatial discretization of the domain. For a standard one-dimensional benchmark problem, we present numerical experiments for which spectral binning exhibits excellent convergence characteristics and outperforms other linear-scaling methods.

Gauge Invariant Spectral Cauchy Characteristic Extraction of Gravitational Waves in Computational General Relativity

We present a complete system for Spectral Cauchy characteristic extraction (Spectral CCE). Implemented in C++ within the Spectral Einstein Code (SpEC), the method employs numerous innovative algorithms to efficiently calculate the Bondi strain, news, and flux.

Spectral CCE was envisioned to ensure physically accurate gravitational wave-forms computed for the Laser Interferometer Gravitational wave Observatory (LIGO) and similar experiments, while working toward a template bank with more than a thousand waveforms to span the binary black hole (BBH) problem’s seven-dimensional parameter space.

The Bondi strain, news, and flux are physical quantities central to efforts to understand and detect astrophysical gravitational wave sources within the Simulations of eXtreme Spacetime (SXS) collaboration, with the ultimate aim of providing the first strong field probe of the Einstein field equation.

In a series of included papers, we demonstrate stability, convergence, and gauge invariance. We also demonstrate agreement between Spectral CCE and the legacy Pitt null code, while achieving a factor of 200 improvement in computational efficiency.

Spectral CCE represents a significant computational advance. It is the foundation upon which further capability will be built, specifically enabling the complete calculation of junk-free, gauge-free, and physically valid waveform data on the fly within SpEC.

Fast numerical methods for mixed, singular Helmholtz boundary value problems and Laplace eigenvalue problems - with applications to antenna design, sloshing, electromagnetic scattering and spectral geometry

This thesis presents a novel class of algorithms for the solution of scattering and eigenvalue problems on general two-dimensional domains under a variety of boundary conditions, including non-smooth domains and certain "Zaremba" boundary conditions - for which Dirichlet and Neumann conditions are specified on various portions of the domain boundary. The theoretical basis of the methods for the Zaremba problems on smooth domains concern detailed information, which is put forth for the first time in this thesis, about the singularity structure of solutions of the Laplace operator under boundary conditions of Zaremba type. The new methods, which are based on use of Green functions and integral equations, incorporate a number of algorithmic innovations, including a fast and robust eigenvalue-search algorithm, use of the Fourier Continuation method for regularization of all smooth-domain Zaremba singularities, and newly derived quadrature rules which give rise to high-order convergence even around singular points for the Zaremba problem. The resulting algorithms enjoy high-order convergence, and they can tackle a variety of elliptic problems under general boundary conditions, including, for example, eigenvalue problems, scattering problems, and, in particular, eigenfunction expansion for time-domain problems in non-separable physical domains with mixed boundary conditions.

Polaron hopping in olivine phosphates studied by nuclear resonant scattering

Valence fluctuations of Fe²⁺ and Fe³⁺ were studied in a solid solution of Li_xFePO₄ by nuclear resonant forward scattering of synchrotron x rays while the sample was heated in a diamond-anvil pressure cell. The spectra acquired at different temperatures and pressures were analyzed for the frequencies of valence changes using the Blume-Tjon model of a system with a fluctuating Hamil- tonian. These frequencies were analyzed to obtain activation energies and an activation volume for polaron hopping. There was a large suppression of hopping frequency with pressure, giving an anomalously large activation volume. This large, positive value is typical of ion diffusion, which indicates correlated motions of polarons, and Li⁺ ions that alter the dynamics of both.

In a parallel study of Na_xFePO₄, the interplay between sodium ordering and electron mobility was investigated using a combination of synchrotron x-ray diffraction and nuclear resonant scattering. Conventional Mossbauer spectra were collected while the sample was heated in a resistive furnace. An analysis of the temperature evolution of the spectral shapes was used to identify the onset of fast electron hopping and determine the polaron hopping rate. Synchrotron x-ray diffraction measurements were carried out in the same temperature range. Reitveld analysis of the diffraction patterns was used to determine the temperature of sodium redistribution on the lattice. The diffraction analysis also provides new information about the phase stability of the system. The temperature evolution of the iron site occupancies from the Mossbauer measurements, combined with the synchrotron diffraction results give strong evidence for a relationship between the onset of fast electron dynamics and the redistribution of sodium in the lattice.

Measurements of activation barriers for polaron hopping gave fundamental insights about the correlation between electronic carriers and mobile ions. This work established that polaron-ion interactions can alter the local dynamics of electron and ion transport. These types of coupled processes may be common in many materials used for battery electrodes, and new details concerning the influence of polaron-ion interactions on the charge dynamics are relevant to optimizing their electrochemical performance.

Spectral density of first order Piecewise linear systems excited by white noise

The Fokker-Planck (FP) equation is used to develop a general method for finding the spectral density for a class of randomly excited first order systems. This class consists of systems satisfying stochastic differential equations of form ẋ + f(x) = m/Ʃ/j = 1 h_j(x)n_j(t) where f and the h_j are piecewise linear functions (not necessarily continuous), and the n_j are stationary Gaussian white noise. For such systems, it is shown how the Laplace-transformed FP equation can be solved for the transformed transition probability density. By manipulation of the FP equation and its adjoint, a formula is derived for the transformed autocorrelation function in terms of the transformed transition density. From this, the spectral density is readily obtained. The method generalizes that of Caughey and Dienes, J. Appl. Phys., 32.11.

This method is applied to 4 subclasses: (1) m = 1, h₁ = const. (forcing function excitation); (2) m = 1, h₁ = f (parametric excitation); (3) m = 2, h₁ = const., h₂ = f, n₁ and n₂ correlated; (4) the same, uncorrelated. Many special cases, especially in subclass (1), are worked through to obtain explicit formulas for the spectral density, most of which have not been obtained before. Some results are graphed.

Dealing with parametrically excited first order systems leads to two complications. There is some controversy concerning the form of the FP equation involved (see Gray and Caughey, J. Math. Phys., 44.3); and the conditions which apply at irregular points, where the second order coefficient of the FP equation vanishes, are not obvious but require use of the mathematical theory of diffusion processes developed by Feller and others. These points are discussed in the first chapter, relevant results from various sources being summarized and applied. Also discussed is the steady-state density (the limit of the transition density as t → ∞).

I. Configurational stability and redistribution equilibria in organomagnesium compounds. II. Redistribution equilibria in organocadmium compounds

I. CONFIGURATIONAL STABILITY AND REDISTRIBUTION EQUILIBRIA IN ORGANOMAGNESIUM COMPOUNDS

The dependence of the rate of inversion of a dialkylmagnesium compound on the solvent has been studied.

Examination of the temperature dependence of the nuclear magnetic resonance spectrum of 1-phenyl-2-propylmagnesium bromide in diethyl ether solution indicates that inversion of configuration at the methylene group of this Grignard reagent occurs with an approximate rate of 2 sec^-1 at room temperature. This is the first example of a rapid inversion rate in a secondary Grignard reagent.

The rates of exchange of alkyl groups between dineopentylmagnesium and di-s-butylmagnesium, bis-(2-methylbutyl)-magnesium and bis-(4, 4-dimethyl-2-pentyl)-magnesium respectively in diethyl ether solution were found to be fast on the nmr time scale. However, the alkyl group exchange rate was found to be slow in a diethyl ether solution of dineopentylmagnesium and bis-(2-methylbutyl)-magnesium containing N, N, N', N'-tetramethylethylenediamine. The unsymmetrical species neopentyl-2-methylbutyl-magnesium was observed at room temperature in the nmr spectrum of the solution containing the diamine.

II. REDISTRIBUTION EQUILIBRIA IN ORGANOCADMIUM COMPOUNDS

The exchange of methyl groups in dimethylcadmium has been studied by nuclear magnetic resonance spectroscopy. Activation parameters for the methyl group exchange have been measured for a neat sample and for a solution in tetrahydrofuran. The exchange is faster in the basic solvent tetrahydrofuran relative to the neat sample and in tetrahydrofuran solution is retarded by the solvating agent N, N, N’, N’-tetramethylethylenediamine and greatly increased by cadmium bromide. The addition of methanol to a solution of dimethylcadmium in tetrahydrofuran appears to have very little effect on the rate of exchange. The exchange was found to proceed with retention of configuration. The rate-limiting step for the exchange of methyl groups in a basic solvent appears to be the dissociation of coordinating solvent from dimethylcadmium.

The equilibrium between methylcadmium bromide, dimethylcadmium and cadmium bromide in tetrahydrofuran solution has also been studied. At room temperature the interconversion of the species is very fast on the nmr time scale but at -100° distinct absorptions for methylcadmium bromide and imethylcadmium are observed.

The species ethylmethylcadmium has been observed in the nmr spectrum.

The rate of exchange of vinyl groups in a solution of divinylcadmium in tetrahydrofuran has been found to be fast on the nmr time scale.

A mathematical and experimental investigation of microcycle spectral estimates of semiconductor flicker noise

The experimental portion of this thesis tries to estimate the density of the power spectrum of very low frequency semiconductor noise, from 10^-6.3 cps to 1. cps with a greater accuracy than that achieved in previous similar attempts: it is concluded that the spectrum is 1/f^α with α approximately 1.3 over most of the frequency range, but appearing to have a value of about 1 in the lowest decade. The noise sources are, among others, the first stage circuits of a grounded input silicon epitaxial operational amplifier. This thesis also investigates a peculiar form of stationarity which seems to distinguish flicker noise from other semiconductor noise.

In order to decrease by an order of magnitude the pernicious effects of temperature drifts, semiconductor "aging", and possible mechanical failures associated with prolonged periods of data taking, 10 independent noise sources were time-multiplexed and their spectral estimates were subsequently averaged. If the sources have similar spectra, it is demonstrated that this reduces the necessary data-taking time by a factor of 10 for a given accuracy.

In view of the measured high temperature sensitivity of the noise sources, it was necessary to combine the passive attenuation of a special-material container with active control. The noise sources were placed in a copper-epoxy container of high heat capacity and medium heat conductivity, and that container was immersed in a temperature controlled circulating ethylene-glycol bath.

Other spectra of interest, estimated from data taken concurrently with the semiconductor noise data were the spectra of the bath's controlled temperature, the semiconductor surface temperature, and the power supply voltage amplitude fluctuations. A brief description of the equipment constructed to obtain the aforementioned data is included.

The analytical portion of this work is concerned with the following questions: what is the best final spectral density estimate given 10 statistically independent ones of varying quality and magnitude? How can the Blackman and Tukey algorithm which is used for spectral estimation in this work be improved upon? How can non-equidistant sampling reduce data processing cost? Should one try to remove common trands shared by supposedly statistically independent noise sources and, if so, what are the mathematical difficulties involved? What is a physically plausible mathematical model that can account for flicker noise and what are the mathematical implications on its statistical properties? Finally, the variance of the spectral estimate obtained through the Blackman/Tukey algorithm is analyzed in greater detail; the variance is shown to diverge for α ≥ 1 in an assumed power spectrum of k/|f|^α, unless the assumed spectrum is "truncated".