Search

KSiH 3 exhibits 4.1 wt% experimental hydrogen storage capacity and shows reversibility under moderate conditions, which provides fresh impetus to the search for other complex hydrides in the K-Si-H system. Here, we reproduce the stable Fm 3̄ m phase of K 2 SiH 6 and uncover two denser phases, space groups P 3̄ m 1 and P 6 3 mc at ambient pressure, by means of first-principles structure searches. We note that P 3̄ m 1-K 2 SiH 6 has a high hydrogen content of 5.4 wt% and a volumetric density of 88.3 g L -1 . Further calculations suggest a favorable dehydrogenation temperature T des of -20.1/55.8 °C with decomposition into KSi + K + H 2 . The higher hydrogen density and appropriate dehydrogenation temperature indicate that K 2 SiH 6 is a promising hydrogen storage material, and our results provide helpful and clear guidance for further experimental studies. We found three further potential hydrogen storage materials stable at high pressure: K 2 SiH 8 , KSiH 7 and KSiH 8 . These results suggest the need for further investigations into hydrogen storage materials among such ternary hydrides at high pressure.

Nuclear Magnetic Resonance (NMR) is one of the most powerful experimental techniques to characterize the structure of molecules and confined liquids. Nevertheless, the complexity of the systems under investigation usually requires complementary computational studies to interpret the NMR results. In this work we focus on polycyclic aromatic hydrocarbons (PAHs), an important class of organic molecules which have been commonly used as simple analogues for the spectroscopic properties of more complex systems, such as porous disordered carbons. We use Density Functional Theory (DFT) to calculate 13C chemical shifts and Nucleus Independent Chemical Shifts (NICS) for 34 PAHs. The results show a clear molecular size dependence of the two quantities, as well as the convergence of the 13C NMR shifts towards the values observed for graphene. We then present two computationally cheap models for the prediction of NICS in simple PAHs. We show that while a simple dipolar model fails to produce accurate values, a perturbative tight-binding approach can be successfully applied for the prediction of NICS in this series of molecules, including some non-planar ones containing 5- and 7-membered rings. This model, one to two orders of magnitude faster than DFT calculations, is very promising and can be further refined in order to study more complex systems.

Predicting phase stabilities of crystal polymorphs is central to computational materials science and chemistry. Such predictions are challenging because they first require searching for potential energy minima and then performing arduous free-energy calculations to account for entropic effects at finite temperatures. Here, we develop a framework that facilitates such predictions by exploiting all the information obtained from random searches of crystal structures. This framework combines automated clustering, classification and visualisation of crystal structures with machine-learning estimation of their enthalpy and entropy. We demonstrate the framework on the technologically important system of TiO 2 , which has many polymorphs, without relying on prior knowledge of known phases. We find a number of new phases and predict the phase diagram and metastabilities of crystal polymorphs at 1600 K, benchmarking the results against full free-energy calculations.

We present computational discoveries of new structural phases of the B2O compound exhibiting novel bonding networks and electronic states at ambient and elevated pressures. Our advanced crystal structure searches in conjunction with density functional theory calculations have identified an orthorhombic phase of B2O that is energetically stable at ambient pressure and contains an intriguing bonding network of icosahedral B12 clusters bridged by oxygen atoms. As pressure increases above 1.9 GPa, a structural transformation takes the orthorhombic B2O into a pseudo-layered trigonal phase. We have performed extensive studies to investigate the evolution of chemical bonds and electronic states associated with the B12 icosahedral unit in the orthorhombic phase and the covalent B-O bonds in the trigonal phase. We have also examined the nature of the charge carriers and their coupling to the lattice vibrations in the newly identified B2O crystals. Interestingly, our results indicate that both B2O phases become superconducting at low temperatures, with transition temperatures of 6.4 K and 5.9 K, respectively, in the ambient and high-pressure phase. The present findings establish new B2O phases and characterize their structural and electronic properties, which offer insights and guidance for exploration toward further fundamental understanding and potential synthesis and application.

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.

We present a systematic study of the stability of nineteen different periodic structures using the finite range Lennard-Jones potential model discussing the effects of pressure, potential truncation, cutoff distance and Lennard-Jones exponents. The structures considered are the hexagonal close packed (hcp), face centred cubic (fcc) and seventeen other polytype stacking sequences, such as dhcp and 9R. We found that at certain pressure and cutoff distance values, neither fcc nor hcp is the ground state structure as previously documented, but different polytypic sequences. This behaviour shows a strong dependence on the way the tail of the potential is truncated.

Solid solution perovskite Pb(Zr 1-x Ti x )O 3 (PZT) is an industrially important material. Despite the long history of experimental and theoretical studies, the structure of this material is still under intensive discussion. In this work, we have applied structure searching coupled with density functional theory methods to provide a multiphase description of this material at x = 0.4. We demonstrate that the permutational freedom of B-site cations leads to the stabilisation of a variety of local phases reflecting a relatively flat energy landscape of PZT. Using a set of predicted local phases we reproduce the experimental pair distribution function (PDF) profile with high accuracy. We introduce a complex multiphase picture of the structure of PZT and show that additional monoclinic and rhombohedral phases account for a better description of the experimental PDF profile. We propose that such a multiphase picture reflects the entropy reached in the sample during the preparation process.

The structural chemistry of materials containing low levels of nonstoichiometric hydrogen is difficult to determine, and producing structural models is challenging where hydrogen has no fixed crystallographic site. Here we demonstrate a computational approach employing ab initio random structure searching (AIRSS) to generate a series of candidate structures for hydrous wadsleyite (β-Mg2SiO4 with 1.6 wt% H2O), a high-pressure mineral proposed as a repository for water in the Earth's transition zone. Aligning with previous experimental work, we solely consider models with Mg3 (over Mg1, Mg2 or Si) vacancies. We adapt the AIRSS method by starting with anhydrous wadsleyite, removing a single Mg(2+) and randomly placing two H(+) in a unit cell model, generating 819 candidate structures. 103 geometries were then subjected to more accurate optimisation under periodic DFT. Using this approach, we find the most favourable hydration mechanism involves protonation of two O1 sites around the Mg3 vacancy. The formation of silanol groups on O3 or O4 sites (with loss of stable O1-H hydroxyls) coincides with an increase in total enthalpy. Importantly, the approach we employ allows observables such as NMR parameters to be computed for each structure. We consider hydrous wadsleyite (∼1.6 wt%) to be dominated by protonated O1 sites, with O3/O4-H silanol groups present as defects, a model that maps well onto experimental studies at higher levels of hydration (J. M. Griffin et al., Chem. Sci., 2013, 4, 1523). The AIRSS approach adopted herein provides the crucial link between atomic-scale structure and experimental studies.

Selenium-rich Ge-Te-Se glasses have been synthesized along the GeSe4-GeTe4 pseudo-composition line and acquired by (77)Se Hahn echo magic-angle spinning NMR. The comparison with the GeSe4 spectrum shows a drastic modification of the typical double-resonance lineshape even at low Te concentrations (<10%). In order to rationalize this feature and to understand the effect of Te on the structure of our glasses, first-principles molecular dynamics simulations and gauge including projector augmented wave NMR parameter calculations have been performed. The distribution of the tellurium atoms in the selenium phase was shown to be mainly responsible for the (77)Se lineshape changes. Another possible factor related to the perturbation of the δiso value due to Te proximity appears to be much more limited in the bulk, while the results obtained using molecular models suggest shifts of several hundreds of ppm.

Structures of calcium peroxide (CaO2) are investigated in the pressure range 0-200 GPa using the ab initio random structure searching (AIRSS) method and density functional theory (DFT) calculations. At 0 GPa, there are several CaO2 structures very close in enthalpy, with the ground-state structure dependent on the choice of exchange-correlation functional. Further stable structures for CaO2 with C2/c, I4/mcm and P21/c symmetries emerge at pressures below 40 GPa. These phases are thermodynamically stable against decomposition into CaO and O2. The stability of CaO2 with respect to decomposition increases with pressure, with peak stability occurring at the CaO B1-B2 phase transition at 65 GPa. Phonon calculations using the quasiharmonic approximation show that CaO2 is a stable oxide of calcium at mantle temperatures and pressures, highlighting a possible role for CaO2 in planetary geochemistry. We sketch the phase diagram for CaO2, and find at least five new stable phases in the pressure-temperature ranges 0 ≤ P ≤ 60 GPa, 0 ≤ T ≤ 600 K, including two new candidates for the zero-pressure ground state structure.

Solid-state (77)Se NMR measurements, first-principles molecular dynamics and DFT calculations of NMR parameters were performed to gain insight into the structure of selenium-rich GexSe(1-x) glasses. We recorded the fully-relaxed NMR spectra on natural abundance and 100% isotopically enriched GeSe4 samples, which led us to reconsider the level of structural heterogeneity in this material. In this paper, we propose an alternative procedure to initialise molecular dynamics runs for the chalcogenide glasses. The (77)Se NMR spectra calculated on the basis of the structural models deduced from these simulations are consistent with the experimental spectrum.

Periodic density functional theory (DFT) calculations have recently emerged as a popular tool for assigning solid-state nuclear magnetic resonance (NMR) spectra. However, in order for the calculations to yield accurate results, accurate structural models are also required. In many cases the structural model (often derived from crystallographic diffraction) must be optimised (i.e., to an energy minimum) using DFT prior to the calculation of NMR parameters. However, DFT does not reproduce weak long-range "dispersion" interactions well, and optimisation using some functionals can expand the crystallographic unit cell, particularly when dispersion interactions are important in defining the structure. Recently, dispersion-corrected DFT (DFT-D) has been extended to periodic calculations, to compensate for these missing interactions. Here, we investigate whether dispersion corrections are important for aluminophosphate zeolites (AlPOs) by comparing the structures optimised by DFT and DFT-D (using the PBE functional). For as-made AlPOs (containing cationic structure-directing agents (SDAs) and framework-bound anions) dispersion interactions appear to be important, with significant changes between the DFT and DFT-D unit cells. However, for calcined AlPOs, where the SDA-anion pairs are removed, dispersion interactions appear much less important, and the DFT and DFT-D unit cells are similar. We show that, while the different optimisation strategies yield similar calculated NMR parameters (providing that the atomic positions are optimised), the DFT-D optimisations provide structures in better agreement with the experimental diffraction measurements. Therefore, it appears that DFT-D calculations can, and should, be used for the optimisation of calcined and as-made AlPOs, in order to provide the closest agreement with all experimental measurements.

(77)Se NMR parameters for three prototypical crystalline compounds (As2Se3, As4Se4 and As4Se3) have been determined from solid-state NMR spectra in the framework of an investigation concerning AsxSe(1-x) glass structure understanding. Density functional NMR calculations using the gauge including projector augmented wave methodology have been performed on X-ray and optimized crystal structures for a set of selenium-based crystals. These theoretical results have been combined with the experimental data in order to achieve a precise assignment of the spectral lines. This work and the high sensitivity of solid-state NMR to local order show that the structure of As4Se3 should be reinvestigated using state-of-the-art diffraction techniques. Calculations performed on several molecules derived from the crystal structures have demonstrated the limited effect of interlayer or intermolecular interactions on the isotropic chemical shifts. These interactions are therefore not responsible for the unexpected large chemical shift difference observed between these three systems that could mostly be attributed to the presence of short rings.

Ab initio DFT calculations of (93)Nb NMR parameters using the NMR-CASTEP code were performed for a series of over fifty individual niobates, and a good agreement has been found with experimental NMR parameters. New experimental and calculated (93)Nb NMR data were obtained for several compounds, AlNbO4, VNb9O25, K8Nb6O19 and Cs3NbO8, which are of particular interest for catalysis. Several interesting trends have been identified between (93)Nb NMR parameters and the specifics of niobium site environments in niobates. These trends may serve as useful guidelines in interpreting (93)Nb NMR spectra of complex niobium oxide systems, including amorphous samples and niobium-based multicomponent heterogeneous catalysts. Potential applications of (93)Nb NMR to study solid polyoxoniobates are discussed.

Solid-state (95)Mo nuclear magnetic resonance (NMR) properties of molybdenum hexacarbonyl have been computed using density functional theory (DFT) based methods. Both quadrupolar coupling and chemical shift parameters were evaluated and compared with parameters of high precision determined using single-crystal (95)Mo NMR experiments. Within a molecular approach, the effects of major computational parameters, i.e. basis set, exchange-correlation functional, treatment of relativity, have been evaluated. Except for the isotropic parameter of both chemical shift and chemical shielding, computed NMR parameters are more sensitive to geometrical variations than computational details. Relativistic effects do not play a crucial part in the calculations of such parameters for the 4d transition metal, in particular isotropic chemical shift. Periodic DFT calculations were tackled to measure the influence of neighbouring molecules on the crystal structure. These effects have to be taken into account to compute accurate solid-state (95)Mo NMR parameters even for such an inorganic molecular compound., (This journal is © the Owner Societies 2011)

The local structure and cation disorder in Y(2)Ti(2-x)Sn(x)O(7) pyrochlores, materials proposed for the encapsulation of lanthanide- and actinide-bearing radioactive waste, is studied using (119)Sn (I = 1/2) NMR spectroscopy. NMR provides an excellent probe of disorder, as it is sensitive to the atomic scale environment without the need for any long-range periodicity. However, the complex and overlapping spectral resonances that often result can be difficult to interpret. Here, we demonstrate how (119)Sn DFT calculations can be used to aid the spectral interpretation and assignment, confirming that Sn occupies only the six-coordinate pyrochlore B site, and that the Sn chemical shift is sensitive to the number of Sn/Ti on the neighbouring B sites. Although distinct resonances are resolved experimentally when the Ti content is low, there is significant spectral overlap for Ti-rich compositions. We establish that this is a result of two competing contributions to the Sn chemical shift; an upfield shift resulting from the incorporation of the more polarizing Ti(4+) cation onto the neighbouring B sites, and a concomitant downfield shift arising from the decrease in unit cell size. Despite the considerably easier spectral acquisition, the lower resolution in the (119)Sn spectra hinders the extraction of the detailed structural information previously obtained using (89)Y NMR. However, the spectra we obtain are consistent with a random distribution of Sn/Ti on the pyrochlore B sites. Finally, we consider whether an equilibrium structure has been achieved by investigating materials that have been annealed for different durations.

A disaccharide is a challenging case for high-resolution (1)H solid-state NMR because of the 24 distinct protons (14 aliphatic and 10 OH) having (1)H chemical shifts that all fall within a narrow range of approximately 3 to 7 ppm. High-resolution (1)H (500 MHz) double-quantum (DQ) combined rotation and multiple pulse sequence (CRAMPS) solid-state NMR spectra of beta-maltose monohydrate are presented. (1)H-(1)H DQ-SQ CRAMPS spectra are presented together with (1)H (DQ)-(13)C correlation spectra obtained with a new pulse sequence that correlates a high-resolution (1)H DQ dimension with a (13)C single quantum (SQ) dimension using the refocused INEPT pulse-sequence element to transfer magnetization via one-bond (13)C-(1)H J couplings. Compared to the observation of only a single broad peak in a (1)H DQ spectrum recorded at 30 kHz magic-angle spinning (MAS), the use of DUMBO (1)H homonuclear decoupling in the (1)H DQ CRAMPS experiment allows the resolution of distinct DQ correlation peaks which, in combination with first-principles chemical shift calculations based on the GIPAW (Gauge Including Projector Augmented Waves) plane-wave pseudopotential approach, enables the assignment of the (1)H resonances to the 24 distinct protons. We believe this to be the first experimental solid-state NMR determination of the hydroxyl OH (1)H chemical shifts for a simple sugar. Variable-temperature (1)H-(1)H DQ CRAMPS spectra reveal small increases in the (1)H chemical shifts of the OH resonances upon decreasing the temperature from 348 K to 248 K.

A protocol for the structure determination of powdered solids at natural abundance by NMR is presented and illustrated for the case of the small drug molecule thymol. The procedure uses proton spin-diffusion data from two-dimensional NMR experiments in combination with periodic DFT refinements incorporating (1)H and (13)C NMR chemical shifts. For thymol, the method yields a crystal structure for the powdered sample, which differs by an atomic root-mean-square-deviation (all atoms except methyl group protons) of only 0.07 A from the single crystal X-ray diffraction structure with DFT-optimized proton positions.

The high-resolution 27Al and 31P NMR spectra of two as-synthesized forms of the microporous aluminophosphate AlPO-14 and the corresponding calcined-dehydrated form were assigned using both "first-principles" calculations of NMR parameters (GIPAW, as implemented in NMR-CASTEP) and a 27Al-31P heteronuclear correlation NMR experiment (MQ-J-HETCOR) that exploits 27Al multiple-quantum coherences and J couplings to identify Al-O-P linkages. NMR parameters calculated from published AlPO-14 crystal structures, which are derived from powder X-ray diffraction (XRD) data, are in poor agreement with experiment and it was necessary to optimize the structure geometry using energy minimization before satisfactory agreement was obtained. Comparison of simulated powder XRD patterns from the experimental and the energy-minimized structures shows that the changes in relative atomic positions in the optimized structure are relatively small and yield only minor adjustments in the Bragg peak intensities. These results indicate that a combination of NMR spectroscopy and first-principles calculation of NMR parameters may soon be considered a generally useful step in the refinement of the structures of microporous materials derived from powder diffraction data.

The nuclear magnetic resonance (NMR) shielding and electric field gradient (EFG) tensors of three polymorphs of Mg(2)SiO(4), forsterite (alpha-Mg(2)SiO(4)), wadsleyite (beta-Mg(2)SiO(4)) and ringwoodite (gamma-Mg(2)SiO(4)), have been calculated using a density functional theory (DFT) approach with a planewave basis set and pseudopotential approximation. These Mg(2)SiO(4) polymorphs are the principal components of the Earth down to depths of 660 km and have been proposed as the hosts of water in the Earth's upper mantle and transition zone. A comparison of our calculations with single-crystal spectroscopic data in the literature for the alpha-polymorph, forsterite, shows that both the magnitude and orientation of the shielding and EFG tensors for O and Si can be obtained with sufficient accuracy to distinguish subtle differences in atomic positions between published structures. We compare calculated (17)O MAS NMR quadrupolar powder lineshapes directly with experimental lineshapes and show that we are able to reproduce them within the precision with which the NMR parameters may be determined from multi-parameter fitting. The relatively small amounts of sample available for the beta- and gamma-polymorphs, arising from the high pressures required for synthesis, has hindered the extraction of NMR parameters in previous work. The application of DFT calculations to these high-pressure polymorphs confirms previous spectral assignments, and provides deeper insight into the empirical correlations and observations reported in the literature. These first-principles methods are highly promising for the determination of local bonding in more complex materials, such as the hydrated forms of Mg(2)SiO(4), by aiding analysis of their multinuclear NMR spectra.

The (13)C CPMAS spectrum is presented for the polymorph of oxybuprocaine hydrochloride which is stable at room temperature, i.e. Mod. II degrees . It shows crystallographic splittings arising from the fact that there are two molecules, with substantially different conformations, in the asymmetric unit. An INADEQUATE two-dimensional experiment was used to link signals for the same independent molecule. The chemical shifts are discussed in relation to the crystal structure. Of the four ethyl groups attached to NH(+) nitrogens, one gives rise to unusually low chemical shifts, very different from those of the other three ethyl groups. This is attributed empirically to gamma-gauche conformational effects, as is confirmed by shielding computations. These considerations allow (13)C signals to be assigned to specific carbons in the two crystallographically inequivalent molecules in the crystal structure. Indeed, information about the conformations is inherent in the NMR spectrum, which thus provides data of crystallographic significance. A (13)C/(1)H HETCOR experiment enabled resolution to be obtained in the (1)H dimension and allowed (1)H and (13)C signals for the same independent molecule to be linked.

We show how powder samples at natural isotopic abundance can be assigned to crystal structures by using high-resolution proton and carbon-13 solid-state NMR spectra in combination with first principles calculations. Homonuclear proton double-quantum spectra in combination with through-bond proton-carbon HSQC spectra are used to assign the NMR spectra. We then show that the proton chemical shifts can be included in the process of assigning the spectra to a crystal structure using first principles calculations. The method is demonstrated on the K salt of penicillin G.

The distorted perovskites NaTaO(3) and NaNbO(3) have been studied using (23)Na multiple-quantum (MQ) MAS NMR. NaTaO(3) was prepared by high temperature solid state synthesis and the NMR spectra are consistent with the expected room temperature structure of the material (space group Pbnm), with a single crystallographic sodium site. Two samples of NaNbO(3) were studied. The first, a commercially available sample which was annealed at 900 degrees C, showed two crystallographic sodium sites, as expected for the room temperature structure of the material (space group Pbcm). The second sample, prepared by a low temperature hydrothermal method, showed the presence of four sodium sites, two of which match the expected room temperature structure and the second pair, another polymorph of the material (space group P21ma). This is consistent with powder X-ray diffraction data which showed weak extra peaks which can be accounted for by the presence of this second polymorph. Density functional theory (DFT) calculations support our conclusions, and aid assignment of the NMR spectra. Finally, we discuss the measured NMR parameters in relation to other studies of sodium in high coordination sites in the solid state.

A (13)C CPMAS NMR experiment at high field (11.7 T) has produced significantly improved dispersion for the alpha form of testosterone, allowing revisions and extensions to be made to the assignments. Correlations shown by an INADEQUATE two-dimensional spectrum, recorded at 16.5 T, have allowed the components of most of the doublet signals to be grouped into two sets (for the two crystallographically independent molecules). First-principles computations, employing a fully solid-state approach, have been used to obtain values for the crystallographic splittings, which are discussed in relation to the experimental values. This procedure enables assignments to the two groups to be suggested for all but one of the remaining doublet signals. It also allows the two sets of signals to be identified specifically to the two independent molecules in the crystal structure. Computations were also carried out for the beta form of testosterone (a dihydrate). The shift differences between the alpha and beta forms were compared with the experimental data, with encouraging results. Comparisons were also made between computed and experimental shielding anisotropies and asymmetries for three of the carbons of the alpha form. The methodology has a high potential for future applications, though more examples need to be evaluated before general conclusions can be drawn.

The 1H, 13C and 19F magic-angle spinning NMR spectra have been recorded for Form 1 of flurbiprofen. In the case of 19F, spinning sideband analysis has produced data for the components of the shielding tensor. The chemical shift of the hydrogen-bonded proton was found to be 14.0 ppm. Shielding parameters for all three nuclei have been calculated using Density Functional Theory (DFT) together with the Gauge Including Projector Augmented Wave (GIPAW) method which takes full allowance for the repetition inherent in crystalline structures. Such computations were made for the reported geometry, for a structure with all the atomic positions relaxed using DFT, and with only the hydrogen positions relaxed. The relationships of the computed shifts to those observed are discussed. In general, the correlations are good.