Conveniencia económica del empleo de seguros para riesgos climáticos en la producción agrícola : un análisis mediante simulación Monte Carlo

p.157-168

Modelación agropecuaria mediante programación lineal a partir de modelos Monte Carlo para el partido de Guaminí (Pcia. de Buenos Aires)

p.147-154

Diseño de agromodelos sostenibles para el partido de Guaminí (pcia. de Buenos Aires) con el método Monte Carlo

p.119-129

Neutron quasi-elastic scattering in disordered solids: a Monte Carlo study of metal-hydrogen systems

The dynamic structure factor of neutron quasi-elastic scattering has been calculated by Monte Carlo methods for atoms diffusing on a disordered lattice. The disorder includes not only variation in the distances between neighbouring atomic sites but also variation in the hopping rate associated with each site. The presence of the disorder, particularly the hopping rate disorder, causes changes in the time-dependent intermediate scattering function which translate into a significant increase in the intensity in the wings of the quasi-elastic spectrum as compared with the Lorentzian form. The effect is particularly marked at high values of the momentum transfer and at site occupancies of the order of unity. The MC calculations demonstrate how the degree of disorder may be derived from experimental measurements of the quasi-elastic scattering. The model structure factors are compared with the experimental quasi-elastic spectrum of an amorphous metal-hydrogen alloy.

Rotating-frame nuclear magnetic relaxation of spins diffusing on a disordered lattice: a Monte Carlo model

The rotating-frame nuclear magnetic relaxation rate of spins diffusing on a disordered lattice has been calculated by Monte Carlo methods. The disorder includes not only variation in the distances between neighbouring spin sites but also variation in the hopping rate associated with each site. The presence of the disorder, particularly the hopping rate disorder, causes changes in the time-dependent spin correlation functions which translate into asymmetry in the characteristic peak in the temperature dependence of the dipolar relaxation rate. The results may be used to deduce the average hopping rate from the relaxation but the effect is not sufficiently marked to enable the distribution of the hopping rates to be evaluated. The distribution, which is a measure of the degree of disorder, is the more interesting feature and it has been possible to show from the calculation that measurements of the relaxation rate as a function of the strength of the radiofrequency spin-locking magnetic field can lead to an evaluation of its width. Some experimental data on an amorphous metal - hydrogen alloy are reported which demonstrate the feasibility of this novel approach to rotating-frame relaxation in disordered materials.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Monte Carlo simulation of cluster growth in surface defects induced by the tip of a scanning tunnelling microscope

First steps are taken to model the electrochemical deposition of metals in nanometer-sized cavities. In the present work, the electrochemical deposition of Cu atoms in nanometer-sized holes dug on Au(111) is investigated through Monte Carlo simulations using the embedded atom method to represent particle interactions. By sweeping the chemical potential of Cu, a cluster is allowed to grow within the hole rising four atomic layers above the surface. Its lateral extension remains confined to the area defined by the borders of the original defect. (C) 2004 Elsevier B.V. All rights reserved.

Monte Carlo simulation of properties of monolayers and metal islands adsorbed on metallic (111) surfaces

To obtain the surface stress changes due to the adsorption of metal monolayers onto metallic surfaces, a new model derived from thermodynamic considerations is presented. Such a model is based on continuum Monte Carlo simulations with embedded atom method potentials in the canonical ensemble, and it is extended to consider the behavior on different islands adsorbed onto (111) substrate surfaces. Homoepitaxial and heteroepitaxial systems are studied. Pseudomorphic growth is not observed for small metal islands with considerable positive misfit with the substrate. Instead, the islands become compressed upon increase of their size. A simple model is proposed to interpolate between the misfits of atoms in small islands and the pseudomorphic behavior of the monolayer.