Estudo teórico da contribuição do espalhamento nuclear não elástico na trajetória de prótons de energias entre 50 e 250 MeV em materiais de interesse em radiologia: aplicação em tomografia com prótons

The goal of this work is to study the process of interaction of protons with matter through Monte Carlo simulation. For this purpose, it was employed the SRIM program (Stopping and Range of Ions in Matter ) and MCNPX (Monte Carlo N-Particle eXtended) v2.50. This work is going to support the development of a tomography system with protons. It was studied the interaction of proton with the follow materials: Polimethyl Mehacralate (PMMA), MS20 Tissue Substitute and water. This work employed energies in range of 50 MeV and 250 MeV, that is the range of clinical interest. The energy loss of proton after cross a material layer, the decreasing of its intensity, the angular and lateral de ection of incident beam, including and excluding nuclear interactions. This work is related with Medical Physics and Material Physics, like interaction of radiation with matter, particle transport phenomena, and the experimental methods in Nuclear Physics like simulation and computational by Monte Carlo method

A Monte Carlo study of the anisotropic N=3 Ashkin-Teller model

The phase diagram of an asymmetric N = 3 Ashkin-Teller model is obtained by a numerical analysis which combines Monte Carlo renormalization group and reweighting techniques. Present results reveal several differences with those obtained by mean-field calculations and a Hamiltonian approach. In particular, we found Ising critical exponents along a line where Goldschmidt has located the Kosterlitz-Thouless multicritical point. On the other hand, we did find nonuniversal exponents along another transition line. Symmetry breaking in this case is very similar to the N = 2 case, since the symmetries associated to only two of the Ising variables are broken. However, for large values of the coupling constant ratio XW = W/K, when the only broken symmetry is of a hidden variable, we detected first-order phase transitions giving evidence supporting the existence of a multicritical point, as suggested by Goldschmidt, but in a different region of the phase diagram. © 2002 Elsevier Science B.V. All rights reserved.

Eficiência energética de sistemas de produção de mamona no Estado de Minas Gerais

The objective of this study is to determine the economical and energetic efficiency of the castor bean culture in "Zona da Mata" and South of Minas Gerais. Through the method of simulation "Monte Carlo", we verified the probabilities of occurrence of the economical, cultural and energetic efficiency indexes. In relation to the production systems of castor bean in Minas Gerais in the season 2005/2006, we established that the variables price and productivity were the most noticeable for the producers from "Zona da Mata", while in the South of Minas it was productivity. We verified that the probability of the economical efficiency index to be lower than one was 43,26% for the producers from "Zona da Mata" and 39,57% for the ones from the South of Minas. The medium price received covered the medium costs of production. However, we observed that the medium costs in these regions of Minas Gerais, were over the minimum price. Regarding the energetic analysis, the results showed that the systems in the regions studied in Minas Gerasi, presented average of cultural efficiency indexes of 8,26 and 18,89. We concluded that despite the result being favorable from the energetic point of view, from the economical sustainability point of view there is the need of a more effective support policy for the castor bean, taking into consideration that the expectations of the producers with PNPB were not confirmed.

Optimization of x-ray spectra in digital mammography through Monte Carlo simulations

In this work, a Monte Carlo code was used to investigate the performance of different x-ray spectra in digital mammography, through a figure of merit (FOM), defined as FOM = CNR2/(D) over bar (g), with CNR being the contrast-to-noise ratio in image and (D) over bar (g) being the average glandular dose. The FOM was studied for breasts with different thicknesses t (2 cm <= t <= 8 cm) and glandular contents (25%, 50% and 75% glandularity). The anode/filter combinations evaluated were those traditionally employed in mammography (Mo/Mo, Mo/Rh, Rh/Rh), and a W anode combined with Al or K-edge filters (Zr, Mo, Rh, Pd, Ag, Cd, Sn), for tube potentials between 22 and 34 kVp. Results show that the W anode combined with K-edge filters provides higher values of FOM for all breast thicknesses investigated. Nevertheless, the most suitable filter and tube potential depend on the breast thickness, and for t >= 6 cm, they also depend on breast glandularity. Particularly for thick and dense breasts, a W anode combined with K-edge filters can greatly improve the digital technique, with the values of FOM up to 200% greater than that obtained with the anode/filter combinations and tube potentials traditionally employed in mammography. For breasts with t < 4 cm, a general good performance was obtained with the W anode combined with 60 mu m of the Mo filter at 24-25 kVp, while 60 mu m of the Pd filter provided a general good performance at 24-26 kVp for t = 4 cm, and at 28-30 and 29-31 kVp for t = 6 and 8 cm, respectively.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Study of surface free energy on a lattice-gas model applied to Liquid bridge

The pulmonary crackling and the formation of liquid bridges are problems that for centuries have been attracting the attention of scientists. In order to study these phenomena, it was developed a canonical cubic lattice-gas­ like model to explain the rupture of liquid bridges in lung airways [A. Alencar et al., 2006, PRE]. Here, we further develop this model and add entropy analysis to study thermodynamic properties, such as free energy and force. The simulations were performed using the Monte Carlo method with Metropolis algorithm. The exchange between gas and liquid particles were performed randomly according to the Kawasaki dynamics and weighted by the Boltzmann factor. Each particle, which can be solid (s), liquid (l) or gas (g), has 26 neighbors: 6 + 12 + 8, with distances 1, √2 and √3, respectively. The energy of a lattice's site m is calculated by the following expression: Em = ∑k=126 Ji(m)j(k) in witch (i, j) = g, l or s. Specifically, it was studied the surface free energy of the liquid bridge, trapped between two planes, when its height is changed. For that, was considered two methods. First, just the internal energy was calculated. Then was considered the entropy. It was fond no difference in the surface free energy between this two methods. We calculate the liquid bridge force between the two planes using the numerical surface free energy. This force is strong for small height, and decreases as the distance between the two planes, height, is increased. The liquid-gas system was also characterized studying the variation of internal energy and heat capacity with the temperature. For that, was performed simulation with the same proportion of liquid and gas particle, but different lattice size. The scale of the liquid-gas system was also studied, for low temperature, using different values to the interaction Jij.

Modelización econométrica a través de cadenas de Markov de Monte Carlo (MCMC) en valoración contingente

Doctorado en Análisis Económico. Programa en Análisis Económico Aplicado

Paralelización del algoritmo Monte Carlo Ray Tracing para Comunicaciones Ópticas Submarinas

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Caracterización computacional de un detector de germanio hiperpuro de rango extendido (HPGe- XtRa) con simulación Montecarlo y optimización con un algoritmo evolutivo

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Quantensimulationen physisorbierter Molekülschichten auf Graphit

Mit Hilfe der Pfadintegral-Monte Carlo-Methode werdenPhasendiagramme von physisorbierten Molekülschichten aufGraphit untersucht. Die Verwendung von realistischen Potenzialen sowie dieBehandlung aller translatorischen und rotatorischenFreiheitsgrade erlaubt einen quantitativen Vergleich mit denExperimenten.Krypton-Atome bilden in der Monolage ein kommensurablesGitter mit den Atomen über der Mitte jeder drittenGraphitwabe.Die Vorgänge am Schmelzübergang werden von der Desorptioneiniger Atome dominiert. Die Argon-Schicht auf Graphit ist dagegen inkommensurabel.Zweiatomigen Stickstoff-Moleküle bilden eineorientierungsgeordnete Tieftemperaturphase(Fischgrät-Struktur). Quantenfluktuationen führen zu einer Erniedrigung der mitklassischen Methoden berechneten Phasenübergangstemperaturum 12%.Damit wird der experimentelle Wert von 28 K erreicht.Die Anisotropie und das Dipolmoment von Kohlenmonoxid führenzu einer dipolar geordneten Tieftemperaturphase.Die experimentell nicht geklärte Struktur kann in derQuantensimulation als antiferroelektrischeFischgrät-Struktur identifiziert werden.Der Phasenübergang liegt mit 6 K sehr nahe am Experiment(5.2 K).Für die Argon-Stickstoff-Mischsysteme wird dasPhasendiagramm in der Konzentrations(x)-Temperatur(T)-Ebeneerstellt. Die Übergangstemperaturen decken sich mit denen desExperiments.In Konfigurationen mit zufälliger Teilchenbesetzung weisen die linearen Moleküle ab Argon-Konzentrationen von10% ein Orientierungsglas-Verhalten auf.Durch einen zusätzlichen Teilchenaustausch wird in denMischsystemen die Bildung einer Windrad-Phase ermöglicht, inder die Argon-Atome eine Überstruktur annehmen.Diese Phase wird experimentell imArgon-Kohlenmonoxid-Mischsystem vorgefunden, dessenx-T-Phasendiagramm in guter Übereinstimmung mit denSimulationsergebnissen steht.Die explizite Berücksichtigung der Quantenmechanik in denComputersimulationen liefert wesentliche Beiträge zurKlärung des Phasenverhaltens und der Bestimmung vonÜbergangstemperaturen der Tieftemperaturstrukturen.