Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temperature heterogeneity”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2021

Challenges of Investigating Redox Reactions in Room Temperature Ionic Liquids [Slides]

The goal of this project is to understand the effects of the morphology and chemical properties of ionic liquids (IL) on the heterogeneous electron transfer processes of simple inorganic ions. This project will seek to understand the fundamental chemical and physical phenomena in ionic liquids through electrochemical, and material characterization techniques.

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

Development and Application of High-Fidelity Models for Heterogeneous CO2 Frost Formation

Carbon America has developed a cryogenic carbon capture technology ("FrostCC") that separates CO2 from point source emissions by solidifying it at cold temperatures through preferential desublimation. Cooling is achieved through a series of interlinked compression, heat exchange, and expansion operations. In the current system, frosting of CO2 happens in heat exchangers, followed by CO2 recovery in a separate extraction step. In this work, multiphysics computational fluid dynamics (CFD) models are developed and validated for compressible and low Mach flows to simulate the formation of solid CO2 in flue gas flowing in a heat exchanger geometry. The models track the mass transfer rate of CO2 from gas phase to solid phase, heat released from desublimation, and the evolution of the solid CO2 layer. Simulations are used to answer scientific questions related to the angle of heat exchanger pipes, where buoyancy effects from flow velocity and pipe orientation influence CO2 frosting. Results show that upwardly angled pipes produce notably different flow structures compared to horizontal or vertical configurations, and that carbon capture efficiency correlates with buoyancy effects for pipe angles within plus or minus 23 degrees of horizontal.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A First Principles Approach to Spectral Phonon Transport in Heterostructures

Understanding thermal transport across interfaces which give rise to a thermal resistance (also known as Kapitza resistance) is a critical issue affecting the development of nanotechnologies. Much modern and emergent nanotechnology consist of adjacent materials, and phonon mediated heat transfer governs thermal behavior across internal interfaces in these devices. The physics of thermal transport in solids are governed both by phenomena occurring at the atomic scale and interactions with the material's microstructure. The forecasting of fundamental quantities such as temperature, heat flux and thermal conductivity typically employs the semi-classical Boltzmann transport equation to predict the macroscopic behavior of materials in terms of the microscopic dynamics of its heat carriers. Kapitza resistance was first discovered in liquid helium experiments and has led to a fundamental research thrust in micro and nano-scale heat transport, the behavior of thermal carriers across internal interfaces. Thermal interfacial resistance (TIR) is a widely studied phenomenon, first engaged by Swartz and Pohl through their development of the acoustic and diffuse mismatch methods, then continued through myriad efforts with varying methods and approaches in an attempt to resolve carrier behavior at thermal interfaces. Many of the fundamental approaches to TIR have been at the nanoscale, and research is conducted with molecular dynamics (MD) and density functional theory (DFT) methods. The limitations of these methods is system size; atomistic methods tend to be limited to system sizes of 100,000 atoms or less. Larger length-scale methods have also been pursued, based on the principles of acoustic or diffuse mismatch, but not all include simulation of TIR using a full phonon band spectrum, or temperature dependent methods. Our approach to enabling phonon transport in layered materials draws upon our previous work of demonstrating spectrally coupled phonon transport in homogeneous and heterogeneous materials. We use a semi-analytical approach in which the Bose-Einstein (B-E) statistics set the strength of the phonon radiance in a frequency group, but the B-E statistics are informed with information from the transport system. The B-E statistics in a single frequency group feels the influence of all the groups through the spatial temperature. We also include a new field term which is an indicator of the amount of non-equilibrium behavior of the phonon spectrum---this is added to the phonon source term in all groups to ensure closure and conservation of energy, as the phonon groups in the transport system and the analytical systems are coupled. This work builds upon our previous approach by adding a phonon coupling term at an internal interface, using the principles of the DMM through transmission and reflection coefficients. In this work, the coefficients are determined through computing a common temperature at the interface, influenced by the phonon band structure of both materials, in effect, providing mixing between the two material systems and using the common temperature to set the strength of the phonon radiance at the boundaries on either side of the interface. Our approach uses material properties computed along various crystallographic orientations, and while some isotropy is built into the interface condition, the material properties weight the phonon distributions in the proper crystalline direction. Greater resolution of phonon behavior in proximity to an interface, and more accurate predictions of TIR are obtained. While it is true the assumption of diffuse mismatch can yield inconsistent results compared to experiment especially at low temperatures, this work focuses on room temperature and beyond effects, for future applications in nuclear fuel, or thermoelectric devices; a modified mismatch approach may be feasible if applied properly. Additionally, our methods focus on bridging mesoscale to engineering scale

36 MATERIALS SCIENCE↗
DOE OSTI2021

Irradiation stability and thermomechanical properties of 3D-printed SiC

Neutron irradiation tests were carried out on 3D-printed SiC derived from binderjet additive manufacturing and chemical vapor infiltration. Irradiation was carried to 2.3 dpa over a temperature range of 400–850 °C. Anisotropy that had been observed in the thermal conductivity of 3D-printed SiC prior to irradiation vanished after irradiation as the irradiation defect thermal resistivity accumulated in the material. No degradation in strength was observed in the material before or after irradiation, at various temperatures, or in different orientations. Electron microscopy of the microstructure after neutron irradiation showed distinct defect morphologies in the heterogenous material, but no evidence for irradiation-induced cracking or degradation in the microstructure was observed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Tendencies of Soil Microbial NO Emissions During HI‐SCALE as Predicted by a Nitrification/Denitrification Scheme

Many atmospheric chemical processes, including the formation of secondary organic aerosol (SOA), are strongly modulated by the reactions of NO and NO 2 (NO x ). Though NO x is controlled by anthropogenic emissions near urban areas, in rural areas soil microbes can be a significant contribution to NO emission globally. The relative rates of emissions of different nitrogen-containing species (e.g., NO, N 2 O, HONO, and N 2 ) are strong functions of soil properties such as temperature, moisture content, pH, and soil carbon and nitrogen pools. However, typical large-scale biogeochemical models either express these emissions simplistically, or not at all. Here we investigate the potential impact of soil NO emissions on atmospheric chemistry and SOA formation over regional and monthly time scales, specifically the 2016 spring and summer Intensive Observational Periods (IOPs) of the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign. We implement the soil NO nitrification/denitrification parameterization of Rasool et al. (2019) into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), supplemented by a 1-km soil moisture analysis. We then simulate both IOPs over the U.S. Great Plains, evaluating against ground stations and flight data. We show that soil NO emissions can account for a large fraction of total NO x and locally increase O 3 concentrations by up to 25%, while alleviating negative biases of gases and aerosols toward observations. Soil moisture and temperature changes between IOP1 and IOP2 lead to overall differences in emissions, but with large regional variability due to heterogeneous surface characteristics.

Atmoshpheric Chemistry↗
DOE OSTI2020

Correlation between the temperature evolution of the interfacial region and the growing dynamic cooperativity length scale

Here, we study experimentally the temperature evolution of the thickness of the interfacial layer, L int (T), between bulk matrices and the surface of nanoparticles in nanocomposites through broadband dielectric spectroscopy. Analyses revealed a power-law dependence between the logarithm of structural relaxation time in the interfacial layer, τ int (T), and the L int (T): ln τ i n t ( T ) / τ 0 ∝ L i n t β ( T ) / T , with τ 0 ~ 10 –12 s, and β index ~0.67 at high temperatures and ~1.7 at temperatures close to the glass transition temperature. In addition, our analysis revealed that the L int (T) is comparable to the length scale of dynamic heterogeneity estimated from previous nonlinear dielectric measurements and the four-point NMR [ξ NMR (T)], with L int (T) ~ ξ NMR (T). These observations may suggest a direct correlation between the L int (T) and the size of the cooperatively rearranging regions and have strong implications for understanding the dynamic heterogeneity and cooperativity in supercool liquids and their role in interfacial dynamics.

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

Nonmonotonic effect of chemical heterogeneity on interfacial crack growth at high-angle grain boundaries in Fe-Ni-Cr alloys

Here, an intermittent pattern is observed in the modeling of interfacial cyclic-loading crack growth at high-angle grain boundaries in ternary Fe-Ni-Cr alloys. Different from conventional wisdom of stress-intensity factor, the abrupt crack advances are found driven by extreme value statistics—namely, the aggregation of atoms with most compressive residual stresses. In addition, inherently non-affine atomic stress fluctuations are discovered, and the fluctuations peak at intermediate level of chemical heterogeneity, causing the fastest crack growth. Implications of such nonmonotonic mechanism in regard to the origin of intermediate-temperature embrittlement phenomena are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Optimizing DOE Opportunities to Research Land–Atmosphere Interactions in the U.S. Southeast (Workshop Report)

The southeastern United States (Southeast), with its complex and varied environments, is an area of tremendous economic, ecological, and societal importance to the country. The region is characterized by heterogeneous landscapes (i.e., geology and soil type) and a long history of human land use coupled with a warm temperature regime and high precipitation. As a result, soil erosion and deposition are pronounced, vegetation recovery is rapid, and human modification is extensive across the region. To better understand land–atmosphere interactions in this important and complex region, research communities supported by the U.S. Department of Energy’s Biological and Environmental Research (BER) program identified the Southeast as a priority region of interest. In fall 2024, the third Atmospheric Radiation Measurement Mobile Facility (AMF3), one of three mobile monitoring facilities designed to collect atmospheric and climate data from undersampled regions around the world, will begin operations in northwestern Alabama’s Bankhead National Forest (BNF). The AMF3-BNF 5-year deployment, from 2024–2029, will monitor the effects of feedbacks among aerosols, clouds, and precipitation on plant physiology and canopy-scale fluxes. It will also focus on scale aggregation to resolve the role of local forcing on larger-scale processes. To enable broader AMF3 involvement by the science community, the BER Environmental System Science (ESS) program organized the Southeast Land– Atmosphere Research Opportunities (SELARO) workshop in August 2023. The purpose was to identify gaps in scientific understanding of terrestrial processes in the Southeast (defined as states bounded by the Gulf of Mexico to the south, the Atlantic Ocean to east, the Mississippi River to the west, and extending through Tennessee and North Carolina to the north) and explore opportunities to use the AMF3-BNF deployment to coordinate and leverage research efforts across the region. Many parts of the Southeast have experienced repeated anthropogenic forcings. Farming, hunting, burning, and settlement of the region by Indigenous Peoples first shaped the distribution of plant communities, which in turn influenced European colonization patterns. Timber harvesting was common during the expansion of European settlements, and production forestry continues today. Agricultural production was extensive and then waned through the 20th century, creating a period of afforestation following agricultural abandonment. Today, many formerly agricultural landscapes are undergoing rapid urbanization and suburbanization. Overlying these patterns of anthropogenic land use are frequent disturbances from hurricanes, tornadoes, wildfires, drought, flooding, ice storms, and the occasional blizzard. An additional characteristic of the Southeast is its overall landscape complexity. Unlike the western United States, where broad expanses may share similar characteristics, Southeast topography, drainage patterns, vegetation, and development patterns vary widely across relatively small spatial scales (<1 km). This is due to the region’s underlying geology and soil development, species biodiversity patterns, and land ownership and use coupled with strong forces of erosion, weathering, and rapid plant growth in the warm, wet climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Heterogeneous Growth of UiO-66-NH 2 on Oxidized Single-Walled Carbon Nanotubes to Form “Beads-on-a-String” Composites

In this work we demonstrate a facile synthesis of UiO-66-NH 2 metal-organic framework (MOF)/oxidized single-walled carbon nanotubes (ox-SWCNTs) composite at room temperature. Acetic acid (HAc) was used as a modulator to manipulate the morphology of the MOF in these composites. With zirconium oxide-cluster (Zr) to 2-aminoteraphthalate linker (ATA) 1:1.42 ratio and acetic acid modulator, we achieved predominately heterogeneous MOF growth on the sidewalls of CNTs. Understanding the growth mechanism of these composites was facilitated by conducting DFT calculations to investigate the interactions between ox-SWCNTs and the MOF precursors. The synthesized composites combine both micro-porosity of the MOF and electrical conductivity of the SWCNTs. Gas sensing tests demonstrated higher response for UiO-66-NH 2 /ox-SWCNT hybrid toward dry air saturated with dimethyl methylphosphonate (DMMP) vapor compared to oxidized single-walled carbon nanotubes (ox-SWCNTs) alone.

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

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Representing surface heterogeneity in land–atmosphere coupling in E3SMv1 single-column model over ARM SGP during summertime

Abstract. The Earth's land surface features spatial and temporal heterogeneity over a wide range of scales below those resolved by current Earth system models (ESMs). State-of-the-art land and atmosphere models employ parameterizations to represent their subgrid heterogeneity, but the land–atmosphere coupling in ESMs typically operates on the grid scale. Communicating the information on the land surface heterogeneity with the overlying atmospheric boundary layer (ABL) remains a challenge in modeling land–atmosphere interactions. In order to account for the subgrid-scale heterogeneity in land–atmosphere coupling, we implement a new coupling scheme in the Energy Exascale Earth system model version 1 (E3SMv1) that uses adjusted surface variances and covariance of potential temperature and specific water content as the lower boundary condition for the atmosphere model. The new lower boundary condition accounts for both the variability of individual subgrid land surface patches and the inter-patch variability. The E3SMv1 single-column model (SCM) simulations over the Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site were performed to assess the impacts. We find that the new coupling parameterization increases the magnitude and diurnal cycle of the temperature variance and humidity variance in the lower ABL on non-precipitating days. The impacts are primarily attributed to subgrid inter-patch variability rather than the variability of individual patches. These effects extend vertically from the surface to several levels in the lower ABL on clear days. We also find that accounting for surface heterogeneity increases low cloud cover and liquid water path (LWP). These cloud changes are associated with the change in cloud regime indicated by the skewness of the probability density function (PDF) of the subgrid vertical velocity. In precipitating days, the inter-patch variability reduces significantly so that the impact of accounting for surface heterogeneity vanishes. These results highlight the importance of accounting for subgrid heterogeneity in land–atmosphere coupling in next-generation ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Kinetically versus thermodynamically controlled factors governing elementary pathways of GaP(111) surface oxidation

GaP and related III-V semiconductors have attracted interest as high-efficiency photoelectrochemical electrodes for solar water splitting. However, their efficacy is linked to the presence, identity, and integrity of native surface oxides, which are structurally and chemically complex and evolve during operation. Using ambient pressure X-ray photoelectron spectroscopy (APXPS) coupled with ab initio simulations, we track key chemical motifs expressed during evolution of GaP(111) surface oxides, their associated reaction kinetics, and their correlations with electronic properties. We identify two distinct thermal regimes corresponding to kinetically and thermodynamically controlled oxidation. Below 600 K, exposure to O2 generates kinetically facile Ga-O-Ga configurations, whereas higher temperatures cause activated oxygen to insert into Ga–P bonds as part of a thermodynamically driven transformation into a complex, heterogeneous 3D network of surface POx (1≤x ≤ 4) groups and Ga2O3 species, the latter of which eventually dominates upon depletion of surface phosphorus. Here our study highlights the critical competition between kinetic and thermodynamic factors during GaP oxidation, yielding insights for fabricating stable III-P-based photoelectrodes with precisely engineered surface properties.

08 HYDROGEN↗
DOE OSTI2026

Quenched disorder in the triangular lattice antiferromagnet YbZn 2⁢ GaO 5

We investigate the crystal electric field (CEF) excitations of Yb 3+ ions in powder samples of the triangular-lattice rare-earth-based antiferromagnet YbZn 2⁢ GaO 5 using inelastic neutron scattering (INS). Three CEF excitations from the ground-state Kramers doublet were observed, each exhibiting significant broadening beyond instrumental resolution. Combining temperature-dependent INS and neutron powder diffraction, we identify a significant static contribution to this broadening and attribute it to heterogeneous coordination of Yb 3+ ions due to Ga 3+ /Zn 2+ site mixing. Rietveld refinement of neutron powder diffraction indicates that 35% of Ga occupies the Zn site and 60% of Zn occupies the Ga site. We show with a point charge model for the CEF Hamiltonian that heterogeneous coordination of Yb3+ ions leads to broadened CEF peaks. First-principles calculations demonstrate that the random Ga 3+ /Zn 2+ distribution can produce the distortions of the YbO 6 octahedra observed from neutron diffraction. Because the documented heterogeneity will extend to exchange interactions, our results suggest that disorder is a significant factor in the unusual magnetism previously reported in YbZn 2 ⁢GaO 5 , including broad low-energy magnetic excitations and the absence of magnetic ordering down to 0.3 K.

Crystal field excitations↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2024

An uncertain future change in aridity over the tropics

Abstract An ensemble of climate models from phase six of the Coupled Model Intercomparison Project shows that temperature and potential evapotranspiration (PET) are projected to increase globally towards the end of the 21st century. However, climate models show a spatially heterogeneous change in precipitation over the tropics. Consequently, future changes in aridity (a measure of water availability) are complex and location-dependent. We assess future changes in aridity using three climate models and several single-forcing experiments. Near-term (2021–2040) changes in aridity are small, and we focus instead on its long-term (2081–2100) changes. We show that the increase in greenhouse gases (GHG) primarily explains the spatial pattern, magnitude and ensemble spread of the long-term future changes in aridity. On this timescale, the effects of changes in emissions of anthropogenic aerosols are moderate compared to the effects of increases in atmospheric GHG concentrations. Model diversity in the responses to GHG concentration is large over northern Africa and North and South America. We suggest the large uncertainty is due to differences between models in simulating the effects of an increase in GHG concentrations on surface air temperature over the North Atlantic Ocean, on the interhemispheric temperature gradient, and on PET over North and South America.

Environmental Sciences & Ecology↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗