Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal properties”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Spontaneous Twirls and Structural Frustration in Moiré Materials

Structural twirls form spontaneously in the domain wall networks of some moiré materials. We show that in heterobilayers, neighboring twirl chiralities tend to antialign, forming staggered patterns that are well described by antiferromagnetic lattice 𝜙 4 theories. In moiré systems with triangular domains, this leads to frustration in the chirality configuration of the twirls and to hysteresis with respect to variation of the average twist angle and possibly other control parameters. As a result, we expect that in typical materials, the ordering temperature of twirls is about 10 3 K, and that thermal fluctuations in individual twirl chiralities freeze below room temperature.

Elastic deformation↗
DOE OSTI2025

TiSe 2 is a band insulator created by lattice fluctuations, not an excitonic insulator

TiSe 2 is a narrow-gap insulator with a rich array of unique properties. In addition to being a superconductor under certain modifications, it is commonly thought to be a rare realisation of an excitonic insulator. Below 200 K, TiSe 2 undergoes a transition from a high-symmetry ($P\bar{3}m1$) phase to a low-symmetry ($P\bar{3}c1$) charge density wave (CDW). Here we establish that it is indeed an insulator in both $P\bar{3}m1$ and $P\bar{3}c1$ phases. However, the insulating state is driven not by excitonic effects but by symmetry-breaking. In the CDW phase it is static. At high temperature, thermally driven instantaneous deviations from $P\bar{3}m1$ break the symmetry on the characteristic time scale of a phonon. Even though the time-averaged lattice structure assumes $P\bar{3}m1$ symmetry, the time-averaged energy band structure is closer to the CDW phase - a rare instance of a metal-insulator transition induced by dynamical symmetry breaking. We establish these conclusions from quasiparticle self-consistent GW (QSGW) and many-body calculations (QS$G\widehat{W}$), in combination with molecular dynamics simulations to capture the effects of thermal disorder. The many-body theory includes explicitly ladder diagrams in the polarizability, which incorporates excitonic effects in an ab initio manner. We find that the excitonic modification to the potential is weak, ruling out the possibility that TiSe 2 is an excitonic insulator.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fatty acid foams for nonselective physical removal of microplastics from aqueous solutions

Microplastics (MPs) are pervasive environmental contaminants whose removal from water remains a major challenge due to their small size, chemical diversity, and dynamic surface properties arising from environmental aging/weathering. Here, we present a concept of foam-based separation method that physically traps MPs in the foam phase using microtubular assemblies of 12-hydroxystearic acid. These foams are stabilized by anisotropic fatty acid microtubules formed in the presence of ethanolamine, which jam within the foam channels and suppress fluid drainage thereby enhancing MP retention and foam stability. MPs of different sizes, polymer compositions (including polystyrene, polypropylene, polyethylene terephthalate, and polytetrafluoroethylene), and weathered states were retained in the foam phase without requiring chemical modification or relying on chemical interactions between the fatty acid and MPs. Thermally induced transition of the fatty acid microtubules into nanomicelles above the characteristic phase transition temperature (∼35 °C) enables controlled foam collapse and recovery of trapped MPs. The cumulative removal efficiency can exceed 85% through multiple foaming cycles, matching predictions from a probabilistic retention model. This work shows that foams can provide a simple platform to trap MPs, thus providing a new physical-removal strategy that does not rely on the particles’ chemistry.

Guillot, Kennedy A. [Louisiana State Univ., Baton ↗
DOE OSTI2025

Investigation of Oxidation Mechanisms in HfS 2 and ZrS 2 via In Situ Electron Microscopy

Transition metal dichalcogenides (TMD) combine semiconducting properties with stability and are therefore widely studied for future microelectronic devices. However, applications require integration with a suitable dielectric at a planar, defect-free interface. For Hf- and Zr-based TMDs, the good dielectric properties of hafnia and zirconia imply that direct oxidation is a promising strategy, provided that a high-quality interface and suitable oxide morphology can be formed. Here, in this study, we investigate HfS 2 and ZrS 2 oxidation, aiming to understand pathways toward planar oxide layer formation and the mechanisms that determine the characteristics of the semiconductor/dielectric interface. We use conventional and environmental transmission electron microscopy to reveal changes in morphology and composition of the TMDs under different gaseous environments. We show that oxidation at ambient conditions causes compositional heterogeneities, with local replacement of S by O in both materials. Thermal oxidation causes desulfurization and formation of a smooth but defective oxide layer. In contrast, plasma oxidation appears the best-suited to forming a controlled TMD/dielectric interface with a smooth oxide layer and without major defects. These results provide insight into the opportunities available from HfS 2 and ZrS 2 through oxidation and suggest a materials processing strategy for electronic device fabrication.

HfS2↗
DOE OSTI2026

FAST-1.2.2: A Computer Code for Thermal-Mechanical Nuclear Fuel Analysis under Steady-state and Transients

Fuel Analysis under Steady-state and Transients (FAST) is the U.S. Nuclear Regulatory Commission (NRC)’s computer code that calculates the steady-state and transient response of nuclear reactor fuel rods during long-term in-reactor burnup, anticipated operational occurrences (AOOs), design basis accidents (DBAs), and dry storage conditions. The code calculates the temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include heat conduction through the fuel and other materials, heat transfer from the cladding-to-coolant, cladding elastic and plastic deformation (including creep), fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, void volume, and cladding oxidation. The code contains necessary material and coolant properties, as well as clad-to-coolant heat transfer correlations, for normal operation through postulated accidents and AOOs for today’s U.S.-based light water reactor (LWR) fuel designs. FAST-1.2.2 also contains preliminary materials and models for new LWR fuel concepts, such as accident tolerant fuel (ATF), and non-LWR fuel concepts such as metallic fuels for sodium fast reactors (SFRs). FAST has been developed for use on Windows and Linux operating systems. This document describes FAST-1.2.2 and is one of a series of documents on the code; the other documents detail the material properties used by FAST as well as its integral assessment to experiments and commercial data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Nondestructive ultrasonic characterization of a triple-weld-bead wire arc additively manufactured ER70S-6 S-curved wall

Understanding build-scale microstructural variation in wire arc additive manufacturing (WAAM) of low-carbon steels is essential for ensuring consistent structure–property relationships throughout large components. Conventional destructive characterization techniques, such as scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), are time-intensive and limited to localized regions, making comprehensive evaluation of large WAAM structures challenging. In this study, a nondestructive ultrasonic approach was employed to characterize a 252 mm tall ER70S-6 S-curved WAAM wall produced using a triple-bead deposition strategy. Optical and SEM analysis revealed a repeating dual-region microstructure consisting of uniform polygonal ferrite at melt pool centers and heterogeneous ferrite with coarse and fine grains near melt pool boundaries, attributed to cyclic thermal conditions. Longitudinal ultrasonic backscatter imaging was used to evaluate the continuity of this periodicity along the full build height. The ultrasonic response exhibited a consistent repeating pattern that correlated with the observed layer-wise microstructural variation. X-ray computed tomography confirmed the absence of detectable porosity, indicating that ultrasonic contrast is primarily governed by grain morphology. Overall, the results indicate that longitudinal backscatter ultrasound is a promising nondestructive characterization technique for validating microstructural variations along the s-curved WAAM wall, with significant potential for microstructure optimization and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molten Salt Corrosion Tests of Additively Manufactured Stainless Steel 316H

Molten salt reactors (MSRs) have drawn considerable interest due to their favorable safety features, high thermal efficiency, and compatibility with different fuel cycles. Yet, the success of MSRs hinges critically on the performance of structural materials to be used in these aggressive molten salt environments, where corrosion and material compatibility remain primary challenges to long-term reliability. Additively manufactured (AM) nuclear structural materials prompt the use of novel geometries and compositions to enhance material performance and reduce costs of constructing MSRs. The rapid solidification conditions inherent to AM processing impart distinctive microstructural features, including cellular sub-structures, dislocation densities, residual stress, and oxide inclusions, which can influence material performance in MSR components. While the mechanical properties of AM stainless steels have been widely studied, their corrosion behavior, particularly in molten salt environments, has received far less attention. Addressing these needs, the Advanced Materials and Manufacturing Technologies (AMMT) program provides a framework for systematically evaluating how unique microstructures produced by AM processes influence the performance of these materials in these demanding environments and for developing reproducible testing workflows that can support future code qualification efforts and standards development. Bridging this knowledge gap is essential for assessing the viability of AM alloys in MSRs and informing qualification strategies. A further challenge is the absence of standardized protocols for molten salt corrosion testing. Accordingly, this report provides an account of the corrosion evaluation of AM 316H stainless steel in NaCl 2 -MgCl 2 molten salt at 550 °C, with exposure times of 100 and 500 hours. It documents the experimental procedures implemented under the AMMT program, including salt preparation, exposure protocols, and post-test characterization methods, to establish reproducibility and transparency. Importantly, the study examines AM 316H samples in the as-fabricated condition, directly reflecting the surface state most relevant to engineering applications, and compares their behavior to machine-cut surfaces. Overall, preliminary evaluations have noted that surface conditions (e.g. morphology, contamination, etc.) have a noticeable impact on the corrosion resiliency. The impact of the corrosion is difficult to detect at 100h, unless, in the case of AM 316H, the specimen surface is decontaminated. After 500 h, as-fabricated surfaces of AM and wrought 316H display evidence of general versus preferential corrosion attack, respectively. Both AM and wrought 316H machine-cut surfaces exhibit a continuous Cr depletion zone, evident of general corrosion. While the estimated extent of corrosion appears within the same order of magnitude regardless of the surface condition, it is apparent that more predictable behavior is observed on machine-cut surfaces. Nonetheless, further investigation is necessary to fully elucidate the corrosion mechanism under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Polymorphs of the n–Type Polymer P(NDI2OD–T2): A Comprehensive Description of the Impact of Processing on Crystalline Morphology and Charge Transport

A systematic study of the polymorphs emerging in P(NDI2OD-T2) (also commercially known as N2200), a prototypical organic semiconducting n-type polymer, is presented. Using a tightly integrated experimental and computational approach, detailed atomistic-level descriptions are provided investigating the three known P(NDI2OD-T2) polymorphs observed at room temperature as a function of thin-film processing. Importantly, over the course of the work, a missing link is uncovered, a fourth polymorph referred to here as Form I-β; this new form is a morphological intermediary observed upon thermal annealing, which evolves from Form I but tends to disappear upon full polymer chain melting. The computationally derived polymorph structures show excellent agreement with experimental X-ray scattering characterization. The relative stabilities of each polymorph are calculated in terms of both the bulk material and the polymorph-air interface. An energy landscape is then constructed to qualitatively compare the thermodynamic versus kinetic origins of each polymorph, and the factors driving (supra)assembly and associated transformations among polymorphs using an approach generalizable to other organic semiconducting polymers. Lastly, the relationships among preferential polymorphic crystallinity, relative chain orientations, and directional charge transport properties in P(NDI2OD-T2) are explored. Altogether, this work provides unprecedented insights into complex structure-processing-transport relationships in a representative semiconducting organic polymer.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis and Characterization of Methacrylamide-Based Block Random Copolymers via Amine Functionalization of Polystyrene- Block -Poly(Pentafluorophenyl Methacrylate) Toward Enhanced Amenability to Manufacturing Criteria in Nanolithography

Recent interest in manipulating the chemistry of block copolymers (BCPs) to manage the covarying properties necessary to meet manufacturing criteria in nanolithography has resulted in the development and expansion of the A-block-(B-random-C) BCP architecture, where the random block allows for the decoupling of thermodynamic and wetting properties. Previous reports of such BCPs have used click chemistry to create the desired random block, but all such instances possess additional functional groups that are susceptible to undesirable side reactions upon annealing, such as cross-linking and surface grafting. This study reports the substitution reaction of polystyrene-block-poly(pentafluorophenyl methacrylate) (PS-b-PPFMA) with primary amines. The resulting methacrylamide structure of the functionalized random block has only an amide linkage between the attached functional group and the polymer backbone, thereby omitting any sources of side reactions within the BCP. The outcomes of this report also demonstrate the enhanced thermal stability of these materials by virtue of their stable amide bonds. A range of random copolymer compositions is assessed to develop BCPs in which the random block has approximately the same surface energy as the PS block. These BCPs can self-assemble into perpendicular lamellae and are therefore promising candidates for directed self-assembly.

block copolymers↗
DOE OSTI2025

A heat and mass transfer model for evaluation of damaged cesium chloride radiological sources

Radiological sources are vital to many applications, and typically contain byproducts from waste streams which decay over time to form complex mixtures of elements. Daughter products resulting from these decay processes can introduce complicating factors to the integrity and safety of the vessel that contains the radiological material. Sources containing cesium 137 ( 137 Cs) are of particular interest, because the decay of this isotope produces barium metal which reacts readily and exothermically with oxygen. This work employs physics-based numerical tools to examine the thermal response of a radiological source containing a mixture of cesium and barium in the event that vessel walls are damaged and atmospheric gases contact source material. A parametric study was conducted to determine the sensitivity of the response to various factors, including vessel geometry, source material age, the degree of vessel damage, and other parameters. Here, it was found that the peak temperatures that occurred within the source material strongly depend on these parameters, particularly vessel geometry and age, which determine whether a breach would be a relatively minor accident or a catastrophic incident. Finally, the model’s sensitivity to uncertain thermophysical properties is discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Multiscale structural analysis of polymorphic phase boundaries in doped antiferroelectric sodium niobate

In the current work, we have performed multiscale structural analysis of a Pb-free sodium niobate-based smart system, i.e., 0.9⁢NaNbO 3 –0.1⁢Ba 0.9 ⁢Ca 0.1 ⁢TiO 3 (NN-10BCT) reported earlier for its high ferroelectric response. We have investigated the temperature-dependent evolution of crystal structure at long, medium, and short ranges using synchrotron x-ray diffraction (SXRD), Raman scattering, and pair distribution function(PDF) techniques in conjunction with dielectric studies. Temperature-dependent synchrotron x-ray diffraction data combined with dielectric analysis suggest two unique polymorphic phase boundaries (PPB) with two coexisting ferroelectric phases stable in the wide temperature ranges. These PPBs are stable in different regions viz. (i) cryogenic temperatures with coexisting R3c and Pmc⁢2 1 phases (ii) vicinity of room temperatures with coexisting Pmc⁢2 1 and Amm2 phases. In contrast to the conclusions drawn from SXRD, PDF reveals structures having lower symmetry (with coexisting Cc+ Pmc2 1 phases at 1.7 ≤r≤ 20 Å) at short ranges for these PPBs. In conclusion, the presence of different long- and short-range symmetries (accommodating tilt-oriented ferroelectric phases) in the unique polymorphic phase boundaries makes them thermally stable and advantageous for technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2025

Deep potential molecular dynamics simulations of ion-enhanced etching of silicon by atomic chlorine

The continued development of plasma-assisted processing techniques requires a fundamental understanding of plasma-surface interactions. Molecular dynamics (MD) simulations have been employed to complement experimental studies and better understand the properties of such systems. Recently, machine learning (ML) methods have enabled the development of ab initio-based interatomic potentials, which can be generalized to complex combinations of multiple atom types. In this work, we use ML potentials developed using the Deep Potential Molecular Dynamics (DeepMD) framework to provide a model of ion-enhanced etching of Si by Cl atoms. We demonstrate the importance of proper selection of the training data set to the accuracy of the DeepMD model and compare our results to MD results using empirical potentials, as well as to experimental measurements. Exposure of undoped Si at 300 K to thermal Cl atoms yields a steady-state Cl coverage of 1.25 monolayers, which is slightly lower than the value obtained in previous experimental studies. Predictions of Si etch yields by simultaneous Cl atom and Ar + ion impacts as a function of ion energy, neutral to ion flux ratio, and angle of incidence of the ions are in reasonably good agreement with classical MD results and experimental measurements. Finally, etch yields and SiCl x mixed layer thicknesses during simultaneous bombardment of the Si(100) surface by Cl atoms and Cl + ions are in good agreement with experimental data. In conclusion, the present work is a necessary condition for the extension of the DeepMD procedure to more complex systems of interest in plasma-surface interactions.

Artificial neural networks↗
DOE OSTI2025

Intrinsic property of defective β-Ga 2 O 3 to self-heal under ionizing irradiation

Damage evolution and phase stability in defective β-Ga 2 O 3 and an irradiation-converted γ-Ga 2 O 3 layer have been studied under ionizing irradiation at 300 K. By exploring athermal nonequilibrium processes in β-Ga 2 O 3 , we succeed in identifying a self-healing mechanism that enables recovery pre-existing damage, characterized by a recovery cross-section of ~0.17 nm 2 . Remarkably, this study further demonstrates that the crystallinity of the irradiation-converted γ-Ga 2 O 3 layer improves under ionizing irradiation. More importantly, X-ray diffraction analysis reveals that the highly-strained 𝛾 -phase transforms into a highly-crystalline structure without film disintegration, contrasting to that reported for isochronal annealing at 1000 K. The inelastic thermal spike calculations provide insights into the important effects of energy transfer to electrons in reordering the local atomic arrangement of both defective β- Ga 2 O 3 and 𝛾-Ga 2 O 3 . This behavior suggests a pathway for low-temperature crystallization, offering a promising strategy for fabricating ultrahigh-speed non-volatile memory devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reducing heat load density with asymmetric and inclined double-crystal monochromators: principles and requirements revisited

Asymmetric double-crystal monochromators (aDCMs) and inclined DCMs (iDCMs) can significantly expand the X-ray beam footprint and consequently reduce the heat load density and gradient. Based on rigorous dynamical theory calculations, the major principles and properties of aDCMs and iDCMs are presented to guide their design and development, particularly for fourth-generation synchrotrons. In addition to the large beam footprint, aDCMs have very large bandwidths (up to ∼10 eV) and angular acceptance, but the narrow angular acceptance of the second crystal requires precise control of the relative orientations and strains. Based on Fourier coupled-wave diffraction theory calculations, it is rigorously proved that the iDCM has almost the same properties as the conventional symmetric DCM, including the efficiency, angular acceptance, bandwidth, tuning energy range and sensitivity to misalignment. The exception is that, for the extremely inclined geometry that can achieve very large footprint expansion, the iDCM has (beneficially) a larger bandwidth and wider angular acceptance. Inclined diffraction has the `rho-kick effect' that can be cancelled by the second reflection of the iDCM (even with misalignment), except that inhomogeneous strains may cause non-uniform rho-kick angles. At present, fabrication/mounting-induced strains pose low risk since they can be controlled to <0.5 µrad over large areas. The only uncertain challenge is the thermally induced strains, yet it is estimated that these strains are naturally lowered by the large footprint and may be further mitigated by optimized cryogenic cooling to the 1–2 µrad level. Overall, aDCMs and iDCMs have more stringent requirements than normal DCMs, but they are feasible schemes in practice.

asymmetric monochromator↗
DOE OSTI2025

In Situ Atomic Tracking on the Interfacial Etching and Reconfiguration of Cu-ReSe 2 Contact during Thermal Annealing

The Schottky barrier height can be greatly affected by the metal diffusion, reaction, and covalent bonding formation at the contact. Exploring novel methods and revealing the fundamental mechanisms for contact engineering are of vital importance for microelectronic devices. Here, in this study, the annealing induced interfacial reactions at Cu-ReSe 2 contact are dynamically revealed from the atomic scale. Accompanied by the diffusion of Se to Cu, ReSe 2 is gradually decomposed to a thin Re interlayer through a “chain-by-chain” manner. Theoretical calculations show that the Cu atoms can facilitate the chemical bond breaking of ReSe 2 , significantly lowering the Se diffusion energy barrier toward Cu. The formed Re/ReSe 2 heterostructure presents a metal-like band structure, which underscores the critical role of Cu in altering the interfacial chemistry and promoting carrier transport across the interface. Our results can provide vital insights into the contact properties of ReSe 2 and provide a possible method for fabricating high-performance ReSe 2 -based devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Plasmonic Ag nanocomposite phosphate glasses produced via γ-ray irradiation as reduction route

This paper reports on the impact of γ-ray irradiation (10, 100 kGy) on melt-quenched Ag + -doped phosphate glass and the effects of subsequent thermal processing leading to the production of plasmonic Ag nanocomposites. The γ-irradiated glasses were characterized alongside the pristine by differential scanning calorimetry (DSC), Raman spectroscopy, electron paramagnetic resonance (EPR) spectroscopy, optical absorption, and photoluminescence (PL) spectroscopy. DSC characterization showed consistent glass transition temperatures (T g ) before and after γ-irradiation whereas the crystallization temperatures tended to decrease with increasing γ-ray dose. However, a lack of alteration of the glass network structure was supported by Raman spectroscopy. Room temperature EPR spectra clearly showed the formation of phosphorus oxygen hole center (POHC) defects in the undoped host, in addition to another doublet likely associated with a P 3 defect. The presence of paramagnetic silver species encompassing Ag 2+ and 107/109 Ag 0 atoms was also indicated in the silver-activated glass together with POHC defects. Optical absorption spectra were also consistent with the presence of various radiation-induced centers. Further analyzing the glass absorption edge via Tauc plots suggested the formation of electron center (EC) defects in γ-irradiated samples wherein the silver-doped glass exhibited decreasing band gap energies with increasing γ-ray dose. The PL characterization showed the silver-related radio-PL effect was induced exhibiting broad band emission with two maxima around 500 and 625 nm stemming from various molecular Ag$^{x+}_{n}$ clusters. Emission decay analyses revealed that the longer wavelength emission exhibited slower decay. The highest radiation dose of 100 kGy however resulted in weaker emission and faster decay kinetics attributed to energy transfer between the luminescent silver species and POHC defects. Finally subjecting the γ-irradiated Ag-doped glasses to heat treatment near the T g at 490 °C led to the development of the surface plasmon resonance of Ag nanoparticles (NPs) and the vanishing of the Ag$^{x+}_{n}$ clusters luminescence. In conclusion, the presence of the matrix-related EC defects was deemed accountable for the thermally induced reduction and consequent precipitation of Ag NPs making the plasmonic glasses attractive for photonic applications such as nonlinear optics.

36 MATERIALS SCIENCE↗
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↗