On the entropic nucleation barrier in a martensitic transformation

Abstract

The nucleation of martensite in alloys is hindered by a free energy nucleation barrier, hence comprising contributions of the potential energy and the entropy. The leading effect is commonly attributed to the potential energy barrier due to strain fields. In this contribution, we investigate the nature of the entropic barrier by means of molecular dynamics (MD) simulations. We study a transformation process of an undercooled single crystal and examine two nucleation events observed under adiabatic conditions using vibrational mode analysis of the atomic trajectories. Our analysis shows that martensitic nucleations are indicated by transit from a state of uncorrelated into a state of correlated atomic motions. This correlation process is built up locally by a small group of atoms even before the product lattice can be recognized morphologically and it produces vibrational ‘soft’ modes along transformation paths. Phase space analyses unveil that the correlation process is characterized by narrow domains – ‘nucleation channels’ – the atomic trajectories have to pass, connecting the phase space domains of the parent and the product lattice. For a successful nucleation event, the nucleus atoms have to pass this channel collectively, which stochastically represents a rare event. Thermal fluctuations prevent finding the channel at elevated temperature and give rise for entropic stabilization of the parent phase. This ‘entropic nucleation barrier’ is reduced in the undercooled state but still effective, thus preventing the parent phase from collapsing into the product. The entropic barrier may be interpreted as the probability of a group of atoms to simultaneously pass the nucleation channel. Such group then represents a nucleus.

Keywords:

• martensitic transformation
• nucleation
• molecular dynamics simulations
• vibrational lattice analysis

Acknowledgments

The authors appreciate funding from the Deutsche Forschungsgemeinschaft (DFG) under contract 446 ISR-113/227/0-1. The simulations were calculated in parallel at the high-performance facility JUMP located at the NIC [101].