I. Singular perturbation problems involving singular points and turning points. II. On the averaged Lagrangian technique for nonlinear dispersive waves

In Part I a class of linear boundary value problems is considered which is a simple model of boundary layer theory. The effect of zeros and singularities of the coefficients of the equations at the point where the boundary layer occurs is considered. The usual boundary layer techniques are still applicable in some cases and are used to derive uniform asymptotic expansions. In other cases it is shown that the inner and outer expansions do not overlap due to the presence of a turning point outside the boundary layer. The region near the turning point is described by a two-variable expansion. In these cases a related initial value problem is solved and then used to show formally that for the boundary value problem either a solution exists, except for a discrete set of eigenvalues, whose asymptotic behaviour is found, or the solution is non-unique. A proof is given of the validity of the two-variable expansion; in a special case this proof also demonstrates the validity of the inner and outer expansions.

Nonlinear dispersive wave equations which are governed by variational principles are considered in Part II. It is shown that the averaged Lagrangian variational principle is in fact exact. This result is used to construct perturbation schemes to enable higher order terms in the equations for the slowly varying quantities to be calculated. A simple scheme applicable to linear or near-linear equations is first derived. The specific form of the first order correction terms is derived for several examples. The stability of constant solutions to these equations is considered and it is shown that the correction terms lead to the instability cut-off found by Benjamin. A general stability criterion is given which explicitly demonstrates the conditions under which this cut-off occurs. The corrected set of equations are nonlinear dispersive equations and their stationary solutions are investigated. A more sophisticated scheme is developed for fully nonlinear equations by using an extension of the Hamiltonian formalism recently introduced by Whitham. Finally the averaged Lagrangian technique is extended to treat slowly varying multiply-periodic solutions. The adiabatic invariants for a separable mechanical system are derived by this method.

A singularly perturbed linear two-point boundary-value problem

We consider the following singularly perturbed linear two-point boundary-value problem:

Ly(x) ≡ Ω(ε)D_xy(x) - A(x,ε)y(x) = f(x,ε) 0≤x≤1 (1a)

By ≡ L(ε)y(0) + R(ε)y(1) = g(ε) ε → 0^+ (1b)

Here Ω(ε) is a diagonal matrix whose first m diagonal elements are 1 and last m elements are ε. Aside from reasonable continuity conditions placed on A, L, R, f, g, we assume the lower right mxm principle submatrix of A has no eigenvalues whose real part is zero. Under these assumptions a constructive technique is used to derive sufficient conditions for the existence of a unique solution of (1). These sufficient conditions are used to define when (1) is a regular problem. It is then shown that as ε → 0^+ the solution of a regular problem exists and converges on every closed subinterval of (0,1) to a solution of the reduced problem. The reduced problem consists of the differential equation obtained by formally setting ε equal to zero in (1a) and initial conditions obtained from the boundary conditions (1b). Several examples of regular problems are also considered.

A similar technique is used to derive the properties of the solution of a particular difference scheme used to approximate (1). Under restrictions on the boundary conditions (1b) it is shown that for the stepsize much larger than ε the solution of the difference scheme, when applied to a regular problem, accurately represents the solution of the reduced problem.

Furthermore, the existence of a similarity transformation which block diagonalizes a matrix is presented as well as exponential bounds on certain fundamental solution matrices associated with the problem (1).

Topics in randomized numerical linear algebra

This thesis studies three classes of randomized numerical linear algebra algorithms, namely: (i) randomized matrix sparsification algorithms, (ii) low-rank approximation algorithms that use randomized unitary transformations, and (iii) low-rank approximation algorithms for positive-semidefinite (PSD) matrices.

Randomized matrix sparsification algorithms set randomly chosen entries of the input matrix to zero. When the approximant is substituted for the original matrix in computations, its sparsity allows one to employ faster sparsity-exploiting algorithms. This thesis contributes bounds on the approximation error of nonuniform randomized sparsification schemes, measured in the spectral norm and two NP-hard norms that are of interest in computational graph theory and subset selection applications.

Low-rank approximations based on randomized unitary transformations have several desirable properties: they have low communication costs, are amenable to parallel implementation, and exploit the existence of fast transform algorithms. This thesis investigates the tradeoff between the accuracy and cost of generating such approximations. State-of-the-art spectral and Frobenius-norm error bounds are provided.

The last class of algorithms considered are SPSD "sketching" algorithms. Such sketches can be computed faster than approximations based on projecting onto mixtures of the columns of the matrix. The performance of several such sketching schemes is empirically evaluated using a suite of canonical matrices drawn from machine learning and data analysis applications, and a framework is developed for establishing theoretical error bounds.

In addition to studying these algorithms, this thesis extends the Matrix Laplace Transform framework to derive Chernoff and Bernstein inequalities that apply to all the eigenvalues of certain classes of random matrices. These inequalities are used to investigate the behavior of the singular values of a matrix under random sampling, and to derive convergence rates for each individual eigenvalue of a sample covariance matrix.

^(18)O/^(16)O and D/H studies on the interactions between heated meteoric ground waters and igneous intrusions: Western San Juan Mountains, Colorado and the Isle of Skye, Scotland

Detailed oxygen, hydrogen and carbon isotope studies have been carried out on igneous and metamorphic rocks of the Stony Mountain complex, Colorado, and the Isle of Skye, Scotland, in order to better understand the problems of hydrothermal meteoric water-rock interaction.

The Tertiary Stony Mountain stock (~1.3 km in diameter), is composed of an outer diorite, a main mass of biotite gabbro, and an inner diorite. The entire complex and most of the surrounding country rocks have experienced various degrees of ¹⁸O depletion (up to 10 per mil) due to interaction with heated meteoric waters. The inner diorite apparently formed from a low-¹⁸O magma with δ¹⁸O ≃ +2.5, but most of the isotopic effects are a result of exchange between H₂O and solidified igneous rocks. The low-¹⁸O inner diorite magma was probably produced by massive assimilation and/or melting of hydrothermally altered country rocks. The δ¹⁸O values of the rocks generally increase with increasing grain size, except that quartz typically has δ¹⁸O = +6 to +8, and is more resistant to hydrothermal exchange than any other mineral studied. Based on atom % oxygen, the outer diorites, gabbros, and volcanic rocks exhibit integrated water/rock ratios of 0.3 ± 0.2, 0.15 ± 0.1, and 0.2 ± 0.1, respectively. Locally, water/rock ratios attain values greater than 1.0. Hydrogen isotopic analyses of sericites, chlorites, biotites, and amphiboles range from -117 to -150. δD in biotites varies inversely with Fe/Fe+Mg, as predicted by Suzuoki and Epstein (1974), and positively with elevation, over a range of 600 m. The calculated δD of the mid-to-late-Tertiary meteoric waters is about -100. Carbonate δ¹³C values average -5.5 (PDB), within the generally accepted range for deep-seated carbon.

Almost all the rocks within 4 km of the central Tertiary intrusive complex of Skye are depleted in ¹⁸O. Whole-rock δ¹⁸O values of basalts (-7. 1 to +8.4), Mesozoic shales (-0.6 to + 12.4), and Precambrian sandstones (-6.2 to + 10.8) systematically decrease inward towards the center of the complex. The Cuillin gabbro may have formed from a ¹⁸O-depleted magma (depleted by about 2 per mil); δ¹⁸O of plagioclase (-7.1 to + 2.5) and pyroxene (-0.5 to + 3.2) decrease outward toward the margins of the pluton. The Red Hills epigranite plutons have δ¹⁸O quartz (-2.7 to + 7.6) and feldspar (-6.7 to + 6.0) that suggest about 3/4 of the exchange took place at subsolidus temperatures; profound disequilibrium quartz-feldspar fractionations (up to 12) are characteristic. The early epigranites were intruded as low-¹⁸O melts (depletions of up to 3 per mil) with δ¹⁸O of the primary, igneous quartz decreasing progressively with time. The Southern Porphyritic Epigranite was apparently intruded as a low-¹⁸O magma with δ¹⁸O ≃ -2.6. A good correlation exists between grain size and δ¹⁸O for the unique, high-¹⁸O Beinn an Dubhaich granite which intrudes limestone having a δ¹⁸O range of +0.5 to +20.8, and δ¹³C of -4.9 to -1.0. The δD values of sericites (-104 to -107), and amphiboles, chlorites, and biotites (-105 to -128) from the igneous rocks , indicate that Eocene surface waters at Skye had δD ≃ -90. The average water/rock ratio for the Skye hydrothermal system is approximately one; at least 2000 km³ of heated meteoric waters were cycled through these rocks.

Thus these detailed isotopic studies of two widely separated areas indicate that (1) ¹⁸O-depleted magmas are commonly produced in volcanic terranes invaded by epizonal intrusions; (2) most of the ¹⁸O-depletion in such areas are a result of subsolidus exchange (particularly of feldspars); however correlation of δ¹⁸O with grain size is generally preserved only for systems that have undergone relatively minor meteoric hydrothermal exchange; (3) feldspar and calcite are the minerals mos t susceptible to oxygen isotopic exchange, whereas quartz is very resistant to oxygen isotope exchange; biotite, magnetite, and pyroxene have intermediate susceptibilities; and (4) basaltic country rocks are much more permeable to the hydrothermal convective system than shale, sandstone, or the crystalline basement complex.

Studies directed toward the synthesis of palau'amine and axinellamines A-D

The use of spiro [2.4]hepta-4,6-diene-1-methanol 7 as a general precursor for the synthesis of highly functionalized cyclopentyl rings is described. Diene 7 was converted to its silyl protected 4-nitrile derivative 24 in 46% overall yield. The cyclopropyl ring of 24 reacted with soft carbanionic nucleophiles to give ring opened homo-conjugate addition products 25a-h in 76-97% yield without loss of optical purity. The addition products could be further manipulated by selective mono-hydrogenation to give 1,2 substituted cyclopentenes 26a-e in 85-96% yield.

Diene 7 was used as a starting material for studies directed toward the synthesis of the stereochemically dense chloro-cyclopentyl core of palau'amine 1. Two advanced intermediates 50 and 72 were synthesized. Attempts to effect intramolecular chlorine transfer with 50 were unsuccessful. Attempted intramolecular chlorine transfer with 72 led, instead, to an oxygenated species resulting from oxygen radical trapping.

The enantioselective synthesis of the stereochemically dense chloro-cyclopenty l core of axinellamines A-D 2-5 starting from 7 is also described. The core is synthesized in 4.6% yield over 24 steps. Nakamura's radical dehalogenative hydroxylation is applied for the first time to a cyclopropyl carbonyl iodide to give the ring-opened product in 86% yield. Bolm's meso-anhydride desymmetrization is used to introduce asymmetry in a norbornene intermediate. The final step is a diastereoselective intermolecular chlorination using Barton's methodology to achieve chlorine transfer in 76% yield.

Biophysics of V(D)J recombination and Genome Packaging: In singulo studies on RAG, HMGB1, and TFAM

The recombination-activating gene products, RAG1 and RAG2, initiate V(D)J recombination during lymphocyte development by cleaving DNA adjacent to conserved recombination signal sequences (RSSs). The reaction involves DNA binding, synapsis, and cleavage at two RSSs located on the same DNA molecule and results in the assembly of antigen receptor genes. Since their discovery full-length, RAG1 and RAG2 have been difficult to purify, and core derivatives are shown to be most active when purified from adherent 293-T cells. However, the protein yield from adherent 293-T cells is limited. Here we develop a human suspension cell purification and change the expression vector to boost RAG production 6-fold. We use these purified RAG proteins to investigate V(D)J recombination on a mechanistic single molecule level. As a result, we are able to measure the binding statistics (dwell times and binding energies) of the initial RAG binding events with or without its co-factor high mobility group box protein 1 (HMGB1), and to characterize synapse formation at the single-molecule level yielding insights into the distribution of dwell times in the paired complex and the propensity for cleavage upon forming the synapse. We then go on to investigate HMGB1 further by measuring it compact single DNA molecules. We observed concentration dependent DNA compaction, differential DNA compaction depending on the divalent cation type, and found that at a particular HMGB1 concentration the percentage of DNA compacted is conserved across DNA lengths. Lastly, we investigate another HMGB protein called TFAM, which is essential for packaging the mitochondrial genome. We present crystal structures of TFAM bound to the heavy strand promoter 1 (HSP1) and to nonspecific DNA. We show TFAM dimerization is dispensable for DNA bending and transcriptional activation, but is required for mtDNA compaction. We propose that TFAM dimerization enhances mtDNA compaction by promoting looping of mtDNA.

Studies in nuclear reactor dynamics : I. The accuracy of point kinetics. II. The effect of delayed neutrons on the spectrum of the group-diffusion operator

This thesis is a theoretical work on the space-time dynamic behavior of a nuclear reactor without feedback. Diffusion theory with G-energy groups is used.

In the first part the accuracy of the point kinetics (lumped-parameter description) model is examined. The fundamental approximation of this model is the splitting of the neutron density into a product of a known function of space and an unknown function of time; then the properties of the system can be averaged in space through the use of appropriate weighting functions; as a result a set of ordinary differential equations is obtained for the description of time behavior. It is clear that changes of the shape of the neutron-density distribution due to space-dependent perturbations are neglected. This results to an error in the eigenvalues and it is to this error that bounds are derived. This is done by using the method of weighted residuals to reduce the original eigenvalue problem to that of a real asymmetric matrix. Then Gershgorin-type theorems .are used to find discs in the complex plane in which the eigenvalues are contained. The radii of the discs depend on the perturbation in a simple manner.

In the second part the effect of delayed neutrons on the eigenvalues of the group-diffusion operator is examined. The delayed neutrons cause a shifting of the prompt-neutron eigenvalue s and the appearance of the delayed eigenvalues. Using a simple perturbation method this shifting is calculated and the delayed eigenvalues are predicted with good accuracy.

La mujer investigadora en el País Vasco: análisis de la presencia de la mujer en actividades de I+D y el gasto en I+D

[ES]La apuesta por la igualdad en cualquier ámbito de nuestro entorno es un aspecto clave para la mejora de él. La igualdad en la investigación científica y desarrollo tecnológico es todavía un reto entre las instituciones que componen el sistema. En el País Vasco, existe un potente sistema y estructura de apoyo institucional a la ciencia, la tecnología y la innovación, y es imprescindible su correcto funcionamiento para un adecuado desarrollo de la I+D en todas sus vertientes. Para conocer la situación actual del País Vasco, este trabajo de investigación tiene como objetivo estudiar la relación existente entre la evolución del gasto en I+D y la evolución de la presencia de la mujer investigadora en actividades de I+D en el País Vasco. En general, el estudio nos permite concluir que la presencia de la mujer investigadora y el gasto en I+D han evolución de manera similar en los últimos veinte años.