Lq estimates for rearrangements of Fourier series

This work is concerned with estimating the upper envelopes S* of the absolute values of the partial sums of rearranged trigonometric sums. A.M. Garsia [Annals of Math. 79 (1964), 634-9] gave an estimate for the L₂ norms of the S*, averaged over all rearrangements of the original (finite) sum. This estimate enabled him to prove that the Fourier series of any function in L₂ can be rearranged so that it converges a.e. The main result of this thesis is a similar estimate of the L_q norms of the S*, for all even integers q. This holds for finite linear combinations of functions which satisfy a condition which is a generalization of orthonormality in the L₂ case. This estimate for finite sums is extended to Fourier series of L_q functions; it is shown that there are functions to which the Men’shov-Paley Theorem does not apply, but whose Fourier series can nevertheless be rearranged so that the S* of the rearranged series is in L_q.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Oxidative chemisorption of SO_2 on γ - Al_2O_3 and deposition of H_2-permselective SiO_2 films

The interaction of SO_2 with γ - Al_2O_3 and the deposition of H_2 permselective SiO_2 films have been investigated. The adsorption and oxidative adsorption of SO_2 on γ - Al_2O_3 have been examined at temperatures 500-700°C by Fourier transform infrared spectroscopy (FTIR) and thermogravimetric analysis (TGA). At temperatures above 500°C most of SO_2 adsorbed on the strong sites on alumina. The adsorbed SO_2 species was characterized by an IR band at 1065 cm^(-1). The equilibrium coverage and initial rate of adsorption decreased with temperature suggesting a two-step adsorption. When γ - Al_2O_3 was contacted with a mixture of SO_2 and O_2, adsorption of SO_2 and oxidation of the adsorbed SO_2 to a surface sulfate characterized by broad IR bands at 1070 cm^(-1), 1390 cm^(-1) took place. The results of a series of TGA experiments under different atmospheres strongly suggest that surface SO_2 and surface sulfate involve the same active sites such that SO_2 adsorption is inhibited by already formed sulfate. The results also indicate a broad range of site strengths.

The desorption of adsorbed SO_2 and the reductive desorption of oxidatively adsorbed SO_2 have been investigated by microreactor experiments and thermogravimetric analysis (TGA). Temperature programmed reduction (TPR) of adsorbed SO_2 showed that SO_2 was desorbed without significant reaction with H_2 when H_2 concentration was low while considerable reaction occurred when 100% H_2 was used. SO_2 adsorbed on the strong sites on alumina was reduced to sulfur and H_2S. The isothermal reduction experiments of oxidatively adsorbed SO_2 reveal that the rate of reduction is very slow below 550°C even with 100% H_2. The reduction product is mainly composed of SO_2. TPR experiments of oxidatively adsorbed SO_2 showed that H_2S arose from a sulfate strongly chemisorbed on the surface.

Films of amorphous SiO_2 were deposited within the walls of porous Vycor tubes by SiH_4 oxidation in an opposing reactants geometry : SiH_4 was passed inside the tube while O_2 was passed outside the tube. The two reactants diffused opposite to each other and reacted within a narrow front inside the tube wall to form a thin SiO_2 film. Once the pores were plugged the reactants could not reach each other and the reaction stopped. At 450°C and 0.1 and 0.33 atm of SiH_4 and O_2, the reaction was complete within 15 minutes. The thickness of the SiO_2 film was estimated to be about 0.1 µm. Measurements of H_2 and N_2 permeation rates showed that the SiO_2 film was highly selective to H_2 permeation. The H_2:N_2 flux at 450°C varied between 2000-3000.

Thin SiO_2 films were heat treated in different gas mixtures to determine their stability in functioning as high-temperature hydrogen-permselective membranes. The films were heat-treated at 450-700°C in dry N_2, dry O_2, N_2-H_2O, and O_2-H_2O mixtures. The permeation rates of H_2 and N_2 changed depending on the original conditions of film formation as well as on the heat treatment. Heating in dry N_2 slowly reduced the permeation rates of both H_2 and N_2. Heating in a N_2-H_2O atmosphere led to a steeper decline of H_2 permeability. But the permeation rate of N_2 increased or decreased according to whether the film deposition had been carried out in the absence or presence of H_2O vapor, respectively. Thermal treatment in O_2 caused rapid decline of the permeation rates of H_2 and N_2 in films that were deposited under dry conditions. The decline was moderate in films deposited under wet conditions.

A new high-order Fourier continuation-based elasticity solver for complex three-dimensional geometries

This thesis presents a new approach for the numerical solution of three-dimensional problems in elastodynamics. The new methodology, which is based on a recently introduced Fourier continuation (FC) algorithm for the solution of Partial Differential Equations on the basis of accurate Fourier expansions of possibly non-periodic functions, enables fast, high-order solutions of the time-dependent elastic wave equation in a nearly dispersionless manner, and it requires use of CFL constraints that scale only linearly with spatial discretizations. A new FC operator is introduced to treat Neumann and traction boundary conditions, and a block-decomposed (sub-patch) overset strategy is presented for implementation of general, complex geometries in distributed-memory parallel computing environments. Our treatment of the elastic wave equation, which is formulated as a complex system of variable-coefficient PDEs that includes possibly heterogeneous and spatially varying material constants, represents the first fully-realized three-dimensional extension of FC-based solvers to date. Challenges for three-dimensional elastodynamics simulations such as treatment of corners and edges in three-dimensional geometries, the existence of variable coefficients arising from physical configurations and/or use of curvilinear coordinate systems and treatment of boundary conditions, are all addressed. The broad applicability of our new FC elasticity solver is demonstrated through application to realistic problems concerning seismic wave motion on three-dimensional topographies as well as applications to non-destructive evaluation where, for the first time, we present three-dimensional simulations for comparison to experimental studies of guided-wave scattering by through-thickness holes in thin plates.

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

The power of quantum Fourier sampling

How powerful are Quantum Computers? Despite the prevailing belief that Quantum Computers are more powerful than their classical counterparts, this remains a conjecture backed by little formal evidence. Shor's famous factoring algorithm [Shor97] gives an example of a problem that can be solved efficiently on a quantum computer with no known efficient classical algorithm. Factoring, however, is unlikely to be NP-Hard, meaning that few unexpected formal consequences would arise, should such a classical algorithm be discovered. Could it then be the case that any quantum algorithm can be simulated efficiently classically? Likewise, could it be the case that Quantum Computers can quickly solve problems much harder than factoring? If so, where does this power come from, and what classical computational resources do we need to solve the hardest problems for which there exist efficient quantum algorithms?

We make progress toward understanding these questions through studying the relationship between classical nondeterminism and quantum computing. In particular, is there a problem that can be solved efficiently on a Quantum Computer that cannot be efficiently solved using nondeterminism? In this thesis we address this problem from the perspective of sampling problems. Namely, we give evidence that approximately sampling the Quantum Fourier Transform of an efficiently computable function, while easy quantumly, is hard for any classical machine in the Polynomial Time Hierarchy. In particular, we prove the existence of a class of distributions that can be sampled efficiently by a Quantum Computer, that likely cannot be approximately sampled in randomized polynomial time with an oracle for the Polynomial Time Hierarchy.

Our work complements and generalizes the evidence given in Aaronson and Arkhipov's work [AA2013] where a different distribution with the same computational properties was given. Our result is more general than theirs, but requires a more powerful quantum sampler.