Gauge invariant calculations of nuclear magnetic shielding constants using the continuous transformation of the origin of the current density approach. II. Density functional and coupled cluster theory

The quantum mechanical current density induced in a molecule by an external magnetic field is invariant to translations of the coordinate system. This fundamental symmetry is exploited to formally annihilate the diamagnetic contribution to the current density via the approach of "continuous transformation of the origin of the current density-diamagnetic zero" (CTOCD-DZ). The relationships obtained by this method for the magnetic shielding at the nuclei are intrinsically independent of the origin of the coordinate system for any approximate computational scheme relying on the algebraic approximation. The authors report for the first time an extended series of origin-independent estimates of nuclear magnetic shielding constants using the CTOCD-DZ approach at the level of density functional theory (DFT) with four different types of functionals and unrelaxed coupled cluster singles and doubles linear response (CCSD-LR) theory. The results obtained indicate that in the case of DFT the procedure employed is competitive with currently adopted computational methods allowing for basis sets of gauge-including atomic orbitals, whereas larger differences between CTOCD-DZ and common origin CCSD-LR results are observed due to the incomplete fulfillment of hypervirial relations in standard CCSD-LR theory. It was found furthermore that the unrelaxed CCSD-LR calculations predict larger correlation corrections for the shielding constants of almost all nonhydrogen atoms in their set of molecules than the usual relaxed energy derivative CCSD calculations. Finally the results confirm the excellent performance of Keal and Tozer's third functional, in particular, for the multiply bonded systems with a lot of electron correlation, but find also that the simple local density functional gives even better results for the few singly bonded molecules in their study where correlation effects are small.     

On CHF calculations of second-order magnetic properties using the method of continuous transformation of origin of the current density

Summary A fully analytical formulation is outlined for computing molecular magnetic susceptibilities and nuclear magnetic shieldings via a continuous change of origin of the electronic current density induced by an external magnetic field. The change of origin is described in terms of a (continuous) arbitrary shift functiond(r). Coupled Hartree-Fock second-order magnetic properties of CH₄ and CO₂ molecules have been computed, using the special choiced(r)=r as generating function. A detailed analysis of results obtained with a variety of basis sets reveals that such a method is not as good as previously suggested. Large basis sets must be used to obtain accurate magnetic properties. On the other hand, all the components of theoretical nuclear magnetic shielding evaluated via this approach are independent of the origin of the vector potential. In general, theoretical magnetic susceptibilities depend linearly on the distance between different coordinate frames, but are origin independent for centre-symmetric molecules.     

Benchmarking the multipole shielding polarizability/reaction field approach to solvation against QM/MM: applications to the shielding constants of N-methylacetamide

We present a benchmark study of a combined multipole shielding polarizability/reaction field (MSP/RF) approach to the calculation of both specific and bulk solvation effects on nuclear magnetic shielding constants of solvated molecules. The MSP/RF scheme is defined by an expansion of the shielding constants of the solvated molecule in terms of electric field and field gradient property derivatives derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics simulations, thereby accounting for solute-solvent dynamical effects. The MSP/RF method is benchmarked against polarizable quantum mechanics/molecular mechanics (QM/MM) calculations. The best agreement between the MSP/RF and QM/MM approaches is found by truncating the electric field expansion in the MSP/RF approach at the linear electric field level which is due to the cancelation of errors. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the property derivatives and find that these derivatives are affected by the basis set in a way similar to the shielding constants themselves.     

The coupling constant polarizability and hyperpolarizabilty of 1J(NH) in N-methylacetamide, and its application for the multipole spin-spin coupling constant polarizability/reaction field approach to solvation

We present a benchmark study of a combined multipole spin-spin coupling constant (SSCC) polarizability/reaction field (MJP/RF) approach to the calculation of both specific and bulk solvation effects on SSCCs of solvated molecules. The MJP/RF scheme is defined by an expansion of the SSCCs of the solvated molecule in terms of coupling constant dipole and quadrupole polarizabilities and hyperpolarizabilities derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics (MD) simulations thereby accounting for solute-solvent dynamical effects. The MJP/RF method is benchmarked against polarizable QM/MM calculations for the one-bond N-H coupling constant in N-methylacetamide. The best agreement between the MJP/RF and QM/MM approaches is found by truncating the electric field expansion in the MJP/RF approach at the linear electric field level. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the coupling constant (hyper-)polarizabilities and find that they are affected by the basis set in a way similar to the coupling constants themselves.     

Density functional theory investigation of hydrogen bonding effects on the oxygen, nitrogen and hydrogen electric field gradient and chemical shielding tensors of anhydrous chitosan crystalline structure

A systematic computational investigation was carried out to characterize the 17O, 14N and 2H electric field gradient, EFG, as well as 17O, 15N, 13C and 1H chemical shielding tensors in the anhydrous chitosan crystalline structure. To include the hydrogen-bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a hexameric cluster. The computations were performed with the B3LYP method and 6-311++G(d,p) and 6-31++G(d,p) standard basis sets using the Gaussian 98 suite of programs. Calculated EFG and chemical shielding tensors were used to evaluate the 17O, 14N and 2H nuclear quadrupole resonance, NQR, and 17O, 15N, 13C and 1H nuclear magnetic resonance, NMR, parameters in the hexameric cluster, which are in good agreement with the available experimental data. The difference between the calculated NQR and NMR parameters of the monomer and hexamer cluster shows how much hydrogen bonding interactions affect the EFG and chemical shielding tensors of each nucleus. These results indicate that both O(3)-H(33)...O(5-3) and N-H(22)...O(6-4) hydrogen bonding have a major influence on NQR and NMR parameters. Also, the quantum chemical calculations indicate that the intra- and intermolecular hydrogen bonding interactions play an essential role in determining the relative orientation of EFG and chemical shielding principal components in the molecular frame axes.     

Quadrupole central transition 17O NMR spectroscopy of biological macromolecules in aqueous solution

We demonstrate a general nuclear magnetic resonance (NMR) spectroscopic approach in obtaining high-resolution (17)O (spin-5/2) NMR spectra for biological macromolecules in aqueous solution. This approach, termed quadrupole central transition (QCT) NMR, is based on the multiexponential relaxation properties of half-integer quadrupolar nuclei in molecules undergoing slow isotropic tumbling motion. Under such a circumstance, Redfield's relaxation theory predicts that the central transition, m(I) = +1/2 ? -1/2, can exhibit relatively long transverse relaxation time constants, thus giving rise to relatively narrow spectral lines. Using three robust protein-ligand complexes of size ranging from 65 to 240 kDa, we have obtained (17)O QCT NMR spectra with unprecedented resolution, allowing the chemical environment around the targeted oxygen atoms to be directly probed for the first time. The new QCT approach increases the size limit of molecular systems previously attainable by solution (17)O NMR by nearly 3 orders of magnitude (1000-fold). We have also shown that, when both quadrupole and shielding anisotropy interactions are operative, (17)O QCT NMR spectra display an analogous transverse relaxation optimized spectroscopy type behavior in that the condition for optimal resolution depends on the applied magnetic field. We conclude that, with the currently available moderate and ultrahigh magnetic fields (14 T and higher), this (17)O QCT NMR approach is applicable to a wide variety of biological macromolecules. The new (17)O NMR parameters so obtained for biological molecules are complementary to those obtained from (1)H, (13)C, and (15)N NMR studies.     

Molecular magnetic properties via formal annihilation of paramagnetic contribution to electronic current density

A novel procedure for calculating magnetic susceptibilities and nuclear magnetic shieldings in molecules is outlined, based on formal annihilation of transverse paramagnetic contribution to quantum mechanical current density induced within the electron cloud by an external homogeneous, static magnetic field. Within this method all the components of nuclear magnetic shielding are independent of a gauge translation, in any calculation relying on the algebraic approximation, irrespective of size and quality of the gaugeless basis set adopted; magnetic susceptibilities are invariant for center-symmetric molecules (virtual invariance is actually observed for molecules of arbitrary symmetry). Large basis set calculations of near-Hartree-Fock magnetic properties and maps describing a current density field of acetylene molecule carried out via the new procedure are compared with corresponding ones adopting the common origin-coupled Hartree-Fock approach and methods formally annihilating the diamagnetic contribution to the current density. © 1996 John Wiley & Sons, Inc.     

Electric field effects on the shielding constants of noble gases: a four-component relativistic Hartree-Fock study

Second derivatives of nuclear shielding constants with respect to an electric field, i.e., shielding polarizabilities, have been calculated for the noble gas atoms from helium to xenon. The calculations have been carried out using the four-component relativistic Hartree-Fock method. In order to assess the importance of the individual relativistic corrections, the shielding polarizabilities have also been calculated at the nonrelativistic Hartree-Fock level, with spin-orbit and scalar (Darwin and mass-velocity) effects having been established by perturbative methods. Electron correlation effects have been estimated using the second-order polarization propagator approach. The relativistic effects on the tensor components of the shielding polarizabilities are found to be larger and changing less regularly with the atomic number than for the shielding constant itself. However, there is a partial cancellation of the contributions to the parallel and perpendicular components of the shielding polarizability and as a consequence the mean shielding polarizability is far less affected than the individual components.     

On the existence of a natural common gauge-origin for the calculation of magnetic properties of atoms and molecules via gaugeless basis sets

It is proven that, within the conventional approach using a common origin and gaugeless basis sets for the calculation of atomic magnetizability and Larmor current density induced by an external magnetic field, the natural gauge origin coincides with the nucleus. Recipes for defining an optimal gauge origin for the calculation of magnetizability and magnetic shielding at the nuclei of a molecule are given. Within the common origin approach, the paramagnetic contributions to the components of magnetic tensors of a molecule are represented by a minimum number of non-vanishing parameters if the gauge origin is chosen at a point characterized by the total molecular symmetry, e.g., the center of electronic charge for magnetizabilities. It is shown that total values of diagonal components of the magnetic shielding tensor σ(I) at a nucleus I in a molecule, as well as separate diamagnetic σ(dI) and paramagnetic σ(pI) contributions, calculated via the common origin method, are origin independent for a number of local point group symmetries. The diagonal components (and the average value) of σ(I) depend on the gauge origin only for nuclear site symmetries C(1), C(s), C(n), C(nv), n = 2, 3.... Group-theoretical methods show interesting features, e.g., for S(4) local symmetry, in a coordinate transformation, the paramagnetic contribution to the zz component and to the trace of the shielding tensor is origin independent, whereas the xx and yy components mix into one another, in such a way that their sum remains constant.     

Chiral discrimination via nuclear magnetic shielding polarisabilities from NMR spectroscopy: theoretical study of (R(a))-1,3-dimethylallene, (2R)-2-methyloxirane, and (2R)-N-methyloxaziridine

Three medium-size optically active molecules have been studied to make a guess at candidates suitable for chiral discrimination in an isotropic medium via nuclear magnetic resonance spectroscopy. The criterion for experimental detection is given by the magnitude of the isotropic part of nuclear magnetic shielding polarisability tensors, related to a pseudoscalar of opposite sign for the two enantiomers. The pseudoscalar shielding polarisability at the (17)O nucleus in N-methyloxaziridine, calculated at the Hartree-Fock level, is approximately 7.8 x10(-)(17) mV(-)(1). To obtain an experimentally observable magnetic field induced at the (17)O nucleus in N-methyloxaziridine, electric fields as large as approximately 10(7) - 10(8) Vm(-)(1) should be applied to the probe. The molecular electric dipole moment induced by precession of the magnetic dipole of the (17)O nucleus in a magnetic field of 10 T is, in absolute value, approximately 8.8 x 10(-)(42) Cm. The estimated rf-voltage at a resonance circuit is approximately 10 nV. Smaller values have been estimated for N, C, and H nuclei in 1,3-dimethylallene and 2-methyloxirane.     

Calculation of dipole-shielding polarizabilities (sigma(alphabetagamma)I): the influence of uniform electric field effects on the shielding of backbone nuclei in proteins

A significant contribution to the chemical shielding of a nucleus can arise from uniform electric fields that act to distort the electronic charge distribution surrounding a nucleus and, hence, affect the nuclear shielding. It has been shown by Buckingham (Buckingham, A. D. Can. J. Chem. 1960, 38, 300) that the nuclear magnetic shielding tensor sigmaalphabetaI of a nucleus in the presence of an external weak static uniform electric field E may be expanded using sigmaalphabetaI(E) = sigmaalphabetaI + sigmaalphabetagammaIEgamma + 1/2sigmaalphabetagammadeltaIEgammaEdelta + sigmaalphabeta,gammadeltaIEgammadelta ... The third rank tensor sigmaalphabetagammaI is referred to as the dipole-shielding polarizability and describes the nonlinear response of the electron cloud to first order in E, muI, and B0. We report calculations of sigmaalphabetagammaI for the N, HN, and C' nuclei in N-methyl acetamide (NMA) and show that these tensors can be used to provide considerable insight into the behavior of uniform electric fields upon the shielding of backbone nuclei in proteins. The sigmaalphabetagammaI values for the N, HN, and C' of NMA were calculated using the continuous transformation of the origin of the current density (CTOCD) scheme with the diamagnetic contribution set to zero (CTOCD-DZ). Values are given for the individual tensor components of sigmaalphabetagammaI for each nucleus. To test that the calculations have provided a reasonable estimate for the sigmaalphabetagammaI of N, HN, and C' nuclei in proteins, a pH titration was performed using Hen Lysozyme (HEWL). The pH-induced isotropic shielding changes for the C', N, and HN nuclei in some peptide bonds close to E35 ( approximately <8 A) were extracted from sets of fitted titration curves. Assuming the experimental shielding changes arise solely from uniform electric field effects caused by the deprotonation of E35, without any other pH-induced structural alterations which might lead to a shielding change, the experimental shielding differences were compared to those calculated via the product Agamma(I).Egamma where Agamma(I) = (1/3)sigmaalphaalphagammaIota. The agreement with the experimental data is in many cases reasonable and suggests that, within the Buckingham formalism, the complete sigmaalphabetagammaI tensors reported here will be helpful to resolve the importance of uniform electric fields upon isotropic and anisotropic shielding in proteins and their complexes.     

Rearrangement correction to higher polarizabilities of two-electron systems

In this paper we include the rearrangement correction (discussed in the preceding paper) in a coupled Hartree–Fock (CHF) calculation of atomic hyperpolarizabilities and other related properties. We have studied the effect of these corrections on properties like electric dipole hyperpolarizabilities, uniform electric field quadrupole polarizabilities and shielding factors in two-electron ions and have noticed significant changes in the computed values over the CHF results.     

Beyond NICS: estimation of the magnetotropicity of inorganic unsaturated planar rings

A simple classical model of magnetic-field induced electron flow is used to evaluate the ring current strength for a few inorganic monocyclic compounds: B(3)H(3)N(3), B(3)H(3)O(3), P(6), N(6), Si(6)H(6), N, Al and H(6). It is shown that, for these neutral and charged systems, sustaining delocalized electron currents in the presence of a magnetic field B(ext) orthogonal to the σ(h) plane, the out-of-plane component of the nuclear magnetic shielding along the central axis is connected to the out-of-plane magnetizability by a simple equation, involving the radius of an average loop of current. A novel estimate of this effective radius is provided. Reliable ring current susceptibilities (that is, current strengths) can be evaluated by a simple relationship, using the out-of-plane components of nuclear shielding and magnetizability tensors. The accuracy of the current susceptibilities calculated by the classical model is established by comparison with corresponding ab initio estimates obtained by integrating the quantum mechanical current-density vector field. The out-of-plane components of nuclear shielding and magnetizability are both strongly biased by the molecular geometry. Their combined use to estimate the ring current susceptibility offers a quantifier of magnetotropicity more reliable than (i) the ξ(∥) out-of-plane component of magnetizability, (ii) the σ(∥)(CM) out-of-plane component of the magnetic shielding at the center of mass, widely reported as NICS(∥)(0) = -σ(∥)(CM). The inadequacy of these commonly adopted magnetotropicity measures is demonstrated by comparing a set of related molecules, C(6)H(6) and Si(6)H(6), N(6) and P(6).     

Calculation of static and dynamic linear magnetic response in approximate time-dependent density functional theory

We report implementations and results of time-dependent density functional calculations (i) of the frequency-dependent magnetic dipole-magnetic dipole polarizability, (ii) of the (observable) translationally invariant linear magnetic response, and (iii) of a linear intensity differential (LID) which includes the dynamic dipole magnetizability. The density functional calculations utilized density fitting. For achieving gauge-origin independence we have employed time-periodic magnetic-field-dependent basis functions as well as the dipole velocity gauge, and have included explicit density-fit related derivatives of the Coulomb potential. We present the results of calculations of static and dynamic magnetic dipole-magnetic dipole polarizabilities for a set of small molecules, the LID for the SF6 molecule, and dispersion curves for M-hexahelicene of the origin invariant linear magnetic response as well as of three dynamic polarizabilities: magnetic dipole-magnetic dipole, electric dipole-electric dipole, and electric dipole-magnetic dipole. We have also performed comparison of the linear magnetic response and magnetic dipole-magnetic dipole polarizability over a wide range of frequencies for H2O and SF6.     

Theoretical predictions of nuclear magnetic resonance parameters in a novel organo-xenon species: chemical shifts and nuclear quadrupole couplings in HXeCCH

We calibrate the methodology for the calculation of nuclear magnetic resonance (NMR) properties in novel organo-xenon compounds. The available state-of-the-art quantum-chemical approaches are combined and applied to the HXeCCH molecule as the model system. The studied properties are (129)Xe, (1)H, and (13)C chemical shifts and shielding anisotropies, as well as (131)Xe and (2)H nuclear quadrupole coupling constants. The aim is to obtain, as accurately as currently possible, converged results with respect to the basis set, electron correlation, and relativistic effects, including the coupling of relativity and correlation. This is done, on one hand, by nonrelativistic correlated ab initio calculations up to the CCSD(T) level and, on the other hand, for chemical shifts and shielding anisotropies by the leading-order relativistic Breit-Pauli perturbation theory (BPPT) with correlated ab initio and density-functional theory (DFT) reference states. BPPT at the uncorrelated Hartree-Fock level as well as the corresponding fully relativistic Dirac-Hartree-Fock method are found to be inapplicable due to a dramatic overestimation of relativistic effects, implying the influence of triplet instability in this multiply bonded system. In contrast, the fully relativistic second-order Moller-Plesset perturbation theory method can be applied for the quadrupole coupling, which is a ground-state electric property. The performance of DFT with various exchange-correlation functionals is found to be inadequate for the nonrelativistic shifts and shielding anisotropies as compared to the CCSD(T) results. The relativistic BPPT corrections to these quantities can, however, be reasonably predicted by DFT, due to the improved triplet excitation spectrum as compared to the Hartree-Fock method, as well as error cancellation within the five main BPPT contributions. We establish three computationally feasible models with characteristic error margins for future calculations of larger organo-xenon compounds to guide forthcoming experimental NMR efforts. The predicted (129)Xe chemical shift in HXeCCH is in a novel range for this nucleus, between weakly bonded or solvated atomic xenon and xenon in the hitherto characterized molecules.     

Ab initio study of the ground state properties of molecular oxygen

The electric and magnetic properties of the ground state of oxygen molecule are calculated by multiconfiguration self-consisted field (MCSCF) method and compared with experimental data: the quadrupole moment, polarizability, the 17O nuclear quadrupole coupling constant, magnetizability tensor, nuclear spin-rotation coupling constant and rotational g factor. The last two constants are calculated for all possible isotope modifications. The rotational, ESR and NMR spectra are discussed. Fermi-contact hyperfine coupling parameter is additionally estimated by different methods. The NMR chemical shielding tensor for 17O16O species at high temperature limit (without electron spin contribution) is predicted. Potential energy curves for 10 excited bound states and the internuclear distance dependence of the studied properties are also presented.     

Effects of adsorbate molecules on the quadrupolar interaction of framework aluminum atoms in dehydrated zeolite H,Na-Y

The effect of adsorbate molecules on the quadrupolar interaction of framework aluminum atoms with the electric field gradient in dehydrated zeolite H,Na-Y has been studied by (27)Al MAS NMR and (27)Al MQMAS NMR spectroscopy at magnetic fields of 9.4 and 17.6 T. Upon adsorption of molecules interacting with bridging OH groups by hydrogen bonds (acetonitrile and acetone), the quadrupole coupling constant of framework aluminum atoms was found to decrease from 16.0 MHz (unloaded zeolite) to 9.4 MHz. Adsorption of molecules, which cause a proton transfer from the zeolite framework to the adsorbates (ammonia and pyridine), reduces the quadrupole coupling constant to 3.8 MHz for coverages of 0.5-2 molecules per bridging OH group. The experiments indicate that the quadrupole coupling constant of framework aluminum atoms in dehydrated zeolite H,Na-Y reflects the chemical state of adsorbate complexes formed at bridging OH groups. In agreement with earlier investigations it was found that a proton affinity of the adsorbate molecules of PA = 812-854 kJ/mol is necessary to induce a proton transfer from the zeolite framework to the adsorbed compounds. This proton transfer is accompanied by a strong improvement of the tetrahedral symmetry of zeolitic framework AlO(4) tetrahedra and a decrease of the electric field gradient.     

Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.     

Characterization of divalent metal metavanadates by 51V magic-angle spinning NMR spectroscopy of the central and satellite transitions

51V quadrupole coupling and chemical shielding tensors have been determined from 51V magic-angle spinning (MAS) NMR spectra at a magnetic field of 14.1 T for nine divalent metal metavanadates: Mg(VO3)2, Ca(VO3)2, Ca(VO3)(2).4H2O, alpha-Sr(VO3)2, Zn(VO3)2, alpha- and beta-Cd(VO3)2. The manifold of spinning sidebands (ssbs) from the central and satellite transitions, observed in the 15V MAS NMR spectra, have been analyzed using least-squares fitting and numerical error analysis. This has led to a precise determination of the eight NMR parameters characterizing the magnitudes and relative orientations of the quadrupole coupling and chemical shielding tensors. The optimized data show strong similarities between the NMR parameters for the isostructural groups of divalent metal metavanadates. This demonstrates that different types of metavanadates can easily be distinguished by their anisotropic NMR parameters. The brannerite type of divalent metal metavanadates exhibits very strong 51V quadrupole couplings (i.e., CQ = 6.46-7.50 MHz), which reflect the highly distorted octahedral environments for the V5+ ion in these phases. Linear correlations between the principal tensor elements for the 51V quadrupole coupling tensors and electric field gradient tensor elements, estimated from point-monopole calculations, are reported for the divalent metal metavanadates. These correlations are used in the assignment of the NMR parameters for the different crystallographic 51V sites of Ca(VO3)(2).4H2O, Pb(VO3)2, and Ba(VO3)2. For alpha-Sr(VO3)2, with an unknown crystal structure, the 51V NMR data strongly suggest that this metavanadate is isostructural with Ba(VO3)2, for which the crystal structure has been reported. Finally, the chemical shielding parameters for orthovanadates and mono- and divalent metal metavanadates are compared.