Search NASA⌕ Search

SEARCH · Search NASA

Results for “near surface”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Independent Review of Groundwater Remediation Strategy for Hexavalent Chromium and RDX Groundwater Plumes at Los Alamos National Laboratory (Rev. 1)

Site operations at the Los Alamos National Laboratory (LANL) resulted in the release of oxidized chromium, Cr(VI), into Sandia Canyon from cooling tower effluent from 1956 until 1972. The chromium traveled with the surface water approximately 3 miles downstream before migrating below ground surface. Chromium concentrations exceed 50 μg/L in the upper portion of the aquifer. Another LANL groundwater plume of concern is associated with RDX (Royal Demolition Explosives, 1,3,5-trinitro-1,3,5-triazine). Between 1951 and 1996, RDX was released to the mesa-top facilities' process water outfall, adjacent and underlying soils, and alluvial sediments, along with surface water in Cañon de Valle. Between 2000 and 2010, two remedial actions were deployed, removing much of the near-surface RDX, however, recharge due to precipitation has transported RDX into the perched-intermediate zone and into the regional aquifer. The report documents an independent technical review by scientists from the Department of Energy (DOE) Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) to provide recommendations for potential near term actions to address and optimize remediation for both the Cr(VI) and RDX plumes. The proposed near-term remedial actions include design of pump and treat systems for Cr(VI) and monitoring and study for natural attenuation for RDX. The review assesses existing data, conceptual and numerical modeling, and it recommends a technical integration process to support identifying and implementing strategic, effective and efficient remedies. The DOE Environmental Management Los Alamos Field Office (EM-LA) and their cleanup contractor Newport News Nuclear-BWTX, LLC Los Alamos (N3B) provided the information required for the review. Interviews were also conducted with regulators to obtain the full spectrum of technical, regulatory and scientific perspectives. The independent review team was impressed by the capabilities, experiences, innovativeness, and insightfulness of the technical representatives from both the regulator, the New Mexico Environment Department (NMED) and N3B. Incorporation of vadose zone flow pathways in the conceptual site model (CSM) and configuring the numerical modeling for the site was generally state-of-the-practice (or better). This could be considered state-of-the-art by addressing uncertainties related to spatial extent of hydraulic windows. The reviews from the regulators were thorough and often provided useful concepts for consideration and future/study resolutions. The overarching consensus recommendation of independent review team is that the LANL groundwater plumes should be addressed in context of the emerging "management of complex sites" paradigm.

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

Probing thermal conductivity of subsurface, amorphous layers in irradiated diamond

In this study, we report on the thermal conductivity of amorphous carbon generated in diamond via nitrogen ion implantation (N 3+ at 16.5 MeV). Here, transmission electron microscopy techniques demonstrate amorphous band formation about the longitudinal projected range, localized approximately 7 μm beneath the sample surface. While high-frequency time-domain thermoreflectance measurements provide insight into the thermal properties of the near-surface preceding the longitudinal projected range depth, a complimentary technique, steady-state thermoreflectance, is used to probe the thermal conductivity at depths which could not otherwise be resolved. Through measurements with a series of focusing objective lenses for the laser spot size, we find the thermal conductivity of the amorphous region to be approximately 1.4 W m –1 K –1 , which is comparable to that measured for amorphous carbon films fabricated through other techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Helium bubble size effects on the surface morphological response of plasma-facing tungsten

Abstract We report a simulation study on the effects of helium (He) bubble size on the morphological evolution and pattern formation on the surface of tungsten used as a plasma-facing component (PFC) in nuclear fusion devices. We have carried out a systematic investigation based on self-consistent dynamical simulations of surface morphological evolution according to an atomistically-informed, 3D continuum-scale model that captures well the relevant length and time scales of surface nanostructure formation in PFC tungsten. The model accounts for PFC surface diffusion, driven by the biaxial compressive stress originating from the over-pressurized He bubbles in the near-surface region of PFC tungsten as a result of He plasma exposure, combined with the formation of self-interstitial atoms in tungsten that diffuse toward the PFC surface and the flux of surface adatoms generated as a result of surface vacancy-adatom pair formation upon He implantation; this transport of surface adatoms contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. Our detailed analysis reveals that varying the average He bubble size plays an important role in the PFC surface growth kinetics as well as the resulting surface topography. Specifically, we find that an increase in the He bubble size leads to a deceleration in the growth rate of the tungsten nanotendrils that emanate from the PFC surface. We also find that the separation distance between the resulting surface features increases with increasing He bubble size, as well as over time. This coarsening effect is a thermally activated process resulting in an accurate description of the temperature dependence of the average surface feature separation by an Arrhenius relation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Seismoacoustic Analysis of the Large Surface Explosion Coupling Experiment Using a Large- N Seismic Array

Seismoacoustic signals at local distance (<~10 km) are widely used as important constraints on source parameters for near-surface events, yet the seismoacoustic wave generation and energy partitioning are not fully understood. Spatially dense sensors could provide observations in high resolution to capture the full wavefield for better understanding wave propagation and improving source estimation. Recently, spatially dense observations of the local seismoacoustic wavefield produced by a pair of 1-ton surface explosions have been recorded using a large-N seismic array. This large-N array consists of 446 geophones and covers an area of approximately 2×2.5 km 2 . The two surface explosions occur at the same location but at different times with different atmospheric conditions. Both seismic and air–ground coupled acoustic waves from the two surface explosions are well observed. Analyses of signals recorded by the large-N seismic array show different acoustic wave speed and amplitude for the two explosions. A strong spatial variability in acoustic wave speed and amplitude for each explosion is also observed. Finally, the observations suggest the important role of local atmosphere state on wave propagation and source estimation and demonstrate how the use of a large-N capability can improve characterization of the propagation medium and source.

58 GEOSCIENCES↗
DOE OSTI2021

Influence of Ligand Structure on Excited State Surface Chemistry of Lead Sulfide Quantum Dots

The ligand-nanocrystal boundaries of colloidal quantum dots (QDs) mediate the primary energy and electron transfer processes that underpin photochemical and photocatalytic transformations at their surfaces. We use mid-infrared transient absorption spectroscopy to reveal the influence that ligand structure and bonding to nanocrystal surfaces have on the changes of the excited state surface chemistry of this boundary in PbS QDs and the corresponding impact on charge transfer processes between nanocrystals. We demonstrate that oleate ligands undergo marked changes in their bonding to surfaces in the excitonic excited states of the nanocrystals, indicating that oleate passivated PbS surfaces undergo significant structural changes following photoexcitation. These changes can impact the surface mobility of the ligands and the ability of redox shuttles to approach the nanocrystal surfaces to undergo charge transfer in photocatalytic reactions. In contrast, markedly different transient vibrational features are observed in iodide/mercaptoproprionic acid passivated PbS QD films that result from charge transfer between neighboring nanocrystals and localization of holes at the nanocrystal surfaces near MPA ligands. This ability to distinguish the influence that excitonic excited states vs charge transfer processes have on the surface chemistry of the ligand-nanocrystal boundary lays the groundwork for exploration of how this boundary can be understood and controlled for the design of nanocrystalline materials tailored for specific applications in solar energy harvesting and photocatalytic reactions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗
DOE OSTI2026

Ionic Precursors Transformed Into Vinyl Acetate Synthesis Catalyst via Reaction-Driven Restructuring

Conventional preparation of supported bimetallic catalysts relies on solution-mediated metal salt immobilization and pre-formation of alloy nanoparticles before reaction. Here, we report a fundamentally different synthesis strategy of using a physical mixture of salt precursors to generate an active catalyst during reaction. The catalytic structure is generated in situ from Pd3(OAc)6, Au(OH)3, and KOAc through H2 treatment and reaction-driven restructuring under vinyl acetate monomer (VAM) synthesis conditions. Ascertained from in situ X-ray diffraction and operando infrared spectroscopy analyses, reduction treatment produces segregated Pd and Au domains, and subsequent exposure to a VAM reaction mixture triggers dynamic extraction of Pd from the metal surface. This latter process, mediated by acetate-assisted redox cycles, facilitates Pd migration toward Au domains to form a near-surface localized Pd50Au50 alloy phase. Monometallic Pd domains serve as a reservoir of Pd to the alloy phase, leading to and sustaining a more Pd-enriched active surface and a higher population of accessible Pd sites, compared to a conventionally prepared K-PdAu/SiO2 catalyst. Consequently, this leads to a twofold increase in the VAM formation rate, demonstrating highly active bimetallic catalysts can be generated through the gas-phase treatment of physically mixed ionic precursors.

Cha, Byeong Jun [Rice University]↗
DOE OSTI2020

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

SEM/EBSD Characterization of Baffle-Former Bolts Designated for High Energy X-ray Microscopy Studies at the Advanced Photon Source

This report documents the microstructural and surface characterization of two high-fluence baffle-former bolts (BFBs) harvested from a commercial pressurized water reactor (PWR) after nearly four decades of service. These bolts experienced extreme conditions including high neutron fluence, elevated temperatures, corrosive primary water, and sustained mechanical loads, all of which contribute to degradation mechanisms such as irradiation-assisted stress corrosion cracking (IASCC). Using scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD), the study assessed bulk and near-surface microstructures, grain boundary characteristics, and in-service oxide layers. The results revealed significant differences between bolts exposed to 41 dpa and 29 dpa, including differences in recrystallization extent, deformation structures, and oxidation morphology. These findings establish critical baseline data for integration with advanced synchrotron-based high-energy X-ray microscopy (HEXM) and computational crystal plasticity models. Ultimately, the work supports improved mechanistic understanding and predictive capability regarding the long-term performance of reactor core internal components.

Lach, Tim [Oak Ridge National Laboratory (ORNL), O↗
DOE OSTI2020

Oxidation State and Surface Reconstruction of Cu under CO 2 Reduction Conditions from In Situ X-ray Characterization

The electrochemical CO 2 reduction reaction (CO 2 RR) using Cu-based catalysts holds great potential for producing valuable multi-carbon products from renewable energy. However, the chemical and structural state of Cu catalyst surfaces during the CO 2 RR remains a matter of debate. In this paper, we show the structural evolution of the near-surface region of polycrystalline Cu electrodes under in situ conditions through a combination of grazing incidence X-ray absorption spectroscopy (GIXAS) and X-ray diffraction (GIXRD). The in situ GIXAS reveals that the surface oxide layer is fully reduced to metallic Cu before the onset potential for CO 2 RR, and the catalyst maintains the metallic state across the potentials relevant to the CO 2 RR. We also find a preferential surface reconstruction of the polycrystalline Cu surface toward (100) facets in the presence of CO 2 . Quantitative analysis of the reconstruction profiles reveals that the degree of reconstruction increases with increasingly negative applied potentials, and it persists when the applied potential returns to more positive values. These findings show that the surface of Cu electrocatalysts is dynamic during the CO 2 RR, and emphasize the importance of in situ characterization to understand the surface structure and its role in electrocatalysis.

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

Surface radiation trends at North Slope of Alaska influenced by large-scale circulation and atmospheric rivers

Arctic amplification manifests as a pervasive warming trend emerging over the past century in near-surface air temperature throughout the Arctic that is double the globally averaged temperature increase throughout most of the year. It results from complex processes involving oceanic, atmospheric and terrestrial components which require detailed study to discern roles of the fundamental processes involved to improve predictions of the Arctic environment. We report on signals that are beginning to emerge, on a timescale predicted by recent satellite remote sensing studies, from the unique 25 year record of detailed surface-based radiometer measurements obtained by the US Department of Energy Atmospheric Radiation Measurement (ARM) Facility North Slope of Alaska (NSA) site at Utqiaġvik, Alaska. Statistically significant warming trends are found at the site in the boreal fall, while a decrease in net radiation occurs in late summer. This decrease is driven primarily by the decrease in shortwave radiation resulting from increasing cloud liquid water path as observed by the microwave radiometer. The magnitude of the liquid water path trends is slightly larger than the root mean square error (RMSE) due to natural variability, while the magnitude of the net radiation trends is smaller than the RMSE natural variability, emphasizing that these signals are only just emerging. Analysis of prevailing meteorological regimes linking NSA with the Arctic Ocean and subarctic latitudes, and atmospheric rivers, suggests that specific changing circulation patterns are the primary driver for these summertime trends.

Lubin, Dan [Univ. of California, San Diego, CA (Un↗
DOE OSTI2025

Mechanistic Insights Into Fatigue Life Enhancement of High‐Strength Steel via Ultrasonic Impact Treatment

In this work we explore the impact of Ultrasonic Impact Treatment (UIT) on the fatigue performance of high‐strength microalloyed steel commonly used in crankshaft applications. Building on prior observations of fatigue life improvement with UIT, this work focuses on unraveling the process–structure–performance relationships underpinning these enhancements. Microstructural analysis revealed significant grain refinement in the near‐surface layers, with deformation depth increasing at higher impact energies. This led to increased surface hardness and the development of deeper compressive residual stresses, particularly in samples treated with higher impact energy. Rotating bending fatigue testing showed a substantial improvement in fatigue life for UIT‐treated samples, with the endurance limit nearly doubling compared with untreated specimens. Fractographic analysis revealed a transition in crack initiation from surface defects in untreated samples to interior regions in UIT‐treated samples, characterized by the formation of noninclusion‐induced granular bright facets (GBFs). The observed fatigue enhancement is attributed to the synergistic effects of strain hardening and compressive residual stresses, which increase the surface crack threshold and promote interior crack initiation. This study provides new mechanistic insight into UIT‐induced fatigue resistance and interior failure behavior in high‐strength steels.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulating Alpha Particles Incident on MKID Chips for Quantum Sensitivity Analysis

Superconducting quantum devices, such as microwave kinetic inductance detectors (MKIDs), are highly sensitive instruments used in quantum computing and advanced sensing technologies. However, their extreme sensitivity also makes them vulnerable to background noise from natural sources like radiation. One significant contributor to this noise is alpha particles emitted by 210Po, a radon decay daughter that accumulates on surfaces near the detector. This project investigates how alpha particles emitted from 210Po interact with MKID chips. These particles can deposit energy on the detector surface, disrupting its operation and generating false signals. Understanding the energy and behavior of these particles is crucial for improving the design and reliability of quantum devices. To explore this, we first modeled the decay chain starting from 210Pb to 210Po using differential equations. This allowed us to predict how the activity of alpha-emitting isotopes changes over time, reaching a steady state after about two years. Next, we simulated alpha particle interactions with the MKID chip using the Geant4 software toolkit. We built a detailed computer model of the detector housing, including the copper lid where alpha particles originate, the silicon chip, and a thin aluminum sensor layer. Alpha particles were emitted isotropically from just beneath the copper lid’s surface, mimicking natural decay conditions. The simulation tracked how these particles deposit energy on the chip, generating electron-hole pairs and phonons. The results provide insight into the behavior of the resultant electron-hole pairs and phonons, giving us a clear understanding of the energy deposition distribution on the chip. This work supports efforts to mitigate background noise in superconducting sensors, advancing their use in quantum computing and sensitive physics experiments.

Hall, Matthew [Fermilab; UCLA]↗
DOE OSTI2022

Energetic particle optimization of quasi-axisymmetric stellarator equilibria

Abstract An important goal of stellarator optimization is to achieve good confinement of energetic particles such as, in the case of a reactor, alphas created by deuterium–tritium fusion. In this work, a fixed-boundary stellarator equilibrium was re-optimized for energetic particle confinement via a two-step process: first, by minimizing deviations from quasi-axisymmetry (QA) on a single flux surface near the mid-radius, and secondly by maintaining this improved QA while minimizing the analytical quantity Γ C , which represents the angle between magnetic flux surfaces and contours of J | | , the second adiabatic invariant. This was performed multiple times, resulting in a group of equilibria with significantly reduced energetic particle losses, as evaluated by Monte Carlo simulations of alpha particles in scaled-up versions of the equilibria. This is the first time that energetic particle losses in a QA stellarator have successfully been reduced by optimizing Γ C . The relationship between energetic particle losses and metrics such as QA error ( E q a ) and Γ C in this set of equilibria were examined via statistical methods and a nearly linear relationship between volume-averaged Γ C and prompt particle losses was found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Operando Unveiling of Hydrogen Spillover Mechanisms on Tungsten Oxide Surfaces

Hydrogen spillover is an important process in catalytic hydrogenation reactions, facilitating H 2 activation and modulating surface chemistry of reducible oxide catalysts. This study focuses on the operando unveiling of platinum-induced hydrogen spillover on monoclinic tungsten trioxide (γ-WO 3 ), employing ambient pressure X-ray photoelectron spectroscopy, density functional theory calculations and microkinetic modeling to investigate the dynamic evolution of surface states at varied temperatures. At room temperature, hydrogen spillover results in the formation of W 5+ and hydrogen intermediates (hydroxyl species and adsorbed water), facilitated by Pt metal clusters. With increasing temperature, water desorption, reverse hydrogen spillover and surface-to-bulk diffusion of hydrogen atoms compete with each other, leading initially to reoxidation and then further reduction of W atoms in the near-surface. The combined experimental results and simulations provide a comprehensive understanding of the mechanisms underlying hydrogen interaction with reducible metal oxides, lending insights of relevance to the design of enhanced hydrogenation catalysts.

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

Correlated metals and unconventional superconductivity in rhombohedral trilayer graphene: A renormalization group analysis

Here, motivated by recent experimental observations of correlated metallic phases and superconductivity in rhombohedral trilayer graphene (RTG), we perform an unbiased study of electronic ordering instabilities in hole-doped RTG. Specifically, we focus on electronic states energetically proximate to Van Hove singularities (VHSs), where a large density of states promotes different interaction-induced symmetry-breaking electronic orders. To resolve the Fermi surface near VHSs, we construct a fermionic hot-spot model and demonstrate that a perpendicular electric field can tune different nesting structures of the Fermi surface. Subsequently, we apply a renormalization group analysis to describe the low-energy phase diagrams of our model under both short-range repulsive interactions as well as realistic (long-range) Coulomb interactions. Our analysis shows instabilities towards either intervalley coherent metallic phases or superconducting phases. The dominant pairing channel depends crucially on the nature of Fermi surface nesting—repulsive Coulomb interaction favors spin-singlet d-wave pairing for relatively small displacement field and spin-singlet i-wave pairing for larger displacement field. We argue that the phase diagram of RTG can be well-understood by modeling the realistic Coulomb interaction as the sum of repulsive density-density interaction and ferromagnetic spin-triplet intervalley coherence (IVC) Hund's coupling, while phonon-mediated electronic interactions have a negligible effect on this system, in sharp contrast to twisted graphene multilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗