Comparing first-principles methods for computing electrical conductivity in warm dense matter
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.
Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.
Silicon carbide has been identified as a useful material for nuclear application due to its heat resistance, low neutron absorbing cross-section, chemical inertness, and its rigidity as a structural material. These properties make it very favorable for high temperature environments such as in TRISO fuel and advanced reactor designs. Treatments such as n-doping, where carrier atoms are added to increase the number of free electrons, are expected to further improve properties such as thermal conductivity and electrical resistivity. Due to the high amount of neutron interactions in a nuclear reactor, it is expected that silicon carbide used in a reactor will be passively n-doped via transmutation. Studying the effects of n-doping in silicon carbide can help us better understand how the material will operate under real world conditions. We have found that n-doping the SiC increases both thermal and electrical conductivity, providing positive feedback under reactor usage. Previously it had been thought that under reactor conditions material properties of SiC could only degrade due to defects caused by neutron interactions. This study shows that transmutation doping can help to counteract and slow this process. This research shows that SiC can be used in nuclear reactor components contrasting with more costly and complicated alternative materials.
Effective lightning strike protection for critical aerospace and wind applications requires high electrical conductivity to dissipate current efficiently. However, polymer matrix composites face a challenge due to their inherently insulating nature. While conventional carbon fiber-reinforced composites (CFRP) exhibit electrical conductivity in the planar direction, achieving through-thickness conductivity remains an ongoing challenge. In this work, we have undertaken the fabrication of CFRP interleaved with vertically oriented carbon fibers (Z-fiber) to impart higher electrical conductivity along the thickness direction. Two Z-fiber composite variations are prepared: Z-1 with a single layer of Z-fiber and Z-5 with five interleaved layers and compared with no Z-fiber layer (Z-0) composite. The composite panels were subjected to lab-scale lightning strike tests with a current magnitude of 100 kA. To emulate real-world service conditions, an aerospace-grade paint coating was applied to the composite laminates. Comparative analysis shows Z-1 reduces damage diameter to ∼22 mm compared to Z-0 (∼26 mm), while Z-5 exhibits the least damage (∼16.7 mm), confirmed by optical microscopy. Z-5 demonstrates nine times higher through-thickness electrical conductivity than Z-0, reducing electrical anisotropy substantially. Thermal-electric finite element damage modeling predicts surface damage within 6% of experimental values for both Z-0 and Z-5 composites. Flexural tests post-lightning reveal Z-5 retains 66% flexural strength and 86% modulus, significantly better than Z-0, which retains less than 40% for both properties. This study highlights the efficacy of Z-fiber composites in lightning strike protection, offering improved through-thickness conductivity and mechanical property retention.
There is a current need for new aluminum alloy design strategies to target applications requiring high strength and conductivity with reductions in mass. A new lightweight Al-2Ni-0.5Zr (wt. %) conductor alloy was fabricated using laser powder bed fusion. A design of experiments probed the alloy's solidification cracking susceptibility. It was observed that solidification cracking was generally reduced with fast scan speeds, above 1500 mm/s, and smaller hatch spacings. The different cooling rates throughout the melt pool produced a heterogeneous distribution of cellular and equiaxed Al 3 Ni precipitates in the as-printed alloy. Additionally, the rapid solidification characteristic of laser powder bed fusion created a super-saturated Zr solid solution. An aging heat treatment at 375 °C for 24 h imparted strengthening through the precipitation of L1 2 -Al 3 Zr nanoprecipitates, which counteracted the softening caused by the fragmentation and coarsening of Al 3 Ni precipitates. The yield strength increased from 138 MPa in the as-printed condition to 168 MPa after aging, while the ductility remained constant at ∼21%. The aging treatment simultaneously increased the electrical conductivity from 40.8% IACS (International Annealed Copper Standard) to 53.5% IACS. Modeling of the strengthening mechanisms and electrical conductivity contributions rationalized the simultaneous increase in strength and conductivity upon aging. Furthermore, the strengthening efficacy of the Al 3 Ni and L1 2 -Al 3 Zr precipitates, combined with the low Ni and Zr solubility in the FCC Al matrix, facilitated both high strength and electrical conductivity. Overall, the combination of strength and electrical conductivity positions this alloy as a suitable choice for additively manufactured lightweight conductors.
Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500 C to 1500 C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 D/m^2-s. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.
Reliable material performance is required for plasma-facing material (PFM) candidates. Previous research has shown that plasma and neutron radiation exposure induces microstructural changes in PFMs; changes in thermal and electrical conductivities and in material hardening and embrittlement were also observed after neutron irradiation. These material property changes will negatively impact the performance of the PFMs in a fusion reactor. Despite the well-known connection between material microstructure, properties, and performance, there is a need for validated modeling capabilities connecting PFM property degradation with microstructural evolution under fusion-relevant conditions. We are developing a simulation capability to couple plasma-induced microstructural evolution to material property degradation. Our approach relies on deliberate mapping between individual simulation models and experimental characterization for validation. The open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) software was used for this simulation capability development. A MOOSE phase-field model was coupled with the cluster dynamics code, Xolotl, to predict microstructural evolution. Microstructure characterization techniques, including scanning electron microscopy (SEM), transmission electron microscopy (TEM), and laser scanning confocal microscopy (LSCM) are used to validate these microstructural evolution simulations. Calculation of thermal and electrical conductivities with first principles simulations was performed for bulk material and for grain boundaries; these results are used within MOOSE models to calculate effective thermal and electrical conductivities as a function of grain characteristics. Thermoreflectance and four-probe techniques were employed to measure the thermal and electrical conductivities, respectively. A MOOSE crystal plasticity model was adapted to predict microstructure-sensitive deformation behavior, and X-ray diffraction (XRD) was used to collect bulk dislocation density data for validation. After individual simulation validation, these models are coupled to predict material property changes resulting from plasma exposure. We focused here on an experimental design to emphasize the separate effects of moderate thermal loads and plasma exposure using tungsten. Annealing of tungsten was performed under a protective environment for temperatures ranging from 500$^o$C to 1500$^o$C. The plasma exposure was completed in the Tritium Plasma Experiment at Idaho National Laboratory under a deuterium flux of 1e22 $\frac{D}{m^2s}$. This incremental approach is employed to build confidence in the modeling capability: separate-effects tests ensure that the models capture key mechanisms from single environmental conditions before predicting PFM property degradation under combined loads. We will show our early results from coupling these simulation models to predict PFM property changes from microstructural evolution. Comparisons of the simulation results with preliminary validation data will be discussed.
Electrically conductive composites are emerging across diverse industries such as electronic, automotive, aerospace, advanced air mobility, biomedical, infrastructure, defense and security offering static charge dissipation, electromagnetic interference shielding, lighting protection, sensing, dicing, corrosion monitoring, etc. Conductivity enhanced composites provide several advantages compared to conventional metals including weight reduction, corrosion resistance, energy efficient processability, tunable properties and multifunctionality. Polymers are typically insulating in nature and require conducting filler for electron transport. However, dispersion and polymer-filler interphases are critical and often disrupt conducting pathways. Besides, conductive fillers such as graphene, carbon nanotube, MXene, silver nanowire, etc. are expensive, limiting their wide adoption in composite industry. On the other hand, a limited number of intrinsically conductive polymers are available among which polyaniline (PANI) has been widely studied due to its high conductivity, thermal and chemical stability. However, PANI is difficult to process and exhibits weak mechanical properties. In brief, there is a significant demand for electrically conductive polymer formulation with cost-effective conducting fillers that offer processability in scale to expand the market of conductivity enhanced materials.
SiGe alloys have a proven track record as robust high-temperature thermoelectric materials, powering NASA missions like SNAP-10A, LES-9, and Voyager 1 and 2. Enhancing thermoelectric efficiency hinges on minimizing thermal conductivity while preserving electrical conductivity. Rapid solidification via LPBF can generate microstructural features such as subgrain cellular boundaries and twinning, which may help reduce thermal conductivity while preserving semiconducting behavior. Here, this study investigates the potential of laser powder bed fusion (LPBF) additive manufacturing to fabricate nanostructured Si and SiGe thermoelectric materials. High cooling rates (10 5 to 10 7 K/s) rates inherent to the LPBF process are conducive to forming such nanostructures. Moreover, this fabrication technique could also be suitable for fabricating complex geometries needed to achieve improved device level performance. Process mapping of commercial Si powder with irregular morphology was first performed to understand the LPBF processing behavior of this semiconductor material. Subsequent studies included in-house synthesized B-doped (p-type) Si7 8 Ge 22 spherical powder that was produced via ultrasonic atomization. Scan strategies involved multiple laser exposures to mitigate solidification cracking: a high density of >98% was achieved, but solidification cracking could not be fully eliminated. Subsequently, a high electrical resistivity (i.e., low conductivity) was observed, but the measured Seebeck coefficient, ∼230 μV/K @ 500 °C, proved that good semiconductor material was being fabricated. A subgrain cellular structure (5–10 μm) was observed as defined by Ge segregation to the intercellular boundaries. The remelting strategies helped lower the cooling rates in processing SiGe, but this still resulted in high residual stresses, which induced a remarkably high density of twins (78–95%) to accommodate the deformation. This unique grain structure offers an avenue for phonon scattering and potential improvements in thermoelectric performance.
Thermoelectric materials with high electrical conductivity and low thermal conductivity (e.g., Bi 2 Te 3 ) can efficiently convert waste heat into electricity. However, despite favorable theoretical predictions, individual Bi 2 Te 3 nanostructures such as two-dimensional (2D) nanoplates tend to underperform bulk Bi 2 Te 3 . We report a novel surface doping technique to synthesize highly n-type Bi 2 Te 3 nanoplates using an external Cr coating followed by a thermal annealing process in a reducing atmosphere, as well as the mechanism by which this surface coating – only a few atoms or less in thickness – can observably impact the thermoelectric performance of 2D Bi 2 Te 3 . The Cr atoms act as n-type carrier donors by directly incorporating into the Bi 2 Te 3 structure during thermal annealing, enhancing electrical conductivity by ∼ 70 % while increasing thermal conductivity by only ∼ 5 % at room temperature. Compared to the uncoated Bi 2 Te 3 nanoplate, the Cr-doped Bi 2 Te 3 nanoplate exhibits a doubled thermoelectric figure of merit (zT), which is still relatively low. Raman spectroscopy and chemical potential simulations further confirm that Cr atoms are incorporated into the Bi 2 Te 3 structure.
Thermal conductivity in nanoparticle solids has been previously reported in the range of 0.1-1 W m -1 K -1 , which is a much smaller variation than the orders of magnitude differences achievable in electrical conductivity of similar systems. Both the low absolute magnitude of thermal conductivity and the relative insensitivity compared to electrical conductivity may be largely attributed to the poor interfacial thermal conductance of the many interfaces of the nanocrystal solid, but a direct experimental study of these interfaces is challenging. Here, in this work, we overcome this challenge via direct spectroscopic observation of heat flow within the components of a nanocrystal solid. These thermal transfers are studied by mixing two distinct types of particles: one that serves as a selectively excited antenna to inject heat and the other as the thermal acceptor to report the time-dependent change in temperature. Using transient spectroscopy, the equilibration between the heat donor and heat acceptor is observed to require similar to ~300 ps at room temperature, speeds up at reduced temperature, and has only weak sensitivity to the relative stoichiometry of the components or the intervening ligands. These results contrast strongly with the 10-20 ps time-scale of through-bond heat transfer at the ligand-particle surface and highlight the substantially lower interfacial thermal conductance of particle-to-particle transport without covalent bonding. It is also found, serendipitously, that the mixed composite films show an unexpected, substantially enhanced nonlinear absorption at the resonant wavelength of the plasmonic particles, which is tentatively attributed to a local field enhancement.
Copper catalyst gas diffusion electrodes (GDEs) have demonstrated unique electrochemical selectivity converting CO 2 to C 2 -hydrocarbons such as ethylene and ethanol but have been challenged by their hydrophobic chemical stability and internal electrical resistance leading to low energy efficiency. Carbon-supported GDEs have low electrical resistance but lack sufficient stability at industrially relevant current densities. While polymer-supported GDEs have improved hydrophobicity, they also display high in-plane electrical resistance, particularly at industrial scales. Here, in this work, we demonstrate a composite gas diffusion layer that combines hydrophobic porous polymers with an electrically conductive backbone addressing these core gas diffusion electrode (GDE) scaling challenges. We investigate the material properties of standalone porous perfluoropolyether (PFPE) polymers, including porosity and surface morphology, under varying processing conditions and then incorporate these polymers into a porous copper foam. This composite enhances the mechanical rigidity necessary for cell assembly and provides a through-plane electrical conduction path to reduce electrical resistive losses. This enhanced PFPE composite GDE displays efficient CO 2 reduction, achieving 15% ethylene energy efficiency at 100 cm 2 . These findings contribute to the development of advanced catalyst materials and electrode architectures and promote scalable strategies for electrochemical conversion of CO 2 into high-value carbon products.
This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.
This dataset contains Electromagnetic Induction (EMI) data collected at Trail Creek, Colorado, in 2024. EMI surveys were conducted to investigate the spatial distribution of electrical conductivity in the subsurface, providing insights into soil moisture and subsurface geological features. The surveys were performed along multiple transects to capture variations in conductivity influenced by changes in soil composition, moisture content, and underlying geological structures. This dataset complements other geophysical data collected in the region, including Electrical Resistivity Tomography (ERT) and Terrestrial LiDAR Scanning (TLS), providing a detailed understanding of the subsurface and its impact on surface vegetation and hydrological processes. The data are valuable for environmental geophysics, ecological research, and hydrological modeling in mountainous ecosystems. The files include: - data.zip: the raw EMI data (.csv) - inversion.zip: the inverted resistivity model (.csv and .kml) - kriging.zip: the kriging resistivity model (.csv, .tif, .kmz) - flmd.csv: file level metadata file describing all files within this dataset - dd.csv: data dictionary file describing the column headers within CSV files This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.
Since the WTP compositional space of HLW glasses is extremely large, development of HLW property-composition models is a multi-year task consisting of multiple phases. The primary focus of this work was to enhance the WTP HLW models of interest including PCT releases, spinel crystallization (T1%), viscosity, electrical conductivity, and TCLP-Cd response. Model development to predict nepheline formation upon CCC is the subject of a separate task. In particular, the earlier work has produced property-composition models for glass melt viscosity and glass melt electrical conductivity that showed good performance, while models for PCT releases and spinel crystallization (T1%) required improvement. Therefore, more efforts were directed in the present work to collect data and improve model performance for HLW glass PCT releases and spinel crystallization than for other properties. The present work is a continuation of earlier development phases and is responsive to the applicable Test Plan.