A New ab-Initio Approach for NMR Chemical Shifts in Periodic Systems

NMR spectroscopy probes microstructure, dynamics and doping of metal halide perovskites

Solid-state magic-angle spinning NMR spectroscopy is a powerful technique to probe atomic-level microstructure and structural dynamics in metal halide perovskites. It can be used to measure dopant incorporation, phase segregation, halide mixing, decomposition pathways, passivation mechanisms, short-range and long-range dynamics, and other local properties. This Review describes practical aspects of recording solid-state NMR data on halide perovskites and how these afford unique insights into new compositions, dopants and passivation agents. We discuss the applicability, feasibility and limitations of ¹H, ¹³C, ¹⁵N, ¹⁴N, ¹³³Cs, ⁸⁷Rb, ³⁹K, ²⁰⁷Pb, ¹¹⁹Sn, ¹¹³Cd, ²⁰⁹Bi, ¹¹⁵In, ¹⁹F and ²H NMR in typical experimental scenarios. We highlight the pivotal complementary role of solid-state mechanosynthesis, which enables highly sensitive NMR studies by providing large quantities of high-purity materials of arbitrary complexity and of chemical shifts calculated using density functional theory. We examine the broader impact of solid-state NMR on materials research and how its evolution over seven decades has benefitted structural studies of contemporary materials such as halide perovskites. Finally, we summarize some of the open questions in perovskite optoelectronics that could be addressed using solid-state NMR. We, thereby, hope to stimulate wider use of this technique in materials and optoelectronics research.

Computation of Solid-State Vibrational Circular Dichroism in the Periodic Gauge

We introduce a new theoretical formalism to compute solid-state vibrational circular dichroism (VCD) spectra from molecular dynamics simulations. Having solved the origin-dependence problem of the periodic magnetic gauge, we present IR and VCD spectra of (1S,2S)-trans-1,2-cyclohexanediol obtained from first-principles molecular dynamics calculations and nuclear velocity perturbation theory, along with the experimental results. Because the structure model imposes periodic boundary conditions, the common origin of the rotational strength has hitherto been ill-defined and was approximated by means of averaging multiple origins. The new formalism reconnects the periodic model with the finite physical system and restores gauge freedom. It nevertheless fully accounts for nonlocal spatial couplings from the gauge transport term. We show that even for small simulation cells the rich nature of solid-state VCD spectra found in experiments can be reproduced to a very satisfactory level.

Structure Solution of Nano-Crystalline Small Molecules Using MicroED and Solid-State NMR Dipolar-Based Experiments

Three-dimensional electron diffraction crystallography (microED) can solve structures of sub-micrometer crystals, which are too small for single crystal X-ray crystallography. However, R factors for the microED-based structures are generally high because of dynamic scattering. That means R factor may not be reliable provided that kinetic analysis is used. Consequently, there remains ambiguity to locate hydrogens and to assign nuclei with close atomic numbers, like carbon, nitrogen, and oxygen. Herein, we employed microED and ssNMR dipolar-based experiments together with spin dynamics numerical simulations. The NMR dipolar-based experiments were ¹H-¹⁴N phase-modulated rotational-echo saturation-pulse double-resonance (PM-S-RESPDOR) and ¹H-¹H selective recoupling of proton (SERP) experiments. The former examined the dephasing effect of a specific ¹H resonance under multiple ¹H-¹⁴N dipolar couplings. The latter examined the selective polarization transfer between a ¹H-¹H pair. The structure was solved by microED and then validated by evaluating the agreement between experimental and calculated dipolar-based NMR results. As the measurements were performed on ¹H and ¹⁴N, the method can be employed for natural abundance samples. Furthermore, the whole validation procedure was conducted at 293 K unlike widely used chemical shift calculation at 0 K using the GIPAW method. This combined method was demonstrated on monoclinic l-histidine.

Machine learning-accelerated quantum mechanics-based atomistic simulations for industrial applications

Atomistic simulations have become an invaluable tool for industrial applications ranging from the optimization of protein-ligand interactions for drug discovery to the design of new materials for energy applications. Here we review recent advances in the use of machine learning (ML) methods for accelerated simulations based on a quantum mechanical (QM) description of the system. We show how recent progress in ML methods has dramatically extended the applicability range of conventional QM-based simulations, allowing to calculate industrially relevant properties with enhanced accuracy, at reduced computational cost, and for length and time scales that would have otherwise not been accessible. We illustrate the benefits of ML-accelerated atomistic simulations for industrial R&D processes by showcasing relevant applications from two very different areas, drug discovery (pharmaceuticals) and energy materials. Writing from the perspective of both a molecular and a materials modeling scientist, this review aims to provide a unified picture of the impact of ML-accelerated atomistic simulations on the pharmaceutical, chemical, and materials industries and gives an outlook on the exciting opportunities that could emerge in the future.

First-Principles Calculation of 1H NMR Chemical Shifts of Complex Metal Polyhydrides: The Essential Inclusion of Relativity and Dynamics

1H NMR spectroscopy has become an important technique for the characterization of transition-metal hydride complexes, whose metal-bound hydrides are often difficult to locate by X-ray diffraction. In this regard, the accurate prediction of 1H NMR chemical shifts provides a useful, but challenging, strategy to help in the interpretation of the experimental spectra. In this work, we establish a density-functional-theory protocol that includes relativistic, solvent, and dynamic effects at a high level of theory, allowing us to report an accurate and reliable interpretation of 1H NMR hydride chemical shifts of iridium polyhydride complexes. In particular, we have studied in detail the hydride chemical shifts of the [Ir6(IMe)8(CO)2H14]2+ complex in order to validate previous assignments. The computed 1H NMR chemical shifts are strongly dependent on the relativistic treatment, the choice of the DFT exchange-correlation functional, and the conformational dynamics. By combining a fully relativistic four-component electronic-structure treatment with ab initio molecular dynamics, we were able to reliably model both the terminal and bridging hydride chemical shifts and to show that two NMR hydride signals were inversely assigned in the experiment.

Computation of NMR Shielding Constants for Solids Using an Embedded Cluster Approach with DFT, Double-Hybrid DFT, and MP2

In this work, we explore the accuracy of post-Hartree–Fock (HF) methods and double-hybrid density functional theory (DFT) for the computation of solid-state NMR chemical shifts. We apply an embedded cluster approach and investigate the convergence with cluster size and embedding for a series of inorganic solids with long-range electrostatic interactions. In a systematic study, we discuss the cluster design, the embedding procedure, and basis set convergence using gauge-including atomic orbital (GIAO) NMR calculations at the DFT and MP2 levels of theory. We demonstrate that the accuracy obtained for the prediction of NMR chemical shifts, which can be achieved for molecular systems, can be carried over to solid systems. An appropriate embedded cluster approach allows one to apply methods beyond standard DFT even for systems for which long-range electrostatic effects are important. We find that an embedded cluster should include at least one sphere of explicit neighbors around the nuclei of interest, given that a sufficiently large point charge and boundary effective potential embedding is applied. Using the pcSseg-3 basis set and GIAOs for the computation of nuclear shielding constants, accuracies of 1.6 ppm for ⁷Li, 1.5 ppm for ²³Na, and 5.1 ppm for ³⁹K as well as 9.3 ppm for ¹⁹F, 6.5 ppm for ³⁵Cl, 7.4 ppm for ⁷⁹Br, and 7.5 ppm for ²⁵Mg as well as 3.8 ppm for ⁶⁷Zn can be achieved with MP2. Comparing various DFT functionals with HF and MP2, we report the superior quality of results for methods that include post-HF correlation like MP2 and double-hybrid DFT.

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

Numerical study of disorder on the orbital magnetization in two dimensions

The modern theory of orbital magnetization (OM) was developed by using Wannier function method, which has a formalism similar with the Berry phase. In this manuscript, we perform a numerical study on the fate of the OM under disorder, by using this method on the Haldane model in two dimensions, which can be tuned between a normal insulator or a Chern insulator at half filling. The effects of increasing disorder on OM for both cases are simulated. Energy renormalization shifts are observed in the weak disorder regime and topologically trivial case, which was predicted by a self-consistentT-matrix approximation. Besides this, two other phenomena can be seen. One is the localization trend of the band orbital magnetization. The other is the remarkable contribution from topological chiral states arising from nonzero Chern number or large value of integrated Berry curvature. If the fermi energy is fixed at the gap center of the clean system, there is an enhancement of |M| at the intermediate disorder, for both cases of normal and Chern insulators, which can be attributed to the disorder induced topological metal state before localization.

Rapid Structure Determination of Molecular Solids Using Chemical Shifts Directed by Unambiguous Prior Constraints

NMR-based crystallography approaches involving the combination of crystal structure prediction methods, ab initio calculated chemical shifts and solid-state NMR experiments are powerful methods for crystal structure determination of microcrystalline powders. However, currently structural information obtained from solid-state NMR is usually included only after a set of candidate crystal structures has already been independently generated, starting from a set of single-molecule conformations. Here, we show with the case of ampicillin that this can lead to failure of structure determination. We propose a crystal structure determination method that includes experimental constraints during conformer selection. In order to overcome the problem that experimental measurements on the crystalline samples are not obviously translatable to restrict the single-molecule conformational space, we propose constraints based on the analysis of absent cross-peaks in solid-state NMR correlation experiments. We show that these absences provide unambiguous structural constraints on both the crystal structure and the gas-phase conformations, and therefore can be used for unambiguous selection. The approach is parametrized on the crystal structure determination of flutamide, flufenamic acid, and cocaine, where we reduce the computational cost by around 50%. Most importantly, the method is then shown to correctly determine the crystal structure of ampicillin, which would have failed using current methods because it adopts a high-energy conformer in its crystal structure. The average positional RMSE on the NMR powder structure is ⟨r_av⟩ = 0.176 Å, which corresponds to an average equivalent displacement parameter U_eq = 0.0103 Ų.

Four-component relativistic 31P NMR calculations for trans-platinum(ii) complexes: importance of the solvent and dynamics in spectral simulations

We report a combined experimental-theoretical study on the 31P NMR chemical shift for a number of trans-platinum(ii) complexes. Validity and reliability of the 31P NMR chemical shift calculations are examined by comparing with the experimental data. A successful computational protocol for the accurate prediction of the 31P NMR chemical shifts was established for trans-[PtCl2(dma)PPh3] (dma = dimethylamine) complexes. The reliability of the computed values is shown to be critically dependent on the level of relativistic effects (two-component vs. four component), choice of density functionals, dynamical averaging, and solvation effects. Snapshots obtained from ab initio molecular dynamics simulations were used to identify those solvent molecules which show the largest interactions with the platinum complex, through inspection by using the non-covalent interaction program. We observe satisfactory accuracy from the full four-component matrix Dirac-Kohn-Sham method (mDKS) based on the Dirac-Coulomb Hamiltonian, in conjunction with the KT2 density functional, and dynamical averaging with explicit solvent molecules.

Chemical shifts in molecular solids by machine learning

Due to their strong dependence on local atonic environments, NMR chemical shifts are among the most powerful tools for strucutre elucidation of powdered solids or amorphous materials. Unfortunately, using them for structure determination depends on the ability to calculate them, which  comes at the cost of high accuracy first-principles calculations. Machine learning has recently emerged as a way to overcome the need for quantum chemical calculations, but for chemical shifts in solids it is hindered by the chemical and combinatorial space spanned by molecular solids, the strong dependency of chemical shifts on their environment, and the lack of an experimental database of shifts. We propose a machine learning method based on local environments to accurately predict chemical shifts of molecular solids and their polymorphs to within DFT accuracy. We also demonstrate that the trained model is able to determine, based on the match between experimentally measured and ML-predicted shifts, the structures of cocaine and the drug 4-[4-(2-adamantylcarbamoyl)-5-tert-butylpyrazol-1-yl]benzoic acid.

Solid-state nuclear magnetic resonance combined with quantum chemical shift predictions is limited by high computational cost. Here, the authors use machine learning based on local atomic environments to predict experimental chemical shifts in molecular solids with accuracy similar to density functional theory.

Electronic orbital response of regular extended and infinite periodic systems to magnetic fields. I. Theoretical foundations for static case

A theoretical treatment for the orbital response of an infinite, periodic system to a static, homogeneous, magnetic field is presented. It is assumed that the system of interest has an energy gap separating occupied and unoccupied orbitals and a zero Chern number. In contrast to earlier studies, we do not utilize a perturbation expansion, although we do assume the field is sufficiently weak that the occurrence of Landau levels can be ignored. The theory is developed by analyzing results for large, finite systems and also by comparing with the analogous treatment of an electrostatic field. The resulting many-electron Hamilton operator is forced to be hermitian, but hermiticity is not preserved, in general, for the subsequently derived single-particle operators that determine the electronic orbitals. However, we demonstrate that when focusing on the canonical solutions to the single-particle equations, hermiticity is preserved. The issue of gauge-origin dependence of approximate solutions is addressed. Our approach is compared with several previously proposed treatments, whereby limitations in some of the latter are identified.

Localized molecular orbitals for calculation and analysis of vibrational Raman optical activity

First calculations of vibrational Raman optical activity based on localized molecular orbitals are presented, which pave the way for novel insight into spectroscopic signatures of chiral systems.

127I NMR calculations in binary metal iodides by PBE-GGA, YS-PBE0 and mBJ exchange correlation potentials

We have calculated Nuclear Magnetic Resonance (NMR) spectroscopy for ¹²⁷I (quadrupolar nuclei I=5/2) in binary metal iodides XI (X=Li, Na, K, Rb and Cs) by using PBE- GGA, YS- PBE0 and mBJ exchange correlation potentials. The results show that the nature of bonds between Iodine and metal atoms are ionic. The main contribution in NMR spectroscopy is related to the induced current inside the atomic sphere and the remainder of the unit cell volume contributes only a few ppm. Obtained NMR shifts are compared with the NMR shielding data and the NMR shielding for metal-p band varies across the series about 221ppm. Density of states results indicate that the largest contribution in the shielding comes from the I-core electrons (1s and 4d). The NMR shielding graphs show that there are negative linear correlation with slope -1.18, -1.16 and -1.01 by PBE- GGA, YS- PBE0 and mBJ, respectively. The computed results by mBJ are in good agreement with the experimental values.

Ab Initio Quality NMR Parameters in Solid-State Materials Using a High-Dimensional Neural-Network Representation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental tools to probe the local atomic order of a wide range of solid-state compounds. However, due to the complexity of the related spectra, in particular for amorphous materials, their interpretation in terms of structural information is often challenging. These difficulties can be overcome by combining molecular dynamics simulations to generate realistic structural models with an ab initio evaluation of the corresponding chemical shift and quadrupolar coupling tensors. However, due to computational constraints, this approach is limited to relatively small system sizes which, for amorphous materials, prevents an adequate statistical sampling of the distribution of the local environments that is required to quantitatively describe the system. In this work, we present an approach to efficiently and accurately predict the NMR parameters of very large systems. This is achieved by using a high-dimensional neural-network representation of NMR parameters that are calculated using an ab initio formalism. To illustrate the potential of this approach, we applied this neural-network NMR (NN-NMR) method on the (17)O and (29)Si quadrupolar coupling and chemical shift parameters of various crystalline silica polymorphs and silica glasses. This approach is, in principal, general and has the potential to be applied to predict the NMR properties of various materials.

Effect of Co-Ordination Chemistry and Oxidation State on the (207)Pb Magnetic-Shielding Tensor: A DFT/ZORA Investigation

The magnetic shielding tensor of (207)Pb is calculated for various solids exhibiting (1) a holodirected lead(II) center containing a stereochemically inactive lone pair, (2) a hemidirected lead(II) center with a stereochemically active lone-pair, or (3) a lead(IV) center. Tensors investigated at the scalar relativistic level are compared with those calculated with the full ZORA/spin-orbit Hamiltonian. The effect of using GGA density functionals is compared to the use of hybrid density functionals.

The Calculation of NMR Chemical Shifts in Periodic Systems Based on Gauge Including Atomic Orbitals and Density Functional Theory

We present here a method that can calculate NMR shielding tensors from first principles for systems with translational invariance. Our approach is based on Kohn-Sham density functional theory and gauge-including atomic orbitals. Our scheme determines the shielding tensor as the second derivative of the total electronic energy with respect to an external magnetic field and a nuclear magnetic moment. The induced current density due to a periodic perturbation from nuclear magnetic moments is obtained through numerical differentiation, whereas the influence of the responding perturbation in terms of the external magnetic field is evaluated analytically. The method is implemented into the periodic program BAND. It employs a Bloch basis set made up of Slater-type or numeric atomic orbitals and represents the Kohn-Sham potential fully without the use of effective core potentials. Results from calculations of NMR shielding constants based on the present approach are presented for isolated molecules as well as systems with one-, two- and three-dimensional periodicity. The reported values are compared to experiment and results from calculations on cluster models.

Covalency of hydrogen bonds in liquid water can be probed by proton nuclear magnetic resonance experiments

The concept of covalency is widely used to describe the nature of intermolecular bonds, to explain their spectroscopic features and to rationalize their chemical behaviour. Unfortunately, the degree of covalency of an intermolecular bond cannot be directly measured in an experiment. Here we established a simple quantitative relationship between the calculated covalency of hydrogen bonds in liquid water and the anisotropy of the proton magnetic shielding tensor that can be measured experimentally. This relationship enabled us to quantify the degree of covalency of hydrogen bonds in liquid water using the experimentally measured anisotropy. We estimated that the amount of electron density transferred between molecules is on the order of 10  m while the stabilization energy due to this charge transfer is ∼15 kJ mol⁻¹. The physical insight into the fundamental nature of hydrogen bonding provided in this work will facilitate new studies of intermolecular bonding in a variety of molecular systems.

Covalency is a fundamental concept in chemical bonding, but experimentally it is not possible to measure the degree of covalency of a particular bond. Here, the authors report a model to link the covalency of hydrogen bonds in water with the anisotropy of the magnetic shielding tensor in the proton NMR.

NMR Shielding in Metals Using the Augmented Plane Wave Method

We present calculations of solid state NMR magnetic shielding in metals, which includes both the orbital and the complete spin response of the system in a consistent way. The latter contains an induced spin-polarization of the core states and needs an all-electron self-consistent treatment. In particular, for transition metals, the spin hyperfine field originates not only from the polarization of the valence s-electrons, but the induced magnetic moment of the d-electrons polarizes the core s-states in opposite direction. The method is based on DFT and the augmented plane wave approach as implemented in the WIEN2k code. A comparison between calculated and measured NMR shifts indicates that first-principle calculations can obtain converged results and are more reliable than initially concluded based on previous publications. Nevertheless large k-meshes (up to 2 000 000 k-points in the full Brillouin-zone) and some Fermi-broadening are necessary. Our results show that, in general, both spin and orbital components of the NMR shielding must be evaluated in order to reproduce experimental shifts, because the orbital part cancels the shift of the usually highly ionic reference compound only for simple sp-elements but not for transition metals. This development paves the way for routine NMR calculations of metallic systems.

Dynamical effects on the magnetic properties of dithiazolyl bistable materials

The magnetic properties of molecule-based magnets are commonly rationalized by considering only a single nuclear configuration of the system under study (usually an X-ray crystal structure). Here, by means of a computational study, we compare the results obtained using such a static approach with those obtained by explicitly accounting for thermal fluctuations, and uncover the serious limitations of the static perspective when dealing with magnetic crystals whose radicals undergo wide-amplitude motions. As a proof of concept, these limitations are illustrated for the magnetically bistable 1,3,5-trithia-2,4,6-triazapentalenyl (TTTA) material. For its high-temperature phase at 300 K, we show that nuclear dynamics induce large fluctuations in the magnetic exchange interactions (JAB) between spins (up to 1000% of the average value). These deviations result in a ∼20% difference between the 300 K magnetic susceptibility computed by explicitly considering the nuclear dynamics and that computed using the X-ray structure, the former being in better agreement with the experimental data. The unveiled strong coupling between JAB interactions and intermolecular vibrations reveals that considering JAB as a constant value at a given temperature (as always done in molecular magnetism) leads to a flawed description of the magnetism of TTTA. Instead, the physically relevant concept in this case is the statistical distribution of JAB values. The discovery that a single X-ray structure is not adequate enough to interpret the magnetic properties of TTTA is also expected to be decisive in other organic magnets with dominant exchange interactions propagating through labile π-π networks.

Gauge-invariant calculation of static and dynamical magnetic properties from the current density

In this work we solve two problems related to the calculation of static and dynamical magnetic properties with ab initio theories. First, we show that the dependence of the dynamical magnetic dipole moment on the reference point of the multipole expansion and on the gauge origin of the vector potential have a clear physical significance. They are due to a dynamical electric dipole moment and an electric field, respectively. Both are fully determined by the experimental setup and do not pose any fundamental problem, contrary to what is commonly assumed. Second, in the static case, any dependence on the gauge origin is an artifact of the computational method. We show that the artificial dependence on the gauge origin can be removed in an elegant way by the introduction of a sum rule that puts the diamagnetic and paramagnetic contributions on equal footing. Our approach can be applied to calculate any magnetic observable that can be derived from the current density, and can be used in combination with any ab initio theory from which it can be obtained. To illustrate our method we apply it here to time-dependent current-density-functional theory for the calculation of static and dynamical magnetizabilities of molecules.

Chemical-shift tensors of heavy nuclei in network solids: a DFT/ZORA investigation of 207Pb chemical-shift tensors using the bond-valence method

Accurate computation of ²⁰⁷Pb magnetic shielding principal components is within the reach of quantum chemistry methods by employing relativistic ZORA/DFT and cluster models adapted from the bond valence model.

NMR characterization of hydrocarbon adsorption on calcite surfaces: a first principles study

The electronic and coordination environment of minerals surfaces, as calcite, are very difficult to characterize experimentally. This is mainly due to the fact that there are relatively few spectroscopic techniques able to detect Ca(2+). Since calcite is a major constituent of sedimentary rocks in oil reservoir, a more detailed characterization of the interaction between hydrocarbon molecules and mineral surfaces is highly desirable. Here we perform a first principles study on the adsorption of hydrocarbon molecules on calcite surface (CaCO3 (101¯4)). The simulations were based on Density Functional Theory with Solid State Nuclear Magnetic Resonance (SS-NMR) calculations. The Gauge-Including Projector Augmented Wave method was used to compute mainly SS-NMR parameters for (43)Ca, (13)C, and (17)O in calcite surface. It was possible to assign the peaks in the theoretical NMR spectra for all structures studied. Besides showing different chemical shifts for atoms located on different environments (bulk and surface) for calcite, the results also display changes on the chemical shift, mainly for Ca sites, when the hydrocarbon molecules are present. Even though the interaction of the benzene molecule with the calcite surface is weak, there is a clearly distinguishable displacement of the signal of the Ca sites over which the hydrocarbon molecule is located. A similar effect is also observed for hexane adsorption. Through NMR spectroscopy, we show that aromatic and alkane hydrocarbon molecules adsorbed on carbonate surfaces can be differentiated.

Aqueous solutions: state of the art in ab initio molecular dynamics

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.

Calculation of chemical-shift tensors of heavy nuclei: a DFT/ZORA investigation of 199Hg chemical-shift tensors in solids, and the effects of cluster size and electronic-state approximations

A ZORA/DFT investigation of the NMR chemical shielding of a suite of ¹⁹⁹Hg-containing solid materials shows the importance of calculations with the molecular-cluster method.

Water and small organic molecules as probes for geometric confinement in well-ordered mesoporous carbon materials

Confinement effects on water, benzene and pyridine in mesoporous carbon materials were probed by ¹H-MAS NMR and chemical shift calculations.

Six-site polarizable model of water based on the classical Drude oscillator

A polarizable water model, SWM6, was developed and optimized for liquid phase simulations under ambient conditions. Building upon the previously developed SWM4-NDP model, additional sites representing oxygen lone-pairs were introduced. The geometry of the sites is assumed to be rigid. Considering the large number of adjustable parameters, simulated annealing together with polynomial fitting was used to facilitate model optimization. The new water model was shown to yield the correct self-diffusion coefficient after taking the system size effect into account, and the dimer geometry is better reproduced than in the SWM4 models. Moreover, the experimental oxygen-oxygen radial distribution is better reproduced, indicating that the new model more accurately describes the local hydrogen bonding structure of bulk phase water. This was further validated by its ability to reproduce the experimental nuclear magnetic shielding and related chemical shift of the water hydrogen in the bulk phase, a property sensitive to the local hydrogen bonding structure. In addition, comparison of the liquid properties of the SWM6 model is made with those of a number of widely used additive and polarizable models. Overall, improved balance between the description of monomer, dimer, clustered, and bulk phase water is obtained with the new model compared to its SWM4-NDP polarizable predecessor, though application of the model requires an approximately twofold increase on computational resources.