Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static analysis”

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 487 records · Page 27

Characterizing the Effect of Capillary Heterogeneity on Multiphase Flow Pulsation in an Intermediate–Scale Beadpack Experiment using Time Series Clustering and Frequency Analysis

An intermediate-scale beadpack drainage experiment was conducted to investigate how simple layered lamination heterogeneity affects CO 2 flow. Two simple layers of capillary barriers are manually packed in the tank and slow drainage was carried out using analog fluids to mimic the capillary- and gravity-dominated CO 2 upward migration process in deep saline aquifers. Nonwetting phase saturation time series clustering analysis and frequency analysis have been conducted on the experimental data. Additionally, modified invasion percolation numerical simulations were done on a digital model of the beadpack to compare to experimental results. Results show that capillary barriers can lead to strong pulsation behavior, which in turn can cause unexpected early breaching through other barriers. The inlet pressure is found to be able to respond to saturation changes in far regions of the domain, indicating that the wetting phase can transmit pressure changes from the other phase. Although static simulations were not able to capture all the dynamic behavior observed in the experiment, Monte Carlo composite simulation results combining many different realizations can better illustrate how the nonwetting phase will behave in the heterogeneous domain. Furthermore, our results suggest the need for CO 2 storage site selection with preference given to aquifers with more capillary barriers with finer grain sizes to avoid flow pulsation and to retard plume upward migration.

58 GEOSCIENCES↗

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗

Convectively coupled Rossby–Gravity waves in a field campaign: How they are captured in reanalysis products

Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.

Wang, Xiaocong↗

Damage Tolerance Assessment of Laser Clad Repairs of Coarse Grain Ti-6Al-4V

Background: Laser cladding offers a promising approach to repairing aviation components, with the ability to restore the geometric shape and the static and fatigue strength. Objective: The damage tolerance of a coarse grain Ti-6Al-4 V material is assessed via three sets of samples: (i) baseline, (ii) laser clad repaired, and (iii) laser cladding followed by heat treatment. Method: The specimens were subjected to fatigue loading using a marker band technique, in which cracks initiated from pre-defined laser slots and the small crack growth rates were measured by quantitative fractography. Additional characterization was performed on these samples, including hardness measurements, microstructure analysis, and multimodal residual stress analyses. Results: The region of the specimen undergoing laser cladding exhibited higher hardness values and a compressive residual stress, based on the slitting method. The samples subjected to heat treatment after cladding displayed negligible differences in grain size and minimal residual stress relaxation. Conclusions: The results demonstrated significant variability, with data demonstrating that cracks in the laser clad region (exhibiting the largest compressive residual stress) grew faster than cracks in the substrate region of the sample or in the baseline cases. Finally, these counterintuitive results are rationalized with additional analyses from synchrotron-based energy dispersive X-ray diffraction and high energy X-ray diffraction, demonstrating significant spatial variability in the residual elastic strains in the millimeter-sized prior β grains within the material.

36 MATERIALS SCIENCE↗

Deformation Mechanism Transition in Additively Manufactured Compositionally Graded Fe-Base Alloys

Microstructure-dependent deformation and fracture behavior was investigated for an additively manufactured compositionally graded alloy (CGA) printed using the laser-directed energy deposition (L-DED) method to explore an alternative approach for dissimilar metal joints in nuclear energy systems. The electron backscatter diffraction (EBSD) maps from scanning electron microscopy (SEM) display a clear microstructural transition with decreasing austenite-forming elements (Ni and Mn), from an austenite (γ) dominant structure, to a complex composite structure containing ferrite (α), martensite (α') and retained austenite, and then to a fully ferritic structure. EBSD data were recorded in situ during tensile testing in SEM, and the evolution of the deformation mechanism and microstructure was characterized using Kikuchi diffraction pattern analysis. Complementary analysis for high-resolution features was also performed using scanning transmission electron microscopy (STEM). The Ni/Mn-rich austenite-dominant microstructures showed a complex deformation mechanism of two-step martensitic transformation (γ→ε→α'), whereas the minor austenite phase retained in the ferrite and/or martensite matrix showed a single transformation route (γ→α'). Ordinary dislocation glide and twinning via partial dislocation glide were observed in the austenite deformation. Meanwhile, the ferrite and martensite grains deformed mainly by ordinary dislocation slips and grain rotation. Furthermore, static tensile fracture was also highly dependent on local composition and phase constituents.

36 MATERIALS SCIENCE↗

CAMFeND: Credibility-Aware Multimodal Fake News Detection with Rotational Attention

In the evolving digital landscape, fake news is a significant challenge, influencing public perception and decision-making. Traditional detection approaches focus on single-modal data or simple multimodal fusion, often overlooking deeper interactions and news credibility. We propose a novel model addressing these limitations by introducing rotational attention and news domain information as a feature. Unlike static attention mechanisms, our rotational attention dynamically shifts query, key, and value roles across text and image inputs, enabling richer cross-modal interaction. Incorporating news domain information further enhances the model’s reliability by associating news posts with top domains extracted from Google search results, reducing false detections. This approach assesses both the content and the broader web context in which the news is discussed. Our model outperforms existing state-of-the-art methods by providing deeper, layered multimodal integration and domain information analysis, resulting in a more robust and adaptive fake news detection system.

Gupta, Nidhi↗

Observation of high-pressure polymorphs in bulk silicon formed at relativistic laser intensities

Silicon polymorphs with exotic electronic and optical properties have recently attracted significant attention due to their wide range of useful band gap characteristics. They are typically formed by static high-pressure techniques, which limits the crystal structures that can be made. This constitutes a major obstacle to study these polymorphs and their incorporation into existing technology. Approaches have attempted to address this shortcoming through using dynamic conditions and chemical precursor materials. Here, we report on an approach to create unusual crystal structures deep in the bulk of a silicon crystal by irradiating it with a laser pulse at ultrarelativistic intensity of up to 7.5 × 10 19 W/ c m 2 . Laser-generated electrons with MeV energy swiftly penetrate the target with speed close to the speed of light and deposit their energy into a large volume across the whole thickness of the sample. The relativistic electron current creates, via branching propagation and ionization, high-energy-density conditions for thermodynamically nonequilibrium phase transformation paths into new crystal polymorphs. X-ray microdiffraction and synchrotron x-ray diffraction analyses indicate, along with conventional dc-Si, the presence of exotic silicon structures in the bulk of the laser intact target volume. These structures are identified as body-centered bc8-Si, rhombohedral r8-Si, hexagonal-diamond hd-Si, and the tetragonal Si-VIII, all phases of Si that have previously been made through static techniques. Additionally, simple-tetragonal st12-Si and body-centered tetragonal bt8-Si were observed along with signatures of not yet identified diffraction spots. Both st12-Si and bt8-Si have only been observed in ultrafast laser microexplosion conditions at much lower laser intensity ∼ 10 14 W/ c m 2 and within a micron-thin surface layer. The findings here are supported by direct observation of nanoparticles with high-resolution transmission electron microscopy and corresponding fast Fourier transform analysis of their interatomic distances. The presented analyses of absorbed laser energy, generation of the MeV electron current, and deposition of energy across the whole target thickness provide a solid basis for drawing the conclusion that the observed silicon polymorphs were produced because of laser-generated high-energy electrons fast-penetrating deeply into the bulk of silicon. In contrast to solid-solid transformations, the plasma-solid transitions offer a paradigm for the creation of exotic, high-energy density materials inside the bulk of the sample by using laser pulses at relativistic intensities. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic infrared gas analysis from longleaf pine fuel beds burned in a wind tunnel: observation of phenol in pyrolysis and combustion phases

Pyrolysis is the first step in a series of chemical and physical processes that produce flammable organic gases from wildland fuels that can result in a wildland fire. We report results using a new time-resolved Fourier transform infrared (FTIR) method that correlates the measured FTIR spectrum with an infrared thermal image sequence, enabling the identification and quantification of gases within different phases of the fire process. The flame from burning fuel beds composed of pine needles (Pinus palustris) and mixtures of sparkleberry, fetterbush, and inkberry plants was the natural heat source for pyrolysis. Extractive gas samples were analyzed and identified in both static and dynamic modes synchronized to thermal infrared imaging: a total of 29 gases were identified including small alkanes, alkenes, aldehydes, nitrogen compounds, and aromatics, most previously measured by FTIR in wildland fires. This study presents one of the first identifications of phenol associated with both pre-combustion and combustion phases using ca. 1 Hz temporal resolution. Preliminary results indicate ~2.5× greater phenol emissions from sparkleberry and inkberry compared to fetterbush, with differing temporal profiles.

09 BIOMASS FUELS↗

Modeling pre-Exascale AMR Parallel I/O Workloads via Proxy Applications

The present work investigates the modeling of preexascale input/output (I/O) workloads of Adaptive Mesh Refinement (AMR) simulations through a simple proxy application. We collect data from the AMReX Castro framework running on the Summit supercomputer for a wide range of scales and mesh partitions for the hydrodynamic Sedov case as a baseline to provide sufficient coverage to the formulated proxy model. The non-linear analysis data production rates are quantified as a function of a set of input parameters such as output frequency, grid size, number of levels, and the Courant-Friedrichs-Lewy (CFL) condition number for each rank, mesh level and simulation time step. Linear regression is then applied to formulate a simple analytical model which allows to translate AMReX inputs into MACSio proxy I/O application parameters, resulting in a simple “kernel” approximation for data production at each time step. Results show that MACSio can simulate actual AMReX nonlinear “static” I/O workloads to a certain degree of confidence on the Summit supercomputer using the present methodology. The goal is to provide an initial level of understanding of AMR I/O workloads via lightweight proxy applications models to facilitate autotune data management strategies in anticipation of exascale systems.

Godoy, William↗

Incorporating Valve Switching Losses Into a Static Optimal Control Algorithm for the Hybrid Hydraulic-Electric Architecture (HHEA)

The Hybrid Hydraulic-Electric Architecture (HHEA) has, in recent years, been proposed as an energy efficient alternative to conventional load-sensing architectures in mobile machines such as excavators and wheel-loaders. HHEA leverages the advantages of hydraulic power and electric power to eliminate throttling valves while also improving the energy and control performance of the system. The architecture utilizes a set of common pressure rails to provide a majority of power and and a small electric motor driven pump to modulate this power to meet the exact demand. Previous work has developed a computationally efficient Lagrange Multiplier approach for determining the optimal pressure rail selections that minimizes the energy losses in the system. The static model used considers only the energy use for each pressure rail selection but not the losses associated with the valves during the transition. This paper presents an approach to include the switching losses in the model and in the optimization procedure. To capture the switching losses, switching events between different rails and at various input and output flow rates were simulated with consideration of valve spool dynamics. A parameterized model that summarizes the losses is then obtained, allowing switching losses to be added to the previous energy analysis. The performance of the switching loss model was compared with reference data obtained from a high-fidelity simulation model. To incorporate the switching losses into optimal control algorithm, an efficient dynamic programming approach that prevents frequent switching is adopted in place of the Lagrange multiplier approach. The overall effect of switching losses on the energy consumption and optimal control decisions is presented. In general, switching losses contribute to about 9–10% of input energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Sweep-tracing algorithm: in silico slip crystallography and tension-compression asymmetry in BCC metals

Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.

36 MATERIALS SCIENCE↗

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

High-pressure phase transitions in a laser directed energy deposited Fe-33Cu Alloy

Here, additively manufactured Fe-Cu alloys contain both equilibrium face-centered cubic (FCC) and metastable body-centered cubic (BCC) crystal structure Cu precipitates depending on their size. However, the stability of these nanoscale precipitates under extreme conditions such as high pressures has not been reported. This study investigates the phase transformations and microstructural stability of laser directed energy deposition (DED-LB) made Fe67Cu33 alloy (nominal composition in at.%) under high static pressure deformation using in-situ synchrotron X-ray diffraction under pressure, postmortem high-resolution scanning transmission electron microscopy (HR-STEM), and molecular dynamics (MD) simulations. In-situ XRD results reveal a reversible phase transformation from BCC to the hexagonal close-packed (HCP) structure in the Fe grains at an onset pressure of 16.4 GPa, significantly higher than reported in the literature for pure Fe. Although no high-pressure phase transition was observed in the FCC Cu grains through XRD, HR-STEM analysis uncovers a phase transition to the HCP structure in nanoscale metastable BCC Cu precipitates within the BCC Fe matrix. After decompression, the Fe matrix reverted back to the BCC structure with periodic lath martensite, while regions of the nanoscale BCC Cu precipitates retained the metastable HCP structure. MD simulations support the BCC → HCP transition in the embedded nanoscale coherent Cu precipitates, consistent with the classical Burgers mechanism. Thus, by leveraging in-situ XRD observation, postmortem high-resolution S/TEM microscopy, and MD simulations, this study offers profound insights into the distinctive phase transformations induced by high pressure in DED-LB Fe67Cu33 alloy, which is distinguished by its hierarchical microstructure.

36 MATERIALS SCIENCE↗

Inelastic relaxation processes in amorphous sodium silicates

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na 2 O-SiO 2 (NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na 2 O-SiO 2 (NS15) composition, the 25% Na 2 O-SiO 2 (NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na 2 O-SiO 2 compositions. Coordination analysis of the Na + ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na + ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

36 MATERIALS SCIENCE↗

Experimental Analysis of Advanced Turbine for Supercritical CO2 Power Cycle

An advanced first stage high pressure turbine blade for a supercritical CO2 power cycle is tested in the Big Rig for Aerothermal Stationary Turbine Analysis (BRASTA) annular cascade at the Purdue Experimental Turbine Aerothermal Laboratory (PETAL) high pressure blow down facility alongside a baseline blade for comparison of aerodynamic performance. Both geometries are tested simultaneously using a novel off-axis rig design to allow for the blade geometries to be scaled up to achieve higher Reynolds numbers. The off-axis design necessitates the use of discrete sectors of airfoils, which are designed and additively manufactured in house using an mSLA printer. Printing the blades allows for unique routing of passage for blade surface static pressure taps that contour to the blade shape instead of requiring a straight line view from surface tap to egress, allowing all instrumentation for 15%, 50%, and 85% span to exist in the same passage. A flow conditioning gauze [1] is placed upstream of the blade passages to impart pressure, Mach, and swirl profiles to mimic the rotor relative frame inlet conditions to the blade. Downstream, a sonic valve is used to alter the backpressure to achieve different pressure ratios for different blowdown setpoints. The design of the rig, sonic wheel, and instrumentation are discussed along with the manufacturing and GD&T of the blades with respect to the build-up of this novel experimental apparatus. Initial commissioning data is obtained from inlet total pressure and temperature rakes, blade static pressure taps, and exit total pressure rakes. Oil visualization is performed using a silicone oil mixture containing titanium dioxide and pigmented for contrast of the blade suction side surface with the hub and shroud endwalls. The viscosity of the oil is tuned so that the oil does not fully thin and blow away during the time of the test, and endoscopic cameras are placed in the rig for live monitoring and recording. The preliminary oil visualization and rake profiles are presented, showcasing the operability of this additively manufactured off-axis flow path design.

Tuite, Logan [Purdue University]↗

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗

Fundamental Studies of the Vibrational, Electronic, and Photophysical Properties of Tetrapyrrolic Architectures

The ability to capture and utilize light in the near-ultraviolet (NUV), visible and near-infrared (NIR-I and NIR-II) spectral regions (i.e., 320–400, 400–700, 700–1000, 1000–1700 nm) is essential for any solar-energy conversion scheme. Nature employs chlorophylls and bacteriochlorophylls in light-harvesting architectures to absorb light in the blue and red/NIR regions. Accessory pigments (carotenoids, bilins) augment absorption of the (bacterio)chlorophylls in the green region. The harvested energy is funneled to a reaction center protein, where charge separation occurs. Subsequent migration of the electron and the hole stabilizes and stores the energy from light via redox chemistry. The long-term objective of the Bocian/Holten&Kirmaier/Lindsey research program under this DOE grant has been to develop tetrapyrrole-based molecular architectures that absorb sunlight, funnel energy and separate charge with high efficiency. Integral to the program has been iterative cycles of design, synthesis and characterization that provided deep insights into the relationships between chemical composition, electronic structure, and key static and dynamic properties (vibrational, redox, photophysical, energy/charge transfer) of tetrapyrrolic systems. Such architectures included monomers, dyads, larger arrays, and complexes with accessory components. The objective was to develop molecular designs and guiding principles to enhance current and future energy-conversion schemes. Molecular arrays targeted to address one or more fundamental questions concerning light harvesting and energy/charge transfer were constructed from analogues of the naturally occurring hemes, chlorophylls and bacteriochlorophylls. Diverse, tunable synthetic building blocks were prepared that spanned the three respective tetrapyrrole families, which are the porphyrins, chlorins and bacteriochlorins. Thus, the research focused on porphyrins as well as synthetic surrogates for chlorophylls (chlorins, 13 1 -oxophorbines and chlorin-imides) and bacteriochlorophylls (bacteriochlorins, bacterio-13 1 -oxophorbines and bacteriochlorin-imides), generically termed hydroporphyrins. Although the three tetrapyrrole classes (porphyrins, chlorins and bacteriochlorins) absorb light strongly in the violet-blue spectral region, the long-wavelength absorption band typically lies in the green-orange, red, and NIR regions, respectively, with increasing intensity. Understanding the spectra, electronic structure, and energy/charge-transfer properties of such tetrapyrrolic macrocycles is of central importance for the rational design of molecular architectures for solar-energy conversion. Our integrated program of molecular design and synthesis coupled with a variety of spectroscopic, electrochemical, and computational studies have probed from first principles how structural and electronic properties of tetrapyrrolic macrocycles dictate spectral properties as well as the rates of ground-state hole/electron transfer and excited-state energy flow in multicomponent architectures. Individual molecules and multicomponent architectures were designed to test ideas of fundamental importance, often requiring the development of new synthetic methodology. The members of the collaborative team had almost daily discussions by phone and/or e-mail concerning design of molecules, flow of compounds between the labs, planning of physical characterization studies, discussing results and analysis and integrating into design of next generation architectures, and the preparation of manuscripts. Furthermore, students and postdocs in the different labs routinely communicated with one another to facilitate the advancement of the research activities. In short, a highly integrated and collaborative research program was well established among the groups. The research effort involved molecular design and synthesis of synthetic molecular architectures by the Lindsey group integrated with physicochemical and photophysical characterization by the Bocian group and the Holten&Kirmaier group (Figure 2). The Bocian group carried out electrochemical, electron paramagnetic resonance (EPR), resonance Raman (RR), and Fourier-transform infrared (FT-IR) studies, as well as density functional theory (DFT) calculations and the time-dependent extension (TDDFT) to gain insight into excited-state properties. The Holten&Kirmaier group carried out static and time-resolved absorption and fluorescence spectroscopy studies and simulated absorption spectra using molecular orbital (MO) energies from DFT as input to the four-orbital model to complement the TDDFT calculations. The combined measurements provided understanding of the vibrational/electronic properties of the individual molecules and the changes that occur upon incorporation into multicomponent architectures. This information underpinned elucidating the mechanisms and timescales of ground-state hole/electron transfer and excited-state energy and charge transfer.

14 SOLAR ENERGY↗