Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction.

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO2 . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO3 ⋅ 3H2 O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO3 )OH ⋅ 2H2 O]. The resolution of this controversy is central to understanding this material's thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance 13 C{1 H} (REDOR), 13 C-1 H distances are determined with precision, and the combination of 13 C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO3 ⋅ 3H2 O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO3 . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

Mechanistic basis of oxygen sensitivity in titanium.

One of the most potent examples of interstitial solute strengthening in metal alloys is the extreme sensitivity of titanium to small amounts of oxygen. Unfortunately, these small amounts of oxygen also lead to a markedly decreased ductility, which in turn drives the increased cost to purify titanium to avoid this oxygen poisoning effect. Here, we report a systematic study on the oxygen sensitivity of titanium that provides a clear mechanistic view of how oxygen impurities affect the mechanical properties of titanium. The increased slip planarity of Ti-O alloys is caused by an interstitial shuffling mechanism, which is sensitive to temperature, strain rate, and oxygen content and leads to the subsequent alteration of deformation twinning behavior. The insights from our experimental and computational work provide a rationale for the design of titanium alloys with increased tolerance to variations in interstitial content, with notable implications for more widespread use of titanium alloys.

NMR Crystallography: Evaluation of Hydrogen Positions in Hydromagnesite by 13 C{1 H} REDOR Solid-State NMR and Density Functional Theory Calculation of Chemical Shielding Tensors.

Solid-state NMR measurements coupled with density functional theory (DFT) calculations demonstrate how hydrogen positions can be refined in a crystalline system. The precision afforded by rotational-echo double-resonance (REDOR) NMR to interrogate 13 C-1 H distances is exploited along with DFT determinations of the 13 C tensor of carbonates (CO3 2- ). Nearby 1 H nuclei perturb the axial symmetry of the carbonate sites in the hydrated carbonate mineral, hydromagnesite [4 MgCO3 ⋅Mg(OH)2 ⋅4 H2 O]. A match between the calculated structure and solid-state NMR was found by testing multiple semi-local and dispersion-corrected DFT functionals and applying them to optimize atom positions, starting from X-ray diffraction (XRD)-determined atomic coordinates. This was validated by comparing calculated to experimental 13 C{1 H} REDOR and 13 C chemical shift anisotropy (CSA) tensor values. The results show that the combination of solid-state NMR, XRD, and DFT can improve structure refinement for hydrated materials.

Electric potential calculation in molecular simulation of electric double layer capacitors.

For the molecular simulation of electric double layer capacitors (EDLCs), a number of methods have been proposed and implemented to determine the one-dimensional electric potential profile between the two electrodes at a fixed potential difference. In this work, we compare several of these methods for a model LiClO4-acetonitrile/graphite EDLC simulated using both the traditional fixed-charged method (FCM), in which a fixed charge is assigned a priori to the electrode atoms, or the recently developed constant potential method (CPM) (2007 J. Chem. Phys. 126 084704), where the electrode charges are allowed to fluctuate to keep the potential fixed. Based on an analysis of the full three-dimensional electric potential field, we suggest a method for determining the averaged one-dimensional electric potential profile that can be applied to both the FCM and CPM simulations. Compared to traditional methods based on numerically solving the one-dimensional Poisson's equation, this method yields better accuracy and no supplemental assumptions.

Evaluation of the constant potential method in simulating electric double-layer capacitors.

A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations in the electrolyte solution. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [S. K. Reed et al., J. Chem. Phys. 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO4-acetonitrile/graphite EDLC. At low potential difference (ΔΨ â©½ 2 V), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher ΔΨ. At ΔΨ â©¾ 4 V, the CPM ion density profiles show significant enhancement (over FCM) of "inner-sphere adsorbed" Li(+) ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li(+) ions to the electrode surface.

Structural phase transformations in metallic grain boundaries.

Structural transformations at interfaces are of profound fundamental interest as complex examples of phase transitions in low-dimensional systems. Despite decades of extensive research, no compelling evidence exists for structural transformations in high-angle grain boundaries in elemental systems. Here we show that the critical impediment to observations of such phase transformations in atomistic modelling has been rooted in inadequate simulation methodology. The proposed new methodology allows variations in atomic density inside the grain boundary and reveals multiple grain boundary phases with different atomic structures. Reversible first-order transformations between such phases are observed by varying temperature or injecting point defects into the boundary region. Owing to the presence of multiple metastable phases, grain boundaries can absorb significant amounts of point defects created inside the material by processes such as irradiation. We propose a novel mechanism of radiation damage healing in metals, which may guide further improvements in radiation resistance of metallic materials through grain boundary engineering.

Dislocation-pairing transitions in hot grain boundaries.

We report the finding of a novel grain-boundary structural phase transition in both molecular-dynamics and phase-field-crystal simulations of classical models of bcc Fe. This transition is characterized by pairing of individual dislocations with mixed screw and edge components. We demonstrate that this type of transition is driven by a combination of factors including elastic softening, core interaction, and core disordering. At high homologous temperatures the occurrence of this transition is shown to prevent premelting at misorientation angles where it would otherwise be expected.

Structural disjoining potential for grain-boundary premelting and grain coalescence from molecular-dynamics simulations.

We describe a molecular-dynamics framework for the direct calculation of the short-ranged structural forces underlying grain-boundary premelting and grain coalescence in solidification. The method is applied in a comparative study of (i) a Sigma9115120 degrees twist and (ii) a Sigma9110{411} symmetric tilt boundary in a classical embedded-atom model of elemental Ni. Although both boundaries feature highly disordered structures near the melting point, the nature of the temperature dependence of the width of the disordered regions in these boundaries is qualitatively different. The former boundary displays behavior consistent with a logarithmically diverging premelted layer thickness as the melting temperature is approached from below, while the latter displays behavior featuring a finite grain-boundary width at the melting point. It is demonstrated that both types of behavior can be quantitatively described within a sharp-interface thermodynamic formalism involving a width-dependent interfacial free energy, referred to as the disjoining potential. The disjoining potential for boundary (i) is calculated to display a monotonic exponential dependence on width, while that of boundary (ii) features a weak attractive minimum. The results of this work are discussed in relation to recent simulation and theoretical studies of the thermodynamic forces underlying grain-boundary premelting.

Method for computing short-range forces between solid-liquid interfaces driving grain boundary premelting.

We present a molecular dynamics based method for accurately computing short-range structural forces resulting from the overlap of spatially diffuse solid-liquid interfaces at wetted grain boundaries close to the melting point. The method is based on monitoring the fluctuations of the liquid layer width at different temperatures to extract the excess interfacial free energy as a function of this width. The method is illustrated for a high-energy Sigma9 twist boundary in pure Ni. The short-range repulsion driving premelting is found to be dominant in comparison to long-range dispersion and entropic forces and consistent with previous experimental findings that nanometer-scale layer widths may be observed only very close to the melting point.

A predictive mechanism for dynamic strain ageing in aluminium-magnesium alloys.

Dynamic strain ageing (DSA) is the phenomenon in which solute atoms diffuse around dislocations and retard dislocation motion, leading to negative strain-rate sensitivity (nSRS) and thus to material instabilities during processing, an important issue in commercial metal alloys. Here, we show the mechanism of DSA and nSRS on experimental strain-rate, temperature and stress scales for Al-Mg to be single-atomic-hop motion of solutes from the compression to the tension side of a dislocation core. We derive an analytic expression for the strengthening versus strain rate and temperature that justifies widely used phenomenological forms, provides specific dependences of the parameters on material properties and is supported by atomistic kinetic Monte Carlo simulations. Using literature material properties, the predicted strengthening quantitatively agrees with the experimentally derived behaviour of Al-2.5% Mg at 300 K, and qualitatively agrees with the strain rate and temperature ranges of DSA and nSRS in Al-Mg alloys. The analyses herein show a clear path for multiscale design, from quantum to continuum mechanics, of solute strengthening in face-centred-cubic metal alloys.

Computing the mobility of grain boundaries.

As current experimental and simulation methods cannot determine the mobility of flat boundaries across the large misorientation phase space, we have developed a computational method for imposing an artificial driving force on boundaries. In a molecular dynamics simulation, this allows us to go beyond the inherent timescale restrictions of the technique and induce non-negligible motion in flat boundaries of arbitrary misorientation. For different series of symmetric boundaries, we find both expected and unexpected results. In general, mobility increases as the grain boundary plane deviates from (111), but high-coincidence and low-angle boundaries represent special cases. These results agree with and enrich experimental observations.