Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometric modelling”

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

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗
DOE OSTI2025

Toward engineering lattice structures with the material point method (MPM)

This study examines the potential of two variants of the material point method—the generalized interpolation material point (GIMP) and dual domain material point (DDMP) methods—in developing a robust computational framework for engineering lattice structures under different loading conditions. The study begins with assessing the ability of the two methods in predicting elastic buckling phenomena using column geometries with and without initial geometric imperfections. The results indicate that both methods effectively capture buckling phenomena when initial geometric imperfections are introduced. After this verification step, we create several models of tetrahedral lattice structures with varying strut diameter and orientation and subject them to quasi-static loading. We then validate the numerical results using laboratory test results. The results show that, while both methods accurately predict load–displacement curves in the pre-buckling regime, their predictive capabilities diminish in the post-buckling regime. Through visual comparison between the numerical and experimental deformed shapes, it appears that the discrepancies between model and experimental results are attributed to initial geometric imperfections in the lattices that occurred during 3D printing. We then establish a second set of lattice models where different types of initial geometric imperfections are considered. The results from these models show that imperfections have a negligible influence in the pre-buckling regime but affect the behavior considerably in the post-buckling regime. As a final step in this work, we subject the lattice models to impact loading and employ hypothetical soft and stiff materials. These results show that the lattice stiffness, which depends on material stiffness, strut diameter, and orientation, significantly influences the ability of a lattice structure to resist impact. In particular, we find that a stiffer lattice (i.e., one made with a stiff material and thicker struts) is capable of absorbing more energy than a softer one during impact. Although material nonlinearities, inelasticity, and detailed contact formulations are not considered in this study, the findings obtained herein lay the groundwork for engineering lattice structures under extreme loading conditions through a simulation-driven framework based on particle-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Ducted Assembly Steady-State Heat Transfer Software (DASSH): Theory Manual

The Ducted Assembly Steady-State Heat transfer software (DASSH) is being developed at Argonne National Laboratory to obtain a steady-state thermal fluids result for the coolant flow and temperature distribution in a hexagonally gridded reactor core with ducted assemblies. In DASSH, each assembly is assumed to contain a hexagonal bundle of wire-wrapped pins. DASSH divides the coolant in the pin bundle into subchannels on which it applies energy balances to determine coolant and duct temperatures. DASSH offers multiple models for inter assembly heat transfer and features the ability to calculate pin temperatures. This document describes the models and equations implemented in DASSH. It includes a description of the geometric system to which the models are applied and a derivation of the subchannel energy balance formulation is presented. It also reviews the available correlations and provides an overview of the process by which power distributions are generated. Finally, it covers the code execution procedure, discussing the main code objects and their interactions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗
DOE OSTI2022

Accurate tight-binding model for twisted bilayer graphene describes topological flat bands without geometric relaxation

A major hurdle in understanding the phase diagram of twisted bilayer graphene is the roles of lattice relaxation and electronic structure on isolated band flattening near magic twist angles. Here in this work, the authors develop an accurate local environment tight-binding model fit to tight-binding parameters computed from ab initio density-functional theory calculations across many atomic configurations. With the accurate parametrization, it is found that the magic angle shifts to slightly lower angles than often quoted, from around 1.05° to around 0.99°, and that isolated flat bands appear for rigidly rotated graphene layers, with enhancement of the flat bands when the layers are allowed to distort. Study of the orbital localization supports the emergence of fragile topology in the isolated flat bands without the need for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
DOE OSTI2026

Novel application of the modal strain energy technique for state-of-the-art damping predictions

Structural damping, which measures the energy dissipation of a vibrating structure, is a key modeling input for lightweight structures but is notoriously hard to predict. This work utilizes vibration-based measurements of centimeter-scale coupons and the modal strain energy approach to predict structural damping of a lightly damped structure. The approach was originally validated with panels shorter than a meter in length. This work extends the validation to a 2.75 m beam made of unidirectional and biaxial glass fiber laminates bonded by adhesive. The comparison between three-dimensional finite element model predictions and full-scale experimental measurements of damping show an average error of 5.2% for the first five modes. Additionally, the modal strain energy approach is newly applied with a one-dimensional geometrically exact beam theory model and a two-dimensional sectional analysis solver. This beam approach accurately predicts the damping behavior of the first bending modes but loses accuracy for higher order modes that are dominated by three-dimensional effects. This novel approach provides faster simulations while allowing arbitrary beam cross sections. The paper also investigates traditional and high-force dynamic mechanical analysis to measure structural damping of coupons. Both alternatives show significant errors in attempted validation against the theoretical thermoelastic damping of aluminum coupons.

17 WIND ENERGY↗
DOE OSTI2024

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2021

Theoretical Analysis of Experimental Data of Sodium Diffusion in Oxidized Molybdenum Thin Films

In this work, the diffusion process of sodium (Na) in molybdenum (Mo) thin films while it was deposited on soda lime glass (SLG) was studied. A small amount of oxygen was present in the chamber while the direct-current (DC) magnetron sputtering was used for the deposition. The substrate temperatures were varied to observe its effect. Such molybdenum films, with or without oxidations, are often used in thin film solar cells, either as back contact or as hole transport layers. Secondary ion mass spectrometry (SIMS) was used to quantify the concentration of the species. A grain diffusion mechanistic model incorporating the effect of grain and grain boundary geometrical shape and size was developed. The model was used to provide an in-depth theoretical analysis of the sodium diffusion in molybdenum thin films that lead to the measured SIMS data. It was observed that not only diffusion coefficients should be considered when analyzing diffusion processes in thin films but also the ratio of grain boundary size to grain size. Both depend on substrate temperature and directly affect the amount of diffused species in the film. The data were analyzed under the light of the film growth speed versus diffusion front speed, the effect of oxygen content, and the effect of substrate temperature on the overall diffusion process. The temperature inversely affects the ratio of grain boundary size and grain size and directly affects the diffusion coefficient, which leads to a preferable temperature at which the highest amount of alkali can be found in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Orders of vanishing and U(1) charges in F-theory

Many interesting questions about F-theory models, including several concerning the F-theory swampland, involve massless matter charged under U(1) gauge symmetries. It is therefore important to better understand the geometric properties of F-theory models realizing various U(1) charges. We propose that, for F-theory models described by elliptic fibrations in Weierstrass form, the U(1) charge of light matter is encoded in the orders of vanishing of the section components corresponding to the U(1) gauge symmetry. We give specific equations relating the U(1) charges to the orders of vanishing that seem to hold for both U(1)-charged singlets and for matter additionally charged under a simply-laced nonabelian gauge algebra. Our formulas correctly describe properties of F-theory models in the prior literature, and we give an argument that they should describe the orders of vanishing for arbitrarily high U(1) charges. They also resemble formulas for the p-adic valuations of elliptic divisibility sequences developed by Stange. These proposals could serve as a U(1) analogue of the Katz-Vafa method, allowing one to determine U(1) charges without resolution. Additionally, they predict geometric information about F-theory models with general U(1) charges, which may be useful for exploring the F-theory landscape and swampland.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Gravo-aeroelastic scaling of a 13-MW downwind rotor for 20% scale blades

A 105-m, 13-MW two-bladed downwind Segmented Ultralight Morphing Rotor (SUMR-13) blade was gravo-aeroelastically scaled by 20% to a 20.87-m-long demonstrator blade and confirmed through structural ground testing. The subscale model was achieved through geometric scaling and by aeroelastic scaling principles based on operational flapwise deflections combined with rotational and structural frequencies while retaining the turbine tip-speed ratio. In particular, the subscale demonstrator was designed to replicate, as closely as possible, the nondimensional geometry, the ratio of centrifugal to gravitational moments, the tip-speed ratio, and the nondimensional rotation rate. The intent for this demonstrator was to achieve the same nondimensional flapwise blade deflections and dynamics of the full-scale 13-MW rotor. The manufactured SUMR-D blade resulted in less than half of the mass of the conventional two-bladed Controls Advanced Research Turbine (CART2) rotor blade based on scaling and a lower power rating, though with some differences in mass and stiffness from the ideal scaled-down design to meet safety requirements at the test site. To achieve proper scaling, operational pitch control set points were altered to account for the differences by evaluating simulated operation of both the SUMR-13 and SUMR-D rotors. Structural testing of the SUMR-D blade investigated the response to well-defined flapwise loads and indicated that the subscale blade had the appropriate elastic properties needed for both scaling and for safe operational field testing.

17 WIND ENERGY↗
DOE OSTI2025

cymyc: $\underline{C}$alabi-$\underline{Y}$au $\underline{M}$etrics, $\underline{Y}$ukawas, and $\underline{C}$urvature

We introduce cymyc, a high-performance Python library for numerical investigation of the geometry of a large class of string compactification manifolds and their associated moduli spaces. We develop a well-defined geometric ansatz to numerically model tensor fields of arbitrary degree on a large class of Calabi-Yau manifolds. cymyc includes a machine learning component which incorporates this ansatz to model tensor fields of interest on these spaces by finding an approximate solution to the system of partial differential equations they should satisfy.

differential and algebraic geometry↗
DOE OSTI2021

Role of a 3D Structure in Energy Transfer in Mixed-Ligand Metal–Organic Frameworks

We present a detailed investigation of the photophysical properties of mixed-ligand pyrene- and porphyrin-based metal–organic frameworks (MOFs) as a function of their 3D structure. Solvothermal reactions between metal salts (InCl 3 , Zr(acac) 4 , and ZrCl 4 ) and suitable ratios of 1,3,6,8-tetrakis(p-benzoic acid)pyrene (TBAPy) and meso-tetrakis- (4-carboxyphenyl)porphyrin (TCPP) were performed to prepare a series of mixed-ligand ROD-7, NU-901, and NU-1000 MOFs. Time-resolved and steady-state fluorescence measurements were conducted on the mixed-ligand MOFs to study their photophysics. Based on the results, we concluded that upon excitation of TBAPy linkers in the MOFs, singlet excitation energy migrates across TBAPy linkers until it finds a TCPP unit. TCPP acts as an energy trap and quenches the excitation. The efficiency of TBAPy-to-TBAPy energy transfer was found to be sensitive to the structural parameters of MOFs. Analysis of steady-state and time-resolved fluorescence data revealed that excitation energy transfer (EET) is most efficient in ROD-7, followed by NU-901 and NU-1000. We propose that topology that invokes the shorter interchromophoric distances between TBAPy linkers in ROD-7 is responsible for its higher EET efficiency. The distance dependence of the EET rate constant (kEET) was investigated to gain insight into the mechanistic aspects of energy transfer in MOFs. This study revealed that (a) energy transfer in MOFs deviates from the classical Förster model and (b) the geometrical arrangement of linkers influences the mechanism of EET in MOFs. A theoretical investigation was also performed to determine energy-transfer rate constants along different directions and assess the directionality of energy transfer in these MOFs. The magnitude of rate constants indicated that energy transfer in ROD-7 should be highly anisotropic along the stacking direction. Furthermore, these findings suggest that ROD-7 is a promising candidate to play the role of the light-harvesting and energy-transfer component in solar energy conversion devices, where directional energy transfer is required.

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

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗