Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Selecting methods to solve multi-well master equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.

An investigation of magnetic antennas for ground penetrating radar

For ground penetrating radar (GPR), smaller antennas would provide considerable practical advantages. Some of which are: portability; ease of use; and higher spatial sampling. A theoretical comparison of the fundamental limits of a small electric field antenna and a small magnetic field antenna shows that the minimum Q constraints are identical. Furthermore, it is shown that only the small magnetic loop antenna can be constructed to approach, arbitrarily closely, the fundamental minimum Q limit. This is achieved with the addition of a high permeability material which reduces energy stored in the magnetic fields. This is of special interest to some GPR applications. For example, applications requiring synthetic aperture data collection would benefit from the increased spatial sampling offered by electrically smaller antennas. Low frequency applications may also benefit, in terms of reduced antenna dimensions, by the use of electrically small antennas. Under these circumstances, a magnetic type antenna should be considered in preference to the typical electric field antenna. Numerical modeling data supports this assertion.

Integrated control of nitrate recirculation and external carbon addition in a predenitrification system

The integrated control of nitrate recirculation and external carbon addition in a predenitrification biological wastewater treatment system is studied. The proposed control structure consists of four feedback control loops, which manipulate the nitrate recirculation and the carbon dosage flows in a highly coordinated manner such that the consumption of external carbon is minimised while the nitrate discharge limits (based on both grab and composite samples) are met. The control system requires the measurement of the nitrate concentrations at the end of both the anoxic and the aerobic zones. Distinct from ordinary control systems, which typically minimise the variation in the controlled variables, the proposed control system essentially maximises the diurnal variation of the effluent nitrate concentration and through this maximises the use of influent COD for denitrification, thus minimising the requirement for external carbon source. Simulation studies using a commonly accepted simulation benchmark show that the controlled system consistently achieves the designated effluent quality with minimum costs.

Antibiotic desensitization in adults with cystic fibrosis

Objective: Allergic reactions to antibiotics occur in up to 30% of patients with cystic fibrosis (CF). Repeated antibiotic exposure and immune hyper-responsiveness increase the risk of allergic reactions and may limit antibiotic choice. Desensitization may allow the successful administration of an antibiotic despite previous allergy. We aimed to determine the success of antibiotic desensitization in patients with CF in an adult CF unit over a 7-year period. Methodology: A retrospective medical record review was performed on the 19 patients who had undergone antibiotic desensitization procedures. Data collected included drug allergy and intolerance profiles, nature of allergies, and the outcome of desensitization procedures. Desensitization procedures were performed in a ward setting according to published methods. Results: Nineteen patients (13 females) reported 62 drug allergies with a mean of 3.3 per patient. Of the 71 desensitization procedures undergone by this group, 54 (76%) were successful. Fifteen of the 19 patients were allergic to two or more beta-lactam antibiotics. Over half of the patients were desensitized to more than one antibiotic. Nine different antibiotics were used in 31 different patient/drug combinations. A successful outcome was achieved in 18/31 (58%) combinations, with three requiring treatment for mild allergic reactions. Allergic reactions caused drug cessation in a total of 19 patient/drug combinations (three after initial successful desensitization and full courses of antibiotics). Over 50% of these reactions occurred on day 1. Desensitization failures were more common in patients with well-documented allergic reactions to a specific drug. Conclusion: This study demonstrates that multiple antibiotic allergies are common in adults with CE Cross-reactivity between beta-lactam antibiotics may limit antibiotic choice for the treatment of pulmonary exacerbations. Antibiotic desensitization allows safe and successful treatment in the ward setting of many patients with previous allergies to an antibiotic. In many patients symptoms of allergy still occur and result in cessation of the antibiotics. Use of corticosteroids and antihistamines may improve the success rate of desensitization procedures.

Gene flow in great bustard populations across the Strait of Gibraltar as elucidated from excremental PCR and mtDNA sequencing

Recent advances in molecular biology have made it possible to use the trace amounts of DNA in faeces to non-invasively sample endangered species for genetic studies. Here we use faeces as a source of DNA and mtDNA sequence data to elucidate the relationship among Spanish and Moroccan populations of great bustards. 834 bp of combined control region and cytochrome-b mtDNA fragments revealed four variable sites that defined seven closely related haplotypes in 54 individuals. Morocco was fixed for a single mtDNA haplotype that occurs at moderate frequency (28%) in Spain. We could not differentiate among the sampled Spanish populations of Caceres and Andalucia but these combined populations were differentiated from the Moroccan population. Estimates of gene flow (Nm = 0.82) are consistent with extensive observations on the southern Iberian peninsular indicating that few individuals fly across the Strait of Gibraltar. We demonstrate that both this sea barrier and mountain barriers in Spain limit dispersal among adjacent great bustard populations to a similar extent. The Moroccan population is of high ornithological significance as it holds the only population of great bustards in Africa. This population is critically small and genetic and observational data indicate that it is unlikely to be recolonised via immigration from Spain should it be extirpated. In light of the evidence presented here it deserves the maximum level of protection.

Surface diffusion of strongly adsorbing vapors in activated carbon by a differential permeation method

Conventional methods to determine surface diffusion of adsorbed molecules are proven to be inadequate for strongly adsorbing vapors on activated carbon. Knudsen diffusion permeability (B-k) for strongly adsorbing vapors cannot be directly estimated from that of inert gases such as helium. In this paper three models are considered to elucidate the mechanism of surface diffusion in activated carbon. The transport mechanism in all three models is a combination of Knudsen diffusion, viscous flow and surface diffusion. The collision reflection factor f (which is the fraction of molecules undergoing collision to the solid surface over reflection from the surface) of the Knudsen diffusivity is assumed to be a function of loading. It was found to be 1.79 in the limit of zero loading, and decreases as loading increases. The surface diffusion permeability increases sharply at very low pressures and then starts to decrease after it has reached a maximum (B(mum)s) at a threshold pressure. The initial rapid increase in the total permeability is mainly attributed to surface diffusion. Interestingly the B(mum)s for all adsorbates appear at the same volumetric adsorbed phase concentration, suggesting that the volume of adsorbed molecules may play an important role in the surface diffusion mechanism in activated carbon. (C) 2003 Elsevier Ltd. All rights reserved.

Experimental investigation of criteria for continuous variable entanglement

We generate a pair of entangled beams from the interference of two amplitude squeezed beams. The entanglement is quantified in terms of EPR paradox and inseparability criteria, with both results clearly beating the standard quantum limit. We experimentally analyze the effect of decoherence on each criterion and demonstrate qualitative differences. We also characterize the number of required and excess photons present in the entangled beams and provide contour plots of the efficacy of quantum information protocols in terms of these variables.

Pseudo-randomness of round-off errors in discretized linear maps on the plane

We analyze the sequences of round-off errors of the orbits of a discretized planar rotation, from a probabilistic angle. It was shown [Bosio & Vivaldi, 2000] that for a dense set of parameters, the discretized map can be embedded into an expanding p-adic dynamical system, which serves as a source of deterministic randomness. For each parameter value, these systems can generate infinitely many distinct pseudo-random sequences over a finite alphabet, whose average period is conjectured to grow exponentially with the bit-length of the initial condition (the seed). We study some properties of these symbolic sequences, deriving a central limit theorem for the deviations between round-off and exact orbits, and obtain bounds concerning repetitions of words. We also explore some asymptotic problems computationally, verifying, among other things, that the occurrence of words of a given length is consistent with that of an abstract Bernoulli sequence.

Inova????es na Gest??o P??blica Federal: balan??o de cinco anos do Concurso organizado pela ENAP

Este livro apresenta uma an??lise de informa????es consolidadas e sistematizadas sobre experi??ncias vencedoras do Concurso de Inova????es na Gest??o P??blica Federal entre os anos de 1996 a 2000, para oferecer algumas sugest??es para o debate aos que est??o engajados no processo de aperfei??oamento do setor p??blico ou se interessam pelo tema

Colet??nea planejamento e or??amento governamental: volume 2: reforma gerencial dos processos de planejamento e or??amento

O prop??sito deste trabalho ?? apresentar os fundamentos, as dificuldades e os desafios de duas reestrutura????es: a da natureza da despesa, vigente desde 1990 e a da classifica????o funcional-program??tica, a partir do ano 2000, ambas inseridas no escopo maior da chamada reforma gerencial dos processos de planejamento e or??amento.

Semin??rio internacional sobre educa????o a dist??ncia

Semin??rio internacional sobre educa????o a dist??ncia, realizado na ENAP em 09 de julho de 2007.

Caracterização morfodinâmica do litoral Norte Fluminense, RJ, Brasil

A morfologia das praias é reflexo da ação hidrodinfunica com o tipo de sedimento disponível. Esta interação se dá a partir da base de ação das ondas, tornando-se mais efetiva na zona de arrebentação onde a energia é dissipada, originando na zona de surfe que se estende em direção à linha de costa até Q limite do espraiamento na face de praia. O objetivo deste trabalho é apresentar as características morfodinâmicas do litoral NE do Estado do Rio de Janeiro, baseando-se na análise da variação espaço Hemporal de perfis de praia, na morfologia da plataforma continental e em fotos aéreas. Os perfis de praia foram levantados durante cinco anos pela PETROBRÁS S/A, em oito estações ao longo do litoral. O comportamento morfodinâmico distinto de setores do litoral estudado é função do padrão de ondas e direção da linha de costa, morfologia da plataforma interna controlando a extensão da zona de surfe e o número de zonas de arrebentação, características dos sedimentos e gradiente da face de praia. Foram reconhecidos quatro compartimentos morfodinâmicos distintos, com base na morfologia dos perfis praiais e plataforma interna, no índice de mobilidade da linha de praia, granulometria x gradiente da face de praia e o parâmetro Q. Os quatro compartimentos são: Atafona/Foz do Rio Paraíba do Sul - estágio morfodinâmico intermediário a dissipativo, com o maior índice de mobilidade e elevada taxa de erosão; Sul de Atafona-Cabo de São Tomé - estágio intermediário a refletivo com baixo índice de mobilidade da praia; Cabo de São Tomé - estágio refletivo a intermediário com alto índice de mobilidade da praia; Cabo de São Tomé-Cabiúnas - estágio refletivo com baixo índice de mobilidade.

Mid- to Late-Holocene estuarine infilling processes studied by radiocarbon dates, high resolution seismic and biofacies at Vitoria Bay, Espirito Santo, Southeastern Brazil

A Baía de Vitória é um estuário com 20 km de comprimento, morfologicamente estreito, com um regime de micromaré e, como outros estuários modernos, formado durante a última transgressão pós-glacial. A morfologia de fundo do estrato estuarino é caracterizada por um canal natural principal limitado por planícies de maré com manguezais desenvolvidos. Datações de radiocarbono originais foram obtidas para a área. Cinco idades de radiocarbono estendendo-se de 1.010 a 7.240 anos AP foram obtidas através de dois testemunhos de sedimento, representando uma sequência estratigráfica de 5 m de espessura. Os resultados indicam que até aproximadamente 4.000 anos cal. AP, as condições ambientais da Baía de Vitória eram ainda de uma baía aberta, com uma conexão livre e aberta com águas marinhas. Durante os últimos 4.000 anos a baía experimentou uma fase de regressão importante, tornando-se mais restrita em termos de circulação da água do mar e provavelmente aumentando a energia de marés. Três superfícies estratigráficas principais foram reconhecidas, limitando fácies transgressiva, transgressiva/nível de mar alto e regressiva. A morfologia do canal atual representa um diastema de maré, mostrando fácies regressivas truncadas e erodidas. Biofácies de foraminíferos, passando de ambiente marinho para ambiente salobro e de manguezais em planície de maré confirmam a interpretação sismoestratigráfica. A ausência de biofácies de mangue em um dos dois testemunhos é tambémuma indicação de ravinamento de maré atual.

Basis set convergence in Hartree-Fock calculations of some diatomic molecules containing first and second-row atoms

Investiga-se a convergência de conjuntos de bases em direção ao limite numérico da energia Hartree-Fock (HF) total para as seqüências hierárquicas dos conjuntos de bases XZP e ccpVXZ. Para as duas hierarquias, melhoramentos significativos são obtidos com cada incremento em X. Para estimar o limite do conjunto de base completo, uma forma exponencial foi usada. Entre as várias aproximações consideradas aqui, uma extrapolação exponencial de três parâmetros aplicada aos resultados TZP, QZP e 5ZP deu os limites do conjunto de bases mais precisos. Em adição, energias HF dos orbitais moleculares ocupados mais altos de algumas moléculas diatômicas foram calculadas com o conjunto 5ZP e comparadas com as correspondentes obtidas com o conjunto cc-pV5Z e com um método numérico HF.