SEARCH · Engineering Papers

Results for “solar constant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Characterizing Current THD’s Dependency on Solar Irradiance and Supraharmonics Profiling for a Grid-Tied Photovoltaic Power Plant

The rapidly increasing distributed energy resources (DERs) in power systems are now getting interconnected to set community grid structures, where power quality will be a major concern. The grid-to-grid (G2G) bidirectional power transfer among the distribution microgrid will not be considered commercially feasible unless the upstream harmonics are under the limits. The aggregation of such harmonics, measured as total harmonic distortion (THD), is feared to be beyond tolerable limits with the progression of rooftop grid-tied PV-like installations. Hence, this THD needs to be characterized with DER generation end variables. In this work, the photovoltaic (PV) DERs’ dependency on environment variables such as irradiance was profiled in the context of generating and injecting harmonics into the grid. A mathematical model of a grid-tied three-phase PV DER was developed as part of this correlation characterization, matching the fundamental unit structure of a 1.4 MW solar canopy located on the Florida International University (FIU) Miami campus. To determine the qualitative association with produced THD patterns, the model was evaluated with various irradiance settings. A real-time digital simulation (RTDS) platform was used to verify it. Following this confirmation, sets of data from power quality meters at the point of common coupling and FIU field sensors were utilized to validate further the correlation model. The results showed that the grid current’s THD exhibited a high correlation with the irradiance profile and its variation over time. The early morning and late afternoon periods of the day, associated with a low irradiance, constantly had higher harmonics generated from the PV DER. The midday THD was rather rational with partial shadings, hence a geolocation-dependent factor. These findings were verified by an RTDS and validated by real field data. In quantifying the THD injected by a single DER at a high-frequency (2–150 kHz) supraharmonics (SH) level, a 3% peak increment in magnitude was observed from the high-fixed to the low-fixed irradiance profile. The correlation characteristics depicted that the hybrid microgrid suffered from a daytime-dependent harmonic insertion from the grid-tied DER. This is a global problem unless specific measures are taken to mitigate the harmonics. The electrically notorious higher-frequency SH was found to increase proportionally. The G2G power transfer can be limited because of the higher THD in the early morning and late afternoon, which will also worsen because the numbers of grid-tied PV DERs (i.e., rooftop solar and industrial solar) are likely to increase rapidly soon. The community grid structure can thus have a controlled harmonics filtration setup purposefully designed to address the findings of this work, which also fall within the scope of our future research.

14 SOLAR ENERGY↗

DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Direct and Indirect Interfacial Electron Transfer at a Plasmonic p-Cu 7 S 4 /CdS Heterojunction

Plasmonic semiconductors exhibit significant potential for harvesting near-IR solar energy, although their mechanisms of plasmon-induced hot electron transfer (HET) are poorly understood. We report a transient absorption study of plasmon-induced HET in p-Cu 7 S 4 /CdS type II heterojunctions. Near-IR excitation of the p-Cu 7 S 4 plasmon band at ~1400 nm leads to ultrafast HET into the CdS conduction band with a time constant of <150 fs and a quantum efficiency of ~0.054%. The injected hot electrons remain in CdS with an amplitude-weighted average lifetime of 1.9 ± 0.5 ns, significantly longer than that in Au/CdS heterostructures, suggesting that plasmonic semiconductors can slow down charge recombination due to the presence of a bandgap. The excited near-IR plasmon does not decay by coupling to the interfacial charge transfer transition, likely due to its energy mismatch. This study provides a detailed mechanistic understanding and possible directions for improving plasmonic HET in plasmonic semiconductor heterojunctions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Asymmetric drift in MaNGA: mass and radially dependent stratification rates in galaxy discs

ABSTRACT We measure the age–velocity relationship from the lag between ionized gas and stellar tangential speeds in ∼500 nearby disc galaxies from MaNGA in Sloan Digital Sky Survey IV (SDSS-IV). Selected galaxies are kinematically axisymmetric. Velocity lags are asymmetric drift, seen in the Milky Way’s (MW) solar neighbourhood and other Local Group galaxies; their amplitude correlates with stellar population age. The trend is qualitatively consistent in rate ($\dot{\sigma }$) with a simple power-law model where σ ∝ tb that explains the dynamical phase-space stratification in the solar neighbourhood. The model is generalized based on disc dynamical times to other radii and other galaxies. We find in-plane radial stratification parameters σ0,r (dispersion of the youngest populations) in the range of 10–40 km s−1 and 0.2 < br < 0.5 for MaNGA galaxies. Overall, brincreases with galaxy mass, decreases with radius for galaxies above 10.4 dex (M⊙) in stellar mass, but is ∼constant with radius at lower mass. The measurement scatter indicates the stratification model is too simple to capture the complexity seen in the data, unsurprising given the many possible astrophysical processes that may lead to stellar population dynamical stratification. None the less, the data show dynamical stratification is broadly present in the galaxy population, with systematic trends in mass and density. The amplitude of the asymmetric drift signal is larger for the MaNGA sample than the MW, and better represented in the mean by what is observed in the discs of M31 and M33. Either typical discs have higher surface-density or, more likely, are dynamically hotter (hence thicker) than the MW.

Bershady, Matthew A. (ORCID:0000000231314374)↗

DOE OSTI2021

Change-point Detection and Image Segmentation for Time Series of Astrophysical Images

Many astrophysical phenomena are time-varying, in the sense that their intensity, energy spectrum, and/or the spatial distribution of the emission suddenly change. This paper develops a method for modeling a time series of images. Under the assumption that the arrival times of the photons follow a Poisson process, the data are binned into 4D grids of voxels (time, energy band, and x-y coordinates), and viewed as a time series of non-homogeneous Poisson images. The method assumes that at each time point, the corresponding multiband image stack is an unknown 3D piecewise constant function including Poisson noise. It also assumes that all image stacks between any two adjacent change points (in time domain) share the same unknown piecewise constant function. The proposed method is designed to estimate the number and the locations of all of the change points (in time domain), as well as all of the unknown piecewise constant functions between any pairs of the change points. The method applies the minimum description length principle to perform this task. A practical algorithm is also developed to solve the corresponding complicated optimization problem. Simulation experiments and applications to real data sets show that the proposed method enjoys very promising empirical properties. Applications to two real data sets, the XMM observation of a flaring star and an emerging solar coronal loop, illustrate the usage of the proposed method and the scientific insight gained from it.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗

DOE OSTI2023

Influence of Correlations on the Thermal Performance Modeling of Parabolic Trough Collectors

The influence of correlations on the thermal performance modeling of parabolic trough collectors was analyzed in this work. A versatile model for a parabolic trough collector was developed that allows one- and two-dimensional analysis and enables the use of correlations to calculate thermophysical properties and convection heat transfer coefficients. The model also allows the use of constant values for properties and/or coefficients obtained from the evaluation correlations at a specific temperature. The effect of each correlation was evaluated independently, and the results were compared with a reference case that considered a two-dimensional approach and used all the correlations. For the analyzed cases, the correlation for the absorber emittance has the strongest impact on the collector efficiency, leading to a lower error when used. Based on the results, a one-dimensional model approach considering a correlation for the absorber emittance leads to efficiency errors below 3% for collector lengths of up to 243.6 m. Compared with the reference case, a one-dimensional approach using all correlations for a collector with a length of 500 m, and operating with an inlet temperature of 773 K, can result in errors around 9%. However, using constant values for properties and heat transfer coefficients could lead to errors of up to 50%. Multiple thermal models for parabolic trough collectors proposed in the literature rely on a one-dimensional approach, estimated values for the heat transfer coefficients, and constant thermophysical properties. The errors associated with those approaches are analyzed and quantified in this work as a function of the collector length and operation temperature.

absorber emittance↗

DOE OSTI2024

Amidination of ligands for chemical and field-effect passivation stabilizes perovskite solar cells

Surface passivation has driven the rapid increase in the power conversion efficiency (PCE) of perovskite solar cells (PSCs). However, state-of-the-art surface passivation techniques rely on ammonium ligands that suffer deprotonation under light and thermal stress. We developed a library of amidinium ligands, of interest for their resonance effect–enhanced N–H bonds that may resist deprotonation, to increase the thermal stability of passivation layers on perovskite surfaces. This strategy resulted in a >10-fold reduction in the ligand deprotonation equilibrium constant and a twofold increase in the maintenance of photoluminescence quantum yield after aging at 85°C under illumination in air. Here, implementing this approach, we achieved a certified quasi–steady-state PCE of 26.3% for inverted PSCs; and we report retention of ≥90% PCE after 1100 hours of continuous 1-sun maximum power point operation at 85°C.

Yang, Yi↗

DOE OSTI2025

A Comparison of Battery Charge Controller Technologies for Wave Energy Converters

Wave energy is a uniquely challenging field for electrical system designers. High peak and low average power potential with a constantly varying energy input is difficult to harness and control through conventional means. To power the blue economy, low-powered wave energy converters (WECs) need batteries for energy storage. Safely and effectively charging batteries from waves requires a charge controller to properly monitor and control voltage and current going to the battery. Currently, off-the-shelf charge controllers exist for other renewable generation such as wind, hydro, and solar. Two topologies were validated: a buck converter and a pulse width modulation (PWM) charge controller. Using an in-lab dry testbed, wave energy power inputs were simulated to properly validate the effectiveness of existing charge controller technologies, identifying the shortcomings and improvements needed to effectively harness wave energy.

battery storage↗

DOE OSTI2024

A Comparison of Battery Charge Controller Technologies for Wave Energy Converters: Preprint

battery storage↗

DOE OSTI2022

Highly efficient organic solar cells with superior deformability enabled by diluting the small molecule acceptor content

Developing efficient organic solar cells (OSCs) with strong mechanical deformability is important to be addressed to ensure their operational reliability in wearable electronics. However, it is challenging to achieve mechanically robust polymer/small molecule OSCs with an efficiency over 17% due to the abundant brittle donor/acceptor (D/A) interface. Decreasing the small molecule content can reduce the brittle D/A interface area to enhance deformability, but it may also cause discontinuous electron-conducting regions and thus deteriorate photovoltaic performance. Here, we incorporate a polymer donor (D18) into a binary PTQ10cm-BTP-PhC6 system at a constant PTQ10:m-BTP-PhC6 ratio of 1 : 1.2 to minimize the D/A interfacial area and modulate phase separation, and subsequently fabricate mechanically reliable OSCs with high efficiency. The large incorporation of D18 molecules increases the crystallinity of m-BTP-PhC6 molecules because of the poor interaction between D18 and the two host materials, affording sufficient pathways for efficient charge percolation at a large D/A ratio. More importantly, highly aggregated m-BTP-PhC6 molecules help to minimize the D/A interfacial area, enabling ternary films with better deformability than their binary counterparts. As a result, ternary OSCs with a large D/A ratio of 1.5 : 1.2 exhibit a high efficiency of 17.3% with a crack-onset strain (COS) of 8.8%, which significantly outperforms their binary counterparts with an efficiency of 13.8% and a COS of 4.4%. Finally, diluting the content of small molecules and controlling intermolecular interaction through the introduction of the second polymer donor represents a promising strategy to fabricate efficient OSCs with superior deformability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Modulation of Photoinduced Iodine Expulsion in Mixed Halide Perovskites with Electrochemical Bias

Hole trapping at iodine (I) sites in MAPbBr 1.5 I 1.5 mixed halide perovskites (MHP) is responsible for iodine migration and its eventual expulsion into solution. We have now modulated the photoinduced iodine expulsion in MHP through an externally applied electrochemical bias. At positive potentials, electron extraction at TiO 2 /MHP interfaces becomes efficient, leading to hole buildup within MHP films. This improved charge separation, in turn, favors iodine migration as evident from the increased apparent rate constant of iodine expulsion (k expulsion = 0.0030 s –1 ). Conversely, at negative potentials (-0.3 V vs Ag/AgCl) electron–hole recombination is facilitated within MHP, slowing down iodine expulsion by an order of magnitude (k expulsion = 0.00018 s –1 ). The tuning of the EFermi level through external bias modulates electron extraction at the TiO 2 /MHP interface and indirectly controls the buildup of holes, ultimately inducing iodine migration/expulsion. Suppressing iodine migration in perovskite solar cells is important for attaining greater stability since they operate under internal electrical bias.

14 SOLAR ENERGY↗

DOE OSTI2021

Conjugated-linker dependence of the photophysical properties and electronic structure of chlorin dyads

The synthesis, photophysical properties and electronic structure of seven new chlorin dyads and associated benchmark monomers are described. Each dyad contains two identical chlorins linked at the macrocycle β-pyrrole 13-position. The extent of electronic communication between chlorin constituents depends on the nature of the conjugated linker. The communication is assessed by modification of prominent ground-state absorption and redox properties, rate constants and yields of excited-state decay processes, and molecular-orbital characteristics. Relative to the benchmark monomers, the chlorin dyads in toluene exhibit a substantial bathochromic shift of the long-wavelength absorption band (30 nm average), two-fold increased radiative rate constant [average (10 ns) -1 vs. (22 ns) -1 ], reduced singlet excited-state lifetimes (average 5.0 ns vs. 8.2 ns), and increased fluorescence quantum yields (average 0.56 vs. 0.42). The excited-state lifetime and fluorescence yield for the chlorin dyad with a benzothiadiazole linker are reduced substantially in benzonitrile vs. toluene due largely to ~25-fold accelerated internal conversion. The results aid design strategies for molecular architectures that may find utility in solar-energy conversion and photomedicine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Reorganization Energy of Electron Transfer in Ionic Liquids

We report bandshape analysis of charge-transfer optical bands in room-temperature ionic liquids (ILs) is performed to extract reorganization energy of electron transfer. Remarkably, the reorganization energies in ILs are close to those in cyclohexane. This result runs against common wisdom in the field since conducting ILs, characterized by an infinite static dielectric constant, and nonpolar cyclohexane fall to the opposite ends of the polarity scale based on their dielectric constants. Theoretical calculations employing structure factors of ILs from molecular dynamics simulations support low values of the reorganization energy. Standard dielectric arguments do not apply to solvation in ILs and nonergodic reorganization energies are required for a quantitative analysis.

14 SOLAR ENERGY↗

DOE OSTI2023

Symmetry-driven persistent spin texture for the two-dimensional nonsymmorphic CdTe and ZnTe crystal structures

In this paper, two nonsymmorphic two-dimensional structures of CdTe and ZnTe are modeled, and using state-of-the-art density functional theory with the group theory of solids, their symmetry-enforced electronic properties are studied. The in-plane ferroelectricity coupled with strong spin-orbit coupling induces a unidirectional out-of-plane Rashba spin-orbit field that can host a momentum-independent uniform spin configuration known as persistent spin texture (PST) at the Brillouin zone center. PST in these structures is found to be robust against external perturbations such as strain, structural distortion, and independent of layer thickness. These unprecedented intrinsic spin transport properties hold utmost importance in spintronics, as the experimental stringent condition of equal Rashba and Dresselhaus constants [Phys. Rev. Lett. 90, 146801 (2003)] is eliminated. The calculated persistent spin helix wavelength of <~5 nm paves the way for developing next-generation nanosized nonballistic spin field-effect transistors compared with micrometer-sized GaAs/AlGaAs quantum wells. Further, these materials exhibit finite spin Hall conductivity at the band edges and hence can be used in ferromagnet-free spin Hall transistors. Although CdTe and ZnTe systems have been widely studied for photocatalysis and solar cell applications over the past few decades, their potential application in spintronic devices has not been explored. Mono/few layers of CdTe and ZnTe synthesized from (110) facets of bulk zinc-blende crystals [Nat. Commun. 3, 1057 (2012)] satisfy all symmetry operations of the nonsymmorphic space group and hence can be considered ideal materials to verify our theoretical results experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Spatial profiles of photon chemical potential in near-field thermophotovoltaic cells

Emitted photons stemming from the radiative recombination of electron–hole pairs carry chemical potential in radiative energy converters. This luminescent effect can substantially alter the local net photogeneration in near-field thermophotovoltaic cells. Several assumptions involving the luminescent effect are commonly made in modeling photovoltaic devices; in particular, the photon chemical potential is assumed to be zero or a constant prescribed by the bias voltage. The significance of photon chemical potential depends upon the emitter temperature, the semiconductor properties, and the injection level. Hence, these assumptions are questionable in near-field thermophotovoltaic devices with nanoscale separation distances between the emitter and the cell. In the present work, an iterative solver that combines fluctuational electrodynamics with the drift-diffusion model is developed to tackle the coupled photon and charge transport problem, enabling the determination of the spatial profile of photon chemical potential beyond the detailed balance approach. Here, the difference between the results obtained by allowing the photon chemical potential to vary spatially and by assuming a constant value demonstrates the limitations of the conventional approaches. This study is critically important for performance evaluation of near-field thermophotovoltaic systems with nanoscale vacuum gaps.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2025

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗