A note on the error analysis of classical Gram-Schmidt

An error analysis result is given for classical Gram--Schmidt factorization of a full rank matrix $A$ into $A=QR$ where $Q$ is left orthogonal (has orthonormal columns) and $R$ is upper triangular. The work presented here shows that the computed $R$ satisfies \normal{R}=\normal{A}+E where $E$ is an appropriately small backward error, but only if the diagonals of $R$ are computed in a manner similar to Cholesky factorization of the normal equations matrix. A similar result is stated in [Giraud at al, Numer. Math. 101(1):87--100,2005]. However, for that result to hold, the diagonals of $R$ must be computed in the manner recommended in this work.Comment: 12 pages This v2. v1 (from 2006) has not the biliographical reference set (at all). This is the only modification between v1 and v2. If you want to quote this paper, please quote the version published in Numerische Mathemati

Multiresolution analysis of active region magnetic structure and its correlation with the Mt. Wilson classification and flaring activity

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.Comment: Accepted for publication in Solar Physic

Multifragmentation of a very heavy nuclear system (I): Selection of single-source events

A sample of `single-source' events, compatible with the multifragmentation of very heavy fused systems, are isolated among well-measured 155Gd+natU 36AMeV reactions by examining the evolution of the kinematics of fragments with Z>=5 as a function of the dissipated energy and loss of memory of the entrance channel. Single-source events are found to be the result of very central collisions. Such central collisions may also lead to multiple fragment emission due to the decay of excited projectile- and target-like nuclei and so-called `neck' emission, and for this reason the isolation of single-source events is very difficult. Event-selection criteria based on centrality of collisions, or on the isotropy of the emitted fragments in each event, are found to be inefficient to separate the two mechanisms, unless they take into account the redistribution of fragments' kinetic energies into directions perpendicular to the beam axis. The selected events are good candidates to look for bulk effects in the multifragmentation process.Comment: 39 pages including 15 figures; submitted to Nucl. Phys.

Genetic, abiotic and social influences on sex differentiation in cichlid fishes and the evolution of sequential hermaphroditism

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/73799/1/j.1467-2979.2005.00184.x.pd

Effects of watershed land use on nitrogen concentrations and δ15 Nitrogen in groundwater

Author Posting. © The Authors, 2005. This is the author's version of the work. It is posted here by permission of Springer for personal use, not for redistribution. The definitive version was published in Biogeochemistry 77 (2006): 199-215, doi:10.1007/s10533-005-1036-2.Eutrophication is a major agent of change affecting freshwater, estuarine, and marine systems. It is largely driven by transportation of nitrogen from natural and anthropogenic sources. Research is needed to quantify this nitrogen delivery and to link the delivery to specific land-derived sources. In this study we measured nitrogen concentrations and δ15N values in seepage water entering three freshwater ponds and six estuaries on Cape Cod, Massachusetts and assessed how they varied with different types of land use. Nitrate concentrations and δ15N values in groundwater reflected land use in developed and pristine watersheds. In particular, watersheds with larger populations delivered larger nitrate loads with higher δ15N values to receiving waters. The enriched δ15N values confirmed nitrogen loading model results identifying wastewater contributions from septic tanks as the major N source. Furthermore, it was apparent that N coastal sources had a relatively larger impact on the N loads and isotopic signatures than did inland N sources further upstream in the watersheds. This finding suggests that management priorities could focus on coastal sources as a first course of action. This would require management constraints on a much smaller population.This work was supported by funds from the Woods Hole Oceanographic Institution Sea Grant Program, from the Cooperative Institute for Coastal and Estuarine Environmental Technology, from Massachusetts Department of Environmental Protection to Applied Science Associates, Narragansett, RI, as well as from Palmer/McLeod and NOAA National Estuarine Research Reserve Fellowships to Kevin Kroeger. This work is the result of research sponsored by NOAA National Sea Grant College Program Office, Department of Commerce, under Grant No. NA86RG0075, Woods Hole Oceanographic Institution Sea Grant Project No. R/M-40

Tomato: a crop species amenable to improvement by cellular and molecular methods

Tomato is a crop plant with a relatively small DNA content per haploid genome and a well developed genetics. Plant regeneration from explants and protoplasts is feasable which led to the development of efficient transformation procedures. In view of the current data, the isolation of useful mutants at the cellular level probably will be of limited value in the genetic improvement of tomato. Protoplast fusion may lead to novel combinations of organelle and nuclear DNA (cybrids), whereas this technique also provides a means of introducing genetic information from alien species into tomato. Important developments have come from molecular approaches. Following the construction of an RFLP map, these RFLP markers can be used in tomato to tag quantitative traits bred in from related species. Both RFLP's and transposons are in the process of being used to clone desired genes for which no gene products are known. Cloned genes can be introduced and potentially improve specific properties of tomato especially those controlled by single genes. Recent results suggest that, in principle, phenotypic mutants can be created for cloned and characterized genes and will prove their value in further improving the cultivated tomato.

Time-integrated luminosity recorded by the BABAR detector at the PEP-II e+e- collider

This article is the Preprint version of the final published artcile which can be accessed at the link below.We describe a measurement of the time-integrated luminosity of the data collected by the BABAR experiment at the PEP-II asymmetric-energy e+e- collider at the ϒ(4S), ϒ(3S), and ϒ(2S) resonances and in a continuum region below each resonance. We measure the time-integrated luminosity by counting e+e-→e+e- and (for the ϒ(4S) only) e+e-→μ+μ- candidate events, allowing additional photons in the final state. We use data-corrected simulation to determine the cross-sections and reconstruction efficiencies for these processes, as well as the major backgrounds. Due to the large cross-sections of e+e-→e+e- and e+e-→μ+μ-, the statistical uncertainties of the measurement are substantially smaller than the systematic uncertainties. The dominant systematic uncertainties are due to observed differences between data and simulation, as well as uncertainties on the cross-sections. For data collected on the ϒ(3S) and ϒ(2S) resonances, an additional uncertainty arises due to ϒ→e+e-X background. For data collected off the ϒ resonances, we estimate an additional uncertainty due to time dependent efficiency variations, which can affect the short off-resonance runs. The relative uncertainties on the luminosities of the on-resonance (off-resonance) samples are 0.43% (0.43%) for the ϒ(4S), 0.58% (0.72%) for the ϒ(3S), and 0.68% (0.88%) for the ϒ(2S).This work is supported by the US Department of Energy and National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), the Commissariat à l’Energie Atomique and Institut National de Physique Nucléaire et de Physiquedes Particules (France), the Bundesministerium für Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale di Fisica Nucleare (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), the Research Council of Norway, the Ministry of Education and Science of the Russian Federation, Ministerio de Ciencia e Innovación (Spain), and the Science and Technology Facilities Council (United Kingdom). Individuals have received support from the Marie-Curie IEF program (European Union) and the A.P. Sloan Foundation (USA)

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio

Observation of a new chi_b state in radiative transitions to Upsilon(1S) and Upsilon(2S) at ATLAS

The chi_b(nP) quarkonium states are produced in proton-proton collisions at the Large Hadron Collider (LHC) at sqrt(s) = 7 TeV and recorded by the ATLAS detector. Using a data sample corresponding to an integrated luminosity of 4.4 fb^-1, these states are reconstructed through their radiative decays to Upsilon(1S,2S) with Upsilon->mu+mu-. In addition to the mass peaks corresponding to the decay modes chi_b(1P,2P)->Upsilon(1S)gamma, a new structure centered at a mass of 10.530+/-0.005 (stat.)+/-0.009 (syst.) GeV is also observed, in both the Upsilon(1S)gamma and Upsilon(2S)gamma decay modes. This is interpreted as the chi_b(3P) system.Comment: 5 pages plus author list (18 pages total), 2 figures, 1 table, corrected author list, matches final version in Physical Review Letter