Efficient simulation of stochastic chemical kinetics with the Stochastic Bulirsch-Stoer extrapolation method

Background Biochemical systems with relatively low numbers of components must be simulated stochastically in order to capture their inherent noise. Although there has recently been considerable work on discrete stochastic solvers, there is still a need for numerical methods that are both fast and accurate. The Bulirsch-Stoer method is an established method for solving ordinary differential equations that possesses both of these qualities. Results In this paper, we present the Stochastic Bulirsch-Stoer method, a new numerical method for simulating discrete chemical reaction systems, inspired by its deterministic counterpart. It is able to achieve an excellent efficiency due to the fact that it is based on an approach with high deterministic order, allowing for larger stepsizes and leading to fast simulations. We compare it to the Euler τ-leap, as well as two more recent τ-leap methods, on a number of example problems, and find that as well as being very accurate, our method is the most robust, in terms of efficiency, of all the methods considered in this paper. The problems it is most suited for are those with increased populations that would be too slow to simulate using Gillespie’s stochastic simulation algorithm. For such problems, it is likely to achieve higher weak order in the moments. Conclusions The Stochastic Bulirsch-Stoer method is a novel stochastic solver that can be used for fast and accurate simulations. Crucially, compared to other similar methods, it better retains its high accuracy when the timesteps are increased. Thus the Stochastic Bulirsch-Stoer method is both computationally efficient and robust. These are key properties for any stochastic numerical method, as they must typically run many thousands of simulations.

Impacts of tillage and application methods on atrazine and alachlor losses from upland fields: Research paper

The effects of tillage practises and the methods of chemical application on atrazine and alachlor losses through run-off were evaluated for five treatments: conservation (untilled) and surface (US), disk and surface, plow and surface, disk and preplant-incorporated, and plow and preplant-incorporated treatments. A rainfall simulator was used to create 63.5 mm h^-1 of rainfall for 60 min and 127 mm h^-1 for 15 min. Rainfall simulation occurred 24-36 h after chemical application. There was no significant difference in the run-off volume among the treatments but the untilled treatment significantly reduced erosion loss. The untilled treatments had the highest herbicide concentration and the disk treatments were higher than the plow treatments. The surface treatments showed a higher concentration than the incorporated treatments. The concentration of herbicides in the water decreased with time. Among the experimental sites, the one with sandy loam soil produced the greatest losses, both in terms of the run-off volume and herbicide loss. The US treatments had the highest loss and the herbicide incorporation treatments had smaller losses through run-off as the residue cover was effective in preventing herbicide losses. Incorporation might be a favorable method of herbicide application to reduce the herbicide losses by run-off.

An overview of land application of pig effluent-P using soil P chemistry and mass balance calculations.

Attention is directed at land application of piggery effluent (containing urine, faeces, water, and wasted feed) as a potential source of water resource contamination with phosphorus (P). This paper summarises P-related properties of soil from 0-0.05 m depth at 11 piggery effluent application sites, in order to explore the impact that effluent application has had on the potential for run-off transport of P. The sites investigated were situated on Alfisol, Mollisol, Vertisol, and Spodosol soils in areas that received effluent for 1.5-30 years (estimated effluent-P applications of 100-310000 kg P/ha in total). Total (PT), bicarbonate extractable (PB), and soluble P forms were determined for the soil (0-0.05 m) at paired effluent and no-effluent sites, as well as texture, oxalate-extractable Fe and Al, organic carbon, and pH. All forms of soil P at 0-0.05 m depth increased with effluent application (PB at effluent sites was 1.7-15 times that at no-effluent sites) at 10 of the 11 sites. Increases in PB were strongly related to net P applications (regression analysis of log values for 7 sites with complete data sets: 82.6 % of variance accounted for, p <0.01). Effluent irrigation tended to increase the proportion of soil PT in dilute CaCl2-extractable forms (PTC: effluent average 2.0 %; no-effluent average 0.6%). The proportion of PTC in non-molybdate reactive forms (centrifuged supernatant) decreased (no-effluent average, 46.4 %; effluent average, 13.7 %). Anaerobic lagoon effluent did not reliably acidify soil, since no consistent relationship was observed for pH with effluent application. Soil organic carbon was increased in most of the effluent areas relative to the no-effluent areas. The four effluent areas where organic carbon was reduced had undergone intensive cultivation and cropping. Current effluent management at many of the piggeries failed to maximise the potential for waste P recapture. Ten of the case-study effluent application areas have received effluent-P in excess of crop uptake. While this may not represent a significant risk of leaching where sorption retains P, it has increased the risk of transport of P by run-off. Where such sites are close to surface water, run-off P loads should be managed.

Simulation of the supply capability of sheep on the Mitchell grass downs of north west Queensland

Sheep in western Queensland have been predominantly reared for wool. When wool prices became depressed interest in the sheep meat industry, increased. For north west Queensland producers, opportunities may exist to participate in live sheep and meat export to Asia. A simulation model was developed to determine whether this sheep producing area has the capability to provide sufficient numbers of sheep under variable climatic conditions while sustaining the land resources. Maximum capacity for sustainability of resources (as described by stock numbers) was derived from an in-depth study of the agricultural and pastoral potential of Queensland. Decades of sheep production and climatic data spanning differing seasonal conditions were collated for analysis. A ruminant biology model adapted from Grazplan was used to simulate pregnancy rate. Empirical equations predict mortalities, marking rates, and weight characteristics of sheep of various ages from simple climatic measures, stocking rate and reproductive status. The initial age structure of flocks was determined by running the model for several years with historical climatic conditions. Drought management strategies such as selling a proportion of wethers progressively down to two-tooth and oldest ewes were incorporated. Management decisions such as time of joining, age at which ewes were cast-for-age, wether turn-off age and turning-off rate of lambs vary with geographical area and can be specified at run time. The model is run for sequences of climatic conditions generated stochastically from distributions based on historical climatic data correlated in some instances. The model highlights the difficulties of sustaining a consistent supply of sheep under variable climatic conditions.

Management of mirids and stinkbugs in Bollgard II

The 'insurance' sprays by the growers run the risk of flaring secondary pests such as whitefly, mites and aphids. Although recent research on mirids has addressed some aspects of this issue, further effort is needed to refine mirid management.

Judah P. Benjamin collection, undated, 1838-1942.

The collection consists of 25 letters written by Benjamin between 1838 and 1881 on a variety of subjects, four Confederate notes and two bonds bearing his picture, miscellaneous items about Benjamin (1893-1942), nine issues of the Congressional globe with speeches by Benjamin, as well as separate copies of his printed speeches, and a photostatic copy of the "Diary of Events" (400 pp.) kept by Benjamin, the original of which is in the Library of Congress (1862-1864).

Co2 emissions and income trajectory in Australia: The role of technological change

This study investigates the relationship between per capita carbon dioxide (CO2) emissions and per capita GDP in Australia, while controlling for technological state as measured by multifactor productivity and export of black coal. Although technological progress seems to play a critical role in achieving long term goals of CO2 reduction and economic growth, empirical studies have often considered time trend to proxy technological change. However, as discoveries and diffusion of new technologies may not progress smoothly with time, the assumption of a deterministic technological progress may be incorrect in the long run. The use of multifactor productivity as a measure of technological state, therefore, overcomes the limitations and provides practical policy directions. This study uses recently developed bound-testing approach, which is complemented by Johansen- Juselius maximum likelihood approach and a reasonably large sample size to investigate the cointegration relationship. Both of the techniques suggest that cointegration relationship exists among the variables. The long-run and short-run coefficients of CO2 emissions function is estimated using ARDL approach. The empirical findings in the study show evidence of the existence of Environmental Kuznets Curve type relationship for per capita CO2 emissions in the Australian context. The technology as measured by the multifactor productivity, however, is not found as an influencing variable in emissionsincome trajectory.

Impact of wildfires and floods on property values: A before and after analysis

One of the most evident casualties of a natural disaster is the property market. The private and social costs from such events run into millions of dollars. In this paper, we use a unique dataset to examine the impact on residential house prices affected by natural disasters using a hedonic property (HP) values approach. For this purpose, we use data before and after a wildfire and floods from Rockhampton in central Queensland, Australia. The data is unique because one suburb was affected by wildfires and another was affected by floods. For the analysis, three suburbs namely Frenchville, Park Avenue and Norman Gardens are used. Frenchville was significantly affected by wildfires in the latter part of 2009 and to a lesser extent in 2012, while Park Avenue was affected by floods at the end of 2010, January 2011–2013. Norman Gardens, which was relatively unaffected, is used as a control site. This enables us to examine the before and after effects on property values in the three suburbs. The results confirm that soon after a natural disaster property prices in affected areas decrease even though the large majority of individual houses remain unaffected. Furthermore, the results indicate that the largely unaffected suburb may gain immediately after a natural disaster but this gain may disappear if natural disasters continue to occur in the area/region due to the stigma created. The results have several important policy decision and welfare implications which are briefly discussed in the paper.

Synthetic investigation on the building of ring A of steroids

Methyl 7-keto-1,2,3,4,4a,5,6,7-octahydronaphthoate (Va) has been prepared by the reduction of 7-methoxy-1,2,3,4-tetrahydronaphthoic acid (III) with lithium and ammonia followed by hydrolysis of the enol ether, esterification and migration of the double bond. Alkylation of Va has led to the substitution at the expected 8-position. Methyl 4-keto-7-methoxy-1,2,3,4-tetrahydronaphthoate (X), an intermediate in the preparation of III, has been converted into methyl 3-methyl-3-cyano-4-keto-7-methoxy-1,2,3,4-tetrahydronaphthoate (XIII).

Synthesis and Identification of a new class of antileukemic agents containing 2-(arylcarboxamide)-(S)-6-amino-4,5,6,7-tetrahydrobenzod]thiazole

Recently we have reported the effect of (S)-6-aryl urea/thiourea substituted-2-amino-4,5,6,7-tetrahydrobenzod]thiazole derivatives as potent anti-leukemic agents. To elucidate further the Structure Activity Relationship (SAR) studies on the anti-leukemic activity of (S)-2,6-diamino-4,5,6,7 tetrahydrobenzod]thiazole moiety, a series of 2-arlycarboxamide substituted-(S)-6-amino-4,5,6,7-tetrahydrobenzod]thiazole were designed, synthesized and evaluated for their anti-leukemic activity by trypan blue exclusion, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), lactate dehydrogenase (LDH) assays and cell cycle analysis. Results suggest that the position, number and bulkiness of the substituent on the phenyl ring of aryl carboxamide moiety at 2nd position of 6-amino-4,5,6,7-tetrhydrobenzod]thiazole play a key role in inhibiting the proliferation of leukemia cells. Compounds with ortho substitution showed poor activity and with meta and para substitution showed good activity. (C) 2010 Elsevier Masson SAS. All rights reserved.

Income inequality in the process of economic development : An empirical approach

This thesis studies the effect of income inequality on economic growth. This is done by analyzing panel data from several countries with both short and long time dimensions of the data. Two of the chapters study the direct effect of inequality on growth, and one chapter also looks at the possible indirect effect of inequality on growth by assessing the effect of inequality on savings. In Chapter two, the effect of inequality on growth is studied by using a panel of 70 countries and a new EHII2008 inequality measure. Chapter contributes on two problems that panel econometric studies on the economic effect of inequality have recently encountered: the comparability problem associated with the commonly used Deininger and Squire s Gini index, and the problem relating to the estimation of group-related elasticities in panel data. In this study, a simple way to 'bypass' vagueness related to the use of parametric methods to estimate group-related parameters is presented. The idea is to estimate the group-related elasticities implicitly using a set of group-related instrumental variables. The estimation results with new data and method indicate that the relationship between income inequality and growth is likely to be non-linear. Chapter three incorporates the EHII2.1 inequality measure and a panel with annual time series observations from 38 countries to test the existence of long-run equilibrium relation(s) between inequality and the level of GDP. Panel unit root tests indicate that both the logarithmic EHII2.1 inequality measure and the logarithmic GDP per capita series are I(1) nonstationary processes. They are also found to be cointegrated of order one, which implies that there is a long-run equilibrium relation between them. The long-run growth elasticity of inequality is found to be negative in the middle-income and rich economies, but the results for poor economies are inconclusive. In the fourth Chapter, macroeconomic data on nine developed economies spanning across four decades starting from the year 1960 is used to study the effect of the changes in the top income share to national and private savings. The income share of the top 1 % of population is used as proxy for the distribution of income. The effect of inequality on private savings is found to be positive in the Nordic and Central-European countries, but for the Anglo-Saxon countries the direction of the effect (positive vs. negative) remains somewhat ambiguous. Inequality is found to have an effect national savings only in the Nordic countries, where it is positive.

Dexamethasone negatively regulates phenobarbitone-activated transcription but synergistically enhances cytoplasmic levels of cytochrome P-450b/e messenger RNA

Dexamethasone has a potentiating effect on phenobarbitone mediated induction of cytochrome P-450b + e mRNAs in adult rat liver. However, the glucocorticoid inhibits phenobarbitone-activated transcription of cytochrome P-450b + e mRNAs by 60-70%. This inhibitory effect is evident in run-off transcription of the endogenous genes as well as in the transcription of an added cloned gene fragment. Dexamethasone inhibits the phenobarbitone-mediated increase in the binding of a transcription factor(s) to the upstream region of the gene as evidenced by gel retardation and Southwestern blot analysis. The glucocorticoid does not stabilize the phenobarbitone-induced polyribosomal cytochrome P-450b + e mRNAs but appears to stabilize the nuclear transcripts. It is proposed that a negative element may mediate the action of dexamethasone at the level of nuclear transcription and stabilization of the nuclear transcript may account for the potentiating effect of the glucocorticoid on phenobarbitone-mediated increase in cytochrome P-450b + e mRNAs in the cytoplasm of the adult rat liver. However, the cytochrome P-450b protein levels are slightly lower in phenobarbitone + dexamethasone treatment than in phenobarbitone-treated liver microsomes.

Influence of metalaxyl on Glomus fasciculatum associated with wheat (Triticum aestivum L)

Soil incorporation of metalaxyl [methyl N-(2-methoxyacetyl)-N-(2,6,xylyl)-DL-alaninate] significantly enhanced root colonization of the vesicular-arbuscular (VA) mycorrhizal fungi Glomus fasciculatum associated with wheat. The stimulatory response of VA mycorrhizal fungi to low concentration of metalaxyl resulted in increased plant biomass production, nutrient uptake and grain yield of wheat. However, higher concentrations of metalaxyl, particularly 2.5 ppm of metalaxyl affected the mycorrhizal infection and seed yield of wheat, Addition of urban compost to an extent ameliorated the toxic effect of fungicide on VA mycorrhizal colonization, plant growth and yield of wheat when compared to unamended soil.

Mechanisms of summer intraseasonal sea surface temperature oscillations in the Bay of Bengal

Intraseasonal variations (ISV) of sea surface temperature (SST) in the Bay of Bengal (BoB) is highest in its northwestern part. An Indian Ocean model forced by QuikSCAT winds and climatological river discharge (QR run) reproduces ISV of SST, albeit with weaker magnitude. Air-sea fluxes, in the presence of a shallow mixed layer, efficiently effect intraseasonal SST fluctuations. Warming during intraseasonal events is smaller (<1°C) for June - July period and larger (1.5° to 2°C) during September, the latter due to a thinner mixed layer. To examine the effect of salinity on ISV, the model was run by artificially increasing the salinity (NORR run) and by decreasing it (MAHA10 run). In NORR, both rainfall and river discharge were switched off and in MAHA10 the discharge by river Mahanadi was increased tenfold. The spatial pattern of ISV as well as its periodicity was similar in QR, NORR and MAHA10. The ISV was stronger in NORR and weaker in MAHA10, compared to QR. In NORR, both intraseasonal warming and cooling were higher than in QR, the former due to reduced air-sea heat loss as the mean SST was lower, and the latter due to enhanced subsurface processes resulting from weaker stratification. In MAHA10, both warming and cooling were lower than in QR, the former due to higher air-sea heat loss owing to higher mean SST, and the latter due to weak subsurface processes resulting from stronger stratification. These model experiments suggest that salinity effects are crucial in determining amplitudes of intraseasonal SST variations in the BoB.

Analysis of grating inscribed micro-cantilever for high resolution AFM probe

We present a mathematical modelling and analysis of reflection grating etched Si AFM cantilever deflections under different loading conditions. A simple analysis of the effect of grating structures on cantilever deflection is carried out with emphasis on optimizing the beam and gratings such that maximum amount of diffracted light remains within the detector area.