Effects of aluminum on hydrogen solubility and diffusion in deformed Fe-Mn alloys

We discuss hydrogen diffusion and solubility in aluminum alloyed Fe-Mn alloys. The systems of interest are subjected to tetragonal and isotropic deformations. Based on ab initio modelling, we calculate solution energies, then employ Oriani's theory which reflects the influence of Al alloying via trap site diffusion. This local equilibrium model is complemented by qualitative considerations of Einstein diffusion. Therefore, we apply the climbing image nudged elastic band method to compute the minimum energy paths and energy barriers for hydrogen diffusion. Both for diffusivity and solubility of hydrogen, we find that the influence of the substitutional Al atom has both local chemical and nonlocal volumetric contributions.Comment: 9 page

The influence of short range forces on melting along grain boundaries

We investigate a model which couples diffusional melting and nanoscale structural forces via a combined nano-mesoscale description. Specifically, we obtain analytic and numerical solutions for melting processes at grain boundaries influenced by structural disjoining forces in the experimentally relevant regime of small deviations from the melting temperature. Though spatially limited to the close vicinity of the tip of the propagating melt finger, the influence of the disjoining forces is remarkable and leads to a strong modification of the penetration velocity. The problem is represented in terms of a sharp interface model to capture the wide range of relevant length scales, predicting the growth velocity and the length scale describing the pattern, depending on temperature, grain boundary energy, strength and length scale of the exponential decay of the disjoining potential. Close to equilibrium the short-range effects near the triple junctions can be expressed through a contact angle renormalisation in a mesoscale formulation. For higher driving forces strong deviations are found, leading to a significantly higher melting velocity than predicted from a purely mesoscopic description.Comment: 10 page

Elastic and plastic effects on heterogeneous nucleation and nanowire formation

We investigate theoretically the effects of elastic and plastic deformations on heterogeneous nucleation and nanowire formation. In the first case, the influence of the confinement of the critical nucleus between two parallel misfitting substrates is investigated using scaling arguments. We present phase diagrams giving the nature of the nucleation regime as a function of the driving force and the degree of confinement. We complement this analytical study by amplitude equations simulations. In the second case, the influence of a screw dislocation inside a nanowire on the development of the morphological surface stability of the wire, related to the Rayleigh-Plateau instability, is examined. Here the screw dislocation provokes a torsion of the wire known as Eshelby twist. Numerical calculations using the finite element method and the amplitude equations are performed to support analytical investigations. It is shown that the screw dislocation promotes the Rayleigh-Plateau instability.Comment: 16 page

Pattern formation during diffusion limited transformations in solids

We develop a description of diffusion limited growth in solid-solid transformations, which are strongly influenced by elastic effects. Density differences and structural transformations provoke stresses at interfaces, which affect the phase equilibrium conditions. We formulate equations for the interface kinetics similar to dendritic growth and study the growth of a stable phase from a metastable solid in both a channel geometry and in free space. We perform sharp interface calculations based on Green's function methods and phase field simulations, supplemented by analytical investigations. For pure dilatational transformations we find a single growing finger with symmetry breaking at higher driving forces, whereas for shear transformations the emergence of twin structures can be favorable. We predict the steady state shapes and propagation velocities, which can be higher than in conventional dendritic growth.Comment: submitted to Philosophical Magazin

Non-linear elastic effects in phase field crystal and amplitude equations: Comparison to ab initio simulations of bcc metals and graphene

We investigate non-linear elastic deformations in the phase field crystal model and derived amplitude equations formulations. Two sources of non-linearity are found, one of them based on geometric non-linearity expressed through a finite strain tensor. It reflects the Eulerian structure of the continuum models and correctly describes the strain dependence of the stiffness. In general, the relevant strain tensor is related to the left Cauchy-Green deformation tensor. In isotropic one- and two-dimensional situations the elastic energy can be expressed equivalently through the right deformation tensor. The predicted isotropic low temperature non-linear elastic effects are directly related to the Birch-Murnaghan equation of state with bulk modulus derivative $K'=4$ for bcc. A two-dimensional generalization suggests $K'_{2D}=5$. These predictions are in agreement with ab initio results for large strain bulk deformations of various bcc elements and graphene. Physical non-linearity arises if the strain dependence of the density wave amplitudes is taken into account and leads to elastic weakening. For anisotropic deformations the magnitudes of the amplitudes depend on their relative orientation to the applied strain.Comment: 16 page

Phase field modelling of grain boundary premelting using obstacle potentials

We investigate the multi-order parameter phase field model of Steinbach and Pezzolla [I. Steinbach, F. Pezzolla, A generalized field method for multiphase transformations using interface fields, Physica D 134 (1999) 385-393] concerning its ability to describe grain boundary premelting. For a single order parameter situation solid-melt interfaces are always attractive, which allows to have (unstable) equilibrium solid-melt-solid coexistence above the bulk melting point. The temperature dependent melt layer thickness and the disjoining potential, which describe the interface interaction, are affected by the choice of the thermal coupling function and the measure to define the amount of the liquid phase. Due to the strictly finite interface thickness also the interaction range is finite. For a multi-order parameter model we find either purely attractive or purely repulsive finite-ranged interactions. The premelting transition is then directly linked to the ratio of the grain boundary and solid-melt interfacial energy.Comment: 12 page

Phase Field Modeling of Fracture and Stress Induced Phase Transitions

We present a continuum theory to describe elastically induced phase transitions between coherent solid phases. In the limit of vanishing elastic constants in one of the phases, the model can be used to describe fracture on the basis of the late stage of the Asaro-Tiller-Grinfeld instability. Starting from a sharp interface formulation we derive the elastic equations and the dissipative interface kinetics. We develop a phase field model to simulate these processes numerically; in the sharp interface limit, it reproduces the desired equations of motion and boundary conditions. We perform large scale simulations of fracture processes to eliminate finite-size effects and compare the results to a recently developed sharp interface method. Details of the numerical simulations are explained, and the generalization to multiphase simulations is presented