Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion current”

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

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

At least 397 records · Page 22

Fluorinated Interphase Enables Reversible Zn 2+ Storage in Aqueous ZnSO 4 Electrolytes

Aqueous zinc-ion batteries (AZIBs) using ZnSO 4 aqueous electrolytes complement Li-ion batteries and offer high safety, low cost, and scalability. However, an inferior cycle life, attributed to the generation of basic layered double hydroxides (LDHs) on the cathode during cycling, hampers the development of AZIBs. Although Zn(CF 3 SO 3 ) 2 (Zn(OTf) 2 ) electrolytes demonstrate exceptional performance, their higher cost in comparison to ZnSO 4 offsets their electrochemical stability benefits. Toward this end, we propose a fluorinated interphase strategy to achieve a stable battery with ZnSO 4 electrolytes by in situ pre-constructing the cathode-electrolyte interphase (CEI) of ZnOTf-LDH on the cathode surface. Unlike ZnSO 4 - LDH, which obstructs the diffusion channel of zinc ions, ZnOTf-LDH facilitates the desolvation of zinc ions, thereby enhancing the cycling stability (over 50 cycles at a low current density of 200 mA g -1 ). Further, this work offers valuable insights into the degradation mechanism of AZIBs and provides an effective approach for developing stable and low-cost AZIBs.

25 ENERGY STORAGE↗

Development of in situ thickness sensor for vapor diffused Nb3Sn films

Fermilab is one of the leaders in development of vapor diffused Nb3Sn films inside niobium cavities. This material has a higher critical temperature (Tc) than niobium, enabling cavity operation at 4.2 K. This higher operational temperature significantly reduces the infrastructure required for cooling compared to 2 K systems, making superconducting radio-frequency (SRF) technology more accessible. Current deposition methods have relied on iterative testing to determine nominal film thickness, a process that can be time-consuming and imprecise. To address this, we are developing a sensor to measure the thickness of Nb3Sn thin film in situ during vapor diffusion. Our design involves the resistance measurement of a thin film of niobium. During coating, the change in resistance reflects the conversion of the film to Nb3Sn, which allows simple integration with the current furnace infrastructure. This sensor would allow real time measurement of the thickness, allowing for increased precision in future depositions.

Howard, K. [Unlisted] (ORCID:0000000325613308)↗

Electrochemical mechanisms of leakage‐current in photovoltaic modules

Abstract This paper analyzes the mechanisms and pathways for leakage current flow observed in Si photovoltaic modules subjected to high temperature and humidity and a large voltage bias with respect to ground. The current inside of the frame is in the form of electron motion, but in the glass and polymer, it is in at least a large part attributable to the movement of ions. When the mode of current flow changes from electronic to ionic conduction, electrochemical reactions will take place at the interface. This can include reactions that produce volatile chemical species like H 2 , CO x , and O 2 , along with ionic species such as OH − and H 3 O + . Here, we see evidence of the importance of different charge carriers with different diffusion rates and the influence of electrochemical processes involved. The application of negative voltage to the cell circuit affects the resistivity of glass producing surfaces with poor conductivity but with some increases in the electrochemical potential producing complicated interactions that are important when the voltage is changed. In the polymer, there is the development of a space charge region and a chemical gradient providing two oppositional forces to current flow, which when released create a complicated discharge process. Here, we give a basic understanding of the charge/discharge of PV cells highlighting how the specific mechanisms are important in understanding some of the degradation processes in PV modules. We find that there is evidence of multiple significant charge carrier species with different diffusion time scales. The glass/polymer interface forms a depleted region of higher resistance after prolonged exposure to current. Charge also builds up at the polymer to cell antireflective coating interface and mostly flows to the gridlines to experience electrochemical reactions. These complexities result in non‐linear behavior where the apparent resistivities of the different layers change during charge/discharging processes, making the modeling of the current flow extremely difficult.

14 SOLAR ENERGY↗

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗

The nanoscale distribution of copper and its influence on charge collection in CdTe solar cells

For decades, Cu has been the primary dopant in CdTe photovoltaic absorbers. Typically, Cu acceptor concentrations in these devices are on the order of 1 1014 cm-3, which has made it notoriously difficult to directly correlate nanoscale Cu distributions to the local charge transport properties of these devices. To measure and correlate these properties, measurement techniques require high sensitivity to elemental concentration, large penetration depth, and operando compatibility. Techniques such as secondary-ion mass spectroscopy and X-ray energy dispersive spectroscopy are widely adopted to measure Cu concentrations, but they are limited by penetration depth, sensitivity, or spatial resolution. Additionally, they lack the operando capabilities required to correlate one-to-one Cu concentrations to electrical performance. In this work, correlative X-ray microscopy is used to investigate the spatial distribution of Cu and its impact on charge collection through the depth and breadth of CdTe photovoltaic devices. Plan-view, nanoscale X-ray fluorescence maps clearly demonstrate the spatial segregation of copper around regions thought to be CdTe grain boundaries. Complementary cross-section imaging unveils the transition of the maximum charge-collection efficiency from the ZnTe-CdTe interface to the CdS-CdTe interface as a function of Cu incorporation. The copper concentration through the depth of the CdTe layer is characterized by slow and fast diffusion components, and cross-section charge-transport modeling shows that the experimentally obtained charge collection can be explained by the modeled acceptor distribution through the depth of the CdTe layer.

14 SOLAR ENERGY↗

Liquid–Solid Boundaries Dominate Activity of CO 2 Reduction on Gas-Diffusion Electrodes

Electrochemical CO 2 electrolysis to produce hydrocarbon fuels or material feedstocks offers a renewable alternative to fossilized carbon sources. Gas-diffusion electrodes (GDEs), composed of solid electrocatalysts on porous supports positioned near the interface of a conducting electrolyte and CO 2 gas, have been able to demonstrate the substantial current densities needed for future commercialization. These higher reaction rates have often been ascribed to the presence of a three-phase interface, where solid, liquid, and gas provide electrons, water, and CO 2 , respectively. Conversely, mechanistic work on electrochemical reactions implicates a fully two-phase reaction interface, where gas molecules reach the electrocatalyst’s surface by dissolution and diffusion through the electrolyte. Because the discrepancy between an atomistic three-phase versus two-phase reaction has substantial implications for the design of catalysts, gas-diffusion layers, and cell architectures, the nuances of nomenclatures and governing phenomena surrounding the three-phase-region require clarification. In this paper we outline the macro, micro, and atomistic phenomena occurring within a gas-diffusion electrode to provide a focused discussion on the architecture of the often-discussed three-phase region for CO 2 electrolysis. From this information, we comment on the outlook for the broader CO 2 electroreduction GDE cell architecture.

CO2 reduction↗

Application of electron backscatter diffraction techniques to quantify effects of aging on sub-grain and spatial heterogeneity in NMC cathodes

Identification and evaluation of structural heterogeneity in spent cathode materials is crucial to developing appropriate remediation strategies for novel recycling processes. Native heterogeneities may be exacerbated during the cell's operational lifetime, as sub-particle-scale variations induce anisotropic expansion and contraction upon cycling. Structural transformations resulting from repeated cycling and calendar aging predominantly occur at the secondary particle surface and at the grain boundaries (GBs) between primary particles. However, the diffusion and stress build up around and across GBs are poorly understood. In this study, electron backscatter diffraction (EBSD) is employed to track sub-grain lattice structure across a statistically relevant number of Li(Ni0.33Mn0.33Co0.33)O2 (NMC-111) particles. Specifically, differences in lattice misorientation – measured as the deviation from the grain's average orientation – are tracked as a function of position within the electrode (near current collector, middle, and near separator) and as a function of electrochemical cycling. Further, a novel method of structural analysis is developed, offering insight into sub-grain diffusion behavior by comparing lattice misorientation near the grain boundary versus in the grain bulk. The present results suggest that electrode-scale spatial heterogeneity in lattice structure is induced by initial manufacturing conditions, and that radial gradients in lattice misorientation evolve at the primary particle scale with repeated electrochemical cycling.

25 ENERGY STORAGE↗

Stability of Cu(In x Ga 1− x )Se 2 Solar Cells Utilizing RbF Postdeposition Treatment under a Sulfur Atmosphere

Alkali halide postdeposition treatments (PDTs) have become a key tool to maximize efficiency in Cu(In x Ga 1− x )Se 2 (CIGS) photovoltaics. RbF PDTs have emerged as an alternative to the more common Na‐ and K‐based techniques. This study utilizes temperature‐dependent current–voltage ( JVT ) measurements to study a unique RbF PDT performed in a S atmosphere. The samples are measured before and after 6 months in a desiccator to study device stability. Both samples contain Na and K which diffuse from the soda–lime glass substrate. A reference sample and a RbF + S PDT sample both show the development of a rear contact barrier after aging. The contact barrier is higher for the RbF + S PDT sample, leading to decreased current in forward bias. Series resistance is also higher in the RbF + S PDT device which leads to lower fill factor. However, after aging the reference sample has a larger decrease in open‐circuit voltage ( V OC ). Ideality factor measurements suggest Shockley–Read–Hall recombination dominates both samples. V OC versus temperature and a temperature‐dependent activation energy model are used to calculate diode activation energies for each sample condition. Both techniques produce similar values that indicate recombination primarily occurs within the bulk absorber.

14 SOLAR ENERGY↗

Development of in situ thickness sensor for vapor diffused Nb3Sn films

Fermilab is one of the leaders in development of vapor diffused Nb3Sn films inside niobium cavities. This material has a higher critical temperature (Tc) than niobium, enabling cavity operation at 4.2 K. This higher operational temperature significantly reduces the infrastructure required for cooling compared to 2 K systems, making superconducting radio-frequency (SRF) technology more accessible. Current deposition methods have relied on iterative testing to determine nominal film thickness, a process that can be time-consuming and imprecise. To address this, we are developing a sensor to measure the thickness of Nb3Sn thin film in situ during vapor diffusion. Our design involves a four-point resistance measurement of a thin film of niobium, inside the coating region. During coating, the change in resistance reflects the conversion of the film from Nb to Nb3Sn, which allows simple integration with the current furnace infrastructure. This sensor would allow real time measurement of the Nb3Sn film thickness, allowing for increased precision in future depositions for cavity applications.

Howard, K. [Chicago U.] (ORCID:0000000325613308)↗

Through‐Plane Conductive Hydrophobic Electrodes for CO 2 Electrolysis to Ethylene

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.

Chemistry↗

Bayesian optimization of dynamic Ce( III ) electrooxidation

Abstract Cerium(III/IV) is a redox couple with promising applications in redox flow batteries and redox‐mediated water electrolyzers, but due to its high oxidation potential, parasitic oxygen evolution can occur and decrease the Faradaic efficiency (FE) of the oxidation reaction. Dynamic potential dosing can help balance the transport limitations and the electrooxidation rates, allowing for control over the cerium(III) concentration near the electrode, enhancing FE. Due to the large number of possible pulsing operating conditions, Bayesian optimization (BO) was used to rapidly identify the optimal potential pulses. With an active pulse voltage of 2.5 V vs. Ag/AgCl, the maximum FE achieved was 0.91 (active pulse time = 5 ms, resting pulse time = 136 ms), underscoring the need for long inactive times for cerium(III) diffusion to the electrode. Furthermore, a multiobjective BO was performed to identify optimal trade‐offs between FE and the partial current density to achieve improvements over constant current operation.

Engineering↗

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On Rate-Limiting Mechanisms in NMC Cathodes: The Interplay of Low and High Current Constraints

In this study, we investigate the rate performance limits in LiNi 0.6 Mn 0.2 Co 0.2 O 2 cathodes for lithium-ion batteries, focusing on how cathode thickness, porosity, and threshold voltage impact discharge capacity. By conducting galvanostatic discharge experiments across a wide range of current densities using cathodes of varying thicknesses and porosities, we identify two distinct rate-limiting mechanisms: ionic liquid-diffusion and Ohmic/charge-transfer limitations. Our findings show that thick, dense cathodes are primarily limited by lithium diffusion in the electrolyte, while thinner and more porous cathodes are dominated by Ohmic/charge-transfer limitations, particularly at higher lower cutoff voltages. Crucially, these results allow us to identify the rate-limiting mechanisms in different cathode configurations, offering clear insights into how cathode design can be optimized for improved performance. Understanding these mechanisms is essential for designing next-generation batteries with enhanced rate performance, which is critical for applications such as electric vehicles and renewable energy storage systems.

Brischetto, Martin (ORCID:0000000339865813)↗

Decoupling Accurate Electrochemical Behaviors for High-Capacity Electrodes via Reviving Three-Electrode Vehicles

Developing high-capacity electrodes requires the evaluation of electrochemical behaviors with an increasing current density. Currently, the current density for evaluation of high-capacity electrodes has reached a new stage where the polarization at the lithium counter electrode has become a technical barrier for the accurate evaluation of battery electrodes, resulting in severe performance and mechanism mischaracterizations. Here, the accurate electrochemical behavior for high-capacity electrodes via a single-channel three-electrode vehicle is decoupled, by which the impact of lithium counter electrode is minimized. The testing high-capacity graphite electrode is capable of delivering an excellent rate capability with 81.7% capacity retention at 0.3 C, as well as stable cycling performance retaining 97.5% practical reversible capacity after 225 cycles, much higher than the graphite electrode tested with traditional half-cell testing vehicle but in close agreement with the results obtained from a well-matched full cell, reflecting accurate electrochemical performance evaluations of high-capacity electrodes. Moreover, detailed electrochemical mechanisms of impedance and diffusion properties for working electrodes are also successfully decoupled individually. Here, this work uncovers the mismatch between traditional evaluation configuration and increasing testing current density and provides a guideline for accurate electrochemical evaluation for ever-increasing high-capacity electrodes, which is of great significance for high-energy lithium or other alkali-metal ion batteries.

25 ENERGY STORAGE↗

MESOSCALE THERMAL TRANSPORT MEASUREMENTS OF MULTI-PHASE AND POROUS NUCLEAR FUELS USING A SQUARE-WAVE PULSE THERMOREFLECTANCE TECHNIQUE

The safe and efficient operation of nuclear reactors require accurate knowledge of peak temperatures in the fuel assemblies. The temperature profiles are governed by the thermal transport properties of the fuel, namely the thermal conductivities (k) and thermal diffusivities (D). These values can be very difficult to measure as they are known to vary considerably from the measured bulk values of the fresh fuel, and quickly degrade with increasing burnup [1 -3] . Laser-based techniques have been effectively used for non-destructive and non-contact thermal transport measurements of a wide variety of materials, including nuclear materials [4-6] that would otherwise prove too hazardous or difficult to measure otherwise. In this study, a new thermoreflectance technique known as square-pulse transient thermoreflectance (SPTR) is described and used to determine the mesoscale thermal diffusivity of both uranium sesquisilicide (U3Si2) and uranium nitride (UN) phases in a composite fuel with micron level spatial resolution [7]. This technique employs a rapid train of square-wave pulses from an excitation laser to create a periodic heat flux on a gold coated sample surface. Surface adsorption results in transient film temperatures and hence rapid fluctuations in thermoreflectance that can be measured via a detection laser coupled with a digital oscilloscope. The lasers are coaxially focused on the sample surface, allowing for a sample measurement area of a single convolved laser spot size (~2 µm). A sensitivity analysis was conducted to identify key measurement parameters of this technique using reference materials with a range of thermal conductivities comparable to those of both ceramic, composite, and metal nuclear fuel types (1.4 – 27.2 W/m-K). The reference materials were measured using the new technique as well as a spatial-domain thermoreflectance technique (SDTR) previously reported for a comparison [8]. Additionally, measurements of several U3Si2 and UN phase regions of a polished UN/U3Si2 (70/30 vol.%) sample were taken, and the resulting calculated D values are reported, with both techniques showing excellent agreement between samples. This technique was used to scan a multiphase region at 5 micron increments to generate a local diffusivity map, demonstrating the utility of the techniques for measuring thermal transport properties in specimens with precipitates and secondary phases. Furthermore, the technique is currently being applied to measure thermal properties of Fast Flux Test Reactor (FFTR) irradiated metal fuel specimens whose porous microstructure make it very difficult to measure using other techniques. Comparison of the local scale measurements are compared with the pre-irradiated fuel samples to show the degradation of thermal transport in fuels due to pores from fission gas bubbles and other irradiation induced defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Towards a neutron multiplicity measurement with the Accelerator Neutrino Neutron Interaction Experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26 ton Gadolinium (Gd)-loaded water Cherenkov detector located on the Booster Neutrino Beam line at Fermilab. Its main goals are the measurement of the neutron multiplicity in neutrino-nucleus interactions as well as the cross-section of Charged Current Quasi-Elastic (CCQE) neutrino interactions on water. Besides the physics goals, the experiment also aims to be a testbed for new technologies such as Large Area Picosecond Photodetectors (LAPPDs) and Water-based Liquid Scintillators (WbLS). This thesis presents a preliminary measurement of the neutron multiplicity with {ANNIE}, using an analysis conducted on a fraction of the 2021 beam year. As preparatory measures, the efficiency of ANNIE's Front Muon Veto (FMV) was determined to be {$\bar{\varepsilon}_{\mathrm{FMV}} = (95.6 \pm 1.6)\%$} while the average efficiency for active scintillator paddles in the Muon Range Detector (MRD) was found to be {$\bar{\varepsilon}_ {\mathrm{MRD}} = (92.1 \pm 7.9)\%$}. Furthermore, the simulation framework used for ANNIE was validated and adapted to reproduce the experimental data by comparing the detector response for samples of Michel electrons, Americium Beryllium neutrons, and through-going muons. The analysis finds average neutron yields of {$\bar{n}_{\mathrm{data}} (\mathrm{beam}) = (0.272 \pm 0.010_{\mathrm{stat}})$} for an inclusive set of all identified muon neutrino candidates and {$\bar{n}_{\mathrm{data}} (\mathrm{beam,FV}) = (0.287 \pm 0.044_{\mathrm{stat}})$ for interactions which happened inside of the Fiducial Volume of ANNIE, which was optimized to increase the neutron detection acceptance. The presented neutron multiplicity values represent the number of detected neutrons after all event selection cuts and are not yet corrected for the neutron detection efficiency. An equivalent analysis on a simulated beam sample predicts neutron yields of $\bar{n}_{\mathrm{MC}}(\mathrm{beam}) = (0.515 \pm 0.0 07_{\mathrm{stat}})$ and $\bar{n}_{\mathrm{MC}}(\mathrm{beam,FV}) = (0.627 \pm 0.031_{\mathrm{stat}})$, indicating that the models tend to overpredict the number of neutrons produced in such interactions. Systematic errors have been briefly considered to contribute {$\sigma_{\mathrm{sys,FMV}} \sim 0.01\,$neutrons/$\nu$-interaction} due to the slight FMV inefficiency and {$\sigma_{\mathrm{sys,n}} \sim 0.05\,$neutrons/$\nu$-interaction} due to the neutron detection efficiency. Simulation studies further highlighted the importance of neutron detection in Diffuse Supernova Background (DSNB) searches. A combination of neutron tagging and Convolutional Neural Networks was found to reduce the most relevant Neutral Current Quasi-Elastic (NCQE) interaction background below the signal level, achieving a Signal-to-Background ratio of 4:1. In a further study, we investigated the positive impact of a deployment of a WbLS target on the energy reconstruction in ANNIE. WbLS provides a scintillation signal from hadronic recoils in addition to the charged lepton that can be included in neutrino energy reconstruction. It was found that a deployed WbLS volume in ANNIE improves the neutrino energy reconstruction from 14\% to 12\%, with the potential of going beyond this if more sophisticated reconstruction algorithms are developed in the future.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhancing carbon dioxide gas-diffusion electrolysis by creating a hydrophobic catalyst microenvironment

Electroreduction of carbon dioxide (CO 2 ) over copper-based catalysts provides an attractive approach for sustainable fuel production. While efforts are focused on developing catalytic materials, it is also critical to understand and control the microenvironment around catalytic sites, which can mediate the transport of reaction species and influence reaction pathways. Here, we show that a hydrophobic microenvironment can significantly enhance CO 2 gas-diffusion electrolysis. For proof-of-concept, we use commercial copper nanoparticles and disperse hydrophobic polytetrafluoroethylene (PTFE) nanoparticles inside the catalyst layer. Consequently, the PTFE-added electrode achieves a greatly improved activity and Faradaic efficiency for CO 2 reduction, with a partial current density >250 mA cm -2 and a single-pass conversion of 14% at moderate potentials, which are around twice that of a regular electrode without added PTFE. The improvement is attributed to a balanced gas/liquid microenvironment that reduces the diffusion layer thickness, accelerates CO 2 mass transport, and increases CO 2 local concentration for the electrolysis.

36 MATERIALS SCIENCE↗

Reexamination of hydrodynamic phonon transport in thin graphite

The recent hydrodynamic phonon transport theory for graphitic materials has been supported by the measurements of the second sound at temperatures up to about 100 K. When boundary scattering becomes comparable to momentum-conserving normal phonon scattering processes that are responsible for phonon hydrodynamics, Poiseuille phonon flow phenomena can emerge to give rise to unique size-dependent thermal conductivity in thin graphite. Here, we examine the thickness range for the Poiseuille phonon flow to become observable in thin graphite with the use of both deviational Monte Carlo simulation of the Peierls-Boltzmann transport equation and four-probe thermal transport measurements. As the basal-plane thermal conductivity calculated by prior first-principles theories saturates to the graphite value when the thickness is increased to five graphene layers, the phonon dispersion of graphite is used in the current calculations of thin graphite of micrometer thickness and a 23-layer thick ultrathin graphite (UTG) sample. The calculations show that diffuse surface scattering by surface defects can lead to Poiseuille phonon flow at 50 K in thin graphite with the thickness close to several micrometers but not in the 65 μm thin graphite and 23-layer UTG, where phonon scattering by the top and bottom surfaces become, respectively, much less and more frequent than the normal processes. In addition, the calculation results with the bulk graphite dispersion and diffuse surface scattering show decreased basal-plane thermal conductivity with decreasing thickness, opposite to recent thermocouple measurements of thin graphite samples. In comparison, the calculation results reveal that partially diffuse surface defect scattering can yield the four-probe measurement results of UTG samples, which are prepared here with an improved process to minimize surface contamination.

36 MATERIALS SCIENCE↗