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.

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.