SEARCH · Search NASA

Results for “Single-Event Effects”

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.

SourceResearch connectionFrom yearThrough year

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

1,098 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS2026

Glenn Research Center Propulsion Systems Laboratory 2026 Customer Guide

This guide describes the Propulsion Systems Laboratory (PSL) at the NASA Glenn Research Center. It was written to help customers understand the various components involved in conducting a test program within the PSL. The PSL complex supports two large-engine test cells that simulate altitude flight conditions for a wide range of research and experimental tests. These test cells operate at altitudes up to 90,000 ft and speeds from subsonic to above supersonic. Test points such as pressure, temperature, and Mach number can be set at the engine or test article inlet based on customer requirements. The facility’s support systems include the heated and cooled combustion air systems; altitude exhaust system; hydraulic system; nitrogen, oxygen, and hydrogen systems; thrust measurement system, which includes the facility’s single- and multi-axis thrust stands; inlet system; and electrical systems. In addition to providing a detailed description of PSL systems and capabilities, this guide discusses the facility’s history and past tests and addresses facility safety procedures, pretest requirements, and test operation standards.

Engine Icing Test Facility↗

NASA NTRS1977

The Effect of Temperature and Pressure on the Distribution of Iron Group Elements Between Metal and Olivine Phases in the Process of Differentiation of Protoplanetary Material

The distribution patterns of Ni, Co, Mn, and Cr were studied in olivines of various origins: from meteorites (chondrites, achondrites, pallasites), which are likely analogs of the protoplanetary material, to peridotite inclusions in kimberlite pipes, which are analogs of mantle material. According to X-ray microanalysis data, each genetic group of olivines is characterized by a specific concentration of these elements. Nickel is concentrated (up to 0.34 percent) in peridotite olivines, while manganese is concentrated in meteoritic olivines. The maximum chromium content (0.2 percent) was found in ureilites, which were formed under reducing conditions. Experiments at pressures of 20 to 70 kbar and temperatures of 1100 to 2000°C have shown that in a mixture of olivine and Ni metal or NiO nickel enters the silicate phase (up to 4 percent), displacing Fe into the metallic phase. Equilibrium temperatures were estimated from the Fe, Ni distribution coefficients between the metal and olivine: 1500 K for pallasites, 1600 K for olivine-bronzite H6 chondrites, 1200 K for olivine-hypersthene L6, 900 K for LL6, and 1900 K for ureilites (at P = 1 atm). The equilibrium conditions of peridotites are close to T = 1800 K and P over 100 kbar. The distribution patterns of the transition elements are explained on the basis of physical-chemical properties. It is concluded that there is a sharp difference between the conditions of differentiation of the protoplanetary material at the time meteorites were formed and the conditions of differentiation of the planets into concentric layers.

A P Vinogradov↗

DOE OSTI2026

Topographic effects on reflected acoustic waves from the OSIRIS-REx reentry observed from stratospheric balloons

During long-distance sound propagation in planetary atmospheres, acoustic waves may reflect off the air/surface interface one or more times. For low sound frequencies and flat interfaces, the incident and reflected wave tend to be nearly identical. However, this may not be the case when the downgoing acoustic wave encounters topography. Here, we describe a set of direct and reflected acoustic signals recorded on free-flying balloons during the hypersonic entry of the OSIRIS-REx sample return capsule (SRC). In two of the three cases presented here, an impulsive reflected arrival similar in form to the direct sonic boom of the SRC is observed, followed by a diffuse coda. In contrast, one of the floating stations lacked an impulsive reflection entirely, with only the coda present. We use full waveform modeling to show how reflection in the presence of complex topography can explain coda seen in all three examples as well as the lack of impulsive arrival on the third. Our results indicate that the complex signals often observed in long range acoustic propagation could be due, in part, to interactions with topography during transmission.

Lees, Jonathan M. [University of North Carolina, C↗

NASA NTRS2026

Moon Base Transportation - Deliveries to the Lunar Surface

Development of the Moon Base will enable a home away from home for astronauts who will live and work at humanity’s first lunar outpost. In this effort, NASA’s Moon Transportation Office is responsible for enabling the transformational missions required to deliver habitats, supplies, science payloads, and all other elements needed to cultivate a permanent presence on the Lunar Surface. The Mission Concept (MC) is characterized through evaluation of an end-to-end architecture that can successfully deliver a generalized heavy large-volume payload, in excess of 4000 kg, to a precision landing and touchdown on the lunar surface. The mission architecture utilizes a single launch configuration of a Lunar Lander (LL) with a unique propellant system. The LL has an integral orbital transfer capability and features jettisonable elements. The design circumvents the need for prop transfer on orbit and multiple launch configurations. The launch vehicle (LV) for this work will assume the capability to deliver a payload in excess of 40,000 kg to orbit, affording multiple LV solutions. Considerations for the LL and payload deployment from the fairing are assumed to be handled through compliance with a launch providers’ Interface Requirements Document (IRD). The MC will span from launch at Kennedy Space Center (KSC) to terminal descent and touchdown on the lunar surface, requiring a total ΔV on the order of 6 km/s beyond what is required to get the vehicle stack to a 200 km circular Low Earth Orbit (LEO). Major mission phases include: launch and launch vehicle separation, transfer operations, pre-landing navigation, and lunar descent and touchdown. A Concept of Operations (ConOps) is used as the primary design driver for defining architecture of the vehicles necessary to achieve final payload delivery. Numerous ground rules and assumptions will be provided for each phase of the mission. Concept designs for the LL is presented. An emphasis of the design maximizes a feasible path for maturation, manufacturing, and operation. A self-imposed practical consideration for this effort is the incorporation of legacy designed hardware to minimize expensive, time-intensive, and high-risk hardware development cycles. The propulsion system of the LL adopts a conventional storable bipropellant configuration of monomethyl hydrazine (MMH) and mixed oxides of nitrogen (MON3). This effort will showcase a unique propellant delivery system to minimize the reliance on propellant management devices (PMDs) during descent. Numerous key constraints have been considered, across the multiple segments of the mission. These include the unique aspects of center of gravity (CG) management, thruster plume effects including self-impingement, propulsion system hardware limitations, navigation during multiple mission phases, and landing gear geometry for uneven terrain. These constraints shape the trades necessary for precision landing of heavy cargo and ensure compatibility with broader Moon Base Transportation concepts. The resulting insights inform future transportation strategies for the Moon and beyond; directly contributing to the development of cargo‑delivery standards that will support the long‑term buildup of a sustainable, continuously inhabited Moon Base.

Lunar Habitat↗

DOE OSTI2026

Effects of Composition and Oxidation States on the Structures of Chromium-Containing Sodium Silicate Glasses: Molecular Dynamics Simulations using Machine Learning Interatomic Potentials

Chromium represents a significant challenge for the vitrification of high-level nuclear waste into silicate and borosilicate glasses due to its low solubility and variable oxidation states, which can limit the waste loading due to promotion of crystallization or phase separation during processing. In this study, we modeled chromium containing silicate glasses using molecular dynamics simulations with three machine learning interatomic potentials (MLIPs), MACE, CHGNet, and PFP were employed, to gain insights on glass composition and oxidation states on the structures of these glasses. One of the goals is to evaluate their ability of these MLIPs to accurately represent the general structure of silicate glasses and chromium local environments as a function of chromium oxidation states. Density Functional Theory (DFT) based calculations and experimental data such as neutron structure factors were used to validate the structural models. It was found that the foundation models of all three MLIPs are able to reproduce general structural features of the sodium silicate glass structure consistent with experimental and DFT data, but only CHGNet and PFP can accurately capture the oxidation states and local environment of chromium: tetrahedral for Cr6+ and octahedral for Cr3+. Furthermore, we studied the effect of varying Cr3+/ Cr6+ (Cr3+/Crtotal) ratio and total chromium content using PFP. Our results show that Cr6+ enhances network polymerization by reducing non-bridging oxygens through Na? charge compensation required due to the formation of chromate (CrO42-) species, while Cr³? acts as a network modifier that disrupts connectivity. System size effects on the structural characteristics and chromium environments were also tested using the PFP potential. This work highlights the importance of careful validation on the precision, transferability, and potential of MLIPs for modeling glasses containing transition metal elements that can exist in multiple oxidation states. It is also encouraging to see the foundational models are all three MLFFs are able to reproduce the basic sodium silicate glass structures, while suggesting additional training or refining is needed to improve the description of more complex systems containing transition metals.

Puga, Christina L.↗

NASA NTRS1985

A Comparison of Doppler Lidar Wind Sensors for Earth-Orbit Global Measurement Applications

Now, there are four Doppler lidar configurations which are being promoted for the measurement of tropospheric winds: (1) the coherent CO2 Lidar, operating in the 9 micrometer region using a pulsed, atmospheric pressure CO2 gas discharge laser transmitter, and heterodyne detection; (2) the coherent Neodymium doped YAG or Glass Lidar, operating at 1.06 micrometers, using flashlamp or diode laser optical pumping of the solid state laser medium, and heterodyne detection; (3) the Neodymium doped YAG/Glass Lidar, operating at the doubled frequency (at 530 nm wavelength), again using flashlamp or diode laser pumping of the laser transmitter, and using a high resolution tandem Fabry-Perot filter and direct detection; and (4) the Raman shifted Xenon Chloride Lidar, operating at 350 nm wavelength, using a pulsed, atmospheric pressure XeCl gas discharge laser transmitter at 308 nm, Raman shifted in a high pressure hydrogen cell to 350 nm in order to avoid strong stratospheric ozone absorption, also using a high resolution tandem Fabry-Perot filter and direct detection. Comparisons of these four systems can include many factors and tradeoffs. The major portion of this comparison is devoted to efficiency. Efficiency comparisons are made by estimating the number of transmitted photons required for a single pulse wind velocity estimate of + or - 1 m/s accuracy in the middle troposphere, from an altitude of 800 km, which is assured to be reasonable for a polar orbiting platform.

Robert T Menzies↗

DOE OSTI2026

Associative polymers with controlled sticker placement: How reversible bond distribution and density govern polymer dynamics

Associative polymers with precisely arranged stickers offer opportunities to program material properties with molecular precision. Yet, it remains unclear how the placement and fraction of stickers dictate structure, dynamics, and macroscopic properties. By developing a model unentangled polymer system with hydrogen-bonding stickers, we show that randomly distributed stickers neither form clusters nor change flow properties, whereas stickers placed at chain ends drive nanocluster formation even at low concentrations. Adding more end stickers produces a rubbery plateau spanning eight decades in frequency with two distinct relaxation timescales, in contrast to the single plateau predicted by the classic sticky Rouse model. These results demonstrate that sticker distribution dictates whether associative polymers undergo nanocluster formation or microphase separation, while substantial alterations in dynamics and viscoelasticity require both sticker aggregation and thermomechanical stability of associated domains. Our findings resolve a longstanding debate on associative polymer dynamics and provide molecular design rules for programmable soft materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Nanostructuring and Underscreening of Bisalt Electrolytes with Dual-Anion Effects: Insights from Small-Angle Scattering Prepeak Analysis

Using small-angle neutron scattering, we studied the nanostructure of a prototypical D(H)-bonded network electrolyte, alkaline sodium aluminate bisalt, at concentrations up to solute volume fraction of ~0.5. Analysis of the structure factor prepeak at 0.1 < Q < 1.0 Å-1 showed that its evolution is associated with nanoscopic species distribution in water-anion network., with differences in characteristic distance d between OD- and Al(OD)4- anions related to their distinct sizes and interactions. When the solute volume fraction approached 0.5, a common/maximum value of ~4.5 Å was found for the correlation length ? that characterizes the electrostatic force in concentrated electrolyte solution before precipitation occurred. This implies that the solutions’ morphology and behavior at high concentrations may be governed by geometric factors, rather than the chemistry of the specific anion. Furthermore, the scaling of ? with concentration yielded an exponent of 1.25(1), suggesting that ? is limited to a few Å. By combining the current observations of structural heterogeneity at the nanoscale, with dynamic heterogeneity at the microscopic scale from our previous quasi-elastic neutron scattering study, we have established a structural origin of local “caging” and restricted structural relaxation processes. These local solvent-solute interactions not only control dynamics heterogeneity in concentrated electrolytes but also are responsible for crystallization processes in industrial setting, such as aluminum production and radioactive waste treatment.

Wang, Hsiu-Wen↗

NASA NTRS1994

Surface Characterization Through Shape Oscillations of Drops in Microgravity and 1-g

The goal of these experiments is to determine the rheological properties of liquid drops of single or multiple components in the presence or absence of surface active materials by exciting drops into their quadrupole resonance and observing their free decay. The resulting data coupled with appropriate theory should give a better description of the physics of the underlying phenomena, providing a better foundation than earlier empirical results could. The space environment makes an idealized geometry available (spherical drops) so that theory and experiment can be properly compared, and allows a 'clean' environment, by which is meant an environment in which no solid surfaces come in contact with the drops during the test period. Moreover, by considering the oscillations of intentionally deformed drops in microgravity, a baseline is established for interpreting surface characterization experiments done on the ground by other groups and ours. Experiments performed on the United States Microgravity Laboratory Laboratory (USML-1) demonstrated that shape oscillation experiments could be performed over a wide parameter range, and with a variety of surfactant materials. Results, however, were compromised by an unexpected, slow drop tumbling, some problems with droplet injection, and the presence of bubbles in the drop samples. Nevertheless, initial data suggests that the space environment will be useful in providing baseline data that can serve to validate theory and permit quantitative materials characterization at 1-g.

Robert E Apfel↗

NASA NTRS2023

Cryogenic Flow Boiling in Microgravity: Effects of Reduced Gravity on Two-Phase Fluid Physics and Heat Transfer

With the growing interest in space exploration, cryogenic technologies involving two-phase flow and heat transfer are in high demand to successfully procure advanced space applications such as fuel depots and nuclear thermal propulsion (NTP) systems for deep space missions. However, the unique and extreme thermal properties of cryogenic fluids introduce distinct flow boiling fluid physics and energy transport phenomena, which differ significantly from those observed with conventional fluids. Understanding the unique two-phase physics in cryogenic flow boiling remains an ongoing challenge. Furthermore, the lack of readily available microgravity cryogenic steady-state heat transfer data hinders the assessment of gravitational effects on cryogenic flow boiling. This study aims to elucidate the gravitational effects on two-phase fluid physics and heat transfer by conducting the first-ever experimental measurement of cryogenic flow boiling performance using a steady-state heated method in a reduced gravity environment. Parabolic flight experiments were performed to acquire both heat transfer measurements and high-speed video of interfacial behaviors, under varying gravity levels (microgravity, hypergravity, Lunar gravity, and Martian gravity). The experiments involved flow boiling of liquid nitrogen (LN 2 ) with a near-saturated inlet along a circular heated tube of dimensions 8.5-mm inner diameter and 680-mm heated length. The operating parameters varied are mass velocity of 398.3 - 1342.8 kg/m2s, inlet quality of -0.08 to -0.01, and inlet pressure of 413.68 - 689.48 kPa. Captured microgravity flow patterns range from bubbly to annular, all having vapor structures that are larger than those under higher gravity levels. Under microgravity, absence of buoyancy yields symmetrical vapor structures without flow stratification, laying a physical foundation for the distinct two-phase heat transfer trends during LN 2 flow boiling in microgravity. Transient data collected during the flight parabolas exhibited decreasing heated wall temperature as the aircraft transitioned from hypergravity to microgravity phases. The temperature variation indicated an enhancement in flow boiling heat transfer with decreasing gravity levels and a reduction with increasing gravity levels. The effect of reduced gravity on cryogenic flow boiling heat transfer coefficient (HTC) is discussed based on steady state heat transfer analysis. Seminal HTC correlations are evaluated against the measured microgravity HTC data, of which one is identified for superior accuracy in predicting microgravity data. Finally, a new HTC correlation is proposed to improve accuracy of microgravity predictions, yet there still exists room for further improvement with future terrestrial flow boiling experiments at different flow orientations relative to Earth gravity.

Microgravity↗

NASA NTRS2004

New Lithium-ion Polymer Battery for the Extravehicular Mobility Unit Suit

The Extravehicular Mobility Unit (EMU) suit currently has a silver-zinc battery that is 20.5 V and 45 Ah capacity. The EMU's portable life support system (PLSS) will draw power from the battery during the entire period of an EVA. Due to the disadvantages of using the silver-zinc battery in terms of cost and performance, a new high energy density battery is being developed for future use, The new battery (Lithium-ion battery or LIB) will consist of Li-ion polymer cells that will provide power to the EMU suit. The battery design consists of five 8 Ah cells in parallel to form a single module of 40 Ah and five such modules will be placed in series to give a 20.5 V, 40 Ah battery. Charging will be accomplished on the Shuttle or Station using the new LIB charger or the existing ALPS (Air Lock Power Supply) charger. The LIB delivers a maximum of 3.8 A on the average, for seven continuous hours, at voltages ranging from 20.5 V to 16.0 V and it should be capable of supporting transient pulses during start up and once every hour to support PLSS fan and pump operation. Figure 1 shows the placement of the battery in the backpack area of the EMU suit. The battery and cells will undergo testing under different conditions to understand its performance and safety characteristics.

Jeevarajan, J. A.↗

DOE OSTI2024

High‐Loading Lithium‐Sulfur Batteries with Solvent‐Free Dry‐Electrode Processing

Abstract Lithium‐sulfur (Li‐S) batteries, with their high energy density, nontoxicity, and the natural abundance of sulfur, hold immense potential as the next‐generation energy storage technology. To maximize the actual energy density of the Li‐S batteries for practical applications, it is crucial to escalate the areal capacity of the sulfur cathode by fabricating an electrode with high sulfur loading. Herein, ultra‐high sulfur loading (up to 12 mg cm −2 ) cathodes are fabricated through an industrially viable and sustainable solvent‐free dry‐processing method that utilizes a polytetrafluoroethylene binder fibrillation. Due to its low porosity cathode architecture formed by the binder fibrillation process, the dry‐processed electrodes exhibit a relatively lower initial capacity compared to the slurry‐processed electrode. However, its mechanical stability is well maintained throughout the cycling without the formation of electrode cracking, demonstrating significantly superior cycling stability. Additionally, through the optimization of the dry‐processing, a single‐layer pouch cell with a loading of 9 mg cm −2 and a novel multi‐layer pouch cell that uses an aluminum mesh as its current collector with a total loading of 14 mg cm −2 are introduced. To address the reduced initial capacity of dry‐processed electrodes, strategies such as incorporating electrocatalysts or employing prelithiated active materials are suggested.

Chemistry↗

DOE OSTI2026

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗

DOE OSTI2026

Design Trade-Offs in Composite Fuel Cell Membranes: Effects of Reinforcement and Chemical Additives

Perfluorosulfonic acid (PFSA) membranes are critical components in proton exchange membrane fuel cells, where performance depends on balancing ionic conductivity, mechanical durability, and chemical stability. This study characterizes a composite membrane (NC700) featuring PFSA-impregnated expanded polytetrafluoroethylene (ePTFE) reinforcement and cerium-based radical scavengers, benchmarked against unreinforced NR211. Complementary techniques, including electron microscopy, X-ray scattering, infrared spectroscopy, thermogravimetric analysis, and dynamic mechanical analysis, identify the structural and compositional strategies employed in NC700. Water sorption isotherms reveal lower water uptake for NC700 across all conditions, attributed to reinforcement and cerium incorporation. Reinforcement reduces in-plane swelling from 11% to 2.1% at 90% RH, confirming strong swelling anisotropy, while maintaining mechanical properties at elevated temperatures. While the ionic conductivity of NC700 is approximately 10% lower than that of NR211, the reduced thickness yields a 40% decrease in calculated area-specific resistance, suggesting the composite architecture can favorably shift the conductivity-stability trade-off. The composite structure also reduces gas permeability, indicating potential for improved separator function alongside favorable transport properties. Systematic deconvolution of reinforcement and additive contributions shows that conductivity losses from cerium incorporation are largely offset by gains from the lower equivalent-weight polymer, providing quantitative relationships that may guide composite membrane design for fuel cells and other electrochemical applications.

25 ENERGY STORAGE↗

NASA NTRS2025

ISS Radiator Face Sheet Anomaly Investigation and Return to Function

Maintaining sufficient heat rejection on the International Space Station (ISS) is critical to the function of the Low-Earth Orbit station. The External Active Thermal Control System (EATCS) rejects the excess heat generated by the US On-Orbit Segment (USOS) modules. The system uses single-phase liquid ammonia to collect the heat and reject it to six radiators (three on the Starboard side, three on the Port side). Each radiator is made up of eight panels. In September 2008, imagery of the Starboard radiators was conducted and showed one of the panels’ face sheets had peeled up from the internal honeycomb core structure. Out of an abundance of caution, the radiator was isolated from the rest of the EATCS and the ammonia was vented to space. Since the radiator was vented, periodic imagery of all radiator panels was taken to monitor any changes in the face sheet and inspect for other radiator anomalies. In January 2025, after years of trending and inspections, the radiator was reintegrated into the system to increase heat rejection capabilities. This paper will document the multi-year investigation that took place to determine root cause and mitigations implemented to reduce risk to the system. A summary of the data since the investigation will show rationale for reintegrating the radiator even with the damaged panel. The goal of the paper is to document the investigation for historical purposes and provide future programs with findings discovered during the investigation.

Aaron Rodriguez↗

NASA NTRS2026

ISS Radiator Face Sheet Anomaly Investigation and Return to Function

Aaron Rodriguez↗

DOE OSTI2026

A Critical Review of Heliostat Design for Concentrating Solar Thermal Technologies

Heliostat designs have undergone a widespread and eclectic development process, with many unique designs demonstrated. However, recent developments do not show the overall cohort converging toward a globally accepted universal design. Here, this study characterizes heliostats by breaking down and evaluating design traits based on emergent patterns from a comprehensive compilation of known heliostat designs spanning several decades. Four main categories for evaluation emerged: heliostat base, heliostat primary axis, heliostat drive, and facet support. Each of these four categories is further defined by four subtypes so that all heliostats fall into a single subtype within each category. The classified heliostats are ranked, yielding several view slices into the heliostat compilation or a breakdown of heliostats by type. An analysis of the breakdown shows several trends: a scale-up of established designs, new approaches at small to medium scales, and a movement toward greater adoption of linear drives. These trends reflect the most meaningful contributing factors to a proposed trajectory for a new era of heliostat designs striving to meet widely considered cost targets of $\$$50/m 2 or $\$$75/m 2 .

14 SOLAR ENERGY↗

DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗