SEARCH · Engineering Papers
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.
We studied the electron dynamic correlation using inelastic X-ray scattering measurements. The data files contain a 2D array of numeric data that is used to produce the figures in the paper, "Observation of Real-space Dynamic Electron Correlation in Beryllium."
Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.
Iron is a widespread structural material in nuclear applications. However some benchmark studies show that for many libraries, including TENDL-2019, producing good iron evaluations is still an open challenge. This work aims at improving TENDL-2019 {sup 56}Fe and {sup 54}Fe evaluations through feedbacks from experimental data and Monte Carlo results of the iron benchmark ASPIS. Part of the T6 code infrastructure was used to produce the evaluations. In the new evaluations, particular attention to resonance parameters and formalism was paid to address the importance of elastic scattering. The total inelastic scattering and inelastic continuum cross sections were shaped by modifying the input parameters of TALYS, one of the T6 codes. Additionally, it was necessary to normalise some reaction channels to those of JENDL-4.0. With the new files ASPIS results, obtained with the Monte Carlo code TRIPOLI-4, improved noticeably compared to TENDL 2019. The evaluations and parameter setting modifications are to be included into TENDL-2021. (authors)
The temperature control of crystal optics is critical for ultrahigh resolution applications such as those used in meV-resolved Inelastic Scattering. Due to the low count rate and long acquisition time of these experiments, for 1-meV energy resolution, the absolute temperature stability of the crystal optics must be maintained below 4 mK to ensure the required stability of lattice constant, thereby ensuring the energy stability of the optics. Furthermore, the temperature control with sub-mK precision enables setting the absolute temperature of individual crystal, making it possible to align the reflection energy of each crystal’s rocking curve in sub-meV resolution thereby maximizing the combined efficiency of the crystal optics. In this contribution, we report the details of an EPICS control system using PT1000 sensors, Keithley 3706A 7.5 digits sensor scanner, and Wiener MPOD LV power supply for the analyzer crystals of the Inelastic X-ray Scattering (IXS) beamline 10-ID at NSLS-II**. We were able to achieve absolute temperature stability below 1 mK and sub-meV energy alignment for several asymmetrically cut analyzer crystals. The EPICS ePID record was used for the control of the power supplies based on the PT1000 sensor input that was read with 7.5 digits accuracy from the Keithley 3706A scanner. The system enhances the performance of the meV-resolved IXS spectrometer with currently a 1.4 meV total energy resolution and unprecedented spectral sharpness for studies of atomic dynamics in a broad range of materials.
Previous work has indicated that inelastic grains undergoing homogeneous cooling may be unstable, giving rise to the formation of velocity vortices and particle clusters for sufficiently large systems. Such instabilities are observed in industrial coal and biomass gasifiers and are known to influence gas-solid contact area, mixing dynamics, and heat/mass transfer rates. However, the driving mechanisms that lead to vortices and clusters are not well understood. Discrete-particle simulations provide a well-established method for understanding such mechanisms but are not a feasible technique for predicting the behavior of large-scale systems. Kinetic-theory-based continuum models offer an effective means of describing such flows, and instabilities present a stringent test of such models due to the transient, three-dimensional nature of instabilities and the large range of time and length scales over which these mechanisms occur.This work begins with the study, via a combination of continuum models and discrete- particle simulations, of a relatively simple flow and includes additional complexities in a stepwise manner to assess various driving mechanisms. Comparisons with discrete-particle simulations, which offer detailed, well-established (but computationally limited) descriptions of particle flows, indicate the ability of continuum models to accurately incorporate each mechanism. Specifically, the critical length scale for velocity vortices and/or particle clusters are studied via direct numerical simulation, molecular dynamics simulations, linear stability analyses of the continuum model, and transient simulations of the continuum model in a range of flow complexities, including moderate dissipation and particle concentration, frictional particles collisions, high gradients, and gas-solid flows. Strong agreement between kinetic-theory-based continuum models and discrete-particle simulations is found for a range for conditions. Furthermore, discrete-particle simulations offer insight into the relative importance of various mechanisms for instability in both granular and gas-solid flows. Normal inelasticity and friction are both found to be important in granular flows. Particle collisions, thermal drag, and mean drag are found to be important in gas-solid flows.
The Main Injector Particle Production (MIPP) experiment is a fixed target hadron production experiment at Fermilab. MIPP is a high acceptance spectrometer which provides excellent charged particle identification using Time Projection Chamber (TPC), Time of Flight (ToF), multicell Cherenkov (Ckov), ring imaging Cherenkov (RICH) detectors, and Calorimeter for neutrons. The MIPP experiment is designed to measure particle production in interactions of 120 GeV/c primary protons from the Main Injector and secondary beams of $\pi^{\pm}, \rm{K}^{\pm}$, p and $\bar{\rm{p}}$ from 5 to 90 GeV/c on nuclear targets which include H, Be, C, Bi and U, and a dedicated run with the NuMI target. The goal of the experiment is to measure hadron production cross sections or yields using these beams and targets. These hadronic interaction data can have a direct impact on the detailed understanding of the neutrino fluxes of several accelerator-based neutrino experiments like MINOS, MINER$\nu$A, NO$\ nu$A and ELBNF. In this thesis, we present inelastic cross section measurements for 58 and 85 GeV/c p+p interactions, and 58 and 120 GeV/c p+Be, p+C and p+Bi interactions. A new method is described to account for the low multiplicity inefficiencies in the interaction trigger using KNO scaling. Inelastic cross sections as a function of multiplicity are also presented for p+H and p+C interactions at different beam momenta. The MIPP data are compared with the Monte Carlo predictions and existing data. We also describe an algorithm to identify charged particles ($\pi^{\pm}/\rm{p}/\bar{\rm{p}}$ etc.), and present the inclusive charged pion production cross sections in bins of true momentum using 58 and 120 GeV/c p+C interactions for both the data and Monte Carlo.
This report describes a modification to the elastic-perfectly plastic (EPP) strain limits design method to account for cyclic softening in Gr. 91 steel. The report demonstrates that the unmodified EPP strain limits method described in current ASME code case is not conservative for materials with substantial cyclic softening behavior like Gr. 91 steel. However, the EPP strain limits method can be modified to be conservative for softening materials by using softened isochronous stress-strain curves in place of the standard curves developed from unsoftened creep experiments. The report provides softened curves derived from inelastic material simulations and factors describing the transformation of unsoftened curves to a softened state. Furthermore, the report outlines a method for deriving these factors directly from creep/fatigue tests. If the material softening saturates the proposed EPP strain limits method can be further simplified, providing a methodology based on temperature-dependent softening factors that could be implemented in an ASME code case allowing the use of the EPP strain limits method with Gr. 91. Finally, the report demonstrates the conservatism of the modified method when applied to inelastic simulation results and two bar experiments.
The Qweak project is seeking to find new physics beyond the Standard Model. It is aimed to measure the weak charge of the proton, which has never been measured, at 4% precision at low momentum transfer. The experiment is performed by scattering electrons from protons and exploiting parity violation in the weak interaction at low four-momentum transfer. In this experiment, two measurements were considered: which are elastic and inelastic. The elastic is to measure the proton's weak charge. In addition, the inelastic asymmetry measurement, which will extract the low energy constant d. That measurement works in the neutral current sector of the weak interaction. Qweak measures the asymmetry in the N → Δ; transition. Because the elastic radiative tail gives a dominant contribution to the uncertainty to the N → Δ; asymmetries, this thesis will discuss the radiative correction. In addition, this thesis will describe in details the extensive simulations preformed to determine the impact of all simulated background processes on extracting the PV N → Δ; asymmetries. In the process of verifying the validity of these background fractions, we determined the best value of a quantity measured during the Qweak experiment: the beam normal single spin asymmetry, Bn, in the N → Δ; transition.
Our group studies inelastic and reactive collisions of small molecules, focusing on radicals important in combustion environments. The goal is the better understanding of kinetic processes that may be difficult to access experimentally. An essential component is the accurate determination and fitting of potential energy surfaces (PESs). After fitting the ab initio points to obtain global PESs, we treat the dynamics using time-independent (close-coupling) methods. Cross sections and rate constants for collisions of are determined with our Hibridon program suite . We have studied energy transfer (rotationally, vibrationally, and/or electronically inelastic) in small hydrocarbon radicals (CH 2 and CH 3 ) and the CN radical. We have made a comparison with experimental measurements of relevant rate constants for collisions of these radicals. Also, we have calculated accurate transport properties using state-of-the-art PESs and to investigate the sensitivity to these parameters in 1-dimensional flame simulations. Of particular interest are collision pairs involving the light H atom.
Techniques were developed to measure properties at sub-grain scales using depleted Uranium Oxide (d-UO2) samples heat-treated to obtain different grain sizes and oxygen stoichiometries, through three main tasks: 1) sample processing and characterization, 2) microscale and conventional testing and 3) modeling. Grain size and crystallography were characterized using Scanning Electron Microscopy (SEM), in conjunction with Electron Backscattering Diffraction (EBSD) and Electron Channeling Contrast Imaging (ECCI). Grains were then carefully selected based on their crystallographic orientations to perform ex-situ micromechanical tests with samples machined via Focused Ion Beam (FIB), with emphasis on micro-cantilever bending. These experiments were performed under controlled atmospheres, to insure stoichiometry control, at temperatures up to 700 °C and allowed measurements involving elastic (effective Young’s modulus), plastic (critical resolved shear stresses) and creep (creep strain rates) behavior. Conventional compression experiments were performed simultaneously to compare with the ex-situ measurements and study potential size effects. Modeling was implemented using anisotropic elasticity and inelastic constitutive relations for plasticity and creep based on kinematics and kinetics of dislocation glide that account for the effects of crystal orientation, and stress. The models will be calibrated and validated using the experimental data. This project provided insight on correlations among stoichiometry, crystallography and mechanical behavior in advanced oxide fuels, provided valuable experimental data to validate and calibrate mesoscale fuel performance codes and also a framework to measure sub-grain scale mechanical properties that should be suitable for use with irradiated samples due to small volumes required. The goals and metrics of the ongoing study of thermo-mechanical behavior in depleted uranium dioxide (d-UO 2 ) outlined in this project have been concluded successfully, resulting in: 1) the successful fabrication, processing, and characterization of large-grained samples with various orientations (up to and including single crystals) having stoichiometric and hyper-stoichiometric O/U ratios; 2) formulation, calibration, and validation of a crystal plasticity constitutive model to describe the creep deformation of UO 2 at the sub-grain length scale (single crystal level) at intermediate temperatures; 3) the successful calibration of a crystal plasticity constitutive model to describe the elasto-plastic deformation of microcantilever beams, also at moderate temperatures. Samples were prepared from natural uranium oxide powder of production-quality provided by Areva. The powder was pressed in a die to a pressure of 100 MPa to produce green pellets with no sintering aids, lubricants, or any other additives. The green pellets were then heated up to 1700 °C under ultra-high purity argon atmosphere (~1 ppm O2). The atmosphere was then changed to 79% Argon, 21% O 2 and the temperature was held at 1700 °C for 2 hours to sinter the pellets under oxidative conditions [1] that are known to increase grain growth kinetics in UO 2 [2]. Samples were then cooled down under Ar-4%H2 atmosphere to reduce the samples back to stoichiometric UO 2 . For macro-scale procedures, testing of UO 2 samples with large grains was performed at 1200 °C using a modified load frame capable of applying dead-weight loads to ensure constant stress conditions, while displacement of the sample produced by the applied load was measured with high precision micrometers to obtain strains. Stress steps were used during testing and the strains were monitored to measured creep strain rates under steady state for each level of stress used, so that stress exponents could be obtained. The results of the mechanical testing, along with sample geometry and crystal orientation of the grains in the samples, as well as post-test sample characterization were used to formulate a viscoplastic model to account for steady state (stage II) creep behavior, along with basic assumptions from crystal plasticity and kinematic constraints due to testing fixtures. In the micro-scale, testing of microcantilever beams at temperatures ranging from 25 to 570 °C was performed in-situ with a scanning electron microscope with a special attachment to apply load and measure displacement while the samples were at temperature. The load-displacement curves obtained showed linear behavior before fracture for all temperatures attempted except 570 °C, where clear deviations from non-linearity were observed before fracture. These deviations were consistently observed for all samples tested for a given orientation. A viscoplastic model was used to account for the presence of inelastic strain, along with basic assumptions from crystal plasticity and beam theory. These models were kept as simple as possible, and results from tests performed in a set of samples with a given crystal orientation were used to calibrate the material constants for the model, while results from a different sample set were then used for validation, thus satisfying the conditions of all main tasks within the parameters of this project. Details of these efforts are outlined in this report.
X-ray Thomson scattering (XRTS) is an experimental technique that can directly probe the complicated physics of warm and hot dense matter. The Compton-shifted profile of inelastically-scattered x-rays reflects the electron velocity distribution and can measure temperature and density. The ratio of elastically versus inelastically scattered x-rays is related to the number of tightly bound versus free electrons, and thus reflects the ionization state of the sample. Recent experiments highlight uncertainties in the ionization models often used to describe warm dense matter. More experimental validation of these models is needed to understand the onset of ionization in dense plasmas.
Over the last several Evaluated Nuclear Data File (ENDF)/B releases, many updates have been introduced to the nuclear data libraries that can significantly change computational results. These changes can be particularly important for advanced reactor concepts, which are not as thoroughly investigated as those used in light water reactor (LWR) systems. Therefore, a performance assessment of the ENDF/B libraries using the SCALE 6.2.3 code package was conducted under the auspices of the U.S., Department of Energy Office of Nuclear Energy (DOE-NE). The following systems relevant to the advanced reactor community were chosen as subject for investigation in this report: sodium cooled fast reactor systems, graphite moderated high temperature gas-cooled reactors, and several molten salt reactor models. At first, the similarity of the integral performance of the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries caused by compensating errors in important differential data was demonstrated. In calculations for a shipping container of high-assay low enriched uranium, the data of individual isotopes were systematically swapped between the two ENDF/B libraries. An eigenvalue difference of up to 450 pcm was found due to the use of 235 U and 238 U from one library and 1 H and 16 O from the other library. This clearly demonstrates a cross-correlation between reaction data sets of different isotopes within a library that should be reported in the evaluations. For sodium-cooled fast reactor (SFR) fuel assemblies, significant eigenvalue (k inf ) differences (200– 450 pcm) were found between calculations using the 2011 ENDF/B-VII.1 data and 2018 ENDF/B-VIII.0 data. These differences were mainly caused by updates of 238 U and 239 Pu neutron cross sections. Results from the multigroup (MG) calculations further revealed that the group structure of the applied MG library strongly influences the MG bias due to the importance of an appropriate energy resolution of the resonances in higher energy ranges. The application of previously used MG libraries requires new verification for advanced reactor simulation in comparisons with reference continuous-energy calculations. Nuclear data uncertainty analyses of these SFR systems resulted in eigenvalue uncertainties between 1,400 and 1,800 pcm, which are three to four times higher than corresponding uncertainties in light water reactor (LWR) systems. The main contributor to this uncertainty was found to be inelastic scattering on 238 U, which shows an uncertainty of up to 50% in the fast energy range. Other important contributors are the scattering reactions of 56 Fe and 23 Na which have so far not appeared in LWR analysis. Since the largest contribution to the eigenvalue uncertainty of SFR systems is coming from scattering reactions, it is expected that uncertainties in the angular scattering distributions also have a significant impact. Those uncertainties are, however, not available for the majority of nuclides, and the capabilities to determine sensitivities to this data are currently underdeveloped. For two graphite moderated high temperature gas-cooled reactor (HTGR) benchmarks, only the ENDF/BVII. 1 calculations resulted in consistent eigenvalues (k eff ) with the corresponding experimental measurements. Eigenvalue differences of about 1,000 pcm were found between the 2006 ENDF/B-VII.0 data and the 2011 ENDF/B-VII.1 data, primarily due to an updated carbon capture cross section in the thermal energy range in the later ENDF/B release. The ENDF/B-VIII.0 eigenvalue was larger than the ENDF/B-VII.1 eigenvalue by about 300 pcm mainly due to updates in multiple 235 U neutron cross sections, with offsetting updates in the 238 U cross sections. With ENDF/B-VIII.0, graphite can be modeled as perfect crystal or with two different porosities. The choice of the graphite evaluation can have a significant influence of the eigenvalue. For a HTGR benchmark, maximum eigenvalue differences as high as 650 pcm due to different porosities were observed. More detailed studies of the impact of the porosity are necessary while considering that the graphite porosity can vary between the used materials and changes as a function of neutron fluence. xii The uncertainty of the HTGR eigenvalues due to nuclear data uncertainties was found between 500 and 600 pcm, with the top contributor being the neutron multiplicity of 235 U. A gap in the form of missing uncertainties in graphite thermal scattering data, that might have significant impact on HTGR reactors, was identified. This gap is being investigated by the DOE-NE Nuclear Data and Benchmarking Program, as documented in a separate report The fast spectrum molten salt reactor calculations revealed similar comparisons as for the SFR assemblies. The eigenvalue uncertainties were significantly influenced by inelastic scattering on 238 U, while the impact of angular scattering distributions is unknown. For the graphite and zirconium hydride moderated system, eigenvalue uncertainties of up to 700 pcm were observed, requiring in depth analyses of the cross sections and corresponding uncertainties of the included materials. In general, it was observed that the flux spectra of the various molten salt reactor systems show significant differences, even between a fresh state and a depleted state. The choice of the energy group structure of MG calculations is therefore highly relevant; the applicability of previously used MG libraries needs to be verified. Due the unavailability of pin power measurements of advanced reactor systems, pin power calculations with both ENDF/B-VII.1 and VIII.0 data were performed for a light water reactor and compared to corresponding measurements. Both calculations show good agreement with the measurements. A comparison between the ENDF/B-VII.1 and VIII.0 results revealed small differences, mostly in the range of about 0.2%.
The Q weak Experiment ran for two and a half years at the Thomas Jefferson National Accelerator Facility in pursuit of Q$^p_w$, the neutral weak analog to the electric charge of the proton. Q weak measured the parity-violating asymmetry in elastic electron-proton scattering at an extreme forward angle (Q 2 = 0.0249 (GeV/c) 2 ). From the data gathered via the 1.16 GeV 180 µA longitudinally polarized electron beam scattering off the unpolarized photons in the liquid hydrogen target, a value of Q$^p_w$ (PVES) = 0.0719 ± 0.0045 was determined. The Standard Model has a definite prediction of Q$^p_w$ (SM) = 0.0708 ± 0.0003, consistent with the value determined by Q weak which sets a limit on possible new physics up to 7.5 TeV. The theory behind the main measurement of the Q weak Experiment is described in this document, along with the apparatus that made the measurement possible. Understanding the kinematics of the apparatus was a vital component to Q weak 's final measurement. An in-depth explanation of the tracking system responsible for benchmarking the momentum transfer and scattering angle simulations is included. The unblinded analysis of Q weak 's final result is outlined, as are additional physics results related to the N → Δ transition. During April 2012, an opportunity was seized to take ancillary data on the inelastic N → Δ transition at a different beam energy(877 MeV) than the nominal Q weak data. This data, combined with the inelastic data taken at nominal beam energy and a previous measurement, determined a constraint on d Δ , a low energy constant related to hadronic parity violation, of (3.8 ± 14.7) g π . It also resulted in a measurement of the beam-normal single-spin asymmetry of the N → Δ transition of 149 ± 3 (stat) ± 72 (syst) ppm. This document includes both the longitudinal and transverse analysis of the 877 MeV data.
This thesis covers two subjects, brie y outlined below, of hadronic physics. The proton radius puzzle and the muon anomalous magnetic moment discrepancy point to possible signs of lepton-universality violation. We introduce the context and background necessary to understand these two problems. And we analyze two suitable beyond-the-standard-model solutions, one vector-based and one scalar-based, which simultaneously solves both issues. Furthermore, we demonstrate that the parameter space of the respective solutions can be chosen so as to not be completely forbidden by considered experimental constraints. Specifically, we show that certain violations, analyzed by other authors in the context of similar solutions to those proposed, are ameliorated if one adds to the proposed solutions particular criteria required for renormalizability. We also explore the field of lattice perturbation theory in the context of deep inelastic scattering. In particular, we study the connection between lattice-based Euclidean parton distributions and their experiment-based Minkowskian equivalents. These parton distributions are extracted from deep inelastic scattering processes in which leptons are scattered of hadronic targets at energy levels high enough to resolve the internal hadronic structure. To synchronize experiment with theory, one must find the proper matching criteria to transform approximations of parton distributions, known as quasidistributions, as they are resolved on a computer (i.e. the lattice), back to the full parton distributions extracted from experiment. We provide the necessary background to understand context of the above and show the procedures necessary to match experiment with the lattice. We demonstrate how one must tread carefully to properly perform such a matching.
The liquid argon time projection chamber (LArTPC) is the detection technology chosen for several Fermilab-based neutrino experiments. This technology will be used in studies of neutrino cross-sections and oscillations, neutrinos from supernovae as well as a variety of studies of beyond the Standard Model physics. This thesis explores the use of these detectors to study MeV-scale activity. MeV-scale electrons arising from Compton scatters of deexcitation photons and photons from inelastic neutron scattering in neutrino-nucleus interactions are reconstructed using novel methods presented here. This work represents the first demonstration of MeV-scale physics capabilities in a LArTPC neutrino experiment as well as the first observation of neutrino-produced photons from nuclear de-excitation and inelastic neutron scattering. A search for millicharged particles, postulated by theories of beyond the standard model physics, is also performed using data from a LArTPC and the low-ener gy reconstruction techniques developed. The results set world-leading bounds on the parameter space of millicharged particles. The work in this thesis demonstrates that studies of MeV-scale activity and new physics are possible with LArTPC technology and provides the foundation for future LArTPC studies of low energy neutrinos and new physics.
A recent journal article by J.A. Gomez et al. measured gamma-ray die-off curves of three subcritical static highly-enriched uranium (HEU) assemblies driven by an external neutron source with a new detector system. This new detector system was developed for being used in dynamically driven subcritical assemblies as part of the Neutron Diagnosed Subcritical Experiments (NDSE) program. Simulations of these die-off curves with various nuclear data and comparison to experimental data indicated that a decrease of the ENDF/B-VII.1 235 U(n,inl) cross section by a factor 0.8 and 0.85 for ENDF/B-VIII.0 would lead to better predictions of experimental data. Here, we test the proposed changes with another type of measurement response, namely neutron-leakage spectra emitted in LLNL pulsed sphere measurements. These spheres were pulsed by 14-MeV neutrons produced via the D+T reaction in their center. The proposed changes in nuclear data have a distinctly smaller impact on predicting pulsed-sphere neutron-leakage spectra than for the die-off curves; they lead to a worsened prediction of the inelastic valley of LLNL pulsed-sphere neutron spectra indicating that the proposed change could constitute a compensating error. Changes in the inelastic angular distributions along with the cross section might be worthwhile to study
The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.
The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.