Understanding the fundamental limits of various energy conversion approaches is essential for assessing their efficiency and power output. In this work, we derive general performance bounds for thermoradiative heat engines that exchange heat radiatively with a cold sink, establishing power-versus-efficiency thermodynamic bounds for several configurations. We find that the performance of these engines is always bounded by that of radiative engines, which harness the thermal radiation emitted by a hot source, making thermoradiative engines inherently less favorable for energy conversion. By unifying the results of radiative and thermoradiative engines within a common thermodynamic framework, which also encompasses dual-engine configurations that combine both features, this work provides a comprehensive understanding of the thermodynamic limits of radiative energy conversion. Our framework offers general metrics for evaluating specific devices and raises critical questions regarding the relevance of thermoradiative cells for energy production.

Concerningly, numerous incorrect expressions for extracting phonon lifetimes from spectroscopic linewidths remain prevalent in the literature, resulting in misinterpretation of experimental data. This has significant consequences in the broader energy materials field, where erroneous lifetimes risk misinforming materials design and device modeling. We outline the origin of this ambiguity, provide historical context on how this error has propagated, and present the correct formula. In addition, we highlight a frequently overlooked source of error: the contribution of phase relaxation to spectral broadening. A unified approach enables fair comparison of values obtained from different spectroscopic methods, accurate benchmarking against first-principles calculations, and, most importantly, precise determination of macroscopic physical properties derived from phonon lifetime measurements.

In recent years, orientationally disordered crystals, or plastic crystals, have transformed the field of solid-state cooling due to the significant latent heat and entropy changes associated with their temperature-induced molecular order-disorder phase transition, which can produce colossal caloric effects under external field stimuli. However, the molecular mechanisms underlying these huge caloric effects remain inadequately understood and general principles for enhancing the performance of caloric plastic crystals are lacking. Previous studies have predominantly focused on molecular rotations, overlooking other potentially critical factors, such as lattice vibrations and molecular conformations. In this study, we employ classical molecular dynamics (MD) simulations to both replicate and elucidate the microscopic origins of the experimentally observed colossal barocaloric (BC) effects—those driven by hydrostatic pressure—in the archetypal plastic crystal neopentyl glycol (NPG). Our MD simulations demonstrate that in NPG, the combined BC response and phase-transition entropy changes arising from lattice vibrations and molecular conformations are nearly equal to those from molecular reorientations. These findings suggest that alongside hydrogen bonding—which directly impacts molecular rotational dynamics—lattice vibrational and molecular structural features, often overlooked, must be integrated into the rational design and modeling of advanced caloric plastic crystals. These insights are not only of significant fundamental interest but also essential for driving the development of next-generation solid-state refrigeration technologies.

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

The migration of mobile ions is one of the leading causes of the degradation of perovskite solar cells. However, quantifying mobile ions in complete perovskite solar cells is challenging due to the complex device stacks and the impact of charge transport layers on the measurement techniques. Here we develop a simple and openly accessible step model that approximates drift-diffusion simulations. The step model is based on the expression of the charge density in the ionic and electronic accumulation and depletion layers as step functions. We can then accurately determine the impact of mobile ions on the dc potential distribution of perovskite solar cells. Furthermore, we can simulate electrical measurement techniques commonly used to quantify mobile ions: capacitance transient, current transient, and capacitance-frequency measurements. By validating the step model with drift-diffusion simulations, we show that an accurate extraction of ion density, diffusion coefficient, and activation energy is possible in an accessible range. We finally apply the step model developed to estimate the ionic conductivity and activation energy of perovskite solar cells.

This study develops a high-fidelity multiphysics coupling framework to systematically investigate the impact of Chalk River unidentified deposits (CRUD) on pressurized-water-reactor (PWR) assemblies, integrating neutronics, thermal hydraulics, and chemical kinetics modeling. CRUD, formed via subcooled nucleate boiling (SNB) and corrosion-product deposition, alters local boron concentrations, thereby inducing a power shift (CIPS) due to CRUD and generating asymmetric axial power distributions. The proposed framework adopts an internally coupled approach, utilizing Picard iteration and a background-mesh technique to resolve multiscale computational challenges while maintaining numerical accuracy and efficiency. Numerical simulations of a $17\times 17$ PWR fuel assembly reveal that CRUD deposition reduces the eigenvalue by 293 percent mille (pcm) and shifts the axial power peak downward by 21 cm, corresponding to a 6.27% axial offset. Radially, the effect of CRUD on power distribution remains minimal (approximately 0.04%). A critical observation is the dominant influence of surface heat flux on CRUD growth, with thicker deposits and elevated boron concentrations observed in high-power regions. CRUD-induced perturbations also exhibit a “neighboring phenomenon,” wherein CRUD-free fuel rods experience comparable axial power shifts due to neutron-flux disturbances. These findings underscore the intricate coupling between power distribution, SNB, and CRUD chemistry, providing valuable insights for enhancing reactor performance and safety.

The p-i-n junction structure that develops dynamically under an applied bias via electrochemical doping (ECD) is of key importance for the performance of light-emitting electrochemical cells (LECs). While the complex electronic and ionic processes that govern its transient formation have been extensively studied by both experiments and drift-diffusion modeling, less attention has been given to the steady-state junction as a function of voltage. Here we study the formed p-i-n structure of a polymer LEC at the steady state by measuring and analyzing its most distinctive feature: the current density–voltage-luminance characteristics. Unexpectedly, we find that the effective conductance of the p-i-n structure exhibits a positive correlation with the applied bias, a behavior not predicted by existing LEC drift-diffusion models. We attribute this discrepancy to the assumption in these models of a constant density of mobile ions. Hence, we present a modified model in which the ECD level—represented by the number of ions with which the organic semiconductor is doped—scales with the applied voltage, implying a voltage-dependent doping efficiency. We validate this hypothesis using electron spin resonance spectroscopy and drift-diffusion modeling, while additionally establishing that only a small fraction of the available ions in our system, which increases from 1% to 3% with increasing bias, contributes to the ECD and is necessary for efficient LEC operation. These findings not only provide fundamental insights into the operational mechanism of LECs but also have direct implications for the broader organic mixed ionic and electronic conductor community.

The development of clean energy solutions requires an extension of electrical grids in size and power. Network evolution modeling is crucial for revealing critical insights into transmission grid growth dynamics. This study proposes an efficient analytical framework based on spatial network physics to enhance previous models and deliver better predictions. Utilizing the French transmission grid as a case study, our methodology demonstrates its aptitude to estimate energy invested and consumed, as well as topological and spatial characteristics. This new approach balances computational efficiency with predictive accuracy, making it particularly valuable for building datasets for research and contributing to well-informed infrastructure development forecast science.

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., $O({10}^2)\mathrm{GPa}$, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the $O({10}^0)\mathrm{GPa}$ and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the $\mathrm{N}$-doped $\mathrm{Lu}$ hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic ${\text{Lu(H,N,Va)}}_3$ system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., $\mathrm{N}$-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the $\mathrm{Lu}$-$\mathrm{H}$-$\mathrm{N}$ system and help rational synthesis of $\mathrm{N}$-doped $\mathrm{Lu}$ hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

We perform a high-throughput computational search for novel phonon-mediated superconductors, starting from the Materials Cloud three-dimensional structure database of experimentally known inorganic stoichiometric compounds. We first compute the Allen-Dynes critical temperature ($T_c^{\mathrm{AD}}$) for $4533$ nonmagnetic metals using a direct and progressively finer sampling of the electron-phonon couplings. For the candidates with the largest $T_c^{\mathrm{AD}}$ value, we use automated Wannierizations and electron-phonon interpolations to obtain a high-quality data set for the most promising $250$ dynamically stable structures, for which we calculate spectral functions, superconducting band gaps, and isotropic Migdal-Eliashberg critical temperatures. For $140$ of these, we also provide anisotropic Migdal-Eliashberg superconducting gaps and critical temperatures. The approach is remarkably successful in finding known superconductors and we find $24$ unknown ones with a predicted anisotropic $T_c$ value above 10 K. Among them, we identify a possible double-gap superconductor ($p$-doped $\mathrm{Ba}\mathrm{B}$${}_2$), a nonmagnetic half-Heusler $\mathrm{Zr}\mathrm{Ru}\mathrm{Sb}$, and the perovskite $\mathrm{Ta}\mathrm{Ru}$${}_3$$\mathrm{C}$, all exhibiting significant $T_c$ values. Finally, we introduce a sensitivity analysis to estimate the robustness of the predictions.

The interplay between photochemical and electrochemical reactions fundamentally influences charge transfer processes at solid-liquid interfaces. Nevertheless, chemical processes at semiconductor surfaces triggered by light excitation under an applied potential remain poorly explored. This work deciphers the synergistic effect of potential and light excitation on boron-doped diamond electrodes in producing either surface redox reactions or emission of solvated electrons in water. The effect of diamond surface termination on electron affinity, band bending, and charge extraction is identified in a photoelectrochemical cell. While photocurrent is observed for excitation as low as 3.5 eV, we show that it is induced mostly by surface redox reactions, whereas solvated electrons are detected only for excitation above the band gap (5.47 eV). Solvated electrons are generated irrespective of band bending, which affects only the emission yield. Depending on the surface band bending, photoreduction of the hydroxylated surface groups and photooxidation of the C—$\mathrm{H}$ surface groups can be induced by direct photoexcitation in the range of 4.2–4.8 eV. The surface of the diamond can be electrochemically reduced when the Fermi level of the oxidized surface decreases below the ${\mathrm{H}}^{+}$/${\mathrm{H}}_2$ redox potential. On the other hand, the hydrogenated surface oxidizes spontaneously for potentials at which the Fermi level drops below the occupied $\text{C---H}$ surface states, depending on both the pH and the electron affinity of the surface. This work provides insights into (photo)redox processes on diamond materials, which may find applications in photoelectrochemical solar fuel generation or energy storage.

Lithium thiophosphate (LPS) is a promising solid electrolyte for next-generation lithium-ion batteries due to its superior energy storage, high ionic conductivity, and low-flammability properties. Despite its potential, the high reactivity of LPS with common contaminants such as atmospheric water, preparation solvents, and electrode materials poses significant challenges for commercialization. The lack of understanding regarding the structure, morphology, and chemical behavior of the surface of LPS slows down the search for solutions to these issues. Here, we utilize a machine-learning interatomic potential to achieve a fundamental atomistic understanding of the mechanical and chemical properties of the $\beta$-$\mathrm{Li}$${}_3$$\mathrm{PS}$${}_4$ surfaces. Employing molecular-dynamics simulations, we identify relevant complexions formed by surface reconstructions, determine their surface energies, and compute the Wulff shape of $\beta$-LPS. The most stable complexions exhibit properties distinctly different from the bulk, including amorphization, increased density, decreased conductivity, and large deformation of the structural building blocks. We demonstrate that these surfaces are not static but undergo significant dynamical activity at room temperature, which is clearly identified by an analysis featuring a time-averaged structural descriptor. Finally, we examine the changes of the electronic structure induced by the surface complexions as well as their effect on the reaction enthalpy with water, which provides us with details on changes in surface reactivity and active sites, underlining the importance of investigating surface complexions under realistic conditions.

Point defects are ubiquitous in materials and play a critical role in determining their electrical, optical, ion-transport, and other properties, making them key to the performance of materials in energy applications. First-principles defect calculations, particularly those based on density functional theory, are powerful tools for predicting defect energetics and associated electronic properties. However, the results are often communicated through defect formation energy diagrams and defect level diagrams, which can be challenging to interpret for researchers who do not specialize in defect physics. This Tutorial provides a practical accessible guide for understanding and interpreting these diagrams, with a focus on materials with a band gap—semiconductors and insulators. Aimed at experimentalists and computationalists new to defect analysis, the Tutorial explains how to extract key physical insights from defect diagrams, such as defect and carrier concentrations, dopability, midgap states, and charge localization. By demystifying the language and visuals of defect calculations, this guide seeks to bridge the gap between theory and experiments, enabling broader use of defect data in the rational design of functional materials for energy applications and beyond.

We investigate the layered magnetic semiconductor ${\mathrm{Eu}\mathrm{Te}}_4$, focusing on its intricate charge-density-wave (CDW) states near and above room temperature through single-crystal x-ray diffraction (XRD), magnetic, and thermoelectric measurements. The XRD measurement revealed that the CDW state inducing the polar lattice distortion persists even at 650 K, demonstrating its remarkable thermal stability. Notably, the Seebeck coefficient near room temperature reaches values exceeding 500 µV K${}^{-1}$. This large Seebeck coefficient, not fully captured by a simple band calculation, is comparable to those observed in heavy-electron semiconducting oxides, suggesting the importance of electron correlation and spin- and lattice instabilities. Furthermore, potentially reflecting the competition of two types of CDW state, the thermal conductivity near room temperature is as low as 0.02 W cm${}^{-1}$ K${}^{-1}$. As a result, the thermoelectric figure of merit zT reaches 0.22 at 460 K. These findings establish ${\mathrm{Eu}\mathrm{Te}}_4$ as a compelling platform to explore novel types of thermoelectric materials with multiple electronic instability.

Magnesium-ion batteries (MIBs) are promising next-generation energy storage devices. Vanadium oxides are promising cathode materials for MIBs due to their rich chemical valences and versatile structures. However, vanadium oxide cathodes face challenges such as low structural stability, slow ${\mathrm{Mg}}^{2+}$ diffusion kinetics, and low electrical conductivity. The work reported herein used hydrothermal synthesis to pre-intercalate ${\mathrm{H}}_2\mathrm{O}$ molecules and ${\mathrm{Mg}}^{2+}$ into ${\mathrm{V}}_2{\mathrm{O}}_5$ interlayers to stabilize the lamellar structure and modulate $\mathrm{V}$ 3d orbital splitting, leading to a ${\mathrm{Mg}}_{0.2}{\mathrm{V}}_2{\mathrm{O}}_5\cdot 1.5{\mathrm{H}}_2\mathrm{O}$ cathode (MVOH1). The $\mathrm{V}$ 3d orbital splitting facilitates electron hopping between ${\mathrm{V}}^{4+}$ and ${\mathrm{V}}^{5+}$ cations, leading to the enhanced conductivity of MVOH1; it also makes it easier for $\mathrm{V}$ atoms to acquire and lose electrons during discharging and charging, leading to high capacity and high rate capability of MVOH1. Consequently, the MVOH1 cathode exhibited a high specific capacity of 300 mA h ${\mathrm{g}}^{-1}$ at a current density of 0.05 A ${\mathrm{g}}^{-1}$ and retained a specific capacity of 85 mA h ${\mathrm{g}}^{-1}$ at 10 A ${\mathrm{g}}^{-1}$; the capacity and rate capability surpass those of previously reported ${\mathrm{V}}_2{\mathrm{O}}_5$ cathode materials. The aqueous MIB based on MVOH1 demonstrated an ultralong lifespan of 11 000 cycles, which is significantly longer than that of previously reported aqueous MIBs. This study highlights the synergistic effect of regulating interlayer spacing and $\mathrm{V}$ 3d orbital splitting on electrochemical performance, which is crucial for developing more efficient and more durable energy storage systems.

Simultaneously achieving interfacial stability and high room-temperature ionic conductivity ($\ge {10}^{-3}$ S/cm) in solid-state electrolytes is one of the most significant challenges in the development of all-solid-state batteries. In this study, by combining first-principles simulations and large-scale machine-learning molecular dynamics modeling, we demonstrate ${\mathrm{K}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$ is such an outstanding electrolyte for $\mathrm{K}$-ion batteries that showcases both intrinsic interfacial stability with metal anode and fast ion diffusivity. Conversely, ${\mathrm{Li}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$ only exhibits mediocre performance for $\mathrm{Li}$-ion batteries. Our results show ${\mathrm{Li}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$ is thermodynamically unstable and will be reduced by the metal anode, whereas thermodynamically stable ${\mathrm{K}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$ presents a wide electrochemical stability window of 0.00–3.27 V versus $\mathrm{K}/{\mathrm{K}}^{+}$, suggesting it is stable in contact with potassium metal. The predicted ionic conductivity is 0.43 and 3.63 mS/cm for ${\mathrm{Li}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$ and ${\mathrm{K}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$, respectively. Moreover, a significant size effect has been observed in ${\mathrm{K}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$, where the accurate ionic conductivity can only be obtained through large-scale simulations. Further analysis reveals low energy barrier and long hopping distance are responsible for the higher ionic conductivity in ${\mathrm{K}}_5{\mathrm{YSi}}_4{\mathrm{O}}_{12}$. This work identifies an excellent electrolyte candidate for $\mathrm{K}$-ion batteries and underscores the importance of scale in the modeling of ionic conductivity.

The radiative open-circuit voltage loss in a solar cell occurs because the absorptance spectrum near the band gap shows a gradual increase rather than a sharp step-function-like transition. This broadening effect has been attributed to band-gap fluctuations and/or to Urbach tails. In this work, we use modeling based on Planck’s generalized law to distinguish between these two effects. Our results demonstrate that Urbach tails have only a minimal effect on the absorptance-edge broadening and clarify that even an ideal direct semiconductor, with no band-gap fluctuations, shows broadening at the absorptance onset. Furthermore, state-of-the-art inorganic thin-film solar cells often incorporate a band-gap gradient across their thickness, which can further contribute to the absorptance broadening. Using $\mathrm{Cu}(\mathrm{In},\mathrm{Ga}){\mathrm{Se}}_2$ (CIGSe) absorbers as a case study, we perform a comprehensive analysis of voltage losses through absolute photoluminescence and electroluminescence spectroscopy, combined with photospectrometry and high-spatial-resolution cathodoluminescence measurements. This study has two main purposes: (i) to present a detailed loss analysis based on Planck’s generalized law and (ii) to quantify the radiative losses arising from the band-gap gradient. We find that the loss analysis based on the combination of radiative, generation, and nonradiative losses covers the entire losses and does not leave out any other significant contribution. Samples with a graded band-gap profile show more pronounced broadening of the absorptance onset and up to 16 mV higher radiative losses compared to the samples with uniform band gap. There is also an indication that band-gap-graded samples have larger lateral band-gap inhomogeneity. Needs text.

Accurate performance calibration of multijunction (MJ) solar cells is critical for advancing this technology toward large-scale terrestrial application, yet existing testing methods developed by photovoltaics calibration laboratories remain prohibitively complex and resource intensive for most research laboratories. Current approaches rely on expensive multisource simulators and/or intricate spectral tuning algorithms, limiting accessibility and hindering standardized comparisons of emerging MJ technologies such as perovskite-based tandems. This paper introduces a streamlined spectral irradiance adjustment method for light-emitting diode (LED)-based solar simulators, which are increasingly adopted in the photovoltaics community due to their versatility and cost-effectiveness. The method we present bins LED channels into groups based on the number of junctions of the MJ photovoltaic device under test (DUT) and their corresponding band gaps and incorporates an automated tuning algorithm that eliminates the need to adjust each channel’s irradiance individually. This tuning algorithm requires the relative spectral responsivities of both the DUT and a broadband reference cell, as well as a calibrated spectroradiometer. Measurement validation across various MJ devices, including III-V and perovskite/silicon tandems, demonstrates excellent agreement within 1% with well-established xenon-tungsten multisource simulators and associated tuning algorithms. By enabling precise spectral tuning with readily available equipment and a simple tuning algorithm, our approach bridges the measurement accuracy gap between research laboratories and accredited testing facilities, fostering more reliable device comparisons and accelerating the translation of MJ technologies into real-world energy systems.

Understanding thermal transport in solid electrolytes is essential for improving the performance, reliability, and safety of all-solid-state batteries. Garnet-type lithium-ion conductors are promising candidates for solid electrolytes, yet their thermal-transport mechanisms remain poorly understood. Here, we connect the lattice and ion dynamics of single-crystal garnet-type ${\mathrm{Li}}_{6.5}$${\mathrm{La}}_3$${\mathrm{Zr}}_{1.5}$${\mathrm{Ta}}_{0.5}$${\mathrm{O}}_{12}$ to its intrinsically low thermal conductivity. Our study reveals that the single crystals grown by the floating-zone method exhibit remarkably low glasslike thermal conductivity. Using first-principles calculations and inelastic-neutron-scattering measurements, we identify both the acoustic and numerous optical phonon modes, which stem from the complex crystal structure of the material. Notably, a low-energy optical branch exhibits an avoided crossing with acoustic phonons near 7 meV. These optical modes can enhance the scattering of heat-carrying acoustic phonons and reduce thermal conductivity. Furthermore, the calculated Grüneisen parameters are large, especially for the vibrational modes around 6 meV, indicating strong anharmonicity, with a noticeable contribution from lithium-ion vibrations. A two-channel thermal-transport model is employed to describe the weak temperature dependence of the thermal conductivity, which can be attributed to the substantial contribution of diffuson transport facilitated by the abundance of optical phonons and intrinsic anharmonicity. These results offer valuable insights into the thermal transport in a broad class of ionic conductors of interest for energy conversion and storage applications.

Recent studies have revealed magnetically controllable thermoelectric effects in superconductor/ferromagnet (S/F) structures. A tunable cryogenic thermoelectric generator needs not only a high conversion factor between electricity and heat, but also a large change in the thermoelectric output when switching the magnetic state of the device. However, the reported modifications in thermoelectric power are either minimal, involve superconductors with relatively low critical temperatures (below 1 K), or do not utilize commercially available spintronic materials. Here, we experimentally measure and numerically model thermoelectric effects in fully epitaxial ferromagnet/superconductor/ferromagnet (F/S/F) junctions based on commercially available easily grown materials, as well as their dependence on the magnetic configuration of the ferromagnetic (F) electrodes. We observe sizeable Seebeck coefficients for the parallel alignment of the F electrodes, reaching values of about 100 $\text{μ}\mathrm{V}/$K. Importantly, we find a decrease of the thermoelectric signal of more than an order of magnitude when switching from a parallel to an antiparallel configuration, constituting a large thermoelectric spin-valve effect. Theoretical modeling based on a self-consistent nonequilibrium Keldysh-Usadel Green’s function theory, combined with micromagnetic simulations, qualitatively reproduce the experimental findings. The thermoelectric effect is optimized when there is a large spin-dependent electron-hole asymmetry in the superconductor combined with spin-dependent transmission through the interfaces. These findings pave the way for the development of efficient and versatile cryogenic thermoelectric heat engines.