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.

Chinese Basic Pension Substitution Rate: A Monte Carlo Demonstration of the Individual Account Model

At the end of 2005, the State Council of China passed ”The Decision on adjusting the Individual Account of Basic Pension System”, which adjusted the individual account in the 1997 basic pension system. In this essay, we will analyze the adjustment above, and use Life Annuity Actuarial Theory to establish the basic pension substitution rate model. Monte Carlo simulation is also used to prove the rationality of the model. Some suggestions are put forward associated with the substitution rate according to the current policy.

Forecasting accuracy and estimation uncertainty using VAR models with short- and long-term economic restrictions: a Monte-Carlo study

Using vector autoregressive (VAR) models and Monte-Carlo simulation methods we investigate the potential gains for forecasting accuracy and estimation uncertainty of two commonly used restrictions arising from economic relationships. The Örst reduces parameter space by imposing long-term restrictions on the behavior of economic variables as discussed by the literature on cointegration, and the second reduces parameter space by imposing short-term restrictions as discussed by the literature on serial-correlation common features (SCCF). Our simulations cover three important issues on model building, estimation, and forecasting. First, we examine the performance of standard and modiÖed information criteria in choosing lag length for cointegrated VARs with SCCF restrictions. Second, we provide a comparison of forecasting accuracy of Ötted VARs when only cointegration restrictions are imposed and when cointegration and SCCF restrictions are jointly imposed. Third, we propose a new estimation algorithm where short- and long-term restrictions interact to estimate the cointegrating and the cofeature spaces respectively. We have three basic results. First, ignoring SCCF restrictions has a high cost in terms of model selection, because standard information criteria chooses too frequently inconsistent models, with too small a lag length. Criteria selecting lag and rank simultaneously have a superior performance in this case. Second, this translates into a superior forecasting performance of the restricted VECM over the VECM, with important improvements in forecasting accuracy ñreaching more than 100% in extreme cases. Third, the new algorithm proposed here fares very well in terms of parameter estimation, even when we consider the estimation of long-term parameters, opening up the discussion of joint estimation of short- and long-term parameters in VAR models.

The importance of common cyclical features in VAR analysis: a Monte-Carlo study

Despite the commonly held belief that aggregate data display short-run comovement, there has been little discussion about the econometric consequences of this feature of the data. We use exhaustive Monte-Carlo simulations to investigate the importance of restrictions implied by common-cyclical features for estimates and forecasts based on vector autoregressive models. First, we show that the ìbestî empirical model developed without common cycle restrictions need not nest the ìbestî model developed with those restrictions. This is due to possible differences in the lag-lengths chosen by model selection criteria for the two alternative models. Second, we show that the costs of ignoring common cyclical features in vector autoregressive modelling can be high, both in terms of forecast accuracy and efficient estimation of variance decomposition coefficients. Third, we find that the Hannan-Quinn criterion performs best among model selection criteria in simultaneously selecting the lag-length and rank of vector autoregressions.

The importance of Common-Cyclical Features in VAR analysis: a Monte-Carlo study (Preliminary Version)

Despite the belief, supported byrecentapplied research, thataggregate datadisplay short-run comovement, there has been little discussion about the econometric consequences ofthese data “features.” W e use exhaustive M onte-Carlo simulations toinvestigate theimportance ofrestrictions implied by common-cyclicalfeatures for estimates and forecasts based on vectorautoregressive and errorcorrection models. First, weshowthatthe“best” empiricalmodeldevelopedwithoutcommoncycles restrictions neednotnestthe“best” modeldevelopedwiththoserestrictions, duetothe use ofinformation criteria forchoosingthe lagorderofthe twoalternative models. Second, weshowthatthecosts ofignoringcommon-cyclicalfeatures inV A R analysis may be high in terms offorecastingaccuracy and e¢ciency ofestimates ofvariance decomposition coe¢cients. A lthough these costs are more pronounced when the lag orderofV A R modelsareknown, theyarealsonon-trivialwhenitis selectedusingthe conventionaltoolsavailabletoappliedresearchers. T hird, we…ndthatifthedatahave common-cyclicalfeatures andtheresearcherwants touseaninformationcriterium to selectthelaglength, theH annan-Q uinn criterium is themostappropriate, sincethe A kaike and theSchwarz criteriahave atendency toover- and under-predictthe lag lengthrespectivelyinoursimulations.

Estimation of DSGE Models: A Monte Carlo Analysis

Neste trabalho investigamos as propriedades em pequena amostra e a robustez das estimativas dos parâmetros de modelos DSGE. Tomamos o modelo de Smets and Wouters (2007) como base e avaliamos a performance de dois procedimentos de estimação: Método dos Momentos Simulados (MMS) e Máxima Verossimilhança (MV). Examinamos a distribuição empírica das estimativas dos parâmetros e sua implicação para as análises de impulso-resposta e decomposição de variância nos casos de especificação correta e má especificação. Nossos resultados apontam para um desempenho ruim de MMS e alguns padrões de viés nas análises de impulso-resposta e decomposição de variância com estimativas de MV nos casos de má especificação considerados.

Utilização do método de Monte Carlo na resolução de problemas de transferência de calor por radiação em cavidades que contém meio participante

Um estudo foi realizado com a finalidade principal de investigar de que maneira a teoria de transferência de calor por radiação em meios participantes pode ser melhor aplicada na resolução de problemas de troca de energia no interior de cavidades que se assemelhem a fornalhas industriais. Verificou-se que as superfícies internas da maioria das fornalhas podem ser consideradas cinzentas e difusoras, o que facilita a aplicação de modelos de soma ponderada de gases cinzas. Sendo assim, as maiores dificuldades aparecem quando a geometria da cavidade é complexa. Com o método de Monte Carlo, é possível resolver problemas com características geométricas bastante complexas, bem como considerar efeitos espectrais e direcionais, quando necessário. Este método foi, então, aplicado, em conjunto com o método da zona, na elaboração de três códigos computacionais, que são capazes de determinar áreas de troca totais em cavidades onde as paredes e o gás são cinzas. Um deles pode ser aplicado a cavidades paralelepipédicas; outro, a cilíndricas. Tanto no primeiro, quanto no segundo, as dimensões podem assumir qualquer valor. Com o terceiro código computacional, é possível obter-se áreas de troca totais entre zonas de cavidades que têm geometria genérica formada a partir de combinações de cubos. Resultados obtidos através dos três códigos computacionais foram comparados com outros disponíveis na literatura.

A new algorithm for Reverse Monte Carlo simulations

We present a new algorithm for Reverse Monte Carlo (RMC) simulations of liquids. During the simulations, we calculate energy, excess chemical potentials, bond-angle distributions and three-body correlations. This allows us to test the quality and physical meaning of RMC-generated results and its limitations. It also indicates the possibility to explore orientational correlations from simple scattering experiments. The new technique has been applied to bulk hard-sphere and Lennard-Jones systems and compared to standard Metropolis Monte Carlo results. (C) 1998 American Institute of Physics.

Monte Carlo investigations of intermolecular interactions in water-amide mixtures

Monte Carlo simulations of water-amides (amide=fonnamide-FOR, methylfonnamide-NMF and dimethylformamide-DMF) solutions have been carried out in the NpT ensemble at 308 K and 1 atm. The structure and excess enthalpy of the mixtures as a function of the composition have been investigated. The TIP4P model was used for simulating water and six-site models previously optimized in this laboratory were used for simulating the liquid amides. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones potential plus a Coulomb term. The interaction energy between solute and solvent has been partitioned what leads to a better understanding of the behavior of the enthalpy of mixture obtained for the three solutions experimentally. Radial distribution functions for the water-amides correlations permit to explore the intermolecular interactions between the molecules. The results show that three, two and one hydrogen bonds between the water and the amide molecules are formed in the FOR, NMF and DMF-water solutions, respectively. These H-bonds are, respectively, stronger for DMF-water, NMF-water and FOR-water. In the NMF-water solution, the interaction between the methyl group of the NMF and the oxygen of the water plays a role in the stabilization of the aqueous solution quite similar to that of an H-bond in the FOR-water solution. (c) 2005 Elsevier B.V. All rights reserved.

Structure of acetone and dimethyl sulfoxide from Monte Carlo simulations and MM2 calculations

The structure of acetone and dimethyl sulfoxide in the liquid phase is investigated using Monte Carlo simulations and MM2 calculations. The principal site - site correlations and degree of structure in both liquids have been investigated. The results showed that dimethyl sulfoxide is more structured than acetone. At short distances the dipoles of neighboring molecules are found to be in antiparallel configurations, but further apart the molecules tend to be aligned predominantly as head to tail. In both liquids there is evidence of strong methyl - oxygen interaction, important to the structure of the liquids. The contacts suggest weak hydrogen bonds between methyl hydrogen and oxygen.

A Monte Carlo Revisiting of N-Methylformamide and Acetone

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Monte Carlo studies of N-methylformamide-dimethyl sulfoxide mixtures

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)