Microbial catabolic activities are naturally selected by metabolic energy harvest rate

The fundamental trade-off between yield and rate of energy harvest per unit of substrate has been largely discussed as a main characteristic for microbial established cooperation or competition. In this study, this point is addressed by developing a generalized model that simulates competition between existing and not experimentally reported microbial catabolic activities defined only based on well-known biochemical pathways. No specific microbial physiological adaptations are considered, growth yield is calculated coupled to catabolism energetics and a common maximum biomass-specific catabolism rate (expressed as electron transfer rate) is assumed for all microbial groups. Under this approach, successful microbial metabolisms are predicted in line with experimental observations under the hypothesis of maximum energy harvest rate. Two microbial ecosystems, typically found in wastewater treatment plants, are simulated, namely: (i) the anaerobic fermentation of glucose and (ii) the oxidation and reduction of nitrogen under aerobic autotrophic (nitrification) and anoxic heterotrophic and autotrophic (denitrification) conditions. The experimentally observed cross feeding in glucose fermentation, through multiple intermediate fermentation pathways, towards ultimately methane and carbon dioxide is predicted. Analogously, two-stage nitrification (by ammonium and nitrite oxidizers) is predicted as prevailing over nitrification in one stage. Conversely, denitrification is predicted in one stage (by denitrifiers) as well as anammox (anaerobic ammonium oxidation). The model results suggest that these observations are a direct consequence of the different energy yields per electron transferred at the different steps of the pathways. Overall, our results theoretically support the hypothesis that successful microbial catabolic activities are selected by an overall maximum energy harvest rate

On identifiability for chemical systems from measurable variables

The final publication is available at Springer via http://dx.doi.org/10.1007/s10910-013-0149-4The dynamics of the composition of chemical species in reacting systems can be characterized by a set of autonomous differential equations derived from mass conservation principles and some elementary hypothesis related to chemical reactivity. These sets of ordinary differential equations are basically non-linear, their complexity grows as much increases the number of substances present in the reacting media an can be characterized by a set of phenomenological constants which contains all the relevant information about the physical system. The determination of these kinetic constants is critical for the design or control of chemical systems from a technological point of view but the non-linear nature of the equations implies that there are hidden correlations between the parameters which maybe can be revealed with a identifiability analysis.This work has been partially supported by MTM2010-18228.Cantó Colomina, B.; Coll, C.; Sánchez, E.; Cardona Navarrete, SC.; Navarro-Laboulais, J. (2014). On identifiability for chemical systems from measurable variables. Journal of Mathematical Chemistry. 52(4):1023-1035. https://doi.org/10.1007/s10910-013-0149-4S10231035524M.J. Almendral, A. Alonso, M.S. Fuentes, Development of new methodologies for on-line determination of the bromate. J. Environ. Monit. 11, 1381–1388 (2009)A. Ben-Zvi, P.J. McLellan, K.B. McAuley, Identifiability of linear time-invariant differential-algebraic systems. I. The generalized Markov parameter approach. Ind. Eng. Chem. Res. 42, 6607–6618 (2003)T.P. Bonacquisti, A drinking water utility’s perspective on bromide, bromate, and ozonation. Toxicology 221, 145–148 (2006)R. Butler, A. Godley, L. Lytton, E. Cartmell, Bromate environmental contamination: review of impact and possible treatment. Crit. Rev. Environ. Sci. Tech. 35, 193–217 (2005)R. Butler, L. Lytton, A.R. Godley, I.E. Tothill, E. Cartmell, Bromate analysis in groundwater and wastewater samples. J. Environ. Monit. 7, 999–1006 (2005)B. Cantó, S.C. Cardona, C. Coll, J. Navarro-Laboulais, E. Sánchez, Dynamic optimization of a gas-liquid reactor. J. Math. Chem. 50, 381–393 (2012)B. Cantó, C. Coll and E. Sánchez, Identifiability of a class of discretized linear partial differential algebraic equations, Math. Problems Eng. 2011, 1–12 (2011)A. Constantinides, N. Mostoufi, Numerical Methods for Chemical Engineers with MATLAB Applications, Alkis Constantinides and Navid Mostoufi, Upper Saddle River (Prentice Hall, New Jersey, 1999)P. Englezos, N. Kalogerakis, Applied Parameter Estimation for Chemical Engineers (Marcel Dekker, New York, 2001)U. von Gunten, Ozonation of drinking water. Part II. Disinfection and by-product formation in presence of bromide, iodide or chlorine. Water Res. 37, 1469–1487 (2003)B. Legube, B. Parinet, K. Gelinet, F. Berne, J-Ph Croue, Modeling of bromate formation by ozonation of surface waters in drinking water treatment. Water Res. 38, 2185–2195 (2004)Q. Liu, L.M. Schurter, C.E. Muller, S. Aloisio, J.S. Francisco, D.W. Margerum, Kinetics and mechanisms of aqueous ozone reactions with bromide, sulfite, hydrogen sulfite, iodide, and nitrite ions. Inorg. Chem. 40, 4436–4442 (2001)J.B. Rawling, J.G. Ekerdt, Chemical Reactor Analysis and Design Fundamentals (Nob Hill Pub, Madison, 2002)W.E. Stewart, M. Caracotsios, Computer Aided Modelling of Reactive Systems (John Wiley and Sons, New York, 2008)P. Westerhoff, R. Song, G. Amy, R. Minear, Numerical kinetic models for bromide oxidation to bromine and bromate. Water Res. 32, 1687–1699 (1998)World Health Organization, Bromate in Drinking-water, Document WHO/SDE/WSH/05.08/78, http://www.who.int/water_sanitation_health/dwq/chemicals/en/ (accesed 26/07/12

Constraining the pˉ/p\bar{p}/p Ratio in TeV Cosmic Rays with Observations of the Moon Shadow by HAWC

An indirect measurement of the antiproton flux in cosmic rays is possible as the particles undergo deflection by the geomagnetic field. This effect can be measured by studying the deficit in the flux, or shadow, created by the Moon as it absorbs cosmic rays that are headed towards the Earth. The shadow is displaced from the actual position of the Moon due to geomagnetic deflection, which is a function of the energy and charge of the cosmic rays. The displacement provides a natural tool for momentum/charge discrimination that can be used to study the composition of cosmic rays. Using 33 months of data comprising more than 80 billion cosmic rays measured by the High Altitude Water Cherenkov (HAWC) observatory, we have analyzed the Moon shadow to search for TeV antiprotons in cosmic rays. We present our first upper limits on the $\bar{p}/p$ fraction, which in the absence of any direct measurements, provide the tightest available constraints of $\sim1\%$ on the antiproton fraction for energies between 1 and 10 TeV.Comment: 10 pages, 5 figures. Accepted by Physical Review

Search for Point Sources of Ultra-High Energy Cosmic Rays Above 40 EeV Using a Maximum Likelihood Ratio Test

We present the results of a search for cosmic ray point sources at energies above 40 EeV in the combined data sets recorded by the AGASA and HiRes stereo experiments. The analysis is based on a maximum likelihood ratio test using the probability density function for each event rather than requiring an a priori choice of a fixed angular bin size. No statistically significant clustering of events consistent with a point source is found.Comment: 7 pages, 7 figures. Accepted for publication in The Astrophysical Journa

Very high energy particle acceleration powered by the jets of the microquasar SS 433

SS 433 is a binary system containing a supergiant star that is overflowing its Roche lobe with matter accreting onto a compact object (either a black hole or neutron star). Two jets of ionized matter with a bulk velocity of $\sim0.26c$ extend from the binary, perpendicular to the line of sight, and terminate inside W50, a supernova remnant that is being distorted by the jets. SS 433 differs from other microquasars in that the accretion is believed to be super-Eddington, and the luminosity of the system is $\sim10^{40}$ erg s$^{-1}$. The lobes of W50 in which the jets terminate, about 40 pc from the central source, are expected to accelerate charged particles, and indeed radio and X-ray emission consistent with electron synchrotron emission in a magnetic field have been observed. At higher energies (>100 GeV), the particle fluxes of $\gamma$ rays from X-ray hotspots around SS 433 have been reported as flux upper limits. In this energy regime, it has been unclear whether the emission is dominated by electrons that are interacting with photons from the cosmic microwave background through inverse-Compton scattering or by protons interacting with the ambient gas. Here we report TeV $\gamma$-ray observations of the SS 433/W50 system where the lobes are spatially resolved. The TeV emission is localized to structures in the lobes, far from the center of the system where the jets are formed. We have measured photon energies of at least 25 TeV, and these are certainly not Doppler boosted, because of the viewing geometry. We conclude that the emission from radio to TeV energies is consistent with a single population of electrons with energies extending to at least hundreds of TeV in a magnetic field of $\sim16$~micro-Gauss.Comment: Preprint version of Nature paper. Contacts: S. BenZvi, B. Dingus, K. Fang, C.D. Rho , H. Zhang, H. Zho

Measurement of the Crab Nebula Spectrum Past 100 TeV with HAWC

We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy-estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much improved energy resolution over the previous analysis and extends HAWC's ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log parabola shape $\left(\frac{dN}{dE} = \phi_0 \left(E/\textrm{7 TeV}\right)^{-\alpha-\beta\ln\left(E/\textrm{7 TeV}\right)}\right)$ with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 meters from the shower axis, the best-fit values are $\phi_o$=(2.35$\pm$0.04$^{+0.20}_{-0.21}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.79$\pm$0.02$^{+0.01}_{-0.03}$, and $\beta$=0.10$\pm$0.01$^{+0.01}_{-0.03}$. For the second estimator, a neural network which uses the charge distribution in annuli around the core and other variables, these values are $\phi_o$=(2.31$\pm$0.02$^{+0.32}_{-0.17}$)$\times$10$^{-13}$ (TeV cm$^2$ s)$^{-1}$, $\alpha$=2.73$\pm$0.02$^{+0.03}_{-0.02}$, and $\beta$=0.06$\pm$0.01$\pm$0.02. The first set of uncertainties are statistical; the second set are systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date.Comment: published in Ap