Combining Legendre's Polynomials and Genetic Algorithm in the Solution of Nonlinear Initial-Value Problems

This work develops a method for solving ordinary differential equations, that is, initial-value problems, with solutions approximated by using Legendre's polynomials. An iterative procedure for the adjustment of the polynomial coefficients is developed, based on the genetic algorithm. This procedure is applied to several examples providing comparisons between its results and the best polynomial fitting when numerical solutions by the traditional Runge-Kutta or Adams methods are available. The resulting algorithm provides reliable solutions even if the numerical solutions are not available, that is, when the mass matrix is singular or the equation produces unstable running processes.

Rumor Propagation Model: An Equilibrium Study

Compartmental epidemiological models have been developed since the 1920s and successfully applied to study the propagation of infectious diseases. Besides, due to their structure, in the 1960s an interesting version of these models was developed to clarify some aspects of rumor propagation, considering that spreading an infectious disease or disseminating information is analogous phenomena. Here, in an analogy with the SIR (Susceptible-Infected-Removed) epidemiological model, the ISS (Ignorant-Spreader-Stifler) rumor spreading model is studied. By using concepts from the Dynamical Systems Theory, stability of equilibrium points is established, according to propagation parameters and initial conditions. Some numerical experiments are conducted in order to validate the model.

The non-equilibrium nature of culinary evolution

Food is an essential part of civilization, with a scope that ranges from the biological to the economic and cultural levels. Here, we study the statistics of ingredients and recipes taken from Brazilian, British, French and Medieval cookery books. We find universal distributions with scale invariant behaviour. We propose a copy-mutate process to model culinary evolution that fits our empirical data very well. We find a cultural 'founder effect' produced by the non-equilibrium dynamics of the model. Both the invariant and idiosyncratic aspects of culture are accounted for by our model, which may have applications in other kinds of evolutionary processes.

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.

Analytical solutions for Tokamak equilibria with reversed toroidal current

In tokamaks, an advanced plasma confinement regime has been investigated with a central hollow electric current with negative density which gives rise to non-nested magnetic surfaces. We present analytical solutions for the magnetohydrodynamic equilibria of this regime in terms of non-orthogonal toroidal polar coordinates. These solutions are obtained for large aspect ratio tokamaks and they are valid for any kind of reversed hollow current density profiles. The zero order solution of the poloidal magnetic flux function describes nested toroidal magnetic surfaces with a magnetic axis displaced due to the toroidal geometry. The first order correction introduces a poloidal field asymmetry and, consequently, magnetic islands arise around the zero order surface with null poloidal magnetic flux gradient. An analytic expression for the magnetic island width is deduced in terms of the equilibrium parameters. We give examples of the equilibrium plasma profiles and islands obtained for a class of current density profile. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624551]

Magnetoconfined levels in a parabolic quantum dot: An analytical solution of a three-dimensional Fock-Darwin problem in a tilted magnetic field

The energy spectrum of an electron confined in a quantum dot (QD) with a three-dimensional anisotropic parabolic potential in a tilted magnetic field was found analytically. The theory describes exactly the mixing of in-plane and out-of-plane motions of an electron caused by a tilted magnetic field, which could be seen, for example, in the level anticrossing. For charged QDs in a tilted magnetic field we predict three strong resonant lines in the far-infrared-absorption spectra.

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

Determination of the molecular weight of proteins in solution from a single small-angle X-ray scattering measurement on a relative scale

This paper describes a new and simple method to determine the molecular weight of proteins in dilute solution, with an error smaller than similar to 10%, by using the experimental data of a single small-angle X-ray scattering (SAXS) curve measured on a relative scale. This procedure does not require the measurement of SAXS intensity on an absolute scale and does not involve a comparison with another SAXS curve determined from a known standard protein. The proposed procedure can be applied to monodisperse systems of proteins in dilute solution, either in monomeric or multimeric state, and it has been successfully tested on SAXS data experimentally determined for proteins with known molecular weights. It is shown here that the molecular weights determined by this procedure deviate from the known values by less than 10% in each case and the average error for the test set of 21 proteins was 5.3%. Importantly, this method allows for an unambiguous determination of the multimeric state of proteins with known molecular weights.

Characterization and greenhouse evaluation of Brazilian calcined nonapatite phosphate rocks for rice

Little information is available on the agronomic effectiveness of calcined nonapatite phosphate rock (PR) sources containing crandallite minerals in the form of Ca-Fe-Al-P for flooded and upland rice (Oryza sativa L.). We conducted laboratory and greenhouse studies to (i) characterize the mineralogical composition, (ii) investigate the solubility and dissolution behavior, and (iii) evaluate the agronomic effectiveness of two nonapatite PR sources (Juquia and Sapucaia) from Brazil and compared them with (i) a highly reactive Gafsa PR (Tunisia) containing apatite in the form of Ca-P and (ii) a reference water-soluble triple superphosphate (TSP) for flooded and upland rice. After calcination at 500 degrees C for 4 h, the solubility of Juquia PR and Sapucaia PR in neutral ammonium citrate (NAC) significantly increased from almost nil to a maximum of 39.3 and 114 g P kg(-1), respectively. X-ray diffraction showed that crystalline crandallite mineral was transformed to an amophorus form after calcination. The solubility behavior of the two calcined PR sources followed the same trend as Gafsa PR, that is, P release decreased with increasing equilibrium pH in the 0.01 M KCl solution (PH 3.0-8.0). At PH 3, the solubility followed: Gafsa PR > calcined Sapucaia PR > calcined Juquia PR. No P release was detected from any of the PR sources at pH >= 5.0 in the solution, indicating the Ca-P characteristic of the Ca-Fe-Al-P mineral controlled P dissolution of the calcined PR. Without calcination, both Juquia PR and Sapucaia PR were totally ineffective for upland rice grown on a Hiwassee clay loam (fine, kaolinitic, thermic Rhodic Kanhapludult) with pH 5.4 whereas a significant P response was observed with the calcined PR samples. For flooded rice grown on Hiwassee soil, the calcined Juquia PR and Sapucaia PR were 66 and 72%, respectively, as effective as TSP in increasing rice grain yield whereas Gafsa PR was ineffective. For upland rice grown on the unlimed soil, Gafsa PR was as effective as TSP in increasing rice grain yield whereas calcined Juquia PR and Sapucaia PR were 89 and 83% of TSP. The effectiveness of Gafsa PR was reduced to 0% after the soil was limed to pH 7.0 whereas the two calcined PR sources were reduced to 49% of TSP. Soil available P extracted by iron oxide impregnated filter paper (Pi test) or anion-exchange resin after rice harvest correlated well with P uptake by rice grain for flooded and upland rice.

Kinetics of the xylitol crystallization in hydro-alcoholic solution

Nyvlt method Was used to determine the kinetic parameters of commercial xylitol in ethanol:water (50:50 %w/w) Solution by batch cooling crystallization. The kinetic exponents (n, g and in) and the system kinetic constant (B(N)) were determined. Model experiments were carried Out in order to verify the combined effects of saturation temperatures (40, 50 and 60 degrees C) and cooling rates (0.10, 0.25 and 0.50 degrees C/min) on these parameters. The fitting between experimental and Calculated crystal sizes has 11.30% mean deviation. (C) 2007 Elsevier B.V. All rights reserved.

Thermodynamic and kinetic investigations of phosphate adsorption onto hydrous niobium oxide prepared by homogeneous solution method

The adsorption kinetics of phosphate onto Nb(2)O(5)center dot nH(2)O was investigated at initial phosphate concentrations 10 and 50 mg L(-1). The kinetic process was described by a pseudo second-order rate model very well. The adsorption thermodynamics was carried out at 298, 308, 318, 328 and 338 K. The positive values of both Delta H and Delta S suggest an endothermic reaction and increase in randomness at the solid-liquid interface during the adsorption. Delta G values obtained were negative indicating a spontaneous adsorption process. The Langmuir model described the data better than the Freundlich isotherm model. The peak appearing at 1050 cm(-1) in IR spectra after adsorption was attributed to the bending vibration of adsorbed phosphate. The effective desorption could be achieved using water at pH 12. (C) 2010 Elsevier B.V. All rights reserved.

Extending the kinetic solution of the classic Michaelis-Menten model of enzyme action

The principal aim of studies of enzyme-mediated reactions has been to provide comparative and quantitative information on enzyme-catalyzed reactions under distinct conditions. The classic Michaelis-Menten model (Biochem Zeit 49:333, 1913) for enzyme kinetic has been widely used to determine important parameters involved in enzyme catalysis, particularly the Michaelis-Menten constant (K (M) ) and the maximum velocity of reaction (V (max) ). Subsequently, a detailed treatment of the mechanisms of enzyme catalysis was undertaken by Briggs-Haldane (Biochem J 19:338, 1925). These authors proposed the steady-state treatment, since its applicability was constrained to this condition. The present work describes an extending solution of the Michaelis-Menten model without the need for such a steady-state restriction. We provide the first analysis of all of the individual reaction constants calculated analytically. Using this approach, it is possible to accurately predict the results under new experimental conditions and to characterize and optimize industrial processes in the fields of chemical and food engineering, pharmaceuticals and biotechnology.

Solution Polymerization of N-vinylcaprolactam in 1,4-dioxane. Kinetic Dependence on Temperature, Monomer, and Initiator Concentrations

The kinetics of the solution free radical polymerization of N-vinylcaprolactam, in 1,4-dioxane and under various polymerization conditions was studied. Azobisisobutyronitrile and 3-mercaptopropionic acid were used as initiator and as chain transfer agent (CTA), respectively. The influence of monomer and initiator concentrations and polymerization temperature on the rate of polymerizations (R(p)) was investigated. In general, high conversions were obtained. The order with respect to initiator was consistent with the classical kinetic rate equation, while the order with respect to the monomer was greater than unity. The overall activation energy of 53.6 kJ mol(-1) was obtained in the temperature range 60-80 degrees C. The decreasing of the absolute molecular weights when increasing the CIA concentration was confirmed by GPC/SEC/LALS analyses. It was confirmed by UV-visible analyses the effect of molecular weights on the lower critical solution temperature of the polymers. It was also verified that the addition of the CTA influenced the kinetic of the polymerizations. (C) 2010 Wiley Periodicals, Inc. J Appl Polym Sci 118: 229-240, 2010

Adsorption of Cr(VI) from aqueous solution by hydrous zirconium oxide

A type of ZrO(2)center dot nH(2)O Was synthesized and its Cr(VI) removal potential was investigated in this study. The kinetic study, adsorption isotherm, pH effect, thermodynamic study and desorption were examined in batch experiments. The kinetic process was described by a pseudo-second-order rate model very well. The Cr(VI) adsorption tended to increase with a decrease of pH. The adsorption data fitted well to the Langmuir model. The adsorption capacity increased from 61 to 66 mg g(-1) when the temperature was increased from 298 to 338 K. The positive values of both Delta H degrees and Delta S degrees suggest an endothermic reaction and increase in randomness at the solid-liquid interface during the adsorption. Delta G degrees values obtained were negative indicating a spontaneous adsorption process. The effective desorption of Cr(VI) on ZrO(2)center dot nH(2)O could be achieved using distilled water at pH 12. (C) 2009 Elsevier B.V. All rights reserved.