Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Probing the Atomic-Scale Structure of Amorphous Aluminum Oxide Grown by Atomic Layer Deposition

Atomic layer deposition (ALD) is a well-established technique for depositing nanoscale coatings with pristine control of film thickness and composition. The trimethylaluminum (TMA) and water (H 2 O) ALD chemistry is inarguably the most widely used and yet to date, we have little information about the atomic-scale structure of the amorphous aluminum oxide (AlO x ) formed by this chemistry. This lack of understanding hinders our ability to establish process-structure-property relationships and ultimately limits technological advancements employing AlOx made via ALD. In this work, we employ synchrotron high-energy X-ray diffraction (HE-XRD) coupled with pair distribution function (PDF) analysis to characterize the atomic structure of amorphous AlO x ALD coatings. We combine ex situ and in operando HE-XRD measurements on ALD AlO x and fit these experimental data using stochastic structural modeling to reveal variations in the Al-O bond length, Al and O coordination environment, and extent of Al vacancies as a function of growth conditions. In particular, the local atomic structure of ALD AlO x is found to change with the substrate and number of ALD cycles. The observed trends are consistent with the formation of bulk Al 2 O 3 surrounded by an O-rich surface layer. We deconvolute these data to reveal atomic-scale structural information for both the bulk and surface phases. Overall, this work demonstrates the usefulness of HE-XRD and PDF analysis in improving our understanding of the structure of amorphous ALD thin films and provides a pathway to evaluate how process changes impact the structure and properties of ALD films.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding Multi-Ion Transport Mechanisms in Bipolar Membranes

Bipolar membranes (BPMs) have the potential to become critical components in electrochemical devices for a variety of electrolysis and electrosynthesis applications. Because they can operate under large pH gradients, BPMs enable favorable environments for electrocatalysis at the individual electrodes. Critical to the implementation of BPMs in these devices is understanding the kinetics of water dissociation that occurs within the BPM as well as the co- and counter-ion crossover through the BPM, which both present significant obstacles to developing efficient and stable BPM-electrolyzers. In this study, a continuum model of multi-ion transport in a BPM is developed and fit to experimental data. Specifically, concentration profiles are determined for all ionic species, and the importance of a water-dissociation catalyst is demonstrated. The model describes internal concentration polarization and co- and counter-ion crossover in BPMs, determining the mode of transport for ions within the BPM and revealing the significance of salt-ion crossover when operated with pH gradients relevant to electrolysis and electrosynthesis. Lastly, a sensitivity analysis reveals that the performance and lifetime of BPMs can be improved substantially by using of thinner dissociation catalysts, managing water transport, modulating the thickness of the individual layers in the BPM to control salt-ion crossover, and increasing the ion-exchange capacity of the ion-exchange layers in order to amplify the water-dissociation kinetics at the interface.

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

Confined Room-Temperature Ferromagnetism in Kagome (Fe 3 Sn 2 /CoSn) Superlattices

We reveal the presence of robust ferromagnetism in topological kagome (Fe 3 Sn 2 /CoSn) ×n superlattices via polarized neutron reflectivity (PNR) and magnetization measurements. Using molecular beam epitaxy, we have synthesized superlattices of alternating ferromagnetic Fe 3 Sn 2 and paramagnetic CoSn layers, the interfacial integrity of which is confirmed using scanning transmission electron microscopy. Magnetization depth profiles, obtained by fitting the PNR data, confirm the presence of ferromagnetism in the confined Fe 3 Sn 2 layers at room temperature, with minimal thickness dependence down to 2 nm. A net magnetization is found in the near-surface region of CoSn, which we attribute to surface oxidation. Our results demonstrate that superlattice engineering is a means to confine ferromagnetism in kagome metals, enabling future studies of thickness-dependent anomalous Hall and Nernst effects as well as spin-dependent transport and tunneling in device structures.

kagome metals↗
DOE OSTI2022

Water-Accelerated Photooxidation of CH 3 NH 3 PbI 3 Perovskite

Optical absorbance is used to study the kinetics of methylammonium lead iodide (MAPbI 3 ) thin film degradation in response to combinations of moisture, oxygen, and illumination over a range of temperatures. 105 degradations were conducted over 41 unique environmental conditions. We discover that water acts synergistically with oxygen in a water-accelerated photo-oxidation (WPO) pathway. This pathway is the dominant pathway at 25 °C and is 10, 100, 1000, and >1000 times faster than dry photooxidation (DPO), degradation via hydrate formation, thermal degradation, and blue light degradation, respectively. We find that the rate determining step for DPO is proton abstraction from methylammonium while for WPO it is proton abstraction from water, which occurs at a faster rate and results in water acting as an accelerant for photooxidation of MAPbI 3 . A full kinetic rate equation is derived and fitted to the data to determine activation energies and rate constants.

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

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2022

Anthropogenic Contributions to the 2021 Pacific Northwest Heatwave

Abstract Daily maximum temperatures during the 2021 heatwave in the Pacific Northwest United States and Canada shattered century old records. Multiple causal factors, including anthropogenic climate change, contributed to these high temperatures, challenging traditional methods of attributing human influence. We demonstrate that the observed 2021 daily maximum temperatures are far above the bounds of Generalized Extreme Value distributions fitted from historical data. Hence, confidence in Granger causal inference statements about the human influence on this heatwave is low. Alternatively, we present a more conditional hindcast attribution study using two regional models. We performed ensembles of simulations of the heatwave to investigate how the event would have changed if it had occurred without anthropogenic climate change and with future warming. We found that global warming caused a ∼0.8°C–1°C increase in heatwave temperatures. Future warming would lead to a ∼5°C increase in heatwave temperature by the end of the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

In Situ Determination of Thermoelastic Properties of Magnesite at High Pressure and Temperature With Implications to Seismic Detectability of Moderately Carbonated Lithologies in the Earth's Mantle

Abstract Subduction of carbonate‐bearing oceanic plates into Earth's interior recycles carbon from the surface to the deep mantle. The subducted carbonates can significantly affect mantle properties and dynamics. Magnesite is recognized as one of the major potential carbon hosts in the deep mantle because of its stability up to deep lower mantle conditions. However, despite many previous studies on the equation of state and elastic properties of magnesite, large discrepancies still exist for its elastic moduli and their pressure and temperature derivatives. Here we report in situ density and elastic wave velocity measurements on a natural magnesite at simultaneous high pressure‐temperature conditions up to ∼8 GPa–1073 K in a multi‐anvil apparatus using ultrasonic and synchrotron X‐ray techniques. Global fitting of the data set to finite strain equations yieldsK S0 = 114.0 ± 1.2 GPa, = 4.9 ± 0.3, = (−0.019 ± 0.002) GPa/K,G 0 = 68.6 ± 0.4 GPa,G′ = 1.6 ± 0.1, and = (−0.018 ± 0.001) GPa/K. Compared to other major upper mantle and transition zone minerals, magnesite has intermediateV P , the lowestV S , and the lowest density. Thus, magnesite possesses higherV P /V S ratio than other mantle minerals in normal mantle regions, whereas this feature is less pronounced in subduction zone environments. Modeling of the velocity profiles of carbonated lithologies along different geotherms suggests that moderately‐enriched magnesite domains are unlikely to be detected seismically in the Earth's mantle.

Geochemistry & Geophysics↗
DOE OSTI2024

Permeability and Induced Polarization of Mudstones

Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.

58 GEOSCIENCES↗
DOE OSTI2025

Sound Velocity of Hematite up to 13 GPa

We performed acoustic measurements on polycrystalline hematite using ultrasonic interferometric techniques in a Kawai-type multi-anvil apparatus. The elastic moduli and their pressure derivatives were determined by fitting the experimental data to third-order finite strain equations, yielding: 𝐾 𝑆,0 = 214.8 (20) GPa, 𝐺 0 = 92.0 (2) GPa, 𝐾$^{′}_{S,0}$ = 5.70 (1), and 𝐺$^{′}_{0}$ = 0.88 (4). Modeling of the velocity profiles of hematite along hot and cold slab geotherms shows that compressional (V P ) and shear (V S ) wave velocities are 22%–24% and 53%–55% lower, respectively, compared to the Preliminary Reference Earth Model at the top of the lower mantle. Furthermore, these findings likely offer a potential explanation for the low-velocity anomalies associated with slab graveyards at the base of the lower mantle.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗
DOE OSTI2025

Néel domain walls with bistable chirality in a perpendicularly magnetized ferrimagnetic insulator

Field-free spin-orbit torque-driven domain wall motion in magnetic thin films with perpendicular magnetic anisotropy (PMA) requires the domain walls to have Néel character. Conventionally, Néel domain walls are stabilized by the Dzyaloshinskii-Moriya interaction (DMI) in ultrathin films. Here, in a europium iron garnet thin film with PMA and an additional uniaxial in-plane anisotropy, we demonstrate two bistable Néel domain wall states in the absence of DMI, and the capability to toggle the wall states with an in-plane field pulse and consequently their directions of motion under a current pulse. We present a phase diagram for the bistable Néel domain wall states as a function of in-plane field pulse width and amplitude. By fitting the experimental data to an analytical model of Néel wall reversal through the nucleation and propagation of Bloch lines, we extract the length of the initial reversed domain wall segment and Bloch line nucleation energy barrier. Current-driven motion of in-plane anisotropy stabilized Néel walls is qualitatively different from that of DMI-stabilized ones owing to the different symmetry of the effective fields that stabilize the Néel configuration. Furthermore, we present a proof of principle demonstration for 2-bit random number generation based on the stochastic reversal of domain wall chirality. These results provide critical insight into the topological energy barrier of Bloch lines and identify paths towards domain wall-based memory and computing devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Catalytic production of tetrahydropyran (THP): a biomass-derived, economically competitive solvent with demonstrated use in plastic dissolution

Tetrahydropyran (THP) is a five-carbon heterocyclic ether that is non-carcinogenic, non-peroxide forming, biodegradable, and economically competitive with tetrahydrofuran (THF) as a solvent. In this work, THP has been synthesized from renewable biomass at >99.8% selectivity and 98% yield via hydrogenation of furfural-derived 3,4-dihydropyran (DHP) over Ni/SiO 2 in a continuous flow reactor at 150–200 °C. The apparent activation energy of THP formation is 31 kJ mol -1 , and the reaction orders with respect to the partial pressures of H 2 , DHP, and THP are: 2, 1, and -0.3. The kinetic data has been fitted to a Hougen–Watson model where the rate limiting step is the hydrogenation of adsorbed DHP. Ni/SiO 2 is shown to have a low deactivation rate constant of 0.012 h -1 over 100 h time on stream and can be regenerated in situ. As a performance advantage, THP is shown to be resistant to ring opening polymerization under strongly acidic conditions that THF is not, revealing it to be a superior solvent. Further, the minimum selling price of THP is competitive with the market price of THF ($$900 – 1400 per ton ) at a DHP feedstock cost of $1000 per ton. Conductor-like Screening Model for Real Solvents (COSMO-RS) and molecular dynamics (MD) simulations with 1008 solvents and 8 common plastics have demonstrated that THP can serve as an alternative solvent to THF, 2-methyltetrahydrofuran (MeTHF), and cyclopentyl methyl ether (CPME) for plastic dissolution, especially low-density polyethylene (LDPE), polypropylene (PP), polystyrene (PS) and polyvinyl chloride (PVC). This work establishes THP as a green solvent with excellent thermal, chemical and peroxidative stability that can be used for numerous applications, including waste plastic recycling.

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

Percolative phase transition in few-layered MoSe 2 field-effect transistors using Co and Cr contacts

The metal-to-insulator phase transition (MIT) in two-dimensional (2D) materials under the influence of a gating electric field has revealed interesting electronic behavior and the need for a deeper fundamental understanding of electron transport processes, while attracting much interest in the development of next-generation electronic and optoelectronic devices. Although the mechanism of the MIT in 2D semiconductors is a topic under debate in condensed matter physics, our work demonstrates the tunable percolative phase transition in few-layered MoSe 2 field-effect transistors (FETs) using different metallic contact materials. Here, we attempted to understand the MIT through temperature-dependent electronic transport measurements by tuning the carrier density in a MoSe 2 channel under the influence of an applied gate voltage. In particular, we have examined this phenomenon using the conventional chromium (Cr) and ferromagnetic cobalt (Co) as two metal contacts. For both Cr and Co, our devices demonstrated n-type behavior with a room-temperature field-effect mobility of 16 cm 2 V −1 s −1 for the device with Cr-contacts and 92 cm 2 V −1 s −1 for the device with Co-contacts, respectively. Further, with low temperature measurements at 50 K, the mobilities increased significantly to 65 cm 2 V −1 s −1 for the device with Cr and 394 cm 2 V −1 s −1 for the device with Co-contacts. By fitting our experimental data to the percolative phase transition theory, the temperature-dependent conductivity data show a transition from an insulating-to-metallic behavior at a bias of ∼28 V for Cr-contacts and ∼20 V for Co-contacts. This cross-over of the conductivity can be attributed to an increase in carrier density as a function of the gate bias in temperature-dependent transfer characteristics. By extracting the critical exponents, we find that the transport behavior in the device with Co-contacts aligns closely with the 2D percolation theory. In contrast, the devices with Cr-contacts deviate significantly from the 2D limit at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ordered magnetic fields around the 3C 84 central black hole

Context. 3C 84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of very-long-baseline interferometry (VLBI) above the hitherto available maximum frequency of 86 GHz.Aims. Using ultrahigh resolution VLBI observations at the currently highest available frequency of 228 GHz, we aim to perform a direct detection of compact structures and understand the physical conditions in the compact region of 3C 84.Methods. We used Event Horizon Telescope (EHT) 228 GHz observations and, given the limited (u, v)-coverage, applied geometric model fitting to the data. Furthermore, we employed quasi-simultaneously observed, ancillary multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure.Results. We report the detection of a highly ordered, strong magnetic field around the central, supermassive black hole of 3C 84. The brightness temperature analysis suggests that the system is in equipartition. We also determined a turnover frequency of νm = (113 ± 4) GHz, a corresponding synchrotron self-absorbed magnetic field of BSSA = (2.9 ± 1.6) G, and an equipartition magnetic field of Beq = (5.2 ± 0.6) G. Three components are resolved with the highest fractional polarisation detected for this object (mnet = (17.0 ± 3.9)%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017–2018. We report a steeply negative slope of the spectrum at 228 GHz. We used these findings to test existing models of jet formation, propagation, and Faraday rotation in 3C 84.Conclusions. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C 84. However, systematic uncertainties due to the limited (u, v)-coverage, however, cannot be ignored. Our upcoming work using new EHT data, which offer full imaging capabilities, will shed more light on the compact region of 3C 84.Key words: techniques: high angular resolution / techniques: interferometric / galaxies: active / galaxies: individual: NGC 1275 / galaxies: jets

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The shock physics of giant impacts: Key requirements for the equations of state

We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg 2 SiO 4 ), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
DOE OSTI2022

Transducerless time domain reflectance measurement of semiconductor thermal properties

The dominant techniques for laser-based thermal measurement, time-domain thermoreflectance (TDTR) and frequency-domain thermoreflectance, generally require that samples be coated with a metal transducer layer. The properties of the transducer and its interface to the sample are then important to the interpretation of results, and heat conduction in the transducer layer can limit measurements of ultra-thin samples or in-plane conductivity. Here, we demonstrate the feasibility of measuring semiconductor samples using standard TDTR instrumentation, without the use of a metal transducer. A mathematical model accounting for non-zero laser penetration depth as well as both thermal and charge carrier effects is presented, and the ability of the model to describe and fit to experimental data for silicon and germanium samples is demonstrated. Though at this stage measurement uncertainties are larger than for traditional TDTR, we anticipate that the use of this technique will expand the range of samples that can be measured with existing TDTR setups.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical readout of singlet fission biexcitons in a heteroacene with photoluminescence detected magnetic resonance

Molecular spin systems based on photoexcited triplet pairs formed via singlet fission (SF) are attractive as carriers of quantum information because of their potentially pure and controllable spin polarization, but developing systems that offer optical routes to readout as well as initialization is challenging. Herein, we characterize the electron spin magnetic resonance change in the photoluminescence intensity for a tailored organic molecular crystal while sweeping a microwave drive up to 10 GHz in a broadband loop structure. We observe resonant transitions for both triplet and quintet spin sublevel populations showing their optical sensitivity and revealing the zero-field parameters for each. We map the evolution of these spectra in both microwave frequency and magnetic field, producing a pattern of optically detected magnetic resonance (ODMR) peaks. Fits to these data using a suitable model suggest significant spin polarization in this system with orientation selectivity. Unusual excitation intensity dependence is also observed, which inverts the sign of the ODMR signal for the triplet features, but not for the quintet. These observations demonstrate optical detection of the spin sublevel population dictated by SF and intermolecular geometry, and highlight anisotropic and multi-scale dynamics of triplet pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗