-
MICROSIM: A high performance phase-field solver based on CPU and GPU implementations
Authors:
Tanmay Dutta,
Dasari Mohan,
Saurav Shenoy,
Nasir Attar,
Abhikshek Kalokhe,
Ajay Sagar,
Swapnil Bhure,
Swaroop . S. Pradhan,
Jitendriya Praharaj,
Subham Mridha,
Anshika Kushwaha,
Vaishali Shah,
M. P. Gururajan,
V. Venkatesh Shenoi,
Gandham Phanikumar,
Saswata Bhattacharyya,
Abhik Choudhury
Abstract:
The phase-field method has become a useful tool for the simulation of classical metallurgical phase transformations as well as other phenomena related to materials science. The thermodynamic consistency that forms the basis of these formulations lends to its strong predictive capabilities and utility. However, a strong impediment to the usage of the method for typical applied problems of industria…
▽ More
The phase-field method has become a useful tool for the simulation of classical metallurgical phase transformations as well as other phenomena related to materials science. The thermodynamic consistency that forms the basis of these formulations lends to its strong predictive capabilities and utility. However, a strong impediment to the usage of the method for typical applied problems of industrial and academic relevance is the significant overhead with regard to the code development and know-how required for quantitative model formulations. In this paper, we report the development of an open-source phase-field software stack that contains generic formulations for the simulation of multi-phase and multi-component phase transformations. The solvers incorporate thermodynamic coupling that allows the realization of simulations with real alloys in scenarios directly relevant to the materials industry. Further, the solvers utilize parallelization strategies using either multiple CPUs or GPUs to provide cross-platform portability and usability on available supercomputing machines. Finally, the solver stack also contains a graphical user interface to gradually introduce the usage of the software. The user interface also provides a collection of post-processing tools that allow the estimation of useful metrics related to microstructural evolution.
△ Less
Submitted 1 April, 2024;
originally announced April 2024.
-
Extracting Rates and Activation Free Energies of Martensitic Transitions Using Nanomechanical Force Statistics: Theory, Models, and Analysis
Authors:
Arijit Maitra,
M. P. Gururajan
Abstract:
Nanomechanical responses (force-time profiles) of crystal lattices under deformation exhibit random critical jumps, reflecting the underlying structural transition processes. Despite extensive data collection, interpreting dynamic critical responses and their underlying mechanisms remains a significant challenge. This study explores a microscopic theoretical approach to analyse critical force fluc…
▽ More
Nanomechanical responses (force-time profiles) of crystal lattices under deformation exhibit random critical jumps, reflecting the underlying structural transition processes. Despite extensive data collection, interpreting dynamic critical responses and their underlying mechanisms remains a significant challenge. This study explores a microscopic theoretical approach to analyse critical force fluctuations in martensitic transitions. Extensive sampling of the critical forces was performed using nonequilibrium molecular dynamics simulations of an atomic model of single-crystalline titanium nickel. We demonstrate that a framework of nonequilibrium statistical mechanics offers a principled explanation of the relationship between strain rate and the critical force distribution as well as its mean. The martensitic transition is represented on a free energy landscape, taking into account the thermally activated evolution of atomic arrangements over a barrier during its time-dependent deformation. The framework enables consistent inference of the relevant fundamental properties (e.g., intrinsic rate, activation free energy) that define the rate process of structural transition. The study demonstrates how the statistical characterisation of nanomechanical response-stimulus patterns can offer microscopic insights into the deformation behaviours of crystalline materials.
△ Less
Submitted 25 July, 2024; v1 submitted 13 February, 2024;
originally announced February 2024.
-
Dislocation assisted coarsening of coherent precipitates: a phase field study
Authors:
Arjun Varma R,
Prita Pant,
M P Gururajan
Abstract:
Coarsening of precipitates in coherent systems is influenced by the elastic fields of the precipitates and the interfacial curvature. It is also known that if precipitates are connected by dislocations, coarsening is affected by the elastic fields of the dislocations and the pipe diffusivity. Although there is experimental evidence of accelerated coarsening in the presence of dislocations, these s…
▽ More
Coarsening of precipitates in coherent systems is influenced by the elastic fields of the precipitates and the interfacial curvature. It is also known that if precipitates are connected by dislocations, coarsening is affected by the elastic fields of the dislocations and the pipe diffusivity. Although there is experimental evidence of accelerated coarsening in the presence of dislocations, these studies do not capture the effect of the elastic fields. There exist generic theoretical models that can predict the average sizes and size distributions of coarsening precipitates considering the coherency related elastic stress fields. In this paper, we use a phase field model to study the coarsening of precipitates connected by dislocations, incorporating its elastic fields and pipe diffusivity. Specifically, we study the effects of misfit strain, elastic moduli mismatch, faster pipe mobility and the elastic fields of the dislocation on the morphology and kinetics of the coarsening precipitates. The dilatational component, associated with the edge character of the dislocations interact with the precipitates in an elastically homogeneous system. In an elastically inhomogeneous system, the deviatoric elastic fields also interact with the precipitates influencing the morphology and kinetics of coarsening precipitates. The kinetics of coarsening is different as well, with hard precipitates coarsening faster as compared to soft precipitates, when connected by dislocations. Further, we note that the modified Gibbs-Thomson equation, which was originally derived for an isolated precipitate in an infinite matrix, is also applicable for coarsening precipitates at close proximity.
△ Less
Submitted 9 January, 2024;
originally announced January 2024.
-
Theoretical issues in the accurate computation of the electron-phonon interaction contribution to the total energy
Authors:
Shilpa Paul,
M. P. Gururajan,
Amrita Bhattacharya,
T. R. S. Prasanna
Abstract:
The total energy is the most fundamental quantity in ab initio studies. To include electron-phonon interaction (EPI) contribution to the total energy, we have recast Allen's equation, for the case of semiconductors and insulators. This equivalent expression can be computed using available software, leading to more accurate total energy. We calculate the total energies and their differences for car…
▽ More
The total energy is the most fundamental quantity in ab initio studies. To include electron-phonon interaction (EPI) contribution to the total energy, we have recast Allen's equation, for the case of semiconductors and insulators. This equivalent expression can be computed using available software, leading to more accurate total energy. We calculate the total energies and their differences for carbon-diamond and carbon-hexagonal polytypes. Contrary to the current understanding, where per-atom quantities are computed from the primitive unit cell, our results, unexpectedly, show that the per-atom total energy (EPI included) depends on the unit cell size and violates the unit cell independence. For example, it differs for carbon-diamond by 1eV/atom between the primitive unit cell and supercells. This norm-violating behaviour arises from the partial Fan-Migdal (FM) self-energy term present in Allen's equation. Consequently, to determine the total energy differences between polytypes, supercells with identical number of atoms are necessary. A crucial inference of general validity is that any equation that contains a partial Fan-Migdal self-energy term violates the unit cell independence, contrary to the current understanding. Further theoretical studies are needed to establish if the total energy (EPI included) is an exception or can be reconciled with the unit cell independence.
△ Less
Submitted 11 October, 2023; v1 submitted 8 September, 2023;
originally announced September 2023.
-
Critical role of electron-phonon interactions in determining the relative stability of Boron Nitride polymorphs
Authors:
Shilpa Paul,
M. P. Gururajan,
Amrita Bhattacharya,
T. R. S. Prasanna
Abstract:
Despite several first principles studies, the relative stability of BN polymorphs remains controversial. The stable polymorph varies between the cubic (c-BN) and hexagonal (h-BN) depending on the van der Waals (vdW) dispersion approximation used. These studies are unable to explain the main experimental results, c-BN is stable, the relative stability order and the large energy difference between h…
▽ More
Despite several first principles studies, the relative stability of BN polymorphs remains controversial. The stable polymorph varies between the cubic (c-BN) and hexagonal (h-BN) depending on the van der Waals (vdW) dispersion approximation used. These studies are unable to explain the main experimental results, c-BN is stable, the relative stability order and the large energy difference between h-BN and c-BN (greater than 150 meV/formula unit). In this study, we introduce contributions from electron-phonon interactions (EPI) to the total energy of BN polymorphs. This clearly establishes c-BN is the stable polymorph irrespective of the vdW approximation. Only by including EPI contributions do the ab initio results match, for the first time, the main experimental results mentioned above. The EPI contribution to the total energy is strongly sensitive to chemical bonding (approximately twice in $sp^2$-bonded layered over $sp^3$-bonded polymorphs) and to crystal structure. The crucial role of EPI contributions is seen in $sp^2$-bonded layered BN polymorphs where it is greater than the vdW contribution. Given that h-BN is a prototype layered material, in bulk or 2D form, our results have a broader relevance, that is, including EPI correction, along with vdW approximation, is vital for the study of energetics in layered materials.
△ Less
Submitted 16 December, 2022;
originally announced December 2022.
-
Electron-phonon interaction contribution to the total energy of group IV semiconductor polymorphs: evaluation and implications
Authors:
R. Arjun Varma,
Shilpa Paul,
Anup Itale,
Pranav Pable,
Radhika Tibrewala,
Samruddhi Dodal,
Harshal Yerunkar,
Saurav Bhaumik,
Vaishali Shah,
M. P. Gururajan,
T. R. S. Prasanna
Abstract:
In density functional theory (DFT) based total energy studies, the van der Waals (vdW) and zero-point vibrational energy (ZPVE) correction terms are included to obtain energy differences between polymorphs. We propose and compute a new correction term to the total energy, due to electron-phonon interactions (EPI). We rely on Allen's general formalism, which goes beyond the Quasi-Harmonic Approxima…
▽ More
In density functional theory (DFT) based total energy studies, the van der Waals (vdW) and zero-point vibrational energy (ZPVE) correction terms are included to obtain energy differences between polymorphs. We propose and compute a new correction term to the total energy, due to electron-phonon interactions (EPI). We rely on Allen's general formalism, which goes beyond the Quasi-Harmonic Approximation (QHA), to include the free energy contributions due to quasiparticle interactions. We show that, for semiconductors and insulators, the EPI contributions to the free energies of electrons and phonons are the corresponding zero-point energy contributions. Using an approximate version of Allen's formalism in combination with the Allen-Heine theory for EPI corrections, we calculate the zero-point EPI corrections to the total energy for cubic and hexagonal polytypes of Carbon, Silicon and Silicon Carbide. The EPI corrections alter the energy differences between polytypes. In SiC polytypes, the EPI correction term is more sensitive to crystal structure than the vdW and ZPVE terms and is thus essential in determining their energy differences. It clearly establishes that the cubic SiC-3C is metastable and hexagonal SiC-4H is the stable polytype. Our results are consistent with the experimental results of Kleykamp. Our study enables the inclusion of EPI corrections as a separate term in the free energy expression. This opens the way to go beyond the QHA by including the contribution of EPI on all thermodynamic properties.
△ Less
Submitted 16 March, 2023; v1 submitted 14 April, 2022;
originally announced April 2022.
-
Dislocation assisted phase separation: a phase field study
Authors:
Arjun Varma R,
Prita Pant,
M P Gururajan
Abstract:
Defects play a key role in deciding the mechanisms and kinetics of phase transformations. In this paper, we show how dislocations influence phase separation in alloys with miscibility gap. Specifically, depending on the ratio of pipe mobility to bulk mobility, it is seen that even in a system with nominal compositions outside the spinodal limit, spinodal phase separation is possible. Surprisingly,…
▽ More
Defects play a key role in deciding the mechanisms and kinetics of phase transformations. In this paper, we show how dislocations influence phase separation in alloys with miscibility gap. Specifically, depending on the ratio of pipe mobility to bulk mobility, it is seen that even in a system with nominal compositions outside the spinodal limit, spinodal phase separation is possible. Surprisingly, phase separation through both nucleation and growth, and spinodal decomposition, is seen concurrently (for the case of intersecting dislocations). Finally, the prominent role played by dislocations in influencing the morphology of precipitates is explored. We show that these results agree qualitatively with recent experimental results in iron based systems obtained using Atom Probe Tomography (APT).
△ Less
Submitted 27 July, 2021;
originally announced July 2021.
-
Phase-field study of surface diffusion enhanced break-ups of nanowire junctions
Authors:
Abhinav Roy,
Arjun Varma R,
M P Gururajan
Abstract:
Using a phase-field model which incorporates enhanced diffusion at the nanowire surfaces, we study the effect of different parameters on the stability of intersecting nanowires. Our study shows that at the intersection of nanowires, sintering (curvature driven material flow) leads to the formation of junctions. These junctions act the initiators of nanowire break-up. The subsequent break-ups take…
▽ More
Using a phase-field model which incorporates enhanced diffusion at the nanowire surfaces, we study the effect of different parameters on the stability of intersecting nanowires. Our study shows that at the intersection of nanowires, sintering (curvature driven material flow) leads to the formation of junctions. These junctions act the initiators of nanowire break-up. The subsequent break-ups take place due to Rayleigh instability at the arms away from these junctions. Finally, at long time scales, the fragments coarsen due to the differences in sizes. The radii of the nanowires that form the junction, the difference in size of the intersecting nanowires and the angle of intersection play a dominant role in determining the kinetics of break-up while the density of intersections has little or no effect on the kinetics. We rationalise our results using maps of (i) mean curvatures (and, hence, chemical potentials), and, (ii) Interfacial Shape Distributions (ISDs) (which are based on probability densities associated with different combinations of the two principal curvatures). Finally, we use the moment of inertia tensor to characterise the (non-spherical) shapes and morphologies of (central) nanowire fragments at the junctions.
△ Less
Submitted 9 October, 2021; v1 submitted 5 July, 2021;
originally announced July 2021.
-
Stacking fault energy induced softening in the nucleation limited plasticity regime: a molecular dynamics study on Cu-Al alloys
Authors:
G Kamalakshi,
Prita Pant,
M P Gururajan
Abstract:
In general, with the addition of solutes, the yield strength of alloys is expected to increase and this phenomenon is known as solid solution strengthening. However, reports on the ``anomalous'' softening with alloying additions are not uncommon; for example, it is known that the heterogeneous nucleation of dislocations at the sites of solute atoms can lead to softening. In this study, using Molec…
▽ More
In general, with the addition of solutes, the yield strength of alloys is expected to increase and this phenomenon is known as solid solution strengthening. However, reports on the ``anomalous'' softening with alloying additions are not uncommon; for example, it is known that the heterogeneous nucleation of dislocations at the sites of solute atoms can lead to softening. In this study, using Molecular Dynamics (MD) simulations, we show anomalous softening in Cu-Al alloys deformed at 300 K; specifically, not only the yield stress but also the magnitude of stress drop at yield decrease with increasing Al content. We calculate the stress needed for the homogeneous nucleation of partial dislocation loops using a continuum model. One of the key inputs to the continuum model is the stacking fault energy (SFE). We carry out the thermodynamic integration to evaluate the free energies in crystals with and without the stacking faults and hence calculate the SFE as a function of Al content at 300 K. Using the SFE values thus obtained we show that softening is a consequence of the reduction of SFE values with the addition of Al in a system where yielding is controlled by the nucleation of dislocation loops. Specifically, the results of continuum model are in good agreement with MD simulation results of pure copper. However, in the Cu-Al alloys, the drop in yield strength can only be explained using a combination of the continuum model (which assumes homogeneous nucleation), and the reduction in heterogeneous nucleation barrier (which is a function of the ratio of the unstable and stable stacking fault energies). We show that the decrease in stress drop can be rationalised in terms of the stored energy available at yielding -- the stored energy at yielding decreases with increasing Al, and, as a result, the maximum dislocation density decreases with increasing Al content.
△ Less
Submitted 11 June, 2021;
originally announced June 2021.
-
Anisotropy in interface stress at the BCC-iron solid-melt interface: molecular dynamics and phase field crystal modelling
Authors:
Sushil Kumar,
Ming-Wei Liu,
Kuo-An Wu,
M P Gururajan
Abstract:
The interface stresses at of the solid-melt interface are, in general, anisotropic. The anisotropy in the interfacial stress can be evaluated using molecular dynamics (MD) and phase field crystal (PFC) models. In this paper, we report our results on the evaluation of the anisotropy in interface stress in a BCC solid with its melt. Specifically, we study Fe using both MD and PFC models. We show tha…
▽ More
The interface stresses at of the solid-melt interface are, in general, anisotropic. The anisotropy in the interfacial stress can be evaluated using molecular dynamics (MD) and phase field crystal (PFC) models. In this paper, we report our results on the evaluation of the anisotropy in interface stress in a BCC solid with its melt. Specifically, we study Fe using both MD and PFC models. We show that while both MD and PFC can be used for the evaluation, and the PFC and the amplitude equations based on PFC give quantitatively consistent results, the MD and PFC results are qualitatively the same but do not match quantitatively. We also find that even though the interfacial free energy is only weakly anisotropic in BCC interfacial stress anisotropy is strong. This strong anisotropy has implications for the equilibrium shapes, growth morphologies and other properties at nano-scale in these materials.
△ Less
Submitted 30 May, 2021;
originally announced May 2021.
-
Finite temperature electronic structure of Diamond and Silicon
Authors:
Vaishali Shah,
Bhavik Sanghavi,
Rahul Ramchandani,
M. P. Gururajan,
T. R. S. Prasanna
Abstract:
The electron-phonon interaction contribution to the electronic energies is included in density functional total energy calculations with ab initio pseudopotentials via the Allen formalism [Phys. Rev. B 18, 5217 (1978)] to obtain temperature dependent electronic structure of diamond and silicon. This method allows us to obtain the thermally-averaged ab initio electronic structure in a straightforwa…
▽ More
The electron-phonon interaction contribution to the electronic energies is included in density functional total energy calculations with ab initio pseudopotentials via the Allen formalism [Phys. Rev. B 18, 5217 (1978)] to obtain temperature dependent electronic structure of diamond and silicon. This method allows us to obtain the thermally-averaged ab initio electronic structure in a straightforward and computationally inexpensive way. Our investigations on the finite temperature electronic structure of diamond and silicon lead to a new criterion, temperature transferability, which is required in the ab initio pseudopotentials for temperature dependent studies. The temperature transferability of the Troullier-Martins pseudopotentials used in this work is strongly dependent on the cut-off radius and the inclusion of the unbound 3d$^0$ state. The finite temperature indirect band gaps are highly sensitive to the choice of cut-off radius used in the pseudopotentials. The finite temperature band structures and density of states show that thermal vibrations affect the electron energies throughout the valence and conduction band. We compare our results on the band gap shifts with that due to the Debye-Waller term in the Allen-Heine theory and discuss the observed differences in the zero point and high temperature band gap shifts. Although, the electron energy shifts in the highest occupied valence band and lowest unoccupied conduction band enable to obtain the changes in the indirect and direct band gaps at finite temperatures, the shifts in other electronic levels with temperature enable investigations into the finite temperature valence charge distribution in the bonding region. Thus, we demonstrate that the Allen theory provides a simple and theoretically justified formalism to obtain finite temperature valence electron charge densities that go beyond the rigid pseudo-atom approximation.
△ Less
Submitted 16 February, 2018;
originally announced February 2018.
-
Phase field modelling of precipitate morphologies in systems with tetragonal interfacial free energy anisotropy
Authors:
Arijit Roy,
M P Gururajan
Abstract:
A wide variety of morphologies arise due to the tetragonal anisotropy in interfacial free energy. In this paper, we report on a family of Extended Cahn-Hilliard (ECH) models for incorporating tetragonal anisotropy in interfacial free energy. We list the non-zero and independent parameters that are introduced in our model and list the constraints on them. For appropriate choice of these parameters,…
▽ More
A wide variety of morphologies arise due to the tetragonal anisotropy in interfacial free energy. In this paper, we report on a family of Extended Cahn-Hilliard (ECH) models for incorporating tetragonal anisotropy in interfacial free energy. We list the non-zero and independent parameters that are introduced in our model and list the constraints on them. For appropriate choice of these parameters, our model can produce a many of the morphologies seen in tetragonal systems such as di-pyramids, rods, plates and their truncated variants. We analyse these morphologies and show that they indeed are equilibrium morphologies consistent with the Wulff construction.
△ Less
Submitted 5 August, 2017;
originally announced August 2017.
-
Effect of stacking fault energy on nucleation limited plasticity in Cu-Al alloys
Authors:
G. Kamalakshi,
Prachi Limaye,
M. P. Gururajan,
Prita Pant
Abstract:
We study the effect of Stacking Fault Energy (SFE) on the deformation behaviour of copper and copper-aluminium alloys using Molecular Dynamics (MD) simulation. We find that both yield stress and the magnitude of stress drop at yield decrease with increasing Al content. This "anomalous" softening behaviour is explained on the basis of nucleation controlled yielding behaviour. Further, the decrease…
▽ More
We study the effect of Stacking Fault Energy (SFE) on the deformation behaviour of copper and copper-aluminium alloys using Molecular Dynamics (MD) simulation. We find that both yield stress and the magnitude of stress drop at yield decrease with increasing Al content. This "anomalous" softening behaviour is explained on the basis of nucleation controlled yielding behaviour. Further, the decrease in stress drop is rationalised in terms of the stored energy available at yielding - we show that this decreases with increasing Al. As a result, the maximum dislocation density is found to decrease with increasing Al content. Finally, we show that the yield stress calculated using a continuum model of homogeneous nucleation of partial loops agrees well with the yield stress seen in the simulations.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Elastic stress effects on microstructural instabilities
Authors:
M P Gururajan,
Arka Lahiri
Abstract:
This review is on the phase field modelling studies in elastic stress effects on microstructural instabilities. We will focus primarily on four elastic stress driven instabilities: (i) Spinodal phase separation; (ii) Particle splitting; (iii) Rafting; and, (iv) Asaro-Tiller-Grinfeld (ATG) instabilities. This review is organised as follows: in Section 1, we briefly describe some of the important an…
▽ More
This review is on the phase field modelling studies in elastic stress effects on microstructural instabilities. We will focus primarily on four elastic stress driven instabilities: (i) Spinodal phase separation; (ii) Particle splitting; (iii) Rafting; and, (iv) Asaro-Tiller-Grinfeld (ATG) instabilities. This review is organised as follows: in Section 1, we briefly describe some of the important and interesting experimental observations on elastic stress effects on microstructural instabilities; in Section 2, we describe, in reasonable detail, the theoretical developments in understanding the effects of elastic stress on microstructural instabilities in solids. Both Section 1 and 2 are neither comprehensive nor complete; however, they are helpful in setting the stage for discussion of (and, in giving a perspective on) phase field modelling studies that will be discussed in Section 3. In Section 4, we conclude with a summary and indication of future directions.
△ Less
Submitted 3 July, 2016;
originally announced July 2016.
-
Phase field modelling of interfacial anisotropy driven faceting of precipitates
Authors:
Arijit Roy,
E S Nani,
Arka Lahiri,
M P Gururajan
Abstract:
We use extended Cahn-Hilliard (ECH) equations to study faceted precipitate morphologies; specifically, we obtain four sided precipitates (in 2-D) and dodecahedron (in 3-D) in a system with cubic anisotropy, and, six-sided precipitates (in 2-D, in the basal plane), hexagonal dipyramids and hexagonal prisms (in 3-D) in systems with hexagonal anisotropy. Our listing of these ECH equations is fairly c…
▽ More
We use extended Cahn-Hilliard (ECH) equations to study faceted precipitate morphologies; specifically, we obtain four sided precipitates (in 2-D) and dodecahedron (in 3-D) in a system with cubic anisotropy, and, six-sided precipitates (in 2-D, in the basal plane), hexagonal dipyramids and hexagonal prisms (in 3-D) in systems with hexagonal anisotropy. Our listing of these ECH equations is fairly comprehensive and complete (upto sixth rank tensor terms of the Taylor expansion of the free energy). We also show how the parameters that enter the model are to be obtained if either the interfacial energy anisotropy or the equilibrium morphology of the precipitate is known.
△ Less
Submitted 18 April, 2016;
originally announced April 2016.
-
Growth kinetics and morphological stability of precipitates in 3-D: a phase field study
Authors:
Arijit Roy,
Mogadalai P Gururajan
Abstract:
We have studied the growth kinetics of isolated precipitates growing from a supersaturated matrix in 3-dimensions (3-D) using phase field models; we assume isotropic interfacial energy consider both constant and variable diffusivity. We report and compare our numerical growth rates with the classic analytical solutions of Zener and Frank (ZF). The numerical results deviate from the analytical ones…
▽ More
We have studied the growth kinetics of isolated precipitates growing from a supersaturated matrix in 3-dimensions (3-D) using phase field models; we assume isotropic interfacial energy consider both constant and variable diffusivity. We report and compare our numerical growth rates with the classic analytical solutions of Zener and Frank (ZF). The numerical results deviate from the analytical ones. These deviations can be understood in terms of the generalised Gibbs-Thomson effect. Specifically, due to the higher capillary contribution in 3-D (curvature is twice for a sphere compared to a circle), the precipitate growth kinetics deviates more from ZF in 3-D as compared to 2-D. In addition, the kinetic parameter associated with the normal velocity of the precipitate-matrix interface also modifies the deviation of the precipitate composition from its equilibrium value and hence its growth kinetics. In phase field models (such as the one used by us) which use a combination of Allen-Cahn and Cahn-Hilliard type equations, we show how to choose the kinetic parameters (namely, mobility and relaxation parameter) so that the kinetic coefficient (in the generalised Gibbs-Thomson equation) is made effectively zero. We also show that the kinetic parameter the precipitate-matrix interface might play a crucial role in making the precipitate undergo morphological instabilities as it grows (leading to "sea-weed"-like structures).
△ Less
Submitted 22 December, 2014;
originally announced December 2014.
-
Finite temperature external potential in crystalline solids
Authors:
T. R. S. Prasanna,
M. P. Gururajan
Abstract:
Thermal vibrations alter the external potential. Allen (Phys. Rev. B 18 (1978) 5217) proved that at finite temperatures the pseudopotential form factors are corrected by a Debye-Waller Factor (DWF). We generalize this result to the crystal potential. (The generalization to the all-electron case of the nuclear potential fails due to the breakdown of the rigid-atom approximation.) This finite temper…
▽ More
Thermal vibrations alter the external potential. Allen (Phys. Rev. B 18 (1978) 5217) proved that at finite temperatures the pseudopotential form factors are corrected by a Debye-Waller Factor (DWF). We generalize this result to the crystal potential. (The generalization to the all-electron case of the nuclear potential fails due to the breakdown of the rigid-atom approximation.) This finite temperature formalism only gives thermal-averaged properties and no dynamical information can be obtained. Hence, it is labeled the Quasi Ab Initio formalism. Analogous to the use of experimental lattice parameters in ab initio studies, experimental DWF can also be used. The justification is identical; the experimental parameters can be validated by separate ab initio studies. Our work transforms, forty years later, Kasowski's empirical study (Phys. Rev. B. 8 (1973) 1378) into the first ab initio finite temperature band structure calculation. This formalism opens the way to obtain ab initio finite temperature thermal-averaged properties from a single calculation.
△ Less
Submitted 21 May, 2015; v1 submitted 21 May, 2014;
originally announced May 2014.
-
On the incorporation of cubic and hexagonal interfacial energy anisotropy in phase field models using higher order tensor terms
Authors:
E. S. Nani,
M. P. Gururajan
Abstract:
In this paper, we show how to incorporate cubic and hexagonal anisotropies in interfacial energies in phase field models; this incorporation is achieved by including upto sixth rank tensor terms in the free energy expansion, assuming that the free energy is only a function of coarse grained composition, its gradient, curvature and aberration. We derive the number of non-zero and independent compon…
▽ More
In this paper, we show how to incorporate cubic and hexagonal anisotropies in interfacial energies in phase field models; this incorporation is achieved by including upto sixth rank tensor terms in the free energy expansion, assuming that the free energy is only a function of coarse grained composition, its gradient, curvature and aberration. We derive the number of non-zero and independent components of these tensors. Further, by demanding that the resultant interfacial energy is positive definite for inclusion of each of the tensor terms individually, we identify the constraints imposed on the independent components of these tensors. The existing results in the invariant group theory literature can be used to simplify the process of construction of some (but not all) of the higher order tensors. Finally, we derive the relevant phase field evolution equations.
△ Less
Submitted 13 April, 2014; v1 submitted 2 April, 2014;
originally announced April 2014.
-
Effect of Epitaxial Strain on Phase Separation in Thin Films
Authors:
Arka Lahiri,
T. A. Abinandanan,
M. P. Gururajan,
Saswata Bhattacharyya
Abstract:
We examine the role of an imposed epitaxial strain e in enhancing or depressing the spinodal instability of an alloy thin film. Since the alloy film starts with an imposed strain, phase separation offers a mechanism to relieve it, but only when the film is elastically inhomogeneous. With composition-dependence of elastic modulus given by y, and that of lattice parameter by η, our simulations using…
▽ More
We examine the role of an imposed epitaxial strain e in enhancing or depressing the spinodal instability of an alloy thin film. Since the alloy film starts with an imposed strain, phase separation offers a mechanism to relieve it, but only when the film is elastically inhomogeneous. With composition-dependence of elastic modulus given by y, and that of lattice parameter by η, our simulations using the Cahn-Hilliard model show (and analytical results for early stages confirm) that, for (ey/η) > 0, the imposed strain adds to the driving force for phase separation, decreases the maximally growing wave length, and expands the coherent spinodal in the phase diagram. Further, when (ey/η) > 0.372, it expands to even outside of chemical spinodal. Phase separation produces islands of elastically softer (harder) phase with (without) a favorable imposed strain. These results are in agreement with experimental results in GeSi thin films on Si and Ge substrates, as well as in InGaAs films on GaAs substrates.
△ Less
Submitted 22 October, 2013;
originally announced October 2013.
