Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spectral simulations”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Spectral–Dynamic Model of the Hot Plasma Layer Expansion

We propose a theoretical model describing the expansion of a plasma layer into vacuum for an arbitrary electron plasma component temperature. Comparison with the known limiting cases of the quasi-neutral outflow and the Coulomb explosion, as well as with the results of the 1D electrostatic simulation, has revealed a high accuracy of the proposed model in the description of the spectral–energy and spatial characteristics of ions accelerated during the plasma expansion. The procedure for obtaining the relations between the characteristics of accelerated ions and the laser pulse and target parameters is described with regard to qualitative predictions and the description of the results of numerical kinetic simulation and experiments on the laser and plasma acceleration of ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Contribution to Solar Brightness of Small-Size Magnetic Elements

Irradiance variability is mostly driven by surface magnetism, each magnetic feature contributing in a different manner. The contribution of small-size magnetic elements observed ubiquitously at high-spatial resolution on the solar photosphere, is still debated, as such features are mostly unresolved on full-disk images employed to model irradiance variability. Understanding the contribution to solar brightness of small-size magnetic elements, especially in quiet regions, will help understanding irradiance variability on the decadal and longer temporal scales, which, in turn, are fundamental to understand the role of the Sun on the Earth’s climate. We present a preliminary study of the brightness of small-size magnetic elements using high-spatial resolution spectro-polarimetric observations acquired with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) during Commissioning Phase 1. We focus on the Intensity contrast vs magnetic field relation, which is a fundamental observable in irradiance studies. Previous studies mostly focused on the FeI 630.1/630.2 nm spectral range and so to contextualize our results we focus here on this spectral region as well. By comparing DKIST results with MURaM magneto-hydrodynamic simulations, and with previous results published in the literature, we conclude that our results are in agreement with studies conducted on HINODE observations. We plan to extend the analysis to UV and IR ViSP observations and VBI blue continuum in order to better constrain the quiet sun contribution to spectral irradiance variations.

SMD↗
NASA NTRS2018

A Geostationary Instrument Simulator for Aerosol Observing System Simulation Experiments

In the near future, there will be several new instruments measuring atmospheric composition from geostationary orbit over North America, East Asia, and Europe. This constellation of satellites will provide high resolution, time resolved measurements of trace gases and aerosols for monitoring air quality and tracking pollution sources. This paper describes a detailed, fast, and accurate (less than 1.0% uncertainty) method for calculating synthetic top of the atmosphere (TOA) radiances from a global simulation with a mesoscale free running model, the GEOS-5 Nature Run, for remote sensing instruments in geostationary orbit that measure in the ultraviolet-visible spectral range (UV-Vis). Generating these synthetic observations is the first step of an Observing System Simulation Experiment (OSSE), a framework for evaluating the impact of a new observation or algorithm. This paper provides details of the model sampling, aerosol and cloud optical properties, surface reflectance modeling, Rayleigh scattering calculations, and a discussion of the uncertainties of the simulated TOA radiance. An application for the simulated TOA radiance observations is demonstrated in the manuscript. Simulated TEMPO (Tropospheric Emissions: Monitoring of Pollution) and GOES-R (Geostationary Operational Environmental Satellites) observations were used to show how observations from the two instruments could be combined to facilitate aerosol type discrimination. The results demonstrate the viability of a detailed instrument simulator for radiance measurements in the UV-Vis that is capable of accurately simulating high resolution, time-resolved measurements with reasonable computational efficiency.

Castellanos, Patricia↗
NASA NTRS1996

The incremental unknowns-a multilevel scheme for the simulation of turbulent channel flows

In numerical simulation of complex flows, it is important to identify different length scales of the flow and treat them differently. In this report, we introduce a new multilevel scheme for simulating turbulent channel flows. Two different versions of the scheme, namely the spectral and finite difference versions, are presented. The spectral version of the scheme is based on a spectral-Galerkin formulation which provides a natural decomposition of the flow into small and large wavelength parts, and which leads to linear systems that can be solved with quasi-optimal computational complexity. In the finite difference version, the Incremental Unknown (IU) is used to separate the length scales. Preliminary numerical results indicate that the scheme is well suited for turbulence computations and provides results which are comparable to that by Direct Numerical Simulation (DNS) but with significantly less CPU time.

Chen, M.↗
DOE OSTI2020

Consideration of temperature-dependent emissivity of selective emitters in thermophotovoltaic systems

Spectral emissivity control is paramount for designing a high-efficiency selective emitter surface required for thermophotovoltaic (TPV) applications. Owing to the temperature dependency of materials optical constants, the spectral properties of a selective emitter surface changes with the emitter temperature. This paper presents the fabrication of a multilayer metal-dielectric ( ${\rm Si}_{3}{\rm N}_{4}/\rm W/Si_{3}{\rm N_{4}}$ ) coated tungsten selective emitter aimed for GaSb-based TPV systems and studies the dependence of its surface spectral emissivity, $\varepsilon (\lambda)$ , upon a temperature ranging from 300 K to 1500 K. Both the simulation and experimental methods were used to characterize $\varepsilon (\lambda)$ as a function of temperature. For wavelengths less than 1.4 µm, $\varepsilon (\lambda)$ was found to have a minimal dependence on temperature. Beyond 1.4 µm, $\varepsilon (\lambda)$ increases with the temperature. At 1.55 µm, the simulation and experimental data estimated a ${\sim}{{4}}\%$ greater emissivity at 1500 K than at room temperature. At 1500 K, the increased $\varepsilon (\lambda)$ at longer wavelengths lowered the spectral conversion efficiency of the selective emitter from 58% to 47%. The output power density, sub-bandgap loss, and TPV conversion efficiency ( ${\eta _{\rm TPV}}$ ) for a GaSb cell illuminated by the selective thermal emitter at 1500 K were estimated. ${\eta _{\rm TPV}}$ drops from 13.7% to 11% due to the increased sub-bandgap emission at 1500 K. Essential approaches for mitigating the sub-bandgap losses to further improve ${\eta _{\rm TPV}}$ are also discussed.

42 ENGINEERING↗
NASA NTRS2020

Optimized Profile Retrievals of Aerosol Microphysical Properties from Simulated Spaceborne Multiwavelength Lidar

This work is an expanded study of one previously published onretrievals of aerosol microphysical properties from space-borne multiwavelengthlidarmeasurements. The earlier studiesand this one weredone in the framework of the NASA Aerosol-Clouds-Ecosystems (now the Aerosol Clouds Convection and Precipitation) NASA mission. The focus here is on the capabilities of a simulated spacebornemultiwavelengthlidar system for retrieving aerosol complex refractive index (m = mr+ imi) and spectral single scattering albedo (SSA(λ)), although other bulk parameters such as effective (reff) radius and particle volume (V) and surface (S) concentrations are also studied. The novelty presented here is the use of recently published, case dependent optimized-constraints on the microphysical retrievals using three backscattering coefficients (β) at 355, 532 and 1064 nm and two extinction coefficients (α) at 355 and 532 nm, typically known as the stand-alone 3β+2α lidar inversion. Case-dependent optimized-constraints (CDOC) limit the ranges of refractive index, both real (mr) and imaginary (mi) parts, and of radii that are permitted in the retrievals. Such constraints are selected directly from the 3β+2α41measurements through an analysis of the relationship between spectral dependence of aerosol extinction-to-backscatter ratios (LR) and the Ångström exponent of extinction. The analyses presented here for different sets ofsize distributions and refractive indices reveal that the direct determination of CDOCareonly feasible for cases where the uncertaintiesin the input optical data areless than 15 %.Forthe same simulated spacebornesystem and yield than in Whiteman et al., (2018), we demonstrated that the use of CDOC as essential for the retrievals of refractive index and also largely improved retrieval of bulk parameters. A discussion of the global representativeness of CDOC is presented using simulated lidar data from a 24-hour satellite track using GEOS model output to initialize the lidar simulator.We found that CDOCare representative of many aerosol mixtures in spite of some outliers (e.g. highly hydrated particles) associatedwith the assumptions of bimodal size distributions and of the same refractive index for fine and coarse modes. Moreover, sensitivity tests performed using synthetic data reveal that retrievals of imaginary refractive index (mi) and SSA are extremely sensitive to β(355).

NASA Aerosol-Clouds-Ecosystems↗
DOE OSTI2021

Inhomogeneous mean-field approach to collective excitations near the superfluid–Mott glass transition

Highlights: • Unconventional collective excitations across the superfluid–Mott glass transition. • Higgs mode localized in superfluid and insulating phases. • Broad uncritical amplitude response function at quantum phase transition. • Goldstone mode delocalizes in the superfluid phase. We develop an inhomogeneous quantum mean-field approach to the behavior of collective excitations across the superfluid–Mott glass quantum phase transition in two dimensions, complementing recent quantum Monte Carlo simulations (Puschmann et al. 2020). In quadratic (Gaussian) approximation, the Goldstone (phase) and Higgs (amplitude) modes completely decouple. Each is described by a disordered Bogoliubov Hamiltonian which can be solved by an inhomogeneous multi-mode Bogoliubov transformation. We find that the Higgs mode is spatially localized in both phases. The corresponding scalar spectral function shows a broad peak that is noncritical in the sense that its peak frequency does not soften but remains nonzero across the quantum phase transition. In contrast, the lowest-energy Goldstone mode delocalizes in the superfluid phase, leading to a zero-frequency spectral peak. We compare these findings to the results of the quantum Monte Carlo simulations. We also relate them to general results on the localization of bosonic excitations, and we discuss the limits and generality of our approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1982

Irradiance measurement errors due to the assumption of a Lambertian reference panel

A technique is presented for determining the error in diurnal irradiance measurements that results from the non-Lambertian behavior of a reference panel under various irradiance conditions. Spectral biconical reflectance factors of a spray-painted barium sulfate panel, along with simulated sky radiance data for clear and hazy skies at six solar zenith angles, were used to calculate the estimated panel irradiances and true irradiances for a nadir-looking sensor in two wavelength bands. The inherent errors in total spectral irradiance (0.68 microns) for a clear sky were 0.60, 6.0, 13.0, and 27.0% for solar zenith angles of 0, 45, 60, and 75 deg, respectively. The technique can be used to characterize the error of a specific panel used in field measurements, and thus eliminate any ambiguity of the effects of the type, preparation, and aging of the paint.

Kimes, D. S.↗
DOE OSTI2019

Adaptive Reweighted Variance Estimation for Monte Carlo Eigenvalue Simulations

Monte Carlo (MC) simulation is used to solve the eigenvalue form of the Boltzmann transport equation to estimate various parameters such as fuel pin flux distributions that are crucial for the safe and efficient operation of nuclear systems (e.g., a nuclear reactor). Monte Carlo eigenvalue simulation uses a sample mean over many stationary cycles (iterations) to estimate various parameters important to nuclear systems. A variance estimate of the sample mean is often used for calculating the confidence intervals. However, MC eigenvalue simulation variance estimators that ignore the intercycle correlation underestimate the true variance of the estimated quantity. This paper presents novel data-adaptive approaches based on a simple autoregressive (AR) model and sigmoid functions to improve MC variance estimation. The standard MC sample-based variance estimator (or naïve estimator) and the spectral density–based MC variance estimator are enhanced by adding data-adaptive components that reduce their bias and improve performance. By investigating the frequency pattern of the AR(1) (order 1) model, two adaptive spectral estimators and one adaptive naïve estimator are proposed. The proposed estimators manifest superior performance when applied to three test problems compared to the standard spectral density–based estimator previously introduced by the authors. In conclusion, these new estimators are straightforward, as they use online algorithms and do not require storage of tallies from all active cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

CFD Evaluation of Pressure Change Along Coolant Passages in Sodium-Cooled Fast Reactor with Nek5000

To support the design efforts of advanced sodium-cooled fast reactors (SFRs), a series of computational fluid dynamics (CFD) simulations are performed to investigate the pressure change along various flow passages in the proposed SFR system. The simulations are carried out with the state-of-the-art spectral element flow solver, Nek5000. Two specific case studies are presented in this paper: the flow exiting the axial neutron reflector channels and the flow entering the fuel pin bundle. Due to the high Reynolds numbers expected, a Reynolds-averaged Navier-Stokes (RANS) approach is necessary to model the turbulence. A newly developed regularized k-omega RANS model is adopted in the related CFD calculations. The first case study explores the effect of Reynolds number on the pressure change when flow exits the reflector channels. The pressure change in this case has two major contributors: the change due to wall friction and the Bernoulli effect. It is noted that the nondimensional pressure loss follows a log-linear trend up to Re = 10 5 , and then the trend is flattened. In the second case study, the advanced NekNek coupling capability is tested where an integral domain can be divided into multiple subdomains with coupling interfaces, which would greatly ease the meshing process of complex engineering geometries and potentially save computational resources. The preliminary results obtained so far confirm the consistency between the NekNek results and those produced by regular Nek5000 simulation. The presented work demonstrates the readiness and flexibility of the related CFD techniques, which is part of the broader effort to leverage cutting-edge CFD to inform the advanced nuclear reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

A NEW MULTI-BAND RADIAL VELOCITY TECHNIQUE FOR DETECTING EXOPLANETS AROUND ACTIVE STARS

The radial velocity (RV) technique is one of the most efficient ways of detecting exoplanets. However, large RV jitters induced by starspots on an active star can inhibit detection of any exoplanet present or even lead to a false positive detection. This paper presents a new multi-band RV technique capable of substantially reducing starspot-induced RV jitters from stellar RV measurements to allow efficient and accurate extraction of RV signals caused by exoplanets. It takes full advantage of the correlation of RV jitters at different spectral bands and the independence of exoplanet signals at the corresponding bands. Simulations with a single-spot model and a multi-spot model have been conducted to investigate the RV jitter reduction capability of this method. The results show that this method can reduce the RV jitter amplitude by at least an order of magnitude, allowing detection of weaker exoplanet signals without significantly increasing RV observation time and cadence. This method can greatly reduce the observation time required to detect Earth-like planets around solar type stars with {approx}0.1 m s{sup -1} long term Doppler precision if spot-induced jitter is the dominant astrophysical noise source for RV measurements. This method can work efficiently for RV jitter removal if: (1) all the spots on a target star have approximately the same temperature during RV observations; (2) the RV jitter amplitude changes with wavelength, i.e., the RV jitter amplitude ratio, {alpha}, between two different spectral bands is not close to one; (3) the spot-induced RV jitter dominates the RV measurement error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SPECTRAL VARIABILITY FROM THE PATCHY ATMOSPHERES OF T AND Y DWARFS

Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present models that include patchy salt and sulfide clouds as well as water clouds for the Y dwarfs. We find that for objects over 375 K, patchy cloud opacity would generate the largest amplitude variability within near-infrared spectral windows. For objects under 375 K, water clouds also become important and generate larger amplitude variability in the mid-infrared. We also present models in which we perturb the temperature structure at different pressure levels of the atmosphere to simulate hot spots. These models show the most variability in the absorption features between spectral windows. The variability is strongest at wavelengths that probe pressure levels at which the heating is the strongest. The most illustrative types of observations for understanding the physical processes underlying brown dwarf variability are simultaneous, multi-wavelength observations that probe both inside and outside of molecular absorption features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Properties of laser-produced GaAs plasmas measured from highly resolved X-ray line shapes and ratios

The properties of hot, dense plasmas generated by the irradiation of GaAs targets by the Titan laser at Lawrence Livermore National Laboratory were determined by the analysis of high resolution K shell spectra in the 9 keV to 11 keV range. The laser parameters, such as relatively long pulse duration and large focal spot, were chosen to produce a steady-state plasma with minimal edge gradients, and the time-integrated spectra were compared to non-LTE steady state spectrum simulations using the FLYCHK and NOMAD codes. The bulk plasma streaming velocity was measured from the energy shifts of the Ga He-like transitions and Li-like dielectronic satellites. The electron density and the electron energy distribution, both the thermal and the hot non-thermal components, were determined from the spectral line ratios. After accounting for the spectral line broadening contributions, the plasma turbulent motion was measured from the residual line widths. The ionization balance was determined from the ratios of the He-like through F-like spectral features. The detailed comparison of the experimental Ga spectrum and the spectrum simulated by the FLYCHK code indicates two significant discrepancies, the transition energy of a Li-like dielectronic satellite (designated t) and the calculated intensity of a He-like line (x), that should lead to improvements in the kinetics codes used to simulate the X-ray spectra from highly-charged ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Self-radiography of imploded shells on OMEGA based on additive-free multi-monochromatic continuum spectral analysis

Radiographs of pure-DT cryogenic imploding shells provide critical validation of progress toward ignition-scalable performance of inertial confinement fusion implosions. Cryogenic implosions on the OMEGA Laser System can be self-radiographed by their own core spectral emission near ≈2 keV. Utilizing the distinct spectral dependences of continuum emissivity and opacity, the projected optical-thickness distribution of imploded shells, i.e., the shell radiograph, can be distinguished from the structure of the core emission distribution in images.Importantly, this can be done without relying on spectral additives (shell dopants), as in previous applications of implosion self-radiography. Furthermore, demonstrations with simulated data show that this technique is remarkably well-suited to cryogenic implosions and can also be applied to self-radiography of imploded room-temperature CH shells at higher spectral energy (hv ≈ 3–5 keV) based on the very similar continuum spectrum of carbon. Experimental demonstration of additive-free self-radiography with warm CH shell implosions on OMEGA will provide an important proof of principle for future applications to cryogenic DT implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inferring the temperature profile of the radiative shock in the COAX experiment with shock radiography, Dante, and spectral temperature diagnostics

Predicting and modeling the behavior of experiments with radiation waves propagating through low-density foams require a detailed quantification of the numerous uncertainties present. In regimes where a prominent radiative shock is produced, key dynamical features include the shock position, temperature, and curvature and the spatial distribution and temperature of the corresponding supersonic radiation wave. The COAX experimental platform is designed to constrain numerical models of such a radiative shock propagating through a low-density foam by employing radiography for spatial and shock information, Dante for characterizing the x-ray flux from the indirectly driven target, and a novel spectral diagnostic designed to probe the temperature profile of the wave. In this work, we model COAX with parameterized 2D simulations and a Hohlraum-laser modeling package to study uncertainties in diagnosing the experiment. The inferred temperature profile of the COAX radiation transport experiments has been shown to differ from simulations more than expected from drive uncertainties that have been constrained by simultaneous soft x-ray flux and radiography measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Simulated Space Environmental Testing on Thin Films

An exploratory program has been conducted, to irradiate some mature commercial and some experimental polymer films with radiation simulating certain Earth orbits, and to obtain data about the response of each test film's reflective and tensile properties. Protocols to conduct optimized tests were considered and developed to a "prototype" level during this program. Fifteen polymer film specimens were arranged on a specially designed test fixture. The fixture featured controlled exposure areas, and protected the ends of the samples for later gripping in tensile tests. The fixture featured controlled exposure areas, and protected the ends of the samples for later gripping in tensile tests. The fixture containing the films was installed in a clean vacuum chamber where protons, electrons and solar ultraviolet (UV) radiation could simultaneously irradiate the specimens. Near realtime UV rates were used, whereas proton and electron rates were accelerated appreciably to simulate 5 years in orbit during a two month test. Periodically, the spectral reflectance of each film was measured in situ. After the end of the irradiation, final reflectance measurements were made in situ, and solar absorptance values were derived for each specimen. These samples were then measured in air for thermal emittance and for tensile strength. Most specimens withstood the irradiation intact, but with reduced reflectance (increased solar absorptance). Thermal emittance changed slightly in several materials, as did their tensile strength and elongation at break. Conclusions are drawn about the performance of the films. Simulated testing to an expected 5 year dose of electrons and protons consistent with those expected at L2 and 0.98 AU orbits and 100 equivalent solar hours exposure.

Russell, Dennis A.↗
DOE OSTI2016

The Symmetrical Quasi-Classical Model for Electronically Non-Adiabatic Processes Applied to Energy Transfer Dynamics in Site-Exciton Models of Light-Harvesting Complexes

In a recent series of papers, it has been illustrated that a symmetrical quasi-classical (SQC) windowing model applied to the Meyer–Miller (MM) classical vibronic Hamiltonian provides an excellent description of a variety of electronically non-adiabatic benchmark model systems for which exact quantum results are available for comparison. In this study, the SQC/MM approach is used to treat energy transfer dynamics in site-exciton models of light-harvesting complexes, and in particular, the well-known 7-state Fenna–Mathews–Olson (FMO) complex. Again, numerically “exact” results are available for comparison, here via the hierarchical equation of motion (HEOM) approach of Ishizaki and Fleming, and it is seen that the simple SQC/MM approach provides very reasonable agreement with the previous HEOM results. Finally, it is noted, however, that unlike most (if not all) simple approaches for treating these systems, because the SQC/MM approach presents a fully atomistic simulation based on classical trajectory simulation, it places no restrictions on the characteristics of the thermal baths coupled to each two-level site, e.g., bath spectral densities (SD) of any analytic functional form may be employed as well as discrete SD determined experimentally or from MD simulation (nor is there any restriction that the baths be harmonic), opening up the possibility of simulating more realistic variations on the basic site-exciton framework for describing the non-adiabatic dynamics of photosynthetic pigment complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗