Estudo de solvatocromismo em líquidos moleculares orgânicos via método seqüencial Monte Carlo/Mecânica Quântica

Utilizou-se o método seqüencial Monte Carlo / Mecânica Quântica para obterem-se os desvios de solvatocromismo e os momentos de dipolo dos sistemas de moléculas orgânicas: Uracil em meio aquoso, -Caroteno em Ácido Oléico, Ácido Ricinoléico em metanol e em Etanol e Ácido Oléico em metanol e em Etanol. As otimizações das geometrias e as distribuições de cargas foram obtidas através da Teoria do Funcional Densidade com o funcional B3LYP e os conjuntos de funções de base 6-31G(d) para todas as moléculas exceto para a água e Uracil, as quais, foram utilizadas o conjunto de funções de base 6-311++G(d,p). No tratamento clássico, Monte Carlo, aplicou-se o algoritmo Metropólis através do programa DICE. A separação de configurações estatisticamente relevantes para os cálculos das propriedades médias foi implementada com a utilização da função de auto-correlação calculada para cada sistema. A função de distribuição radial dos líquidos moleculares foi utilizada para a separação da primeira camada de solvatação, a qual, estabelece a principal interação entre soluto-solvente. As configurações relevantes da primeira camada de solvatação de cada sistema foram submetidas a cálculos quânticos a nível semi-empírico com o método ZINDO/S-CI. Os espectros de absorção foram obtidos para os solutos em fase gasosa e para os sistemas de líquidos moleculares comentados. Os momentos de dipolo elétrico dos mesmos também foram obtidos. Todas as bandas dos espectros de absorção dos sistemas tiveram um desvio para o azul, exceto a segunda banda do sistema de Beta-Caroteno em Ácido Oléico que apresentou um desvio para o vermelho. Os resultados encontrados apresentam-se em excelente concordância com os valores experimentais encontrados na literatura. Todos os sistemas tiveram aumento no momento de dipolo elétrico devido às moléculas dos solventes serem moléculas polares. Os sistemas de ácidos graxos em álcoois apresentaram resultados muito semelhantes, ou seja, os ácidos graxos mencionados possuem comportamentos espectroscópicos semelhantes submetidos aos mesmos solventes. As simulações através do método seqüencial Monte Carlo / Mecânica Quântica estudadas demonstraram que a metodologia é eficaz para a obtenção das propriedades espectroscópicas dos líquidos moleculares analisados.

Investigação da dinâmica sequêncial Monte Carlo/Mecânica Quântica para sitemas moleculares orgânicos

O principal objetivo deste trabalho é investigar teóricamente as propriedades eletrônicas e óticas de compostos orgânicos em ambiente líquido. A compreensão das interações em sistemas líquidos é muito importante para a descrição de fenômenos em muitas partes da ciência, como Física, Química, Biologia e Materiais, bem como no desenvolvimento de materiais óticos. As flutuações térmicas fazem que inumaram das configurações para os líquidos possam ser gerados. Esta é a razão do comportamento estatístico observado em sistemas líquidos. Considerando estas dificuldades, o tratamento sequencial Monte Carlo / Mecânica Quântica (SMC/QM) é usado neste trabalho. Neste procedimento, a estrutura líquida é gerada primeiramente por simulações clássicas de MC e mais tarde, somente a parte mais importante do sistema é tratada com mecânica quântica. Usando o procedimento acima, os propriedades do quantum o pirazine dimethyl do thiene Methyl da laranja (MO) e do O 2,3-dimetil tieno[3,4-b] pirazina (DTP) foram investigados. O MO é um conhecido indicador de pH e pode ser encontrado sob circunstâncias básicas e acidas. Suas geometrias de mínima energia foram obtidas mediante a Teoria do Funcional da Densidade pelo funcional B3LYP, sendo o sistema descrito pelas bases de Pople com uma função de polarização (6-31G*). Para obter as propriedades médias dos observaveis, cálculos de química quântica foram executados dentro da aproximação semi-empírica de INDO/S-CI. Com respeito ao espectro de absorção, os dados experimentais existentes na literatura científica reportam a existência de uma larga banda localizada na região de baixas energias, mais precisamente entre 400 e 600 nm. Nossos resultados teóricos para a forma alcalina mostram uma transição intensa transição π → π* aproximadamente à 432.4 ± 0.03 nm e, sob condições ácidas, esta transição aparece aproximadamente à 507.4 ± 0.12, 496.4 ± 0.28 ou 545.3 ± 0.10 nm, dependendo da estrutura, mostrando bom acordo com resultados experimentais. O DTP é um sistema particular usado na produção de polímeros de baixo gap. Suas propriedades elétricas e óticas foram obtidas através de um novo procedimento conhecido por Configuração Eletrostática Média do Solvent (ASEC). O procedimento ASEC inclui moléculas do solvente como cargas pontuais e permite o obtenção das quantidades quânticas executando somente poucos cálculos de mecânica quântica. Para o DTP, usando a teoria das perturbações de segunda ordem Mφller-Plesset (MP2) e o conjunto de bases aug-cc-pVDZ, a convergência do momento de dipolo foi alcançada com apenas quatro cálculos de mecânica quântica à 1.16 D, apresentando um aumento de 42% quando comparado ao dipolo isolado. O polarizabilidade corresponde à outra característica elétrica que pôde ser medida. Considerando o mesmo nível empregado ao cálculo do dipolo, o valor médio 132.7 a30 foi observado. A região de mais baixas da energias do espectro de absorption foi investigada também atravé de procedimento de ASEC usando ambos as aproximações, semi-empírico e DFT. Esta região de absorção é motivo de conclusões contraditórias com relação à natureza das transições n → π* e π → π*. Nossos resultados mostram que realmente que essas excitações são realmente observadas simultanemente podendo sobrepôr-se. Como exemplo, nossos resultados para DFT, encontrados usando o funcional B3LYP nos mostra que estas transições aparecem aproximadamente à 360.6 e 351.1 nm.

Cálculo das funções de anisotropia e radial de dose para braquiterapia com semente de 192Ir em diferentes materiais de interesse em radiologia com o auxílio do programa baseado no método de Monte Carlo MCNPX

The brachytherapy braquiterapia is an x-ray modality radiotherapy in which stamped or halfstamped radioactive sources in format of seeds are used, wires or to one short distance, in contact or implanted to the fabric to be treated. All the treatment modalities require a previous planning. The formalism recommended for calculation of dose was considered by the AAPM for the first TG-43 Report. In it distributions of dose of isolated seeds measured and calculated by Monte Carlo method in water instead of with models half-empiricists. In this work we in accordance with present some preliminary results of the calculation of functions of radial anisotropy and of dose in the distance for seed of 192Ir, wide used in brachytherapy treatments of high tax of dose (HDR), with the aid of the program based on the Monte Carlo method MCNPX v2.50 (Mount Carlo N ParticleeXtended). The materials chosen in the simulation beyond water, had been MS20 and estriado muscle

A simulação de Monte Carlo como instrumento fundamental para a análise de viabilidade econômica: a implementação de uma fundição de ferro em uma empresa metalúrgica

In uncertainly economic scenarios, an economic feasibility analysis must be done to accept a project based on investment criteria, such as NPV and IRR, mainly because the shareholders tend to invest their budget in a project if it has a great chance to return their investments. The studied company outsources all of its foundry services, what makes it very dependent of its suppliers, because the products have a lower quality level, long delivery periods and high prices. Therefore, this work will analyze the project of building an iron-casting foundry to decrease the dependence of its suppliers. In order to develop this research, all needed data related to the construction of the foundry and sales were collected to create deterministic and probabilistic (Monte Carlo Simulation) cash flows using MsExcel® and Oracle's Crystal Ball® software. As a result, the project was found to be risky by the NPV and IRR in the case of this new production line supplying only the internal needs. However, when the company offers its services to the foundry market, the project turns to be feasible

Estudo do método Monte Carlo e sua aplicação no modelo SIR

In this work we present a discussion and the results of the simulation of disease spread using the Monte Carlo method. The dissemination model is the SIR model and presents as main characteristic the disease evolution among individuals of the population subdivided into three groups: susceptible (S), infected (I) and recovered (R). The technique used is based on the introduction of transition probabilities S-> I and I->R to do the spread of the disease, they are governed by a Poisson distribution. The simulation of the spread of disease was based on the randomness introduced, taking into account two basic parameters of the model, the power of infection and average time of the disease. Considering appropriate values of these parameters, the results are presented graphically and analysis of these results gives information on a group of individuals react to the changes of these parameters and what are the chances of a disease becoming a pandemic

Application of reverse Monte Carlo simulations to diatomic molecules. Part 1. Noncomplete radial distribution functions

The reverse Monte Carlo (RMC) method generates sets of points in space which yield radial distribution functions (RDFS) that approximate those of the system of interest. Such sets of configurations should, in principle, be sufficient to determine the structural properties of the system. In this work we apply the RMC technique to fluids of hard diatomic molecules. The experimental RDFs of the hard-dimer fluid were generated by the conventional MC method and used as input in the RMC simulations. Our results indicate that the RMC method is only satisfactory in determining the local structure of the fluid studied by means of only mono-variable RDF. Also we suggest that the use of multi-variable RDFs would improve the technique significantly. However, the accuracy of the method turned out to be very sensitive to the variance of the input experimental RDF. © 1995.

Monte Carlo simulation of water-tetrahydrofuran mixtures

Monte Carlo simulations of water-tetrahydrofuran (THF) mixtures were performed in the isothermal-isobaric ensemble (NPT) at T = 298 K and p = 1 atm. The interaction energy was calculated using the TIP4P model for water and a five-site united atom representation for the THF molecule. The potential energy surfaces for water-THF interactions were obtained by using combining rules and the original potential functions used for pure liquids. Theoretical values obtained for the average interaction energy as a function of concentration are in good agreement with available experimental data. Results from the partitioning of the total interaction energy into water-water, water-THF and THF-THF contributions are presented. These results are useful to distinguish between the quantitative contributions of these molecular interactions to the energetic behavior of the water-THF mixing process. The radial distribution functions for HW-OTHF and OW-OTHF site-site interactions show the salient features of hydrogen-bonded liquids. Comparison of the average number of water-water complexes interacting through hydrogen bonding in water-THF and water-methanol mixtures shows an enhancement of the water-water coordination number in a THF rich environment. © 1995.

Using the Ferrenberg-Swendsen Monte Carlo method for nonuniversal exponents

We use the method of Monte Carlo evolution in the coupling constant space of Ferrenberg and Swendsen to evaluate the nonuniversal exponent η* associated to a linear defect in a 2d Ising model. © 1989.

Simulação de Monte Carlo: de modelos de spin à teoria de campos na rede

Pós-graduação em Física - IFT

Aplicação da Simulação de Monte Carlo no estudo de testes de fadiga em rodas agrícolas

Due to the increasing demand from clients and the search for better performances in the heavy vehicles industry, a progressive evolution in technology in a general way was needed. This paper uses a scientific method to validate, prior to its manufacture, the project of an agricultural wheel for sugar cane harvesters. Monte Carlo Simulation is used in conjunction with Finite Elements Method, in order to simulate the wheel's behavior in a cornering test, identify possible failure regions and get an estimate for its life under fatigue. To this end, test conditions according to EUWA Standards were simulated and obeyed, relevant to fatigue. Simulation results were interesting, according to industry experts involved in the project and manufacture of the product in question, and have provided important elements for the decision making regarding improvements that could be made on the product project before its execution

Evaluating Measurement Invariance with Censored Ordinal Data: A Monte Carlo Comparison of Alternative Model Estimators and Scales of Measurement

Evaluations of measurement invariance provide essential construct validity evidence. However, the quality of such evidence is partly dependent upon the validity of the resulting statistical conclusions. The presence of Type I or Type II errors can render measurement invariance conclusions meaningless. The purpose of this study was to determine the effects of categorization and censoring on the behavior of the chi-square/likelihood ratio test statistic and two alternative fit indices (CFI and RMSEA) under the context of evaluating measurement invariance. Monte Carlo simulation was used to examine Type I error and power rates for the (a) overall test statistic/fit indices, and (b) change in test statistic/fit indices. Data were generated according to a multiple-group single-factor CFA model across 40 conditions that varied by sample size, strength of item factor loadings, and categorization thresholds. Seven different combinations of model estimators (ML, Yuan-Bentler scaled ML, and WLSMV) and specified measurement scales (continuous, censored, and categorical) were used to analyze each of the simulation conditions. As hypothesized, non-normality increased Type I error rates for the continuous scale of measurement and did not affect error rates for the categorical scale of measurement. Maximum likelihood estimation combined with a categorical scale of measurement resulted in more correct statistical conclusions than the other analysis combinations. For the continuous and censored scales of measurement, the Yuan-Bentler scaled ML resulted in more correct conclusions than normal-theory ML. The censored measurement scale did not offer any advantages over the continuous measurement scale. Comparing across fit statistics and indices, the chi-square-based test statistics were preferred over the alternative fit indices, and ΔRMSEA was preferred over ΔCFI. Results from this study should be used to inform the modeling decisions of applied researchers. However, no single analysis combination can be recommended for all situations. Therefore, it is essential that researchers consider the context and purpose of their analyses.

Aplicação da Simulação de Monte Carlo no estudo de testes de fadiga em rodas agrícolas

Due to the increasing demand from clients and the search for better performances in the heavy vehicles industry, a progressive evolution in technology in a general way was needed. This paper uses a scientific method to validate, prior to its manufacture, the project of an agricultural wheel for sugar cane harvesters. Monte Carlo Simulation is used in conjunction with Finite Elements Method, in order to simulate the wheel's behavior in a cornering test, identify possible failure regions and get an estimate for its life under fatigue. To this end, test conditions according to EUWA Standards were simulated and obeyed, relevant to fatigue. Simulation results were interesting, according to industry experts involved in the project and manufacture of the product in question, and have provided important elements for the decision making regarding improvements that could be made on the product project before its execution

Characterization of an extrapolation chamber for low-energy X-rays: Experimental and Monte Carlo preliminary results

The extrapolation chamber is a parallel-plate ionization chamber that allows variation of its air-cavity volume. In this work, an experimental study and MCNP-4C Monte Carlo code simulations of an ionization chamber designed and constructed at the Calibration Laboratory at IFEN to be used as a secondary dosimetry standard for low-energy X-rays are reported. The results obtained were within the international recommendations, and the simulations showed that the components of the extrapolation chamber may influence its response up to 11.0%. (C) 2011 Elsevier Ltd. All rights reserved.

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.