Monte-Carlo simulation of the phase diagrams in Ising model of [Sn(Pb)]_2p_2[S(Se)]_6 crystals

摘要

In Sn_2P_2S_6 crystals ferroelectric second order phase transition (PT) occurs. In this case the number of formula units in elementary cell stays unchanged. In Sn_2P_2Se_6 crystals the intermediate incommensurate phase exists. A dipole Ising model of Sn_2P_2[S(Se)]_6 ferroelectrics has been developed. To consider the dipole long-range interaction around the dipole an area of simple form is singled out. Within this area the interaction between the dipoles is calculated precisely. The interaction with the other dipoles is considered by introduction of a local field. For better accordance with the real experiment an amendment on the short-range interaction with two nearest neighbours has been introduced. As a result of calculations made by Monte-Carlo method it has been found that when a certain short-range interaction coefficient in Sn_2P_2S_6 is selected the intermediate modulated phase disappears. In this case the temperatures of PTs are as follows: for Sn_2P_2S_6 model T_0* approx 10.3;;for Sn_2P_2Se_6 model T_c* approx 8.7, T_i* approx 9.7 (in the effective temperatures). These data are in satisfactory agreement with the results of the real experiment. The concentration dependences of temperatures of PTs in Sn_2P_2(S_xSe_(1-x)0_6 crystals are simulated and the Lifshitz point was observed. In mixed (Pb_ySn_(1-y))_2P_2S_6 and (Pb_ySn_(1-y))_2P_2S_6 crystals increase of Pb concentration causes reduction of temperatures of PTs to OK. In this case the character of PTs stays unchanged. By the method of diluted Ising lattice the concentration dependences of temperatures of PTs in this crystals are simulated.