Efficient Exploitation of Numerical Quadrature with Distance-Dependent Integral Screening in Explicitly Correlated F12 Theory: Linear Scaling Evaluation of the Most Expensive RI-MP2-F12 Term.

We present a linear scaling atomic orbital based algorithm for the computation of the most expensive exchange-type RI-MP2-F12 term by employing numerical quadrature in combination with CABS-RI to avoid six-center-three-electron integrals. Furthermore, a robust distance-dependent integral screening scheme, based on integral partition bounds [Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101], is used to drastically reduce the number of the required three-center-one-electron integrals substantially. The accuracy of our numerical quadrature/CABS-RI approach and the corresponding integral screening is thoroughly assessed for interaction and isomerization energies across a variety of numerical integration grids. Our method outperforms the standard density fitting/CABS-RI approach with errors below 1 µEh even for small grid sizes and moderate screening thresholds. The choice of the grid size and screening threshold allows us to tailor our ansatz to a desired accuracy and computational efficiency. We showcase the approach's effectiveness for the chemically relevant system valinomycin, employing a triple-ζ F12 basis set combination (C54H90N6O18, 5757 AO basis functions, 10,266 CABS basis functions, 735,783 grid points). In this context, our ansatz achieves higher accuracy combined with a 135× speedup compared to the classical density fitting based variant, requiring notably less computation time than the corresponding RI-MP2 calculation. Additionally, we demonstrate near-linear scaling through calculations on linear alkanes. We achieved an 817-fold acceleration for C80H162 and an extrapolated 28,765-fold acceleration for C200H402, resulting in a substantially reduced computational time for the latter─from 229 days to just 11.5 min. Our ansatz may also be adapted to the remaining MP2-F12 terms, which will be the subject of future work.

Polymeric Backbone Eutectogel Electrolytes for High-Energy Lithium-Ion Batteries.

This work introduces a polymeric backbone eutectogel (P-ETG) hybrid solid-state electrolyte with an N-isopropylacrylamide (NIPAM) backbone for high-energy lithium-ion batteries (LIBs). The NIPAM-based P-ETG is (electro)chemically compatible with commercially relevant positive electrode materials such as the nickel-rich layered oxide LiNi0.6Mn0.2Co0.2O2 (NMC622). The chemical compatibility was demonstrated through (physico)chemical characterization methods. The nonexistence (within detection limits) of interfacial reactions between the electrolyte and the positive electrode, the unchanged bulk crystallographic composition, and the absence of transition metal ions leaching from the positive electrode in contact with the electrolyte were demonstrated by Fourier transform infrared spectroscopy, powder X-ray diffraction, and elemental analysis, respectively. Moreover, the NIPAM-based P-ETG demonstrates a wide electrochemical stability window (1.5-5.0 V vs Li+/Li) and a reasonably high ionic conductivity at room temperature (0.82 mS cm-1). The electrochemical compatibility of a high-potential NMC622-containing positive electrode and the P-ETG is further demonstrated in Li|P-ETG|NMC622 cells, which deliver a discharge capacity of 134, 110, and 97 mAh g-1 at C/5, C/2, and 1C, respectively, after 90 cycles. The Coulombic efficiency is >95% at C/5, C/2, and 1C. Hence, gaining scientific insights into the compatibility of the electrolytes with positive electrode materials that are relevant to the commercial market, like NMC622, is important because this requires going beyond the electrolyte design itself, which is essential to their practical applications.

Effects of the energy balance transition on bone mass and strength.

Chronic positive energy balance has surged among societies worldwide due to increasing dietary energy intake and decreasing physical activity, a phenomenon called the energy balance transition. Here, we investigate the effects of this transition on bone mass and strength. We focus on the Indigenous peoples of New Mexico in the United States, a rare case of a group for which data can be compared between individuals living before and after the start of the transition. We show that since the transition began, bone strength in the leg has markedly decreased, even though bone mass has apparently increased. Decreased bone strength, coupled with a high prevalence of obesity, has resulted in many people today having weaker bones that must sustain excessively heavy loads, potentially heightening their risk of a bone fracture. These findings may provide insight into more widespread upward trends in bone fragility and fracture risk among societies undergoing the energy balance transition.

Ab Initio Simulation of the Ultrafast Circular Dichroism Spectrum of Provitamin D Ring-Opening.

We present a method to simulate ultrafast pump-probe time-resolved circular dichroism (TRCD) spectra based on time-dependent density functional theory trajectory surface hopping. The method is applied to simulate the TRCD spectrum along the photoinduced ring-opening of provitamin D. Simulations reveal that the initial decay of the signal is due to excited state relaxation, forming the rotationally flexible previtamin D. We further show that oscillations in the experimental TRCD spectrum arise from isomerizations between previtamin D rotamers with different chirality, which are associated with the helical conformation of the triene unit. We give a detailed description of the formation dynamics of different rotamers, playing a key role in the natural regulation of vitamin D photosynthesis. Going beyond the sole extraction of decay rates, simulations greatly increase the amount of information that can be retrieved from ultrafast TRCD, making it a sensitive tool to unravel details in the subpicosecond dynamics of photoinduced chirality changes.

Ab initio simulation of the ultrafast circular dichroism spectrum of provitamin D ring-opening.

We present a method to simulate ultrafast pump-probe time-resolved circular dichroism (TRCD) spectra based on time-dependent density functional theory trajectory surface hopping. The method is applied to simulate the TRCD spectrum along the photoinduced ring-opening of provitamin D. Simulations reveal that the initial decay of the signal is due to excited state relaxation, forming the rotationally flexible previtamin D. We further show that oscillations in the experimental TRCD spectrum arise from isomerizations between previtamin D rotamers with different chirality, which are associated with the helical conformation of the triene unit. We give a detailed description of the formation dynamics of different rotamers, playing a key role in the natural regulation vitamin D photosynthesis. Going beyond the sole extraction of decay rates, simulations greatly increase the amount of information that can be retrieved from ultrafast TRCD, making it a sensitive tool to unravel details in the sub-picosecond dynamics of photoinduced chirality changes.

High-Temperature Thermal Reactivity and Interface Evolution of the NMC-LATP-Carbon Composite Cathode.

Temperature-assisted densification methods are typically used in oxide-based solid-state batteries to suppress resistive interfaces. However, chemical reactivity among the different cathode components (which include a catholyte, the conducting additive, and the electroactive material) still represents a major challenge and processing parameters need thus to be carefully selected. In this study, we evaluate the impact of temperature and heating atmosphere in the LiNi0.6Mn0.2Co0.2O2 (NMC), Li1+xAlxTi2-xP3O12 (LATP), and Ketjenblack (KB) system. A rationale of the chemical reactions between components is proposed from the combination of bulk and surface techniques and overall involves a cation redistribution in the NMC cathode material that is accompanied by the loss of lithium and oxygen from the lattice enhanced by LATP and KB, which act as lithium and oxygen sinks. The final result is the formation of several degradation products, starting at the surface, that lead to a rapid capacity decay above 400 °C. Both the reaction mechanism and threshold temperature depend on the heating atmosphere, with the air atmosphere being more favorable compared to oxygen or any other inert gases.

Braiding Braak and Braak: Staging patterns and model selection in network neurodegeneration.

A hallmark of Alzheimer's disease is the aggregation of insoluble amyloid-beta plaques and tau protein neurofibrillary tangles. A key histopathological observation is that tau protein aggregates follow a structured progression pattern through the brain. Mathematical network models of prion-like propagation have the ability to capture such patterns, but a number of factors impact the observed staging result, thus introducing questions regarding model selection. Here, we introduce a novel approach, based on braid diagrams, for studying the structured progression of a marker evolving on a network. We apply this approach to a six-stage 'Braak pattern' of tau proteins, in Alzheimer's disease, motivated by a recent observation that seed-competent tau precedes tau aggregation. We show that the different modeling choices, from the model parameters to the connectome resolution, play a significant role in the landscape of observable staging patterns. Our approach provides a systematic way to approach model selection for network propagation of neurodegenerative diseases that ensures both reproducibility and optimal parameter fitting.

Protein-protein interactions in neurodegenerative diseases: A conspiracy theory.

Neurodegenerative diseases such as Alzheimer's or Parkinson's are associated with the prion-like propagation and aggregation of toxic proteins. A long standing hypothesis that amyloid-beta drives Alzheimer's disease has proven the subject of contemporary controversy; leading to new research in both the role of tau protein and its interaction with amyloid-beta. Conversely, recent work in mathematical modeling has demonstrated the relevance of nonlinear reaction-diffusion type equations to capture essential features of the disease. Such approaches have been further simplified, to network-based models, and offer researchers a powerful set of computationally tractable tools with which to investigate neurodegenerative disease dynamics. Here, we propose a novel, coupled network-based model for a two-protein system that includes an enzymatic interaction term alongside a simple model of aggregate transneuronal damage. We apply this theoretical model to test the possible interactions between tau proteins and amyloid-beta and study the resulting coupled behavior between toxic protein clearance and proteopathic phenomenology. Our analysis reveals ways in which amyloid-beta and tau proteins may conspire with each other to enhance the nucleation and propagation of different diseases, thus shedding new light on the importance of protein clearance and protein interaction mechanisms in prion-like models of neurodegenerative disease.

Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing Units.

We present a highly efficient and asymptotically linear-scaling graphic processing unit accelerated seminumerical exact-exchange method (sn-LinK). We go beyond our previous central processing unit-based method (Laqua, H.; Kussmann, J.; Ochsenfeld, C. J. Chem. Theory Comput. 2018, 14, 3451-3458) by employing our recently developed integral bounds (Thompson, T. H.; Ochsenfeld, C. J. Chem. Phys. 2019, 150, 044101) and high-accuracy numerical integration grid (Laqua, H.; Kussmann, J.; Ochsenfeld, C. J. Chem. Phys. 2018, 149, 204111). The accuracy is assessed for several established test sets, providing errors significantly below 1mEh for the smallest grid. Moreover, a comprehensive performance analysis for large molecules between 62 and 1347 atoms is provided, revealing the outstanding performance of our method, in particular, for large basis sets such as the polarized quadruple-zeta level with diffuse functions.

A Schwarz inequality for complex basis function methods in non-Hermitian quantum chemistry.

A generalization of the Schwarz bound employed to reduce the scaling of quantum-chemical calculations is introduced in the context of non-Hermitian methods employing complex-scaled basis functions. Non-Hermitian methods offer a treatment of molecular metastable states in terms of L2-integrable wave functions with complex energies, but until now, an efficient upper bound for the resulting electron-repulsion integrals has been unavailable due to the complications from non-Hermiticity. Our newly formulated bound allows us to inexpensively and rigorously estimate the sparsity in the complex-scaled two-electron integral tensor, providing the basis for efficient integral screening procedures. We have incorporated a screening algorithm based on the new Schwarz bound into the state-of-the-art complex basis function integral code by White, Head-Gordon, and McCurdy [J. Chem. Phys. 142, 054103 (2015)]. The effectiveness of the screening is demonstrated through non-Hermitian Hartree-Fock calculations of the static field ionization of the 2-pyridoxine 2-aminopyridine molecular complex.

Impact of genetic subtypes of Prader-Willi syndrome with growth hormone therapy on intelligence and body mass index.

Prader-Willi syndrome (PWS) is a genomic imprinting disorder characterized by infantile hypotonia with a poor suck and failure to thrive, hypogenitalism/hypogonadism, behavior and cognitive problems, hormone deficiencies, hyperphagia, and obesity. The Stanford Binet and Wechsler (WAIS-R; WISC-III) intelligence (IQ) tests were administered on 103 individuals with PWS from two separate cohorts [University of California, Irvine (UCI) (N = 56) and Vanderbilt University (N = 47)] and clinical information obtained including growth hormone (GH) treatment, PWS molecular classes, weight and height. Significantly higher IQ scores (p < .02) were found representing the vocabulary section of the Stanford Binet test in the growth hormone (GH) treated group when compared with non-GH treatment in the pediatric-based UCI PWS cohort with a trend for stabilization of vocabulary IQ scores with age in the GH treated maternal disomy (UPD) 15 subject group. Significant differences (p = .05) were also found in the adult-based Vanderbilt PWS cohort with 15q11-q13 deletion subjects having lower Verbal IQ scores compared with UPD 15. No difference in body mass index was identified based on the PWS molecular class or genetic subtype. Medical care and response to treatment with growth hormone may influence intelligence impacted by PWS genetic subtypes and possibly age, but more studies are needed.

Integral partition bounds for fast and effective screening of general one-, two-, and many-electron integrals.

We introduce tight upper bounds for a variety of integrals appearing in electronic structure theories. These include electronic interaction integrals involving any number of electrons and various integral kernels such as the ubiquitous electron repulsion integrals and the three- and four-electron integrals found in explicitly correlated methods. Our bounds are also applicable to the one-electron potential integrals that appear in great number in quantum mechanical (QM), mixed quantum and molecular mechanical (QM/MM), and semi-numerical methods. The bounds are based on a partitioning of the integration space into balls centered around electronic distributions and their complements. Such a partitioning leads directly to equations for rigorous extents, which we solve for shell pair distributions containing shells of Gaussian basis functions of arbitrary angular momentum. The extents are the first general rigorous formulation we are aware of, as previous definitions are based on the inverse distance operator 1/r12 and typically only rigorous for simple spherical Gaussians. We test our bounds for six different integral kernels found throughout quantum chemistry, including exponential, Gaussian, and complementary error function based forms. We compare to previously developed estimates on the basis of significant integral counts and their usage in both explicitly correlated second-order Møller-Plesset theory (MP2-F12) and density functional theory calculations employing screened Hartree-Fock exchange.

An implicit discontinuous Galerkin method for modeling acute edema and resuscitation in the small intestine.

Edema, also termed oedema, is a generalized medical condition associated with an abnormal aggregation of fluid in a tissue matrix. In the intestine, excessive edema can lead to serious health complications associated with reduced motility. A $7.5\%$ solution of hypertonic saline (HS) has been hypothesized as an effective means to reduce the effects of edema following surgery or injury. However, detailed clinical edema experiments can be difficult to implement, or costly, in practice. In this manuscript we introduce an implicit in time discontinuous Galerkin method with novel adaptations for modeling edema in the 3D layered physiology of the intestine. The model improves over early work via inclusion of the tissue intrinsic storage coefficient, and the effects of Starling overestimation for high venous pressures. Validation against a recent clinical experiment in HS resuscitation of acute edema is presented; the results support the clinical hypothesis that 7.5% HS solution may be effective in the resuscitation of acute edema formation. New results include an improved view into the effects of resuscitation on the hydrostatic pressure profile of edematous rats, effects on lumenal volume attenuation, relative fluid gain and an estimation of the impacts of both acute edema and resuscitation on intestinal motility.

Tuning the photoreactivity of Z-hexatriene photoswitches by substituents - a non-adiabatic molecular dynamics study.

To understand how substituents can be used to increase the quantum yield of photochemical electrocyclic ring-closing of the Z-hexa-1,3,5-triene (HT) photoswitch forming cyclohexadiene (CHD), we investigate the S1 photo dynamics of HT and its derivatives 2,5-dimethyl-HT (DMHT), 2-isopropyl-5-methyl-HT (2,5-IMHT), 1-isopropyl-4-methyl-HT (1,4-IMHT), and 2,5-diisopropyl-HT (DIHT) using time-dependent density functional theory surface hopping dynamics. We report detailed photoproduct distributions, formation mechanisms, branching ratios, and wavelength-dependent product quantum yields. Most products have been confirmed experimentally and include all-trans HT derivatives, cyclopropanes, cyclobutenes, cyclopentene, cyclohexadienes, and bicyclic compounds. Regarding CHD formation, we find that for the 2,5-substituted derivatives DMHT, 2,5-IMHT, and DIHT, the branching ratios increase with increasing size of the substituents. In contrast the branching ratios of the E/Z-isomerization decrease with increasing size of the substituents. Due to steric interactions, increasing the size of the substituents increases the amount of gZg rotamers in the ground state, which are prone to CHD formation and have lower E/Z-isomerization probability. Furthermore, we find [1,4], [1,5], and [1,6]-sigmatropic hydrogen shift reactions occurring at large percentages (5% to 15%); for sterical reasons these reactions stem from tZg conformers. DIHT shows the lowest percentage of side product formation among the 2,5-substituted molecules and highest CHD branching ratio; its CHD quantum yield can be increased up to more than 64%, by excitation of DIHT on the red tail of its absorption spectrum, whereas the Z/E-isomerization is reduced below 5% and side reactions practically vanish. This makes DIHT the best candidate for applications in molecular switches.

First-Principles Prediction of Wavelength-Dependent Product Quantum Yields.

We present a method to predict wavelength-dependent product quantum yields (PQYs) for photochemical reactions and applied it to Z/E-isomerization and several ring-closing reactions of Z-2,5-dimethyl-1,3,5-hexatriene and truncated previtamin D. Using branching ratios from surface hopping molecular dynamics, individual trajectories are correlated with the absorption spectra of their initial structures. The wavelength-dependent PQYs are computed by dividing the average spectrum of the initial structures of the product-forming trajectories by the average spectrum of all initial structures. Accurate absorption spectra are calculated using the correlated ADC(2) method with an implicit solvent. Calculations reproduce the experimentally found trend of increasing six-ring formation and decreasing Z/E-isomerization on the red side of the spectrum. Over all seven reactions studied, the mean absolute error (MAE) between experimental and calculated PQYs (MAE) amounts to 8.1%, and the largest MAE is 18.6%. For four reactions, predicted values agree quantitatively with experiments within 5.6%.

Distance-including rigorous upper bounds and tight estimates for two-electron integrals over long- and short-range operators.

We introduce both rigorous and non-rigorous distance-dependent integral estimates for four-center two-electron integrals derived from a distance-including Schwarz-type inequality. The estimates are even easier to implement than our so far most efficient distance-dependent estimates [S. A. Maurer et al., J. Chem. Phys. 136, 144107 (2012)] and, in addition, do not require well-separated charge-distributions. They are also applicable to a wide range of two-electron operators such as those found in explicitly correlated theories and in short-range hybrid density functionals. For two such operators with exponential distance decay [e-r12 and erfc(0.11⋅r12)/r12], the rigorous bound is shown to be much tighter than the standard Schwarz estimate with virtually no error penalty. The non-rigorous estimate gives results very close to an exact screening for these operators and for the long-range 1/r12 operator, with errors that are completely controllable through the integral screening threshold. In addition, we present an alternative form of our non-rigorous bound that is particularly well-suited for improving the PreLinK method [J. Kussmann and C. Ochsenfeld, J. Chem. Phys. 138, 134114 (2013)] in the context of short-range exchange calculations.

The role of tachysterol in vitamin D photosynthesis - a non-adiabatic molecular dynamics study.

To investigate the role of tachysterol in the photophysical/photochemical regulation of vitamin D photosynthesis, we studied its electronic absorption properties and excited state dynamics using time-dependent density functional theory (TDDFT), second-order approximate coupled cluster theory (CC2), and non-adiabatic surface hopping molecular dynamics in the gas phase. In excellent agreement with experiments, the simulated electronic spectrum shows a broad absorption band with a remarkably higher extinction coefficient than the other vitamin D photoisomers provitamin D, lumisterol, and previtamin D. The broad band arises from the spectral overlap of four different ground state rotamers. After photoexcitation, the first excited singlet state (S1) decays with a lifetime of 882 fs. The S1 dynamics is characterized by a strong twisting of the central double bond. In 96% of all trajectories this is followed by unreactive relaxation to the ground state near a conical intersection. The double-bond twisting in the chemically unreactive trajectories induces a strong interconversion between the different rotamers. In 2.3% of the trajectories we observed [1,5]-sigmatropic hydrogen shift forming the partly deconjugated toxisterol D1. 1.4% previtamin D formation is observed via hula-twist double bond isomerization. In both reaction channels, we find a strong dependence between photoreactivity and dihedral angle conformation: hydrogen shift only occurs in cEc and cEt rotamers and double bond isomerization occurs mainly in cEc rotamers. Hence, our study confirms the previously formed hypothesis that cEc rotamers are more prone to previtamin D formation than other isomers. In addition, we also observe the formation of a cyclobutene-toxisterol in the hot ground state in 3 trajectories (0.7%). Due to its large extinction coefficient and mostly unreactive behavior, tachysterol acts mainly as a Sun shield suppressing previtamin D formation. Tachysterol shows stronger toxisterol formation than previtamin D and can thus be seen as the major degradation route of vitamin D. Absorption of low energy ultraviolet light by the cEc rotamer can lead to previtamin D formation. In addition, the cyclobutene-toxisterol, which possibly reacts thermally to previtamin D, is also preferably formed at long wavelengths. These two mechanisms are consistent with the wavelength dependent photochemistry found in experiments. Our study reinforces a recent hypothesis that tachysterol constitutes a source of previtamin D when only low energy ultraviolet light is available, as it is the case in winter or in the morning and evening hours of the day.

Cyclohexadiene Revisited: A Time-Resolved Photoelectron Spectroscopy and ab Initio Study.

We have reinvestigated the excited state dynamics of cyclohexa-1,3-diene (CHD) with time-resolved photoelectron spectroscopy and fewest switches surface hopping molecular dynamics based on linear response time-dependent density functional theory after excitation to the lowest lying ππ* (1B) state. The combination of both theory and experiment revealed several new results: First, the dynamics progress on one single excited state surface. After an incubation time of 35 ± 10 fs on the excited state, the dynamics proceed to the ground state in an additional 60 ± 10 fs, either via a conrotatory ring-opening to hexatriene or back to the CHD ground state. Moreover, ring-opening predominantly occurs when the wavepacket crosses the region of strong nonadiabatic coupling with a positive velocity in the bond alternation coordinate. After 100 fs, trajectories remaining in the excited state must return to the CHD ground state. This extra time delay induces a revival of the photoelectron signal and is an experimental confirmation of the previously formulated model of two parallel reaction channels with distinct time constants. Finally, our simulations suggest that after the initially formed cis-Z-cis HT rotamer the trans-Z-trans isomer is formed, before the thermodynamical equilibrium of three possible rotamers is reached after 1 ps.

Analysis of human rotaviruses from a single location over an 18-year time span suggests that protein coadaption influences gene constellations.

UNLABELLED: Rotaviruses (RVs) are 11-segmented, double-stranded RNA viruses that cause severe gastroenteritis in children. In addition to an error-prone genome replication mechanism, RVs can increase their genetic diversity by reassorting genes during host coinfection. Such exchanges allow RVs to acquire advantageous genes and adapt in the face of selective pressures. However, reassortment may also impose fitness costs if it unlinks genes/proteins that have accumulated compensatory, coadaptive mutations and that operate best when kept together. To better understand human RV evolutionary dynamics, we analyzed the genome sequences of 135 strains (genotype G1/G3/G4-P[8]-I1-C1-R1-A1-N1-T1-E1-H1) that were collected at a single location in Washington, DC, during the years 1974 to 1991. Intragenotypic phylogenetic trees were constructed for each viral gene using the nucleotide sequences, thereby defining novel allele level gene constellations (GCs) and illuminating putative reassortment events. The results showed that RVs with distinct GCs cocirculated during the vast majority of the collection years and that some of these GCs persisted in the community unchanged by reassortment. To investigate the influence of protein coadaptation on GC maintenance, we performed a mutual information-based analysis of the concatenated amino acid sequences and identified an extensive covariance network. Unexpectedly, amino acid covariation was highest between VP4 and VP2, which are structural components of the RV virion that are not thought to directly interact. These results suggest that GCs may be influenced by the selective constraints placed on functionally coadapted, albeit noninteracting, viral proteins. This work raises important questions about mutation-reassortment interplay and its impact on human RV evolution. IMPORTANCE: Rotaviruses are devastating human pathogens that cause severe diarrhea and kill >450,000 children each year. The virus can evolve by accumulating mutations and by acquiring new genes from other strains via a process called reassortment. However, little is known about the relationship between mutation accumulation and gene reassortment for rotaviruses and how it impacts viral evolution. In this study, we analyzed the genome sequences of human strains found in clinical fecal specimens that were collected at a single hospital over an 18-year time span. We found that many rotaviruses did not reassort their genes but instead maintained them as specific sets (i.e., constellations). By analyzing the encoded proteins, we discovered concurrent amino acid changes among them, which suggests that they are functionally coadapted to operate best when kept together. This study increases our understanding of how rotaviruses evolve over time in the human population.

Conversion efficiency of skutterudite-based thermoelectric modules.

Presently, the only commercially available power generating thermoelectric (TE) modules are based on bismuth telluride (Bi2Te3) alloys and are limited to a hot side temperature of 250 °C due to the melting point of the solder interconnects and/or generally poor power generation performance above this point. For the purposes of demonstrating a TE generator or TEG with higher temperature capability, we selected skutterudite based materials to carry forward with module fabrication because these materials have adequate TE performance and are mechanically robust. We have previously reported the electrical power output for a 32 couple skutterudite TE module, a module that is type identical to ones used in a high temperature capable TEG prototype. The purpose of this previous work was to establish the expected power output of the modules as a function of varying hot and cold side temperatures. Recent upgrades to the TE module measurement system built at the Fraunhofer Institute for Physical Measurement Techniques allow for the assessment of not only the power output, as previously described, but also the thermal to electrical energy conversion efficiency. Here we report the power output and conversion efficiency of a 32 couple, high temperature skutterudite module at varying applied loading pressures and with different interface materials between the module and the heat source and sink of the test system. We demonstrate a 7% conversion efficiency at the module level when a temperature difference of 460 °C is established. Extrapolated values indicate that 7.5% is achievable when proper thermal interfaces and loading pressures are used.