Engineering topics

Singer, Philip M.

Publications and source records attributed to Singer, Philip M..

8 records

DOE OSTI2024

Molecular Modes Elucidate the Nuclear Magnetic Resonance Relaxation of Viscous Fluids

The Bloembergen, Purcell, and Pound (BPP) theory of nuclear magnetic resonance (NMR) relaxation in fluids dating back to 1948 continues to be the linchpin in interpreting NMR relaxation data in applications ranging from characterizing fluids in porous media to medical imaging (MRI). The BPP theory is founded on assuming molecules are hard spheres with 1 H– 1 H dipole pairs reorienting randomly; assumptions that are severe in light of modern understanding of liquids. Nevertheless, it is intriguing to this day that the BPP theory was consistent with the original experimental data for glycerol, a hydrogen-bonding molecular fluid for which the hard-sphere-rigid-dipole assumption is inapplicable. To better understand this incongruity, atomistic molecular simulations are used to compute 1 H NMR T 1 relaxation dispersion (i.e., frequency dependence) in two contrasting cases: glycerol, and a (non hydrogen-bonding) viscosity standard. At high viscosities, simulations predict distinct functional forms of T 1 for glycerol compared to the viscosity standard, in agreement with modern measurements, yet both in contrast to BPP theory. The cause of these departures from BPP theory is elucidated, without assuming any relaxation models and without any free parameters, by decomposing the simulated T 1 response into dynamic molecular modes for both intramolecular and intermolecular interactions. The decomposition into dynamic molecular modes provides an alternative framework to understand the physics of NMR relaxation for viscous fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Theory and modeling of molecular modes in the NMR relaxation of fluids

Traditional theories of the nuclear magnetic resonance (NMR) autocorrelation function for intra-molecular dipole pairs assume a single-exponential decay, yet the calculated autocorrelation of realistic systems displays a rich, multi-exponential behavior, resulting in anomalous NMR relaxation dispersion (i.e., frequency dependence). We develop an approach to model and interpret the multi-exponential intra-molecular autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational diffusion and translational diffusion in the same framework. Here, we recast the problem of evaluating the autocorrelation in terms of averaging over a diffusion propagator whose evolution is described by a Fokker–Planck equation. The time-independent part admits an eigenfunction expansion, allowing us to write the propagator as a sum over modes. Each mode has a spatial part that depends on the specified eigenfunction and a temporal part that depends on the corresponding eigenvalue (i.e., correlation time) with a simple, exponential decay. The spatial part is a probability distribution of the dipole pair, analogous to the stationary states of a quantum harmonic oscillator. Drawing inspiration from the idea of inherent structures in liquids, we interpret each of the spatial contributions as a specific molecular mode. These modes can be used to model and predict the NMR dipole–dipole relaxation dispersion of fluids by incorporating phenomena on the molecular level. We validate our statistical mechanical description of the distribution in molecular modes with molecular dynamics simulations interpreted without any relaxation models or adjustable parameters: the most important poles in the Padé–Laplace transform of the simulated autocorrelation agree with the eigenvalues predicted by the theory

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effect of Nanoconfinement on NMR Relaxation of Heptane in Kerogen from Molecular Simulations and Measurements

Kerogen-rich shale reservoirs will play a key role during the energy transition, yet the effects of nanoconfinement on the NMR relaxation of hydrocarbons in kerogen are poorly understood. Here, we use atomistic MD simulations to investigate the effects of nanoconfinement on the 1H NMR relaxation times T 1 and T 2 of heptane in kerogen. In the case of T 1 , we discover the important role of confinement in reducing T 1 by ~3 orders of magnitude from that of bulk heptane, in agreement with measurements of heptane dissolved in kerogen from the Kimmeridge Shale, without any models or free parameters. In the case of T 2 , we discover that confinement breaks spatial isotropy and gives rise to residual dipolar coupling which reduces T 2 by ~5 orders of magnitude from the value for bulk heptane. We use the simulated T 2 to calibrate the surface relaxivity and thence predict the pore-size distribution of the organic nanopores in kerogen, without additional experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Emergence of the spin polarized domains in the kagome lattice Heisenberg antiferromagnet Zn-barlowite (Zn 0.95 Cu 0.05 )Cu 3 (OD) 6 FBr

Kagome lattice Heisenberg antiferromagnets are known to be highly sensitive to perturbations caused by the structural disorder. NMR is a local probe ideally suited for investigating such disorder-induced effects, but in practice, large distributions in the conventional one-dimensional NMR data make it difficult to distinguish the intrinsic behavior expected for pristine kagome quantum spin liquids from disorder-induced effects. Here we report the development of a two-dimensional NMR data acquisition scheme applied to Zn-barlowite (Zn 0.95 Cu 0.05 )Cu 3 (OD) 6 FBr kagome lattice, and successfully correlate the distribution of the low energy spin excitations with that of the local spin susceptibility. We present evidence for the gradual growth of domains with a local spin polarization induced by 5% Cu 2+ defect spins occupying the interlayer non-magnetic Zn 2+ sites. These spin-polarized domains account for ~60% of the sample volume at 2 K, where gapless excitations induced by interlayer defects dominate the low-energy sector of spin excitations within the kagome planes.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Thermal and concentration effects on 1 H NMR relaxation of Gd 3+ -aqua using MD simulations and measurements

We report gadolinium-based contrast agents are key in clinical MRI for enhancing the longitudinal NMR relativity (r 1 ) of hydrogen nuclei ( 1 H) in water and improving the contrast among different tissues. The importance of MRI in clinical practice cannot be gainsaid, yet the interpretation of MRI relies on models with severe assumptions, reflecting a poor understanding of the molecular-scale relaxation processes. In a step towards building a clearer understanding of the relaxation processes, here we investigate thermal and concentration effects on r 1 of the Gd 3+ -aqua complex using both semi-classical molecular dynamics (MD) simulations and measurements. We follow the MD simulation approach recently introduced by [Singer et al., Phys. Chem. Chem. Phys., 2021, 23, 20974], in which no NMR relaxation model or free-parameter is assumed to predict r 1 , thereby bringing new insights into the physics of r 1 on a molecular scale. We expand the autocorrelation function G(t) in terms of molecular modes and determine the thermal activation energies of the two largest modes, both of which are consistent with the range of literature values for rotational diffusion. We also determine the activation energies for translational diffusion and low-field electron-spin relaxation, both of which are consistent with the literature. Furthermore, we validate the MD simulations at human body temperature and concentrations of the paramagnetic ion used in clinical MRI, and we quantify the uncertainties in both simulations and measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Freezing of the Lattice in the Kagome Lattice Heisenberg Antiferromagnet Zn-Barlowite ZnCu 3 (OD) 6 FBr

Here, we use 79 Br nuclear quadrupole resonance (NQR) to demonstrate that ultraslow lattice dynamics set in below the temperature scale set by the Cu-Cu superexchange interaction J (≃160 K) in the kagome lattice Heisenberg antiferromagnet Zn-barlowite. The lattice completely freezes below 50 K, and 79 Br NQR line shapes become twice broader due to increased lattice distortions. Moreover, the frozen lattice exhibits an oscillatory component in the transverse spin echo decay, a typical signature of pairing of nuclear spins by indirect nuclear spin-spin interaction. This indicates that some Br sites form structural dimers via a pair of kagome Cu sites prior to the gradual emergence of spin singlets below ~30 K. Our findings underscore the significant roles played by subtle structural distortions in determining the nature of the disordered magnetic ground state of the kagome lattice.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Predicting 1 H NMR relaxation in Gd 3+ -aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict 1 H NMR T 1 relaxation of water from paramagnetic Gd 3+ ions in solution at 25 °C. Simulations of the T 1 relaxivity dispersion function r 1 computed from the Gd 3+ – 1 H dipole–dipole autocorrelation function agree within ≃8% of measurements in the range f 0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. Here the simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f 0 ≲ 5 MHz, the simulation overestimates r 1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r 1 at high frequencies in chelated Gd 3+ contrast-agents used for clinical MRI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Emergence of spin singlets with inhomogeneous gaps in the kagome lattice Heisenberg antiferromagnets Zn-barlowite and herbertsmithite

The kagome Heisenberg antiferromagnet formed by frustrated spins arranged in a lattice of corner-sharing triangles is a prime candidate for hosting a quantum spin liquid (QSL) ground state consisting of entangled spin singlets. However, the existence of various competing states makes a convincing theoretical prediction of the QSL ground state difficult, calling for experimental clues from model materials. The kagome lattice materials Zn-barlowite (ZnCu 3 (OD) 6 FBr) and herbertsmithite (ZnCu 3 (OD) 6 Cl 2 ) do not exhibit long-range order and are considered the best realizations of the kagome Heisenberg antiferromagnet known so far. In this study, we use 63 Cu nuclear quadrupole resonance combined with the inverse Laplace transform to locally probe the inhomogeneity of delicate quantum ground states affected by disorder. We present direct evidence for the gradual emergence of spin singlets with spatially varying excitation gaps, but even at temperatures far below the super-exchange energy scale their fraction is limited to ~60% of the total spins. Theoretical models need to incorporate the role of disorder to account for the observed inhomogeneously gapped behaviour.

36 MATERIALS SCIENCE↗