Search NASA⌕ Search

SEARCH · Search NASA

Results for “Aqua”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Americium Oxalate: An Experimental and Computational Investigation of Metal–Ligand Bonding

A novel actinide-containing coordination polymer, [Am(C 2 O 4 )(H 2 O) 3 Cl] (Am-1), has been synthesized and structurally characterized. The crystallographic analysis reveals that the structure is two-dimensional and comprised of pseudo-dimeric Am 3+ nodes that are bridged by oxalate ligands to form sheets. Each metal center is nine-coordinate, forming a distorted capped square antiprism geometry with a C 1 symmetry, and features bound oxalate, aqua, and chloro ligands. The Am 3+ –ligand bonds were probed computationally using the quantum theory of atoms in molecules nd natural localized molecular orbital approaches to investigate the underlying mechanisms and hybrid atomic orbital contributions therein. The analyses indicate that the bonds within Am-1 are predominantly ionic and the 5f shell of the Am 3+ metal centers does not add a significant covalent contribution to the bonds. Further, our bonding assessment is supported by measurements on the optical properties of Am-1 using diffuse reflectance and photoluminescence spectroscopies. The position of the principal absorption band at 507 nm ( 5 L 6' ← 7 F 0' ) is notable because it is consistent with previously reported americium oxalate complexes in solution, indicating similarities in the electronic structure and ionic bonding. Compound Am-1 is an active phosphor, featuring strong bright-blue oxalate-based luminescence with no evidence of metal-centered emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2022

Proton-Coupled Electron Transfer in a Ruthenium(II) Bipyrimidine Complex in Its Ground and Excited Electronic States

In this work, proton-coupled electron transfer (PCET) was studied for the ground and excited electronic states of a [Ru(terpy)(bpm)(OH 2 )(PF 6 ) 2 ] complex, Ru-bpm. Cyclic voltammetry measurements show that the Ru(II)-aqua moiety undergoes PCET to form a Ru(IV)-oxo moiety in the anodic region, while the bpm ligand undergoes PCET to form bpmH 2 in the cathodic region. The photophysical behavior of Ru-bpm was studied using steady-state and femtosecond transient UV–vis absorption spectroscopy, coupled with density functional theory (DFT) calculations. The lowest-lying excited state of Ru-bpm is described as a (Ru → bpm) metal-to-ligand charge-transfer (MLCT) state, while the metal-centered (MC) excited state was found computationally to be close in energy to the lowest-energy bright MLCT state (MC state was 0.16 eV above the MLCT state). The excited-state kinetics of Ru-bpm were found via transient absorption spectroscopy to be short-lived and were fit well to a biexponential function with lifetimes $τ_1$ = 4 ps and $τ_2$ = 65 ps in aqueous solution. Kinetic isotope effects of 1.75 ($τ_1$) and 1.61 ($τ_2$) were observed for both decay components, indicating that the solvent plays an important role in the excited-state dynamics of Ru-bpm. Based on the pH-dependent studies and the results from prior studies of similar Ru-complexes, we hypothesize that the 3 MLCT state forms an excited-state hydrogen-bond adduct with the solvent molecules and that this process occurs with a 4 ps lifetime. The formation of such a hydrogen-bond complex is consistent with the electronic density accumulation at the peripheral N atoms of the bpm moiety in the 3 MLCT state. The hydrogen-bonded state 3 MLCT decays to the ground state with a 65 ps lifetime. Such a short lifetime is likely associated with the efficient vibrational energy transfer from the 3 MLCT state to the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Anatase TiO 2 Surface Hydration: Infrared Vibrational Assignment and Facet-Specific Stability

Anatase TiO 2 nanoparticles support a subtle diversity of surface aqua and hydroxy coordination environments that are characteristic of distinct facets and defect sites. The disparate properties and stability of each molecularly or dissociatively absorbed water environment may provide unique reactivity for facet- and site-selective atomic layer deposition (ALD). Here we probe the temperature-dependent and facet-specific stretching frequency and stability of surface hydroxyls on faceted anatase TiO 2 nanoparticles via in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). Well-defined truncated octahedra, disks, and rods fabricated through hydrothermal synthesis present dominant (101), (001), and (010) facets, respectively. The selective dehydration of discrete local coordination environments as a function of temperature provides a potential pathway to site-selective and facet-selective ALD on metal oxide nanomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

Understanding the solvation structure of transition metal ions is important due to its broad impact on the kinetics, stability, and reactivity in a variety of applications in geochemistry, biochemistry, energy storage, and environmental chemistry. Using water as the common ligand, we study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a near-complete series of hydrated first-row transition metal ions with d-orbital occupancy ranging from d 2 to d 10 . Starting with optimized structures that were derived from an explicit solvation treatment at the density functional theory (DFT) level, we then compute the pre-edge X-ray absorption spectra at the metal ion K-edge at the time dependent density functional theory (TDDFT) and restricted active-space second-order perturbation theory (RASPT2) levels of theory. TDDFT calculations provide accurate results for spectra that are dominated by single excitations, while significant improvements were obtained with RASPT2 calculations that correctly distinguish between singly and doubly excited states, where relevant, with quantitative accuracy compared with experiment. We analyze and assign the dierent pre-edge features for each metal ion in order to reveal the impact of the variations in the d orbital occupancy on the rst-shell coordination environment. For completeness, we also report the lowest energy ligand-field d-d transitions for all the transition metal aqua ions considered in this study using complete active space second order perturbation theory (CASPT2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Evaluating the Representations of Atmospheric Rivers and Their Associated Precipitation in Reanalyses With Satellite Observations

Atmospheric rivers (ARs) are filaments of enhanced horizontal moisture transport in the atmosphere. Due to their prominent role in the meridional moisture transport and regional weather extremes, ARs have been studied extensively in recent years. Yet, the representations of ARs and their associated precipitation on a global scale remains largely unknown. In this study, we developed an AR detection algorithm specifically for satellite observations using moisture and the geostrophic winds derived from 3D geopotential height field from the combined retrievals of the Atmospheric Infrared Sounder and the Advanced Microwave Sounding Unit on NASA Aqua satellite. This algorithm enables us to develop the first global AR catalog based solely on satellite observations. The satellite–based AR catalog is then combined with the satellite–based precipitation (Integrated Muti–SatellitE Retrievals for GPM) to evaluate the representations of ARs and AR–induced precipitation in reanalysis products. Here our results show that the spreads in AR frequency and AR length distribution are generally small across data sets, while the spread in AR width is relatively larger. Reanalysis products are found to consistently underestimate both mean and extreme AR–related precipitation. However, all reanalyses tend to precipitate too often under AR conditions, especially over low latitude regions. This finding is consistent with the “drizzling” bias which has plagued generations of climate models. Overall, the findings of this study can help to improve the representations of ARs and associated precipitation in reanalyses and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Comparing Gravity Waves in a Kilometer‐Scale Run of the IFS to AIRS Satellite Observations and ERA5

Abstract Atmospheric gravity waves (GWs) impact the circulation and variability of the atmosphere. Sub‐grid scale GWs, which are too small to be resolved, are parameterized in weather and climate models. However, some models are now available at resolutions at which these waves become resolved and it is important to test whether these models do this correctly. In this study, a GW resolving run of the European Center for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS), run with a 1.4 km average grid spacing (TCo7999 resolution), is compared to observations from the Atmospheric Infrared Sounder (AIRS) instrument, on NASA's Aqua satellite, to test how well the model resolves GWs that AIRS can observe. In this analysis, nighttime data are used from the first 10 days of November 2018 over part of Asia and surrounding regions. The IFS run is resampled with AIRS's observational filter using two different methods for comparison. The ECMWF ERA5 reanalysis is also resampled as AIRS, to allow for comparison of how the high resolution IFS run resolves GWs compared to a lower resolution model that uses GW drag parametrizations. Wave properties are found in AIRS and the resampled models using a multi‐dimensional S‐Transform method. Orographic GWs can be seen in similar locations at similar times in all three data sets. However, wave amplitudes and momentum fluxes in the resampled IFS run are found to be significantly lower than in the observations. This could be a result of horizontal and vertical wavelengths in the IFS run being underestimated.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Efficient and interpretable graph network representation for angle-dependent properties applied to optical spectroscopy

Abstract Graph neural networks are attractive for learning properties of atomic structures thanks to their intuitive graph encoding of atoms and bonds. However, conventional encoding does not include angular information, which is critical for describing atomic arrangements in disordered systems. In this work, we extend the recently proposed ALIGNN (Atomistic Line Graph Neural Network) encoding, which incorporates bond angles, to also include dihedral angles (ALIGNN-d). This simple extension leads to a memory-efficient graph representation that captures the complete geometry of atomic structures. ALIGNN-d is applied to predict the infrared optical response of dynamically disordered Cu(II) aqua complexes, leveraging the intrinsic interpretability to elucidate the relative contributions of individual structural components. Bond and dihedral angles are found to be critical contributors to the fine structure of the absorption response, with distortions that represent transitions between more common geometries exhibiting the strongest absorption intensity. Future directions for further development of ALIGNN-d are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Image masks of global ship tracks for NASA MODIS data products

Ship tracks, long thin artificial cloud features formed from the pollutants in ship exhaust, are satellite-observable examples of aerosol-cloud interactions (ACI) that can lead to increased cloud albedo and thus increased solar reflectivity, phenomena of interest in solar radiation management. In addition to ship tracks being of interest to meteorologists and policy makers, their observed cloud perturbations provide benchmark evidence of ACI that remain poorly captured by climate models. To broadly analyze the effects of ship tracks, high-resolution satellite imagery data highlighting their presence are required. To support this, we provide a hand labelled dataset to serve as a benchmark for a variety of subsequent analyses. Established from a previous dataset that identified ship track presence using NASA’s MODIS Aqua satellite imager, our first-of-its-kind dataset is comprised of image masks: capturing full ship track regions, including their contours, emission points and dispersive patterns. In total, 300 images, or around 2,500 masked ship tracks, observed under varying conditions are provided, and may facilitate training of machine learning algorithms to automate extraction.

Atmospheric dynamics↗
DOE OSTI2020

Concerted proton–electron transfer reactions of manganese–hydroxo and manganese–oxo complexes

The enzymes manganese superoxide dismutase and manganese lipoxygenase use MnIII–hydroxo centres to mediate proton-coupled electron transfer (PCET) reactions with substrate. As manganese is earth-abundant and inexpensive, manganese catalysts are of interest for synthetic applications. Recent years have seen exciting reports of enantioselective C–H bond oxidation by Mn catalysts supported by aminopyridyl ligands. Such catalysts offer economic and environmentally-friendly alternatives to conventional reagents and catalysts. Mechanistic studies of synthetic catalysts highlight the role of Mn–oxo motifs in attacking substrate C–H bonds, presumably by a concerted proton–electron transfer (CPET) step. (CPET is a sub-class of PCET, where the proton and electron are transferred in the same step.) Knowledge of geometric and electronic influences for CPET reactions of Mn–hydroxo and Mn–oxo adducts enhances our understanding of biological and synthetic manganese centers and informs the design of new catalysts. In this Feature article, we describe kinetic, spectroscopic, and computational studies of MnIII–hydroxo and MnIV–oxo complexes that provide insight into the basis for the CPET reactivity of these species. Systematic perturbations of the ligand environment around MnIII–hydroxo and MnIV–oxo motifs permit elucidation of structure–activity relationships. For MnIII–hydroxo centers, electron-deficient ligands enhance oxidative reactivity. However, ligand perturbations have competing consequences, as changes in the MnIII/II potential, which represents the electron-transfer component for CPET, is offset by compensating changes in the pKa of the MnII–aqua product, which represents the proton-transfer component for CPET. For MnIV–oxo systems, a multi-state reactivity model inspired the development of significantly more reactive complexes. Weakened equatorial donation to the MnIV–oxo unit results in large rate enhancements for C–H bond oxidation and oxygen-atom transfer reactions. These results demonstrate that the local coordination environment can be rationally changed to enhance reactivity of MnIII–hydroxo and MnIV–oxo adducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The effect of sodium thiosulfate on cytotoxicity of a diimine Re(i) tricarbonyl complex

Recently, diimine Re( I ) tricarbonyl complexes have attracted great interest due to their promising cytotoxic effects. Here, we compare the cytotoxicity and cellular uptake of two Re( I ) compounds fac -[(Re(CO) 3 (bpy)(H 2 O)](CF 3 SO 3 ) ( 1 ) and Na( fac -[(Re(CO) 3 (bpy)(S 2 O 3 )])·H 2 O (bpy = 2,2'-bipyridine) ( 2 ). The Re-thiosulfate complex in 2 was characterized in two solvated crystal structures {Na( fac -[Re(CO) 3 (bpy)(S 2 O 3 )])·1.75H 2 O·C 2 H 5 OH} 4 ( 2 + 0.75H 2 O + C 2 H 5 OH) 4 and ( fac -[Re(CO) 3 (bpy)(H 2 O)]) ( fac -[Re(CO) 3 (bpy)(S 2 O 3 )])·4H 2 O ( 3 ). The cytotoxicity of 1 and 2 was tested in the MDA-MB-231 breast cancer cell line and compared with that of cisplatin. The cellular localization of the Re( I ) complexes was investigated using synchrotron-based X-ray fluorescence microscopy (XFM). Overall, the results show that replacement of the aqua ligand with thiosulfate renders the complex less toxic most likely by distrupting its cellular entry. As a result, thiosulfate could potentially have a similar chemoprotective effect against diimine fac -Re(CO) 3 complexes as it has against cisplatin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Incorporation of free halide ions stabilizes metal–organic frameworks (MOFs) against pore collapse and renders large-pore Zr-MOFs functional for water harvesting

Chemically and hydrolytically stable MOFs have shown promising water-vapor adsorption properties. However, MOFs that can simultaneously satisfy the following three requirements for effective water harvesting from low-humidity air are quite rare: (1) high water-uptake capacity; (2) hydrolytic and mechanical stability; (3) complete uptake at ~20–30% relative humidity (RH). Here we show that incorporating free halide ions is effective for enabling a representative Zr-MOF to meet these requirements for water harvesting. As-synthesized MOF-808 initially exhibits very good capacity at RH ≥ 30%, but quickly suffers large capacity losses due to water-evacuation-induced pore collapse. Via a framework-charging and free counter-ion inclusion approach, we were able to replace node-ligated formate anions with charge-neutral aqua ligands and site desired water-sorbing free-halide ions within the large pores of MOF-808 . In this study, altered samples show increased gravimetric water uptake, show beneficial shifts of water sorption isotherms toward lower water-vapor partial pressure, eliminate undesirable sorption/desorption isotherm hysteresis, and render MOF-808-Br indefinitely recyclable for ambient-temperature uptake of water vapor and lower-temperature liquid-water release.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Insights into dual-functional modification for water stability enhancement of mesoporous zirconium metal–organic frameworks

The stability of metal–organic frameworks (MOFs) in water affects their ability to function as chemical catalysts, their capacity as adsorbents for separations in water vapor presence, and their usefulness as recyclable water harvesters. Here, we have examined water stability of four node-modified variants of the mesoporous MOF, NU-1000, namely formate-, Acac-, TFacac-, and Facac-NU-1000, comparing these with node-accessible NU-1000. These NU-1000 variants present ligands grafted to NU-1000's hexa-Zr(IV)-oxy nodes by displacing terminal aqua and hydroxo ligands. Facac-NU-1000, containing the most hydrophobic ligands, showed the greatest water stability, being able to undergo at least 20 water adsorption/desorption cycles without loss of water uptake capacity. Computational studies revealed dual salutary functions of installed Facac ligands: (1) enhancement of framework mechanical stability due to electrostatic interactions; and (2) transformation and shielding of the otherwise highly hydrophilic nodes from H-bonding interactions with free water, presumably leading to weaker channel-stressing capillary forces during water evacuation – consistent with trends in free energies of dehydration across the NU-1000 variants. Water harvesting and hydrolysis of chemical warfare agent simulants were examined to gauge the functional consequences of modification and mechanical stabilization of NU-1000 by Facac ligands. The studies revealed a harvesting capacity of ~1.1 L of water vapor per gram of Facac-NU-1000 per sorption cycle. They also revealed retention of catalytic MOF activity following 20 water uptake and release cycles. This study provides insights into the basis for node-ligand-engendered stabilization of wide-channel MOFs against collapse during water removal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

In situ Raman investigation of Dy complexation in Cl-bearing aqueous solutions at 20–300 °C

Raman spectroscopy provides a versatile tool for in situ characterization of aqueous rare earth elements (REE) speciation at the molecular level. Complexation of REE with ligands such as Cl – and OH – is of particular interest for understanding the mobility of REE in NaCl-bearing hydrothermal fluids responsible for enriching REE to economic levels in nature. Raman spectroscopic studies of REE speciation in Cl-bearing aqueous fluids are primarily conducted at ambient temperature, whereas natural systems indicate temperatures of >100–600 °C. In this study, the speciation of Dy in acidic chloride-bearing hydrothermal solutions was investigated using confocal Raman spectroscopy with a new capillary Raman heating stage at 20–300 °C. Background solutions (pure water, NaCl-solutions) and solutions with 0.14–1.8 mol kg –1 dissolved DyCl 3 were sealed in quartz capillary cells. Comparison of the spectra for Dy chloride solutions with those for background solutions and the spectra for reference Dy-bearing solids was used to identify Raman bands specific to Dy–O and Dy–Cl bonds. The Raman band for the Dy–O stretching mode of hydrated Dy 3+ aqua ions was measured at 365–384 cm –1 and a Raman band for the Dy–Cl stretching modes of Dy chloride complexes was measured near 240 cm –1 . The Dy–O band decreases systematically with temperature, whereas the Dy–Cl band systematically increases, indicating a systematic increase in the stability of Dy chloride complexes with temperature. Here, this study provides the framework for expanding the use of in situ Raman spectroscopy to investigate the speciation of REE in aqueous solutions to hydrothermal conditions.

58 GEOSCIENCES↗
DOE OSTI2024

Strategically designed metal-free deep-ultraviolet birefringent crystals with superior optical properties

Finding new birefringent materials with deep-ultraviolet (DUV, λ < 200 nm) transparency is urgent, as current commercial materials cannot meet the rapidly growing demands in related application fields. Herein, three guanidinium-based compounds, C(NH 2 ) 3 CH 3 SO 3 , β-C(NH 2 ) 3 Cl, and γ-C(NH 2 ) 3 Cl, all featuring [C(NH 2 ) 3 ·X] ∞ (X = CH 3 SO 3 and Cl) pseudo layers, were designed through structural motif tailoring. Theoretical calculations indicate that these metal-free compounds all possess broad bandgaps (6.49–6.71 eV, HSE06) and remarkable birefringence (cal. 0.166–0.211 @ 1064 nm). Centimeter-sized C(NH 2 ) 3 CH 3 SO 3 crystals have been grown using a feasible aqua-solution method. Subsequently, to further optimize the properties, β/γ-C(NH 2 ) 3 Cl was remolded by further tailoring the [C(NH 2 ) 3 ] + cationic unit and the acceptor Cl - anion, and then the fourth compound NH 2 COF was theoretically constructed. Interestingly, NH 2 COF exhibits the desired coexistence of a wider bandgap (7.87 eV, HSE06) and giant birefringence (cal. 0.241 @ 1064 nm) attributed to its higher density of well-aligned birefringence-active groups (BAGs). Furthermore, among these four designed compounds, C(NH 2 ) 3 CH 3 SO 3 has been experimentally synthesized and exhibits a short UV cutoff edge. Centimeter-sized crystals have been grown using a feasible aqueous solution method. This study provides an effective strategy to optimize the density of BAGs for large birefringence and offers valuable insights into the strategic design of metal-free DUV birefringent crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗