Simplified Mathematical Model for an Anaerobic Sequencing Batch Biofilm Reactor Treating Lipid-Rich Wastewater Subject to Rising Organic Loading Rates

This study proposes a simplified mathematical model to describe the processes occurring in an anaerobic sequencing batch biofilm reactor (ASBBR) treating lipid-rich wastewater. The reactor, subjected to rising organic loading rates, contained biomass immobilized cubic polyurethane foam matrices, and was operated at 32 degrees C +/- 2 degrees C, using 24-h batch cycles. In the adaptation period, the reactor was fed with synthetic substrate for 46 days and was operated without agitation. Whereas agitation was raised to 500 rpm, the organic loading rate (OLR) rose from 0.3 g chemical oxygen demand (COD) . L(-1) . day(-1) to 1.2 g COD . L(-1) . day(-1). The ASBBR was fed fat-rich wastewater (dairy wastewater), in an operation period lasting for 116 days, during which four operational conditions (OCs) were tested: 1.1 +/- 0.2 g COD . L(-1) . day(-1) (OC1), 4.5 +/- 0.4 g COD . L(-1) . day(-1) (OC2), 8.0 +/- 0.8 g COD . L(-1) . day(-1) (OC3), and 12.1 +/- 2.4 g COD . L(-1) . day(-1) (OC4). The bicarbonate alkalinity (BA)/COD supplementation ratio was 1:1 at OC1, 1:2 at OC2, and 1:3 at OC3 and OC4. Total COD removal efficiencies were higher than 90%, with a constant production of bicarbonate alkalinity, in all OCs tested. After the process reached stability, temporal profiles of substrate consumption were obtained. Based on these experimental data a simplified first-order model was fit, making possible the inference of kinetic parameters. A simplified mathematical model correlating soluble COD with volatile fatty acids (VFA) was also proposed, and through it the consumption rates of intermediate products as propionic and acetic acid were inferred. Results showed that the microbial consortium worked properly and high efficiencies were obtained, even with high initial substrate concentrations, which led to the accumulation of intermediate metabolites and caused low specific consumption rates.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Turbulent magnetic pumping in a Babcock-Leighton solar dynamo model

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. Aims. We explore the effects of turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. Methods. As a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction equation. In the second step, a latitudinal pumping term was included and then, a near-surface shear was included. Results. The results reveal the importance of the pumping mechanism in solving current limitations in mean field dynamo modeling, such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. If a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspot activity. In models that consider shear near the surface, a second shell of toroidal field is generated above r = 0.95 R(circle dot) at all latitudes. If the full pumping is also included, the polar toroidal fields are efficiently advected inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint.

uvby-beta photometry of solar twins The solar colors, model atmospheres, and the T(eff) and metallicity scales

Aims. Solar colors have been determined on the uvby-beta photometric system to test absolute solar fluxes, to examine colors predicted by model atmospheres as a function of stellar parameters (T(eff), log g, [Fe/H]), and to probe zero-points of T(eff) and metallicity scales. Methods. New uvby-beta photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our data with photometry available in the literature, the solar colors on the uvby-beta system have been inferred. Our solar colors are compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model atmospheres. The zero-points of different T(eff) and metallicity scales are verified and corrections are proposed. Results. Our solar colors are (b - y)(circle dot) = 0.4105 +/- 0.0015, m(1,circle dot) = 0.2122 +/- 0.0018, c(1,circle dot) = 0.3319 +/- 0.0054, and beta(circle dot) = 2.5915 +/- 0.0024. The (b - y)(circle dot) and m(1,circle dot) colors obtained from absolute spectrophotometry of the Sun agree within 3-sigma with the solar colors derived here when the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the c(1,circle dot) and beta(circle dot) synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values computed from the Kurucz models, (b - y)(circle dot) and m(1,circle dot) are in excellent agreement with our solar colors independently of the adopted zero-points, but for c(1,circle dot) and beta circle dot agreement is found only when adopting the ATLAS9 zero-points. The c(1,circle dot) color computed from both the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u band. The T(eff) calibration of Alonso and collaborators has the poorest performance (similar to 140 K off), while the relation of Casagrande and collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (similar to 10%) zero-point correction to place the stars and the Sun on the same metallicity scale. Finally, we confirm that the c(1) index in solar analogs has a strong metallicity sensitivity.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

Human Dental Pulp Cells: A New Source of Cell Therapy in a Mouse Model of Compressive Spinal Cord Injury

Strategies aimed at improving spinal cord regeneration after trauma are still challenging neurologists and neuroscientists throughout the world. Many cell-based therapies have been tested, with limited success in terms of functional outcome. In this study, we investigated the effects of human dental pulp cells (HDPCs) in a mouse model of compressive spinal cord injury (SCI). These cells present some advantages, such as the ease of the extraction process, and expression of trophic factors and embryonic markers from both ecto-mesenchymal and mesenchymal components. Young adult female C57/BL6 mice were subjected to laminectomy at T9 and compression of the spinal cord with a vascular clip for 1 min. The cells were transplanted 7 days or 28 days after the lesion, in order to compare the recovery when treatment is applied in a subacute or chronic phase. We performed quantitative analyses of white-matter preservation, trophic-factor expression and quantification, and ultrastructural and functional analysis. Our results for the HDPC-transplanted animals showed better white-matter preservation than the DMEM groups, higher levels of trophic-factor expression in the tissue, better tissue organization, and the presence of many axons being myelinated by either Schwann cells or oligodendrocytes, in addition to the presence of some healthy-appearing intact neurons with synapse contacts on their cell bodies. We also demonstrated that HDPCs were able to express some glial markers such as GFAP and S-100. The functional analysis also showed locomotor improvement in these animals. Based on these findings, we propose that HDPCs may be feasible candidates for therapeutic intervention after SCI and central nervous system disorders in humans.

Renal denervation in an animal model of diabetes and hypertension: Impact on the autonomic nervous system and nephropathy

Background: The effects of renal denervation on cardiovascular reflexes and markers of nephropathy in diabetic-hypertensive rats have not yet been explored. Methods: Aim: To evaluate the effects of renal denervation on nephropathy development mechanisms (blood pressure, cardiovascular autonomic changes, renal GLUT2) in diabetic-hypertensive rats. Forty-one male spontaneously hypertensive rats (SHR) similar to 250 g were injected with STZ or not; 30 days later, surgical renal denervation (RD) or sham procedure was performed; 15 days later, glycemia and albuminuria (ELISA) were evaluated. Catheters were implanted into the femoral artery to evaluate arterial pressure (AP) and heart rate variability (spectral analysis) one day later in conscious animals. Animals were killed, kidneys removed, and cortical renal GLUT2 quantified (Western blotting). Results: Higher glycemia (p < 0.05) and lower mean AP were observed in diabetics vs. nondiabetics (p < 0.05). Heart rate was higher in renal-denervated hypertensive and lower in diabetic-hypertensive rats (384.8 +/- 37, 431.3 +/- 36, 316.2 +/- 5, 363.8 +/- 12 bpm in SHR, RD-SHR, STZ-SHR and RD-STZ-SHR, respectively). Heart rate variability was higher in renal-denervated diabetic-hypertensive rats (55.75 +/- 25.21, 73.40 +/- 53.30, 148.4 +/- 93 in RD-SHR, STZ-SHR-and RD-STZ-SHR, respectively, p < 0.05), as well as the LF component of AP variability (1.62 +/- 0.9, 2.12 +/- 0.9, 7.38 +/- 6.5 in RD-SHR, STZ-SHR and RD-STZ-SHR, respectively, p < 0.05). GLUT2 renal content was higher in all groups vs. SHR. Conclusions: Renal denervation in diabetic-hypertensive rats improved previously reduced heart rate variability. The GLUT2 equally overexpressed by diabetes and renal denervation may represent a maximal derangement effect of each condition.

Melatonin treatment decreases c-fos expression in a headache model induced by capsaicin

The aim of the present work was to analyze c-fos response within the trigeminal nucleus caudalis (TNC) of pinealectomized rats and animals that received intraperitoneal melatonin, after intracisternal infusion of capsaicin, used to induce intracranial trigeminovascular stimulation. Experimental groups consisted of animals that received vehicle solution (saline-ethanol-Tween 80, 8:1:1, diluted 1:50) only (VEI, n = 5); animals that received capsaicin solution (200 nM) only (CAP, n = 6); animals submitted to pinealectomy (PX, n = 5); sham-operated animals (SH, n = 5); animals submitted to pinealectomy followed by capsaicin stimulation (200 nM) after 15 days (PX + CAP, n = 7); and animals that received capsaicin solution (200 nM) and intraperitoneal melatonin (10 mg/kg) (CAP + MEL, n = 5). Control rats, receiving vehicle in the cisterna magna, showed a small number of c-fos-positive cells in the TNC (layer I/II) as well as the sham-operated and pinealectomized rats, when compared to animals stimulated by capsaicin. On the other hand, pinealectomized rats, which received capsaicin, presented the highest number of c-fos-positive cells. Animals receiving capsaicin and melatonin treatment had similar expression of the vehicle group. Our data provide experimental evidence to support the role of melatonin and pineal gland in the pathophysiology of neurovascular headaches.

Respiratory allergy to Blomia tropicalis: Immune response in four syngeneic mouse strains and assessment of a low allergen-dose, short-term experimental model

Background: The dust mite Blomia tropicalis is an important source of aeroallergens in tropical areas. Although a mouse model for B. tropicalis extract (BtE)-induced asthma has been described, no study comparing different mouse strains in this asthma model has been reported. The relevance and reproducibility of experimental animal models of allergy depends on the genetic background of the animal, the molecular composition of the allergen and the experimental protocol. Objectives: This work had two objectives. The first was to study the anti-B. tropicalis allergic responses in different mouse strains using a short-term model of respiratory allergy to BtE. This study included the comparison of the allergic responses elicited by BtE with those elicited by ovalbumin in mice of the strain that responded better to BtE sensitization. The second objective was to investigate whether the best responder mouse strain could be used in an experimental model of allergy employing relatively low BtE doses. Methods: Groups of mice of four different syngeneic strains were sensitized subcutaneously with 100 mu g of BtE on days 0 and 7 and challenged four times intranasally, at days 8, 10, 12, and 14, with 10 mu g of BtE. A/J mice, that were the best responders to BtE sensitization, were used to compare the B. tropicalis-specific asthma experimental model with the conventional experimental model of ovalbumin (OVA)-specific asthma. A/J mice were also sensitized with a lower dose of BtE. Results: Mice of all strains had lung inflammatory-cell infiltration and increased levels of anti-BtE IgE antibodies, but these responses were significantly more intense in A/J mice than in CBA/J, BALB/c or C57BL/6J mice. Immunization of A/J mice with BtE induced a more intense airway eosinophil influx, higher levels of total IgE, similar airway hyperreactivity to methacholine but less intense mucous production, and lower levels of specific IgE, IgG1 and IgG2 antibodies than sensitization with OVA. Finally, immunization with a relatively low BtE dose (10 mu g per subcutaneous injection per mouse) was able to sensitize A/J mice, which were the best responders to high-dose BtE immunization, for the development of allergy-associated immune and lung inflammatory responses. Conclusions: The described short-term model of BtE-induced allergic lung disease is reproducible in different syngeneic mouse strains, and mice of the A/J strain was the most responsive to it. In addition, it was shown that OVA and BtE induce quantitatively different immune responses in A/J mice and that the experimental model can be set up with low amounts of BtE.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 2: Model sensitivity to the biomass burning

We describe an estimation technique for biomass burning emissions in South America based on a combination of remote-sensing fire products and field observations, the Brazilian Biomass Burning Emission Model (3BEM). For each fire pixel detected by remote sensing, the mass of the emitted tracer is calculated based on field observations of fire properties related to the type of vegetation burning. The burnt area is estimated from the instantaneous fire size retrieved by remote sensing, when available, or from statistical properties of the burn scars. The sources are then spatially and temporally distributed and assimilated daily by the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) in order to perform the prognosis of related tracer concentrations. Three other biomass burning inventories, including GFEDv2 and EDGAR, are simultaneously used to compare the emission strength in terms of the resultant tracer distribution. We also assess the effect of using the daily time resolution of fire emissions by including runs with monthly-averaged emissions. We evaluate the performance of the model using the different emission estimation techniques by comparing the model results with direct measurements of carbon monoxide both near-surface and airborne, as well as remote sensing derived products. The model results obtained using the 3BEM methodology of estimation introduced in this paper show relatively good agreement with the direct measurements and MOPITT data product, suggesting the reliability of the model at local to regional scales.

Model Independent Tests of Skyrmions and Their Holographic Cousins

We describe a new exact relation for large N(c) QCD for the long-distance behavior of baryon form factors in the chiral limit. This model-independent relation is used to test the consistency of the structure of several baryon models. All 4D semiclassical chiral soliton models satisfy the relation, as does the Pomarol-Wulzer holographic model of baryons as 5D Skyrmions. However, remarkably, we find that the holographic model treating baryons as instantons in the Sakai-Sugimoto model does not satisfy the relation.

Charm and longitudinal structure functions within the Kharzeev-Levin-Nardi model

We use the Kharzeev-Levin-Nardi (KLN) model of the low x gluon distributions to fit recent HERA data on F(L) and F(2)(c)(F(2)(b)). Having checked that this model gives a good description of the data, we use it to predict F(L) and F(2)(c) to be measured in a future electron-ion collider. The results are similar to those obtained with the de Florian-Sassot and Eskola-Paukkunen-Salgado nuclear gluon distributions. The conclusion of this exercise is that the KLN model, simple as it is, may still be used as an auxiliary tool to make estimates for both heavy-ion and electron-ion collisions.

Phase diagram of a model for a binary mixture of nematic molecules on a Bethe lattice

We investigate the phase diagram of a discrete version of the Maier-Saupe model with the inclusion of additional degrees of freedom to mimic a distribution of rodlike and disklike molecules. Solutions of this problem on a Bethe lattice come from the analysis of the fixed points of a set of nonlinear recursion relations. Besides the fixed points associated with isotropic and uniaxial nematic structures, there is also a fixed point associated with a biaxial nematic structure. Due to the existence of large overlaps of the stability regions, we resorted to a scheme to calculate the free energy of these structures deep in the interior of a large Cayley tree. Both thermodynamic and dynamic-stability analyses rule out the presence of a biaxial phase, in qualitative agreement with previous mean-field results.

Critical behavior of the susceptible-infected-recovered model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c(0)=0.176 500 5(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda(c)=(1-c(0))/c(0)=4.665 71(3) and a net transmissibility of (1-c(0))/(1+3c(0))=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.