Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

A hypothesis for explaining single outbreaks (like the Black Death in European cities) of vector-borne infections

We propose a mechanism by which single outbreaks of vector-borne infections can happen even when the value of the basic reproduction number, R(o), of the infection is below one. With this hypothesis we have shown that dynamical models simulations demonstrate that the arrival of a relatively small (with respect to the host population) number of infected vectors can trigger a short-lived epidemic but with a huge number of cases. These episodes are characterized by a sudden outbreak in a previously virgin area that last from weeks to a few months, and then disappear without leaving vestiges. The hypothesis proposed in this paper to explain those single outbreaks of vector-borne infections, even when total basic reproduction number, Ro, is less than one (which explain the fact that those infections fail to establish themselves at endemic levels), is that the vector-to-host component of Ro is greater than one and that a sufficient amount of infected vectors are imported to the vulnerable area, triggering the outbreak. We tested the hypothesis by performing numerical simulations that reproduce the observed outbreaks of chikungunya in Italy in 2007 and the plague in Florence in 1348. The theory proposed provides an explanation for isolated outbreaks of vector-borne infections, ways to calculate the size of those outbreaks from the number of infected vectors arriving in the affected areas. Given the ever-increasing worldwide transportation network, providing a high degree of mobility from endemic to virgin areas, the proposed mechanism may have important implications for public health planning. (C) 2009 Elsevier Ltd. All rights reserved.

Cost-effectiveness analysis of a hypothetical hepatitis C vaccine compared to antiviral therapy

We propose a mathematical model to simulate the dynamics of hepatitis C virus (HCV) infection in the state of Sao Paulo, Brazil. We assumed that a hypothetical vaccine, which cost was taken to be the initial cost of the vaccine against hepatitis B exists and it is introduced in the model. We computed its cost-effectiveness compared with the anti-HCV therapy. The calculated basic reproduction number was 1.20. The model predicts that without intervention a steady state exists with an HCV prevalence of 3%, in agreement with the Current epidemiological data. Starting from this steady state three interventions were simulated: indiscriminate vaccination, selective vaccination and anti-HCV therapy. Selective vaccination proved to be the strategy with the best cost-effectiveness ratio, followed by indiscriminate vaccination and anti-HCV therapy.

The risk of chikungunya fever in a dengue-endemic area

Background. Chikungunya, an alphavirus of the Togaviridae family, causes a febrile disease transmitted to humans by the bite of infected Aedes mosquitoes. This infection is reaching endemic levels in many Southeast Asian countries. Symptoms include sudden onset of fever, chills, headache, nausea, vomiting, joint pain with or without swelling, low back pain, and rash. According to the World Health Organization, there are 2 billion people living in Aedes-infested areas. In addition, traveling to these areas is popular, making the potential risk of infections transmitted by the bite of infected Aedes mosquitoes very high. Methods. We proposed a mathematical model to estimate the risk of acquiring chikungunya fever in an Aedes-infested area by taking the prevalence of dengue fever into account. The basic reproduction number for chikungunya fever R-0chik can be written as a function of the basic reproduction number of dengue R-0dengue by calculating the ratio R-0chik/R-0dengue. From R-0chik, we estimated the force of infection and the risk of acquiring the disease both for local residents of a dengue-endemic area and for travelers to this area. Results. We calculated that R-0chik is 64.4% that of R-0dengue. The model was applied to a hypothetical situation, namely, estimating the individual risk of acquiring chikungunya fever in a dengue-endemic area, both for local inhabitants (22% in steady state) and for visiting travelers (from 0.31% to 1.23% depending on the time spent in the area). Conclusions. The method proposed based on the output of a dynamical model is innovative and provided an estimation of the risk of infection, both for local inhabitants and for visiting travelers.

FOR DEBATE A hypothesis for the 2007 dengue outbreak in Singapore

A previous mathematical model explaining dengue in Singapore predicted a reasonable outbreak of about 6500 cases for 2006 and a very mild outbreak with about 2000 cases for 2007. However, only 3051 cases were reported in 2006 while more than 7800 were reported in the first 44 weeks of 2007. We hypothesized that the combination of haze with other local sources of particulate matter had a significant impact on mosquito life expectancy, significantly increasing their mortality rate. To test the hypothesis a mathematical model based on the reproduction number of dengue fever and aimed at comparing the impact of several possible alternative control strategies was proposed. This model also aimed at contributing to the understanding of the causes of dengue resurgence in Singapore in the last decade. The model`s simulation demonstrated that an increase in mosquito mortality in 2006 and either a reduction in mortality or an increase in the carrying capacity of mosquitoes in 2007 explained the patterned observed in Singapore. Based on the model`s simulation we concluded that the fewer than expected number of dengue cases in Singapore in 2006 was caused by an increase in mosquito mortality due to the disproportionate haze affecting the country that year and that particularly favourable environmental conditions in 2007 propitiated mosquitoes with a lower mortality rate, which explains the greater than expected number of dengue cases in 2007. Whether our hypothesis is plausible or not should be debated further.

Modeling the Dynamics of Viral Evolution Considering Competition Within Individual Hosts and at Population Level: The Effects of Treatment

We consider two viral strains competing against each other within individual hosts (at cellular level) and at population level (for infecting hosts) by studying two cases. In the first case, the strains do not mutate into each other. In this case, we found that each individual in the population can be infected by only one strain and that co-existence in the population is possible only when the strain that has the greater basic intracellular reproduction number, R (0c) , has the smaller population number R (0p) . Treatment against the one strain shifts the population equilibrium toward the other strain in a complicated way (see Appendix B). In the second case, we assume that the strain that has the greater intracellular number R (0c) can mutate into the other strain. In this case, individual hosts can be simultaneously infected by both strains (co-existence within the host). Treatment shifts the prevalence of the two strains within the hosts, depending on the mortality induced by the treatment, which is, in turn, dependent upon the doses given to each individual. The relative proportions of the strains at the population level, under treatment, depend both on the relative proportions within the hosts (which is determined by the dosage of treatment) and on the number of individuals treated per unit time, that is, the rate of treatment. Implications for cases of real diseases are briefly discussed.

Estimation of R(0) from the initial phase of an outbreak of a vector-borne infection

The magnitude of the basic reproduction ratio R(0) of an epidemic can be estimated in several ways, namely, from the final size of the epidemic, from the average age at first infection, or from the initial growth phase of the outbreak. In this paper, we discuss this last method for estimating R(0) for vector-borne infections. Implicit in these models is the assumption that there is an exponential phase of the outbreaks, which implies that in all cases R(0) > 1. We demonstrate that an outbreak is possible, even in cases where R(0) is less than one, provided that the vector-to-human component of R(0) is greater than one and that a certain number of infected vectors are introduced into the affected population. This theory is applied to two real epidemiological dengue situations in the southeastern part of Brazil, one where R(0) is less than one, and other one where R(0) is greater than one. In both cases, the model mirrors the real situations with reasonable accuracy.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

Mass transfer correlations for rotating drum bioreactors

Evaporative cooling is extremely important for large-scale operation of rotating drum bioreactors (RDBs). Outlet water vapour concentrations were measured for a RDB containing wet wheat bran with the aim of determining the mass transfer coefficient for evaporation from the bran bed to the headspace. Mass transfer was expressed as the mass transfer coefficient times the area for transfer per unit volume of void space in the drum. Values of ka' were determined under combinations of aeration superficial velocities ranging from 0.006 to 0.017 ms(-1) and rotation rates ranging from 0 to 9 rpm. Mass transfer coefficients were evaluated using a variety of residence time distributions (RTDs) for flow in the gas phase including plug flow and well-mixed and a Central Jet RTD based on RTD studies. If plug flow is assumed, the degree of holdup at low effective Peclet (Pe(eff)) numbers gives an apparent under-estimate of ka' compared with empirical correlations. Values of ka' calculated using the Central Jet RTD agree well with values of ka' from literature correlations. There was a linear relationship between ka' and effective Peclet number: ka' = 2.32 x 10(-3) Pe(eff). (C) 2002 Elsevier Science B.V. All rights reserved.

Dengue and the risk of urban yellow fever reintroduction in São Paulo State, Brazil

OBJECTIVE: To propose a mathematical method for the estimation of the Basic Reproduction Number, R0, of urban yellow fever in a dengue-infested area. METHODS: The method is based on the assumption that, as the same vector (Aedes aegypti) causes both infections, all the quantities related to the mosquito, estimated from the initial phase of dengue epidemic, could be applied to yellow fever dynamics. It is demonstrated that R0 for yellow fever is, on average, 43% lower than that for dengue. This difference is due to the longer dengue viremia and its shorter extrinsic incubation period. RESULTS: In this study the analysis was expanded to the epidemiological situation of dengue in São Paulo in the year 2001. The total number of dengue cases increased from 3,582 in 2000 to 51,348 in 2001. It was then calculated R0 for yellow fever for every city which have shown R0 of dengue greater than 1. It was also estimated the total number of unprotected people living in highly risky areas for urban yellow fever. CONCLUSIONS: Currently there is a great number of non-vaccinated people living in Aedes aegypti infested area in the state of São Paulo.

Imported and autochthonous cases in the dynamics of dengue epidemics in Brazil

OBJECTIVE: To estimate the basic reproduction number (R0) of dengue fever including both imported and autochthonous cases. METHODS: The study was conducted based on epidemiological data of the 2003 dengue epidemic in Brasília, Brazil. The basic reproduction number is estimated from the epidemic curve, fitting linearly the increase of initial cases. Aiming at simulating an epidemic with both autochthonous and imported cases, a "susceptible-infectious-resistant" compartmental model was designed, in which the imported cases were considered as an external forcing. The ratio between R0 of imported versus autochthonous cases was used as an estimator of real R0. RESULTS: The comparison of both reproduction numbers (only autochthonous versus all cases) showed that considering all cases as autochthonous yielded a R0 above one, although the real R0 was below one. The same results were seen when the method was applied on simulated epidemics with fixed R0. This method was also compared to some previous proposed methods by other authors and showed that the latter underestimated R0 values. CONCLUSIONS: It was shown that the inclusion of both imported and autochthonous cases is crucial for the modeling of the epidemic dynamics, and thus provides critical information for decision makers in charge of prevention and control of this disease.

Mathematical model for HIV dynamics in HIV-specific helper cells

In this paper we study a delay mathematical model for the dynamics of HIV in HIV-specific CD4 + T helper cells. We modify the model presented by Roy and Wodarz in 2012, where the HIV dynamics is studied, considering a single CD4 + T cell population. Non-specific helper cells are included as alternative target cell population, to account for macrophages and dendritic cells. In this paper, we include two types of delay: (1) a latent period between the time target cells are contacted by the virus particles and the time the virions enter the cells and; (2) virus production period for new virions to be produced within and released from the infected cells. We compute the reproduction number of the model, R0, and the local stability of the disease free equilibrium and of the endemic equilibrium. We find that for values of R0<1, the model approaches asymptotically the disease free equilibrium. For values of R0>1, the model approximates asymptotically the endemic equilibrium. We observe numerically the phenomenon of backward bifurcation for values of R0⪅1. This statement will be proved in future work. We also vary the values of the latent period and the production period of infected cells and free virus. We conclude that increasing these values translates in a decrease of the reproduction number. Thus, a good strategy to control the HIV virus should focus on drugs to prolong the latent period and/or slow down the virus production. These results suggest that the model is mathematically and epidemiologically well-posed.

New findings on the dynamics of HIV and TB coinfection models

In this paper we study a model for HIV and TB coinfection. We consider the integer order and the fractional order versions of the model. Let α∈[0.78,1.0] be the order of the fractional derivative, then the integer order model is obtained for α=1.0. The model includes vertical transmission for HIV and treatment for both diseases. We compute the reproduction number of the integer order model and HIV and TB submodels, and the stability of the disease free equilibrium. We sketch the bifurcation diagrams of the integer order model, for variation of the average number of sexual partners per person and per unit time, and the tuberculosis transmission rate. We analyze numerical results of the fractional order model for different values of α, including α=1. The results show distinct types of transients, for variation of α. Moreover, we speculate, from observation of the numerical results, that the order of the fractional derivative may behave as a bifurcation parameter for the model. We conclude that the dynamics of the integer and the fractional order versions of the model are very rich and that together these versions may provide a better understanding of the dynamics of HIV and TB coinfection.

Dynamics of the 2006/2007 dengue outbreak in Brazil

We analyzed dengue incidence in the period between October 2006-July 2007 of 146 cities around the country were Larval Index Rapid Assay (LIRA) surveillance was carried out in October 2006. Of these, we chosen 61 cities that had 500 or more cases reported during this period. We calculated the incidence coefficient, the force of infection (») and the basic reproduction number (R0) of dengue in those 61 cities and correlated those variables with the LIRA. We concluded that » and R0 are more associated with the number of cases than LIRA. In addition, the average R0 for the 2006/2007 dengue season was almost as high as that calculated for the 2001/2002 season, the worst in Brazilian history.

Statistical Analysis of an SEIR Epidemic Model

This thesis was focussed on statistical analysis methods and proposes the use of Bayesian inference to extract information contained in experimental data by estimating Ebola model parameters. The model is a system of differential equations expressing the behavior and dynamics of Ebola. Two sets of data (onset and death data) were both used to estimate parameters, which has not been done by previous researchers in (Chowell, 2004). To be able to use both data, a new version of the model has been built. Model parameters have been estimated and then used to calculate the basic reproduction number and to study the disease-free equilibrium. Estimates of the parameters were useful to determine how well the model fits the data and how good estimates were, in terms of the information they provided about the possible relationship between variables. The solution showed that Ebola model fits the observed onset data at 98.95% and the observed death data at 93.6%. Since Bayesian inference can not be performed analytically, the Markov chain Monte Carlo approach has been used to generate samples from the posterior distribution over parameters. Samples have been used to check the accuracy of the model and other characteristics of the target posteriors.