Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural model decomposition”

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 451 records · Page 25

Nutrient Dynamics in a Coupled Terrestrial Biosphere and Land Model (ELM-FATES-CNP)

We present a representation of nitrogen and phosphorus cycling in the Functionally Assembled Terrestrial Ecosystem Simulator, a demographic vegetation model within the Energy Exascale Earth System land model. This representation is modular, and designed to allow testing of multiple hypothetical approaches for carbon-nutrient coupling in plants. Novel model hypotheses introduced in this work include, (a) the controls on plant acquisition of aqueous mineralized nutrients in the soil and (b) fairly straight forward methods of allocating nutrients to specific plant organs and their losses through live plant turnover as well as litter fluxes generated through plant mortality. This combines the new with pre-existing hypotheses (such as nitrogen fixation and soil decomposition) into a system that can accommodate plant-soil dynamics for a large number of size- and functional-type-resolved plant cohorts within a time-since-disturbance-resolved ecosystem. Root uptake of nutrients is governed by fine root biomass, and plants vary in their fine root biomass allocation in order to balance carbon and nutrient limitations to growth. We test the sensitivity of the model to a wide range of parameter variations and structural representations, and in the context of observations at Barro Colorado Island, Panama. A key model prediction is that plants in the high-light-availability canopy positions allocate more carbon to fine roots than plants in low-light understory environments, given the widely different carbon versus nutrient constraints of these two niches within a given ecosystem. This model provides a basis for exploring carbon-nutrient coupling with vegetation demography within Earth system models.

54 ENVIRONMENTAL SCIENCES↗

A Programming Environment Evaluation Methodology for Object-Oriented Systems

The object-oriented design strategy as both a problem decomposition and system development paradigm has made impressive inroads into the various areas of the computing sciences. Substantial development productivity improvements have been demonstrated in areas ranging from artificial intelligence to user interface design. However, there has been very little progress in the formal characterization of these productivity improvements and in the identification of the underlying cognitive mechanisms. The development and validation of models and metrics of this sort require large amounts of systematically-gathered structural and productivity data. There has, however, been a notable lack of systematically-gathered information on these development environments. A large part of this problem is attributable to the lack of a systematic programming environment evaluation methodology that is appropriate to the evaluation of object-oriented systems.

Dominick, Wayne D.↗

The response of soil carbon storage and microbially mediated carbon turnover to simulated climatic disturbance in a northern peatland forest. Revisiting the concept of soil organic matter recalcitrance

The goal of this project was to investigate changes in the structure of dissolved and solid phase organic matter, the production of CO 2 and CH 4 , and the composition of decomposer microbial communities in response to the climatic forcing of environmental processes that determine the balance between carbon gas production versus storage and sequestration in peatlands. Cutting-edge analytical chemistry and next generation sequencing of microbial genes were been applied to habitats at the Marcell Experimental Forest (MEF), where the US DOE’s Oak Ridge National Laboratory and the USDA Forest Service are constructing a large-scale ecosystem study entitled, “Spruce and Peatland Responses Under Climatic and Environmental Change”(SPRUCE). Our study represented a comprehensive characterization of the sources, transformation, and decomposition of organic matter in the S1 bog at MEF. Multiple lines of evidence point to distinct, vertical zones of organic matter transformation: 1) the acrotelm consisting of living mosses, root material, and newly formed litter (0-30 cm), 2) the mesotelm, a mid-depth transition zone (30-75 cm) characterized by labile organic C compounds and intense decomposition, and 3) the underlying catotelm (below 75cm) characterized by refractory organic compounds as well as relatively low decomposition rates. These zones are in part defined by physical changes in hydraulic conductivity and water table depth. O-alkyl-C, which represents the carbohydrate fraction in the peat, was shown to be an excellent proxy for soil decomposition rates. The carbon cycle in deep peat was shown to be fueled by modern carbon sources further indicating that hydrology and surface vegetation play a role in belowground carbon cycling. We provide the first metagenomic study of an ombrotrophic peat bog, with novel insights into microbial specialization and functions in this unique terrestrial ecosystem. Vertical structuring of microbial communities closely paralleled the chemical evolution of peat, with large shifts in microbial populations occurring in the biogeochemical hotspot, the mesotelm, where the highest rates of decomposition were detected. Stable isotope geochemistry and potential rates of methane production paralleled vertical changes in methanogen community composition to indicate a predominance of acetoclastic methanogenesis mediated by the Methanosarcinales in the mesotelm, while hydrogen-utilizing methanogens dominated in the deeper catotelm. Evidence pointed to the availability of phosphorus as well as nitrogen limiting the microbially-mediated turnover of organic carbon at MEF. Prior to initiation of the experimental treatments, our study provided key baseline data for the SPRUCE site on the vertical stratification of peat decomposition, key enzymatic pathways, and microbial taxa containing these pathways. The sensitivity of soil carbon turnover to climate change is strongly linked to recalcitrant carbon stocks and the temperature sensitivity of decomposition is thought to increase with increasing molecular complexity of carbon substrates. This project delivered results on how climate change perturbations impact the microbially-mediated turnover of recalcitrant organic matter in peatland forest soils, both under controlled conditions in the laboratory and at the ecosystem-scale in the field. This project revisited the concept of “recalcitrance” in the regulation of soil carbon turnover using a combination of natural abundance radiocarbon and optical spectroscopic measurements on bulk DOM, and high resolution molecular characterization of DOM. The project elucidated how organic matter reactivity and decomposition will respond to climate change in a both a qualitative (organic matter lability) and quantitiative (increased rates) manner. An Aromaticity Index was developed to represent a more direct and accurate parameter for modeling of DOM reactivity in peatlands. The abundance and community composition of soil microorganisms that mediate C cycling were interrogated with depth in the peat, with season, and in manipulated climate enclosures at unprecedented resolution. Therefore this project delivered strategic new insights on the functioning of peatland ecosystems that collectively store approximately one-third of the world's soil carbon. Furthermore, results from the detailed characterization of DOM lability and microbial community structure/ function will be employed to further develop biogeochemical models to include microbial respiration pathways as well as to track carbon flow with a term that incorporates relative reactivity based on aromaticity index. As it stands now, detailed soil organic matter structure and microbial parameters are not included in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

Ion Pairing and Redissociaton in Low-Permittivity Electrolytes for Multivalent Battery Applications

Detailed speciation of electrolytes as a function of chemical system and concentration provides the foundation for understanding bulk transport as well as possible decomposition mechanisms. In particular, multivalent electrolytes have shown a strong coupling between anodic stability and solvation structure. Furthermore, solvents that are found to exhibit reasonable stability against alkaline-earth metals generally exhibit low permittivity, which typically increases the complexity of the electrolyte species. To improve our understanding of ionic population and associated transport in these important classes of electrolytes, the speciation of Mg(TFSI)2 in monoglyme and diglyme systems is studied via a multiscale thermodynamic model using first-principles calculations for ion association and molecular dynamics simulations for dielectric properties. The results are then compared to Raman and dielectric relaxation spectroscopies, which independently confirm the modeling insights. We find that the significant presence of free ions in the low-permittivity glymes in the concentration range from 0.02 to 0.6 M is well-explained by the low-permittivity redissociation hypothesis. Here, salt speciation is largely dictated by long-range electrostatics, which includes permittivity increases due to polar contact ion pairs. The present results suggest that other low-permittivity multivalent electrolytes may also reach high conductivities as a result of redissociation.

25 ENERGY STORAGE↗

The phase and turbulent properties of Cattail

ABSTRACT A recent publication discovered one of the largest filamentary neutral hydrogen features dubbed Cattail from high-resolution Five-hundred-meter Aperture Spherical radio Telescope observations that might be a new galactic arm of the Milky Way. We evaluate the turbulent and phase properties of Cattail via the newly developed Velocity Decomposition Algorithm and Force Balancing Model. We discover that if there exists a phase transition, then Cattail is unlikely in the cold neutral media phase. We also show that the Cattail is two disjoint features in caustics space, suggesting that the Cattail has two different turbulent systems. We check the spectra of the individual system separated via VDA to confirm this argument. We do not exclude the existence of smaller scale cold media being embedded within this structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Zero to Integration in Eight Months, the Dawn Ground Data System Engineering Challange

The Dawn Project has presented the Ground Data System (GDS) with technical challenges driven by cost and schedule constraints commonly associated with National Aeronautics and Space Administration (NASA) Discovery Projects. The Dawn mission consists of a new and exciting Deep Space partnership among: the Jet Propulsion Laboratory (JPL), responsible for project management and flight operations; Orbital Sciences Corporation (OSC), spacecraft builder and responsible for flight system test and integration; and the University of California, at Los Angeles (UCLA), responsible for science planning and operations. As a cost-capped mission, one of Dawn s implementation strategies is to leverage from both flight and ground heritage. OSC's ground data system is used for flight system test and integration as part of the flight heritage strategy. Mission operations, however, are to be conducted with JPL s ground system. The system engineering challenge of dealing with two heterogeneous ground systems emerged immediately. During the first technical interchange meeting between the JPL s GDS Team and OSC's Flight Software Team, August 2003, the need to integrate the ground system with the flight software was brought to the table. This need was driven by the project s commitment to enable instrument engineering model integration in a spacecraft simulator environment, for both demonstration and risk mitigation purposes, by April 2004. This paper will describe the system engineering approach that was undertaken by JPL's GDS Team in order to meet the technical challenge within a non-negotiable eight-month schedule. Key to the success was adherence to an overall systems engineering process and fundamental systems engineering practices: decomposition of the project request into manageable requirements; definition of a structured yet flexible development process; integration of multiple ground disciplines and experts into a focused team effort; in-process risk management; and aggregation of the intermediate products to an integrated final product. In addition, this paper will highlight the role of lessons learned from the integration experience. The lessons learned from an early GDS deployment have served as the foundation for the design and implementation of the Dawn Ground Data System.

systems engineering↗

Structural Analysis of Novel Lignin-derived Carbon Composite Anodes

The development of novel lignin-based carbon composite anodes consisting of nanocrystalline and amorphous domains motivates the understanding of a relationship of the structural properties characterizing these materials, such as crystallite size, intracrystallite dspacing, crystalline volume fraction and composite density, with their pair distribution functions (PDF), obtained from both molecular dynamics simulation and neutron scattering. A model for these composite materials is developed as a function of experimentally measurable parameters and realized in fifteen composite systems, three of which directly match all parameters of their experimental counterparts. The accurate reproduction of the experimental PDFs using the model systems validates the model. The decomposition of the simulated PDFs provides an understanding of each feature in the PDF and allows for the development of a mapping between the defining characteristics of the PDF and the material properties of interest.

McNutt, Nicholas W↗

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗

Testing the thermal Sunyaev-Zel’dovich power spectrum of a halo model using hydrodynamical simulations

Statistical properties of large-scale cosmological structures serve as powerful tools for constraining the cosmological properties of our Universe. Tracing the gas pressure, the thermal Sunyaev-Zel’dovich (tSZ) effect is a biased probe of mass distribution and, hence, can be used to test the physics of feedback or cosmological models. Therefore, it is crucial to develop robust modelling of hot gas pressure for applications to tSZ surveys. Since gas collapses into bound structures, it is expected that most of the tSZ signal is within halos produced by cosmic accretion shocks. Hence, simple empirical halo models can be used to predict the tSZ power spectra. In this study, we employed the HMx halo model to compare the tSZ power spectra with those of several hydrodynamical simulations: the Horizon suite and the Magneticum simulation. We examine various contributions to the tSZ power spectrum across different redshifts, including the one- and two-halo term decomposition, the amount of bound gas, the importance of different masses, and the electron pressure profiles. Our comparison of the tSZ power spectrum reveals discrepancies between the halo model and cosmological simulations that increase with redshift. We find a 20% to 50% difference between the measured and predicted tSZ angular power spectrum over the multipole range ℓ = 10 3 − 10 4 . Our analysis reveals that these differences are driven by the excess of power in the predicted two-halo term at low k and in the one-halo term at high k . At higher redshifts ( z ∼ 3), simulations indicate that more power comes from outside the virial radius than from inside, suggesting a limitation in the applicability of the halo model. We also observe differences in the pressure profiles, despite the fair level of agreement on the tSZ power spectrum at low redshift with the default calibration of the halo model. In conclusion, our study suggests that the properties of the halo model need to be carefully controlled against real or mock data to be proven useful for cosmological purposes.

Ayçoberry, Emma (ORCID:0000000292351195)↗

Spin-wave multiple excitations in nanoscale classical Heisenberg antiferromagnets

Monte Carlo and spin dynamics techniques have been used to perform large-scale simulations of the dynamic behavior of a nanoscale, classical, Heisenberg antiferromagnet on a simple-cubic lattice with linear sizes L≤ 40 at a temperature below the Neel temperature. In this study, nanoparticles are modeled with completely free boundary conditions, i.e., six free surfaces, and nanofilms are modeled with two free surfaces in the spatial z direction and periodic boundaries parallel to the surfaces in the xy direction, which are compared to the infinite system with periodic boundary conditions. The temporal evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using a fast spin dynamics algorithm with the fourth-order Suzuki-Trotter decomposition of exponential operators, with initial spin configurations generated by Monte Carlo simulations. The local dynamic structure factor S(q,ω) was calculated from the local space- and time-displaced spin-spin correlation function. Multiple excitation peaks for wave vectors within the first Brillouin zone appear in the spin-wave spectra of the transverse component of dynamic structure factor S T (q,ω) in the nanoscale classical Heisenberg antiferromagnet, which are lacking if periodic boundary conditions are used. With the assumption of q-space spin-wave reflections with broken momentum conservation due to free-surface confinements, we successfully explained those spectra quantitatively in the linear dispersion region. Meanwhile, we also observed two unexpected quantized spin-wave excitation modes in the spatial z direction in nanofilms for S T (q,ω) not expected in bulk systems. In conclusion, the results of this study indicate the presence of unexpected forms of spin-wave excitation behavior that have yet to be observed experimentally but could be directly tested through neutron scattering experiments on nanoscale RbMnF 3 particles or films.

36 MATERIALS SCIENCE↗

Sensitivity method for integrated structure/active control law design

The development is described of an integrated structure/active control law design methodology for aeroelastic aircraft applications. A short motivating introduction to aeroservoelasticity is given along with the need for integrated structures/controls design algorithms. Three alternative approaches to development of an integrated design method are briefly discussed with regards to complexity, coordination and tradeoff strategies, and the nature of the resulting solutions. This leads to the formulation of the proposed approach which is based on the concepts of sensitivity of optimum solutions and multi-level decompositions. The concept of sensitivity of optimum is explained in more detail and compared with traditional sensitivity concepts of classical control theory. The analytical sensitivity expressions for the solution of the linear, quadratic cost, Gaussian (LQG) control problem are summarized in terms of the linear regulator solution and the Kalman Filter solution. Numerical results for a state space aeroelastic model of the DAST ARW-II vehicle are given, showing the changes in aircraft responses to variations of a structural parameter, in this case first wing bending natural frequency.

Gilbert, Michael G.↗

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗

The microbe-mediated mechanisms affecting topsoil carbon stock in Tibetan grasslands

Warming has been shown to cause soil carbon (C) loss in northern grasslands owing to accelerated microbial decomposition that offsets increased grass productivity. Yet, a multi-decadal survey indicated that the surface soil C stock in Tibetan alpine grasslands remained relatively stable. To investigate this inconsistency, we analyzed the feedback responses of soil microbial communities to simulated warming by soil transplant in Tibetan grasslands. Microbial functional diversity decreased in response to warming, whereas microbial community structure did not correlate with changes in temperature. The relative abundance of catabolic genes associated with nitrogen (N) and C cycling decreased with warming, most notably in genes encoding enzymes associated with more recalcitrant C substrates. By contrast, genes associated with C fixation increased in relative abundance. The relative abundance of genes associated with urease, glutamate dehydrogenase and ammonia monoxygenase ( ureC , gdh and amoA ) were significantly correlated with N 2 O efflux. These results suggest that unlike arid/semiarid grasslands, Tibetan grasslands maintain negative feedback mechanisms that preserve terrestrial C and N pools. To examine whether these trends were applicable to the whole plateau, we included these measurements in a model and verified that topsoil C stocks remained relatively stable. Thus, by establishing linkages between microbial metabolic potential and soil biogeochemical processes, we conclude that long-term C loss in Tibetan grasslands is ameliorated by a reduction in microbial decomposition of recalcitrant C substrates.

54 ENVIRONMENTAL SCIENCES↗

Nonlinear Ultrasonic Diagnosis and Prognosis of ASR Damage in Dry Cask Storage

Alkali-silica reaction (ASR) is a deleterious chemical process that may occur in cement-based materials such as mortars and concretes, where the hydroxyl ions in the highly alkaline pore solution attack the siloxane groups in the siliceous minerals in the aggregates. The reaction produces a cross-linked alkali-silica gel. The ASR gel swells in the presence of water. Expansion of the gel results in cracking when the swelling-induced stress exceeds the fracture toughness of the concrete. As the ASR continues, cracks may grow and eventually coalesce, which results in reduced service life and a decrease safety of concrete structures. Since concrete is widely used as a critical structural component in dry cask storage of used nuclear fuels, ASR damage poses a significant threat to the sustainability of long term dry cask storage systems. Therefore, techniques for effectively detecting, managing and mitigating ASR damage are needed. Currently, there are no nondestructive methods to accurately detect ASR damage in existing concrete structures. The only current way of accurately assessing ASR damage is to drill a core from an existing structure, and conduct microscopy on this drilled cylindrical core. Clearly, such a practice is not applicable to dry cask storage systems. To meet these needs, this research is aimed at developing (1) a suite of nonlinear ultrasonic quantitative nondestructive evaluation (QNDE) techniques to characterize ASR damage, and (2) a physics-based model for ASR damage evolution using the QNDE data. Outcomes of this research will provide a nondestructive diagnostic tool to evaluate the extent of the ASR damage, and a prognostic tool to estimate the future reliability and safety of the concrete structures in dry cask storage systems

42 ENGINEERING↗

Experimental and computing strategies in advanced material characterization problems

The mechanical characterization of materials relies more and more often on sophisticated experimental methods that permit to acquire a large amount of data and, contemporarily, to reduce the invasiveness of the tests. This evolution accompanies the growing demand of non-destructive diagnostic tools that assess the safety level of components in use in structures and infrastructures, for instance in the strategic energy sector. Advanced material systems and properties that are not amenable to traditional techniques, for instance thin layered structures and their adhesion on the relevant substrates, can be also characterized by means of combined experimental-numerical tools elaborating data acquired by full-field measurement techniques. In this context, parameter identification procedures involve the repeated simulation of the laboratory or in situ tests by sophisticated and usually expensive non-linear analyses while, in some situation, reliable and accurate results would be required in real time. The effectiveness and the filtering capabilities of reduced models based on decomposition and interpolation techniques can be profitably used to meet these conflicting requirements. This communication intends to summarize some results recently achieved in this field by the author and her co-workers. The aim is to foster further interaction between engineering and mathematical communities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real Time Observation of Lithium Insertion into Pre-Cycled Conversion-Type Materials

Conversion-type electrode materials for lithium-ion batteries experience significant structural changes during the first discharge–charge cycle, where a single particle is taken apart into a number of nanoparticles. This structural evolution may affect the following lithium insertion reactions; however, how lithiation occurs in pre-cycled electrode materials is elusive. In this work, in situ transmission electron microscopy was employed to see the lithium-induced structural and chemical evolutions in pre-cycled nickel oxide as a model system. The introduction of lithium ions induced the evolution of metallic nickel, with volume expansion as a result of a conversion reaction. After pre-cycling, the phase evolutions occurred in two separate areas almost at the same time. This is different from the first lithiation, where the phase change takes place successively, with a boundary dividing the reacted and unreacted areas. Structural changes were restricted at the areas having large amount of fluorine, implying the residuals from the decomposition of electrolytes may have hindered the electrochemical reactions. This work provides insights into phase and chemical evolutions in pre-cycled conversion-type materials, which govern electrochemical properties during operation.

25 ENERGY STORAGE↗

Coarsening of complex microstructures following spinodal decomposition

Coarsening plays a pivotal role in materials engineering, but our understanding of the dynamics of coarsening in morphologically complex systems is still limited. In this paper, we examine the correlations between the interfacial velocity and interfacial morphologies, and then predict the evolution of mean curvature based on the correlations. Three simulated structures with varying volume fractions, two bicontinuous and one nonbicontinuous, are generated using the Cahn-Hilliard equation. We find general correlations between interfacial velocity and mean curvature, as well as between interfacial velocity and the surface Laplacian of the mean curvature. Furthermore, we find that the probability of finding a patch of interface with a given normal velocity and the same local principal curvatures is described well by a Gaussian distribution, independent of the principal curvature values and the volume fractions of the structures. We also find that average interfacial velocity is described by a polynomial of the mean curvature and the net curvature. Based on this finding, we develop a semi-analytical approach to predicting the rate of change of the mean curvature, which determines the morphological evolution of complex microstructures.

36 MATERIALS SCIENCE↗

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↗