Parallel performance and scalability experiments with the Danish Eulerian Model on the EPCC supercomputers

The Danish Eulerian Model (DEM) is a powerful air pollution model, designed to calculate the concentrations of various dangerous species over a large geographical region (e.g. Europe). It takes into account the main physical and chemical processes between these species, the actual meteorological conditions, emissions, etc.. This is a huge computational task and requires significant resources of storage and CPU time. Parallel computing is essential for the efficient practical use of the model. Some efficient parallel versions of the model were created over the past several years. A suitable parallel version of DEM by using the Message Passing Interface library (AIPI) was implemented on two powerful supercomputers of the EPCC - Edinburgh, available via the HPC-Europa programme for transnational access to research infrastructures in EC: a Sun Fire E15K and an IBM HPCx cluster. Although the implementation is in principal, the same for both supercomputers, few modifications had to be done for successful porting of the code on the IBM HPCx cluster. Performance analysis and parallel optimization was done next. Results from bench marking experiments will be presented in this paper. Another set of experiments was carried out in order to investigate the sensitivity of the model to variation of some chemical rate constants in the chemical submodel. Certain modifications of the code were necessary to be done in accordance with this task. The obtained results will be used for further sensitivity analysis Studies by using Monte Carlo simulation.

Estimating relationships between air mass origin and chemical composition

Observations of a chemical at a point in the atmosphere typically show sudden transitions between episodes of high and low concentration. Often these are associated with a rapid change in the origin of air arriving at the site. Lagrangian chemical models riding along trajectories can reproduce such transitions, but small timing errors from trajectory phase errors dramatically reduce the correlation between modeled concentrations and observations. Here the origin averaging technique is introduced to obtain maps of average concentration as a function of air mass origin for the East Atlantic Summer Experiment 1996 (EASE96, a ground-based chemistry campaign). These maps are used to construct origin averaged time series which enable comparison between a chemistry model and observations with phase errors factored out. The amount of the observed signal explained by trajectory changes can be quantified, as can the systematic model errors as a function of air mass origin. The Cambridge Tropospheric Trajectory model of Chemistry and Transport (CiTTyCAT) can account for over 70% of the observed ozone signal variance during EASE96 when phase errors are side-stepped by origin averaging. The dramatic increase in correlation (from 23% without averaging) cannot be achieved by time averaging. The success of the model is attributed to the strong relationship between changes in ozone along trajectories and their origin and its ability to simulate those changes. The model performs less well for longer-lived chemical constituents because the initial conditions 5 days before arrival are insufficiently well known.

Interaction of anthocyanins with human serum albumin: influence of pH and chemical structure on binding

The affinity of anthocyanins for human serum albumin (HSA) was determined by a fluorescence quenching method. The effects of pH and structure of anthocyanins on the binding constants were studied. The constants for binding of anthocyanins to HSA ranged from 1.08 x 10^5 M-1 to 13.16 x 10^5 M-1. A hydrophobic effect at acidic pH was shown by the relatively high positive entropy values under the conditions studied. Electrostatic interactions including hydrogen bonding contributed to the binding at pH 7.4. The effect of structure of anthocyanins on the affinity was pH dependent, particularly the effect of additional hydroxyl substituents. Hydroxyl substituents and glycosylation of anthocyanins decreased the affinity for binding to HSA at lower pH (especially pH 4), but increased the strength of binding at pH 7.4. In contrast, methylation of a hydroxyl group enhanced the binding at acidic pH, while this substitution reduced the strength of binding at pH 7.4. This paper has shown that changes in anthocyanin structure or reductions in pH, which may occur in the region of inflammatory sites, have an effect of the binding of anthocyanins to HSA.

Reconstruction and prediction of variations in the open solar magnetic flux and interplanetary conditions

Historic geomagnetic activity observations have been used to reveal centennial variations in the open solar flux and the near-Earth heliospheric conditions (the interplanetary magnetic field and the solar wind speed). The various methods are in very good agreement for the past 135 years when there were sufficient reliable magnetic observatories in operation to eliminate problems due to site-specific errors and calibration drifts. This review underlines the physical principles that allow these reconstructions to be made, as well as the details of the various algorithms employed and the results obtained. Discussion is included of: the importance of the averaging timescale; the key differences between “range” and “interdiurnal variability” geomagnetic data; the need to distinguish source field sector structure from heliospherically-imposed field structure; the importance of ensuring that regressions used are statistically robust; and uncertainty analysis. The reconstructions are exceedingly useful as they provide calibration between the in-situ spacecraft measurements from the past five decades and the millennial records of heliospheric behaviour deduced from measured abundances of cosmogenic radionuclides found in terrestrial reservoirs. Continuity of open solar flux, using sunspot number to quantify the emergence rate, is the basis of a number of models that have been very successful in reproducing the variation derived from geomagnetic activity. These models allow us to extend the reconstructions back to before the development of the magnetometer and to cover the Maunder minimum. Allied to the radionuclide data, the models are revealing much about how the Sun and heliosphere behaved outside of grand solar maxima and are providing a means of predicting how solar activity is likely to evolve now that the recent grand maximum (that had prevailed throughout the space age) has come to an end.

A review of measurement-based assessments of the aerosol direct radiative effect and forcing

Aerosols affect the Earth's energy budget directly by scattering and absorbing radiation and indirectly by acting as cloud condensation nuclei and, thereby, affecting cloud properties. However, large uncertainties exist in current estimates of aerosol forcing because of incomplete knowledge concerning the distribution and the physical and chemical properties of aerosols as well as aerosol-cloud interactions. In recent years, a great deal of effort has gone into improving measurements and datasets. It is thus feasible to shift the estimates of aerosol forcing from largely model-based to increasingly measurement-based. Our goal is to assess current observational capabilities and identify uncertainties in the aerosol direct forcing through comparisons of different methods with independent sources of uncertainties. Here we assess the aerosol optical depth (τ), direct radiative effect (DRE) by natural and anthropogenic aerosols, and direct climate forcing (DCF) by anthropogenic aerosols, focusing on satellite and ground-based measurements supplemented by global chemical transport model (CTM) simulations. The multi-spectral MODIS measures global distributions of aerosol optical depth (τ) on a daily scale, with a high accuracy of ±0.03±0.05τ over ocean. The annual average τ is about 0.14 over global ocean, of which about 21%±7% is contributed by human activities, as estimated by MODIS fine-mode fraction. The multi-angle MISR derives an annual average AOD of 0.23 over global land with an uncertainty of ~20% or ±0.05. These high-accuracy aerosol products and broadband flux measurements from CERES make it feasible to obtain observational constraints for the aerosol direct effect, especially over global the ocean. A number of measurement-based approaches estimate the clear-sky DRE (on solar radiation) at the top-of-atmosphere (TOA) to be about -5.5±0.2 Wm-2 (median ± standard error from various methods) over the global ocean. Accounting for thin cirrus contamination of the satellite derived aerosol field will reduce the TOA DRE to -5.0 Wm-2. Because of a lack of measurements of aerosol absorption and difficulty in characterizing land surface reflection, estimates of DRE over land and at the ocean surface are currently realized through a combination of satellite retrievals, surface measurements, and model simulations, and are less constrained. Over the oceans the surface DRE is estimated to be -8.8±0.7 Wm-2. Over land, an integration of satellite retrievals and model simulations derives a DRE of -4.9±0.7 Wm-2 and -11.8±1.9 Wm-2 at the TOA and surface, respectively. CTM simulations derive a wide range of DRE estimates that on average are smaller than the measurement-based DRE by about 30-40%, even after accounting for thin cirrus and cloud contamination. A number of issues remain. Current estimates of the aerosol direct effect over land are poorly constrained. Uncertainties of DRE estimates are also larger on regional scales than on a global scale and large discrepancies exist between different approaches. The characterization of aerosol absorption and vertical distribution remains challenging. The aerosol direct effect in the thermal infrared range and in cloudy conditions remains relatively unexplored and quite uncertain, because of a lack of global systematic aerosol vertical profile measurements. A coordinated research strategy needs to be developed for integration and assimilation of satellite measurements into models to constrain model simulations. Enhanced measurement capabilities in the next few years and high-level scientific cooperation will further advance our knowledge.

Physical and sensory properties of regular and reduced-fat pound cakes with added amaranth flour

The sum of wheat flour and corn starch was replaced by 10, 20, or 30% whole amaranth flour in both conventional (C) and reduced fat (RF) pound cakes. and the effects on physical and sensory properties of the cakes were investigated. RF presented 33% fat reduction. The increasing amaranth levels darkened crust and crumb of cakes, which decreased color acceptability. Fresh amaranth-containing cakes had similar texture characteristics to (he controls, evaluated both instrumentally and sensorially. Sensory evaluation revealed that replacement by 30% amaranth flour decreased C cakes overall acceptability scores, clue to its lower specific volume and darker color. Amaranth flour levels had no significant effect on overall acceptability of RF cakes. Hence, the sum of wheat flour and corn starch could be successfully replaced by up to 20% amaranth flour in C and up to 30% in RF pound cakes without negatively affecting sensory quality in fresh cakes. Moisture losses for all the cakes were similar, approximate to 1% per day during storage. After six days of storage, both C and RF amaranth-containing cakes had higher hardness and chewiness values than control cakes. Further experiments involving sensory evaluation during storage are necessary to determine the exact limit of amaranth flour replacement.

Polymethine cyanine dyes in beta-cyclodextrin solution: multiple equilibria and chemical oxidation

Absorption and fluorescence spectroscopy, electrochemical techniques, and semiempirical calculations were employed to characterize the multiple complexation equilibria between two polymethine cyanine dyes (IR-786 and Indocyanine green-ICG, 5) and beta-cyclodextrin (beta-CD, L), as well as the chemical reactivity of the complexed and uncomplexed species against the oxidizing agents hypochlorite (HC) and hydrogen peroxide (HP). IR-786 dimerization is favored with the increase in beta-CD concentration in the form of (SL)(2) complexes. In the case of ICG, free dimers (D) and SL complexes are favored. Both IR-786 and ICG react and discolor in the presence of HC and HP. For IR-786, the reaction with HP and HC proceeds with observed rate constants of 10(-3) and 0.28 s(-1) and second-order rate constants (k(2)) of similar to 10(-3) and 10(4) M(-1) s(-1), respectively. The intermediate species observed in the bleaching reactions of IR-786 and ICG were shown, by cyclic voltammetry and VIS absorption, to result from one electron oxidation. IR-786 complexed with beta-CD is protected against bleaching in the presence of HP and HC by factors of 20 and 4, respectively. This protection was not observed in ICG complexes. Superdelocalizability profile of both dyes and frontier orbital analysis indicates that beta-CD does not protect ICG from oxidation by HP or HC, whereas the 2:2 IR-786/beta-Cd complex is able to avoid the oxidation of IR-786. We concluded that the decrease in the chemical reactivity of the dyes against oxidant agents in the presence of beta-CD is due to the formation of (SL)(2) complexes. Copyright (C) 2010 John Wiley & Sons, Ltd.

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Electronic Structure and Chemical Bonding in the Lowest Electronic States of TcN

Multiconfiguration second-order perturbation theory, with the inclusion of relativistic effects and spin-orbit Coupling, was employed to investigate the nature of the ground and low-lying Lambda-S and Omega states of the TcN molecule. Spectroscopic constants, effective bond order, and potential energy curves for 13 low-lying Lambda-S states and 5 Omega states are given, The computed ground state of TcN is of Omega = 3 symmetry (R(e) = 1.605 angstrom and omega(e) = 1085 cm(-1)), originating mainly from the (3)Delta Lambda-S ground state. This result is contrasted with the nature of the ground state for other VIIB transtion-metal mononitrides, including X(3)Sigma(-) symmetry for MnN and Omega = 0(+) symmetry for ReN, derived also from a X(3)Sigma(-) state.

Peat as a natural solid-phase for copper preconcentration and determination in a multicommuted flow system coupled to flame atomic absorption spectrometry

The physical and chemical characteristics of peat were assessed through measurement of pH, percentage of organic matter, cationic exchange capacity (CEC), elemental analysis, infrared spectroscopy and quantitative analysis of metals by ICP OES. Despite the material showed to be very acid in view of the percentage of organic matter, its CEC was significant, showing potential for retention of metal ions. This characteristic was exploited by coupling a peat mini-column to a flow system based on the multicommutation approach for the in-line copper concentration prior to flame atomic absorption spectrometric determination. Cu(II) ions were adsorbed at pH 4.5 and eluted with 0.50 mol L(-1) HNO(3). The influence of chemical and hydrodynamic parameters, such as sample pH, buffer concentration, eluent type and concentration, sample flow-rate and preconcentration time were investigated. Under the optimized conditions, a linear response was observed between 16 and 100 mu g L(-1), with a detection limit estimated as 3 mu g L(-1) at the 99.7% confidence level and an enrichment factor of 16. The relative standard deviation was estimated as 3.3% (n = 20). The mini-column was used for at least 100 sampling cycles without significant variation in the analytical response. Recoveries from copper spiked to lake water or groundwater as well as concentrates used in hemodialysis were in the 97.3-111 % range. The results obtained for copper determination in these samples agreed with those achieved by graphite furnace atomic absorption spectrometry (GFAAS) at the 95% confidence level. (C) 2009 Elsevier B.V. All rights reserved.

Electronic Structure and Chemical Bonding in the Ground and Low-Lying Electronic States of Ta(2)

The electronic structure and chemical bonding of the ground and low-lying Lambda - S and Omega states of Ta(2) were investigated at the multiconfiguration second-order perturbation theory (CASSCF//CASPT2) level. The ground state of Ta(2) is computed to be a X(3)Sigma(-)(g) state (R(e) = 2.120 angstrom, omega(e) = 323 cm(-1), and D(e) = 4.65 eV), with two low-lying singlet states close to it (a(1) Sigma(+)(g) : T(e) = 409 cm(-1), R(e) = 2.131 angstrom, and omega(e) = 313 cm(-1); b(1) Gamma(g): T(e) = 1, 038 cm(-1), R(e) = 2.127 angstrom, and omega(e) = 316 cm(-1)). These electronic states are derived from the same electronic configuration: vertical bar 13 sigma(2)(g)14 sigma(2)(g)7 delta(2)(g)13 pi(4)(u)>. The effective bond order of the X(3) Sigma(-)(g) state is 4.52, which indicates that the Ta atoms are bound by a quintuple chemical bond. The a(1) Sigma(+)(g) state interacts strongly with the X(3)Sigma(-)(g) g ground state by a second-order spin-orbit interaction, giving rise to the (1)0(g)(+) (ground state) (dominated by the X(3)Sigma(-)(g) Lambda - S ground state) and (9)0(g)(+) (dominated by the a(1) Sigma(+)(g) Lambda - S state) Omega states. These results are in line with those reported for the group 5B homonuclear transition metal diatomics. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1306-1315, 2011

Rotifer abundance, biomass, and secondary production after the recovery of hydrologic connectivity between a river and two marginal lakes (São Paulo, Brazil)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Maturation of Eugenia pyriformis seeds under different hydric and thermal conditions

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Reproductive physiology, and physical and sexual development of female offspring born to diabetic dams

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Power law distribution in education: Effect of economical, teaching, and study conditions in university entrance examination

We studied the statistical distribution of student's performance, which is measured through their marks, in university entrance examination (Vestibular) of UNESP (Universidade Estadual Paulista) with respect to (i) period of study - day versus night period (ii) teaching conditions - private versus public school (iii) economical conditions - high versus low family income. We observed long ubiquitous power law tails in physical and biological sciences in all cases. The mean value increases with better study conditions followed by better teaching and economical conditions. In humanities, the distribution is close to normal distribution with very small tail. This indicates that these power law tails in science subjects axe due to the nature of the subjects themselves. Further and better study, teaching and economical conditions axe more important for physical and biological sciences in comparison to humanities at this level of study. We explain these statistical distributions through Gradually Truncated Power law distributions. We discuss the possible reason for this peculiar behavior.