A numerical model for multigroup radiation hydrodynamics

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well established kinetic model through Marshak wave tests with frequency dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.Comment: 27 pages, 11 figures, Accepted for publication in JQSR

ARMC 5 Variants and Risk of Hypertension in Blacks: MH- GRID Study.

Background We recently found that ARMC 5 variants may be associated with primary aldosteronism in blacks. We investigated a cohort from the MH - GRID (Minority Health Genomics and Translational Research Bio-Repository Database) and tested the association between ARMC 5 variants and blood pressure in black s. Methods and Results Whole exome sequencing data of 1377 black s were analyzed. Target single-variant and gene-based association analyses of hypertension were performed for ARMC 5, and replicated in a subset of 3015 individuals of African descent from the UK Biobank cohort. Sixteen rare variants were significantly associated with hypertension ( P=0.0402) in the gene-based (optimized sequenced kernel association test) analysis; the 16 and one other, rs116201073, together, showed a strong association ( P=0.0003) with blood pressure in this data set. The presence of the rs116201073 variant was associated with lower blood pressure. We then used human embryonic kidney 293 and adrenocortical H295R cells transfected with an ARMC 5 construct containing rs116201073 (c.*920T>C). The latter was common in both the discovery ( MH - GRID ) and replication ( UK Biobank) data and reached statistical significance ( P=0.044 [odds ratio, 0.7] and P=0.007 [odds ratio, 0.76], respectively). The allele carrying rs116201073 increased levels of ARMC5 mRNA , consistent with its protective effect in the epidemiological data. Conclusions ARMC 5 shows an association with hypertension in black s when rare variants within the gene are considered. We also identified a protective variant of the ARMC 5 gene with an effect on ARMC 5 expression confirmed in vitro. These results extend our previous report of ARMC 5's possible involvement in the determination of blood pressure in blacks

Head and neck target delineation using a novel PET automatic segmentation algorithm

Purpose To evaluate the feasibility and impact of using a novel advanced PET auto-segmentation method in Head and Neck (H&N) radiotherapy treatment (RT) planning. Methods ATLAAS, Automatic decision Tree-based Learning Algorithm for Advanced Segmentation, previously developed and validated on pre-clinical data, was applied to 18F-FDG-PET/CT scans of 20 H&N patients undergoing Intensity Modulated Radiation Therapy. Primary Gross Tumour Volumes (GTVs) manually delineated on CT/MRI scans (GTVpCT/MRI), together with ATLAAS-generated contours (GTVpATLAAS) were used to derive the RT planning GTV (GTVpfinal). ATLAAS outlines were compared to CT/MRI and final GTVs qualitatively and quantitatively using a conformity metric. Results The ATLAAS contours were found to be reliable and useful. The volume of GTVpATLAAS was smaller than GTVpCT/MRI in 70% of the cases, with an average conformity index of 0.70. The information provided by ATLAAS was used to grow the GTVpCT/MRI in 10 cases (up to 10.6 mL) and to shrink the GTVpCT/MRI in 7 cases (up to 12.3 mL). ATLAAS provided complementary information to CT/MRI and GTVpATLAAS contributed to up to 33% of the final GTV volume across the patient cohort. Conclusions ATLAAS can deliver operator independent PET segmentation to augment clinical outlining using CT and MRI and could have utility in future clinical studies

Simulations of protostellar collapse using multigroup radiation hydrodynamics. I. The first collapse

Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larson's first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larson's first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.Comment: 9 pages, 5 figures, accepted for publication in A&

Numerical simulation of strongly nonlinear and dispersive waves using a Green-Naghdi model

We investigate here the ability of a Green-Naghdi model to reproduce strongly nonlinear and dispersive wave propagation. We test in particular the behavior of the new hybrid finite-volume and finite-difference splitting approach recently developed by the authors and collaborators on the challenging benchmark of waves propagating over a submerged bar. Such a configuration requires a model with very good dispersive properties, because of the high-order harmonics generated by topography-induced nonlinear interactions. We thus depart from the aforementioned work and choose to use a new Green-Naghdi system with improved frequency dispersion characteristics. The absence of dry areas also allows us to improve the treatment of the hyperbolic part of the equations. This leads to very satisfying results for the demanding benchmarks under consideration

Technical note: development of a 3D printed subresolution sandwich phantom for validation of brain SPECT analysis

Purpose: To make an adaptable, head shaped radionuclide phantom to simulate molecular imaging of the brain using clinical acquisition and reconstruction protocols. This will allow the characterization and correction of scanner characteristics, and improve the accuracy of clinical image analysis, including the application of databases of normal subjects. Methods: A fused deposition modeling 3D printer was used to create a head shaped phantom made up of transaxial slabs, derived from a simulated MRI dataset. The attenuation of the printed polylactide (PLA), measured by means of the Hounsfield unit on CT scanning, was set to match that of the brain by adjusting the proportion of plastic filament and air (fill ratio). Transmission measurements were made to verify the attenuation of the printed slabs. The radionuclide distribution within the phantom was created by adding 99mTc pertechnetate to the ink cartridge of a paper printer and printing images of gray and white matter anatomy, segmented from the same MRI data. The complete subresolution sandwich phantom was assembled from alternate 3D printed slabs and radioactive paper sheets, and then imaged on a dual headed gamma camera to simulate an HMPAO SPECT scan. Results: Reconstructions of phantom scans successfully used automated ellipse fitting to apply attenuation correction. This removed the variability inherent in manual application of attenuation correction and registration inherent in existing cylindrical phantom designs. The resulting images were assessed visually and by count profiles and found to be similar to those from an existing elliptical PMMA phantom. Conclusions: The authors have demonstrated the ability to create physically realistic HMPAO SPECT simulations using a novel head-shaped 3D printed subresolution sandwich method phantom. The phantom can be used to validate all neurological SPECT imaging applications. A simple modification of the phantom design to use thinner slabs would make it suitable for use in PET

The Fourth SeaWiFS HPLC Analysis Round-Robin Experiment (SeaHARRE-4)

Ten international laboratories specializing in the determination of marine pigment concentrations using high performance liquid chromatography (HPLC) were intercompared using in situ samples and a mixed pigment sample. Although prior Sea-viewing Wide Field-of-view Sensor (SeaWiFS) High Performance Liquid Chromatography (HPLC) Round-Robin Experiment (SeaHARRE) activities conducted in open-ocean waters covered a wide dynamic range in productivity, and some of the samples were collected in the coastal zone, none of the activities involved exclusively coastal samples. Consequently, SeaHARRE-4 was organized and executed as a strictly coastal activity and the field samples were collected from primarily eutrophic waters within the coastal zone of Denmark. The more restrictive perspective limited the dynamic range in chlorophyll concentration to approximately one and a half orders of magnitude (previous activities covered more than two orders of magnitude). The method intercomparisons were used for the following objectives: a) estimate the uncertainties in quantitating individual pigments and higher-order variables formed from sums and ratios; b) confirm if the chlorophyll a accuracy requirements for ocean color validation activities (approximately 25%, although 15% would allow for algorithm refinement) can be met in coastal waters; c) establish the reduction in uncertainties as a result of applying QA procedures; d) show the importance of establishing a properly defined referencing system in the computation of uncertainties; e) quantify the analytical benefits of performance metrics, and f) demonstrate the utility of a laboratory mix in understanding method performance. In addition, the remote sensing requirements for the in situ determination of total chlorophyll a were investigated to determine whether or not the average uncertainty for this measurement is being satisfied

Intercalibration of the barrel electromagnetic calorimeter of the CMS experiment at start-up

Calibration of the relative response of the individual channels of the barrel electromagnetic calorimeter of the CMS detector was accomplished, before installation, with cosmic ray muons and test beams. One fourth of the calorimeter was exposed to a beam of high energy electrons and the relative calibration of the channels, the intercalibration, was found to be reproducible to a precision of about 0.3%. Additionally, data were collected with cosmic rays for the entire ECAL barrel during the commissioning phase. By comparing the intercalibration constants obtained with the electron beam data with those from the cosmic ray data, it is demonstrated that the latter provide an intercalibration precision of 1.5% over most of the barrel ECAL. The best intercalibration precision is expected to come from the analysis of events collected in situ during the LHC operation. Using data collected with both electrons and pion beams, several aspects of the intercalibration procedures based on electrons or neutral pions were investigated

Positioning in cyberspace: evaluating telecom web sites using correspondence analysis

The growth of presence in the marketspace of the Web has been exponential, both in general and within specific industries. While the academic literature on the phenomena is still in its infancy, there exists a pressing need to establish methodologies to evaluate and map Web sites within and between industries. With the exception of a few notable papers little or no theoretical, descriptive, or normative research has been conducted into Web site positioning and evaluation. This article goes a little towards addressing this lacuna. Specifically, Web sites from the Telecom Industry are evaluated and mapped using correspondence analysis. The technique of correspondence analysis and the interpretation of the maps produced are described in detail. The implications for management are also discussed.</p

Dilepton production in heavy ion collisions at intermediate energies

We present a unified description of the vector meson and dilepton production in elementary and in heavy ion reactions. The production of vector mesons ($\rho,\omega$) is described via the excitation of nuclear resonances ($R$). The theoretical framework is an extended vector meson dominance model (eVMD). The treatment of the resonance decays $R\longmapsto NV$ with arbitrary spin is covariant and kinematically complete. The eVMD includes thereby excited vector meson states in the transition form factors. This ensures correct asymptotics and provides a unified description of photonic and mesonic decays. The resonance model is successfully applied to the $\omega$ production in $p+p$ reactions. The same model is applied to the dilepton production in elementary reactions ($p+p, p+d$). Corresponding data are well reproduced. However, when the model is applied to heavy ion reactions in the BEVALAC/SIS energy range the experimental dilepton spectra measured by the DLS Collaboration are significantly underestimated at small invariant masses. As a possible solution of this problem the destruction of quantum interference in a dense medium is discussed. A decoherent emission through vector mesons decays enhances the corresponding dilepton yield in heavy ion reactions. In the vicinity of the $\rho/\omega$-peak the reproduction of the data requires further a substantial collisional broadening of the $\rho$ and in particular of the $\omega$ meson.Comment: 32 pages revtex, 19 figures, to appear in PR