Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuse interface methods”

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

Deep learning of interface structures from simulated 4D STEM data: cation intermixing vs. roughening ∗

Abstract Interface structures in complex oxides remain an active area of condensed matter physics research, largely enabled by recent advances in scanning transmission electron microscopy (STEM). Yet the nature of the STEM contrast in which the structure is projected along the given direction precludes separation of possible structural models. Here, we utilize deep convolutional neural networks (DCNN) trained on simulated 4D STEM datasets to predict structural descriptors of interfaces. We focus on the widely studied interface between LaAlO 3 and SrTiO 3 , using dynamical diffraction theory and leveraging high performance computing to simulate thousands of possible 4D STEM datasets to train the DCNN to learn properties of the underlying structures on which the simulations are based. We test the DCNN on simulated data and show that it is possible (with >95% accuracy) to identify a physically rough from a chemically diffuse interface and create a DCNN regression model to predict step positions. We quantify the applicability of the model to different thicknesses and the transferability of the approach. The method shown here is general and can be applied for any inverse imaging problem where forward models are present.

42 ENGINEERING↗
DOE OSTI2026

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2022

Coupling TRACE with nodal neutronics code Ants using the Exterior Communications Interface and VTT's multi-physics driver Cerberus

In this paper a recently developed coupling between system code TRACE and nodal neutronics code Ants is presented. The TRAC/RELAP Advanced Computational Engine (TRACE) is a reactor system code developed by U.S. Nuclear Regulatory Commission (NRC). The code has been designed to analyse loss-of-coolant accidents (LOCAs), operational transients and other accident scenarios in light water reactors. Ants is the reduced order nodal neutronics code of the Kraken framework. The diffusion solution method in Ants is based on the Analytic Function Expansion Nodal method (AFEN) and Flux Expansion Nodal Method (FENM). Rectangular, hexagonal and triangular geometries are supported. Both steady state and burnup problems including micro-depletion can be simulated. The coupling is implemented using VTT's multi-physics driver Cerberus and the Exterior Communications Interface (ECI) which is supplied with TRACE and enables the coupling of TRACE without source code modifications. The main focus in the paper is on describing the TRACE-Ants coupling in detail in order to share to other code developers what we have learned from the ECI based coupling process. The coupled code system is briefly tested by simulating a simple transient in a SMR core initiated by control rod movement. The results of the test calculation look reasonable and more complicated test cases will be modelled in the near future.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Multidimensional visualization of the dynamic evolution of Li metal via in situ/operando methods

The growing demands for high-energy density electrical energy storage devices stimulate the coupling of conversion-type cathodes and lithium (Li) metal anodes. While promising, the use of these “Li-free” cathodes brings new challenges to the Li anode interface, as Li needs to be dissolved first during cell operation. In this study, we have achieved a direct visualization and comprehensive analysis of the dynamic evolution of the Li interface. The critical metrics of the interfacial resistance, Li growth, and solid electrolyte interface (SEI) distribution during the initial dissolution/deposition processes were systematically investigated by employing multidimensional analysis methods. They include three-electrode impedance tests, in situ atomic force microscopy, scanning electrochemical microscopy, and cryogenic scanning transmission electron microscopy. The high-resolution imaging and real-time observations show that a loose, diffuse, and unevenly distributed SEI is formed during the initial dissolution process. This leads to the dramatically fast growth of Li during the subsequent deposition, deviating from Fick’s law, which exacerbates the interfacial impedance. The compactness of the interfacial structure and enrichment of electrolyte species at the surface during the initial deposition play critical roles in the long-term stability of Li anodes, as revealed by operando confocal Raman spectroscopic mapping. In conclusion, our observations relate to ion transfer, morphological and structural evolution, and Li (de)solvation at Li interfaces, revealing the underlying pathways influenced by the initial dissolution process, which promotes a reconsideration of anode investigations and effective protection strategies.

25 ENERGY STORAGE↗
DOE OSTI2021

Hydrogen Permeation Analysis Tool

There is a rising need for the measurement and understanding of hydrogen transport properties in materials for the energy and transportation industries. The Hydrogen Permeation Analysis Tool (HyPAT) provides a user friendly interface for batch analysis of dynamic hydrogen permeation data to calculate the transport properties permeability, solubility, and diffusivity. A built in literature database provides an easy comparison to known values. A common method for measuring these hydrogen transport properties is dynamic pressure rise permeation experiments, for example, those done in the Static Gas Absorption and Permeation (SGAP) experiment at INL. Most laboratories develop individual, in-house codes for data analysis. HyPAT provides a platform that can be adopted by other user groups for such experimental analysis and significantly reduce personnel time spent performing data analysis.

Fuerst, ThomasF↗
DOE OSTI2023

A novel design of transitional layer structure between reduced activation ferritic martensitic steels and tungsten for plasma facing materials

Plasma-facing components (PFCs) are among the most critical gaps for fusion energy to establish technical and economic feasibility. Tungsten as a first wall/blanket material in PFCs requires to be integrating with reduced activation ferritic martensitic (RAFM) steels as a structural component. Currently, major drawbacks are the requirement of brazing, the formation of a brittle interface, and a large difference between the coefficients of thermal expansion of tungsten and steel. Here, a novel transitional multilayer structure was designed and investigated to join tungsten and RAFM steels using three interlayers. The composition of each interlayer was selected based on computational thermodynamics and diffusion kinetics to ensure a body-centered cubic (bcc) single-phase structure and prevent the formation of a brittle intermetallic phase region in the temperature range of 600–1150 °C. Although the transitional layer structure was designed for additive manufacturing, spark plasma sintering (SPS) as proof of concept was used to bond the individual layers. Interfaces were investigated using scanning and transmission electron microscopy methods but no layered intermetallic phase was observed. Nanoindentation maps across the interface suggest major hardness differences at the interface between tungsten and the vanadium interlayer, as well as the interface between RAFM steel and the FeCrAl interlayer.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2025

An end-to-end deep learning method for solving nonlocal Allen–Cahn and Cahn–Hilliard phase-field models

Here, we propose an efficient end-to-end deep learning method for solving nonlocal Allen–Cahn (AC) and Cahn–Hilliard (CH) phase-field models. One motivation for this effort emanates from the fact that discretized partial differential equation-based AC or CH phase-field models result in diffuse interfaces between phases, with the only recourse for remediation is to severely refine the spatial grids in the vicinity of the true moving sharp interface whose width is determined by a grid-independent parameter that is substantially larger than the local grid size. In this work, we introduce non-mass conserving nonlocal AC or CH phase-field models with regular, logarithmic, or obstacle double-well potentials. Because of non-locality, some of these models feature totally sharp interfaces separating phases. The discretization of such models can lead to a transition between phases whose width is only a single grid cell wide. Another motivation is to use deep learning approaches to ameliorate the otherwise high cost of solving discretized nonlocal phase-field models. To this end, loss functions of the customized neural networks are defined using the residual of the fully discrete approximations of the AC or CH models, which results from applying a Fourier collocation method and a temporal semi-implicit approximation. To address the long-range interactions in the models, we tailor the architecture of the neural network by incorporating a nonlocal kernel as an input channel to the neural network model. We then provide the results of extensive computational experiments to illustrate the accuracy, predictive capabilities, and cost reductions of the proposed method.

42 ENGINEERING↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2022

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li 7 La 3 Zr 2 O 12 solid electrolytes

Batteries based on solid-state electrolytes, including Li 7 La 3 Zr 2 O 12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. Here these examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

25 ENERGY STORAGE↗
DOE OSTI2025

Multimaterial 3D Printing in Activating Bath Enables In Situ Polymerization of Thermosets with Intricate Geometries and Diverse Elastic Behaviors

Polydicyclopentadiene, p(DCPD), is a high‐performance thermoset valued for its exceptional toughness, strength, and stiffness. When copolymerized with 1,5‐cyclooctadiene (COD), its mechanical properties can be tuned from glassy to rubbery at room temperature. While frontal polymerization enables a rapid and energy‐efficient route to 3D print DCPD‐based materials, challenges such as ink shelf life and gravitational distortion, especially in direct ink writing of soft COD‐rich formulations, must be considered. Here, a complementary chemical strategy is presented, embedded 3D printing, that enables localized in situ polymerization of printed DCPD/COD inks within a reactive support matrix. The matrix provides both physical support and a reservoir of chemical activator, which diffuses into the ink, activates a latent bis(N‐heterocyclic carbene) Ru precatalyst, and initiates ring‐opening metathesis polymerization. Curing begins at the ink–matrix interface and propagates inward via diffusion, stabilizing the interface and preventing capillary‐driven deformation regardless of the matrix yield stress. This approach eliminates the need for cold storage, external curing, or photoinitiation, significantly expanding the processing window. Using this method, diverse thermosetting and elastomeric architectures are fabricated with features as small as 5 µm and aspect ratios of 100, including interlinked chains, shallow spherical shells exhibiting snap‐through buckling, and hair‐like fin arrays inaccessible through traditional techniques.

chemical activation↗
DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermally Stable Ruthenium Contact for Robust p‑Type Tellurium Transistors

Tellurium (Te) is attractive for p-channel transistors due to its high hole mobility. Despite having a low thermal budget suitable for back-end-of-line (BEOL) monolithic integration, the practical realization of Te transistors is hindered by its thermal stability. In this work, we investigate thermal stability for Te thin films grown via scalable thermal evaporation. Our findings identify ruthenium as a more thermally stable contact for p-type Te transistors, capable of withstanding temperatures up to 250 °C. Ruthenium exhibits significantly lower diffusivity in Te compared to other contact metals commonly used such as nickel and palladium. Using the transfer-length method, we measured a contact resistance of 1.25 kΩ·μm at the ruthenium-tellurium interface. Additionally, the incorporation of high-κ ZrO2 encapsulation not only suppresses the sublimation of the Te channel at elevated temperatures but also serves as the gate dielectric in top-gate devices operating at 1 V, achieving an on/off current ratio of 105.

Rahman, I K M Reaz↗
DOE OSTI2023

Negligible magnetic losses at low temperatures in liquid phase epitaxy grown Y 3 Fe 5 O 12 films

Yttrium iron garnet (Y 3 Fe 5 O 12 ; YIG) has a unique combination of low magnetic damping, high spin-wave conductivity, and insulating properties that make it a highly attractive material for a variety of applications in the fields of magnetics and spintronics. While the room-temperature magnetization dynamics of YIG have been extensively studied, there are limited reports correlating the low-temperature magnetization dynamics to the material structure or growth method. Here, in this study, we investigate liquid phase epitaxy grown YIG films and their magnetization dynamics at temperatures down to 10 K. We show there is a negligible increase in the ferromagnetic resonance linewidth down to 10 K, which is unique when compared with YIG films grown by other deposition methods. From the broadband ferromagnetic resonance measurements, polarized neutron reflectivity, and scanning transmission electron microscopy, we conclude that these liquid phase epitaxy grown films have negligible rare-earth impurities present, specifically the suppression of Gd diffusion from the Gd 3 Ga 5 O 12 (GGG) substrate into the Y 3 Fe 5 O 12 film, and therefore negligible magnetic losses attributed to the slow-relaxation mechanism. Overall, liquid phase epitaxy YIG films have a YIG/GGG interface that is five times sharper and have ten times lower ferromagnetic resonance linewidths below 50 K than comparable YIG films by other deposition methods. Thus, liquid phase epitaxy grown YIG films are ideal for low-temperature experiments/applications that require low magnetic losses, such as quantum transduction and manipulation via magnon coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electrical Analysis of Pulsed Laser Annealed Poly-Si:Ga/SiOx Passivating Contacts

High-efficiency single crystalline silicon (c-Si) solar cells require precise control of dopant diffusion profiles and highly active doping concentrations through thermal annealing. Conventional furnace annealing has been successfully employed for dopant diffusion in aluminum back-surface field (Al-BSF), passivated emitter and rear contact (PERC), and Topcon cells. However, furnace annealing limits some next-generation polycrystalline silicon (poly-Si) on SiOx passivating contacts. This limitation largely affects p-type passivating contacts by not being able to: 1) control the dopant diffusion accurately to prevent B segregation at the SiOx and c-Si interface; and 2) provide highly activated doping concentrations for low contact resistivity. In this contribution, we use a nanosecond excimer laser to examine the passivation quality and electrical performance of both B- and Ga-doped poly-Si/SiOx passivating contacts. The core of this method is to take advantage of the nonequilibrium nature of the anneal through rapid melting and recrystallizing the poly-Si in a short timescale and to achieve doping concentrations above the solid solubility limit. Simulations were performed on polished surfaces to visualize the doping diffusion profiles under different laser conditions, and secondary ion mass spectrometry (SIMS) was used to verify the experimental diffusion profiles post pulsed laser melting (PLM). The results show the dopant diffusion profiles can be tuned precisely through PLM. The electrical analysis using VdP-Hall measurements reveals that the active doping concentration reached 10^21 cm-3 for B and ~2×10^21 cm-3 for Ga in poly-Si, far exceeding their solid solubility limit in Si, with a nearly 100% dopant activation achieved for B, and 20% for Ga. This highly activated dopant profile results in a low contact resistivity of <15 mO·cm^2 for B-doped contacts and <40 mO·cm^2 for Ga-doped contacts. We demonstrate that the passivation quality of the laser annealed contacts shows a linear dependency with the laser energy density and the number of pulses. PLM was also performed on random pyramid textured samples, and we found that preferential melting exists for the tips of pyramids, which can be pursued further for selective pinhole opening. Additionally, we perform thermal stability tests of these hyperdoped passivating contacts and show that under certain laser conditions, the B-doped samples exhibit good thermal stability with near full activation for temperatures up to 800°C. Lastly, we transfer this technique to a larger area by overlapping laser spots and show uniform passivation quality and electrical performance. We are currently implementing this PLM annealed p-type passivating contacts into double-side back-junction poly-Si/SiOx passivating contact devices. Our results show that PLM can not only enable the next generation of high-efficiency c-Si photovoltaics technology but can also benefit the c-Si integrated circuit industry through hyperdoping other semiconductors with atoms that may otherwise have a low solubility.

doping diffusion profile↗
DOE OSTI2023

Chemical and Electronic Structure at the Interface between a Sputter-Deposited Zn(O,S) Buffer and a Cu(In,Ga)(S,Se) 2 Solar Cell Absorber

The chemical and electronic structure of the interface between a sputter-deposited Zn(O,S) buffer layer and an industrial Cu(In,Ga)(S,Se) 2 (CIGSSe) absorber for thin-film solar cells is investigated with X-ray and UV photoelectron spectroscopy, inverse photoemission spectroscopy, and X-ray emission spectroscopy. We find a CIGSSe absorber surface band gap of 1.61 (±0.14) eV, which is significantly increased as compared to the minimal value derived with bulk-sensitive methods (≈1.1 eV). We find no indication for diffusion of absorber elements into the buffer layer. Surface- and bulk-sensitive measurements of the buffer layer suggest the presence of S-Zn and S-O bonds in the Zn(O,S) layer. We find that the naturally existing downward band bending toward the CIGSSe absorber surface is increased by the formation of the interface, likely enhancing carrier separation under illumination. We also derive a flat conduction band alignment, in line with the reported high conversion efficiencies of corresponding large-area solar cells.

14 SOLAR ENERGY↗
DOE OSTI2025

In Situ Atomic Tracking on the Interfacial Etching and Reconfiguration of Cu-ReSe 2 Contact during Thermal Annealing

The Schottky barrier height can be greatly affected by the metal diffusion, reaction, and covalent bonding formation at the contact. Exploring novel methods and revealing the fundamental mechanisms for contact engineering are of vital importance for microelectronic devices. Here, in this study, the annealing induced interfacial reactions at Cu-ReSe 2 contact are dynamically revealed from the atomic scale. Accompanied by the diffusion of Se to Cu, ReSe 2 is gradually decomposed to a thin Re interlayer through a “chain-by-chain” manner. Theoretical calculations show that the Cu atoms can facilitate the chemical bond breaking of ReSe 2 , significantly lowering the Se diffusion energy barrier toward Cu. The formed Re/ReSe 2 heterostructure presents a metal-like band structure, which underscores the critical role of Cu in altering the interfacial chemistry and promoting carrier transport across the interface. Our results can provide vital insights into the contact properties of ReSe 2 and provide a possible method for fabricating high-performance ReSe 2 -based devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗