1:45 PM - 2:15 PM
[SMP33-12] Molecular Dynamics Study of Soret Effect in Silicate Liquids
★Invited Papers
Keywords:Silicate Liquids, Molecular Dynamics Simulation, Soret Effect
I applied the molecular dynamics (MD) simulations for pseudo binary and aluminosilicate liquids to investigate the mechanism of fractionation by temperature gradient. The MD simulation appropriate method for this study because it gives trajectories of each atoms in the simulation cells and potentials of each atom at any point in simulation time. MD simulations of sodium/calcium-silicate/aluminosilicate liquids were performed using MXDORTO code (Sakuma & Kawamura, 2009). The simulation conditions are as follows: The inter-atomic potential model was taken from Feuston & Garofalini (1988). Systems of approximately 30000 particles in rod-shaped (approximately 5 x 5 x 17 nm) simulation cell in periodic boundary condition were firstly annealed for 2 ns at 2000 K from randomly generated structure. Then we started the simulations in temperature gradient. The temperatures in sliced region (approximately 0.35 nm thickness) perpendicular to the longest axis at the end and the center of simulation cell were maintained 2000 and 3500 K using the scaling procedure. After several tens nano-second simulations, we confirmed the changes in distribution of composition. The concentration of SiO2 in the high-temperature center part slightly increases as simulation proceeds. In contrast the concentration of CaO and Na2O in that part slightly decreases. Soret coefficients for each cation has linear relation with the half width of potential energy distribution in binary silicates. The existence of Al2O3 content impedes fractionation in sodium-aluminosilicate whereas it enhances in calcium-aluminosilicate. This behavior might be caused by the charge compensation of 4-fold aluminum cation.