SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Understanding Novel Lewis Acid Doping Mechanisms in Organic Semiconductors

Recent studies supported by this DOE grant have provided a more detailed understanding of the doping mechanism of organic semiconductors by the Lewis acid, tris(pentafluorophenyl)borane (BCF). Our studies have examined the intermolecular semiconductor-dopant interactions as well as the thermodynamic feasibility of different doping mechanisms. Most notably, we found a different proposed doping scheme is much more thermodynamically favoured (exergonic) than our proposed two-step BCF doping mechanism. The state-of-the-art experimental and theoretical techniques of High-Field Multi-Dimensional Solid State NMR spectroscopy and Time-Dependent Density Functional theory (TD-DFT) and DFT were employed in this work. The data obtained helped identify the structural details and intermolecular interactions which impact the organic semiconductor doping efficiency for BCF and how these compare to the case of doping by the standard integer charge transfer dopant molecule F 4 TCNQ. We have, for the first time, assigned specific chemical shift values for the F 4 TCNQ molecules clustered outside of polymer chains as well as those intercalated between the polymer backbones. Further, we provide theoretical support for a novel doping mechanism which combines those proposed in previous literature involving the evolution of hydrogen gas along with a doping species more complex than the simple BCF·H 2 O. While our previously proposed doping mechanism was found to be highly thermodynamically unfavorable, the formation of large anion [BCF(OH)(OH 2 )BCF]- is the key change making the proposed mechanism exergonic while the elimination of gaseous H 2 helps drive the doping reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

A hybrid Monte Carlo, discontinuous Galerkin method for linear kinetic transport equations

Here we present a hybrid method for time-dependent particle transport problems that combines Monte Carlo (MC) estimation with deterministic solutions based on discrete ordinates. For spatial discretizations, the MC algorithm computes a piecewise constant solution and the discrete ordinates use bilinear discontinuous finite elements. From the hybridization of the problem, the resulting problem solved by Monte Carlo is scattering free, resulting in a simple, efficient solution procedure. Between time steps, we use a projection approach to “relabel” collided particles as uncollided particles. In conclusion, from a series of standard 2-D Cartesian test problems we observe that our hybrid method has improved accuracy and reduction in computational complexity of approximately an order of magnitude relative to standard discrete ordinates solutions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Orbitally selective resonant photodoping to enhance superconductivity

Signatures of superconductivity at elevated temperatures above T c in high-temperature superconductors have been observed near 1/8 hole doping for photoexcitation with infrared or optical light polarized either in the CuO 2 plane or along the c axis. While the use of in-plane polarization has been effective for incident energies aligned to specific phonons, c-axis laser excitation in a broad range between 5 μm and 400 nm was found to affect the superconducting dynamics in striped La 1.88 5Ba 0.115 CuO 4 , with a maximum enhancement in the 1/ω dependence to the conductivity observed at 800 nm. This broad energy range, specifically 800 nm, is not resonant with any phonon modes, yet induced electronic excitations appear to be connected to superconductivity at energy scales well above the typical gap energies in the cuprates. A critical question is, What can be responsible for such an effect at 800 nm? Using time-dependent exact diagonalization, we demonstrate that the holes in the CuO 2 plane can be photoexcited into the charge reservoir layers at resonant wavelengths within a multiband Hubbard model. This orbitally selective photoinduced charge transfer effectively changes the in-plane doping level, which can lead to an enhancement of T c near the 1/8 anomaly.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Picosecond volume expansion drives a later-time insulator–metal transition in a nano-textured Mott insulator

There is significant technological interest in developing ever faster switching between different electronic and magnetic states of matter. Manipulating properties at terahertz rates requires accessing the intrinsic timescales of both electrons and associated phonons, which is possible with short-pulse photoexcitation. However, in many Mott insulators, the electronic transition is accompanied by the nucleation and growth of percolating domains of the changed lattice structure, leading to empirical timescales dominated by slowly coarsening dynamics. Here, in this study, we use time-resolved X-ray diffraction and reflectivity measurements to show that the photoinduced insulator-to-metal transition in an epitaxially strained Mott insulating thin film occurs without observable domain formation and coarsening effects, allowing the study of the intrinsic electronic and lattice dynamics. Above a fluence threshold, the initial electronic excitation drives a fast lattice rearrangement, which is followed by a slower electronic evolution into a metastable nonequilibrium state. Microscopic model calculations based on time-dependent dynamical mean-field theory and semiclassical lattice dynamics explain the threshold behaviour and elucidate the delayed onset of the electronic phase transition. This work highlights the importance of combined electronic and structural studies in unravelling the physics of dynamic transitions and the timescales of photoinduced processes. During a photoinduced phase transition, electronic rearrangements are usually faster than lattice ones. Time-resolved measurements now show that the insulator-to-metal transition in a thin-film Mott insulator is preceded by lattice reconfiguration.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Return on Investment and Sustainability of HVDC Links: Role of Diagnostics, Condition Monitoring, and Material Innovations

HVDC cable systems are becoming an upscaled technical option, compared to AC, because of various factors, including easier interconnections, lower losses, and longer transmission distances. In addition, renewables providing direct DC energy, electrified transportation, and aerospace where DC can be favored because of higher carried specific power all point in the direction of broad future usage of HV and MV DC links. However, contrary to AC, there is little return from on-field installation as regards long-term cable reliability and aging processes. This gap must be covered by intensive research, and contributing to this research is the purpose of this paper. The focus is on key points for HVDC (and MVDC) cable reliability and sustainability, from design modeling able to account for voltage transients and extrinsic aging (such as that caused by partial discharges) to the impact of aging on insulation conductivity (which rules the electric field distribution, thus aging rate). Also, recyclable and nanostructured materials, as well as health conditions, are considered. It is shown how cable design can account for accelerated aging due to voltage transients, as well as for aging-time dependence of conductivity, and how design can be free of extrinsic aging caused by PDs. Algorithms for health condition evaluations, which have additional value in a relatively new technology such as HVDC polymeric cables, are applied to insulation system aging under partial discharges, showing how they can provide an indication of insulation degradation globally or locally (weak spots) and of possible maintenance times. All of this can effectively contribute to reducing the risk of major cable breakdown and damage under operation, which would significantly affect the return on investment (ROI).

Montanari, Gian Carlo (ORCID:0000000320258693)↗

DOE OSTI2024

Pre-Transient Characterization of Historic EBR-II Pins for Transient Testing

Current interest in sodium-cooled fast reactor (SFR) designs, such as TerraPower’s Natrium Reactor, has highlighted the need for advanced reactor fuel technology development. Modern U-Zr and U- Pu-Zr pin designs are primary candidates to fuel SFRs and boast high fuel utilization capacity, increased fuel-cladding compatibility, and improved safety through inherent feedback mechanisms. Despite over 60 years of metallic fuel irradiation, uncertainties exist in the performance of the fuel system, particularly under transient overpower (TOP) and loss of flow (LOF) scenarios. Throughout historical testing within the Experimental Breeder Reactor II (EBR-II) and the Fast Flux Test Facility (FFTF), fuel behavior has demonstrated benign response to transient reactor conditions; however, accurate predictions of failure thresholds to inform operational limitations rely heavily on fuel composition, burnup, and irradiation history. In expanding TOP and LOF testing, the Transient Heat sink Overpower Response (THOR) Capsule will be used to test modern fuel technologies in a static sodium environment in the Transient Reactor Test (TREAT) Facility. The THOR capsule is highly instrumented and will provide time-dependent thermal behavior of SFR fuel pins subjected to accident conditions within TREAT. The THOR-Metallic (THOR- M) campaign aims to validate and expand historical TOP and LOF testing on high burnup U-Zr and U-Pu- Zr fuel alloys previously irradiated in EBR-II by running the rods to failure. This contribution focuses primarily on the pre-transient engineering-scale destructive and non- destructive characterization that has been conducted on both the test and sibling pins used for the TOP and LOF tests. All pins underwent visual examination, neutron radiography, element contact profilometry, and precise gamma scan. The sibling pins used for each test were further analyzed using gas assay, sampling, and recharge analysis (GASR), and optical microscopy. The results from each technique confirmed that the fuel pins were intact and devoid of any atypical developments when compared to historical data. Additionally, the analyzed measurements establish a baseline for comparison to post-transient analysis. Key fuel behaviors quanitifed include axial elongation of the fuel column, diametral strain of the pin, patterns in fluff structure geometry, changes in axial isotope distribution, evolution of constituent redistribution, porosity, and fission gas release. The pre-transient measurements and changes attributed to transient behavior from post-transient measurement will be compared to historical data to capture the behavioral dependence on composition, burnup, and irradiation history. Results from this work advance the initiatives of the THOR-M campaign, which aid in informing fuel performance models and establishing safety criteria for SFR operational limits. The novel combination of test environment, in-situ instrumentation, and comprehensive suite of characterization methods provides greater understanding of transient fuel behavior. Overall, information on the time and condition of pin failure for high burnup U-Pu-Zr will greatly expand the limited existing TOP and LOF test data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Pre-Transient Characterization of Historic EBR-II Pins for Transient Testing

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

An introduction to Spent Nuclear Fuel decay heat for Light Water Reactors: a review from the NEA WPNCS

This paper summarized the efforts performed to understand decay heat estimation from existing spent nuclear fuel (SNF), under the auspices of the Working Party on Nuclear Criticality Safety (WPNCS) of the OECD Nuclear Energy Agency. Needs for precise estimations are related to safety, cost, and optimization of SNF handling, storage, and repository. The physical origins of decay heat (a more correct denomination would be decay power) are then introduced, to identify its main contributors (fission products and actinides) and time-dependent evolution. Due to limited absolute prediction capabilities, experimental information is crucial; measurement facilities and methods are then presented, highlighting both their relevance and our need for maintaining the unique current full-scale facility and developing new ones. The third part of this report is dedicated to the computational aspect of the decay heat estimation: calculation methods, codes, and validation. Different approaches and implementations currently exist for these three aspects, directly impacting our capabilities to predict decay heat and to inform decision-makers. Finally, recommendations from the expert community are proposed, potentially guiding future experimental and computational developments. One of the most important outcomes of this work is the consensus among participants on the need to reduce biases and uncertainties for the estimated SNF decay heat. If it is agreed that uncertainties (being one standard deviation) are on average small (less than a few percent), they still substantially impact various applications when one needs to consider up to three standard deviations, thus covering more than 95% of cases. The second main finding is the need of new decay heat measurements and validation for cases corresponding to more modern fuel characteristics: higher initial enrichment, higher average burnup, as well as shorter and longer cooling time. Similar needs exist for fuel types without public experimental data, such as MOX, VVER, or CANDU fuels. A third outcome is related to SNF assemblies for which no direct validation can be performed, representing the vast majority of cases (due to the large number of SNF assemblies currently stored, or too short or too long cooling periods of interest). A few solutions are possible, depending on the application. For the final repository, systematic measurements of quantities related to decay heat can be performed, such as neutron or gamma emission. This would provide indications of the SNF decay heat at the time of encapsulation. For other applications (short- or long-term cooling), the community would benefit from applying consistent and accepted recommendations on calculation methods, for both decay heat and uncertainties. This would improve the understanding of the results and make comparisons easier.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗

DOE OSTI2020

Individual doses for super cohort members exposed to atmospheric radioiodine from the Mayak releases with an emphasis on prenatal doses

Time-dependent thyroid doses were reconstructed for 45,837 members of the Southern Urals Population Exposed to Radiation Cohort (SUPER-C) living in the region around the Mayak Production Association facilities in Russia from 131I released to the atmosphere from all relevant exposure pathways. The dose calculations are implemented in a Monte Carlo framework that produces best estimates and stochastic realizations of dose time-histories. The mean thyroid dose from 131I for SUPER-C members was 195 mGy; the median was 65 mGy. Overall, 131I-thyroid doses for about 3.6% of SUPER-C members were larger than 1 Gy. For children born in 1940-1950, the dose was about 10% higher than in previous studies because doses during the prenatal period for 9,117 individuals are included in the current work. Half of the individuals born in the region in 1950–1960 who remained in the study domain through 1972 received 9.4% or more of their total dose during the prenatal period. SUPER-C members residing in areas contaminated by discharges of liquid radioactive releases into the Techa River and the Kyshtym Accident in 1957 received 80% of their thyroid dose from airborne 131I emissions. The uncertainty in the 131I dose estimates is low enough for this approach to be used in regional epidemiological studies.

Dose assessment, 131I, airborne radioactivity, thy↗

DOE OSTI2020

A novel method to measure ion density in ICF experiments using x-ray spectroscopy of cylindrical tracers

The indirect drive approach to inertial confinement fusion has undergone important advances in the past few years. Improvements in temperature and density diagnostic methods are leading to more accurate measurements of the plasma conditions inside the Hohlraum and therefore to more efficient experimental designs. The implementation of dot spectroscopy has proven to be a versatile approach to extracting space- and time-dependent electron temperatures. In this method, a microdot of a mid-Z material is placed inside the Hohlraum and its K-shell emission spectrum is used to determine the plasma temperature. However, radiation transport of optically thick lines acting within the cylindrical dot geometry influences the outgoing spectral distribution in a manner that depends on the viewing angle. This angular dependence has recently been studied in the high energy density regime at the OMEGA laser facility, which allowed us to design and benchmark appropriate radiative transfer models that can replicate these geometric effects. By combining these models with the measurements from the dot spectroscopy experiments at the National Ignition Facility, we demonstrate here a novel technique that exploits the transport effects to obtain time-resolved measurements of the ion density of the tracer dots, without the need for additional diagnostics. Furthermore, we find excellent agreement between experiment and simulation, opening the possibility of using these geometric effects as a density diagnostic in future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Three-dimensional simulation of the University of Hawai‘i FEL oscillator with cavity desynchronization

In this paper, we present three-dimensional, time-dependent simulations of the University of Hawai‘i (UH) at Mānoa free-electron laser (FEL) oscillator. Using beam parameters from the UH facility, we study the pass-by-pass evolution of the radiation field, including its temporal, spectral, and transverse properties. At nominal bunch length, the radiation pulse develops temporal spiking near saturation, together with sideband formation and increased sensitivity to machine timing jitter. Our results show that modest cavity desynchronization can enhance the radiation energy by 63%. Large cavity desynchronization, on the other hand, can effectively suppress the spiking instabilities and improve robustness to timing fluctuations. Finally, we simulate a short-bunch operational mode with a bunch length comparable to the slippage length, which accelerates saturation and further amplifies the FEL power. Overall, these results provide a quantitative foundation for pulse control studies in the UH FEL oscillator and a critical benchmark for future experimental validation and machine optimization.

free-electron lasers↗

DOE OSTI2023

Excitonic Interactions and Mechanism for Ultrafast Interlayer Photoexcited Response in van der Waals Heterostructures

Optical dynamics in van der Waals heterobilayers is of fundamental scientific and practical interest. Based on a time-dependent adiabatic GW approach, we discover a new many-electron (excitonic) channel for converting photoexcited intralayer to interlayer excitations and the associated ultrafast optical responses in heterobilayers, which is conceptually different from the conventional single-particle picture. We find strong electron-hole interactions drive the dynamics and enhance the pump-probe optical responses by an order of magnitude with a rise time of ~3⁢0⁢0 fs in MoSe 2 /WSe 2 heterobilayers, in agreement with experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Effects of pore size on water dynamics in mesoporous silica

Water confined in mesoporous silica plays a central role in its many uses ranging from gas sorption to nanoconfined chemical reactions. Here in this paper, the influence of pore diameter (2.5–5.4 nm) on water hydrogen bond (H-bond) dynamics in MCM41 and SBA15 mesoporous silicas is investigated using femtosecond infrared vibrational spectroscopy and molecular dynamics simulations on selenocyanate (SeCN-) anions dissolved in the pores. As shown recently, SeCN- spectral diffusion is a reliable probe of surrounding water H-bond structural motions. Additionally, the long CN stretch vibrational lifetime facilitates measurement of the full range of confined dynamics, which are much slower than in bulk water. The simulations shed light on quantitative details that are inaccessible from the spatially averaged observables. The dependence of SeCN- orientational relaxation and that of spectral diffusion on the distance from the silica interface are quantitatively described with an exponential decay and a smoothed step-function, respectively. The distance-dependence of both quantities is found to be independent of the diameter of the pores, and the spatial distribution of SeCN- is markedly non-uniform, reaching a maximum between the interface and the pore center. The results indicate that the commonly invoked two-state, or core–shell, model is a more appropriate description of spectral diffusion. Using these insights, we model the full time-dependence of the measured dynamics for all pore sizes and extract the “core” and “shell” dynamical correlation functions and SeCN- spatial probability distributions. The results are critically compared to those for water confined in reverse micelles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Sensor Development for Liquid Water Detection in Dry Storage Casks: FY 2023 Status

Modeling efforts were undertaken in fiscal year (FY) 2023 to evaluate the feasibility and capability for sensing of liquid water inside canister-based dry cask storage systems (DCSSs). The focus was on the development of full-scale finite element models (FEMs) of ultrasound propagation through DCSS canister components. For these initial investigations, simulation of the baseplate component was targeted, envisioning water collection at the bottom of a vertically oriented canister. The effects of internal canister components (i.e., fuel basket, fuel assemblies) on ultrasound propagation in the baseplate component were evaluated in these efforts. The environment inside a DCSS confinement is intended to be inert and free of water to prevent potential corrosion of used fuel cladding or other internal hardware. Spent fuel assemblies are dried, after storage in water pools, to make sure water has been removed from assembly cavities. However, there is some uncertainty about the amount of residual water potentially left behind in a DCSS after drying processes, because water can become trapped in cavities or other small crevices in the surfaces formed by the fuel cladding, fuel assemblies, and other internal hardware components. Considering the complex space- and time-dependent temperature profiles in DCSSs, water may be in a liquid or gas phase depending on its location in the cask and how long the cask has been in storage. Evacuating most water and oxidizing agents contained within a canister is recommended by NUREG-1536 (NRC 2010), which covers DCSSs. As summarized by Salazar et al. (2020), existing guidance typically relies on achievable vacuum pumping pressures sustained over a hold period as a signal of dryness and water removal. However, time series data about water removal from full-scale commercial drying procedures are lacking (Hanson and Alsaed 2019). A review of drying specifications from several vendors led to the conclusion that if the specifications are followed correctly, the residual moisture in DCSSs should present an insignificant risk of cladding degradation (Knoll and Gilbert 1987). A more recent analysis concluded that much larger quantities of residual water could remain in DCSSs, but the amount would still not be expected to lead to significant corrosion of fuel cladding or other internal components (Jung et al. 2013). Industry drying procedures are mostly prescriptive in nature, and operational issues arising during the process could result in incomplete drying. A summary of operational issues and potential negative effects of residual water is provided by Salazar et al. (2019). Experimental efforts are ongoing to validate the extent of water removal in a DCSS based on drying procedures used at nuclear power plants through well-designed investigations of drying process efficacy and water retention (Durbin et al. 2021; Pulido et al. 2022a, Pulido et al. 2022b). This has been approached by simulating limited portions of a DCSS internal volume and fuel assemblies. So far, these efforts have focused on the effects of potential water trapping in the dashpot region of control rod guide tubes but are expected to be expanded to include the effects of other internal hardware features and failed fuel rod cladding. A method for detection and measurement of liquid water (Meyer et al. 2022; Meyer et al. 2021), in tandem with the drying information gained through experimental investigations, provides comprehensive bases for understanding the internal conditions of DCSSs and provides technical information to support licensing decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Multiphysics Time-Integration for Turbulent Combustion at the Exascale

Turbulent reacting flow systems are often modeled with coupled time-dependent partial differential equations (PDEs). Solving such equations can easily tax the world's largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This generic operator-splitting strategy leads to a set of ordinary differential equations (ODEs) that need to be solved as part of an "outer-loop" time-stepping approach. In many combustion applications, the ODEs to be solved can be very stiff, exhibiting timescales that span many orders of magnitude. The SUNDIALS library provides a plethora of robust time integration algorithms for solving these ODEs on exascale-capable computing hardware, yet for many complex applications (such multicomponent fuels or emissions predictions), the chemical models remain too complex to solve using reasonable resources. The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of the simulations. In this talk, I will discuss the use of the SUDIALS library of ODE solvers together with automatic code generation tools to solve complex turbulent reacting flow problems using QSSA models.

chemistry↗

DOE OSTI2021

How the Size and Density of Charge-Transfer Excitons Depend on Heterojunction’s Architecture

We have characterized the size, intensity, density, and distribution of charge-transfer (CT) excitons as a function of the acceptor–donor architecture of prototypical organic interfaces. This characterization was done by computational analysis of 17 models of varying numbers, positions, and orientations of the donor and acceptor molecules. The models’ building blocks were phenyl-C61-butyric acid methyl ester (PCBM) fullerene acceptors and dual-band donor polymers composed of thiophene, benzothiadiazole, and benzotriazole subunits. The electronic structure of the donor–acceptor complexes was computed with the time-dependent long-range-corrected density-functional tight-binding method and analyzed with the fragment-based one-electron transition density matrix. In all models, the complexes with edge-on orientation have denser spectra of low-energy CT states lying below the absorption bands compared to the complexes with face-on orientation. This CT-state distribution in edge-on complexes provides a gate to efficiently populate cold CT excitons. Moreover, the cold CT excitons have a higher degree of charge separation in the edge-on than in the face-on complexes. The CT amount and the CT exciton size generally increase with the energy of the CT states, although the electron remains localized on a single molecule in cold CT states. Delocalization over two PCBM molecules was observed for high-energy CT states. The exciton size also depends on the orientation. Larger excitons are produced by the delocalization of the electrons perpendicularly to the interface. When the delocalization is parallel, the smaller electron–hole distances yield moderately sized CT excitons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗