Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spin relaxation”

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 577 records · Page 32

Non-Synchronous Rotation on Europa Driven by Ocean Currents

It has been suggested that the ice shell of Jupiter's moon Europa may drift non-synchronously due to tidal torques. Here we argue that torques applied by the underlying ocean are also important and can result in non-synchronous rotation. The resulting spin rate can be slightly slower than the synchronous angular rate that would have kept the same point of the ice shell facing Jupiter. We develop an ice shell rotation model, driven by ocean stress calculated using a high-resolution state-of-the-art ocean general circulation model, and take into account the viscoelastic deformation of the ice shell. We use the ice shell model results together with observed limits on the ice shell drift speed to constrain ice shell parameters such as effective viscosity, which is currently uncertain by at least four orders of magnitude. Our results suggest, at best, sluggish ice shell convection. Depending on the relaxation time scale of the ice shell and on the ocean currents, the ice shell may exhibit negligible drift, constant drift, or oscillatory drift superimposed on random fluctuations. The expected rotation rate exceeds ~30 m/yr; future spacecraft observations can be used to test these predictions and yield insight into the properties of the ice shell and underlying ocean.

15 GEOTHERMAL ENERGY↗

Understanding Interfacial Block Copolymer Structure and Dynamics

Here, block copolymer (BCP) structure and dynamics were studied using small-angle neutron scattering (SANS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations to obtain a fundamental understanding of the impact of an interfacial block on chain dynamics. A glassy block acted as the interface, and the dynamics of a rubbery block was studied. The rubbery block was protonated near the interface in one sample and near the chain end in another sample to observe the interfacial effect on the rubbery polymer. Analysis of the structure and dynamics revealed that the interfacial rubbery block was confined in layered morphologies and exhibited much slower dynamics than the chain-end rubbery block that was dispersed in the rubbery matrix. The interfacial rubbery block showed weaker dynamical relaxation than that at the chain end, and it also had critically important length scale dependence. Dynamical slowing was only observed at length scales significantly larger than the characteristic segmental length, and the disparity between interfacial and chain-end dynamics increased with increasing length.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Highly Ordered Nitroxide Side Chain for Distance Mapping and Monitoring Slow Structural Fluctuations in Proteins

Abstract Site-directed spin labeling electron paramagnetic resonance (SDSL-EPR) is an established tool for exploring protein structure and dynamics. Although nitroxide side chains attached to a single cysteine via a disulfide linkage are commonly employed in SDSL-EPR, their internal flexibility complicates applications to monitor slow internal motions in proteins and to structure determination by distance mapping. Moreover, the labile disulfide linkage prohibits the use of reducing agents often needed for protein stability. To enable the application of SDSL-EPR to the measurement of slow internal dynamics, new spin labels with hindered internal motion are desired. Here, we introduce a highly ordered nitroxide side chain, designated R9, attached at a single cysteine residue via a non-reducible thioether linkage. The reaction to introduce R9 is highly selective for solvent-exposed cysteine residues. Structures of R9 at two helical sites in T4 Lysozyme were determined by X-ray crystallography and the mobility in helical sequences was characterized by EPR spectral lineshape analysis, Saturation Transfer EPR, and Saturation Recovery EPR. In addition, interspin distance measurements between pairs of R9 residues are reported. Collectively, all data indicate that R9 will be useful for monitoring slow internal structural fluctuations, and applications to distance mapping via dipolar spectroscopy and relaxation enhancement methods are anticipated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decorrelated singlet and triplet exciton delocalization in acetylene-bridged Zn-porphyrin dimers

The controlled delocalization of molecular excitons remains an important goal towards the application of organic chromophores in processes ranging from light-initiated chemical transformations to classical and quantum information processing. In this study, we present a methodology to couple optical and magnetic spectroscopic techniques and assess the delocalization of singlet and triplet excitons in model molecular chromophores. By comparing the steady-state and time-resolved optical spectra of Zn-porphyrin monomers and weakly coupled dimers, we show that we can use the identity of substituents bound at specific positions of the macromolecules' rings to control the inter-ring delocalization of singlet excitons stemming from their B states through acetylene bridges. While broadened steady-state absorption spectra suggest the presence of delocalized B state excitons in mesityl-substituted Zn-tetraphenyl porphyrin dimers (Zn 2 U-D), we confirm this conclusion by measuring an enhanced ultrafast non-radiative relaxation from these inter-ring excitonic states to lower lying electronic states relative to their monomer. In contrast to the delocalized nature of singlet excitons, we use time-resolved EPR and ENDOR spectroscopies to show that the triplet states of the Zn-porphyrin dimers remain localized on one of the two macrocyclic sub-units. We use the analysis of EPR and ENDOR measurements on unmetallated model porphyrin monomers and dimers to support this conclusion. The results of DFT calculations also support the interpretation of localized triplet states. These results demonstrate researchers cannot conclude triplet excitons delocalize in macromolecular based on the presence of spatially extended singlet excitons, which can help in the design of chromophores for application in spin conversion and information processing technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unifying femtosecond and picosecond single-pulse magnetic switching in Gd-Fe-Co

Many questions are still open regarding the physical mechanisms behind the magnetic switching in Gd-Fe-Co alloys by single optical pulses. Phenomenological models suggest a femtosecond scale exchange relaxation between sublattice magnetization as the driving mechanism for switching. The recent observation of thermally induced switching in Gd-Fe-Co by using both several picosecond optical laser pulse as well as electric current pulses has questioned this previous understanding. This has raised the question of whether or not the same switching mechanics are acting at the femtosecond and picosecond scales. In this work, we aim at filling this gap in the understanding of the switching mechanisms behind thermal single-pulse switching. To that end, we have studied experimentally thermal single-pulse switching in Gd-Fe-Co alloys, for a wide range of system parameters, such as composition, laser power, and pulse duration. We provide a quantitative description of the switching dynamics using atomistic spin dynamics methods with excellent agreement between the model and our experiments across a wide range of parameters and timescales, ranging from femtoseconds to picoseconds. Furthermore, we find distinct element-specific damping parameters as a key ingredient for switching with long picosecond pulses and argue that switching with pulse durations as long as 15 ps is possible due to a low damping constant of Gd. Our findings can be easily extended to speed up dynamics in other contexts where ferrimagnetic Gd-Fe-Co alloys have been already demonstrated to show fast and energy-efficient processes, e.g., domain-wall motion in a track and spin-orbit torque switching in spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evolution of the Cellulose Microfibril through Gamma-Valerolactone-Assisted Co-Solvent and Enzymatic Hydrolysis

Biomass recalcitrance during deconstruction remains a key bottleneck to affordable biomass processing technologies. A clear connection between the cell wall structure and biomass deconstruction is necessary to understand how lignocellulosic material is broken down to valuable monomeric components. Here, we monitor changes in the cellulose microfibril domains of poplar, sorghum, and switchgrass throughout gamma-valerolactone (GVL)–water co-solvent pretreatment and enzymatic hydrolysis using solid-state 13 C cross-polarization magic angle spinning nuclear magnetic resonance spectroscopy (CP/MAS 13 C-NMR) and wide-angle X-ray scattering (WAXS). Spectral fitting of NMR peaks corresponding to different cellulose microenvironments at the C 4 carbon center suggests that a mildly acidic GVL–water co-solvent pretreatment of poplar leads to nearly full removal of xylan–cellulose linkages, which primes the cellulose for enzymatic attack. The spectral fitting also suggests that the pretreatment causes significant depletion of the inaccessible fibril surface domains with an increase in more thermally stable crystalline resonances (I β ). WAXS confirmed a decrease in the lattice spacing between (200) crystalline planes with increasing co-solvent pretreatment severity. These results are interpreted as an opening of bound microfibril surfaces previously inaccessible to the co-solvent system, which leaves behind a more thermally stable, crystalline domain that is potentially prone to relaxation and recrystallization. Full conversion of residual GVL-pretreated biomass was achieved after the GVL co-solvent pretreatment at 140 °C using a commercial enzyme cocktail, CTec2, which contains different cellulases and other enzymes. Spectral fitting of enzymatically hydrolyzed samples by a single engineered cellulase, CelR, suggests that the residual cellulose recalcitrance is mainly due to the inability of CelR to digest the I β crystalline domain present in pretreated samples. This work helps to provide new information regarding the structure of the cell wall and recalcitrance throughout GVL–water mild acidolysis and CelR enzymatic biomass deconstruction by tracking the evolution of structural domains within the cellulose microfibril. This work further directs recommendations for improving the conversion and sugar yields in future studies. Finally, our findings inform inquiry into larger questions of cellulose recalcitrance through GVL pretreatment and CelR enzymatic hydrolysis and give insight into subsequent required steps for full cellulose conversion with attention to the most recalcitrant cellulose structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear Nano‐Imaging of Interlayer Coupling in 2D Graphene‐Semiconductor Heterostructures

The emergent electronic, spin, and other quantum properties of 2D heterostructures of graphene and transition metal dichalcogenides are controlled by the underlying interlayer coupling and associated charge and energy transfer dynamics. However, these processes are sensitive to interlayer distance and crystallographic orientation, which are in turn affected by defects, grain boundaries, or other nanoscale heterogeneities. This obfuscates the distinction between interlayer charge and energy transfer. Here, nanoscale imaging in coherent four‐wave mixing (FWM) and incoherent two‐photon photoluminescence (2PPL) is combined with a tip distance‐dependent coupled rate equation model to resolve the underlying intra‐ and inter‐layer dynamics while avoiding the influence of structural heterogeneities in mono‐ to multi‐layer graphene/WSe 2 heterostructures. With selective insertion of hBN spacer layers, it is shown that energy, as opposed to charge transfer, dominates the interlayer‐coupled optical response. From the distinct nano‐FWM and ‐2PPL tip‐sample distance‐dependent modification of interlayer and intralayer relaxation by tip‐induced enhancement and quenching, an interlayer energy transfer time of τ ET ≈ ( 0 . 35 − 0.15 + 0.65 ) $\tau _{\rm ET} \approx (0.35^{+0.65}_{-0.15})$ ps consistent with recent reports is derived. As a local probe technique, this approach highlights the ability to determine intrinsic sample properties even in the presence of large sample heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗

Ferromagnetic Exchange and Slow Magnetic Relaxation in Cobalt Bis(1,2-dithiolene)-Bridged Dilanthanide Complexes

The construction of multinuclear lanthanide-based molecules with significant magnetic exchange interactions represents a key challenge in the realization of single-molecule magnets with high operating temperatures. Here, we report the synthesis and magnetic characterization of two series of heterobimetallic compounds, (Cp* 2 Ln) 2 (μ-Co(pdt) 2 ) (Ln = Y 3+ , Gd 3+ , Dy 3+ ; pdt 2– = 1,2-diphenylethylenedithiolate) and [K(18-crown-6)][(Cp* 2 Ln) 2 (μ-Co(pdt) 2 )] (Ln = Y 3+ , Gd 3+ ), featuring two lanthanide centers bridged by a cobalt bis(1,2-dithiolene) complex. Dc magnetic susceptibility data collected for the Gd congeners indicate significant Gd–Co ferromagnetic exchange interactions with fits affording J = +11.5 and +7.33 cm –1 , respectively. Magnetization decay and ac magnetic susceptibility measurements carried out on the single-molecule magnet (Cp* 2 Dy) 2 (μ-Co(pdt) 2 ) reveal full suppression of quantum tunneling and open-loop hysteresis persisting up to 3.5 K. These results, along with those of high-field EPR spectroscopy, suggest that transition metalloligands can enforce strong exchange interactions with adjacent lanthanide centers while maintaining a geometry that preserves molecular anisotropy. Furthermore, the magnetic properties of [K(18-crown-6)][(Cp* 2 Gd) 2 (μ-Co(pdt) 2 )] show that increasing the spin of the ground state of the bridging complex may be a viable alternative to increasing J in obtaining well-isolated, strongly coupled magnetic ground states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗

The Advanced X-ray Spectroscopy and Imaging Observatory (AXSIO)

Following recommendations from the 2010 "New Worlds, New Horizons" (NWNH) report, the Advanced X-ray Spectroscopy and Imaging Observatory (AXSIO) concept streamlines the International X-ray Observatory (IXO) mission to concentrate on the science objectives that are enabled by high-resolution spectroscopic capabilities. AXSIO will trace orbits close to the event horizon of black holes, measure black hole spin for tens of supermassive black holes (SMBH), use spectroscopy to characterize outflows and the environment of AGN during their peak activity, observe 5MBH out to redshift z=6, map bulk motions and turbulence in galaxy clusters, find the missing baryons in the cosmic web using background quasars, and observe the process of cosmic feedback where black holes and supernovae inject energy on galactic and intergalactic scales. These measurements are enabled by a 0.9 sq m collecting area at 1.25 keV, a micro calorimeter array providing high-resolution spectroscopic imaging and a deployable high efficiency grating spectrometer. AXSIO delivers a 30-fold increase in effective area for high resolution spectroscopy. The key simplifications are guided by recommendations in the NWNH panel report include a reduction in focal length from 20m to 10m, eliminating the extendable optical bench, and a reduction in the instrument complement from six to two, avoiding a movable instrument platform. A focus on spectroscopic science allows the spatial resolution requirement to be relaxed to 10 arc sec (with a 5 arc sec goal). These simplifications decrease the total mission cost to under the $2B cost to NASA recommended by NWNH. AXSIO will be available to the entire astronomical community with observing allocations based on peer-review.

White, Nicholas E.↗

Negligible magnetic losses at low temperatures in liquid phase epitaxy grown Y 3 Fe 5 O 12 films

Yttrium iron garnet (Y 3 Fe 5 O 12 ; YIG) has a unique combination of low magnetic damping, high spin-wave conductivity, and insulating properties that make it a highly attractive material for a variety of applications in the fields of magnetics and spintronics. While the room-temperature magnetization dynamics of YIG have been extensively studied, there are limited reports correlating the low-temperature magnetization dynamics to the material structure or growth method. Here, in this study, we investigate liquid phase epitaxy grown YIG films and their magnetization dynamics at temperatures down to 10 K. We show there is a negligible increase in the ferromagnetic resonance linewidth down to 10 K, which is unique when compared with YIG films grown by other deposition methods. From the broadband ferromagnetic resonance measurements, polarized neutron reflectivity, and scanning transmission electron microscopy, we conclude that these liquid phase epitaxy grown films have negligible rare-earth impurities present, specifically the suppression of Gd diffusion from the Gd 3 Ga 5 O 12 (GGG) substrate into the Y 3 Fe 5 O 12 film, and therefore negligible magnetic losses attributed to the slow-relaxation mechanism. Overall, liquid phase epitaxy YIG films have a YIG/GGG interface that is five times sharper and have ten times lower ferromagnetic resonance linewidths below 50 K than comparable YIG films by other deposition methods. Thus, liquid phase epitaxy grown YIG films are ideal for low-temperature experiments/applications that require low magnetic losses, such as quantum transduction and manipulation via magnon coupling.

36 MATERIALS SCIENCE↗

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗

Real- and momentum-space description of the excitons in bulk and monolayer chromium tri-halides

Excitons with large binding energies ~2–3 eV in CrX 3 have been characterized as being localized (Frenkel) excitons that emerge from the atomic d – d transitions between the Cr-3d-t 2g and e g orbitals. The argument has gathered strength in recent years as the excitons in recently made monolayers are found at almost the same energies as the bulk. The Laporte rule, which restricts such parity forbidden atomic transitions, can relax if a symmetry-breaking mechanism is present. While what can be classified as a purely Frenkel exciton is a matter of definition, we show using an advanced first principles parameter-free approach that these excitons in CrX 3 , in both its bulk and monolayer variants, have band origin and it is the dp hybridization between Cr and X that primarily acts as the symmetry-breaking mechanism that relaxes the Laporte rule. We show that the character of these excitons is mostly determined by the Cr-d orbital manifold, nevertheless, the fractions of the spectral weight shared with the ligand halogen states increases as the dp hybridization enhances. The hybridization enhances as the halogen atom becomes heavier, bringing the X-p states closer to the Cr-d states in the sequence Cl → Br → I, with an attendant increase in exciton intensity and a decrease in binding energy. By applying a range of different kinds of perturbations that qualitatively mimics the effects originating from the missing vertex in self-energy, we show that moderate changes to the two-particle Hamiltonian that essentially modifies the Cr-d-X-p hybridization, can alter both the intensities and positions of the exciton peaks. A detailed analysis of several deep-lying excitons, with and without strain, elucidates the fact that the exciton is most Frenkel-like in CrCl 3 and CrBr 3 and acquires mixed Frenkel–Wannier character in CrI 3 , making the excitons in CrI 3 most susceptible to environmental screening and spin–orbit coupling.

36 MATERIALS SCIENCE↗

Electronic structure of boron and aluminum δ -doped layers in silicon

Recent work on atomic-precision dopant incorporation technologies has led to the creation of both boron and aluminum δ-doped layers in silicon with densities above the solid solubility limit. For this study, we use density functional theory to predict the band structure and effective mass values of such δ layers, first modeling them as ordered supercells. Structural relaxation is found to have a significant impact on the impurity band energies and effective masses of the boron layers, but not the aluminum layers. However, disorder in the δ layers is found to lead to a significant flattening of the bands in both cases. We calculate the local density of states and doping potential for these δ-doped layers, demonstrating that their influence is highly localized with spatial extents at most 4 nm. We conclude that acceptor δ-doped layers exhibit different electronic structure features dependent on both the dopant atom and spatial ordering. This suggests prospects for controlling the electronic properties of these layers if the local details of the incorporation chemistry can be fine-tuned.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning Organic Semiconductor Packing and Morphology through Non-equilibrium Solution Processing (Final Report)

Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.

36 MATERIALS SCIENCE↗

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗

Quasi-Lindblad pseudomode theory for open quantum systems

Here, we introduce a new framework to study the dynamics of open quantum systems with linearly coupled Gaussian baths. Our approach replaces the continuous bath with an auxiliary discrete set of pseudomodes with dissipative dynamics, but we further relax the complete positivity requirement in the Lindblad master equation and formulate a quasi-Lindblad pseudomode theory. We show that this quasi-Lindblad pseudomode formulation directly leads to a representation of the bath correlation function in terms of a complex weighted sum of complex exponentials, an expansion that is known to be rapidly convergent in practice and thus leads to a compact set of pseudomodes. The pseudomode representation is not unique and can differ by a gauge choice. When the global dynamics can be simulated exactly, the system dynamics is unique and independent of the specific pseudomode representation. However, the gauge choice may affect the stability of the global dynamics, and we provide an analysis of why and when the global dynamics can retain stability despite losing positivity. We showcase the performance of this formulation across various spectral densities in both bosonic and fermionic problems, finding significant improvements over conventional pseudomode formulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗