Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “intermolecular interaction”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Ion‐Specific Interactions Engender Dynamic and Tailorable Properties in Biomimetic Cationic Polyelectrolytes

Abstract Biomaterials such as spider silk and mussel byssi are fabricated by the dynamic manipulation of intra‐ and intermolecular biopolymer interactions. Organisms modulate solution parameters, such as pH and ion co‐solute concentration, to effect these processes. These biofabrication schemes provide a conceptual framework to develop new dynamic and responsive abiotic soft material systems. Towards these ends, the chemical diversity of readily available ionic compounds offers a broad palette to manipulate the physicochemical properties of polyelectrolytes via ion‐specific interactions. In this study, we show for the first time that the ion‐specific interactions of biomimetic polyelectrolytes engenders a variety of phase separation behaviors, creating dynamic thermal‐ and ion‐responsive soft matter that exhibits a spectrum of physical properties, spanning viscous fluids to viscoelastic and viscoplastic solids. These ion‐dependent characteristics are further rendered general by the merger of lysine and phenylalanine into a single, amphiphilic vinyl monomer. The unprecedented breadth, precision, and dynamicity in the reported ion‐dependent phase behaviors thus introduce a broad array of opportunities for the future development of responsive soft matter; properties that are poised to drive developments in critical areas such as chemical sensing, soft robotics, and additive manufacturing.

Aubrecht, Filip J.↗
DOE OSTI2025

Gold-Standard Chemical Database 137 (GSCDB137): A Diverse Set of Accurate Energy Differences for Assessing and Developing Density Functionals

We present GSCDB137, a rigorously curated benchmark library of 137 data sets (8377 entries) covering main-group and transition-metal reaction energies and barrier heights, (intra- and intermolecular) noncovalent interactions, dipole moments, polarizabilities, electric-field response energies, and vibrational frequencies. Legacy data from GMTKN55 and MGCDB84 have been updated to today's best reference values; redundant or low-quality points were removed, and many new, property-focused sets were added. Testing 29 popular density functional approximations (DFAs) confirms the expected Jacob's-ladder hierarchy overall but also reveals notable exceptions: functional performance for frequencies and electric-field properties correlates poorly with that for other ground-state energetics. ωB97M-V and ωB97X-V are the most balanced hybrid meta-GGA and hybrid GGA, respectively; B97M-V and revPBE-D4 lead the meta-GGA and GGA classes. Double hybrids lower mean errors by about 30% versus their hybrid analogues but demand careful frozen-core, basis set, and spin contamination treatment. GSCDB137 offers a comprehensive, openly documented platform for rigorous validation of DFA and universal machine learning potentials, and training of the next generation of exchange-correlation functionals.

Liang, Jiashu [University of California, Berkeley,↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2025

Influence of counterion substitution on the properties of imidazolium-based ionic liquid clusters

Due to their unique physiochemical properties that may be tailored for specific purposes, ionic liquids (ILs) have been investigated for various applications, including chemical separations, catalysis, energy storage, and space propulsion. The different cations and anions comprising ILs may be selected to optimize a range of desired properties, such as thermal stability, ionic conductivity, and volatility, leading to the designation of certain ILs as designer “green” solvents. The effect of counterions on the properties of ILs is of both fundamental scientific interest and technological importance. Herein, we report a systematic experimental and theoretical investigation of the size, charge, stability toward dissociation, and geometric/electronic structure of 1-ethyl-3-methyl imidazolium (EMIM)-based IL clusters containing two different atomic counterions (i.e., bromide [Br − ] and iodide [I − ]). This work extends our studies of EMIM + cations with atomic chloride (Cl − ) and molecular tetrafluoroborate (BF 4 − ) anions reported previously by Baxter et al. [Chem. Mater. 34, 2612 (2022)] and Zhang et al . [J. Phys. Chem. Lett. 11, 6844 (2020)], respectively. Distributions of anionic IL clusters were generated in the gas phase using electrospray ionization and characterized by high mass resolution mass spectrometry, energy-resolved collision-induced dissociation, and negative ion photoelectron spectroscopy experiments. The experimental results reveal anion-dependent trends in the size distribution, relative abundance, ionic charge state, stability toward dissociation, and electron binding energies of the IL clusters. Complementary global optimization theory provides molecular-level insights into the bonding and electronic structure of a selected subset of clusters, including their low energy structures and electrostatic potential maps, and how these fundamental characteristics are influenced by anion substitution. Collectively, our findings demonstrate how the fundamental properties of ILs, which determine their suitability for many applications, may be tuned by substituting counterions. These observations are critical in the sub-nanometer cluster size regime where phenomena do not scale predictably to the bulk phase, and each atom counts toward determining behavior.

cluster↗
DOE OSTI2024

Kerfeld - Dwyer Data 20240614.sf3

Mass spectrometry data for the DERA-RPI complex submission related to Supplemental Figure S3 in Dwyer et. al., Communications Biology.

Analytical biochemistry↗
DOE OSTI2026

A Solid State Zwitterionic Plastic Crystal With High Static Dielectric Constant

Solid materials with a high dielectric constant have a wide range of applications in the energy storage field. In this research, an imidazolium-based zwitterion is designed, synthesized, and confirmed to have a plastic crystal phase based on the following experimental and computational evidence: (i) the presence of long-range order with weak intermolecular forces and competing attractive-repulsive interactions along different crystallographic directions; (ii) the observation of more than one endotherm on heating including a solid-solid phase transition at T s-s = −26 °C and melting of the plastic crystal at T m = 72°C; (iii) a low entropy of fusion at melting (2.1 JK −1 mol −1 ); (iv) a strongly anisotropic morphology; (v) relatively fast dynamics originating from short-range degrees of freedom. Furthermore, it exhibits a very high dielectric constant in the plastic crystal solid state (147 at −10°C and 103 at 70°C) due to the rotational degrees of freedom of plastic crystals that arise from weak net intermolecular interactions of zwitterions due to only having two carbons between the anion and cation. This material conveniently remains in the plastic crystal phase within 50 K of ambient temperature. This discovery opens new opportunities in the search for solid-state high dielectric constant materials.

dielectric relaxation↗
DOE OSTI2026

Molecular To Mesoscale Targeting of Oxoanions with Multi-Tasking Hosts

Achieving a better understanding of anion interactions both in solution and crystalline state was the overarching goal of this project. Anions are everywhere throughout Nature and play important roles in biological and environmental processes. They can be beneficial or deleterious or both in different situations and concentrations. For either reason it is important to have molecules that can bind anions for key needs that benefit society. However, recognition of specific anions is challenging due to the diffuse nature of their negative charge(s) as well as their various shapes and sizes. Understanding the basic properties of anions and how they interact with other molecules and ions in surrounding environments is key to selective recognition. In this project multi-tasking molecules for selective binding of targeted anions were designed to achieve cooperativity and synergism in one rather than multiple host molecules, including (1) cation:anion pair hosts for anions with charges of -2 or greater; (2) pH and redox activated hosts for on-off binding and release; and (3) multiple anion capture in extended host networks. Our design strategy was to combine the use of simple inexpensive building blocks and high yield synthetic pathways to provide economically feasible scale-up for applications. Oxoanions representing multiple shapes and charges were chosen based on having the potential for significant impact on DOE separations needs. Amide/amine-based macrocycles and urea/amine-based chelates and macrocycles with multiple hydrogen bonding sites provided the basic anion-binding frameworks. Successful multi-tasking outcomes were forthcoming in all three tasks. In Task 1, successful ion pair binding for anions with multiple charges was achieved. Furthermore, the ion pair molecules were capable of extended interactions through supramolecular intertwining, like fishing nets for capturing pools of fish (also fitting with Task 3). In Task 2, molecules were synthesized possessing on-off switches. These included a pH sensitive sensor for on-off binding of anions in general, as well as an electrochemical sensor selective for sulfate capture. Three new classes of extended anion host networks capable of binding multiple ions was a major outcome of Task 3. These systems included: anion sensitive, fluorescent organogels; channel-forming macrocycles for studying anion-water including larger macrocyclic cluster sandwiches; and, the offshoot of Task 1, fishing net ion-pair networks for higher valent anions. These strategies can be expanded in the future to other ions and molecules for a better understanding of intermolecular and interionic interactions.

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

Atomic-Scale Imaging Reveals Polar-π Interactions in Two-Dimensional Molecular Superlattices

Controlling coassembly of synthetic oligomers into binary superlattices at the atomic level is challenging. Here, we report a strategy for programming polar-π interactions in oligomeric peptoids, a class of sequence-defined peptidomimetics, facilitating the formation of homogeneous two-dimensional (2D) superlattices. N-2-phenylethyl and N-(2-perfluorophenyl)ethyl side chains, similar in size, but with contrasting electrostatic characteristics, were introduced at defined sequence positions to generate favorable dipolar aromatic interactions. The resulting nanosheets exhibit different crystal motifs depending on the side chain interactions: systems containing only one type of aromatic side chain form a parallel V-shaped motif driven by π-π interactions, whereas a combination of both types of aromatic side chains, either within one backbone or through the coassembly of two distinct peptoids, adopt an antiparallel V-shaped superlattice with higher thermal stability, driven by polar-π interactions. Cryogenic transmission electron microscopy directly resolved the packing arrangement of perfluorophenyl and phenyl rings in individual nanosheet superlattices, confirming that intermolecular polar-π interaction dominates the superlattice motifs and increases lattice stability. Molecular dynamics simulations and density functional theory calculations further substantiate the energetic favorability of polar-π interactions over π-π interactions, rationalizing the formation of homogeneous superlattices with enhanced thermal stability. Our discoveries establish a design principle for binary coassembly using sequence-defined oligomers, which enables control over unit cell geometry, lattice stability, and molecular registration through aromatic side chain polarization and sequence control. This ability to program atomic-scale binary superlattices opens new avenues for designing functional 2D soft materials.

Lee, Yen Jea [Lawrence Berkeley National Laborator↗
DOE OSTI2024

Precision Synthesis of Conjugated Polymer Films by Surface-Confined Stepwise Sonogashira Cross-Coupling

Thin films of poly(arylene ethynylene)-conjugated polymers, including low-energy-gap donor–acceptor polymers, can be prepared via stepwise polymerization utilizing surface-confined Sonogashira cross-coupling. This robust and efficient polymerization protocol yields conjugated polymers with a precise molecular structure and with nanometer-level control of the organization and the uniform alignment of the macromolecular chains in the densely packed film. In addition to high stability and predictable and well-defined molecular organization and morphology, the surface-confined conjugated polymer chains experience significant interchain electronic interactions, resulting in dominating intermolecular π-electron delocalization which is primarily responsible for the electronic and spectroscopic properties of polymer films. The fluorescent films demonstrate remarkable performance in chemosensing applications, showing a turn-off fluorescent response on the sub-ppt (part per trillion) level of nitroaromatic explosives in water. This unique sensitivity is likely related to the enhanced exciton mobility in the uniformly aligned and structurally monodisperse polymer films.

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

Hydrophobic Metal–Organic Frameworks Enable Superior High-Pressure Ammonia Storage through Geometric Design

Hydrophobic metal–organic frameworks (MOFs) are typically overlooked for ammonia storage due to weak host–guest interactions. Here, we demonstrate that four structurally analogous aluminum-based MOFs exhibit a counterintuitive behavior whereby framework geometry, rather than ligand hydrophilicity, determines high-pressure NH 3 adsorption performance. The hydrophobic CAU-23 achieved an exceptional capacity matching hydrophilic analogs despite its poor low-pressure uptake. This pressure-dependent enhancement stems from the unique 4-cis-4-trans geometry of CAU-23 compared to the purely cis arrangement of MIL-160 and KMF-1 and the alternating cis-trans configuration of MOF-303. Critically, CAU-23 retained 95% capacity over three high-pressure cycles, whereas hydrophilic MOFs suffered 39–46% irreversible losses due to strong NH 3 -framework interactions that compromise structural integrity. Grand canonical Monte Carlo simulations reveal that high pressure enables NH 3 clustering through intermolecular hydrogen bonding, bypassing the need for strong host–guest interactions. High-pressure powder X-ray diffraction measurements confirm the exceptional mechanical resilience of CAU-23, showing complete structural recovery upon decompression despite exhibiting the highest pressure sensitivity among the studied MOFs. An extended analog, HE-CAU-23, validates this design principle with further enhanced capacity. Furthermore, these findings reveal a paradigm shift toward hydrophobic MOFs with optimized geometry for high-performance and regenerable gas storage applications.

MOFs↗
DOE OSTI2024

Self-assembly and condensation of intermolecular poly(UG) RNA quadruplexes

Abstract Poly(UG) or ‘pUG’ dinucleotide repeats are highly abundant sequences in eukaryotic RNAs. In Caenorhabditis elegans, pUGs are added to RNA 3′ ends to direct gene silencing within Mutator foci, a germ granule condensate. Here, we show that pUG RNAs efficiently self-assemble into gel condensates through quadruplex (G4) interactions. Short pUG sequences form right-handed intermolecular G4s (pUG G4s), while longer pUGs form left-handed intramolecular G4s (pUG folds). We determined a 1.05 Å crystal structure of an intermolecular pUG G4, which reveals an eight stranded G4 dimer involving 48 nucleotides, 7 different G and U quartet conformations, 7 coordinated potassium ions, 8 sodium ions and a buried water molecule. A comparison of the intermolecular pUG G4 and intramolecular pUG fold structures provides insights into the molecular basis for G4 handedness and illustrates how a simple dinucleotide repeat sequence can form complex structures with diverse topologies.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2024

Structural basis for inhibition of coagulation factor VIII reveals a shared antigenic hotspot on the C1 domain

Hemophilia A arises from dysfunctional or deficient coagulation factor (F)VIII and leads to inefficient fibrin clot formation and uncontrolled bleeding events. The development of antibody inhibitors is a clinical complication in hemophilia A patients receiving FVIII replacement therapy. LE2E9 is an anti-C1 domain inhibitor previously isolated from a mild/moderate hemophilia A patient and disrupts FVIII interactions with von Willebrand factor and FIXa, though the intermolecular contacts that underpin LE2E9-mediated FVIII neutralization are undefined. To determine the structure of the complex between FVIII and LE2E9 and characterize its mechanism of inhibition. FVIII was bound to the antigen binding fragment (Fab) of NB2E9, a recombinant construct of LE2E9, and its structure was determined by cryogenic electron microscopy. Here, this report communicates the 3.46 Å structure of FVIII bound to NB2E9, with its epitope comprising FVIII residues S2040 to Y2043, K2065 to W2070, and R2150 to H2155. Structural analysis reveals that the LE2E9 epitope overlaps with portions of the epitope for 2A9, a murine-derived inhibitor, suggesting that these residues represent a shared antigenic region on the C1 domain between FVIII –/– mice and hemophilia A patients. Furthermore, the FVIII:NB2E9 structure elucidates the orientation of the LE2E9 glycan, illustrating how the glycan sterically blocks interactions between the FVIII C1 domain and the von Willebrand factor D' domain. A putative model of the FVIIIa:FIXa complex suggests potential clashing between the NB2E9 glycan and FIXa light chain. These results describe an antigenic “hotspot” on the FVIII C1 domain and provide a structural basis for engineering FVIII replacement therapeutics with reduced antigenicity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

How alkyl branching shapes structure in imidazolium and pyrrolidinium NTf 2 ionic liquids

High-energy X-ray scattering experiments and molecular dynamics (MD) simulations were carried out on ionic liquids (ILs) consisting of 1-alkyl-3-methylimidazolium and 1-alkyl-1-methylpyrrolidinium cations. These cations were paired with bis(trifluoromethylsulfonyl)amide anions and identical alkyl tails were used for both cationic species. The goal of this work is to investigate how the nanoscale structure of the ionic liquid changes with the length and with the degree of branching of the alkyl tail, for ILs having a common anion. We investigate spatial correlations in the intermolecular region, focusing on the intrinsic charge-charge interactions that characterize all ionic liquids, as well as the nanoscale domain segregation that is present in IL species with significant nonpolar components.

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

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Predicting compatibility between ferredoxins and the Fe protein of nitrogenase using in silico protein modeling

Biological nitrogen fixation is the process by which certain bacteria and archaea use the enzyme nitrogenase to reduce atmospheric nitrogen into bioavailable ammonium. Engineering non‐nitrogen‐fixing organisms, like plants, to use nitrogenase could reduce dependency on synthetic fertilizer and mitigate the environmental impacts of industrial fertilizer production. However, nitrogenase activity requires delivery of reducing power by small electron carrying proteins known as ferredoxins and flavodoxins, and successfully engineering nitrogenase into new systems will require a mechanistic understanding of electron delivery by these proteins. Most organisms often have multiple ferredoxins, raising the question of which ferredoxin can support nitrogenase activity. The purpose of this study is to gain insight into how we can predict which ferredoxin is compatible with the Fe protein, the component of nitrogenase that interacts with ferredoxin or flavodoxin. Our in silico protein–protein docking simulations reveal that most ferredoxins and flavodoxins involved in nitrogen fixation have the shortest distance (≤10 Å) between their redox cofactor and the [4Fe‐4S] cluster of the Fe protein. We found shorter cofactor distance contributes to faster intermolecular electron tunneling rates. Bacterial ferredoxins that play a role in nitrogen fixation also exhibit more complementary interactions with the Fe protein than bacterial and plant ferredoxins not involved in this process. Heterologous expression of a set of ferredoxins from both nitrogen‐fixing and non‐nitrogen‐fixing bacteria in the diazotroph Rhodopseudomonas palustris supports our model‐derived prediction that shorter distances between the electron‐carrying cofactors favor nitrogenase compatibility. These findings offer a framework to predict and potentially enhance ferredoxin–nitrogenase compatibility, which will help to improve our ability to engineer nitrogen fixation into non‐nitrogen‐fixing organisms like plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Distinct Mechanisms of Triplet Pair Decay in Amorphous and Crystalline Heteroacene Thin Films

Triplet pairs (TT) in crystalline molecular semiconductors have unique spin properties of interest for quantum information or enhancing solar photoconversion. The population and diffusion dynamics of TT have been the subject of recent studies, both in covalent dimers and in crystalline systems. Here, in this study, we monitor the triplet population in neat polycrystalline and amorphous films of a heteroacene with known TT spectral properties and tunable spin polarization depending on the intermolecular geometry. Transient measurements reveal an anomalous power dependence in polycrystalline films that we attribute to the fast diffusion and interaction of dissociated triplet pairs confined to one-dimensional stacks of strongly coupled molecules. The nongeminate triplet interaction after dephasing facilitates conversion to the triplet 3 TT and eventually T 1 +S 0 . Amorphous films have no power dependence and proceed directly from 1 TT to 3 TT and subsequently T 1 +S 0 via state mixing facilitated by nonparallel geometries and weak exchange coupling.

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

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗