Accelerating BST Methods for Model Reduction with Graphics Processors

Model order reduction of dynamical linear time-invariant system appears in many scientific and engineering applications. Numerically reliable SVD-based methods for this task require O(n3) floating-point arithmetic operations, with n being in the range 103 − 105 for many practical applications. In this paper we investigate the use of graphics processors (GPUs) to accelerate model reduction of large-scale linear systems via Balanced Stochastic Truncation, by off-loading the computationally intensive tasks to this device. Experiments on a hybrid platform consisting of state-of-the-art general-purpose multi-core processors and a GPU illustrate the potential of this approach

The pursuit of isotopic and molecular fire tracers in the polar atmosphere and cryosphere

We present an overview of recent multidisciplinary, multi-institutional efforts to identify and date major sources of combustion aerosol in the current and paleoatmospheres. The work was stimulated, in part, by an atmospheric particle \u27sample of opportunity\u27 collected at Summit, Greenland in August 1994, that bore the 14C imprint of biomass burning. During the summer field seasons of 1995 and 1996, we collected air filter, surface snow and snowpit samples to investigate chemical and isotopic evidence of combustion particles that had been transported from distant fires. Among the chemical tracers employed for source identification are organic acids, potassium and ammonium ions, and elemental and organic components of carbonaceous particles. Ion chromatography, performed by members of the Climate Change Research Center (University of New Hampshire), has been especially valuable in indicating periods at Summit that were likely to have been affected by the long range transport of biomass burning aerosol. Univariate and multivariate patterns of the ion concentrations in the snow and ice pinpointed surface and snowpit samples for the direct analysis of particulate (soot) carbon and carbon isotopes. The research at NIST is focusing on graphitic and polycyclic aromatic carbon, which serve as almost certain indicators of fire, and measurements of carbon isotopes, especially 14C, to distinguish fossil and biomass combustion sources. Complementing the chemical and isotopic record, are direct \u27visual\u27 (satellite imagery) records and less direct backtrajectory records, to indicate geographic source regions and transport paths. In this paper we illustrate the unique way in which the synthesis of the chemical, isotopic, satellite and trajectory data enhances our ability to develop the recent history of the formation and transport of soot deposited in the polar snow and ice

Evidence of distributed subpial T2* signal changes at 7T in multiple sclerosis : an histogram based approach

Subpial lesions are the most frequent type of cortical lesion in multiple sclerosis (MS), and are thought to be closely associated with poor clinical outcome. Neuropathological studies report that subpial lesions may come in two major types: they may appear as circumscribed, focal lesions, or extend across multiple adjacent gyri leading to a phenomenon termed “general subpial demyelination” [1]. The in vivo evaluation of diffuse subpial disease is challenging – signal changes may be subtle, and extend across large regions where signal inhomogeneities due to B1 and RF receive coil non-uniformities become more pronounced. Here, we investigate whether a histogram-based analysis of T2* signal intensity in the cortex, at 7T MRI, can show evidence of distributed subpial cortical changes in patients with MS, as described histopathologically. We hypothesized that this phenomenon would be associated with significantly increased T2* signal intensity in patients compared to age-matched controls.Center Algoritm

Constraints on the perturbed mutual motion in Didymos due to impact-induced deformation of its primary after the DART impact

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.Comment: 8 pages, 7 figures, 2 tables, accepted for publication in MNRA

Universality of rain event size distributions

We compare rain event size distributions derived from measurements in climatically different regions, which we find to be well approximated by power laws of similar exponents over broad ranges. Differences can be seen in the large-scale cutoffs of the distributions. Event duration distributions suggest that the scale-free aspects are related to the absence of characteristic scales in the meteorological mesoscale.Comment: 16 pages, 10 figure

Effects of High-Molecular-Weight Dissolved Organic Matter on Nitrogen Dynamics in the Mississippi River Plume

The dynamics of N and its interactions with labile dissolved organic C (DOC), bacteria, and phytoplankton were studied to determine potential effects of dissolved organic matter (DOM) and light on N dynamics in surface waters of the Mississippi River (USA) plume in the Gulf of Mexico. Bacterial uptake of added labeled N compounds ( 15NH4+ or 15N-labeled dissolved free amino acids. DFAA) was stimulated more by high-molecular-weight (HMW, \u3el kDa) DOM than by low-molecular-weight (LMW, \u3c l kDa) DOM. An index that inversely indicated the presence of labile DOC was defined as the fraction of assimilated Amino acid-15N that was Recovered as 15N -Ammonium (ANRA), following the additions of high-levels (4 µM) of 15N -DFAA. ANRA ratios were high in the absence of other available carbon sources because heterotrophic bacteria were forced to use the added amino acids as a carbon source for respiration rather than as a nutrient source for biomass formation. In dynamic light/dark experiments, conducted with in situ populations of organisms, uptake rates of added 15NH4+ were significantly enhanced both by the presence of light and by the addition of HMW DOM. Uptake rates of added 15N -labeled DFAA were increased by the addition of HMW DOM but not by light. ANRA ratios were consistently lower in the presence of added HMW DOM than in controls. Added HMW DOM thus appeared to stimulate the incorporation of assimilated DFAA into bacterial biomass. Bacterial growth rates were relatively high in both light and dark bottles with DFAA additions and in light bottles with HMW DOM plus NH4+ additions, but they remained comparatively low in dark bottles with added NH4+ These results are consistent with the idea that bacterial N dynamics in these euphotic waters may be tightly coupled to photosynthetic activities over short time scales

Band gap renormalization in photoexcited semiconductor quantum wire structures in the GW approximation

We investigate the dynamical self-energy corrections of the electron-hole plasma due to electron-electron and electron-phonon interactions at the band edges of a quasi-one dimensional (1D) photoexcited electron-hole plasma. The leading-order $GW$ dynamical screening approximation is used in the calculation by treating electron-electron Coulomb interaction and electron-optical phonon Fr\"{o}hlich interaction on an equal footing. We calculate the exchange-correlation induced band gap renormalization (BGR) as a function of the electron-hole plasma density and the quantum wire width. The calculated BGR shows good agreement with existing experimental results, and the BGR normalized by the effective quasi-1D excitonic Rydberg exhibits an approximate one-parameter universality.Comment: 11 pages, 3 figure

Hominids adapted to metabolize ethanol long before human-directed fermentation

Paleogenetics is an emerging field that resurrects ancestral proteins from now-extinct organisms to test, in the laboratory, models of protein function based on natural history and Darwinian evolution. Here, we resurrect digestive alcohol dehydrogenases (ADH4) from our primate ancestors to explore the history of primate-ethanol interactions. The evolving catalytic properties of these resurrected enzymes show that our ape ancestors gained a digestive dehydrogenase enzyme capable of metabolizing ethanol near the time that they began using the forest floor, about 10 million y ago. The ADH4 enzyme in our more ancient and arboreal ancestors did not efficiently oxidize ethanol. This change suggests that exposure to dietary sources of ethanol increased in hominids during the early stages of our adaptation to a terrestrial lifestyle. Because fruit collected from the forest floor is expected to contain higher concentrations of fermenting yeast and ethanol than similar fruits hanging on trees, this transition may also be the first time our ancestors were exposed to (and adapted to) substantial amounts of dietary ethanol

Strong anisotropy of superexchange in the copper-oxygen chains of La_{14-x}Ca_{x}Cu_{24}O_{41}

Electron spin resonance data of Cu^{2+} ions in La_{14-x}Ca_{x}Cu_{24}O_{41} crystals (x=9,11,12) reveal a very large width of the resonance line in the paramagnetic state. This signals an unusually strong anisotropy of ~10% of the isotropic Heisenberg superexchange in the Cu-O chains of this compound. The strong anisotropy can be explained by the specific geometry of two symmetrical 90 degree Cu-O-Cu bonds, which boosts the importance of orbital degrees of freedom. Our data show the apparent limitations of the applicability of an isotropic Heisenberg model to the low dimensional cuprates.Comment: 14 pages, 3 figures included, to be published in Phys. Rev. Let

On the Mechanism of Time--Delayed Feedback Control

The Pyragas method for controlling chaos is investigated in detail from the experimental as well as theoretical point of view. We show by an analytical stability analysis that the revolution around an unstable periodic orbit governs the success of the control scheme. Our predictions concerning the transient behaviour of the control signal are confirmed by numerical simulations and an electronic circuit experiment.Comment: 4 pages, REVTeX, 4 eps-figures included Phys. Rev. Lett., in press also available at http://athene.fkp.physik.th-darmstadt.de/public/wolfram.htm