Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Contact thermodynamics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Undercooling minimization in ultrasound coupled DTA measurements of molten salts

Molten salts are ionic liquids that are used for the electrolytic pyroprocessing of metals and as heat transfer fluids in very high temperature processes. Recently, halide-based molten salt reactors (MSR) have gained momentum for high density and environmentally responsible electricity generation. The function of these reactors and their fuel cycle depend on a knowledge of the halide salt’s thermodynamic properties [1]. Therefore, high accuracy phase equilibria of MSR relevant base halide and actinide containing salts are needed. Modern thermal analysis of the phase transitions of halide salts is usually done during heating at relatively high scan rates with commercial devices. Normally such measurements are adequate for pure or pseudo-binary salts. However, as the number of components in the mixture increases, accurately resolving the liquidus becomes increasingly difficult. Reversing the scanning mode greatly increases the sensitivity of phase transition measurements but can decrease accuracy due to undercooling [2] — a common occurrence in molten halide systems [3]. This work presents the development of a differential thermal analysis (DTA) cell intended for use in radiological gloveboxes. Non-contact ultrasonic agitation of the halide salt is implemented to limit kinetic limitations on crystallization during cooling to minimize undercooling [4]. Measurements on halide salts are also presented to elucidate the effect of non-contact mixing on undercooling. Reference [1] S. Boyd and C. Taylor, “3 - Chemical fundamentals and applications of molten salts,” in Molten Salt Reactors and Thorium Energy, T. J. Dolan, Ed., Woodhead Publishing, 2017, pp. 29–91. doi: 10.1016/B978-0-08-101126-3.00003-8. [2] K. Nitsch, A. Cihlár, and M. Rodová, “Molten state and supercooling of lead halides,” J. Cryst. Growth, vol. 264, no. 1, pp. 492–498, Mar. 2004, doi: 10.1016/j.jcrysgro.2004.01.011. [3] L. Rycerz, “Practical remarks concerning phase diagrams determination on the basis of differential scanning calorimetry measurements,” J. Therm. Anal. Calorim., vol. 113, no. 1, pp. 231–238, Jul. 2013, doi: 10.1007/s10973-013-3097-0. [4] Md. H. Zahir, S. A. Mohamed, R. Saidur, and F. A. Al-Sulaiman, “Supercooling of phase-change materials and the techniques used to mitigate the phenomenon,” Appl. Energy, vol. 240, pp. 793–817, Apr. 2019, doi: 10.1016/j.apenergy.2019.02.045.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗
DOE OSTI2023

Multiscale molecular simulations for the solvation of lignin in ionic liquids

Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (H E ) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and H E ) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet–Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA] + -based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA] + -based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.

09 BIOMASS FUELS↗
NASA NTRS1983

Limits to solar power conversion efficiency with applications to quantum and thermal systems

An analytical framework is presented that permits examination of the limit to the efficiency of various solar power conversion devices. Thermodynamic limits to solar power efficiency are determined for both quantum and thermal systems, and the results are applied to a variety of devices currently considered for use in space systems. The power conversion efficiency for single-threshold energy quantum systems receiving unconcentrated air mass zero solar radiation is limited to 31 percent. This limit applies to photovoltaic cells directly converting solar radiation, or indirectly, as in the case of a thermophotovoltaic system. Photoelectrochemical cells rely on an additional chemical reaction at the semiconductor-electrolyte interface, which introduces additional second-law demands and a reduction of the solar conversion efficiency. Photochemical systems exhibit even lower possible efficiencies because of their relatively narrow absorption bands. Solar-powered thermal engines in contact with an ambient reservoir at 300 K and operating at maximum power have a peak conversion efficiency of 64 percent, and this occurs for a thermal reservoir at a temperature of 2900 K. The power conversion efficiency of a solar-powered liquid metal magnetohydrodydnamic generator, a solar-powered steam turbine electric generator, and an alkali metal thermoelectric converter is discussed.

Byvik, C. E.↗
DOE OSTI2020

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗
DOE OSTI2022

Investigating Abiotic and Biotic Mechanisms of Pyrite Reduction

Pyrite (FeS 2 ) has a very low solubility and therefore has historically been considered a sink for iron (Fe) and sulfur (S) and unavailable to biology in the absence of oxygen and oxidative weathering. Anaerobic methanogens were recently shown to reduce FeS 2 and assimilate Fe and S reduction products to meet nutrient demands. However, the mechanism of FeS 2 mineral reduction and the forms of Fe and S assimilated by methanogens remained unclear. Thermodynamic calculations described herein indicate that H 2 at aqueous concentrations as low as 10 –10 M favors the reduction of FeS 2 , with sulfide (HS – ) and pyrrhotite (Fe 1–x S) as products; abiotic laboratory experiments confirmed the reduction of FeS 2 with dissolved H 2 concentrations greater than 1.98 × 10 –4 M H 2 . Growth studies of Methanosarcina barkeri provided with FeS 2 as the sole source of Fe and S resulted in H 2 production but at concentrations too low to drive abiotic FeS 2 reduction, based on abiotic laboratory experimental data. A strain of M. barkeri with deletions in all [NiFe]-hydrogenases maintained the ability to reduce FeS 2 during growth, providing further evidence that extracellular electron transport (EET) to FeS 2 does not involve H 2 or [NiFe]-hydrogenases. Physical contact between cells and FeS 2 was required for mineral reduction but was not required to obtain Fe and S from dissolution products. The addition of a synthetic electron shuttle, anthraquinone-2,6-disulfonate, allowed for biological reduction of FeS 2 when physical contact between cells and FeS 2 was prohibited, indicating that exogenous electron shuttles can mediate FeS 2 reduction. Here, transcriptomics experiments revealed upregulation of several cytoplasmic oxidoreductases during growth of M. barkeri on FeS 2 , which may indicate involvement in provisioning low potential electrons for EET to FeS 2 . Collectively, the data presented herein indicate that reduction of insoluble FeS 2 by M. barkeri occurred via electron transfer from the cell surface to the mineral surface resulting in the generation of soluble HS – and mineral-associated Fe 1–x S. Solubilized Fe(II), but not HS – , from mineral-associated Fe 1–x S reacts with aqueous HS – yielding aqueous iron sulfur clusters (FeS aq ) that likely serve as the Fe and S source for methanogen growth and activity. FeS aq nucleation and subsequent precipitation on the surface of cells may result in accelerated EET to FeS 2 , resulting in positive feedback between cell activity and FeS 2 reduction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Mirostructure Characterization of Friction Consolidated Copper-Nickel using a Machine Learning Approach: Developing Process to Microstructure Associations

Friction consolidation (FC) is a solid phase processing approach where discrete material forms such as powders, chips, nuggets, etc. are densified via shear deformation. The precursors are placed in a billet container and brought in contact with a rotating tool that applying the desirable amount of normal force. Under the combined action of the rotation and normal pressure, the discrete precursor is consolidated through porosity reduction and shear deformation. FC is increasingly being studied as an attractive approach to manufacturing fully dense parts from powder forms owing to its ability to mix, alloy and consolidate difficult-to-process precursors in minimal number of process steps. Material consolidation and deformation in shear consolidation processes have been studied extensively previously for different material combinations previously. However, despite the extensive research in this area, understanding of the mechanistic processes in pore consolidation, deformation-induced mixing and material solubility during FC is still evolving. Material development using solid phase processing approaches such as FC is often performed based on research experience/education, which can be biased. Conventional analysis and simulation tools in this area tend to be successful only when material thermodynamic pathways and microstructural evolution sequences resulting from processing are clearly defined or known. They are not as effective for emerging advanced manufacturing technologies where material evolution pathways are not well established. The ability to predict optimal process parameters based on material chemistry and bulk properties is essential to accelerate materials design and processing, as are an understanding of the relevant structure-processing-property relationships. These structure-processing-property-performance relationships are at the core of materials science research. Microstructure characterization provides the link to these four core areas, often through visualizing material microstructure using imaging techniques. However, linking microstructure image data (i.e., micrographs) to variables of interest (e.g., processing parameters, material chemistry) in a reproducible, generalizable, and quantitative manner is a significant challenge. Typically, quantitatively linking image data to processing history relies on significant domain knowledge and manual or subject matter expert (SME)-heuristic based image analysis. Such an approach to image analysis has the potential to be biased, inefficient, and difficult to replicate.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CO2 Capture Strategies via Mineralization with Industrial Waste Brines

Large coal-fired power plants (>500 MW) account for 30% of global CO2 emissions, and long-term management of this CO2 to is urgently needed mitigate global temperature increases. Sequestration of CO2 within stable mineral carbonates (e.g., CaCO3) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO2 in an environmentally friendly form. We have previously described a mineralization process in which divalent cations are sourced from various waste streams (e.g., produced water and brackish water) and alkalinity is induced via regenerable ion-exchange materials (Bustillos et. al. Frontiers in Energy Research. 2020, 8, 352). In our process, aqueous carbonate-bearing streams with pH > 8 are produced by contacting fresh water and carbon dioxide with various ion-exchange materials (e.g., Na form zeolites or ion exchange resins). These streams are mixed with produced water containing varying concentrations (~0.01 – 1.0 M) of Ca2+ leading to the precipitation of solid calcium carbonate (PCC). This process has the advantages of using regenerable solids in a simple and continuous process to increase the pH of water by ion exchange instead of relying on the consumption of costly and unsustainable sources of alkalinity (e.g., sodium hydroxide). While once-through column experiments showed the above benefits, the same were yet to established in a steady-state process with recycle streams. In this work, we set up a process simulation to quantify the energy requirements and CO2 emissions associated with the process and seek optimal produced water compositions and CO2 concentrations (5 – 20 vol%). The process simulation was set up in ASPEN Plus using eRNTL as the thermodynamic property method and sequential modular strategy. Ion exchange alkaline solution was simulated using sodium hydroxide and validated against the experimental data obtained from once-through kinetic experiments. Nanofiltration and reverse osmosis membrane steps were also implemented for the separation of divalent cations and production of fresh water and a regeneration stream following mineralization. Sensitivity analysis was carried out using a range of produced water compositions (0.01 – 1.0 M Ca2+, 0.001 – 0.15 M Mg2+, 0.5 – 3.5 M Na+ and 0.0004 – 0.002 M Fe2+) according to the United States Geological Survey (USGS) database. Calcium carbonate yields increased with increasing CO2 concentrations and were maximized using produced water compositions with larger Ca2+ concentrations. Maximum calcium carbonate yields produced at 5 vol%, 12 vol% and 20 vol% CO2 were 2.3 mmol/L, 5.5 mmol/L, and 9.3 mmol/L, respectively, with the formation of brucite (a magnesium hydroxide phase, Mg(OH)2) and goethite (an iron hydroxide phase, FeOOH) as the primary contaminant phases (99% calcite, 0.6% brucite, 0.4% goethite), which agree with phases detected by XRD experimentally. These results indicate high purity calcium carbonate can be precipitated using industrial waste streams. Consequentially, energy consumption and net CO2 emissions were minimized where precipitated calcium carbonate was maximized for all produced water compositions and CO2 concentrations. Minimum energy consumptions were 0.21 kWh/ton CO2 processed, with 98% of the energy input required coming from the membrane filtration steps. Produced water compositions with large Na+ concentrations (> 0.5 M) were effective at reducing energy consumptions due to faster regeneration time of ion exchange materials. Additionally, calculated net CO2 emissions were negative for the process and ranged from -0.02 kg/ton CO2 to -0.15 kg/ton CO2 processed, indicating a low emission process. We will also present techno-economic assessment showing the economic benefits of the current process as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO3.

Simonetti, Dante↗
DOE OSTI2025

Vacuum-assisted carbon molecular sieve membrane reactor for non-oxidative ethane dehydrogenation

Non-oxidative ethane dehydrogenation (EDH) is equilibrium-limited and endothermic. Selective hydrogen removal using a gas-permeable membrane within the EDH reaction zone can overcome the thermodynamic equilibrium, enabling higher ethane conversions. Employing vacuum as the permeation driving force, rather than a sweep gas, enhances the industrial viability of membrane reactors by eliminating additional post-reaction separation units. This study presents a membrane reactor that integrates H 2 -permeable carbon molecular sieve (CMS) hollow fiber membranes embedded in a fixed bed of cobalt in a dealuminated beta zeolite (Co@DeAl-BEA) catalyst, utilizing a vacuum to remove hydrogen efficiently. The CMS membrane exhibits high hydrogen permeance and an excellent H 2 /C 2 H 6 separation factor. The membrane reactor significantly enhanced the ethane conversion under reaction conditions comparable to those reported in the literature. A Langmuir-Hinshelwood kinetic rate expression was developed and incorporated into a one-dimensional steady-state reactor model. The experimentally validated model indicates that increasing the number of hollow fibers improves ethane conversion, although ethane loss to the permeate limits the benefit. The contact area between the catalyst and the membrane limits the reactor performance more than the catalytic throughput. Furthermore, we find that the location of the catalyst packing relative to the hollow fiber membranes influences ethane loss and conversion. Higher reactor pressures and inlet ethane flow rates improve space-time yield at the expense of lower ethane conversion. Increasing reactor temperature or packing length promotes both performance metrics. The EDH membrane reactor demonstrated durability over 200 h of continuous operation, maintaining record-low deactivation rates and high ethylene selectivity. Protocols for catalyst regeneration were developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Molecular Dynamics Simulations of Calcite Fracture in Water

Calcite (CaCO 3 ) is one of the most common minerals in geologic and engineered systems. It is often in contact with aqueous solutions, causing chemically assisted fracture that is critical to understanding the stability of subsurface systems and manmade structures. Calcite fracture was evaluated with reactive molecular dynamics simulations, including the impacts of crack tip geometry (notch), the presence of water, and surface hydroxyl groups. Chemo-mechanical weakening was assessed by comparing the loads where fracture began to propagate. Our analyses show that in the presence of a notch, the load at which crack growth begins is lower, compared to the effect of water or surface hydroxyls. Additionally, the breaking of two adjacent Ca–O bonds is the kinetic limitation for crack initiation, since transiently broken bonds can reform, not resulting in crack growth. In aqueous environments, fresh (not hydroxylated) calcite surfaces exhibited water strengthening. Manual addition of H + and/or OH – species on the (104) calcite surface resulted in chemo-mechanical weakening of calcite by 9%. Achieving full hydroxylation of the calcite surface was thermodynamically and kinetically limited, with only 0.17–0.01 OH/nm 2 surface hydroxylation observed on the (104) surface at the end of the simulations. In conclusion, the limited reactivity of pure water with the calcite surface restricts the chemo-mechanical effects and suggests that reactions between physiosorbed water and localized structural defects may be dominating the chemo-mechanical process in the studies where water weakening has been reported.

Calcite↗
DOE OSTI2021

Iron speciation in blast furnace slag cements

Slag-containing pastes and concretes were analysed by element-specific synchrotron-based techniques to determine the speciation of iron on crushed materials through spatially resolved micro-spectroscopic studies. The investigated cement samples were hydrated either in the laboratory, or exposed to river or sea water. Metallic iron, along with minor proportions of iron sulphide and magnetite was detected in the laboratory sample. Iron sulphide, goethite, and siliceous hydrogarnet were discovered in the blended slag cements hydrated in contact with river water for up to 7 years. In contrast, no Fe(0) was observed in blended concretes exposed to sea water. Instead, iron sulphide, iron(II)-hydroxide and -oxide, hematite, magnetite, siliceous hydrogarnet, and goethite were detected as well as ilmenite (FeTiO{sub 3}) in the aggregates. The strong acceleration of Fe oxidation in samples exposed to sea water and the long-term passivation observed in the other samples indicate comparable processes as those occurring on steel bars.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

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

Electron Spin Polarization and Rectification Driven by Chiral Perylene Diimide-Based Nanodonuts

The chirality-induced spin selectivity (CISS) effect allows thin-film layers of chiral conjugated molecules to function as spin filters at ambient temperature. Through solvent-modulated dropcasting of chiral l- and d-perylene diimide (PDI) monomeric building blocks, two types of aggregate morphologies, nanofibers and nanodonuts, may be realized. Spin-diode behavior is evidenced in the nanodonut structures. Stacked PDI units, which form the conjugated core of these nanostructures, dominate the nanodonut–Au electrode contact; in contrast, the AFM tip contacts largely the high-resistance solubilizing alkyl chains of the chiral monomers that form these nanodonuts. Current–voltage responses of the nanodonuts, measured by magnetic conductive AFM (mC-AFM), demonstrate substantial spin polarizations as well as spin current rectification ratios (>10) that exceed the magnitudes of those determined to date for other chiral nanoscale systems. Furthermore, these results underscore the potential for chiral nanostructures, featuring asymmetric molecular junctions, to enable CISS-based nanoscale spin current rectifiers.

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

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
NASA NTRS1983

Improved fabrication techniques for infrared bolometers

Ion implantation and sputter metallization are used to produce ohmic electrical contacts to Ge:Ga chips. The method is shown to give a high yield of small monolithic bolometers with very little low-frequency noise. It is noted that when one of the chips is used as the thermometric element of a composite bolometer it must be bonded to a dielectric substrate. The thermal resistance of the conventional epoxy bond is measured and found to be undesirably large. A procedure for soldering the chip to a metallized portion of the substrate in such a way as to reduce this resistance is outlined. An evaluation is made of the contribution of the metal film absorber to the heat capacity of a composite bolometer. It is found that the heat capacity of a NiCr absorber at 1.3 K can dominate the bolometer performance. A Bi absorber possesses significantly lower heat capacity. A low-temperature blackbody calibrator is built to measure the optical responsivity of bolometers. A composite bolometer system with a throughput of approximately 0.1 sr sq cm is constructed using the new techniques. The noise in this bolometer is white above 2.5 Hz and is slightly below the value predicted by thermodynamic equilibrium theory.

Lange, A. E.↗
NASA NTRS1994

Troilite formation kinetics and growth mechanism in the solar nebula

Troilite formation via the reaction Fe(s) + H2S(g) + H2(g) is the major mechanism for S retention in grains in the solar nebula. Thermodynamic calculations predict that troilite condenses from a solar composition gas. We present experimental results on the kinetics and growth of troilite crystals on Fe metal at temperature (450-650 C) and composition (50-1000 ppm H2S in H2) conditions similar to those in the solar nebula. The fraction of Fe reacted (based on gravimetric data) is plotted at 450, 505, 575, and 650 C. The thickness change of unreacted iron (measured by optical microscopy) is plotted at 575 and 650 C vs. time. the weight change per unit area varies as the square root of time at the lower temperatures and varies linearly with time at the highest temperature. The growth behavior along the lower isotherms is due to diffusion. This behavior suggests sulfide diffusion to the metal-sulfide interface and suggests Fe(2+) diffusion to the sulfide-gas interface. The reaction along the highest isotherm appears to be interface controlled. The formation of troilite crystals is a rapid process forming measurable layers in a few hours. The crystalgrowth is complicated. Initially there are intergrowths of troilite into the pure Fe metal. As the reaction progresses two distinct layers of troilite crystals form. One is in contact with the Fe metal and consists of small randomly oriented crystals with pore space between them. The outermost layer contains large crystals that are all oriented in the same direction. The intergrowth layer is much smaller at 650 C than at 575 C. This suggest that FeS nucleation is inhibited at the higher temperature, accounting for the initially slower reaction rate. Once nucleated, the reaction kinetics are apparently controlled by the growth of the crystals at the interface.

Lauretta, D. S.↗
DOE OSTI2026

The Complexation Properties of Self-Defensive Microgel-Modified Antimicrobial Surfaces

The complexation of cationic antimicrobials with polyanionic microgels on a biomaterial surface can render that surface self-defensive against bacteria by killing those bacteria which physically contact the antimicrobial-loaded microgels. This killing has been attributed to the contact-driven transfer of antimicrobial from a microgel to a challenging bacterium, though much remains unknown about this process. Here, in this study, we use a combination of experiments and computational modeling to identify key aspects of the complexation phenomena which influence the self-defensive properties. We synthesize poly(acrylic acid) (PAA) microgels (∼2–5 μm diameter) via membrane emulsification and electrostatically deposit them onto polycaprolactone (PCL) coupons or onto glass to form a discontinuous submonolayer. Subsequent microgel loading with colistin or with Sub5 antimicrobial peptide (AMP) causes microgel deswelling. Under physiological conditions Sub5 remains stably sequestered whereas colistin is quickly released. Coarse-grained molecular dynamics (CGMD) simulations confirm stronger Sub5/PAA complexation. CGMD calculations also indicate that Sub5 forms dimers and higher-order structures, a prediction confirmed experimentally by Small-Angle X-ray Scattering (SAXS). Supramolecular structure entropically enhances the complexation strength because of enhanced counterion release per complexation event, and this finding can help identify other antimicrobials well suited for such a nonelutive yet self-defensive strategy. CGMD simulations also show that Sub5 has a higher complexation strength with the Staphylococcus aureus membrane than it does with PAA, confirming that there is a thermodynamic driving force for antimicrobial transfer. Such self-defensive surfaces significantly reduce S. aureus colonization (over 90% reduction relative to unmodified controls) in an in vitro hematogenous contamination model and remain cyto-compatible as evidenced by mesenchymal stem cell spreading and proliferation.

36 MATERIALS SCIENCE↗
NASA NTRS1983

Improved fabrication techniques for infrared bolometers

Techniques are described for producing improved infrared bolometers from doped germanium. Ion implantation and sputter metalization have been used to make ohmic electrical contacts to Ge:Ga chips. This method results in a high yield of small monolithic bolometers with very little low-frequency noise. When one of these chips is used as the thermometric element of a composite bolometer, it must be bonded to a dielectric substrate. The thermal resistance of the conventional epoxy bond has been measured and found to be undesirably large. A procedure for soldering the chip to a metalized portion of the substrate is described which reduced this resistance. The contribution of the metal film absorber to the heat capacity of a composite bolometer has been measured. The heat capacity of a NiCr absorber at 1.3 K can dominate the bolometer performance. A Bi absorber has significantly lower heat capacity. A low temperature blackbody calibrator has been built to measure the optical responsivity of bolometers. A composite bolometer system with a throughput of approx. 0.1 sr sq cm was constructed using the new techniques. In negligible background it has an optical NEP of 3.6 10((exp -15) W/sq root of Hz at 1.0 K with a time constant of 20 ms. The noise in this bolometer is white above 2.5 Hz and is somewhat below the value predicted by thermodynamic equilibrium theory. It is in agreement with calculations based on a recent nonequilibrium theory.

Lange, A. E.↗