Search NASA⌕ Search

SEARCH · Search NASA

Results for “oscillator spectroscopy”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Vibronic Coupling and Exciton Chirality: Electronic and Structural Rearrangement between Helical to Zero Momentum Molecular Exciton States

Helical porphyrin aggregates are attractive macromolecular nanostructures for biomimetic light-harvesting and energy transfer in solar energy technology due to their large light absorbing cross sections and promising energy conservation benefits from chiral-induced spin selectivity. However, with these soft materials, it is imperative that we understand how the nuclear degrees of freedom impact the excitonic energy landscape and dynamics, particularly pertaining to how it influences the chiral angular momentum of the excitons. To this end, we have measured time-resolved depolarization ratios using femtosecond stimulated Raman spectroscopy to uncover the vibronic coupling guided by molecular vibrations between excitons of different angular momentum in helical tetra(sulfonatophenyl)porphyrin aggregates. We find that while transient absorption anisotropy reveals rapid (1 ps lifetime) exciton rotation from helical to achiral states in the Q-band, time-resolved vibrational depolarization ratios evolve on remarkably faster time scales (<500 fs lifetime) for the totally symmetric 1530 cm –1 vibration and the nontotally symmetric 1540 cm –1 vibration. The time-resolved depolarization ratios of the 1540 cm –1 vibration, in particular, show strong evidence of vibronic coupling and nonadiabatic exciton behavior that mediates the nonradiative transition between chiral and achiral excitons. Furthermore, our study shows that for rational design of helical molecular aggregates for exciton transport, the vibronic coupling between excitons of different angular momenta driven by molecular vibrations must be considered especially in the case where chiral-induced spin selectivity is desired.

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

Deep Learning of Dark Energy Spectroscopic Instrument Mock Spectra to Find Damped Lyα Systems

We have updated and applied a convolutional neural network (CNN) machine-learning model to discover and characterize damped Ly α systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99% for spectra that have signal-to-noise ratios (S/N) above 5 per pixel. The classification accuracy is the rate of correct classifications. This accuracy remains above 97% for lower S/N ≈1 spectra. This CNN model provides estimations for redshift and H i column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 pixel -1 . Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of baryon acoustic oscillations (BAO) is investigated. The cosmological fitting parameter result for BAO has less than 0.61% difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above 1.7%. We also compared the performances of the CNN and Gaussian Process (GP) models. Our improved CNN model has moderately 14% higher purity and 7% higher completeness than an older version of the GP code, for S/N > 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by 24% less standard deviation. A credible DLA catalog for the DESI main survey can be provided by combining these two algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Completed Sloan Digital Sky Survey IV Extended Baryon Oscillation Spectroscopic Survey: The Damped Lyα Systems Catalog

We present the characteristics of the damped Lyα (DLA) systems found in data release DR16 of the extended Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey. The DLAs were identified using the convolutional neural network of Parks et al. (2018). A total of 117,458 absorber candidates were found with 2 ≤ $z_{DLA}$ ≤ 5.5 and 19.7 ≤ log(N(H I )/cm -2 ) ≤ 22, including 57,136 DLA candidates with log(N(H I )/cm -2 ) ≥ 20.3. Mock quasar spectra were used to estimate the DLA detection efficiency and the purity of the resulting catalog. Restricting the quasar sample to bright forests, i.e., those with mean forest fluxes $\bar{f_{λ}}$ > 2 × 10 -19 W m -2 nm -1 , the efficiency and purity are greater than 90% for DLAs with column densities in the range 20.1 ≤ log(N(H I )/cm -2 ) ≤ 22.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Comprehensive anisotropic linear optical properties of the Weyl semimetals TaAs and NbAs

We report that TaAs and NbAs are two of the earliest identified Weyl semimetals that possess many intriguing optical properties, such as chirality-dependent optical excitations and giant second harmonic generation (SHG). Linear and nonlinear optics have been employed as tools to probe the Weyl physics in these crystals. Here, we extend these studies to address two important aspects: determining the complete anisotropic dielectric response and exploring if and how they can reveal essential Weyl physics. We determine the complete anisotropic dielectric functions of TaAs and NbAs by combining spectroscopic ellipsometry and density functional theory (DFT). Parameterized Lorentz oscillators are reported from 1.2 - 6 eV (experiment) and 0 - 6 eV (DFT), and good agreement is shown between them. Both linear and nonlinear optical properties have been previously reported to reveal Weyl physics. We suggest that strong optical resonances from trivial bands are the likely origin of the large optical SHG previously reported at these energies. Furthermore, by comparing the contribution of a small k -space centered around the Weyl cones with the total linear dielectric function, we find that these contributions are highly anisotropic and are < 25% of the total dielectric function below 0.5 eV; abov1 eV, these contributions are negligible. Thus, the study of Weyl physics using optical techniques requires very low energies, and even there, a careful assessment is required in distinguishing the much smaller contributions of the Weyl bands from the dominant contributions of the trivial bands and Drude response to the total dielectric function.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Effects of Shockwave Pressures on Ultrafast Vibrational Energy Transfer in BNFF, a Hydrogen-Free Energetic Material

Energy conversion in energetic materials from shock-wave-induced lattice compression to bond breaking critically depends on vibrational coupling and energy transfer between intra- and intermolecular vibrations, though the details of the mechanisms remain unknown. Herein, we indirectly tune the strength of intermolecular interactions in 3,4-bis(3-nitrofurazan-4-yl)furoxan (BNFF), a hydrogen-free energetic material characterized by van der Waals interactions, by applying high static pressure using a diamond anvil cell and monitoring vibrational energy transfer (VET) with ultrafast broadband infrared pump–probe spectroscopy. As BNFF is compressed from ambient pressure to 9 GPa, we find that VET accelerates by ∼ 0.9 ps/GPa. Density functional theory is applied in tandem with experiments to assign mode character and elucidate VET pathways. In conclusion, we find that furazan ring O–N–O vibrations, which are high-frequency detonation-relevant vibrational modes, experience increased sensitivity to lattice compression under shockwave pressures. These findings provide new mechanistic insight into how intermolecular interactions govern the rate and selectivity of VET.

Energy transfer↗
DOE OSTI2021

Multi-mode lasing in terahertz metasurface quantum-cascade VECSELs

To date, terahertz quantum-cascade vertical-external-cavity surface-emitting lasers (QC-VECSELs) have tended to oscillate in only one or two lasing modes at a time. This is due to the fact that the interaction between all of the longitudinal external cavity modes and the QC gain material is mediated through a single metasurface resonance, whose spatial overlap changes little with frequency; this suppresses spatial-hole-burning induced multi-mode operation. In this Letter, a VECSEL external cavity is demonstrated using an output coupler based upon a high-resistivity silicon etalon, which presents a periodic reflectance spectrum that is nearly matched with the external cavity mode spectrum. As the cavity length is varied, a systematic transition between a single/double-mode lasing regime and a multi-mode lasing regime is realized due to the Vernier effect. Up to nine modes lasing simultaneously with a free-spectral-range of approximately 21 GHz is demonstrated. Furthermore, this result provides a path toward the multi-mode operation necessary for eventual frequency comb operation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Surface Carboxylate Sensitivity to Electron and Hole Relaxation in Photoexcited Cadmium Sulfide Nanocrystals

Understanding how passivating surface ligands couple to excitonic states in nanocrystal photocatalysts is crucial for controlling nonradiative relaxation pathways which compete with interfacial charge transfer. Here, we report femtosecond transient infrared (IR) spectroscopy to resolve ∼100 fs ligand-exciton coupling between 1S exciton states in oleate-capped cadmium sulfide (CdS) nanocrystals and vibrational modes of surface carboxylates. Differential mid-IR spectra show distinct negative amplitude and positive photoinduced absorption signals at ∼1540 cm –1 (carboxylate asymmetric stretch) and ∼1440 cm –1 (carboxylate symmetric stretch), respectively. Fluence-dependent transient IR measurements reveal that the symmetric stretch is uniquely sensitive to picosecond Auger recombination, while the asymmetric stretch shows no analogous decay. Our results provide direct measurement of femtosecond ligand-exciton coupling in CdS nanocrystals and demonstrate how surface-bound carboxylate ligands serve as carrier-specific reporters of nanocrystal photophysics. Furthermore, these findings offer critical insights for designing and developing predictive models for ligand-mediated strategies in next-generation nanocrystal photocatalysts.

Cadmium sulfide↗
DOE OSTI2025

Cerium Dimer Anion and the Contribution of 4f Electrons to Lanthanide Metal–Metal Bonds

Direct metal−metal bonding between lanthanide atoms has been challenging to observe. We report on the first spectroscopic characterization of the cerium dimer anion (Ce 2 − ) and its neutral analog (Ce 2 ) in the gas phase, achieved using photoelectron and ultrafast spectroscopy combined with high-level quantum chemistry calculations. The electron affinity of Ce 2 is 0.24 eV, from which a dissociation energy of 2.21 eV is derived for Ce 2 − . The wave-packet dynamics upon photodetachment are studied and yield vibrational frequencies for electronically excited Ce 2 . Ce 2 − exhibits a conventional metal−metal triple bond with minimal contribution from 4f electrons. However, evidence of 4f-electron participation in bonding is identified for the low-energy excited states only 0.1 eV higher. The results challenge the assumption of inert 4f electrons in metal−metal bonding, and we propose a promising strategy for forming stable lanthanide−lanthanide bonds involving significant 4f-electron contributions.

Anions↗
DOE OSTI2021

HJ-aggregates of donor–acceptor–donor oligomers and polymers

A vibronic exciton model is developed to account for the spectral signatures of HJ-aggregates of oligomers and polymers containing donor–acceptor–donor (DAD) repeat units. In (DAD) N π-stacks, J-aggregate-promoting intrachain interactions compete with H-aggregate promoting interchain interactions. The latter includes Coulombic coupling, which arises from “side-by-side” fragment transition dipole moments as well as intermolecular charge transfer (ICT), which is enhanced in geometries with substantial overlap between donors on one chain and acceptors on a neighboring chain. J-behavior is dominant in single (DAD)N chains with enhanced intrachain order as evidenced by an increased red-shift in the low-energy absorption band along with a heightened A 1 /A 2 peak ratio, where A 1 and A 2 are the oscillator strengths of the first two vibronic peaks in the progression sourced by the symmetric quinoidal–aromatic vibration. By contrast, the positive H-promoting interchain Coulomb interactions operative in aggregates cause the vibronic ratio to attenuate, similar to what has been established in H-aggregates of homopolymers such as P3HT. An attenuated A 1 /A 2 ratio can also be caused by H-promoting ICT which occurs when the electron and hole transfer integrals are out-of-phase. In this case, the A 1 peak is red-shifted, in contrast to conventional Kasha H-aggregates. With slight modifications, the ratio formula derived previously for P3HT aggregates is shown to apply to (DAD) N aggregates as well, allowing one to determine the effective free-exciton interchain coupling from the A 1 /A 2 ratio. Applications are made to polymers based on 2T-DPP-2T and 2T-BT-2T repeat units, where the importance of the admixture of the excited acceptor state in the lowest energy band is emphasized.

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

Laser-Induced Enhancement of Acoustic Mode Quality Factor Revealed by Correlated Single-Particle Optical and Electron Microscopy

Acoustic modes in plasmonic nanostructures provide fundamental insights into their optomechanical behavior at the nanoscale, enabling emerging applications in plasmon-enhanced optomechanics, ultrasensitive sensing, and nanoscale energy transduction. Here we explore the modulation of acoustic phonon dynamics in lithographically fabricated gold nanodisks via laser-induced photothermal annealing. Using a correlated approach that utilizes both single-particle transient extinction spectroscopy and advanced electron microscopy, we directly link nanoscale structural transformations to changes in mechanical properties as probed through the coherence of the excited acoustic modes. Specifically, ultrafast pump–probe microscopy reveals an enhancement in the acoustic mode quality factor of annealed gold nanodisks, indicative of reduced damping and improved vibrational coherence. Structural characterization via scanning electron microscopy and electron backscatter diffraction confirms that photoinduced annealing results in smoother surface morphology and overall enhanced crystallinity. The improved crystalline order reduces defect and crystal boundary scattering, which we suggest as the reason underlying the lower quality factor before annealing. Furthermore, these findings demonstrate that targeted structural engineering at the nanoscale offers a powerful strategy for optimizing the optomechanical performance of plasmonic nanostructures, with broad implications for the design of next-generation nanophotonic and optomechanical systems.

Annealing (metallurgy)↗
DOE OSTI2023

The DESI Bright Galaxy Survey: Final Target Selection, Design, and Validation

Over the next 5 yr, the Dark Energy Spectroscopic Instrument (DESI) will use 10 spectrographs with 5000 fibers on the 4 m Mayall Telescope at Kitt Peak National Observatory to conduct the first Stage IV dark energy galaxy survey. At z < 0.6, the DESI Bright Galaxy Survey (BGS) will produce the most detailed map of the universe during the dark-energy-dominated epoch with redshifts of >10 million galaxies spanning 14,000 deg 2 . In this work, we present and validate the final BGS target selection and survey design. From the Legacy Surveys, BGS will target an r < 19.5 mag limited sample (BGS Bright), a fainter 19.5 < r < 20.175 color-selected sample (BGS Faint), and a smaller low-z quasar sample. BGS will observe these targets using exposure times scaled to achieve homogeneous completeness and cover the footprint three times. We use observations from the Survey Validation programs conducted prior to the main survey along with simulations to show that BGS can complete its strategy and make optimal use of "bright" time. BGS targets have stellar contamination <1%, and their densities do not depend strongly on imaging properties. BGS Bright will achieve >80% fiber assignment efficiency. Finally, BGS Bright and BGS Faint will achieve >95% redshift success over any observing condition. BGS meets the requirements for an extensive range of scientific applications. BGS will yield the most precise baryon acoustic oscillation and redshift-space distortion measurements at z < 0.4. It presents opportunities for new methods that require highly complete and dense samples (e.g., N-point statistics, multitracers). BGS further provides a powerful tool to study galaxy populations and the relations between galaxies and dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗
DOE OSTI2023

High-Pressure Structural and Thermodynamic Properties of Cerium Orthosilicates (CeSiO 4 )

Pressure-induced phase transitions from the zircon structure-type (I4 1 /amd) to the scheelite structure type (I4 1 /a) are known for many ternary oxides systems (ABO 4 ). In this work, we present the first high-pressure study on synthetic stetindite (CeSiO 4 ) by a combination of in situ high-pressure synchrotron powder X-ray diffraction up to 36 GPa, implemented with and without dual sided laser heating, and in situ high-pressure Raman spectroscopy up to 43 GPa. Two phase transitions were identified: zircon to a high-pressure low-symmetry (HPLS) phase at 15 GPa and then to a scheelite at 18 GPa. The latter from HPLS scheelite phase was found irreversible; i.e., scheelite is fully quenchable at ambient conditions, as in other zircon-type phases. The bulk moduli (K 0 ) of stetindite, HPLS, and high-pressure scheelite phases were determined, respectively, as 171(5), 105(4), and 221(40) GPa by fitting to a second-order Birch-Murnaghan equation of state. The pressure derivatives of vibrational modes and Gru''neisen parameters of the zircon-structured polymorph are similar to those of other orthosilicate minerals. In conclusion, due to the larger ionic radii of Ce 4+ , with respect to Zr 4+ , stetindite was found to possess a softer bulk modulus and undergo the phase transitions at a lower pressure than zircon (ZrSiO 4 ), such observations are consistent with what were found in coffinite (USiO 4 ).

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

Origin of mechanical and dielectric losses from two-level systems in amorphous silicon

Amorphous silicon contains tunneling two-level systems, which are the dominant energy loss mechanisms for amorphous solids at low temperatures. These two-level systems affect both mechanical and electromagnetic oscillators and are believed to produce thermal and electromagnetic noise and energy loss. However, it is unclear whether the two-level systems that dominate mechanical and dielectric losses are the same; the former relies on the coupling between phonons and two-level systems, with an elastic field coupling constant γ while the latter depends on a two-level systems dipole moment p0, which couples to the electromagnetic field. Mechanical and dielectric loss measurements as well as structural characterization were performed on amorphous silicon thin films grown by electron beam deposition with a range of growth parameters. Samples grown at 425 °C show a large reduction of mechanical loss (34 times) and a far smaller reduction of dielectric loss (2.3 times) compared to those grown at room temperature. Additionally, mechanical loss shows lower loss for thicker films, while dielectric loss shows lower loss for thinner films. Overall, analysis of these results indicate that mechanical loss correlates with atomic density, while dielectric loss correlates with dangling-bond density, suggesting a different origin for these two energy dissipation processes in amorphous silicon.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Conformational Dynamics of Bacteriochlorophyll c in Chlorosomes from the bchQ Mutant of Chlorobaculum tepidum

In contrast to the common viewpoint that bacteriochlorophyll (BChl) motion is largely absent within the chlorosome assembly, physics-based modeling points to a crucial role of the nanoscale librational motion of the macrocycle for the transfer of excitons. To elucidate this motion experimentally, compositional uniformity and high sensitivity are required. We focused on uniformly 13 C labeled chlorosome preparations from the bchQ mutant Chlorobaculum tepidum with significantly enhanced structural homogeneity. The librational motion is characterized using Rotational Echo DOuble Resonance (REDOR), and in addition, the impact of temperature on specific functionalities within BChl molecules is studied with 1-dimensional and 2-dimensional dipolar and scalar-based MAS NMR measurements. Results show the gradual freezing of the tails and side chains of the BChls with decreasing temperature. However, the librational motion analyzed by measuring the 5C–H dipolar coupling strength obtained from REDOR data sets persists at different temperatures. REDOR simulations show a close match to the experimental dephasing frequency of oscillation for a dipolar coupling strength of 17.5 ± 0.5 kHz which is considerably less than the dipolar coupling strength of 22.7 kHz in the rigid limit. Following a two-site jump model, we arrive at an estimate for BChl libration sampling at an angle of θ = 48 ± 4°, corroborating that the macrocycle indeed experiences significant librational motion on a time scale that is short compared to the NMR measurement time. This finding is in full quantitative support of the dominant rotational motion exhibited by the BChl macrocycle estimated from early MD simulations.

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

Mixing effects on spectroscopy and partonic observables of heavy mesons with logarithmic confining potential in a light-front quark model

Using the variational principle, we systematically investigate the mass spectra and wave functions of both 1⁢𝑆 and 2⁢𝑆 state heavy pseudoscalar (𝑃) and vector (𝑉) mesons within the light-front quark model. This approach incorporates a Coulomb plus logarithmic confinement potential to accurately describe the constituent quark and antiquark dynamics. Additionally, spin hyperfine interactions are introduced perturbatively to compute the masses of pseudoscalar and vector mesons. The present analyses of the 1⁢𝑆 and 2⁢𝑆 states require the consideration of mixing between them to account for empirical constraints. These constraints include the mass gap Δ⁢𝑀 𝑃 >Δ⁢𝑀 𝑉 , where Δ⁢𝑀 𝑃⁡(𝑉) =𝑀$^{2⁢𝑆}_{𝑃⁡(𝑉)}$−𝑀$^{1⁢𝑆}_{𝑃⁡(𝑉)}$ and the hierarchy of the decay constants 𝑓 1⁢𝑆 >𝑓 2⁢𝑆 . We find the optimal value of the mixing angle to be 𝜃 =1⁢8°, significantly enhancing the consistency between our spectroscopic predictions and the experimental data compiled by the Particle Data Group. Furthermore, based on the predicted mass, the newly observed resonance 𝐵 𝐽⁡ (5840) could be assigned as a 2 1⁢ 𝑆 0 state in the 𝐵 meson family. The study also reports various pertinent observables, including twist-two distribution amplitudes, electromagnetic form factors, charge radii, 𝜉 moments, and transition form factors that are found to be consistent with both available lattice simulations and experimental data. In addition, our predicted branching ratios for the channels of 𝐵 + →𝜏 + ⁢𝜈 𝜏 as well as rare decays of 𝐵 0 and 𝐵$^0_𝑠$ appear in accordance with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗