Point defects induced by doping rare earth elements RE (Nd, Er) into the magnesium oxide host were investigated by classical atomistic simulations utilising the General Utility Lattice Program (GULP). Formation and association energies were calculated for the potential defect structures. Both isolated defects and defect complexes were considered. The most energetically favourable structures of defect complexes were found for rare earth doped and Li co-doped systems. The correlation between the association energy and the structure of the defect complex was investigated. The influences of Li were revealed with respect to energy and structure. The simulation results contribute to the understanding of the point defects of doped MgO and how Li influences the doping of rare earth elements in the MgO host.

Meta-aminobenzoic acid, an important model system in the study of polymorphism and crystallization of active pharmaceutical ingredients, exist in water in both the nonionic (mABA) and zwitterionic (mABA^±) forms. However, the constituent molecules of the polymorph that crystallizes from aqueous solutions are zwitterionic. This study reports atomistic simulations of the events surrounding the early stage of crystal nucleation of meta-aminobenzoic acid from aqueous solutions. Ab initio molecular dynamics was used to simulate the hydration of mABA± and mABA, and to quantify the interaction of these molecules with the surrounding water molecules. Density functional theory calculations were conducted to determine the low-lying energy conformers of meta-aminobenzoic acid dimers and compute the Gibbs free energies in water of nonionic, (mABA)₂, zwitterionic, (mABA^±)₂ and nonionic-zwitterionic, (mABA)(mABA^±), species. Classical molecular dynamics simulations of mixed mABA–mABA^± aqueous solutions were carried out to examine the aggregation of meta-aminobenzoic acid. According to these simulations the selective crystallization of the polymorph which constituent molecules are zwitterionic is driven by the formation of zwitterionic dimers in solution, which are thermodynamically more stable than (mABA)₂ and (mABA)(mABA^±) pairs. This work represents a paradigm of the role of molecular processes during the early stages of crystal nucleation in affecting polymorph selection during crystallization from solution.

The microstructure of oil pipeline steels is polycrystalline due to the presence of a complex grain boundary (GB) network in them. The performance of them when exposed to stresses in corrosive environments is largely dictated by the interplay between the hydrogen (H) diffusion, the dislocation slip, and their interactions with the GBs. However, the slip-GB interaction in the presence of H and the subsequent cracking along the GBs remains not fully understood up to date. It is challenging to simultaneously resolve the collective dislocation motions away from the GBs together with the H diffusion near the atomically structured GBs using single-scale techniques. To address this challenge, taking bcc iron as a model material, here we present a a concurrent atomistic-continuum (CAC) computational analysis of the interaction between dislocation slip and a H-charged GB. With the large number of dislocations-mediated plastic flow away from the GB together with the atomistic structure evolution at the GB being simultaneously retained, several main findings from our simulations are: (i) the large number of dislocation-mediated pileup-induced internal stress nearby the H-charged GB can be as high as 3GPa; (ii) the more dislocations accumulated nearby the GB, the slower H diffusion ahead of the slip-GB intersection; (iii) H atoms diffuse fast behind the pileup tip, get trapped within the GB, and diffuse slow ahead of the pileup tip; (iv) the local stresses nearby the pileup tip exhibits a strong heterogeneity along the dislocation line direction. This differs from the common wisdom in many existing 2D theories/models. The buildup of such a high local stress heterogeneity leads to inhomogeneous H diffusion within the GB plane. The CAC simulation-predicted local H diffusivity, Dpileup−tip, and local stresses, σ, nearby the GB are then correlated with each other and consolidated into a mechanics model by considering the dislocation pileup as an Eshelby inclusion. These findings will provide researchers with opportunities to: (a) characterize the coupled dynamics between plasticity, H diffusion, and crack initiation underlying the hydrogen-induced cracking (HIC); (b) develop mechanism-based constitutive rules to be used in diffusion-plasticity coupling models for understanding the interplay between mechanical and mass transport in materials at the continuum level; and (c) connect the atomistic deformation physics of polycrystalline materials with their overall performance in aqueous environments, which is currently difficult to achieve in laboratory experiments.

Since its introduction, back in the late 1970s, the non-affine or slip parameter, ξ, has been routinely employed by numerous constitutive models as a constant parameter. However, the evidence seems to imply that it should be a function of polymer deformation. In the present work, we phenomenologically modify a constitutive model for the rheology of unentangled polymer melts [P. S. Stephanou et al. J. Rheol. 53, 309 (2009)] to account for a non-constant slip parameter. The revised model predictions are compared against newly accumulated rheological data for a C48 polyethylene melt obtained via direct non-equilibrium molecular dynamics simulations in shear. We find that the conformation tensor data are very well predicted; however, the predictions of the material functions are noted to deviate from the NEMD data, especially at large shear rates.