A critical investigation of the Tanford-Kirkwood scheme by means of Monte Carlo simulations

Monte Carlo simulations are used to assess the adequacy of the Tanford-Kirkwood prescription for electrostatic interactions in macromolecules. Within a continuum dielectric framework, the approach accurately describes salt screening of electrostatic interactions for moderately charged systems consistent with common proteins at physiological conditions. The limitations of the Debye-Huckel theory, which forms the statistical mechanical basis for the Tanford-Kirkwood result, become apparent for highly charged systems. It is shown, both by an analysis of the Debye-Huckel theory and by numerical simulations, that the difference in dielectric permittivity between macromolecule and surrounding solvent does not play a significant role for salt effects if the macromolecule is highly charged. By comparison to experimental data, the continuum dielectric model (combined with either an approximate effective Hamiltonian as in the Tanford-Kirkwood treatment or with exact Monte Carlo simulations) satisfactorily predicts the effects of charge mutation on metal ion binding constants, but only if the macromolecule and solvent are assigned the same or similar permittivities.

Application of a new reverse Monte Carlo algorithm to polyatomic molecular systems. I. Liquid water

Using a new reverse Monte Carlo algorithm, we present simulations that reproduce very well several structural and thermodynamic properties of liquid water. Both Monte Carlo, molecular dynamics simulations and experimental radial distribution functions used as input are accurately reproduced using a small number of molecules and no external constraints. Ad hoc energy and hydrogen bond analysis show the physical consistency and limitations of the generated RMC configurations. (C) 2001 American Institute of Physics.

MONTE-CARLO SIMULATION OF POLYAMPHOLYTE CHAINS

Polyampholyte copolymers containing both positive and negative monomers regularly dispersed along the chain were studied. The Monte Carlo method was used to simulate chains with charged monomers interacting by screened Coulomb potential. The neutral polyampholyte chains collapse due to the attractive electrostatic interactions. The nonneutral chains are in extended conformations due to the repulsive polyelectrolyte effects that dominate the attractive polyampholyte interactions. The results are in good agreement with experiment.

C - H ⋯ 0 and N - H ⋯ 0 hydrogen bonds in liquid amides investigated by monte carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N, N-dimethytformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ⋯ O and N-H ⋯ O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ⋯ O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This interaction is particularly important in the structure of MF. The intensity of the N - H ⋯ O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. © 1997 John Wiley & Sons, Inc.

Coil to extended coil transition in polygalacturonic acid: Conductometric titration and monte carlo simulations

The conformational transition from coil to extended coil for polygalacturonic acid has been studied by conductometric titrations and Monte Carlo simulations. The results of conductometric titrations at different polymer concentrations have been analyzed using the model proposed by Manning,1 which describes the conductivity of polyelectrolitic solutions. This experimental approach provides the transport factor and the average distance between charged groups at different degrees of ionization (α). The mean distances between charged groups have been compared with the values obtained by Monte Carlo simulations. In these simulations the polymer chain is modeled as a self-avoiding random walk in a cubic lattice. The monomers interact through the unscreened Coulombic potential. The ratio between the end-to-end distance and the number of ionized beads provides the average distance between charged monomers. The experimental and theoretical values are in good agreement for the whole range of ionization degrees accessed by conductometric titrations. These results suggest that the electrostatic interactions seem to be the major contribution for the coil to extended coil conformational change. The small deviations for α ≤ 0.5 suggests that the stiffness of the chain, associated with local interactions, becomes increasingly significant as the fraction of charged groups is decreased. © 2000 American Chemical Society.

Monte Carlo and analytical calculation of lateral deflection of proton beams in homogeneous targets

Proton radiation therapy is a precise form of radiation therapy, but the avoidance of damage to critical normal tissues and the prevention of geographical tumor misses require accurate knowledge of the dose delivered to the patient and the verification of his position demand a precise imaging technique. In proton therapy facilities, the X-ray Computed Tomography (xCT) is the preferred technique for the planning treatment of patients. This situation has been changing nowadays with the development of proton accelerators for health care and the increase in the number of treated patients. In fact, protons could be more efficient than xCT for this task. One essential difficulty in pCT image reconstruction systems came from the scattering of the protons inside the target due to the numerous small-angle deflections by nuclear Coulomb fields. The purpose of this study is the comparison of an analytical formulation for the determination of beam lateral deflection, based on Molière's theory and Rutherford scattering with Monte Carlo calculations by SRIM 2008 and MCNPX codes. © 2010 American Institute of Physics.

Monte Carlo simulation applied in total reflection X-ray fluorescence: Preliminary results

The X-ray Fluorescence (XRF) analysis is a technique for the qualitative and quantitative determination of chemical constituents in a sample. This method is based on detection of the characteristic radiation intensities emitted by the elements of the sample, when properly excited. A variant of this technique is the Total Reflection X-ray Fluorescence (TXRF) that utilizes electromagnetic radiation as excitation source. In total reflection of X-ray, the angle of refraction of the incident beam tends to zero and the refracted beam is tangent to the sample support interface. Thus, there is a minimum angle of incidence at which no refracted beam exists and all incident radiation undergoes total reflection. In this study, we evaluated the influence of the energy variation of the beam of incident x-rays, using the MCNPX code (Monte Carlo NParticle) based on Monte Carlo method. © 2013 AIP Publishing LLC.

Estudo de percolação de clusters de Monte Carlo para o modelo de Ising bidimensional

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

Estudo de Monte Carlo da produção de nova ressonância neutra decaindo em dois bósons vetoriais no experimento CMS

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

Teoria quântica do campo escalar real com autoacoplamento quártico - simulações de Monte Carlo na rede com um algoritmo worm

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Regressão Simplex aplicada a delineamentos de mistura e utilização do Algoritmo Boosting

In the composition of this work are present two parts. The first part contains the theory used. The second part contains the two articles. The first article examines two models of the class of generalized linear models for analyzing a mixture experiment, which studied the effect of different diets consist of fat, carbohydrate, and fiber on tumor expression in mammary glands of female rats, given by the ratio mice that had tumor expression in a particular diet. Mixture experiments are characterized by having the effect of collinearity and smaller sample size. In this sense, assuming normality for the answer to be maximized or minimized may be inadequate. Given this fact, the main characteristics of logistic regression and simplex models are addressed. The models were compared by the criteria of selection of models AIC, BIC and ICOMP, simulated envelope charts for residuals of adjusted models, odds ratios graphics and their respective confidence intervals for each mixture component. It was concluded that first article that the simplex regression model showed better quality of fit and narrowest confidence intervals for odds ratio. The second article presents the model Boosted Simplex Regression, the boosting version of the simplex regression model, as an alternative to increase the precision of confidence intervals for the odds ratio for each mixture component. For this, we used the Monte Carlo method for the construction of confidence intervals. Moreover, it is presented in an innovative way the envelope simulated chart for residuals of the adjusted model via boosting algorithm. It was concluded that the Boosted Simplex Regression model was adjusted successfully and confidence intervals for the odds ratio were accurate and lightly more precise than the its maximum likelihood version.

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.

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 teórico da contribuição a dose de partículas secundárias geradas por prótons de energias entre 100-200 MeV num phantom de músculo estriado

In the treatment plans in conventional Proton therapy are considered only the elastic interactions of protons with electrons and/or nuclei, it means, mainly ionization and coulomb excitation processes. As the energy needed to reach the deep tumors should be of several hundred of MeVs, certainly the nuclear inelastic channels are open. Only some previous studies of the contribution of these processes in the full dose have been made towards targets composed of water. In this study will be presented the results of the simulation of the processes of interaction of beams of protons in the range of 100-200 MeV of energy with a cylindrical phantom composed by striated muscle (ICRU), emphasizing in the contribution to total dose due to the deposition of energy by secondary particles alpha (α), deuterium (2H), tritium (3H), neutron (n) and hélio3 (3He), originated by nuclear inelastic processes. The simulations were performed by using the method of Monte Carlo, via the computer code MCNPX v2.50 (Monte Carlo N-Particle eXtended). The results will be shown demonstrated through the graphics of the deposited dose with or without nuclear interaction, the percentual of dose deposited by secondary particles, the radial dispersion of neutrons, as well as the multiplicity of secondary particles