The STATFLUX code: a statistical method for calculation of flow and set of parameters, based on the Multiple-Compartment Biokinetical Model

The code STATFLUX, implementing a new and simple statistical procedure for the calculation of transfer coefficients in radionuclide transport to animals and plants, is proposed. The method is based on the general multiple-compartment model, which uses a system of linear equations involving geometrical volume considerations. Flow parameters were estimated by employing two different least-squares procedures: Derivative and Gauss-Marquardt methods, with the available experimental data of radionuclide concentrations as the input functions of time. The solution of the inverse problem, which relates a given set of flow parameter with the time evolution of concentration functions, is achieved via a Monte Carlo Simulation procedure.Program summaryTitle of program: STATFLUXCatalogue identifier: ADYS_v1_0Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADYS_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions: noneComputer for which the program is designed and others on which it has been tested: Micro-computer with Intel Pentium III, 3.0 GHzInstallation: Laboratory of Linear Accelerator, Department of Experimental Physics, University of São Paulo, BrazilOperating system: Windows 2000 and Windows XPProgramming language used: Fortran-77 as implemented in Microsoft Fortran 4.0. NOTE: Microsoft Fortran includes non-standard features which are used in this program. Standard Fortran compilers such as, g77, f77, ifort and NAG95, are not able to compile the code and therefore it has not been possible for the CPC Program Library to test the program.Memory, required to execute with typical data: 8 Mbytes of RAM memory and 100 MB of Hard disk memoryNo. of bits in a word: 16No. of lines in distributed program, including test data, etc.: 6912No. of bytes in distributed Program, including test data, etc.: 229 541Distribution format: tar.gzNature of the physical problem: the investigation of transport mechanisms for radioactive substances, through environmental pathways, is very important for radiological protection of populations. One such pathway, associated with the food chain, is the grass-animal-man sequence. The distribution of trace elements in humans and laboratory animals has been intensively studied over the past 60 years [R.C. Pendlenton, C.W. Mays, R.D. Lloyd, A.L. Brooks, Differential accumulation of iodine-131 from local fallout in people and milk, Health Phys. 9 (1963) 1253-1262]. In addition, investigations on the incidence of cancer in humans, and a possible causal relationship to radioactive fallout, have been undertaken [E.S. Weiss, M.L. Rallison, W.T. London, W.T. Carlyle Thompson, Thyroid nodularity in southwestern Utah school children exposed to fallout radiation, Amer. J. Public Health 61 (1971) 241-249; M.L. Rallison, B.M. Dobyns, F.R. Keating, J.E. Rall, F.H. Tyler, Thyroid diseases in children, Amer. J. Med. 56 (1974) 457-463; J.L. Lyon, M.R. Klauber, J.W. Gardner, K.S. Udall, Childhood leukemia associated with fallout from nuclear testing, N. Engl. J. Med. 300 (1979) 397-402]. From the pathways of entry of radionuclides in the human (or animal) body, ingestion is the most important because it is closely related to life-long alimentary (or dietary) habits. Those radionuclides which are able to enter the living cells by either metabolic or other processes give rise to localized doses which can be very high. The evaluation of these internally localized doses is of paramount importance for the assessment of radiobiological risks and radiological protection. The time behavior of trace concentration in organs is the principal input for prediction of internal doses after acute or chronic exposure. The General Multiple-Compartment Model (GMCM) is the powerful and more accepted method for biokinetical studies, which allows the calculation of concentration of trace elements in organs as a function of time, when the flow parameters of the model are known. However, few biokinetics data exist in the literature, and the determination of flow and transfer parameters by statistical fitting for each system is an open problem.Restriction on the complexity of the problem: This version of the code works with the constant volume approximation, which is valid for many situations where the biological half-live of a trace is lower than the volume rise time. Another restriction is related to the central flux model. The model considered in the code assumes that exist one central compartment (e.g., blood), that connect the flow with all compartments, and the flow between other compartments is not included.Typical running time: Depends on the choice for calculations. Using the Derivative Method the time is very short (a few minutes) for any number of compartments considered. When the Gauss-Marquardt iterative method is used the calculation time can be approximately 5-6 hours when similar to 15 compartments are considered. (C) 2006 Elsevier B.V. All rights reserved.

Modelling frequency distributions of 5 minute-averaged solar radiation indexes using Beta probability functions

Five minute-averaged values of sky clearness, direct and diffuse indices, were used to model the frequency distributions of these variables in terms of optical air mass. From more than four years of solar radiation observations it was found that variations in the frequency distributions of the three indices of optical air mass for Botucatu, Brazil, are similar to those in other places, as published in the literature. The proposed models were obtained by linear combination of normalized Beta probability functions, using the observed distributions derived from three years of data. The versatility of these functions allows modelling of all three irradiance indexes to similar levels of accuracy. A comparison with the observed distributions obtained from one year of observations indicate that the models are able to reproduce the observed frequency distributions of all three indices at the 95% confidence level.

Friction losses in valves and fittings for viscoplastic fluids

Data on pressure drop were obtained in stainless steel, sanitary fittings and valves during laminar and turbulent flow of aqueous suspensions of sucrose and bentonite. The rheological properties of these suspensions were determined and the Bingham model provided the best fitting with the experimental data. Friction losses were measured in fully- and partially-open butterfly and plug valves, bends and union. Values of loss coefficients (k(f)) were calculated and correlated as functions of the classical Reynolds number and the Reynolds number proposed by Govier and Aziz (1972) for viscoplastic fluids. The two-k method and a new proposed model presented the best adjustments for the Govier and Aziz Reynolds number, and Hedstrom and classical Reynolds numbers, respectively.

Equivalence of many-gluon Green's functions in the Duffin-Kemmer-Petieu and Klein-Gordon-Fock statistical quantum field theories

We prove the equivalence of many-gluon Green's functions in the Duffin-Kemmer-Petieu and Klein-Gordon-Fock statistical quantum field theories. The proof is based on the functional integral formulation for the statistical generating functional in a finite-temperature quantum field theory. As an illustration, we calculate one-loop polarization operators in both theories and show that their expressions indeed coincide.

MAPPING CORRELATION-FUNCTIONS OF SCHWINGER AND AXIAL MODELS

We show that all Green's functions of the Schwinger and axial models can be obtained one from the other. In particular, we show that the two models have the same chiral anomaly. Finally it is demonstrated that the Schwinger model can keep gauge invariance for an arbitrary mass, dispensing with an additional gauge group integration.

Fitting hyperbolic universes to Cayon-Smoot spots in COBE maps

On the possibility that the universe's matter density is low (Ohm(0) < 1), cosmologies can be considered with the metric of Friedmann's open universe but with closed hyperbolic manifolds as the physical three-space. These models have nontrivial spatial topology, with the property of producing multiple images of cosmic sources. Here a fit is attempted of 10 of these models to the physical cold and hot spots found by Cayon & Smoot in the COBE/DMR maps. These spots are interpreted as early, distant images of much nearer sources of inhomogeneity. The source for one of the cold spots is seen as the seed of a known supercluster.

Conjectures and theorems in the theory of entire functions

Motivated by the recent solution of Karlin's conjecture, properties of functions in the Laguerre-Polya class are investigated. The main result of this paper establishes new moment inequalities fur a class of entire functions represented by Fourier transforms. The paper concludes with several conjectures and open problems involving the Laguerre-Polya class and the Riemann xi -function.

Perturbative bosonization from two-point correlation functions

Here we address the problem of bosonizing massive fermions without making expansions in the fermion masses in both massive QED(2) and QED(3) with N fermion flavors including also a Thirring coupling. We start from two-point correlators involving the U(1) fermionic current and the gauge field. From the tensor structure of those correlators we prove that the U(1) current must be identically conserved (topological) in the corresponding bosonized theory in both D=2 and D=3 dimensions. We find an effective generating functional in terms of bosonic fields which reproduces these two-point correlators and from that we obtain a map of the Lagrangian density (ψ) over bar (r)(ipartial derivative-m)psi(r) into a bosonic one in both dimensions. This map is nonlocal but it is independent of the electromagnetic and Thirring couplings, at least in the quadratic approximation for the fermionic determinant.

Kinetic aspects of ion exchange in KhFek[Fe(CN)(6)](l)center dot mH(2)O compounds: A combined electrical and mass transfer functions approach

The present paper quantifies and develops the kinetic aspects involved in the mechanism of interplay between electron and ions presented elsewhere(1) for KhFek[Fe(CN)(6)](l)center dot mH(2)O (Prussian Blue) host materials. Accordingly, there are three different electrochemical processes involved in the PB host materials: H3O+, K+, and H+ insertion/extraction mechanisms which here were fully kinetically studied by means of the use of combined electronic and mass transfer functions as a tool to separate all the processes. The use of combined electronic and mass transfer functions was very important to validate and confirm the proposed mechanism. This mechanism allows the electrochemical and chemical processes involved in the KhFek[Fe(CN)(6)](l)center dot mH(2)O host and Prussian Blue derivatives to be understood. In addition, a formalism was also developed to consider superficial oxygen reduction. From the analysis of the kinetic processes involved in the model, it was possible to demonstrate that the processes associated with K+ and H+ exchanges are reversible whereas the H3O+ insertion process was shown not to present a reversible pattern. This irreversible pattern is very peculiar and was shown to be related to the catalytic proton reduction reaction. Furthermore, from the model, it was possible to calculate the number density of available sites for each intercalation/deintercalation processes and infer that they are very similar for K+ and H+. Hence, the high prominence of the K+ exchange observed in the voltammetric responses has a kinetic origin and is not related to the amount of sites available for intercalation/deintercalation of the ions.

HELICITY WAVE-FUNCTIONS FOR MASSLESS AND MASSIVE SPIN-2 PARTICLES

Starting from general properties of a spin-2 field, we construct helicity wave functions in the framework of the Weyl-van der Waerden spinor formalism. We discuss here the cases of massless and massive spin-2 particles.

Effects of guanethidine-induced sympathectomy on the spermatogenic and steroidogenic testicular functions of prepubertal to mature rats

Selective chemical sympathectomy of the internal genital organs of prepubertal to mature male Wistar rats was performed by chronic treatment with low doses of guanethidine. Sympathetic denervation caused an increase in intratesticular progesterone levels in prepubertal and early pubertal rats in addition to a decrease in androstenedione and testosterone levels in prepubertal animals, thus indicating a decrease in the conversion of progesterone into androgen, probably by blocking the steroidogenic enzymatic pathway at the 17 alpha-hydroxylase/17,20 desmolase level. A lower degree of testicular maturation, probably related to reduced androgen activity, was observed in prepubertal and early pubertal sympathectomized rats. Concentration of spermatozoa, on the other hand, was increased in the enlarged cauda epididymidis of late pubertal and mature denervated animals. This result is discussed in terms of the impairment of epididymal mechanisms of seminal emission, fluid resorption and spermatozoal disposal.

Integration of polyharmonic functions

The results in this paper are motivated by two analogies. First, m-harmonic functions in R(n) are extensions of the univariate algebraic polynomials of odd degree 2m-1. Second, Gauss' and Pizzetti's mean value formulae are natural multivariate analogues of the rectangular and Taylor's quadrature formulae, respectively. This point of view suggests that some theorems concerning quadrature rules could be generalized to results about integration of polyharmonic functions. This is done for the Tchakaloff-Obrechkoff quadrature formula and for the Gaussian quadrature with two nodes.

Temperature and concentration dependence of density of model liquid foods

Density of binary solutions and combinations of sucrose, glucose, fructose, citric acid, malic acid, pectin, and inorganic salts were measured with an oscillating tube density meter in the temperature range from 10degrees to 60degreesC, at varying concentrations. Density can be predicted with accuracy better than 5 x 10(-5) g cm(-3) using predictive equations obtained by fitting the experimental data. Available literature values agreed well with experimental data. Relations for the excess molar volume of these solutions were derived in terms of mole fraction and temperature. A thermodynamic model for the volumetric analysis of multicomponent aqueous solutions containing electrolyte and non-electrolyte compounds was also proposed. These models can be used for prediction of density of liquid food systems, specially fruit juices and beverages, based on composition and temperature, with high accuracy and without elaborate experimental work.