Source author record

Philip M. Singer

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

9works
5topics
4close collaborators

Published work

9 publications

preprint2024arXiv

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

Traditional theories of the NMR autocorrelation function for intramolecular dipole pairs assume single-exponential decay, yet the calculated autocorrelation of realistic systems display 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 autocorrelation using simple, physical models within a rigorous statistical mechanical development that encompasses both rotational and translational diffusion in the same framework. 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 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.

preprint2022arXiv

Emergence of spin singlets with inhomogeneous gaps in the kagome 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. But 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 ZnCu3_{3}(OD)6_{6}FBr and herbertsmithite ZnCu3_{3}(OD)6_{6}Cl2_2 do not exhibit long range order, and they are considered the best realizations of the kagome Heisenberg antiferromagnet known to date. Here we use 63^{63}Cu nuclear quadrupole resonance combined with the inverse Laplace transform (ILT) to probe locally 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 approximately 60\% of the total spins. Theoretical models need to incorporate the role of disorder to account for the observed inhomogeneously gapped behaviour.

preprint2022arXiv

Emergence of the spin polarized domains in the kagome lattice Heisenberg antiferromagnet Zn-barlowite (Zn0.95_{0.95}Cu0.05_{0.05})Cu3_{3}(OD)6_{6}FBr

Kagome lattice Heisenberg antiferromagnets are known to be highly sensitive to perturbations caused by 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 (Zn0.95_{0.95}Cu0.05_{0.05})Cu3_{3}(OD)6_{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\% Cu2+^{2+} defect spins occupying the interlayer non-magnetic Zn2+^{2+} sites. These spin polarized domains account for ∼60\sim60\% 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.

preprint2022arXiv

Freezing of the Lattice in the Kagome Lattice Heisenberg Antiferromagnet Zn-barlowite ZnCu3_3(OD)6_6FBr

We use 79^{79}Br nuclear quadrupole resonance (NQR) to demonstrate that ultra slow lattice dynamics set in below the temperature scale set by the Cu-Cu super-exchange interaction JJ~(≃160\simeq160~K) in the kagome lattice Heisenberg antiferromagnet Zn-barlowite. The lattice completely freezes below 50~K, and 79^{79}Br NQR lineshapes 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\sim30~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.

preprint2022arXiv

Spin excitations of a proximate Kitaev quantum spin liquid realized in Cu2_2IrO3_3

Magnetic moments arranged at the corners of a honeycomb lattice are predicted to form a novel state of matter, Kitaev quantum spin liquid, under the influence of frustration effects between bond-dependent Ising interactions. Some layered honeycomb iridates and related materials, such as Na2_{2}IrO3_{3} and αα-RuCl3_{3}, are proximate to Kitaev quantum spin liquid, but bosonic spin-wave excitations associated with undesirable antiferromagnetic long-range order mask the inherent properties of Kitaev Hamiltonian. Here, we use 63^{63}Cu nuclear quadrupole resonance to uncover the low energy spin excitations in the nearly ideal honeycomb lattice of effective spin S=1/2S = 1/2 at the Ir4+^{4+} sites in Cu2_{2}IrO3_{3}. We demonstrate that, unlike Na2_{2}IrO3_{3}, Ir spin fluctuations exhibit no evidence for critical slowing down toward magnetic long range order in zero external magnetic field. Moreover, the low energy spin excitation spectrum is dominated by a mode that has a large excitation gap comparable to the Ising interactions, a signature expected for Majorana fermions of Kitaev quantum spin liquid.

preprint2020arXiv

1^1H-NMR Dipole-Dipole Relaxation in Fluids: Relaxation of Individual 1^1H-1^1H Pairs versus Relaxation of Molecular Modes

The intra-molecular 1^1H-NMR dipole-dipole relaxation of molecular fluids has traditionally been interpreted within the Bloembergen-Purcell-Pound (BPP) theory of NMR intra-molecular relaxation. The BPP theory draws upon Debye's theory for describing the rotational diffusion of the 1^1H-1^1H pair and predicts a mono-exponential decay of the 1^1H-1^1H dipole-dipole autocorrelation function between distinct spin pairs. Using molecular dynamics (MD) simulations, we show that for both nn-heptane and water this is not the case. In particular, the autocorrelation function of individual 1^1H-1^1H intra-molecular pairs itself evinces a rich stretched-exponential behavior, implying a distribution in rotational correlation times. However for the high-symmetry molecule neopentane, the individual 1^1H-1^1H intra-molecular pairs do conform to the BPP description, suggesting an important role of molecular symmetry in aiding agreement with the BPP model. The inter-molecular autocorrelation functions for nn-heptane, water, and neopentane also do not admit a mono-exponential behavior of individual 1^1H-1^1H inter-molecular pairs at distinct initial separations. We suggest expanding the auto-correlation function in terms of molecular modes, where the molecular modes do have an exponential relaxation behavior. With care, the resulting Fredholm integral equation of the first kind can be inverted to recover the probability distribution of the molecular modes. The advantages and limitations of this approach are noted.

preprint2020arXiv

Elucidating the 1^1H NMR relaxation mechanism in polydisperse polymers and bitumen using measurements, MD simulations, and models

The mechanism behind the 1^1H NMR frequency dependence of T1T_1 and the viscosity dependence of T2T_2 for polydisperse polymers and bitumen remains elusive. We elucidate the matter through NMR relaxation measurements of polydisperse polymers over an extended range of frequencies (f0=0.01↔f_0 = 0.01 \leftrightarrow 400 MHz) and viscosities (η=385↔102,000η= 385 \leftrightarrow 102,000 cP) using T1T_{1} and T2T_2 in static fields, T1T_{1} field-cycling relaxometry, and T1ρT_{1ρ} in the rotating frame. We account for the anomalous behavior of the log-mean relaxation times T1LM∝f0T_{1LM} \propto f_0 and T2LM∝(η/T)−1/2T_{2LM} \propto (η/T)^{-1/2} with a phenomenological model of 1^1H-1^1H dipole-dipole relaxation which includes a distribution in molecular correlation times and internal motions of the non-rigid polymer branches. We show that the model also accounts for the anomalous T1LMT_{1LM} and T2LMT_{2LM} in previously reported bitumen measurements. We find that molecular dynamics (MD) simulations of the T1∝f0T_{1} \propto f_0 dispersion and T2T_2 of similar polymers simulated over a range of viscosities (η=1↔1,000η= 1 \leftrightarrow 1,000 cP) are in good agreement with measurements and the model. The T1∝f0T_{1} \propto f_0 dispersion at high viscosities agrees with previously reported MD simulations of heptane confined in a polymer matrix, which suggests a common NMR relaxation mechanism between viscous polydisperse fluids and fluids under confinement, without the need to invoke paramagnetism.

preprint2020arXiv

Molecular Dynamics Simulations of NMR Relaxation and Diffusion of Heptane Confined in a Polymer Matrix

The mechanism behind the NMR surface relaxation and the large T1T_1/T2T_2 ratio of light hydrocarbons confined in the nano-pores of kerogen remains poorly understood, and consequently has engendered much debate. Towards bringing a molecular-scale resolution to this problem, we present molecular dynamics (MD) simulations of 1^1H NMR relaxation and diffusion of heptane in a polymer matrix, where the high-viscosity polymer is a model for kerogen and bitumen that provides an organic "surface" for heptane. We calculate the autocorrelation function G(t)G(t) for 1^1H-1^1H dipole-dipole interactions of heptane in the polymer matrix and use this to generate the NMR frequency (f0f_0) dependence of T1T_1 and T2T_2 relaxation times as a function of ϕC7ϕ_{C7}. We find that increasing molecular confinement increases the correlation time of the heptane molecule, which decreases the surface relaxation times for heptane in the polymer matrix. For weak confinement (ϕC7>50ϕ_{C7} > 50 vol%), we find that T1S/T2S≃1T_{1S}/T_{2S} \simeq 1. Under strong confinement (ϕC7≲ϕ_{C7} \lesssim 50 vol%), we find that the ratio T1S/T2S≳4T_{1S}/T_{2S} \gtrsim 4 increases with decreasing ϕC7ϕ_{C7}, and that the dispersion relation T1S∝f0T_{1S} \propto f_0 is consistent with previously reported measurements of polymers and bitumen. Such frequency dependence in bitumen has been previously attributed to paramagnetism, but our studies suggests that 1^1H-1^1H dipole-dipole interactions enhanced by organic nano-pore confinement dominates the NMR response in saturated organic-rich shales, without the need to invoke paramagnetism.

preprint2009arXiv

Numerical Modeling of Complex Porous Media For Borehole Applications

The diffusion/relaxation behavior of polarized spins of pore filling fluid, as often probed by NMR relaxometry, is widely used to extract information on the pore-geometry. Such information is further interpreted as an indicator of the key transport property of the formation in the oil industry. As the importance of reservoirs with complex pore geometry grows, so does the need for deeper understanding of how these properties are inter-related. Numerical modeling of relevant physical processes using a known pore geometry promises to be an effective tool in such endeavor. Using a suite of numerical techniques based on random-walk (RW) and Lattice-Boltzmann (LB) algorithms, we compare sandstone and carbonate pore geometries in their impact on NMR and flow properties. For NMR relaxometry, both laboratory measurement and simulation were done on the same source to address some of the long-standing issues in its borehole applications. Through a series of "numerical experiments" in which the interfacial relaxation properties of the pore matrix is varied systematically, we study the effect of a variable surface relaxivity while fully incorporating the complexity of the pore geometry. From combined RW and LB simulations, we also obtain diffusion-convection propagator and compare the result with experimental and network-simulation counterparts.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.