Absolute value preconditioning for symmetric indefinite linear systems

We introduce a novel strategy for constructing symmetric positive definite (SPD) preconditioners for linear systems with symmetric indefinite matrices. The strategy, called absolute value preconditioning, is motivated by the observation that the preconditioned minimal residual method with the inverse of the absolute value of the matrix as a preconditioner converges to the exact solution of the system in at most two steps. Neither the exact absolute value of the matrix nor its exact inverse are computationally feasible to construct in general. However, we provide a practical example of an SPD preconditioner that is based on the suggested approach. In this example we consider a model problem with a shifted discrete negative Laplacian, and suggest a geometric multigrid (MG) preconditioner, where the inverse of the matrix absolute value appears only on the coarse grid, while operations on finer grids are based on the Laplacian. Our numerical tests demonstrate practical effectiveness of the new MG preconditioner, which leads to a robust iterative scheme with minimalist memory requirements

The non-uniform, dynamic atmosphere of Betelgeuse observed at mid-infrared wavelengths

We present an interferometric study of the continuum surface of the red supergiant star Betelgeuse at 11.15 microns wavelength, using data obtained with the Berkeley Infrared Spatial Interferometer each year between 2006 and 2010. These data allow an investigation of an optically thick layer within 1.4 stellar radii of the photosphere. The layer has an optical depth of ~1 at 11.15 microns, and varies in temperature between 1900 K and 2800 K and in outer radius between 1.16 and 1.36 stellar radii. Electron-hydrogen atom collisions contribute significantly to the opacity of the layer. The layer has a non-uniform intensity distribution that changes between observing epochs. These results indicate that large-scale surface convective activity strongly influences the dynamics of the inner atmosphere of Betelgeuse, and mass-loss processes.Comment: 13 pages, 5 figures, in press (ApJ

Building A High Performance Parallel File System Using Grid Datafarm and ROOT I/O

Sheer amount of petabyte scale data foreseen in the LHC experiments require a careful consideration of the persistency design and the system design in the world-wide distributed computing. Event parallelism of the HENP data analysis enables us to take maximum advantage of the high performance cluster computing and networking when we keep the parallelism both in the data processing phase, in the data management phase, and in the data transfer phase. A modular architecture of FADS/ Goofy, a versatile detector simulation framework for Geant4, enables an easy choice of plug-in facilities for persistency technologies such as Objectivity/DB and ROOT I/O. The framework is designed to work naturally with the parallel file system of Grid Datafarm (Gfarm). FADS/Goofy is proven to generate 10^6 Geant4-simulated Atlas Mockup events using a 512 CPU PC cluster. The data in ROOT I/O files is replicated using Gfarm file system. The histogram information is collected from the distributed ROOT files. During the data replication it has been demonstrated to achieve more than 2.3 Gbps data transfer rate between the PC clusters over seven participating PC clusters in the United States and in Japan.Comment: Talk from the 2003 Computing in High Energy and Nuclear Physics (CHEP03), La Jolla, Ca, USA, March 2003, 4 pages, PDF. PSN TUDT01

Imaging the dynamical atmosphere of the red supergiant Betelgeuse in the CO first overtone lines with VLTI/AMBER

We present the first 1-D aperture synthesis imaging of the red supergiant Betelgeuse in the individual CO first overtone lines with VLTI/AMBER. The reconstructed 1-D projection images reveal that the star appears differently in the blue wing, line center, and red wing of the individual CO lines. The 1-D projection images in the blue wing and line center show a pronounced, asymmetrically extended component up to ~1.3 stellar radii, while those in the red wing do not show such a component. The observed 1-D projection images in the lines can be reasonably explained by a model in which the CO gas within a region more than half as large as the stellar size is moving slightly outward with 0--5 km s^-1, while the gas in the remaining region is infalling fast with 20--30 km s^-1. A comparison between the CO line AMBER data taken in 2008 and 2009 shows a significant time variation in the dynamics of the CO line-forming region in the photosphere and the outer atmosphere. In contrast to the line data, the reconstructed 1-D projection images in the continuum show only a slight deviation from a uniform disk or limb-darkened disk. We derive a uniform-disk diameter of 42.05 +/- 0.05 mas and a power-law-type limb-darkened disk diameter of 42.49 +/- 0.06 mas and a limb-darkening parameter of (9.7 +/- 0.5) x 10^{-2}. This latter angular diameter leads to an effective temperature of 3690 +/- 54 K for the continuum-forming layer. These diameters confirm that the near-IR size of Betelgeuse was nearly constant over the last 18 years, in marked contrast to the recently reported noticeable decrease in the mid-IR size. The continuum data taken in 2008 and 2009 reveal no or only marginal time variations, much smaller than the maximum variation predicted by the current 3-D convection simulations.Comment: 21 pages, 12 figures, accepted for publication in Astronomy and Astrophysic

De Novo Growth Zone Formation from Fission Yeast Spheroplasts

Eukaryotic cells often form polarized growth zones in response to internal or external cues. To understand the establishment of growth zones with specific dimensions we used fission yeast, which grows as a rod-shaped cell of near-constant width from growth zones located at the cell tips. Removing the cell wall creates a round spheroplast with a disorganized cytoskeleton and depolarized growth proteins. As spheroplasts recover, new growth zones form that resemble normal growing cell tips in shape and width, and polarized growth resumes. Regulators of the GTPase Cdc42, which control width in exponentially growing cells, also control spheroplast growth zone width. During recovery the Cdc42 scaffold Scd2 forms a polarized patch in the rounded spheroplast, demonstrating that a growth zone protein can organize independent of cell shape. Rga4, a Cdc42 GTPase activating protein (GAP) that is excluded from cell tips, is initially distributed throughout the spheroplast membrane, but is excluded from the growth zone after a stable patch of Scd2 forms. These results provide evidence that growth zones with normal width and protein localization can form de novo through sequential organization of cellular domains, and that the size of these growth zones is genetically controlled, independent of preexisting cell shape

Identification of Genes Affecting the Toxicity of Anti-Cancer Drug Bortezomib by Genome-Wide Screening in S. pombe

Bortezomib/PS-341/Velcade, a proteasome inhibitor, is widely used to treat multiple myeloma. While several mechanisms of the cytotoxicity of the drug were proposed, the actual mechanism remains elusive. We aimed to identify genes affecting the cytotoxicity of Bortezomib in the fission yeast S.pombe as the drug inhibits this organism's cell division cycle like proteasome mutants. Among the 2815 genes screened (covering 56% of total ORFs), 19 genes, whose deletions induce strong synthetic lethality with Bortezomib, were identified. The products of the 19 genes included four ubiquitin enzymes and one nuclear proteasome factor, and 13 of them are conserved in humans. Our results will provide useful information for understanding the actions of Bortezomib within cells

Trends in incidence and mortality of tuberculosis in Japan : a population-based study, 1997–2016

Japan is still a medium-burden tuberculosis (TB) country. We aimed to examine trends in newly notified active TB incidence and TB-related mortality in the last two decades in Japan. This is a population-based study using Japanese Vital Statistics and Japan Tuberculosis Surveillance from 1997 to 2016. We determined active TB incidence and mortality rates (per 100 000 population) by sex, age and disease categories. Joinpoint regression was applied to calculate the annual percentage change (APC) in age-adjusted mortality rates and to identify the years showing significant trend changes. Crude and age-adjusted incidence rates reduced from 33.9 to 13.9 and 37.3 to 11.3 per 100 000 population, respectively. Also, crude and age-adjusted mortality rates reduced from 2.2 to 1.5 and 2.8 to 1.0 per 100 000 population, respectively. Average APC in the incidence and mortality rates showed significant decline both in men (−6.2% and −5.4%, respectively) and women (−5.7% and −4.6%, respectively). Age-specific analysis demonstrated decreases in incidence and mortality rates for every age category, except for the incidence trend in the younger population. Although trends in active TB incidence and mortality rates in Japan have favourably decreased, the rate of decline is far from achieving TB elimination by 2035

The Arctic predictability and prediction on seasonal-to-interannual timescales (APPOSITE) data set version 1

This is the final version of the article. Available from the publisher via the DOI in this record. Discussion paper (published on 15 Oct 2015)Recent decades have seen significant developments in seasonal-to-interannual timescale climate prediction capabilities. However, until recently the potential of such systems to predict Arctic climate had not been assessed. This paper describes a multi- 5 model predictability experiment which was run as part of the Arctic Predictability and Prediction On Seasonal to Inter-annual Timescales (APPOSITE) project. The main goal of APPOSITE was to quantify the timescales on which Arctic climate is predictable. In order to achieve this, a coordinated set of idealised initial-value predictability experiments, with seven general circulation models, was conducted. This was the first model 10 intercomparison project designed to quantify the predictability of Arctic climate on seasonal to inter-annual timescales. Here we present a description of the archived data set (which is available at the British Atmospheric Data Centre) and an update of the project's results. Although designed to address Arctic predictability, this data set could also be used to assess the predictability of other regions and modes of climate vari15 ability on these timescales, such as the El Niño Southern Oscillation.This work was supported by the Natural Environment Research Council (grant NE/I029447/1). Helge Goessling was supported by a fellowship of the German Research Foundation (DFG grant GO 2464/1-1). Data storage and processing capacity was kindly provided by the British Atmospheric Data Centre (BADC). Thanks to Yanjun Jiao (CCCma) for his assistance with the CanCM4 simulations and to Bill Merryfield for his comments on a draft of the pape

The coordination of cell growth during fission yeast mating requires Ras1-GTP hydrolysis

The spatial and temporal control of polarity is fundamental to the survival of all organisms. Cells define their polarity using highly conserved mechanisms that frequently rely upon the action of small GTPases, such as Ras and Cdc42. Schizosaccharomyces pombe is an ideal system with which to study the control of cell polarity since it grows from defined tips using Cdc42-mediated actin remodeling. Here we have investigated the importance of Ras1-GTPase activity for the coordination of polarized cell growth during fission yeast mating. Following pheromone stimulation, Ras1 regulates both a MAPK cascade and the activity of Cdc42 to enable uni-directional cell growth towards a potential mating partner. Like all GTPases, when bound to GTP, Ras1 adopts an active conformation returning to an inactive state upon GTP-hydrolysis, a process accelerated through interaction with negative regulators such as GAPs. Here we show that, at low levels of pheromone stimulation, loss of negative regulation of Ras1 increases signal transduction via the MAPK cascade. However, at the higher concentrations observed during mating, hyperactive Ras1 mutations promote cell death. We demonstrate that these cells die due to their failure to coordinate active Cdc42 into a single growth zone resulting in disorganized actin deposition and unsustainable elongation from multiple tips. These results provide a striking demonstration that the deactivation stage of Ras signaling is fundamentally important in modulating cell polarity

Parameter identification problems in the modelling of cell motility

We present a novel parameter identification algorithm for the estimation of parameters in models of cell motility using imaging data of migrating cells. Two alternative formulations of the objective functional that measures the difference between the computed and observed data are proposed and the parameter identification problem is formulated as a minimisation problem of nonlinear least squares type. A Levenberg–Marquardt based optimisation method is applied to the solution of the minimisation problem and the details of the implementation are discussed. A number of numerical experiments are presented which illustrate the robustness of the algorithm to parameter identification in the presence of large deformations and noisy data and parameter identification in three dimensional models of cell motility. An application to experimental data is also presented in which we seek to identify parameters in a model for the monopolar growth of fission yeast cells using experimental imaging data. Our numerical tests allow us to compare the method with the two different formulations of the objective functional and we conclude that the results with both objective functionals seem to agree