Simulação via dinâmica molecular do Fluoreto de Bário (BaF2)

Abstract

The system studied is considered a superionic material type AB2. This shows a variety of interesting properties, not only in academic context, but also in technological applications. However, there are not still many studies about its behavior and dynamics. The computational method used in this work is the method of Molecular Dynamics (MD), a technique able to generate the trajectories in phase space of N particles that interacting with each other. The system in question is which is composed of 1500 particles, with 1000 of Barium and 500 of Fluoride. The …rst step is to provide simulation system its initial conditions and from there use an algorithm, appropriate to the system under study, to make the numerical integrations of the forces and the calculations of the positions and velocities of the particles. Complied all these steps is then started the simulation in fact. The system under consideration in steps of illness started time _t as the addition initially temperature while the system pressure was kept constant and equal to zero. The ensembles isoentalpic-isobaric (N, P, H), in which the volume ensemble acts as a dynamic variable, and ensemble microcanonic were used (N, V, E), where the calculations of statistical averages were performed. The potential chosen to describe our system has the potential pairs proposed by Vashishta-Rahman. From there we started the simulation via DM Barium ‡uoride BaF2. The temperature of the system was to vary from an initial temperature of 300K to a …nal temperature of 3000K in increments of 50K. From there, they were possible to plot graphs with T dependence of this parameter as well as the graphs of the pair correlation g(r) and coordination number c(r), the density function of vibrational states G(w) functions, and peaks selected from G(w) function that shows the relation of the frequency temperature T. The next step in the simulation was to investigate, from the application of pressure, as the system would behave. The system su¤ered compression 0GPa 3GPa up in increments of 0.5 GPa and temperature of 300K. The superionic phase of the compound is con…rmed and its melting temperature.