Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “type transport”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

P-Type Electrochemical Doping Can Occur by Cation Expulsion in a High Performing Polymer for Organic Electrochemical Transistors

We investigate the mechanism of ion-dependent charge compensation during electrochemical oxidation (doping) of the model mixed ionic/electronic transporting polythiophene derivative poly(3-{[2-(2-methoxyethoxy)ethoxy]methyl}thiophene-2,5-diyl) (P3MEEMT). Using a combination of electrochemical quartz microbalance gravimetry and glow discharge optical emission spectroscopy, we show that charge compensation during polymer redox processes proceeds via a cation-dependent mechanism. For p-type polymer oxidation in certain electrolytes, charge compensation is achieved by both eventual injection of anions into the film, as well as initial expulsion of cations from the film. We compare doping mechanisms for a variety of electrolyte salts including potassium chloride, tetrabutylammonium chloride, potassium hexafluorophosphate (KPF6), and tetrabutylammonium hexafluorophosphate. For the electrolyte KPF6, both the cations and anions coexist in the water-swelled polymer even prior to application of electrical bias. Our data indicate that electrochemical doping (hole injection into the polymer and ionic charge compensation) proceeds via the following mechanism: (1) hydration of the neutral film by electrolyte (water, cations, anions), (2) cation (K+) expulsion from the film upon initial application of an oxidative bias, and (3) anion injection into the film at higher oxidation/doping levels (>~2 × 1020/cm3). Understanding the mechanism of charge compensation during the doping process should allow for the design of improved mixed ionic/electronic conductors for use in applications ranging from organic supercapacitors and redox flow batteries to bioelectronic sensors, thermoelectrics, and devices for neuromorphic computing.

Flagg, Lucas Q.↗
DOE OSTI2020

Transport in magnetically doped one-dimensional wires: can the helical protection emerge without the global helicity?

We study the phase diagram and transport properties of arbitrarily doped quantum wires functionalized by magnetic adatoms. The appropriate theoretical model for these systems is a dense one-dimensional Kondo lattice (KL) which consists of itinerant electrons interacting with localized quantum magnetic moments. We discover the novel phase of the locally helical metal where transport is protected from a destructive influence of material imperfections. Paradoxically, such a protection emerges without a need of the global helicity, which is inherent in all previously studied helical systems and requires breaking the spin-rotation symmetry. We explain the physics of this protection of the new type, find conditions, under which it emerges, and discuss possible experimental tests. Our results pave the way to the straightforward realization of the protected ballistic transport in quantum wires made of various materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Understanding the p-doping of spiroOMeTAD by tris(pentafluorophenyl)borane

The solid-state organization of photoabsorber, hole and electron transporting layers, and interfaces between them plays an important role in governing the performance and stability of emerging optoelectronic devices such as perovskite solar cells (PSCs). The molecular organic semiconductor (OSC) 2,2' ,7,7' -tetrakis [N,N-di(4-methoxyphenyl)amino]-9,9'-spirobifluorene (spiroOMeTAD) is a promising hole-transporting material (HTM) for PSCs, which is p-doped by molecular dopants to augment the charge carrier mobility. Here, the p-type doping of spiroOMeTAD by tris(pentafluorophenyl)borane (BCF) is investigated by a combination of techniques including optical spectroscopy, X-ray diffraction, Fourier transform infrared (FTIR), solid-state (ss)NMR, and electron paramagnetic resonance (EPR) spectroscopy. BCF molecules interact with traces of water molecules to form BCF-water complexes. Optical spectroscopy analysis suggests that the BCF/BCF-water complexes oxidize spiroOMeTAD molecules and facilitate p-type doping of spiroOMeTAD molecules. The different distributions of BCF and BCF-water molecules in doped spiroOMeTAD are characterized by FTIR and 11B NMR spectroscopy. An NMR crystallography approach which combines two-dimensional (2D) ssNMR and crystallography modeling is employed to unravel the packing interactions in spiroOMeTAD, and this analysis is extended to probe the morphological and structural changes in spiroOMeTAD:BCF blends. The hyperfine interactions are characterized by 2D hyperfine sub-level correlation (HYSCORE) spectroscopy. In this way, insight into the complex spiroOMeTAD:BCF blend morphology is obtained and compared for different dopant concentrations. Further, molecular-level analysis of doped HTMs enabled by this study has much wider relevance for further investigation, for example, chemical design and interfacial engineering of p-type doped HTMs for stable and efficient hybrid perovskite photovoltaics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Preliminary Analysis of Advanced Reactor Spent Nuclear Fuel Storage, Transportation, and Disposal

Due to increased interest in advanced reactor deployment and their associated potential new fuel cycles, the U.S. Department of Energy (DOE) Spent Fuel and Waste Science and Technology (SFWST) program has begun to evaluate the possible implications of long term management and final disposition of the spent nuclear fuel (SNF) generated. Safely managing and dispositioning this SNF, along with any other associated radioactive waste, is the primary focus of this initial preliminary assessment. This paper summarizes efforts to evaluate the characteristics and packaging options for three types of advanced reactor SNF forms: (1) tristructural isotropic (TRISO), (2) metallic, and (3) irradiated fuel salt presented in the report titled “Storage, Transportation, and Disposal of Advanced Reactor Spent Nuclear Fuel and High-Level Waste”. TRISO and metallic SNF and their associated waste streams were emphasized because of the near-term anticipated demonstrations of X-energy’s Xe-100 and TerraPower and GE Hitachi’s Natrium advanced reactors. Preliminary information on spent fuel salts discharged from molten-salt reactors (MSRs) was also examined to provide a baseline for future efforts. All calculations and assumptions were based on publicly available information. This paper identifies several different reactors that produce either TRISO or metallic SNF as well as a few of the reactor and fuel characteristics used for the assessments. Based on these characteristics, calculations were performed to determine the applicability of packaging SNF into existing canister designs. The evaluations included geometric (e.g., dimension, volume) and mass/weight considerations, known operational approaches and loading procedures, physical and chemical considerations/conditions for storage environments, and as-loaded radiation, thermal, and criticality analyses to identify constraints on storage, transportation, and disposal. Gaps in publicly available data pertaining to reactor operation and/or fuel composition provide increased uncertainty in some evaluations. Additionally, uncertainty in packaging and SNF management operations provide additional uncertainty. However, preliminary conclusions can still be assessed through this work and are presented in this paper.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

The Zintl pnictides Yb10CdSb9 and Yb14CdSb11: New candidate thermoelectric materials

The synthesis of new materials is the lifeline of solid-state science, and it continues to offer us unique opportunities for testing various theoretical formulations and models on a practical material. Such an avenue, therefore, provides a breeding ground for technological innovations and advancements that can completely revolutionize our world. Here, we report the results of our exploratory syntheses in the Yb–Cd–Sb compositional space that lead to the identification of two new Zintl antimonides, namely, Yb10CdSb9 and Yb14CdSb11. Their crystal structures were established via single-crystal X-ray diffraction methods; the basic electronic and transport properties of the new materials were also characterized. Yb10CdSb9 crystallizes in a disordered variant of the tetragonal Ca10LiMgSb9 structure type with unit cell parameters a = 11.8473(8) Å and c = 17.1302(12) Å (space group P42/mnm). Yb14CdSb11 crystallizes in the tetragonal Ca14AlSb11 structure type with unit cell parameters: a = 16.605(3) Å and c = 12.144(7) Å (space group I41/acd). Although the structures of both compounds can be rationalized within the framework of the Zintl formalism, based on the partitioning of the valence electrons in the much disordered Yb10CdSb9 phase, the charge is indicative of a slightly electron-rich composition. Electronic structure calculations in both cases support the notion of intrinsic semiconductor behavior, as expected for a Zintl phase. The temperature dependence of the electrical resistivity of a single crystal of Yb10CdSb9 is in line with that, and the evolution of the Seebeck coefficient indicates an electron-dominated transport mechanism, and a respectable power factor of 0.71 μW/cm K2 at 460 K can be calculated for Yb10CdSb9. The electrical resistivity of Yb14CdSb11, however, evolves in a semimetallic manner, which could suggest an overdoped sample or degenerate semiconducting behavior.

Ogunbunmi, Michael O. (ORCID:0000000253409198)↗
DOE OSTI2021

Enhancement of anticancer drug sensitivity in multidrug resistance cells overexpressing ATP-binding cassette (ABC) transporter ABCC10 by CP55, a synthetic derivative of 5-cyano-6-phenylpyrimidin

Highlights: • Multidrug resistance (MDR) in cancer is one of the major obstacles in chemotherapy. • ABCC10 mediates resistance to wide spectrum of chemotherapeutic agents, but only limited number of studies are available. • Overexpression of ABCC10 has been linked to various types of cancers including colorectal, gastric and lung cancer. • CP55 is a potent ABCC10 inhibitor which enhances the anticancer drug sensitivity in ABCC10-overexpression MDR cells. ATP-binding cassette (ABC) transporter C10 (ABCC10), also named multidrug resistance protein 7 (MRP7), is a member of ABC transporter superfamily and has been revealed to transport a wide range of chemotherapeutic agents including taxanes, epothilone B, Vinca alkaloids, and anthracyclines. In our previous study, a 5-cyano-6-phenylpyrimidin derivative CP55 was synthesized and found significantly reversal effect of multidrug resistance (MDR) mediated by ABCB1. In this study, we found CP55 also efficiently reversed MDR mediated by ABCC10. Our in vitro study showed that co-treatment with CP55 significantly increased the efficacy of ABCC10-substrate anticancer drugs in MDR cells overexpressing ABCC10. Furthermore, we showed that treatment with CP55 increased the intracellular accumulation of [{sup 3}H]-labeled anticancer drugs and in-turn decreasing drug efflux by inhibiting the transport activity, without altering ABCC10 protein ex-pression level or cellular localization. Potential CP55-ABCC10 interactions were predicted via docking analysis using human ABCC10 homology model and obtained high docking score. Therefore, CP55 represents a promising therapeutic agent in the combinational treatment of chemo-resistant cancer related to ABCC10.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Interplay between topological valley and quantum Hall edge transport

An established way of realising topologically protected states in a two-dimensional electron gas is by applying a perpendicular magnetic field thus creating quantum Hall edge channels. In electrostatically gapped bilayer graphene intriguingly, even in the absence of a magnetic field, topologically protected electronic states can emerge at naturally occurring stacking domain walls. While individually both types of topologically protected states have been investigated, their intriguing interplay remains poorly understood. Here, we focus on the interplay between topological domain wall states and quantum Hall edge transport within the eight-fold degenerate zeroth Landau level of high-quality suspended bilayer graphene. We find that the two-terminal conductance remains approximately constant for low magnetic fields throughout the distinct quantum Hall states since the conduction channels are traded between domain wall and device edges. For high magnetic fields, however, we observe evidence of transport suppression at the domain wall, which can be attributed to the emergence of spectral minigaps. This indicates that stacking domain walls potentially do not correspond to a topological domain wall in the order parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Effect of Aromatic Solvents Residuals on Electron Mobility of Organic Single Crystals

Balanced electron and hole transport properties are essential for various organic electrical/optoelectrical applications such as organic solar cells, complementary circuits, and light-emitting transistors. However, the electron transport in organic semiconductor lags far behind the hole side, making it with vital significance to seek the factors that limit the electron mobility. Here, the authors demonstrate that the π-conjugated solvents, as essential components in the widely-used solution processing techniques, can significantly suppress the electron transport if they are not properly removed and trapped as residuals. Single crystals of typical p-type materials exhibit a transition from p-type behavior to ambipolar conductance via reducing the amount of solvent residuals, and the same effect is confirmed by varied p-type materials and solvents. The highest electron mobility reaches 0.027 cm 2 V –1 s –1 and 0.029 cm 2 V –1 s –1 for 6,13-bis(triisopropylsilylethynyl)pentacene and 5,15-bis(triisopropylsilylethynyl)tetrabenzoporphyrin, respectively. Finally, this work discloses the non-negligible existence of solvent residuals even in high quality, long-range ordered organic single crystals, and further provides efficient strategies to harvest the n-type behaviors of organic semiconductor materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Layer-by-Layer-Processed Organic Solar Cells with 18.02% Efficiency Enabled by Regulating the Aggregation of Bottom Polymers

The fabrication of organic solar cells (OSCs) by a layer-by-layer (LBL) method has attracted growing attention in recent years. As already known, the pre-aggregates of conjugated polymers in solution have a profound impact on their microstructure morphology in films. Herein, by simply controlling the solution temperature and annealing processes, the pre-aggregation behavior of D18 polymer in solution can be fine-tuned and the microstructure of D18 bottom layer is well manipulated. The optimized D18 bottom layer can effectively regulate L8-BO upper-layer-forming suitable networks for efficient charge transportation. In addition, a vertical phase separation with a special D/D:A/A structure (P-i-N-type component distribution) is also formed. As a result, compared to the 16.43% power conversion efficiency (PCE) of the bulk heterojunction devices, such control enables bilayer OSC devices based on the polymer D18 and L8-BO to deliver an enhanced PCE of 18.02% with simultaneously improved short-circuit current density, open-circuit voltage, and fill factor. Importantly, it is also demonstrated in these results that the LBL deposition process utilizing the pre-aggregation of polymer and its fiber-network-forming ability is a very promising approach to improve charge dynamics, suppress carrier recombination, and fabricate highly efficient OSCs.

14 SOLAR ENERGY↗
DOE OSTI2021

Recombination rate analysis in long minority carrier lifetime mid-wave infrared InGaAs/InAsSb superlattices

Gallium is incorporated into the strain-balanced In(Ga)As/InAsSb superlattice system to achieve the same mid-wave infrared cutoff tunability as conventional Ga-free InAs/InAsSb type-II superlattices, but with an additional degree of design freedom to enable optimization of absorption and transport properties. Time-resolved photoluminescence measurements of InGaAs/InAsSb superlattice characterization- and doped device structures are reported from 77 to 300 K and compared to InAs/InAsSb. The low-injection photoluminescence decay yields the minority carrier lifetime, which is analyzed with a recombination rate model, enabling the determination of the temperature-dependent Shockley–Read–Hall, radiative, and Auger recombination lifetimes and extraction of defect energy levels and capture cross section defect concentration products. The Shockley–Read–Hall-limited lifetime of undoped InGaAs/InAsSb is marginally reduced from 2.3 to 1.4 μs due to the inclusion of Ga; however, given that Ga improves the vertical hole mobility by a factor of >10×, a diffusion-limited InGaAs/InAsSb superlattice nBn could expect a lower bound of 2.5× improvement in diffusion length with significant impact on photodetector quantum efficiency and radiation hardness. At temperatures below 120 K, the doped device structures are Shockley–Read–Hall limited at 0.5 μs, which shows promise for detector applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Unveiling defect-mediated carrier dynamics in monolayer semiconductors by spatiotemporal microwave imaging

The optoelectronic properties of atomically thin transition-metal dichalcogenides are strongly correlated with the presence of defects in the materials, which are not necessarily detrimental for certain applications. For instance, defects can lead to an enhanced photoconduction, a complicated process involving charge generation and recombination in the time domain and carrier transport in the spatial domain. Here, we report the simultaneous spatial and temporal photoconductivity imaging in two types of WS 2 monolayers by laser-illuminated microwave impedance microscopy. The diffusion length and carrier lifetime were directly extracted from the spatial profile and temporal relaxation of microwave signals, respectively. Time-resolved experiments indicate that the critical process for photoexcited carriers is the escape of holes from trap states, which prolongs the apparent lifetime of mobile electrons in the conduction band. As a result, counterintuitively, the long-lived photoconductivity signal is higher in chemical-vapor deposited (CVD) samples than exfoliated monolayers due to the presence of traps that inhibits recombination. Overall, our work reveals the intrinsic time and length scales of electrical response to photoexcitation in van der Waals materials, which is essential for their applications in optoelectronic devices.

42 ENGINEERING↗
DOE OSTI2025

Existence of a flat band in highly correlated metal BaCo 2 P 2

We report on the observation of a flat band situated at the Fermi level E F along with the structural, electrical transport, and magnetic properties of BaCo 2 P 2 that crystallizes in the ThCr 2 Si 2 -type body-center tetragonal structure. This compound has the largest inter-layer pnictide (Pn) distance d Pn–Pn as well as the largest c/a ratio among all the known ACo 2 Pn 2 (A = alkaline earth metal) compounds, where a and c are the tetragonal lattice parameters. Hence, the magnetic and electronic properties of this compound are expected to have a quasi-two-dimensional character. Despite the evidence of the presence of sizable magnetic interactions, magnetic susceptibility χ(T) of BaCo 2 P 2 does not show magnetic ordering down to 1.8 K. The material shows good metallic conduction with a large residual resistivity ratio ρ 300 K /ρ 1.8 K ≈ 70 and a Fermi liquid behavior at low temperature. Kadowaki–Woods ratio R KW of BaCo 2 P 2 suggests the presence of sizable electronic correlations within this system. Additionally, a large many-body enhancement of 2.3 of the experimental density of states D γ (E F ) over the band-structure D band (E F ) is inferred to arise from sizable electron-electron and/or electron-phonon interactions leading to a substantial deviation from the free-electron behavior.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Development of Design for the STS Extraction Magnet System

The Oak Ridge National Laboratory Second Target Station Project will enhance the Spallation Neutron Source by adding a new neutron source. The upgrade includes a 30% increase in beam energy and a 50% boost in beam current, doubling the accelerator's power capability to 2.8 MW. The Ring-to-Second-Target Beam Transport (RTST) system is vital in directing high-energy proton beams to the target. A key element of the RTST system is the extraction magnets, which are tasked with precise beam extraction and transport. Within the framework of this work, Fermilab is responsible for carrying out the development of 3 types of magnets: Pulsed Dipole, Large Aperture Quadrupole and Narrow Quadrupole.

Chemenok, Vitalii [Unlisted, US]↗
DOE OSTI2024

Recent Advanced Reactor Multiphysics Model Highlights in the Virtual Test Bed (VTB)

The Virtual Test Bed (VTB) host over 30 distinct simulations that showcase state-of-the art capabilities across the national lab complex. An update on the status of models on the VTB is summarized here, along with a more detailed overview of select recent new capabilities to showcase. All of the major advanced reactor types are represented in the VTB. The first example consists of a multiphysics simulation to track the transport of species in Molten Salt Reactors using depletion, advection, and thermochemical calculations. The second consists of a coupled neutronic and thermal hydraulic simulation to validate a gas cooled reactor. The third consist of pebble-bed equilibrium model for a fluoride high-temperature reactor. The fourth is a high-fidelity neutronic and thermal hydraulic model of a liquid metal reactor assembly. And lastly the fifth consists of transient multiphysics simulations of heat pipe microreactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Mitigation of the internal p-n junction in CoS 2 -contacted FeS 2 single crystals: Accessing bulk semiconducting transport

Pyrite FeS 2 is an outstanding candidate for a low-cost, nontoxic, sustainable photovoltaic material, but efficient pyrite-based solar cells are yet to materialize. Recent studies of single crystals have shed much light on this by uncovering a p-type surface inversion layer on n-type (S-vacancy doped) crystals, and the resulting internal p-n junction. This leaky internal junction likely plays a key role in limiting efficiency in pyrite-based photovoltaic devices, also obscuring the true bulk semiconducting transport properties of pyrite crystals. Here, we demonstrate complete mitigation of the internal p-n junction in FeS 2 crystals by fabricating metallic CoS 2 contacts via a process that simultaneously diffuses Co (a shallow donor) into the crystal, the resulting heavy n doping yielding direct Ohmic contact to the interior. Low-temperature bulk transport studies of controllably Co- and S-vacancy doped semiconducting crystals then enable a host of previously inaccessible observations and measurements, including determination of donor activation energies (which are as low as 5 meV for Co), observation of an unexpected second activated transport regime, realization of electron mobility up to 2100 cm 2 V –1 s –1 , elucidation of very different mobilities in Co- and S-vacancy-doped cases, and observation of an abrupt temperaturedependent crossover to bulk Efros-Shklovskii variable-range hopping, accompanied by an unusual form of nonlinear Hall effect. Aspects of the results are interpreted with the aid of first-principles electronic structure calculations on both Co- and S-vacancy-doped FeS 2 . Furthermore, this work thus demonstrates unequivocal mitigation of the internal p-n junction in pyrite single crystals, with important implications for both future fundamental studies and photovoltaic devices

36 MATERIALS SCIENCE↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗