Redução de incertezas em radioterapia utilizando simulação Monte Carlo: análise espectral aplicada à correção de dose absorvida

OBJETIVO: Determinar, por simulação Monte Carlo, os espectros de feixes de cobaltoterapia em profundidade na água e fatores de correção para doses absorvidas em dosímetros termoluminescentes de fluoreto de lítio. MATERIAIS E MÉTODOS: As simulações dos espectros secundários da fonte clínica de cobalto-60 foram realizadas com o código Monte Carlo PENELOPE, em diversas profundidades na água. Medidas experimentais de dose profunda foram obtidas com dosímetros termoluminescentes e câmara de ionização em condições de referência em radioterapia. Os fatores de correção para os dosímetros termoluminescentes foram obtidos através da razão entre as absorções relativas ao espectro de baixa energia e ao espectro total. RESULTADOS: A análise espectral em profundidade revelou a existência de espectros secundários de baixa energia responsáveis por uma parcela significativa da deposição de dose. Foram observadas discrepâncias de 3,2% nas doses medidas experimentalmente com a câmara de ionização e com os dosímetros termoluminescentes. O uso dos fatores de correção nessas medidas permitiu diminuir a discrepância entre as doses absorvidas para, no máximo, 0,3%. CONCLUSÃO: Os espectros simulados permitem o cálculo de fatores de correção para as leituras de dosímetros termoluminescentes utilizados em medidas de dose profunda, contribuindo para a redução das incertezas associadas ao controle de qualidade de feixes clínicos em radioterapia.

Simulação e análise dosimétrica de protonterapia e íons de carbono no tratamento do melanoma uveal

OBJETIVO: Este artigo apresenta a avaliação dosimétrica da radioterapia por íons de carbono em comparação à protonterapia. MATERIAIS E MÉTODOS: As simulações computacionais foram elaboradas no código Geant4 (GEometry ANd Tracking). Um modelo de olho discretizado em voxels implementado no sistema Siscodes (sistema computacional para dosimetria em radioterapia) foi empregado, em que perfis de dose em profundidade e curvas de isodose foram gerados e superpostos. Nas simulações com feixe de íons de carbono, distintos valores de energia do feixe foram adotados, enquanto nas simulações com feixe de prótons os dispositivos da linha de irradiação foram incluídos e diferentes espessuras do material absorvedor foram aplicadas. RESULTADOS: As saídas das simulações foram processadas e integradas ao Siscodes para gerar as distribuições espaciais de dose no modelo ocular, considerando alterações do posicionamento de entrada do feixe. Os percentuais de dose foram normalizados em função da dose máxima para um feixe em posição de entrada específica, energia da partícula incidente e número de íons de carbono e de prótons incidentes. CONCLUSÃO: Os benefícios descritos e os resultados apresentados contribuem para o desenvolvimento das aplicações clínicas e das pesquisas em radioterapia ocular por íons de carbono e prótons.

Simulação da Química da Atmosfera Poluída por Automóveis Movidos a Álcool

Numerical modeling studies of urban atmospheres were performed with a photochemical box submodel which incorporates primary emission rates, temperature and solar irradiance rate coefficients. Actinic flux calculations were made for an appropriate latitude for Rio de Janeiro and initial concentrations and emission rates were taken from experimental data. The relative importance of individual reactions is discussed as well as the effect of enhanced aldehyde emissions from vehicles using ethanol and ethanol-containing fuel.

Simulação do processo Solvay no laboratório didático

The Solvay Process is an important route for the industrial preparation of sodium carbonate from brine and CO2. However, experiments illustrating such process are performed only in a few chemistry courses, mainly because of the difficulties to precipitate out the intermediate product NaHCO3 within the 3 or 4 hours generally available for laboratory classes. In this article we describe a method that properly simulates the Solvay process using glass apparatus that can be found in any chemistry laboratory. The good yields and purity attainable, associated with the short time (~3 hours) required to get the sodium bicarbonate, account for the suitability of our method for chemistry or chemical engineering undergraduate students laboratory classes.

Termodinâmica estatística de líquidos com o método de Monte Carlo. I. Metodologia

Statistical mechanics Monte Carlo simulation is reviewed as a formalism to study thermodynamic properties of liquids. Considering the importance of free energy changes in chemical processes, the thermodynamic perturbation theory implemented in the Monte Carlo method is discussed. The representation of molecular interaction by the Lennard-Jones and Coulomb potential functions is also discussed. Charges derived from quantum molecular electrostatic potential are also discussed as an useful methodology to generate an adequate set of partial charges to be used in liquid simulation.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Otimização de um mecanismo fotoquímico para a simulação da atmosfera urbana brasileira

A photochemical reaction mechanism for the description of air quality in Brazilian urban regions is described and evaluated by comparison with chamber experiments. The mechanism was developed for use in ozone modeling and application of control strategies. The oxidation of ethanol and methyl-ter-butyl-ether is also considered. Using this chemical model, a trajectory simulation of Brazil Avenue, Rio de Janeiro, was performed. The model predicts that ozone should reach a maximum of 22.4 ppb at 14:57 h. This value is in good agreement with the experimental measurements of 22.5 ppb for 14:00 h and 22.3 ppb for 15:00 h.

Estudo da geometria da uréia por métodos ab initio e simulação computacional de líquidos

A study was carried out on the urea geometries using ab initio calculation and Monte Carlo computational simulation of liquids. The ab initio calculated results showed that urea has a non-planar conformation in the gas phase in which the hydrogen atoms are out of the plane formed by the heavy atoms. Free energies associated to the rotation of the amino groups of urea in water were obtained using the Monte Carlo method in which the thermodynamic perturbation theory is implemented. The magnitude of the free energy obtained from this simulation did not permit us to conclude that urea is non-planar in water.

Comparação de duas metodologias de amostragem atmosférica com ferramenta estatística não paramétrica

In atmospheric aerosol sampling, it is inevitable that the air that carries particles is in motion, as a result of both externally driven wind and the sucking action of the sampler itself. High or low air flow sampling speeds may lead to significant particle size bias. The objective of this work is the validation of measurements enabling the comparison of species concentration from both air flow sampling techniques. The presence of several outliers and increase of residuals with concentration becomes obvious, requiring non-parametric methods, recommended for the handling of data which may not be normally distributed. This way, conversion factors are obtained for each of the various species under study using Kendall regression.

Simulação de desdobramentos complexos de espectros de RMN de ¹H

Complex ¹H NMR spectra multiplets that cannot be easily understood by simple inspection are rather frequent in the daily work of the organic chemistry analyst. The multiple and excellent new techniques available from modern instruments usually provide satisfactory solutions, but there are still many cases where a simulation is necessary, at least to obtain a final confirmation. It is extremely convincing to see that a graph, obtained by calculations with chemical-shift and coupling-constant values only, can be virtually identical to the experimental spectrum. This paper describes a computer program to make such calculations. The program is free and can be downloaded from http://artemis.ffclrp.usp.br/NMR.htm (click on SimEsp_NMR_Compil.zip). All routines are also available and may be used without any restrictions. The paper includes a fairly detailed discussion about how the calculations are made.

Análise estatística e optimização de perfis de redução termoprogramada (TPR)

The effect of operational variables and their interaction in TPR profiles was studied using a fractional factorial experimental design. The heating rate and the reducing agent concentration were found to be the most important variables determining the resolution and sensitivity of the technique. They showed opposite effects. Therefore, they should be manipulated preferentially in order to obtain optimized TPR profiles. The effect of sample particle size was also investigated. The tests were carried out within a Cu/Zn/Al2O3 catalyst used for the water-gas shift reaction that presented two distinct species of Cu2+ in TPR profiles.

TitGer: uma planilha eletrônica para simulação de titulação de mistura de compostos polipróticos

The aim of this paper is to present a simple way of treating the general equation for acid-base titrations based on the concept of degree of dissociation, and to propose a new spreadsheet approach for simulating the titration of mixtures of polyprotic compounds. The general expression, without any approximation, is calculated a simple iteration method, making number manipulation easy and painless. The user-friendly spreadsheet was developed by using MS-Excel and Visual-Basic-for-Excel. Several graphs are drawn for helping visualizing the titration behavior. A Monte Carlo function for error simulation was also implemented. Two examples for titration of alkalinity and McIlvaine buffer are presented.

Avaliação de contaminação ambiental causada por poluentes orgânicos persistentes utilizando simulação computacional

Environmental concern is growing in the current days and there is global agreement to banish production and use of persistent organic pollutants (POP). The synthetic insecticides chlordecone and mirex, classified as POPs, have similar structures and they are potentially toxic. This work uses properties and physicochemical constants related to the pesticides and computational simulation to evaluate the leach phenomenon and persistency in soil. The largest tendency of persistence of the compound is found to be in the surface of soil, but even low concentration in water represents a high risk due to bioaccumulation in adipose tissue.

Hidrocarbonetos policíclicos aromáticos como traçadores da queima de cana-de-açúcar: uma abordagem estatística

In this study atmospheric particulates of PAHs were measured in Araraquara, Piracicaba and São Paulo in July 2003 (sugarcane harvest season in Araraquara and Piracicaba) and in Araraquara in March of 2003. The results were normalized to the total PAH concentrations. Comparison among the sites and principal component analysis (PCA) were used to investigate possible tracers of emission. Fluoranthene and pyrene concentrations were higher in Piracicaba and Araraquara samples. These PAH were also responsible for the largest negative loadings on the second principal component and account for the negative scores and for the formation of the Araraquara and Piracicaba group.

Simulação de espectros de ressonância paramagnética eletrônica (RPE) através do programa NLSL

EPR users often face the problem of extracting information from frequently low-resolution and complex EPR spectra. Simulation programs that provide a series of parameters, characteristic of the investigated system, have been used to achieve this goal. This work describes the general aspects of one of those programs, the NLSL program, used to fit EPR spectra applying a nonlinear least squares method. Several motion regimes of the probes are included in this computational tool, covering a broad range of spectral changes. The meanings of the different parameters and rotational diffusion models are discussed. The anisotropic case is also treated by including an orienting potential and order parameters. Some examples are presented in order to show its applicability in different systems.