Hyperbolic reaction-diffusion equations for a forest fire model

Forest fire models have been widely studied from the context of self-organized criticality and from the ecological properties of the forest and combustion. On the other hand, reaction-diffusion equations have interesting applications in biology and physics. We propose here a model for fire propagation in a forest by using hyperbolic reaction-diffusion equations. The dynamical and thermodynamical aspects of the model are analyzed in detail

Reaction-diffusion wave fronts: Multigeneration biological species under climate change

A generalization of reaction-diffusion models to multigeneration biological species is presented. It is based on more complex random walks than those in previous approaches. The new model is developed analytically up to infinite order. Our predictions for the speed agree to experimental data for several butterfly species better than existing models. The predicted dependence for the speed on the number of generations per year allows us to explain the change in speed observed for a specific invasion

Speed of reaction-transport processes

We present an approach to determining the speed of wave-front solutions to reaction-transport processes. This method is more accurate than previous ones. This is explicitly shown for several cases of practical interest: (i) the anomalous diffusion reaction, (ii) reaction diffusion in an advective field, and (iii) time-delayed reaction diffusion. There is good agreement with the results of numerical simulations

Speed of wave-front solutions to hyperbolic reaction-diffusion equations

The asymptotic speed problem of front solutions to hyperbolic reaction-diffusion (HRD) equations is studied in detail. We perform linear and variational analyses to obtain bounds for the speed. In contrast to what has been done in previous work, here we derive upper bounds in addition to lower ones in such a way that we can obtain improved bounds. For some functions it is possible to determine the speed without any uncertainty. This is also achieved for some systems of HRD (i.e., time-delayed Lotka-Volterra) equations that take into account the interaction among different species. An analytical analysis is performed for several systems of biological interest, and we find good agreement with the results of numerical simulations as well as with available observations for a system discussed recently

Comparison of Work Scheduling Using Constraint Programming or Auctions

Business processes designers take into account the resources that the processes would need, but, due to the variable cost of certain parameters (like energy) or other circumstances, this scheduling must be done when business process enactment. In this report we formalize the energy aware resource cost, including time and usage dependent rates. We also present a constraint programming approach and an auction-based approach to solve the mentioned problem including a comparison of them and a comparison of the proposed algorithms for solving them

Second-order Møller-Plesset perturbation theory without basis set superposition error : II : open-shell systems

The basis set superposition error-free second-order MØller-Plesset perturbation theory of intermolecular interactions was studied. The difficulties of the counterpoise (CP) correction in open-shell systems were also discussed. The calculations were performed by a program which was used for testing the new variants of the theory. It was shown that the CP correction for the diabatic surfaces should be preferred to the adiabatic ones

The hardness profile as a tool to detect spurious stationary points in the potential energy surface

The energy and hardness profile for a series of inter and intramolecular conformational changes at several levels of calculation were computed. The hardness profiles were found to be calculated as the difference between the vertical ionization potential and electron affinity. The hardness profile shows the correct number of stationary points independently of the basis set and methodology used. It was found that the hardness profiles can be used to check the reliability of the energy profiles for those chemical system

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Predicted effects of small decreases in energy expenditure on weight gain in adult women.

Small daily positive energy imbalances of 200 to 800 kJ (about 50 to 200 kcal) due to reduced resting energy expenditure (REE), reduced diet-induced thermogenesis, or physical inactivity are believed to predispose to obesity. However, estimates of the magnitude of the weight gain often fail to account for concurrent changes in body composition and increases in maintenance energy requirements as weight increases and energy equilibrium is re-established. Using previously reported data on body composition and REE in women and the energy cost of tissue deposition, we used mathematical models to predict the theoretical effect of a persistent reduction in energy expenditure on long-term weight gain, assuming no adaptation in energy intake. The analyses indicate the following effects of a reduced level of energy expenditure in lean and obese women: (i) REE rises more slowly with increasing degrees of obesity due to a declining proportion of the more metabolically active fat-free mass; so, for the same positive energy balance, a significantly greater weight gain is expected for obese than for lean women before energy equilibrium is re-established; (ii) due to the greater energy density of adipose tissue, the time course of weight gain to achieve energy balance is longer for obese subjects: in general, this is approximately five years for lean and ten years for obese women; (iii) the magnitude of weight gain of lean women in response to a reduced energy expenditure of 200 to 800 kJ/day is only about 3 to 15 kg, amounts insufficient to explain severe obesity.

Formation of the long range Dpp morphogen gradient.

The TGF-β homolog Decapentaplegic (Dpp) acts as a secreted morphogen in the Drosophila wing disc, and spreads through the target tissue in order to form a long range concentration gradient. Despite extensive studies, the mechanism by which the Dpp gradient is formed remains controversial. Two opposing mechanisms have been proposed: receptor-mediated transcytosis (RMT) and restricted extracellular diffusion (RED). In these scenarios the receptor for Dpp plays different roles. In the RMT model it is essential for endocytosis, re-secretion, and thus transport of Dpp, whereas in the RED model it merely modulates Dpp distribution by binding it at the cell surface for internalization and subsequent degradation. Here we analyzed the effect of receptor mutant clones on the Dpp profile in quantitative mathematical models representing transport by either RMT or RED. We then, using novel genetic tools, experimentally monitored the actual Dpp gradient in wing discs containing receptor gain-of-function and loss-of-function clones. Gain-of-function clones reveal that Dpp binds in vivo strongly to the type I receptor Thick veins, but not to the type II receptor Punt. Importantly, results with the loss-of-function clones then refute the RMT model for Dpp gradient formation, while supporting the RED model in which the majority of Dpp is not bound to Thick veins. Together our results show that receptor-mediated transcytosis cannot account for Dpp gradient formation, and support restricted extracellular diffusion as the main mechanism for Dpp dispersal. The properties of this mechanism, in which only a minority of Dpp is receptor-bound, may facilitate long-range distribution.

What other sciences look like

In order to have references for discussing mathematical menus in political science, Ireview the most common types of mathematical formulae used in physics andchemistry, as well as some mathematical advances in economics. Several issues appearrelevant: variables should be well defined and measurable; the relationships betweenvariables may be non-linear; the direction of causality should be clearly identified andnot assumed on a priori grounds. On these bases, theoretically-driven equations onpolitical matters can be validated by empirical tests and can predict observablephenomena.

Driver scheduling problem modelling

The Drivers Scheduling Problem (DSP) consists of selecting a set of duties for vehicle drivers, for example buses, trains, plane or boat drivers or pilots, for the transportation of passengers or goods. This is a complex problem because it involves several constraints related to labour and company rules and can also present different evaluation criteria and objectives. Being able to develop an adequate model for this problem that can represent the real problem as close as possible is an important research area.The main objective of this research work is to present new mathematical models to the DSP problem that represent all the complexity of the drivers scheduling problem, and also demonstrate that the solutions of these models can be easily implemented in real situations. This issue has been recognized by several authors and as important problem in Public Transportation. The most well-known and general formulation for the DSP is a Set Partition/Set Covering Model (SPP/SCP). However, to a large extend these models simplify some of the specific business aspects and issues of real problems. This makes it difficult to use these models as automatic planning systems because the schedules obtained must be modified manually to be implemented in real situations. Based on extensive passenger transportation experience in bus companies in Portugal, we propose new alternative models to formulate the DSP problem. These models are also based on Set Partitioning/Covering Models; however, they take into account the bus operator issues and the perspective opinions and environment of the user.We follow the steps of the Operations Research Methodology which consist of: Identify the Problem; Understand the System; Formulate a Mathematical Model; Verify the Model; Select the Best Alternative; Present the Results of theAnalysis and Implement and Evaluate. All the processes are done with close participation and involvement of the final users from different transportation companies. The planner s opinion and main criticisms are used to improve the proposed model in a continuous enrichment process. The final objective is to have a model that can be incorporated into an information system to be used as an automatic tool to produce driver schedules. Therefore, the criteria for evaluating the models is the capacity to generate real and useful schedules that can be implemented without many manual adjustments or modifications. We have considered the following as measures of the quality of the model: simplicity, solution quality and applicability. We tested the alternative models with a set of real data obtained from several different transportation companies and analyzed the optimal schedules obtained with respect to the applicability of the solution to the real situation. To do this, the schedules were analyzed by the planners to determine their quality and applicability. The main result of this work is the proposition of new mathematical models for the DSP that better represent the realities of the passenger transportation operators and lead to better schedules that can be implemented directly in real situations.

Accelerated tumor invasion under non-isotropic cell dispersal in glioblastomas

Glioblastomas are highly diffuse, malignant tumors that have so far evaded clinical treatment. The strongly invasive behavior of cells in these tumors makes them very resistant to treatment, and for this reason both experimental and theoretical efforts have been directed toward understanding the spatiotemporal pattern of tumor spreading. Although usual models assume a standard diffusion behavior, recent experiments with cell cultures indicate that cells tend to move in directions close to that of glioblastoma invasion, thus indicating that a biasedrandom walk model may be much more appropriate. Here we show analytically that, for realistic parameter values, the speeds predicted by biased dispersal are consistent with experimentally measured data. We also find that models beyond reaction–diffusion–advection equations are necessary to capture this substantial effect of biased dispersal on glioblastoma spread

Dispersal probability distributions and the wave-front speed problem

The speed and width of front solutions to reaction-dispersal models are analyzed both analytically and numerically. We perform our analysis for Laplace and Gaussian distribution kernels, both for delayed and nondelayed models. The results are discussed in terms of the characteristic parameters of the models