Example 04: Automated calculation of melting temperature#
In Example 03, we calculated the melting temperature for Cu. The same can be done in a fully automated way using Calphy.
The EAM potential we will use is : Mishin, Y., M. J. Mehl, D. A. Papaconstantopoulos, A. F. Voter, and J. D. Kress. “Structural Stability and Lattice Defects in Copper: Ab Initio , Tight-Binding, and Embedded-Atom Calculations.” Physical Review B 63, no. 22 (May 21, 2001): 224106.
The calculation block gives the input conditions at which the calculation is carried out. First of all, the
melting_temperature. This mode is special, and needs much less input information than the other common modes. For example, options such as
temperature are (generally; see below for specific cases) not needed.
The input file in this case is shown below:
element: Cu mass: 63.546 calculations: - mode: melting_temperature pressure:  repeat: [5, 5, 5] nsims: 1 md: pair_style: eam/alloy pair_coeff: "* * ../potentials/Cu01.eam.alloy Cu" timestep: 0.001 nsmall: 10000 tdamp: 0.1 pdamp: 0.1 te: 10000 ts: 15000 queue: scheduler: local cores: 4 commands: - conda activate calphy
Once the input file is set up, the calculation can be run using:
calphy -i input.yaml
A log file called
calphy.log is also be produced. The log file contains a lot of information about the calculation. Here the most important ones are discussed. The important log file contents are prefixed with
STATE. On running
grep STATE calphy.log, the following output is produced.
calphy.helpers INFO STATE: Temperature range of 1156.600000-1556.600000 K calphy.helpers INFO STATE: Tm = 1340.39 K +/- 0.00 K, Exp. Tm = 1356.60 K
The calculated melting temperature for this interatomic potential is 1340 K.
How does it work?#
A flowchart of the automated melting temperature calculation is shown below:
- How can I tune the initial guess temperature?The initial guess temperature can be tuned using the keyword
- How can I tune the width of the temperature range?The width of the temperature scan range can be tuned using the keyword
- What if the system undergoes a solid-solid phase transition before melting?The lattice that calphy automatically chooses for the solid is the ground state. For some elements, for example Ti, a solid-solid phase transition occurs before melting. HCP Ti transforms to BCC Ti and then to liquid. Therefore to use the automated method, the solid lattice has to be specified using the
- How can I calculate melting temperature at non zero pressure?The required pressure can be specified using the
pressurekeyword. Note that specifying
tguesscould help speed up the calculation.