Precise analysis of array usage in scientific programs

The automatic transformation of sequential programs for efficient execution on parallel computers involves a number of analyses and restructurings of the input. Some of these analyses are based on computing array sections, a compact description of a range of array elements. Array sections describe the set of array elements that are either read or written by program statements. These sections can be compactly represented using shape descriptors such as regular sections, simple sections, or generalized convex regions. However, binary operations such as Union performed on these representations do not satisfy a straightforward closure property, e.g., if the operands to Union are convex, the result may be nonconvex. Approximations are resorted to in order to satisfy this closure property. These approximations introduce imprecision in the analyses and, furthermore, the imprecisions resulting from successive operations have a cumulative effect. Delayed merging is a technique suggested and used in some of the existing analyses to minimize the effects of approximation. However, this technique does not guarantee an exact solution in a general setting. This article presents a generalized technique to precisely compute Union which can overcome these imprecisions.

Spatial and temporal variability of the concentration field from localized releases in a regular building array

Spatial and temporal fluctuations in the concentration field from an ensemble of continuous point-source releases in a regular building array are analyzed from data generated by direct numerical simulations. The release is of a passive scalar under conditions of neutral stability. Results are related to the underlying flow structure by contrasting data for an imposed wind direction of 0 deg and 45 deg relative to the buildings. Furthermore, the effects of distance from the source and vicinity to the plume centreline on the spatial and temporal variability are documented. The general picture that emerges is that this particular geometry splits the flow domain into segments (e.g. “streets” and “intersections”) in each of which the air is, to a first approximation, well mixed. Notable exceptions to this general rule include regions close to the source, near the plume edge, and in unobstructed channels when the flow is aligned. In the oblique (45 deg) case the strongly three-dimensional nature of the flow enhances mixing of a scalar within the canopy leading to reduced temporal and spatial concentration fluctuations within the plume core. These fluctuations are in general larger for the parallel flow (0 deg) case, especially so in the long unobstructed channels. Due to the more complex flow structure in the canyon-type streets behind buildings, fluctuations are lower than in the open channels, though still substantially larger than for oblique flow. These results are relevant to the formulation of simple models for dispersion in urban areas and to the quantification of the uncertainties in their predictions.

Using NMR-based metabolomics to evaluate postprandial urinary responses following consumption of minimally processed wheat bran or wheat aleurone by men and women

Wheat bran, and especially wheat aleurone fraction, are concentrated sources of a wide range of components which may contribute to the health benefits associated with higher consumption of whole-grain foods. This study used NMR metabolomics to evaluate urine samples from baseline at one and two hours postprandially, following the consumption of minimally processed bran, aleurone or control by 14 participants (7 Females; 7 Males) in a randomized crossover trial. The methodology discriminated between the urinary responses of control, and bran and aleurone, but not between the two fractions. Compared to control, consumption of aleurone or bran led to significantly and substantially higher urinary concentrations of lactate, alanine, N-acetylaspartate acid and N-acetylaspartylglutamate and significantly and substantially lower urinary betaine concentrations at one and two hours postprandially. There were sex related differences in urinary metabolite profiles with generally higher hippurate and citrate and lower betaine in females compared to males. Overall, this postprandial study suggests that acute consumption of bran or aleurone is associated with a number of physiological effects that may impact on energy metabolism and which are consistent with longer term human and animal metabolomic studies that used whole-grain wheat diets or wheat fractions.

Internal nanoparticle structure of temperature-responsive self-assembled PNIPAM-b-PEG-b-PNIPAM triblock copolymers in aqueous solutions: NMR, SANS and light scattering studies

In this study we report detailed information on the internal structure of PNIPAM-b-PEG-b-PNIPAM nanoparticles formed from self-assembly in aqueous solutions upon increase in temperature. NMR spectroscopy, light scattering and small-angle neutron scattering (SANS) were used to monitor different stages of nanoparticle formation as a function of temperature, providing insight into the fundamental processes involved. The presence of PEG in a copolymer structure significantly affects the formation of nanoparticles, making their transition to occur over a broader temperature range. The crucial parameter that controls the transition is the ratio of PEG/PNIPAM. For pure PNIPAM, the transition is sharp; the higher the PEG/PNIPAM ratio results in a broader transition. This behavior is explained by different mechanisms of PNIPAM block incorporation during nanoparticle formation at different PEG/PNIPAM ratios. Contrast variation experiments using SANS show that the structure of nanoparticles above cloud point temperatures for PNIPAM-b-PEG-b-PNIPAM copolymers is drastically different from the structure of PNIPAM mesoglobules. In contrast with pure PNIPAM mesoglobules, where solid-like particles and chain network with a mesh size of 1-3 nm are present; nanoparticles formed from PNIPAM-b-PEG-b-PNIPAM copolymers have non-uniform structure with “frozen” areas interconnected by single chains in Gaussian conformation. SANS data with deuterated “invisible” PEG blocks imply that PEG is uniformly distributed inside of a nanoparticle. It is kinetically flexible PEG blocks which affect the nanoparticle formation by prevention of PNIPAM microphase separation.

Determination of Very Slow Diffusion of Paramagnetic Ions by NMR Spectroscopy

In this paper, we propose a new method of measuring the very slow paramagnetic ion diffusion coefficient using a commercial high-resolution spectrometer. If there are distinct paramagnetic ions influencing the hydrogen nuclear magnetic relaxation time differently, their diffusion coefficients can be measured separately. A cylindrical phantom filled with Fricke xylenol gel solution and irradiated with gamma rays was used to validate the method. The Fricke xylenol gel solution was prepared with 270 Bloom porcine gelatin, the phantom was irradiated with gamma rays originated from a (60)Co source and a high-resolution 200 MHz nuclear magnetic resonance (NMR) spectrometer was used to obtain the phantom (1)H profile in the presence of a linear magnetic field gradient. By observing the temporal evolution of the phantom NMR profile, an apparent ferric ion diffusion coefficient of 0.50 mu m(2)/ms due to ferric ions diffusion was obtained. In any medical process where the ionizing radiation is used, the dose planning and the dose delivery are the key elements for the patient safety and success of treatment. These points become even more important in modern conformal radio therapy techniques, such as stereotactic radiosurgery, where the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Several methods have been proposed to obtain the three-dimensional (3-D) dose distribution. Recently, we proposed an alternative method for the 3-D radiation dose mapping, where the ionizing radiation modifies the local relative concentration of Fe(2+)/Fe(3+) in a phantom containing Fricke gel and this variation is associated to the MR image intensity. The smearing of the intensity gradient is proportional to the diffusion coefficient of the Fe(3+) and Fe(2+) in the phantom. There are several methods for measurement of the ionic diffusion using NMR, however, they are applicable when the diffusion is not very slow.

Thermal Phase Behavior of DMPG Bilayers in Aqueous Dispersions as Revealed by (2)H- and (31)P-NMR

The synthetic lipid 1,2-dimyristoyl-sn-3-phosphoglycerol (DMPG), when dispersed in water/NaCl exhibits a complex phase behavior caused by its almost unlimited swelling in excess water. Using deuterium ((2)H)- and phosphorus ((31)P)-NMR we have studied the molecular properties of DMPG/water/NaCl dispersions as a function of lipid and NaCl concentration. We have measured the order profile of the hydrophobic part of the lipid bilayer with deuterated DMPG while the orientation of the phosphoglycerol headgroup was deduced from the (31)P NMR chemical shielding anisotropy. At temperatures > 30 degrees C we observe well-resolved (2)H- and (31)P NMR spectra not much different from other liquid crystalline bilayers. From the order profiles it is possible to deduce the average length of the flexible fatty acyl chain. Unusual spectra are obtained in the temperature interval of 20-25 degrees C, indicating one or several phase transitions. The most dramatic changes are seen at low lipid concentration and low ionic strength. Under these conditions and at 25 degrees C, the phosphoglycerol headgroup rotates into the hydrocarbon layer and the hydrocarbon chains show larger flexing motions than at higher temperatures. The orientation of the phosphoglycerol headgroup depends on the bilayer surface charge and correlates with the degree of dissociation of DMPG-Na(+). The larger the negative surface charge, the more the headgroup rotates toward the nonpolar region.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Trigger and aperture of the surface detector array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above 3 x 10(18) eV, for all zenith angles between 0 degrees and 60 degrees, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance. (C) 2009 Elsevier B.V. All rights reserved.

Investigation by Combined Solid-State NMR and SAXS Methods of the Morphology and Domain Size in Polystyrene-b-Polyethylene Oxide-b-Polystyrene Triblock Copolymers

The microphase structure of a series of polystyrene-b-polyethylene oxide-b-polystyrene (SEOS) triblock copolymers with different compositions and molecular weights has been studied by solid-state NMR, DSC, wide and small angle X-ray scattering (WAXS and SAXS). WAXS and DSC measurements were used to detect the presence of crystalline domains of polyethyleneoxide (PEO) blocks at room temperature as a function of the copolymer chemical composition. Furthermore, DSC experiments allowed the determination of the melting temperatures of the crystalline part of the PEO blocks. SAXS measurements, performed above and below the melting temperature of the PEO blocks, revealed the formation of periodic structures, but the absence or the weakness of high order reflections peaks did not allow a clear assessment of the morphological structure of the copolymers. This information was inferred by combining the results obtained by SAXS and (1)H NMR spin diffusion experiments, which also provided an estimation of the size of the dispersed phases of the nanostructured copolymers. (C) 2009 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 48:55-64,2010

Structural Characterization of Rare-Earth Doped Yttrium Aluminoborate Laser Glasses Using Solid State NMR

The structure of laser glasses in the system (B(2)O(3))(0.6){(Al(2)O(3))(0.4-x)(Y(2)O(3))(x)} (0.1 <= x <= 0.25) has been investigated by means of (11)B, (27)Al, and (89)Y solid state NMR as well as Y-3d core-level X-ray photoelectron spectroscopy, (11)B magic-angle spinning (MAS) NMR spectra reveal that the majority of the boron atoms are three-coordinated, and a slight increase of four-coordinated boron content with increasing x can be noticed. (27)Al MAS NMR spectra show that the alumina species are present in the coordination states four, five and six. All of them are in intimate contact with both the three- and the four-coordinate boron species and vice versa, as indicated by (11)B/(27)Al rotational echo double resonance (REDOR) data. These results are consistent with the formation of a homogeneous, nonsegregated glass structure. For the first time, (89)Y solid state NMR has been used to probe the local environment of Y(3+) ions in a glass-forming system. The intrinsic sensitivity problem associated with (89)Y NMR has been overcome by combining the benefits of paramagnetic doping with those of signal accumulation via Carr-Purcell spin echo trains. Both the (89)Y chemical shifts and the Y-3d core level binding energies are found to be rather sensitive to the yttrium bonding state and reveal that the bonding properties of the yttrium atoms in these glasses are similar to those found in the model compounds YBO(3) and YAl(3)(BO(3))(4), Based on charge balance considerations as well as (11)B NMR line shape analyses, the dominant borate species are concluded to be meta- and pyroborate anions.

Direct Observation of Millisecond to Second Motions in Proteins by Dipolar CODEX NMR Spectroscopy

We present a site-resolved study of stow (ms to s) motions in a protein in the solid (microcrystalline) state performed with the use of a modified version of the centerband-only detection of exchange (CODEX) NMR experiment. CODEX was originally based on measuring changes in molecular orientation by means of the chemical shift anisotropy (CSA) tensor, and in our modification, angular reorientations of internuclear vectors are observed. The experiment was applied to the study of stow (15)N-(1)H motions of the SH3 domain of chicken a-spectrin. The protein was perdeuterated with partial back-exchange of protons at labile sites. This allowed indirect (proton) detection of (15)N nuclei and thus a significant enhancement of sensitivity. The diluted proton system also made negligible proton-driven spin diffusion between (15)N nuclei, which interferes with the molecular exchange (motion) and hampers the acquisition of dynamic parameters. The experiment has shown that approximately half of the peaks in the 2D (15)N-(1)H correlation spectrum exhibit exchange in a different extent. The correlation time of the slow motion for most peaks is 1 to 3 s. This is the first NMR study of the internal dynamics of proteins in the solid state on the millisecond to second time scale with site-specific spectral resolution that provides both time-scale and geometry information about molecular motions.

Temperature dependence of molecular dynamics and supramolecular aggregation in MEH-PPV films: A solid-state NMR, X-ray and fluorescence spectroscopy study

This article presents an investigation of the temperature induced modification in the microstructure and dynamics of poly[2-methoxy-5-(2`-ethylhexyloxy)-1,4-phenylenevinylene] (MEH-PPV) cast films using Wide-Angle X-ray Scattering (WAXS), solid-state Nuclear Magnetic Resonance (NMR), and Fluorescence Spectroscopy (PL). MEH-PPV chain motions were characterized as a function of temperature by NMR. The results indicated that the solvent used to cast the films influences the activation energy of the side-chain motions. This was concluded from the comparison of the activation energy of the toluene cast film, E(a) = (54 +/- 8) kJ/mol, and chloroform cast film, E(a) = (69 +/- 5) kJ/mol, and could be attributed to the higher side-chain packing provided by chloroform, that preferentially solvates the side chain in contrast to toluene that solvates mainly the backbone. Concerning the backbone mobility, it was observed that the torsional motions in the MEH-PPV have average amplitude of similar to 10 degrees at 300 K, which was found to be independent of the solvent used to cast the films. In order to correlate the molecular dynamics processes with the changes in the microstructure of the polymer, in situ WAXS experiments as a function of temperature were performed and revealed that the interchain spacing in the MEH-PPV molecular aggregates increases as a function of temperature, particularly at temperatures where molecular relaxations occur. It was also observed that the WAXS peak associated with the bilayer spacing becomes narrower and its intensity increases whereas the peak associated with the inter-backbone planes reduces its intensity for higher temperatures. This last result Could be interpreted as a decrease in the number of aggregates and the reduction of the interchain species during the MEH-PPV relaxation processes. These WAXS results were correlated with PL spectra modifications observed upon temperature treatments. (C) 2008 Elsevier Ltd. All rights reserved.

Solid state NMR as a new approach for the structural characterization of rare-earth doped lead lanthanum zirconate titanate laser ceramics

To facilitate the design of laser host materials with optimized emission properties, detailed structural information at the atomic level is essential, regarding the local bonding environment of the active ions (distribution over distinct lattice sites) and their extent of local clustering as well as their population distribution over separate micro- or nanophases. The present study explores the potential of solid state NMR spectroscopy to provide such understanding for rare-earth doped lead lanthanum zirconate titanate (PLZT) ceramics. As the NMR signals of the paramagnetic dopant species cannot be observed directly, two complementary approaches are utilized: (1) direct observation of diamagnetic mimics using Sc-45 NMR and (2) study of the paramagnetic interaction of the constituent host lattice nuclei with the rare-earth dopant, using Pb-207 NMR lineshape analysis. Sc-45 MAS NMR spectra of scandium-doped PLZT samples unambiguously reveal scandium to be six-coordinated, suggesting that this rare-earth ion substitutes in the B site. Static Pb-207 spin echo NMR spectra of a series of Tm-doped PLZT samples reveal a clear influence of paramagnetic rare-earth dopant concentration on the NMR lineshape. In the latter case high-fidelity spectra can be obtained by spin echo mapping under systematic incrementation of the excitation frequency, benefiting from the signal-to-noise enhancement afforded by spin echo train Fourier transforms. Consistent with XRD data, the Pb-207 NMR lineshape analysis suggests that statistical incorporation into the PLZT lattice occurs at dopant levels of up to 1 wt.% Tm3+, while at higher levels the solubility limit is reached. (C) 2008 Elsevier Masson SAS. All rights reserved.