Computational simulations of nuclear magnetic resonance (NMR) experiments are essential for extracting information about molecular structure and dynamics, but are often intractable on classical computers for large molecules such as proteins and protocols such as zero-field NMR. We demonstrate the first quantum simulation of a NMR spectrum, computing the zero-field spectrum of the methyl group of acetonitrile on a trapped-ion quantum computer. We reduce the sampling cost of the quantum simulation by an order of magnitude using compressed sensing techniques. Our work opens a new practical application for quantum computation, and we show how the inherent decoherence of NMR systems may enable the simulation of classically hard molecules on near-term quantum hardware.
more »
« less
This content will become publicly available on November 15, 2024
Digital quantum simulation of NMR experiments
Simulations of nuclear magnetic resonance (NMR) experiments can be an important tool for extracting information about molecular structure and optimizing experimental protocols but are often intractable on classical computers for large molecules such as proteins and for protocols such as zero-field NMR. We demonstrate the first quantum simulation of an NMR spectrum, computing the zero-field spectrum of the methyl group of acetonitrile using four qubits of a trapped-ion quantum computer. We reduce the sampling cost of the quantum simulation by an order of magnitude using compressed sensing techniques. We show how the intrinsic decoherence of NMR systems may enable the zero-field simulation of classically hard molecules on relatively near-term quantum hardware and discuss how the experimentally demonstrated quantum algorithm can be used to efficiently simulate scientifically and technologically relevant solid-state NMR experiments on more mature devices. Our work opens a practical application for quantum computation.
more » « less- NSF-PAR ID:
- 10482946
- Publisher / Repository:
- American Association for the Advancement of Science
- Date Published:
- Journal Name:
- Science Advances
- Volume:
- 9
- Issue:
- 46
- ISSN:
- 2375-2548
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Zero- to ultralow-field (ZULF) nuclear magnetic resonance (NMR) is a version of NMR that allows studying molecules and their transformations in the regime dominated by intrinsic spin-spin interactions. While spin dynamics at zero magnetic field can be probed indirectly, J-spectra can also be measured at zero field by using non-inductive sensors, for example, optically-pumped magnetometers (OPMs). A J-spectrum can be detected when a molecule contains at least two different types of magnetic nuclei (i.e., nuclei with different gyromagnetic ratios) that are coupled via J-coupling. Up to date, no pure J-spectra of molecules featuring the coupling to quadrupolar nuclei were reported. Here we show that zero-field J-spectra can be collected from molecules containing quadrupolar nuclei with I = 1 and demonstrate this for solutions containing various isotopologues of ammonium cations. Lower ZULF NMR signals are observed for molecules containing larger numbers of deuterons compared to protons; this is attributed to less overall magnetization and not to the scalar relaxation of the second kind. We analyze the energy structure and allowed transitions for the studied molecular cations in detail using perturbation theory and demonstrate that in the studied systems, different lines in J-spectra have different dependencies on the magnetic pulse length allowing for unique on-demand zero-field spectral editing. Precise values for the 15N-1H, 14N-1H, and D-1H coupling constants are extracted from the spectra and the difference in the reduced coupling constants is explained by the secondary isotope effect. Simple symmetric cations such as ammonium do not require expensive isotopic labeling for the observation of J-spectra and, thus, may expand the applicability of ZULF NMR spectroscopy in biomedicine and energy storage.more » « less
-
more » « less
Abstract This paper reports the principal values of the13C chemical shift tensors for five nitrogen‐dense compounds (i.e., cytosine, uracil, imidazole, guanidine hydrochloride, and aminoguanidine hydrochloride). Although these are all fundamentally important compounds, the majority do not have13C chemical shift tensors reported in the literature. The chemical shift tensors are obtained from1H→13C cross‐polarization magic‐angle spinning (CP/MAS) experiments that were conducted at a high field of 18.8 T to suppress the effects of14N‐13C residual dipolar coupling. Quantum chemical calculations using density functional theory are used to obtain the13C magnetic shielding tensors for these compounds. The best agreement with experiment arises from calculations using the hybrid functional PBE0 or the double‐hybrid functional PBE0‐DH, along with the triple‐zeta basis sets TZ2P or pc‐3, respectively, and intermolecular effects modeled using large clusters of molecules with electrostatic embedding through the COSMO approach. These measurements are part of an ongoing effort to expand the catalog of accurate13C chemical shift tensor measurements, with the aim of creating a database that may be useful for benchmarking the accuracy of quantum chemical calculations, developing nuclear magnetic resonance (NMR) crystallography protocols, or aiding in applications involving machine learning or data mining. This work was conducted at the National High Magnetic Field Laboratory as part of a 2‐week school for introducing undergraduate students to practical laboratory experience that will prepare them for scientific careers or postgraduate studies.
-
Molybdenum disulfide (MoS 2 ) is a layered material with outstanding electrical and optical properties. Numerous studies evaluate the performance in sensors, catalysts, batteries, and composites that can benefit from guidance by simulations in all-atom resolution. However, molecular simulations remain difficult due to lack of reliable models. We introduce an interpretable force field for MoS 2 with record performance that reproduces structural, interfacial, and mechanical properties in 0.1% to 5% agreement with experiments. The model overcomes structural instability, deviations in interfacial and mechanical properties by several 100%, and empirical fitting protocols in earlier models. It is compatible with several force fields for molecular dynamics simulation, including the interface force field (IFF), CVFF, DREIDING, PCFF, COMPASS, CHARMM, AMBER, and OPLS-AA. The parameters capture polar covalent bonding, X-ray structure, cleavage energy, infrared spectra, bending stability, bulk modulus, Young's modulus, and contact angles with polar and nonpolar solvents. We utilized the models to uncover the binding mechanism of peptides to the MoS 2 basal plane. The binding strength of several 7mer and 8mer peptides scales linearly with surface contact and replacement of surface-bound water molecules, and is tunable in a wide range from −86 to −6 kcal mol −1 . The binding selectivity is multifactorial, including major contributions by van-der-Waals coordination and charge matching of certain side groups, orientation of hydrophilic side chains towards water, and conformation flexibility. We explain the relative attraction and role of the 20 amino acids using computational and experimental data. The force field can be used to screen and interpret the assembly of MoS 2 -based nanomaterials and electrolyte interfaces up to a billion atoms with high accuracy, including multiscale simulations from the quantum scale to the microscale.more » « less
-
We demonstrate that natural isotopic abundance 2D heteronuclear correlation (HETCOR) solid-state NMR spectra can be used to significantly reduce or eliminate the broadening of 1 H and 13 C solid-state NMR spectra of organic solids due to anisotropic bulk magnetic susceptibility (ABMS). ABMS often manifests in solids with aromatic groups, such as active pharmaceutical ingredients (APIs), and inhomogeneously broadens the NMR peaks of all nuclei in the sample. Inhomogeneous peaks with full widths at half maximum (FWHM) of ∼1 ppm typically result from ABMS broadening and the low spectral resolution impedes the analysis of solid-state NMR spectra. ABMS broadening of solid-state NMR spectra has previously been eliminated using 2D multiple-quantum correlation experiments, or by performing NMR experiments on diluted materials or single crystals. However, these experiments are often infeasible due to their poor sensitivity and/or provide limited gains in resolution. 2D 1 H– 13 C HETCOR experiments have previously been applied to reduce susceptibility broadening in paramagnetic solids and we show that this strategy can significantly reduce ABMS broadening in diamagnetic organic solids. Comparisons of 1D solid-state NMR spectra and 1 H and 13 C solid-state NMR spectra obtained from 2D 1 H– 13 C HETCOR NMR spectra show that the HETCOR spectrum directly increases resolution by a factor of 1.5 to 8. The direct gain in resolution is determined by the ratio of the inhomogeneous 13 C/ 1 H linewidth to the homogeneous 1 H linewidth, with the former depending on the magnitude of the ABMS broadening and the strength of the applied field and the latter on the efficiency of homonuclear decoupling. The direct gains in resolution obtained using the 2D HETCOR experiments are better than that obtained by dilution. For solids with long proton longitudinal relaxation times, dynamic nuclear polarization (DNP) was applied to enhance sensitivity and enable the acquisition of 2D 1 H– 13 C HETCOR NMR spectra. 2D 1 H– 13 C HETCOR experiments were applied to resolve and partially assign the NMR signals of the form I and form II polymorphs of aspirin in a sample containing both forms. These findings have important implications for ultra-high field NMR experiments, optimization of decoupling schemes and assessment of the fundamental limits on the resolution of solid-state NMR spectra.more » « less
