Three parallel programming paradigms: Comparisons on an archetypal PDE computation

Three paradigms for distributed-memory parallel computation that free the application programmer from the details of message passing are compared for an archetypal structured scientific computation -- a nonlinear, structured-grid partial differential equation boundary value problem -- using the same algorithm on the same hardware. All of the paradigms -- parallel languages represented by the Portland Group's HPF, (semi-)automated serial-to-parallel source-to-source translation represented by CAP-Tools from the University of Greenwich, and parallel libraries represented by Argonne's PETSc -- are found to be easy to use for this problem class, and all are reasonably effective in exploiting concurrency after a short learning curve. The level of involvement required by the application programmer under any paradigm includes specification of the data partitioning, corresponding to a geometrically simple decomposition of the domain of the PDE. Programming in SPMD style for the PETSc library requires writing only the routines that discretize the PDE and its Jacobian, managing subdomain-to-processor mappings (affine global-to-local index mappings), and interfacing to library solver routines. Programming for HPF requires a complete sequential implementation of the same algorithm as a starting point, introduction of concurrency through subdomain blocking (a task similar to the index mapping), and modest experimentation with rewriting loops to elucidate to the compiler the latent concurrency. Programming with CAPTools involves feeding the same sequential implementation to the CAPTools interactive parallelization system, and guiding the source-to-source code transformation by responding to various queries about quantities knowable only at runtime. Results representative of "the state of the practice" for a scaled sequence of structured grid problems are given on three of the most important contemporary high-performance platforms: the IBM SP, the SGI Origin 2000, and the CRAYY T3E.

Application of accelerated carbonation with a combination of Na2CO3 and CO2 in cement-based solidification/stabilization of heavy metal-bearing sediment

The efficient remediation of heavy metal-bearing sediment has been one of top priorities of ecosystem protection. Cement-based solidification/stabilization (s/s) is an option for reducing the mobility of heavy metals in the sediment and the subsequent hazard for human beings and animals. This work uses sodium carbonate as an internal carbon source of accelerated carbonation and gaseous CO2 as an external carbon source to overcome deleterious effects of heavy metals on strength development and improve the effectiveness of s/s of heavy metal-bearing sediment. In addition to the compressive strength and porosity measurements, leaching tests followed the Chinese solid waste extraction procedure for leaching toxicity - sulfuric acid and nitric acid method (HJ/T299-2007), German leaching procedure (DIN38414-S4) and US toxicity characteristic leaching procedures (TCLP) have been conducted. The experimental results indicated that the solidified sediment by accelerated carbonation was capable of reaching all performance criteria for the disposal at a Portland cement dosage of 10 wt.% and a solid/water ratio of 1: 1. The concentrations of mercury and other heavy metals in the leachates were below 0.10 mg/L and 5 mg/L, respectively, complying with Chinese regulatory level (GB5085-2007). Compared to the hydration, accelerated carbonation improved the compressive strength of the solidified sediment by more than 100% and reduced leaching concentrations of heavy metals significantly. It is considered that accelerated carbonation technology with a combination of Na2CO3 and CO2 may practically apply to cement-based s/s of heavy metal-bearing sediment. (C) 2008 Elsevier B.V. All rights reserved.

Seasonal ITCZ migration dynamically controls the location of the (sub)tropical Atlantic biogeochemical divide

Inorganic nitrogen depletion restricts productivity in much of the low-latitude oceans, generating a selective advantage for diazotrophic organisms capable of fixing atmospheric dinitrogen (N2). However, the abundance and activity of diazotrophs can in turn be controlled by the availability of other potentially limiting nutrients, including phosphorus (P) and iron (Fe). Here we present high-resolution data (∼0.3°) for dissolved iron, aluminum, and inorganic phosphorus that confirm the existence of a sharp north–south biogeochemical boundary in the surface nutrient concentrations of the (sub)tropical Atlantic Ocean. Combining satellite-based precipitation data with results from a previous study, we here demonstrate that wet deposition in the region of the intertropical convergence zone acts as the major dissolved iron source to surface waters. Moreover, corresponding observations of N2 fixation and the distribution of diazotrophic Trichodesmium spp. indicate that movement in the region of elevated dissolved iron as a result of the seasonal migration of the intertropical convergence zone drives a shift in the latitudinal distribution of diazotrophy and corresponding dissolved inorganic phosphorus depletion. These conclusions are consistent with the results of an idealized numerical model of the system. The boundary between the distinct biogeochemical systems of the (sub)tropical Atlantic thus appears to be defined by the diazotrophic response to spatial–temporal variability in external Fe inputs. Consequently, in addition to demonstrating a unique seasonal cycle forced by atmospheric nutrient inputs, we suggest that the underlying biogeochemical mechanisms would likely characterize the response of oligotrophic systems to altered environmental forcing over longer timescales.

Satellite estimates of net community production indicate predominance of net autotrophy in the Atlantic Ocean

There is ongoing debate as to whether the oligotrophic ocean is predominantly net autotrophic and acts as a CO2 sink, or net heterotrophic and therefore acts as a CO2 source to the atmosphere. This quantification is challenging, both spatially and temporally, due to the sparseness of measurements. There has been a concerted effort to derive accurate estimates of phytoplankton photosynthesis and primary production from satellite data to fill these gaps; however there have been few satellite estimates of net community production (NCP). In this paper, we compare a number of empirical approaches to estimate NCP from satellite data with in vitro measurements of changes in dissolved O2 concentration at 295 stations in the N and S Atlantic Ocean (including the Antarctic), Greenland and Mediterranean Seas. Algorithms based on power laws between NCP and particulate organic carbon production (POC) derived from 14C uptake tend to overestimate NCP at negative values and underestimate at positive values. An algorithm that includes sea surface temperature (SST) in the power function of NCP and 14C POC has the lowest bias and root-mean square error compared with in vitro measured NCP and is the most accurate algorithm for the Atlantic Ocean. Nearly a 13 year time series of NCP was generated using this algorithm with SeaWiFS data to assess changes over time in different regions and in relation to climate variability. The North Atlantic subtropical and tropical Gyres (NATL) were predominantly net autotrophic from 1998 to 2010 except for boreal autumn/winter, suggesting that the northern hemisphere has remained a net sink for CO2 during this period. The South Atlantic subtropical Gyre (SATL) fluctuated from being net autotrophic in austral spring-summer, to net heterotrophic in austral autumn–winter. Recent decadal trends suggest that the SATL is becoming more of a CO2 source. Over the Atlantic basin, the percentage of satellite pixels with negative NCP was 27%, with the largest contributions from the NATL and SATL during boreal and austral autumn–winter, respectively. Variations in NCP in the northern and southern hemispheres were correlated with climate indices. Negative correlations between NCP and the multivariate ENSO index (MEI) occurred in the SATL, which explained up to 60% of the variability in NCP. Similarly there was a negative correlation between NCP and the North Atlantic Oscillation (NAO) in the Southern Sub-Tropical Convergence Zone (SSTC),which explained 90% of the variability. There were also positive correlations with NAO in the Canary Current Coastal Upwelling (CNRY) and Western Tropical Atlantic (WTRA)which explained 80% and 60% of the variability in each province, respectively. MEI and NAO seem to play a role in modifying phases of net autotrophy and heterotrophy in the Atlantic Ocean.

Conventional 3D staging PET/CT in CT simulation for lung cancer: impact of rigid and deformable target volume alignments for radiotherapy treatment planning

Objective: Positron emission tomography (PET)/CT scans can improve target definition in radiotherapy for non-small cell lung cancer (NSCLC). As staging PET/CT scans are increasingly available, we evaluated different methods for co-registration of staging PET/CT data to radiotherapy simulation (RTP) scans.

Methods: 10 patients underwent staging PET/CT followed by RTP PET/CT. On both scans, gross tumour volumes (GTVs) were delineated using CT (GTVCT) and PET display settings. Four PET-based contours (manual delineation, two threshold methods and a source-to-background ratio method) were delineated. The CT component of the staging scan was co-registered using both rigid and deformable techniques to the CT component of RTP PET/CT. Subsequently rigid registration and deformation warps were used to transfer PET and CT contours from the staging scan to the RTP scan. Dice’s similarity coefficient (DSC) was used to assess the registration accuracy of staging-based GTVs following both registration methods with the GTVs delineated on the RTP PET/CT scan.

Results: When the GTVCT delineated on the staging scan after both rigid registration and deformation was compared with the GTVCT on the RTP scan, a significant improvement in overlap (registration) using deformation was observed (mean DSC 0.66 for rigid registration and 0.82 for deformable registration, p50.008). A similar comparison for PET contours revealed no significant improvement in overlap with the use of deformable registration.

Conclusions: No consistent improvements in similarity measures were observed when deformable registration was used for transferring PET-based contours from a staging PET/CT. This suggests that currently the use of rigid registration remains the most appropriate method for RTP in NSCLC.

Computer simulating a trial of a load-bearing implant: example of an intramedullary prosthesis

Computational modelling is becoming ever more important for obtaining regulatory approval for new medical devices. An accepted approach is to infer performance in a population from an analysis conducted for an idealised or ‘average’ patient; we present here a method for predicting the performance of an orthopaedic implant when released into a population—effectively simulating a clinical trial. Specifically we hypothesise that an analysis based on a method for predicting the performance in a population will lead to different conclusions than an analysis based on an idealised or ‘average’ patient. To test this hypothesis we use a finite element model of an intramedullary implant in a bone whose size and remodelling activity is different for each individual in the population. We compare the performance of a low Young’s modulus implant (View the MathML source) to one with a higher Young’s modulus (200 GPa). Cyclic loading is applied and failure is assumed when the migration of the implant relative to the bone exceeds a threshold magnitude. The analysis for an idealised of ‘average’ patient predicts that the lower modulus device survives longer whereas the analysis simulating a clinical trial predicts no statistically-significant tendency (p=0.77) for the low modulus device to perform better. It is concluded that population-based simulations of implant performance–simulating a clinical trial–present a very valuable opportunity for more realistic computational pre-clinical testing of medical devices.

Sulphur chemistry and evolution in hot cores

We compare the results of our JCMT spectral line survey of molecular gas towards ultracompact HII regions with the predictions of models of sulphur chemistry in hot cores. We investigate the range of evolutionary models that are consistent with the observed physical conditions and chemical abundances, and see to what extent it is possible to constrain core ages by comparing abundances with the predictions of chemical models. The observed abundance ratios vary little from source to source, suggesting that all the sources are at a similar evolutionary stage. The models are capable of predicting the observed abundances of H2S, SO, SO2, and CS. The models fail to predict the amount of OCS observed, suggesting that an alternative formation route is required. An initial H2S abundance from grain mantle evaporation of similar to 10(-7) is preferred.

Modelling the cycling of persistent organic pollutants (POPs) in the North Sea system: Fluxes, loading, seasonality, trends

The fate and cycling of two selected POPs is investigated for the North Sea system with an improved version of a fate and transport ocean model (FANTOM). The model uses atmospheric data from the EMEP MSC East POP model (Gusev et al., 2009), giving reasonable concentrations and seasonal distributions for the entire region, as opposed to the three observation stations that Ilyina et al. (2006) were limited to. Other model improvements include changes in the calculation of POP exchange between the water column and sediment.

We chose to simulate the fate of two POPs with very different properties, ?-HCH and PCB 153. Since the fate and cycling of POPs are strongly affected by hydrodynamic processes, a high resolution version of the Hamburg Shelf Ocean Model (HAMSOM) was developed and utilised. Simulations were made for the period 1996–2005. Both models were validated by comparing results with available data, which showed that the simulations were of very satisfactory quality.

Model results show that the North Sea is a net sink for ?-HCH and a net source to the atmosphere of PCB 153. Total masses of ?-HCH and PCB 153 in 2005 are reduced to 30% and 50%, respectively, of 1996 values.

Storms resuspending bottom sediments into the water column mobilise POPs into the atmosphere and have the potential to deliver substantial loads of these POPs into Europe.

Food-poisoning and commercial air travel

With the introduction of budget airlines and greater competitiveness amongst all airlines, air travel has now become an extremely popular form of travel, presenting its own unique set of risks from food poisoning. Foodborne illness associated with air travel is quite uncommon in the modern era. However, when it occurs, it may have serious implications for passengers and when crew are affected, has the potential to threaten safety. Quality, safe, in-flight catering relies on high standards of food preparation and storage; this applies at the airport kitchens (or at subcontractors' facilities), on the aircraft and in the transportation vehicles which carry the food from the ground source to the aircraft. This is especially challenging in certain countries. Several foodborne outbreaks have been recorded by the airline industry as a result of a number of different failures of these systems. These have provided an opportunity to learn from past mistakes and current practice has, therefore, reached such a standard so as to minimise risk of failures of this kind. This review examines: (i) the origin of food safety in modern commercial aviation; (ii) outbreaks which have occurred previously relating to aviation travel; (iii) the microbiological quality of food and water on board commercial aircraft; and (iv) how Hazard Analysis Critical Control Points may be employed to maintain food safety in aviation travel.

Nanoporous metal organic framework materials for hydrogen storage

Hydrogen is expected to play an important role in future transportation as a promising alternative clean energy source to carbon-based fuels. One of the key challenges to commercialize hydrogen energy is to develop appropriate onboard hydrogen storage systems, capable of charging and discharging large quantities of hydrogen with fast enough kinetics to meet commercial requirements. Metal organic framework (MOF) is a new type of inorganic and organic hybrid nanoporous particulate materials. Its diverse networks can enhance hydrogen storage through tuning the structure and property of MOFs. The MOF materials so far developed adsorb hydrogen through weak dispersion interactions, which allow significant quantity of hydrogen to be stored at cryogenic temperatures with fast kinetics. Novel MOFs are being developed to strengthen the interactions between hydrogen and MOFs in order to store hydrogen under ambient conditions. This review surveys the development of such candidate materials, their performance and future research needs. (C) 2009 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

Driven optomechanical systems for mechanical entanglement distribution

We consider the distribution of entanglement from a multimode optical driving source to a network of remote and independent optomechanical systems. By focusing on the tripartite case, we analyse the effects that the features of the optical input states have on the degree and sharing structure of the distributed, fully mechanical, entanglement. This study, which is conducted looking at the mechanical steady state, highlights the structure of the entanglement distributed among the nodes and determines the relative efficiency between bipartite and tripartite entanglement transfer. We discuss a few open points, some of which are directed towards the bypassing of such limitations.

Physician recommendation for invasive prenatal testing: the case of the 'precious baby'

Do clinicians manage pregnancies conceived by assisted reproductive technologies (ART) differently from spontaneous pregnancies?

Clinicians decisions about prenatal testing during pregnancy depend, at least partially, on the method of conception.

Research thus far has shown that patients decisions regarding prenatal screening are different in ART pregnancies compared with spontaneous ones, such that ART pregnancies may be considered more valuable or precious than pregnancies conceived without treatment.

In this cross-sectional study, preformed during the year 2011, 163 obstetricians and gynecologists in Israel completed an anonymous online questionnaire.

Clinicians were randomly assigned to read one of two versions of a vignette describing the case of a pregnant woman. The two versions differed only with regard to the method of conception (ART; n 78 versus spontaneous; n 85). Clinicians were asked to provide their recommendations regarding amniocentesis.

The response rate among all clinicians invited to complete the questionnaire was 16.7. Of the 85 clinicians presented with the spontaneous pregnancy scenario, 37 (43.5) recommended amniocentesis. In contrast, of the 78 clinicians presented with the ART pregnancy scenario, only 15 (19.2) recommended the test. Clinicians were 3.2 (95 confidence interval [CI]: 1.66.6) times more likely to recommend amniocentesis for a spontaneous pregnancy than for an ART pregnancy.

The study is limited by a low response rate, the relatively small sample and the hypothetical nature of the decision, as clinician recommendations may have differed in an actual clinical setting.

Our findings show that fertility history and use of ART may affect clinicians recommendations regarding amniocentesis following receipt of screening test results. This raises the question of how subjective factors influence clinicians decisions regarding other aspects of pregnancy management.

There was no funding source to this study. The authors declare no conflicts of interest.

Defining Target Volumes for Stereotactic Ablative Radiotherapy of Early-stage Lung Tumours:A Comparison of Three-dimensional F-fluorodeoxyglucose Positron Emission Tomography and Four-dimensional Computed Tomography

Aims: High local control rates are achieved in stage I lung cancer using stereotactic ablative radiotherapy. Target delineation is commonly based on four-dimensional computed tomography (CT) scans. Target volumes defined by positron emission tomography/computed tomography (PET/CT) are compared with those defined by four-dimensional CT and conventional ('three-dimensional') F-fluorodeoxyglucose (F-FDG) PET/CT. Materials and methods: For 16 stage I non-small cell lung cancer tumours, six approaches for deriving PET target volumes were evaluated: manual contouring, standardised uptake value (SUV) absolute threshold of 2.5, 35% of maximum SUV (35%SUV), 41% of SUV (41%SUV) and two different source to background ratio techniques (SBR-1 and SBR-2). PET-derived target volumes were compared with the internal target volume (ITV) from the modified maximum intensity projection (MIP ITV). Volumetric and positional correlation was assessed using the Dice similarity coefficient (DSC). Results: PET-based target volumes did not correspond to four-dimensional CT-based target volumes. The mean DSC relative to MIP ITV were: PET manual = 0.64, SUV2.5 = 0.64, 35%SUV = 0.63, 41%SUV = 0.57. SBR-1 = 0.52, SBR-2 = 0.49. PET-based target volumes were smaller than corresponding MIP ITVs. Conclusions: Conventional three-dimensional F-FDG PET-derived target volumes for lung stereotactic ablative radiotherapy did not correspond well with those derived from four-dimensional CT, including those in routine clinical use (MIP ITV). Caution is required in using three-dimensional PET for motion encompassing target volume delineation. © 2012 The Royal College of Radiologists.

Defining target volumes for stereotactic ablative radiotherapy of early-stage lung tumours: a comparison of three-dimensional 18F-fluorodeoxyglucose positron emission tomography and four-dimensional computed tomography.:PETCT versus 4DCT for SABR in Early Stage NSCLC

AIMS: High local control rates are achieved in stage I lung cancer using
stereotactic ablative radiotherapy. Target delineation is commonly based on
four-dimensional computed tomography (CT) scans. Target volumes defined by
positron emission tomography/computed tomography (PET/CT) are compared with those defined by four-dimensional CT and conventional ('three-dimensional')
(18)F-fluorodeoxyglucose ((18)F-FDG) PET/CT.

MATERIALS AND METHODS: For 16 stage I non-small cell lung cancer tumours, six
approaches for deriving PET target volumes were evaluated: manual contouring,
standardised uptake value (SUV) absolute threshold of 2.5, 35% of maximum SUV
(35%SUV(MAX)), 41% of SUV(MAX) (41%SUV(MAX)) and two different source to
background ratio techniques (SBR-1 and SBR-2). PET-derived target volumes were compared with the internal target volume (ITV) from the modified maximum
intensity projection (MIP(MOD) ITV). Volumetric and positional correlation was
assessed using the Dice similarity coefficient (DSC).

RESULTS: PET-based target volumes did not correspond to four-dimensional CT-based target volumes. The mean DSC relative to MIP(MOD) ITV were: PET manual = 0.64, SUV2.5 = 0.64, 35%SUV(MAX) = 0.63, 41%SUV(MAX) = 0.57. SBR-1 = 0.52, SBR-2 =0.49. PET-based target volumes were smaller than corresponding MIP ITVs.

CONCLUSIONS: Conventional three-dimensional (18)F-FDG PET-derived target volumes for lung stereotactic ablative radiotherapy did not correspond well with those derived from four-dimensional CT, including those in routine clinical use
(MIP(MOD) ITV). Caution is required in using three-dimensional PET for motion
encompassing target volume delineation.

In situ electron holography of the dynamic magnetic field emanating from a hard-disk drive writer

The proliferation of mobile devices in society accessing data via the ‘cloud’ is imposing a dramatic increase in the amount of information to be stored on hard disk drives (HDD) used in servers. Forecasts are that areal densities will need to increase by as much as 35% compound per annum and by 2020 cloud storage capacity will be around 7 zettabytes corresponding to areal densities of 2 Tb/in2. This requires increased performance from the magnetic pole of the electromagnetic writer in the read/write head in the HDD. Current state-of-art writing is undertaken by morphologically complex magnetic pole of sub 100 nm dimensions, in an environment of engineered magnetic shields and it needs to deliver strong directional magnetic field to areas on the recording media around 50 nm x 13 nm. This points to the need for a method to perform direct quantitative measurements of the magnetic field generated by the write pole at the nanometer scale. Here we report on the complete in situ quantitative mapping of the magnetic field generated by a functioning write pole in operation using electron holography. Opportunistically, it points the way towards a new nanoscale magnetic field source to further develop in situ Transmission Electron Microscopy.