Output Files -

*mcphas.qvc*- the set of test q-vectors used for calculation of free energy. Components of these q vectors refer to the reciprocal lattice .
*mcphas.phs**spin-configuration table of different types of spin-configurations.**Note*: there is no natural criteria for deciding, if one spin-configuration is different from another one. Therefore the list of ''different'' spin-configurations is dependent on the meaning of ''different''.*The program**McPhase*decides whether a spin-configuration is different from another by a simple criteria, namely by the angle between the spins. Comparing two spin configurations it calculates the angle between corresponding spins and if for one spin the angle is not small, the configuration is treated as a different configuration. Therefore for example a ferromagnet with moments in has a different spin configuration than a ferromagnet with moments in direction.*mcphas.sps**dependence of spin-configuration. The spin configurations stored in this file may be displayed using the program**spins*, an example is given in figure 18. The expectation values ( , note that for the single ion module*so1ion*, for the complete list see appendix G) are stored in this file for i=1,2,...,number of spins in magnetic unit cell.*mcphas.mf**dependence of exchange field configuration, stored as (unit is in meV) for i=1,2,...,number of spins in magnetic unit cell.**mcphas.fum**free energy according to equation (18), magnetic energy per site*

(20) *Note, that form the numerical derivative of with respect to temperature the specific heat per site in mean field approximation can be calculated . Finally,**mcphas.fum*contains the calculated magnetisation (only if a*mcalc*function is present in the single ion module). As an example for the information contained in this file the calculated magnetisation of NdCu is shown in fig 9. Note, in case of non-orthogonal axes the convention for applied field and also for the magnetisation components in these tables is , and perpendicular to and .*mcphas1.j1 .j1 .j2 ...**spin-spin correlation functions for sub-lattice 1 neighbour 1 2 ... (linear combination is proportional to magnetostriction) The spin-spin correlation functions for neighbour are defined by the following sum of dyadic products:*

(21) *with being the number of magnetic atoms in the magnetic unit cell. For example for the single ion module**so1ion*, for the complete list see appendix G. Single ion and two-ion magnetostriction can be calculated using the from*mcphas.xyt*and the spin-spin correlation functions from*mcphas1.j1 .j1 .j2 ...*. As an example the magnetostriction analysis of NdCu is shown in figure 10. For details please refer to [16]. Note, in case of non-orthogonal axes the convention for applied field and also for the moment components in these tables is , and perpendicular to and .*mcphas.xyt*- phase diagram as x,y,T, H, phase-number j according to spin-configuration table given in mcphas.phs, a periodicity key, thermal expectationvalues (averaged over the unit cell). Figure 21 shows the phase diagram of NdCu for magnetic fields parallel to the orthorhombic -direction. Note, in case of non-orthogonal axes the convention for applied field in these tables is , and perpendicular to and . Note that the phase-number is not necessarily a unique number for each phase in the phase diagram, because the program numbers phases according to their similarity in the moments. For example, any tetragonal easy plane magnet will have domains with moments parallel to and . The program will generate both domains and find equal energy for both. It will give each domain a separate phase number, because spin directions are completely different. Because both domains have the same free energy at every temperature, the phase-index will vary in a random manner. In order to generate a beautiful phase diagram figure, such equivalent domains have to be identified and phase numbers have to be substituted.
*mcphas.hkl**calculated (unpolarised) neutron diffraction data (the calculated magnetic intensities correspond to the magnetic structure + Polarisation factor + magnetic Formfactor (dipole approx). The Lorentz-factor and instrumental corrections are not calculated.). The single ion module function**mcalc*is used to evaluate the magnetic moment on each site. As an example figure 20 shows the calculated temperature dependence of magnetic amplitudes for NdCu. refer to the reciprocal lattice . Note, in case of non-orthogonal axes the convention for applied field in these tables is , and perpendicular to and .*mcphasa.hkl*- Fourier Transform of the -component of the magnetic Moments.
refer to the reciprocal lattice
. The single ion module function
*mcalc*is used to evaluate the magnetic moment on each site. Note, in case of non-orthogonal axes the convention for applied field and the magnetic moment component in these tables is , and perpendicular to and . *mcphasb.hkl*- Fourier Transform of the -component of the magnetic Moments. The single ion module function
*mcalc*is used to evaluate the magnetic moment on each site. refer to the reciprocal lattice . Note, in case of non-orthogonal axes the convention for applied field and the magnetic moment component in these tables is , and perpendicular to and . *mcphasc.hkl*- Fourier Transform of the -component of the magnetic Moments. The single ion module function
*mcalc*is used to evaluate the magnetic moment on each site. refer to the reciprocal lattice . Note, in case of non-orthogonal axes the convention for applied field and the magnetic moment component in these tables is , and perpendicular to and .

*Exercises:*

- Look at the output files of
*McPhase*in the directory*examples/ndcu2b_new/results*. At which magnetic field the ferromagnetically aligned state is achieved (at 2 K)? - What is the propagation vector in the different antiferromagnetic phases at 2 K ?

Martin Rotter 2017-01-10