Molecular simulation of gas adsorption equilibria in nanoporous materials

This work is divided into two distinct parts. The first part consists of the study of the metal organic framework UiO-66Zr, where the aim was to determine the force field that best describes the adsorption equilibrium properties of two different gases, methane and carbon dioxide. The other part of the work focuses on the study of the single wall carbon nanotube topology for ethane adsorption; the aim was to simplify as much as possible the solid-fluid force field model to increase the computational efficiency of the Monte Carlo simulations. The choice of both adsorbents relies on their potential use in adsorption processes, such as the capture and storage of carbon dioxide, natural gas storage, separation of components of biogas, and olefin/paraffin separations. The adsorption studies on the two porous materials were performed by molecular simulation using the grand canonical Monte Carlo (μ,V,T) method, over the temperature range of 298-343 K and pressure range 0.06-70 bar. The calibration curves of pressure and density as a function of chemical potential and temperature for the three adsorbates under study, were obtained Monte Carlo simulation in the canonical ensemble (N,V,T); polynomial fit and interpolation of the obtained data allowed to determine the pressure and gas density at any chemical potential. The adsorption equilibria of methane and carbon dioxide in UiO-66Zr were simulated and compared with the experimental data obtained by Jasmina H. Cavka et al. The results show that the best force field for both gases is a chargeless united-atom force field based on the TraPPE model. Using this validated force field it was possible to estimate the isosteric heats of adsorption and the Henry constants. In the Grand-Canonical Monte Carlo simulations of carbon nanotubes, we conclude that the fastest type of run is obtained with a force field that approximates the nanotube as a smooth cylinder; this approximation gives execution times that are 1.6 times faster than the typical atomistic runs.

Use of Poisson spatiotemporal regression models for the Brazilian Amazon Forest: malaria count data

INTRODUCTION: Malaria is a serious problem in the Brazilian Amazon region, and the detection of possible risk factors could be of great interest for public health authorities. The objective of this article was to investigate the association between environmental variables and the yearly registers of malaria in the Amazon region using Bayesian spatiotemporal methods. METHODS: We used Poisson spatiotemporal regression models to analyze the Brazilian Amazon forest malaria count for the period from 1999 to 2008. In this study, we included some covariates that could be important in the yearly prediction of malaria, such as deforestation rate. We obtained the inferences using a Bayesian approach and Markov Chain Monte Carlo (MCMC) methods to simulate samples for the joint posterior distribution of interest. The discrimination of different models was also discussed. RESULTS: The model proposed here suggests that deforestation rate, the number of inhabitants per km², and the human development index (HDI) are important in the prediction of malaria cases. CONCLUSIONS: It is possible to conclude that human development, population growth, deforestation, and their associated ecological alterations are conducive to increasing malaria risk. We conclude that the use of Poisson regression models that capture the spatial and temporal effects under the Bayesian paradigm is a good strategy for modeling malaria counts.

Pharmacokinetic/pharmacodynamic target attainment of intravenous β-lactam regimens against Gram-negative bacteria isolated in a Brazilian teaching hospital

ABSTRACTINTRODUCTION: Monte Carlo simulations have been used for selecting optimal antibiotic regimens for treatment of bacterial infections. The aim of this study was to assess the pharmacokinetic and pharmacodynamic target attainment of intravenous β-lactam regimens commonly used to treat bloodstream infections (BSIs) caused by Gram-negative rod-shaped organisms in a Brazilian teaching hospital.METHODS: In total, 5,000 patients were included in the Monte Carlo simulations of distinct antimicrobial regimens to estimate the likelihood of achieving free drug concentrations above the minimum inhibitory concentration (MIC; fT > MIC) for the requisite periods to clear distinct target organisms. Microbiological data were obtained from blood culture isolates harvested in our hospital from 2008 to 2010.RESULTS: In total, 614 bacterial isolates, including Escherichia coli, Enterobacterspp., Klebsiella pneumoniae, Acinetobacter baumannii, and Pseudomonas aeruginosa, were analyzed Piperacillin/tazobactam failed to achieve a cumulative fraction of response (CFR) > 90% for any of the isolates. While standard dosing (short infusion) of β-lactams achieved target attainment for BSIs caused by E. coliand Enterobacterspp., pharmacodynamic target attainment against K. pneumoniaeisolates was only achieved with ceftazidime and meropenem (prolonged infusion). Lastly, only prolonged infusion of high-dose meropenem approached an ideal CFR against P. aeruginosa; however, no antimicrobial regimen achieved an ideal CFR against A. baumannii.CONCLUSIONS:These data reinforce the use of prolonged infusions of high-dose β-lactam antimicrobials as a reasonable strategy for the treatment of BSIs caused by multidrug resistant Gram-negative bacteria in Brazil.

Otimização de esquemas terapêuticos do carcinoma da próstata envolvendo braquiterapia de baixa taxa de dose e radioterapia externa : abordagens física e radiobiológica

RESUMO: Este trabalho teve como objetivo a determinação de esquemas de tratamento alternativos para o carcinoma da próstata com radioterapia externa (EBRT) e braquiterapia de baixa taxa de dose (LDRBT) com implantes permanentes de Iodo-125, biologicamente equivalentes aos convencionalmente usados na prática clínica, com recurso a modelos teóricos e a métodos de Monte Carlo (MC). Os conceitos de dose biológica efetiva (BED) e de dose uniforme equivalente (EUD) foram utilizados, com o modelo linear-quadrático (LQ), para a determinação de regimes de tratamento equivalentes. Numa primeira abordagem, utilizou-se a BED para determinar: 1) esquemas hipofracionados de EBRT mantendo as complicações retais tardias de regimes convencionais com doses totais de 75,6 Gy, 77,4 Gy, 79,2 Gy e 81,0 Gy; e 2) a relação entre as doses totais de EBRT e LDRBT de modo a manter a BED do regime convencional de 45 Gy de EBRT e 110 Gy de LDRBT. Numa segunda abordagem, recorreu-se ao código de MC MCNPX para a simulação de distribuições de dose de EBRT e LDRBT em dois fantomas de voxel segmentados a partir das imagens de tomografia computorizada de pacientes com carcinoma da próstata. Os resultados das simulações de EBRT e LDRBT foram somados e determinada uma EUD total de forma a obterem-se: 1) esquemas equivalentes ao tratamento convencional de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT; e 2) esquemas equivalentes a EUD na próstata de 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy e 110 Gy. Em todos os resultados nota-se um ganho terapêutico teórico na utilização de esquemas hipofracionados de EBRT. Para uma BED no reto equivalente ao esquema convencional, tem-se um aumento de 2% na BED da próstata com menos 5 frações. Este incremento dá-se de forma cada vez mais visível à medida que se reduz o número de frações, sendo da ordem dos 10-11% com menos 20 frações e dos 35-45% com menos 40 frações. Considerando os resultados das simulações de EBRT, obteve-se uma EUD média de 107 Gy para a próstata e de 42 Gy para o reto, com o esquema convencional de 110 Gy de LDRBT, seguidos de 25 frações de 1,8 Gy de EBRT. Em termos de probabilidade de controlo tumoral (igual EUD), é equivalente a este tratamento a administração de EBRT em 66 frações de 1,8 Gy, 56 de 2 Gy, 40 de 2,5 Gy, 31 de 3 Gy, 20 de 4 Gy ou 13 de 5 Gy. Relativamente à administração de 66 frações de 1,8 Gy, a EUD generalizada no reto reduz em 6% com o recurso a frações de 2,5 Gy e em 10% com frações de 4 Gy. Determinou-se uma BED total de 162 Gy para a administração de 25 frações de 1,8 Gy de EBRT em combinação com 110 Gy de LDRBT. Variando-se a dose total de LDRBT (TDLDRBT) em função da dose total de EBRT (TDEBRT), de modo a garantir uma BED de 162 Gy, obteve-se a seguinte relação:.......... Os resultados das simulações mostram que a EUD no reto diminui com o aumento da dose total de LDRBT para dose por fração de EBRT (dEBRT) inferiores a 2, Gy e aumenta para dEBRT a partir dos 3 Gy. Para quantidades de TDLDRBT mais baixas (<50 Gy), o reto beneficia de frações maiores de EBRT. À medida que se aumenta a TDLDRBT, a EUD generalizada no reto torna-se menos dependente da dEBRT. Este trabalho mostra que é possível a utilização de diferentes regimes de tratamento para o carcinoma da próstata com radioterapia que possibilitem um ganho terapêutico, quer seja administrando uma maior dose biológica com efeitos tardios constantes, quer mantendo a dose no tumor e diminuindo a toxicidade retal. A utilização com precaução de esquemas hipofracionados de EBRT, para além do benefício terapêutico, pode trazer vantagens ao nível da conveniência para o paciente e economia de custos. Os resultados das simulações deste estudo e conversão para doses de efeito biológico para o tratamento do carcinoma da próstata apresentam linhas de orientação teórica de interesse para novos ensaios clínicos. --------------------------------------------------ABSTRACT: The purpose of this work was to determine alternative radiotherapy regimens for the treatment of prostate cancer using external beam radiotherapy (EBRT) and low dose-rate brachytherapy (LDRBT) with Iodine-125 permanent implants which are biologically equivalent to conventional clinical treatments, by the use of theoretical models and Monte Carlo techniques. The concepts of biological effective dose (BED) and equivalent uniform dose (EUD), together with the linear-quadratic model (LQ), were used for determining equivalent treatment regimens. In a first approach, the BED concept was used to determine: 1) hypofractionated schemes of EBRT maintaining late rectal complications as with the conventional regimens with total doses of 75.6 Gy, 77.4 Gy, 79.2 Gy and 81.0 Gy; and 2) the relationship between total doses of EBRT and LDRBT in order to keep the BED of the conventional treatment of 45 Gy of EBRT and 110 Gy of LDRBT. In a second approach, the MC code MCNPX was used for simulating dose distributions of EBRT and LDRBT in two voxel phantoms segmented from the computed tomography of patients with prostate cancer. The results of the simulations of EBRT and LDRBT were added up and given an overall EUD in order to obtain: 1) equivalent to conventional treatment regimens of 25 fraction of 1.8 Gy of EBRT in combination with 110Gy of LDRBT; and 2) equivalent schemes of EUD of 67 Gy, 72 Gy, 80 Gy, 90 Gy, 100 Gy, and 110Gy to the prostate. In all the results it is noted a therapeutic gain using hypofractionated EBRT schemes. For a rectal BED equivalent to the conventional regimen, an increment of 2% in the prostate BED was achieved with less 5 fractions. This increase is visibly higher as the number of fractions decrease, amounting 10-11% with less 20 fractions and 35-45% with less 20 fractions. Considering the results of the EBRT simulations an average EUD of 107 Gy was achieved for the prostate and of 42 Gy for the rectum with the conventional scheme of 110 Gy of LDRBT followed by 25 fractions of 1.8 Gy of EBRT. In terms of tumor control probability (same EUD) it is equivalent to this treatment, for example, delivering the EBRT in 66 fractions of 1.8 Gy, 56 fractions of 2 Gy, 40 fractions of 2.5 Gy, 31 fractions of 3 Gy, 20 fractions of 4 Gy or 13 fractions of 5 Gy. Regarding the use of 66 fractions of 1.8 Gy, the rectum EUD is reduced to 6% with 2.5 Gy per fraction and to 10% with 4 Gy. A total BED of 162 Gy was achieved for the delivery of 25 fractions of 1.8 Gy of EBRT in combination with 110 Gy of LDRBT. By varying the total dose of LDRBT (TDLDRBT) with the total dose of EBRT (TDEBRT) so as to ensure a BED of 162 Gy, the following relationship was obtained: ....... The simulation results show that the rectum EUD decreases with the increase of the TDLDRBT, for EBRT dose per fracion (dEBRT) less than 2.5 Gy and increases for dEBRT above 3 Gy. For lower amounts of TDLDRBT (< 50Gy), the rectum benefits of larger EBRT fractions. As the TDLDRBT increases, the rectum gEUD becomes less dependent on the dEBRT. The use of different regimens which enable a therapeutic gain, whether deivering a higher dose with the same late biological effects or maintaining the dose to the tumor and reducing rectal toxicity is possible. The use with precaution of hypofractionated regimens, in addition to the therapeutic benefit, can bring advantages in terms of convenience for the patient and cost savings. The simulation results of this study together with the biological dose conversion for the treatment of prostate cancer serve as guidelines of interest for new clinical trials.

Structure-property relationships in polymer nanocomposites

Tese de Doutoramento em Ciência e Engenharia de Polímeros e Compósitos

Measurement of the tt¯ production cross-section as a function of jet multiplicity and jet transverse momentum in 7 TeV proton--proton collisions with the ATLAS detector

The tt¯ production cross-section dependence on jet multiplicity and jet transverse momentum is reported for proton--proton collisions at a centre-of-mass energy of 7 TeV in the single-lepton channel. The data were collected with the ATLAS detector at the CERN Large Hadron Collider and comprise the full 2011 data sample corresponding to an integrated luminosity of 4.6 fb−1. Differential cross-sections are presented as a function of the jet multiplicity for up to eight jets using jet transverse momentum thresholds of 25, 40, 60, and 80 GeV, and as a function of jet transverse momentum up to the fifth jet. The results are shown after background subtraction and corrections for all detector effects, within a kinematic range closely matched to the experimental acceptance. Several QCD-based Monte Carlo models are compared with the results. Sensitivity to the parton shower modelling is found at the higher jet multiplicities, at high transverse momentum of the leading jet and in the transverse momentum spectrum of the fifth leading jet. The MC@NLO+HERWIG MC is found to predict too few events at higher jet multiplicities.

On the Hill estimator: a comparison of methods

Extreme value theory (EVT) deals with the occurrence of extreme phenomena. The tail index is a very important parameter appearing in the estimation of the probability of rare events. Under a semiparametric framework, inference requires the choice of a number k of upper order statistics to be considered. This is the crux of the matter and there is no definite formula to do it, since a small k leads to high variance and large values of k tend to increase the bias. Several methodologies have emerged in literature, specially concerning the most popular Hill estimator (Hill, 1975). In this work we compare through simulation well-known procedures presented in Drees and Kaufmann (1998), Matthys and Beirlant (2000), Beirlant et al. (2002) and de Sousa and Michailidis (2004), with a heuristic scheme considered in Frahm et al. (2005) within the estimation of a different tail measure but with a similar context. We will see that the new method may be an interesting alternative.

Extreme value Birnbaum-Saunders regression models applied to environmental data

Extreme value models are widely used in different areas. The Birnbaum–Saunders distribution is receiving considerable attention due to its physical arguments and its good properties. We propose a methodology based on extreme value Birnbaum–Saunders regression models, which includes model formulation, estimation, inference and checking. We further conduct a simulation study for evaluating its performance. A statistical analysis with real-world extreme value environmental data using the methodology is provided as illustration.

Strain induced edge magnetism at zigzag edge of a graphene quantum dot

We study the temperature dependent magnetic susceptibility of a strained graphene quantum dot by using the determinant quantum Monte Carlo method. Within the Hubbard model on a honeycomb lattice, our unbiased numerical results show that a relative small interaction $U$ may lead to a edge ferromagnetic like behavior in the strained graphene quantum dot, and a possible room temperature transition is suggested. Around half filling, the ferromagnetic fluctuations at the zigzag edge is strengthened both markedly by the on-site Coulomb interaction and the strain, especially in low temperature region. The resultant strongly enhanced ferromagnetic like behavior may be important for the development of many applications.

Estimativa de área de vegetação secundária na Amazônia Legal Brasileira

A vegetação secundária tem funções relevantes para os ecossistemas, tais como a fixação de carbono atmosférico, a manutenção da biodiversidade, o estabelecimento da conectividade entre remanescentes florestais, manutenção dos regime hidrológico e a recuperação da fertilidade do solo. O objetivo deste trabalho é, através de uma abordagem amostral, estimar a área ocupada por vegetação secundária na Amazônia Legal Brasileira (AML) em 2006. A amostragem se baseia em uma abordagem estratificada pelo grau de desflorestamento das cenas LANDSAT-TM que recobrem a AML. Foram selecionadas 26 cenas para o ano de 2006, distribuídas em sete estratos conforme o percentual de desflorestamento, nas quais foram mapeadas as áreas de vegetação secundária a partir de técnicas de classificação de imagens. Foi desenvolvido um modelo multivariado de regressão para estimar a área de vegetação secundária utilizando como variáveis independentes a área de desflorestamento, a área de hidrografia, a estrutura agrária, e área das unidades de conservação. A análise de regressão encontrou um R2 ajustado de 0,84 , e coeficientes positivos para a proporção de hidrografia na imagem (2,055) e para a estrutura agrária (0,197), e coeficientes negativos para o grau de desflorestamento na imagem (-0,232) e para a proporção de Unidades de Conservação na imagem (-0,262). O modelo de regressão estimou uma área de 131.873 km² de vegetação secundária para o ano de 2006. Aplicando uma simulação Monte Carlo foi estimada uma incerteza de aproximadamente 12.445 km² para a área.

Development of a methodology to incorporate risk and uncertainty in electricity power planning

Dissertação de mestrado em Engenharia Industrial

Conformal symmetry of the critical 3D Ising model inside a sphere

We perform Monte-Carlo simulations of the three-dimensional Ising model at the critical temperature and zero magnetic field. We simulate the system in a ball with free boundary conditions on the two dimensional spherical boundary. Our results for one and two point functions in this geometry are consistent with the predictions from the conjectured conformal symmetry of the critical Ising model.

Programa com abordagem probabilística do sistema de avaliação de qualidade de taludes Slope Quality Index (SQI)

Dissertação de mestrado integrado em Engenharia Civil

Procesos físicos en relatividad general

Para 1997 se propone realizar las siguientes tareas: Teoría de Superficies Nulas en Relatividad General. * Estudiar la formación de cáusticas y singularidades de las superficies nulas en forma cinemática y dinámica. * Analizar la dinámica de superficies características hasta segundo orden en un esquema perturbativo. * Formalizar la teoría de Superficies Nulas usando espacios fibrados. Cuantificación asintótica del campo gravitatorio. Continuar el estudio del espacio de Hilbert asintótico para el campo gravitatorio. Es conveniente recordar que la cuantifiación de este campo es uno de los problemas teóricos más importantes sin resolver. Cuantificación de Superficies Nulas. Continuar con el estudio de la cuantización del campo gravitatorio usando el formalismo de superficies nulas. Hasta el momento el principal resultado obtenido fue probar que no sólo el campo gravitatorio sino también los puntos del espacio-tiempo se vuelven operadores con relaciones de conmutación. Esto demuestra de manera concreta la especulación de Wheeler que en la gravedad cuántica los puntos del espacio tiempo no están bien definidos. Simulación de un detector de rayos Gamma. (...) Para este año se prevé cumplir las siguientes etapas: * Se estudiarán los efectos de la polarización del fotón incidente en la sección eficaz y su posible método de detección. * Se implementarán las subrutinas correspondientes en las simulaciones Monte Carlo para incorporar secciones eficaces polarizadas en el Código GEANT. Asimismo, se piensa comenzar con el diseño de un nuevo detector de rayos Gamma en un rango de energías (0.1 / 1 MeV) muy inferior al anterior para ser usado como cámara Gamma ultra sensible.

Espectrometría de emisiones en la interacción de radiaciones ionizantes con la materia

En el área de los estudios básicos de la física de radiaciones, el Laboratorio de Espectrometría de Radiaciones (L.E.R.) de la U.N.C. desarrolla líneas de investigación cuyos resultados son de uso directo en física de radiaciones aplicadas. (...) Objetivos generales y específicos: * Estudios de Detectores de Radiaciones Ionizantes: se continúan los estudios en Argón, teniendo como objetivo específico completar y comprobar técnicas de simulación con el método Monte Carlo para detectores con simetría cilíndrica y gradiente de campo eléctrico. * Análisis por fluorescencia de Rayos X: los objetivos específicos están dedicados a la evaluación de diseños para aumentar la eficiencia de polarización de Rayos X. En cuanto al análisis de elementos, se considerarán las posibles mejoras experimentales mediante el aumento de la relación señal-ruido en la fluorescencia con Rayos X polarizados, disminuyendo los umbrales de detección de concentración de elementos, en muestras compuestas por muchos elementos. * Dosimetría de radiaciones ionizantes: En cuanto a la dosimetría mediante LiF se planea seguir estudiando la difusión de los positrones con bajos contenidos de centros de color y la linealidad de respuesta de la vida media en función de la dosis de radiación para diferentes clases de cristales termoluminiscentes. Para la dosimetría con centelladores plásticos, se pretende en el siguiente paso, la puesta a punto final del dosímetro y electrómetro desarrollado como así también su calibración para su posterior prueba en centros de radioterapia. * Análisis de composición de elementos químicos en el cuerpo humano: En un futuro mediato, se prevé la culminación de los estudios de una nueva técnica para determinación de concentración de elementos minoritarios y mayoritarios en órganos y en extremidades por dispersión incoherente de Rayos X. Esta técnica contempla la ulterior determinación in vivo de la composición del cuerpo humano.