Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution. (C) 2009 American Institute of Physics. [doi:10.1063/1.3253297]

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Monte Carlo investigation of the critical behavior of Stavskaya's probabilistic cellular automaton

Stavskaya's model is a one-dimensional probabilistic cellular automaton (PCA) introduced in the end of the 1960s as an example of a model displaying a nonequilibrium phase transition. Although its absorbing state phase transition is well understood nowadays, the model never received a full numerical treatment to investigate its critical behavior. In this Brief Report we characterize the critical behavior of Stavskaya's PCA by means of Monte Carlo simulations and finite-size scaling analysis. The critical exponents of the model are calculated and indicate that its phase transition belongs to the directed percolation universality class of critical behavior, as would be expected on the basis of the directed percolation conjecture. We also explicitly establish the relationship of the model with the Domany-Kinzel PCA on its directed site percolation line, a connection that seems to have gone unnoticed in the literature so far.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

Triple- and quadruple-escape peaks in HPGe detectors: Experimental observation and Monte Carlo simulation

The triple- and quadruple-escape peaks of 6.128 MeV photons from the (19)F(p,alpha gamma)(16)O nuclear reaction were observed in an HPGe detector. The experimental peak areas, measured in spectra projected with a restriction function that allows quantitative comparison of data from different multiplicities, are in reasonably good agreement with those predicted by Monte Carlo simulations done with the general-purpose radiation-transport code PENELOPE. The behaviour of the escape intensities was simulated for some gamma-ray energies and detector dimensions; the results obtained can be extended to other energies using an empirical function and statistical properties related to the phenomenon. (C) 2010 Elsevier B.V. All rights reserved.

Extending the use of canonical and microcanonical Monte Carlo algorithms to spin models

We describe the canonical and microcanonical Monte Carlo algorithms for different systems that can be described by spin models. Sites of the lattice, chosen at random, interchange their spin values, provided they are different. The canonical ensemble is generated by performing exchanges according to the Metropolis prescription whereas in the microcanonical ensemble, exchanges are performed as long as the total energy remains constant. A systematic finite size analysis of intensive quantities and a comparison with results obtained from distinct ensembles are performed and the quality of results reveal that the present approach may be an useful tool for the study of phase transitions, specially first-order transitions. (C) 2009 Elsevier B.V. All rights reserved.

Modelling water adsorption on Au(210) surfaces: II. Monte Carlo simulations

Canonical Monte Carlo simulations for the Au(210)/H(2)O interface, using a force field recently proposed by us, are reported. The results exhibit the main features normally observed in simulations of water molecules in contact with different noble metal surfaces. The calculations also assess the influence of the surface topography on the structural aspects of the adsorbed water and on the distribution of the water molecules in the direction normal to the metal surface plane. The adsorption process is preferential at sites in the first layer of the metal. The analysis of the density profiles and dipole moment distributions points to two predominant orientations. Most of the molecules are adsorbed with the molecular plane parallel to surface, while others adsorb with one of the O-H bonds parallel to the surface and the other bond pointing towards the bulk liquid phase. There is also evidence of hydrogen bond formation between the first and second solvent layers at the interface. (c) 2007 Elsevier B.V. All rights reserved.

Modelos de volatilidade estocástica em séries financeiras: uma aplicação para o IBOVESPA

Neste artigo apresentamos uma análise Bayesiana para o modelo de volatilidade estocástica (SV) e uma forma generalizada deste, cujo objetivo é estimar a volatilidade de séries temporais financeiras. Considerando alguns casos especiais dos modelos SV usamos algoritmos de Monte Carlo em Cadeias de Markov e o software WinBugs para obter sumários a posteriori para as diferentes formas de modelos SV. Introduzimos algumas técnicas Bayesianas de discriminação para a escolha do melhor modelo a ser usado para estimar as volatilidades e fazer previsões de séries financeiras. Um exemplo empírico de aplicação da metodologia é introduzido com a série financeira do IBOVESPA.

Obtaining high purity silica from rice hulls

Many routes for extracting silica from rice hulls are based on direct calcining. These methods, though, often produce silica contaminated with inorganic impurities. This work presents the study of a strategy for obtaining silica from rice hulls with a purity level adequate for applications in electronics. The technique is based on two leaching steps, using respectively aqua regia and Piranha solutions, which extract the organic matrix and inorganic impurities. The material was characterized by Fourier-transform infrared spectroscopy (FTIR), powder x-ray diffraction (XRD), x-ray fluorescence (XRF), scanning electron microscopy (SEM), particle size analysis by laser diffraction (LPSA) and thermal analysis.

A new species of Acteon (Opisthobranchia: Acteonidae) from Northeast Brazil

A new species of Acteon Montfort 1810, Acteon mirim sp. nov., from Canopus Bank, state of Ceará, Brazil is described based on shell morphology. The new species is compared with other species of the genus reported from Brazil. It differs from other Brazilian species in having a whitish color with dark orange-brown spiral bands and a shell surface covered with small spiral grooves, regularly rectangular.

Prevalência de anticorpos anti-Chlamydophila spp. em propriedades rurais com histórico de aborto bovino no estado do Paraná

Chlamydophila abortus é o agente etiológico do aborto epizoótico bovino, cujas manifestações clínicas mais freqüentes são aborto, nascimento de bezerros prematuros e de animais fracos, natimorto e repetição de cio em intervalos irregulares. O objetivo deste trabalho foi estimar a prevalência de anticorpos anti-Chlamydophila spp. em fêmeas bovinas de propriedades rurais com histórico de aborto, selecionadas dentro do delineamento amostral do Plano Nacional de Controle e Erradicação da Brucelose e Tuberculose no estado do Paraná. Foram testadas pela prova de fixação de complemento 3.102 amostras de soro de fêmeas bovinas (idade > 24 meses), provenientes de 373 propriedades. Ao total, 44 (1,42%) animais foram positivos com títulos > 32. A prevalência de focos foi de 8,82% (6,15%-12,17%). Animais confinados ou semi-confinados (OR=3.339, P=0.004), propriedade com menos de 35 matrizes (OR=3.339, P=0.017), presença de produtos do aborto na pastagem (OR=2.372, P=0.037) e aluguel de pasto (OR=3.398, P=0.006) foram considerados fatores de risco para Chlamydophila spp. A infecção por Chlamydophila spp. acometeu um número pequeno de animais, oriundos de propriedades com histórico de aborto. A importância deste agente como causa de aborto em bovinos no estado do Paraná, se existir, é muito pequena.

Troponina cardíaca T para estratificação de risco na insuficiência cardíaca crônica descompensada

FUNDAMENTO: As troponinas cardíacas são marcadores altamente sensíveis e específicos de lesão miocárdica. Esses marcadores foram detectados na insuficiência cardíaca (IC) e estão associadas com mau prognóstico. OBJETIVO: Avaliar a relação da troponina T (cTnT) e suas faixas de valores com o prognóstico na IC descompensada. MÉTODOS: Estudaram-se 70 pacientes com piora da IC crônica que necessitaram de hospitalização. Na admissão, o modelo de Cox foi utilizado para avaliar as variáveis capazes de predizer o desfecho composto por morte ou re-hospitalização em razão de piora da IC durante um ano. RESULTADOS: Durante o seguimento, ocorreram 44 mortes, 36 re-hospitalizações por IC e 56 desfechos compostos. Na análise multivariada, os preditores de eventos clínicos foram: cTnT (cTnT > 0,100 ng/ml; hazard ratio (HR) 3,95 intervalo de confiança (IC) 95%: 1,64-9,49, p = 0,002), diâmetro diastólico final do ventrículo esquerdo (DDVE >70 mm; HR 1,92, IC95%: 1,06-3,47, p = 0,031) e sódio sérico (Na <135 mEq/l; HR 1,79, IC95%: 1,02-3,15, p = 0,044). Para avaliar a relação entre a elevação da cTnT e o prognóstico na IC descompensada, os pacientes foram estratificados em três grupos: cTnT-baixo (cTnT < 0,020 ng/ml, n = 22), cTnT-intermediário (cTnT > 0,020 e < 0,100 ng/ml, n = 36) e cTnT-alto (cTnT > 0,100 ng/ml, n = 12). As probabilidades de sobrevida e sobrevida livre de eventos foram: 54,2%, 31,5%, 16,7% (p = 0,020), e 36,4%, 11,5%, 8,3% (p = 0,005), respectivamente. CONCLUSÃO: A elevação da cTnT está associada com mau prognóstico na IC descompensada, e o grau dessa elevação pode facilitar a estratificação de risco