Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural dynamics and control”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

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

Cations Control Lipid Bilayer Memcapacitance Associated with Long-Term Potentiation

Phospholipid bilayers can be described as capacitors whose capacitance per unit area (specific capacitance, C m ) is determined by their thickness and dielectric constant–independent of applied voltage. It is also widely assumed that the C m of membranes can be treated as a “biological constant”. Recently, using droplet interface bilayers (DIBs), it was shown that zwitterionic phosphatidylcholine (PC) lipid bilayers can act as voltage-dependent, nonlinear memory capacitors, or memcapacitors. When exposed to an electrical “training” stimulation protocol, capacitive energy storage in lipid membranes was enhanced in the form of long-term potentiation (LTP), which enables biological learning and long-term memory. LTP was the result of membrane restructuring and the progressive asymmetric distribution of ions across the lipid bilayer during training, which is analogous, for example, to exponential capacitive energy harvesting from self-powered nanogenerators. Here, we describe how LTP could be produced from a membrane that is continuously pumped into a nonequilibrium steady state, altering its dielectric properties. During this time, the membrane undergoes static and dynamic changes that are fed back to the system’s potential energy, ultimately resulting in a membrane whose modified molecular structure supports long-term memory storage and LTP. Here, we also show that LTP is very sensitive to different salts (KCl, NaCl, LiCl, and TmCl 3 ), with LiCl and TmCl 3 having the most profound effect in depressing LTP, relative to KCl. This effect is related to how the different cations interact with the bilayer zwitterionic PC lipid headgroups primarily through electric-field-induced changes to the statistically averaged orientations of water dipoles at the bilayer headgroup interface.

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

Dislocation-induced stop-and-go kinetics of interfacial transformations

Most engineering materials are based on multiphase microstructures produced either through the control of phase equilibria or by the fabrication of different materials as in thin-film processing. In both processes, the microstructure relaxes towards equilibrium by mismatch dislocations (or geometric misfit dislocations) across the heterophase interfaces (1-5) . Despite their ubiquitous presence, directly probing the dynamic action of mismatch dislocations has been unachievable owing to their buried nature. In this work, using the interfacial transformation of copper oxide to copper as an example, we demonstrate the role of mismatch dislocations in modulating oxide-to-metal interfacial transformations in an intermittent manner, by which the lateral flow of interfacial ledges is pinned at the core of mismatch dislocations until the dislocation climbs to the new oxide/metal interface location. Together with atomistic calculations, we identify that the pinning effect is associated with the non-local transport of metal atoms to fill vacancies at the dislocation core. These results provide mechanistic insight into solid-solid interfacial transformations and have substantial implications for utilizing structural defects at buried interfaces to modulate mass transport and transformation kinetics.

25 ENERGY STORAGE↗
DOE OSTI2024

Electromagnetic Energy-Assisted Thermal Conversion of Fossil-Based Hydrocarbons to Low-Cost Hydrogen

The goal of this project was to develop and optimize catalysts for methane decomposition, particularly focusing on regeneration via an electromagnetic energy-assisted mechanism, to produce hydrogen more cost-effectively compared to electrolysis routes. To achieve this goal, the project pursued several key objectives. The project began with the preparation and testing of various catalysts. A nickel-silica based catalyst was identified as the most promising material for the pyrolysis of methane into carbon and hydrogen. Kinetic parameters for methane decomposition were determined, aiding in computational modeling efforts. Structured catalysts were investigated, highlighting the need for frequent cleaning or regeneration to maintain performance, with methane conversion rates exceeding 70% in tube furnace tests. Computational fluid dynamics modeling was employed to optimize reactor designs and electrode angles, leading the project team to propose a multi-compartment thermal conversion system for larger setups. This modeling work was important in understanding reaction characteristics, carbon deposition rates, and temperature profiles under various conditions. A bench-scale reactor system was assembled to evaluate catalyst regeneration using electromagnetic energy-assisted mechanisms. Experiments demonstrated the potential for carbon removal, though further optimization is needed. The carbon produced from the methane conversion process was evaluated for potential use in lithium-ion battery electrodes. The carbon exhibited properties similar to commercially available high-purity multi-walled carbon nanotubes and nanofibers, with a carbon content greater than 95%. Coin cell batteries assembled with this carbon showed that lower replacement levels (10% to 33%) outperformed the control group, improving specific capacity density and stability. However, higher replacement levels (100%) demonstrated poorer performance, suggesting that excessive carbon substitution negatively impacts battery performance. These findings indicate the potential marketability of the produced carbon as a component in lithium-ion batteries, though further testing is necessary to confirm long-term advantages and disadvantages associated with the use of the carbon product. These results however justified further technoeconomic assessments to determine if the process can provide low-cost hydrogen. The economic feasibility and technical performance of methane decomposition for hydrogen production were assessed, focusing on three plant configurations: 100E (electrically heated), 100C (combustion heated using produced hydrogen), and PE-Hybrid (a combination of pyrolysis (indicating decomposition) and electrolysis). The Levelized Cost of Hydrogen (LCOH) for the pyrolysis configurations was found to be approximately 25% lower than that of electrolysis. The 100E configuration had the lowest LCOH at $\$$3.12/kg. Including carbon product sales significantly improved the economics, with the 100C configuration achieving a negative LCOH of -$\$$0.35/kg. The PE-Hybrid configuration was not economically advantageous compared to pure pyrolysis plants due to its complexity and additional equipment requirements. Ultimately, methane pyrolysis presents a viable method for near carbon dioxide-free hydrogen production, with significant economic advantages over electrolysis, especially when considering the sale of carbon byproducts. The 100E and 100C configurations showed the most promise, with the choice between them ultimately depending on the prices of power and natural gas. In conclusion, this technology has the potential to lower hydrogen production costs by leveraging the methane decomposition process with the sale of valuable carbon byproducts. By optimizing catalyst performance and integrating electromagnetic energy-assisted regeneration, the process can achieve higher efficiency and economic viability, making it a competitive alternative to traditional hydrogen production methods.

08 HYDROGEN↗
DOE OSTI2021

Bridging the gap between NMR measured mean silicate chain length and nano-scale silicate polymorphism of calcium silicate hydrates

Variable stoichiometry and silicate polymorphism in calcium-silicate-hydrates (C-S-H) has impeded the revelation of point defects distribution in the silicate tetrahedral network of C-S-H, which resembles tobermorite crystal structure with some bridging tetrahedra (BT) and paired tetrahedra (PT) vacancies in dreierketten chains. Here we use a computational approach to characterize silicate polymorphism by introducing the vacancy ratio of BT to PT (α) and establishing a three-term empirical mean chain length (MCL) formula for different calcium-to‑silicon ratios (Ca/Si) ranging from 1.2 to 2.3. The formula identifies BT and PT controls at low and high α respectively and allows an inverse mapping of polymorphic range based on NMR experiments. The proposed computational framework quantitatively describes silicate polymorphs and links NMR-measured MCL to C-S-H atomistic configurations at the molecular level.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Bioproduction, bioprotection, and biocontainment in multi-kingdom microbial systems with 3D spatial control

Engineered living materials (ELMs) are a class of hybrid materials that include engineered microbes encapsulated by a polymer matrix. The biotic and abiotic components define the ELMs design space and can be altered to improve performance and function. While current synthetic materials in the field display robust biocompatibility with both native and engineered living systems, we have a limited understanding of how to leverage three-dimensional (3D) form factors to spatially organize and control microbial dynamics within the material. Motivated by this knowledge gap, we employed extrusion-based 3D printing to fabricate multi-kingdom hydrogel constructs for the encapsulation of both single and multi-kingdom microbial systems. Core–shell cubic constructs enabled the spatial organization of a constitutive multi-kingdom system of levodopa (L-DOPA)-producing E. coli and betaxanthins (BXN)-producing S. cerevisiae. This spatial organization in 3D materials can introduce precise control over bioproduction, bioprotection, and biocontainment features that are critical to the efficacy of current ELMs. The relative spatial organization of the organisms, as well as the surface area-to-volume ratio were investigated to determine how these design elements impact microbial behavior (metabolite production, growth, expression, and cell distribution) over time. We demonstrated that F127-bis-urethane methacrylate (F127-BUM) core–shell geometries enable the hierarchical 3D printing of multi-kingdom constructs, offering customizable control over bioproduction, bioprotection, and biocontainment. With the optimization of these core–shell structures for continuous bioproduction, these ELMs could be deployed as compact and sustainable bioreactors in remote environments.

additive manufacturing↗
DOE OSTI2025

Defect-Mediated Diffusion Pathways in Spodumene Accelerate Lithium Transport

Lithium extraction from naturally occurring α-spodumene is hindered by poor lithium diffusivity, necessitating high-temperature phase transformation to a low-density β polymorph. Although β spodumene exhibits up to 5 orders of magnitude higher lithium-ion diffusivity, both phases have diffusion activation energies between 0.8 and 1 eV, indicating that polymorph density is not the controlling factor over diffusivity. We show that aluminum vacancies facilitate lithium-ion diffusion in α-spodumene by reducing the migration barrier from 2.4 to 0.9 eV. Bond valence site energy and nudged elastic band calculations show a new lithium local minimum site which promotes a one-dimensional percolation network by reducing the lithium intersite distance from 4.5 Å to 2.9 Å. However, aluminum vacancies are energetically unfavorable to percolate through the whole structure, resulting in very low net lithium diffusivity and highlighting the critical role of nonstoichiometric defects in facilitating lithium transport in rigid aluminosilicate structures.

Chemical structure↗
DOE OSTI2024

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Coarse-grained molecular dynamics integrated with convolutional neural network for comparing shapes of temperature sensitive bottlebrushes

Quantification of shape changes in nature-inspired soft material architectures of stimuli-sensitive polymers is critical for controlling their properties but is challenging due to their softness and flexibility. Here, we have computationally designed uniquely shaped bottlebrushes of a thermosensitive polymer, poly(N-isopropylacrylamide) (PNIPAM), by controlling the length of side chains along the backbone. Coarse-grained molecular dynamics simulations of solvated bottlebrushes were performed below and above the lower critical solution temperature of PNIPAM. Conventional analyses (free volume, asphericity, etc.) show that lengths of side chains and their immediate environments dictate the compactness and bending in these architectures. We further developed 100 unique convolutional neural network models that captured molecular-level features and generated a statistically significant quantification of the similarity between different shapes. Thus, our study provides insights into the shapes of complex architectures as well as a general method to analyze them. The shapes presented here may inspire the synthesis of new bottlebrushes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical readout of singlet fission biexcitons in a heteroacene with photoluminescence detected magnetic resonance

Molecular spin systems based on photoexcited triplet pairs formed via singlet fission (SF) are attractive as carriers of quantum information because of their potentially pure and controllable spin polarization, but developing systems that offer optical routes to readout as well as initialization is challenging. Herein, we characterize the electron spin magnetic resonance change in the photoluminescence intensity for a tailored organic molecular crystal while sweeping a microwave drive up to 10 GHz in a broadband loop structure. We observe resonant transitions for both triplet and quintet spin sublevel populations showing their optical sensitivity and revealing the zero-field parameters for each. We map the evolution of these spectra in both microwave frequency and magnetic field, producing a pattern of optically detected magnetic resonance (ODMR) peaks. Fits to these data using a suitable model suggest significant spin polarization in this system with orientation selectivity. Unusual excitation intensity dependence is also observed, which inverts the sign of the ODMR signal for the triplet features, but not for the quintet. These observations demonstrate optical detection of the spin sublevel population dictated by SF and intermolecular geometry, and highlight anisotropic and multi-scale dynamics of triplet pairs.

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

The Relationship between Precipitation and Precipitable Water in CMIP6 Simulations and Implications for Tropical Climatology and Change

It is well documented that over the tropical oceans, column-integrated precipitable water (pw) and precipitation (P) have a nonlinear relationship. In this study moisture budget analysis is used to examine this P–pw relationship in a normalized precipitable water framework. It is shown that the parameters of the nonlinear relationship depend on the vertical structure of moisture convergence. Specifically, the precipitable water values at which precipitation is balanced independently by evaporation versus by moisture convergence define a critical normalized precipitable water, pw nc . This is a measure of convective inhibition that separates tropical precipitation into two regimes: a local evaporation-controlled regime with widespread drizzle and a precipitable water–controlled regime. Most of the 17 CMIP6 historical simulations examined here have higher pw nc compared to ERA5, and more frequently they operate in the drizzle regime. When compared to observations, they overestimate precipitation over the high-evaporation oceanic regions off the equator, thereby producing a ‘‘double ITCZ’’ feature, while underestimating precipitation over the large tropical landmasses and over the climatologically moist oceanic regions near the equator. The responses to warming under the SSP585 scenario are also examined using the normalized precipitable water framework. It is shown that the critical normalized precipitable water value at which evaporation versus moisture convergence balance precipitation decreases as a result of the competing dynamic and thermodynamic responses to warming, resulting in an increase in drizzle and total precipitation. Statistically significant historical trends corresponding to the thermodynamic and dynamic changes are detected in ERA5 and in lowintensity drizzle precipitation in the PERSIANN precipitation dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Optical and photoelectric properties of surface barrier diodes based on n- Hg3In2Te6 crystals

The high radiation resistance of Hg3In2Te6 (MIT) single crystals to large doses of ionizing radiation makes devices quite attractive. Such a feature of the crystal structure has not only important advantages, but also causes serious technological difficulties in doping and controlling the specific resistance within the required limits and type of conductivity of single crystals. These technological difficulties create serious problems when creating structures with a potential barrier, especially with a p-n junction. At the moment, there are no publications on the creation of a p-n junction based on MIT. The only alternative is to create a potential barrier such as a Schottky diode or a metal-dielectric-semiconductor structure. This limits the use of MIT to produce, for example, rectifier diodes that can operate at relatively high reverse bias voltages. The typical resistivity of single crystals of MIT of n-type conductivity, which are used for the manufacture of diodes and photodiodes, is approximately equal to ρ≈103-104 Ω×cm. Such a specific resistance leads to a relatively large resistance of the neutral part of the single crystal, which in turn limits the direct currents of the rectifier diodes and significantly limits an important parameter of photodiodes - the dynamic range of the watt-ampere characteristic. This problem arises from the relatively high series resistance of the substrate. For rectifier diodes, this relatively large series resistance leads to a significant forward bias voltage drop, which reduces both the rectification factor and the amount of rectified current. We have already created diodes and photodiodes on Cr/MIT/Cr single crystals with a specific resistance of ρ≈103-104 Ω×cm. The diodes had relatively small dark currents but had the disadvantage of a relatively large specific resistance of the substrate, which we have already discussed above [1,2,3].

Sklyarchuk, V.↗
DOE OSTI2020

Dynamically Reconfigurable Data Readout of Pixel Detectors for Automatic Synchronization with Data Acquisition Systems

Reconfigurable detectors with dynamically selectable sensing and readout modes are highly desirable for implementing edge computing as well as enabling advanced imaging techniques such as foveation. The concept of a camera system capable of simultaneous passive imaging and dynamic ranging in different regions of the detector is presented. Such an adaptive-autonomous detector with both spatial and temporal control requires programmable window of exposure (time frames), ability to switch between readout modes such as full-frame imaging and zero-suppressed data, modification of the number of pixel data bits and independent programmability for distinct detector regions. In this work, a method is presented for seamlessly changing time frames and readout modes without data corruption while still ensuring that the data acquisition system (DAQ) does not need to stop and resynchronize at each change of setting, thus avoiding significant dead time. Data throughput is maximized by using a minimum unique data format, rather than lengthy frame headers, to differentiate between consecutive frames. A data control and transmitter (DCT) synchronizes data transfer from the pixel to the periphery, reconfigures the data to transmit it serially off-chip, while providing optimized decision support based on a DAQ definable mode. Measurements on a test structure demonstrate that the DCT can operate at 1 GHz in a 65 nm LP CMOS process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A novel random forest approach to revealing interactions and controls on chlorophyll concentration and bacterial communities during coastal phytoplankton blooms

Abstract Increasing occurrence of harmful algal blooms across the land–water interface poses significant risks to coastal ecosystem structure and human health. Defining significant drivers and their interactive impacts on blooms allows for more effective analysis and identification of specific conditions supporting phytoplankton growth. A novel iterative Random Forests (iRF) machine-learning model was developed and applied to two example cases along the California coast to identify key stable interactions: (1) phytoplankton abundance in response to various drivers due to coastal conditions and land-sea nutrient fluxes, (2) microbial community structure during algal blooms. In Example 1, watershed derived nutrients were identified as the least significant interacting variable associated with Monterey Bay phytoplankton abundance. In Example 2, through iRF analysis of field-based 16S OTU bacterial community and algae datasets, we independently found stable interactions of prokaryote abundance patterns associated with phytoplankton abundance that have been previously identified in laboratory-based studies. Our study represents the first iRF application to marine algal blooms that helps to identify ocean, microbial, and terrestrial conditions that are considered dominant causal factors on bloom dynamics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Cations impact radical reaction dynamics in concentrated multicomponent aqueous solutions

Ultraviolet (UV) photolysis of nitrite ions (NO 2 – ) in aqueous solutions produces a suite of radicals, viz., NO · , O – , · OH, and · NO 2 . The O – and NO · radicals are initially formed from the dissociation of photoexcited NO 2 – . The O – radical undergoes reversible proton transfer with water to generate · OH. Both · OH and O – oxidize the NO 2 – to · NO 2 radicals. The reactions of · OH occur at solution diffusion limits, which are influenced by the nature of the dissolved cations and anions. Here, we systematically varied the alkali metal cation, spanning the range from strongly to weakly hydrating ions, and measured the production of NO · , · OH, and · NO 2 radicals during UV photolysis of alkaline nitrite solutions using electron paramagnetic resonance spectroscopy with nitromethane spin trapping. Comparing the data for the different alkali cations revealed that the nature of the cation had a significant effect on production of all three radical species. Radical production was inhibited in solutions with high charge density cations, e.g., lithium, and promoted in solutions containing low charge density cations, e.g., cesium. Through complementary investigations with multinuclear single pulse direct excitation nuclear magnetic resonance (NMR) spectroscopy and pulsed field gradient NMR diffusometry, cation-controlled solution structures and extent of NO 2 – solvation were determined to alter the initial yields of · NO and · OH radicals as well as alter the reactivity of NO 2 – toward · OH, impacting the production of · NO 2 . The implications of these results for the retrieval and processing of low-water, highly alkaline solutions that comprise legacy radioactive waste are discussed.

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

Engineering Relaxor Behavior in (BaTiO 3 ) n /(SrTiO 3 ) n Superlattices

Complex-oxide superlattices provide a pathway to numerous emergent phenomena because of the juxtaposition of disparate properties and the strong interfacial interactions in these unit-cell-precise structures. This is particularly true in superlattices of ferroelectric and dielectric materials, wherein new forms of ferroelectricity, exotic dipolar textures, and distinctive domain structures can be produced. Here, in this study, relaxor-like behavior, typically associated with the chemical inhomogeneity and complexity of solid solutions, is observed in (BaTiO 3 ) n /(SrTiO 3 ) n (n = 4–20 unit cells) superlattices. Dielectric studies and subsequent Vogel–Fulcher analysis show significant frequency dispersion of the dielectric maximum across a range of periodicities, with enhanced dielectric constant and more robust relaxor behavior for smaller period n. Bond-valence molecular-dynamics simulations predict the relaxor-like behavior observed experimentally, and interpretations of the polar patterns via 2D discrete-wavelet transforms in shorter-period superlattices suggest that the relaxor behavior arises from shape variations of the dipolar configurations, in contrast to frozen antipolar stripe domains in longer-period superlattices (n = 16). Moreover, the size and shape of the dipolar configurations are tuned by superlattice periodicity, thus providing a definitive design strategy to use superlattice layering to create relaxor-like behavior which may expand the ability to control desired properties in these complex systems.

36 MATERIALS SCIENCE↗
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↗