Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Experimental and theoretical study of three interacting, closely-spaced, sharp-edged 60 deg delta wings at low speeds

Wind tunnel tests were conducted to determine the subsonic longitudinal aerodynamic characteristics of lifting configuration consisting of a 60 deg delta main wing with two smaller 60 deg delta wings (called sub-wings) attached underneath. The test was designed to determine the effects on lift, drag, and pitching moment due to various placement of the subwings in relation to the main wing. Test results indicate the increasing vertical separation between the main wing and the sub-wings produced the most significant results; a 23.1% increase in maximum lift coefficient, a reduction in drag coefficient at high lift coefficients, and an increase in longitudinal stability. Lateral separation of the sub-wings produced no significant changes. Placement of the sub-wings rearward increases the initial lift curve slope and maximum lift coefficient and also increase the longitudinal stability. Results of a computer study using a vortex lattice code supported the experimental conclusions.

Faery, H. F., Jr.↗
DOE OSTI2022

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Multigrid methods for numerical simulation of laminar diffusion flames

This paper documents the result of a computational study of multigrid methods for numerical simulation of 2D diffusion flames. The focus is on a simplified combustion model, which is assumed to be a single step, infinitely fast and irreversible chemical reaction with five species (C3H8, O2, N2, CO2 and H2O). A fully-implicit second-order hybrid scheme is developed on a staggered grid, which is stretched in the streamwise coordinate direction. A full approximation multigrid scheme (FAS) based on line distributive relaxation is developed as a fast solver for the algebraic equations arising at each time step. Convergence of the process for the simplified model problem is more than two-orders of magnitude faster than other iterative methods, and the computational results show good grid convergence, with second-order accuracy, as well as qualitatively agreement with the results of other researchers.

Liu, C.↗
DOE OSTI2024

Synthesis, Isolation, and Study of Heterobimetallic Uranyl Crown Ether Complexes

Although crown ethers can selectively bind many metal cations, little is known regarding the solution properties of crown ether complexes of the uranyl dication, UO 2 2+ . Here, in this work, the synthesis and characterization of isolable complexes in which the uranyl dication is bound in an 18-crown-6-like moiety are reported. A tailored macrocyclic ligand, templated with a Pt(II) center, captures UO 2 2+ in the crown moiety, as demonstrated by results from single-crystal X-ray diffraction analysis. The U(V) oxidation state becomes accessible at a quite positive potential (E 1/2 ) of –0.18 V vs Fc +/0 upon complexation, representing the most positive U VI /U V potential yet reported for the UO 2 n+ core. Isolation and characterization of the U(V) form of the crown complex are also reported here; there are no prior reports of reduced uranyl crown ether complexes, but U(V) is clearly stabilized by crown chelation. Joint computational studies show that the electronic structure of the U(V) form results in significant weakening of U–O oxo bonding despite the quite positive reduction potential at which this species can be accessed, underscoring that crown-ligated uranyl species could demonstrate unique reactivity under only modestly reducing conditions.

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

Tuning the Functionalities of Porous Liquids for Emergent Gas-Capture Properties

Type 3 Porous Liquids (PLs) are a class of materials with the potential to revolutionize gas capture, storage, and utilization. These PLs are formed by suspending sorbent nanoparticles (e.g., metal−organic frameworks) in sterically excluded solvents, creating permanent porosity for gas capture in a processable, low-viscosity phase. Herein, a computational study revealed sorbent surface functionalization strategies to enhance CO 2 sorption, and the molecular structural signatures underpinning the enhancements in gas uptake. PLs composed of a ZIF-8 surface functionalized with 3-amino-1,2,4-triazole (Atz) in glyceryl triacetate were targeted for emergent CO 2 capture, exceeding that of the unfunctionalized ZIF-8 PL. ZIF-8 was surface functionalized with Atz at various surface coverage fractions (f), and classical molecular dynamics simulations predicted an increase in CO 2 sorption capacity with increasing f, up to f = 0.75. Additionally, detailed structural analyses revealed that solvent orientational order, derived from the solvent triplet-angle distribution, can identify the gas-capture potential of a PL without requiring computationally expensive direct modeling of the CO 2 sorption. Combined with experimental validation, initial computational screening of PL compositions promises to accelerate the discovery of PL compositions for novel gas separation materials platforms.

gas capture↗
NASA NTRS2004

Evaluation of CO2, N2 and He as Fire Suppression Agents in Microgravity

The U.S. modules of the International Space Station use gaseous CO2 as the fire extinguishing agent. This was selected as a result of extensive experience with CO2 as a fire suppressant in terrestrial applications, trade studies on various suppressants, and experiments. The selection of fire suppressants and suppression strategies for NASA s Lunar and Martian exploration missions will be based on the same studies and normal-gravity data unless reduced gravity fire suppression data is obtained. In this study, the suppressant agent concentrations required to extinguish a flame in low velocity convective flows within the 20-sec of low gravity on the KC-135 aircraft were investigated. Suppressant gas mixtures of CO2, N2, and He with the balance being oxygen/nitrogen mixtures with either 21% or 25% O2 were used to suppress flames on a 19-mm diameter PMMA cylinder in reduced gravity. For each of the suppressant mixtures, limiting concentrations were established that would extinguish the flame at any velocity. Similarly, concentrations were established that would not extinguish the flame. The limiting concentrations were generally consistent with previous studies but did suggest that geometry had an effect on the limiting conditions. Between the extinction and non-extinction limits, the suppression characteristics depended on the extinguishing agent, flow velocity, and O2 concentration. The limiting velocity data from the CO2, He, and N2 suppressants were well correlated using an effective mixture enthalpy per mole of O2, indicating that all act via O2 displacement and cooling mechanisms. In reduced gravity, the agent concentration required to suppress the flames increased as the velocity increased, up to approximately 10 cm/s (the maximum velocity evaluated in this experiment). The effective enthalpy required to extinguish flames at velocities of 10 cm/s is approximately the same as the concentrations in normal gravity. A computational study is underway to further evaluate these findings.

Ruff, Gary A.↗
NASA NTRS2001

Thermal Effects on the Bearing Behavior of Composite Joints

Thermal effects on the pin-bearing behavior of an IM7/PET15 composite laminate are studied comprehensively. A hypothesis presents factors influencing a change in pin-bearing strength with a change in temperature for a given joint design. The factors include the change in the state of residual cure stress, the material properties, and the fastener fit with a change in temperature. Experiments are conducted to determine necessary lamina and laminate material property data for the IM7/PET15 being utilized in this study. Lamina material properties are determined between the assumed stress free temperature of 460 F down to -200 F. Laminate strength properties are determined for several lay-ups at the operating temperatures of 350 F, 70 F, and -200 F. A three-dimensional finite element analysis model of a composite laminate subject to compressive loading is developed. Both the resin rich layer located between lamina and the thermal residual stresses present in the laminate due to curing are determined to influence the state of stress significantly. Pin-bearing tests of several lay-ups were conducted to develop an understanding on the effect of temperature changes on the pin-bearing behavior of the material. A computational study investigating the factors influencing pin-bearing strength was performed. A finite element model was developed and used to determine the residual thermal cure stresses in the laminate containing a hole. Very high interlaminar stress concentrations were observed two elements away from the hole boundary at all three operating temperatures. The pin-bearing problem was modeled assuming a rigid frictionless pin and restraining only radial displacements at the hole boundary. A uniform negative pressure load was then applied to the straight end of the model. A solution, where thermal residual stresses were combined with the state of stress due to pin-bearing loads was evaluated. The presence of thermal residual stresses intensified the interlaminar stresses predicted at the hole boundary in the pin-bearing problem. This dissertation shows that changes in material properties drives pin-bearing strength degradation with increasing temperature.

Walker, Sandra Polesky↗
NASA NTRS2009

Experimental Studies of the Aerothermal Characteristics of the Project Orion CEV heat Shield in High Speed Transitional and Turbulent Flows

An experimental program has been completed by CUBRC exploring laminar, transitional, and turbulent flows over a 7.0% scale model of the Project ORION CEV geometry. This program was executed primarily to answer questions concerning the increase in heat transfer on the windward, or "hot shoulder" of the CEV heat shield from laminar to turbulent flow. To answer these questions CUBRC constructed and instrumented a 14.0 inch diameter Project ORION CEV model and ran a range of Reynolds numbers based on diameter from 1.0 to over 40 million at a Mach number of 8.0. These Reynolds numbers were selected to cover laminar to turbulent heating data on the "hot shoulder". Data obtained during these runs will be used to guide design decisions as they apply to heat shield thickness and extent. Several experiments at higher enthalpies were achieved to obtain data for code validation with real gas effects and transition. CUBRC also performed computation studies of these experiments to aid in the data reduction process and study turbulence modeling.

Wadhams, T.P.↗
NASA NTRS1996

An experimental and computational investigation of a translating throat single expansion-ramp nozzle

A translating throat single expansion-ramp nozzle (SERN) concept was designed to improve the off-design performance of a SERN with a large, fixed expansion ratio. The concept of translating the nozzle throat provides the SERN with a variable expansion ratio. An experimental and computational study was conducted to predict and verify the internal performance of this concept. Three nozzles with expansion ratios designed for low, intermediate, and high Mach number operating conditions were tested in the Jet-Exit Test Facility at the NASA Langley Research Center. Each nozzle was tested with a concave and a convex geometric expansion ramp surface design. Internal nozzle performance, paint-oil flow and focusing Schlieren flow visualization were obtained for nozzle pressure ratios (NPR's) up to 13. The Navier-Stokes code, PAB3D, with a k-epsilon turbulence model was utilized to verify experimental results at selected NPRs and to predict the performance at conditions unattainable in the test facility. Two-dimensional simulations were computed with near static free-stream conditions and at nozzle pressure ratios of 5, 9, and 13 for the concave ramp, low Mach number configuration and at the design NPR of 102 for the concave ramp, high Mach number configuration. Remarkable similarities between predicted and experimental flow characteristics, as well as performance quantities, were obtained.

Deere, Karen A.↗
NASA NTRS1996

An Experimental and Computational Investigation of a Translating Throat Single Expansion-Ramp Nozzle

A translating throat single expansion-ramp nozzle (SERN) concept was designed to improve the off-design performance of a SERN with a large, fixed expansion ratio. The concept of translating the nozzle throat provides the SERN with a variable expansion ratio. An experimental and computational study was conducted to predict and verify the internal performance of this concept. Three nozzles with expansion ratios designed for low, intermediate, and high Mach number operating conditions were tested in the Jet-Exit Test Facility at the NASA Langley Research Center. Each nozzle was tested with a concave and a convex geometric expansion ramp surface design. Internal nozzle performance, paint-oil flow and focusing Schlieren flow visualization were obtained for nozzle pressure ratios (NPR's) up to 13. The Navier-Stokes code, PAB3D, with a k-epsilon turbulence model was utilized to verify experimental results at selected NPR's and to predict the performance at conditions unattainable in the test facility. Two-dimensional simulations were computed with near static free-stream conditions and at nozzle pressure ratios of 5, 9, and 13 for the concave ramp, low Mach number configuration and at the design NPR of 102 for the concave ramp, high Mach number configuration. Remarkable similarities between predicted and experimental flow characteristics, as well as performance quantities, were obtained.

Deere, Karen A.↗
NASA NTRS1987

Mach number effects on conical surface features of swept shock boundary-layer interactions

A joint experimental and computational study is made of the shock-wave turbulent boundary-layer interaction generated by sharp fins, with emphasis on Mach-number effects. The Mach number range is from 2 to 4 and the unit Reynolds number is from 50 to 80 million per meter. Fin angles are varied from 4 to 22 deg. Surface-flow patterns are obtained using a color surface-flow-visualization technique. The results show that the upstream-influence response in the conical far-field region is a function of the freestream Mach number and the shock strength. A new interpretation of the behavior of the upstream influence with changes of the inviscid shock angle is given. Agreement between the experimental and the computed upstream-influence lines becomes poorer for stronger interactions, with the computations underpredicting the upstream-influence line.

Lu, F. K.↗
DOE OSTI2021

Nickel–Sulfonate Mode of Substrate Binding for Forward and Reverse Reactions of Methyl-SCoM Reductase Suggest a Radical Mechanism Involving Long-Range Electron Transfer

Methyl-coenzyme M reductase (MCR) catalyzes both the synthesis and the anaerobic oxidation of methane (AOM). Its catalytic site contains Ni at the core of cofactor F 430 . The Ni ion, in its low-valent Ni(I) state, lights the fuse leading to homolysis of the C–S bond of methyl-coenzyme M (methyl-SCoM) to generate a methyl radical, which abstracts a hydrogen atom from coenzyme B (HSCoB) to generate methane and the mixed disulfide CoMSSCoB. Direct reversal of this reaction activates methane to initiate anaerobic methane oxidation. On the basis of the crystal structures, which reveal a Ni–thiol interaction between Ni(II)–MCR and inhibitor CoMSH, a Ni(I)–thioether complex with substrate methyl-SCoM has been transposed to canonical MCR mechanisms. Similarly, a Ni(I)–disulfide with CoMSSCoB is proposed for the reverse reaction. However, this Ni(I)–sulfur interaction poses a conundrum for the proposed hydrogen-atom abstraction reaction because the >6 Å distance between the thiol group of SCoB and the thiol of SCoM observed in the structures appears to be too long for such a reaction. The spectroscopic, kinetic, structural, and computational studies described here establish that both methyl-SCoM and CoMSSCoB bind to the active Ni(I) state of MCR through their sulfonate groups, forming a hexacoordinate Ni(I)–N/O complex, not Ni(I)–S. These studies rule out direct Ni(I)–sulfur interactions in both substrate-bound states. As a solution to the mechanistic conundrum, we propose that both the forward and the reverse MCR reactions emanate through long-range electron transfer from the Ni(I)–sulfonate complexes with methyl-SCoM and CoMSSCoB, respectively.

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

A High‐Rate Li–CO 2 Battery Enabled by 2D Medium‐Entropy Catalyst

Abstract Lithium‐air batteries based on CO 2 reactant (Li–CO 2 ) have recently been of interest because it has been found that reversible Li/CO 2 electrochemistry is feasible. In this study, a new medium‐entropy cathode catalyst, (NbTa) 0.5 BiS 3 , that enables the reversible electrochemistry to operate at high rates is presented. This medium entropy cathode catalyst is combined with an ionic liquid‐based electrolyte blend to give a Li–CO 2 battery that operates at high current density of 5000 mA g −1 and capacity of 5000 mAh g −1 for up to 125 cycles, far exceeding reported values in the literature for this type of battery. The higher rate performance is believed to be due to the greater stability of the multi‐element (NbTa) 0.5 BiS 3 catalyst because of its higher entropy compared to previously used catalysts with a smaller number of elements with lower entropies. Evidence for this comes from computational studies giving very low surface energies (high surface stability) for (NbTa) 0.5 BiS 3 and transmission electron microscopystudies showing the structure being retained after cycling. In addition, the calculations indicate that Nb‐terminated surface promotes Li–CO 2 electrochemistry resulting in Li 2 CO 3 and carbon formation, consistent with the products found in the cell. These results open new direction to design and develop high‐performance Li–CO 2 batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Benchmarking nitrous oxide adsorption and activation in metal–organic frameworks bearing coordinatively unsaturated metal centers

Anthropogenic emissions of N 2 O, the third most abundant greenhouse gas after CO 2 and CH 4 , are contributing to global climate change. Although metal–organic frameworks (MOFs) have been widely studied as adsorbents for CO 2 and CH 4 , less effort has focused on the use of MOFs to remove N 2 O from emission streams or from air. Further, N 2 O activation would enable its use as an inexpensive oxidant for fine chemical synthesis. Here, in this paper, we identify features that contribute to strong binding and high uptake of N 2 O at coordinatively unsaturated metal sites in the M 2 Cl 2 (btdd) (M = Mn, Co, Ni, Cu; btdd 2– = bis(1,2,3-triazolo[4,5-b],[4',5'-i])dibenzo[1,4]dioxin) and M 2 (dobdc) (M = Mg, Mn, Fe, Co, Ni, Cu, Zn; dobdc 4– = 2,5-dioxido-1,4-benzenedicarboxylate) series of MOFs. Combined experimental and computational studies suggest that N 2 O adsorption at open-metal-sites is primarily based on electrostatic interactions, rather than π-backbonding, causing MOFs with more Lewis acidic metal centers to be superior N 2 O adsorbents. As a result, Mg 2 (dobdc) demonstrates strong binding and record-setting N 2 O uptake (8.75 mmol g –1 at 1 bar and 298 K). Using density functional theory (DFT) to characterize reactive intermediates and transition states, we demonstrate that N 2 O activation to form a M(IV)–oxo species and N 2 is thermodynamically favorable in Mn 2 (dobdc) and Fe 2 (dobdc) but appears to be kinetically limited in Mn 2 (dobdc). Our work lays a foundation for understanding N 2 O adsorption and activation in MOFs, paving the way for the design of promising next-generation materials for N 2 O capture and utilization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Distinguishing Desirable and Undesirable Reactions in Multicomponent Systems for Redox Activation of the Uranyl Ion

Although it has been established that covalent functionalization of the U–O bonds in the uranyl dication (UO 2 2+ ) generally requires use of strong reductants and electrophiles, little work has examined how interactions between the individual reaction components could affect final outcomes in solution. Here, the patterns of such reactivity have been studied in a UO 2 2+ -containing model system supported by a workhorse pentadentate ligand, 2,2′-[(methylimino)bis(2,1-ethanediylnitrilomethylidyne)]bis-phenol. Oxo activation and functionalization have been tested with (i) electrochemical and chemical reduction, and (ii) coordinating and noncoordinating solvents. In acetonitrile, uranyl reduction was achieved cleanly, but treatment of the reduced species with tris(pentafluorophenyl)borane (BCF) resulted in a mixture of products arising from direct electron transfer to BCF. In dichloromethane (CH 2 Cl 2 ), electrochemical reduction of uranyl was achieved cleanly, but clean chemical reactivity was inaccessible. Despite these challenges, one trinuclear and oxo-deficient uranium-containing product was crystallized from CH 2 Cl 2 solution and characterized; thus, desirable electrophilic reactivity can proceed to some degree in CH 2 Cl 2 with BCF. Computational studies were used to investigate the properties of the trinuclear uranium product and the changes that could be inducible by further reduction. Here, taken together, the reactivity patterns identified here could inform design of improved systems for actinyl oxo functionalization.

cations↗
DOE OSTI2022

Symmetry-Breaking and Self-Sorting in Block Copolymer-Based Multicomponent Nanocomposites

Co-assembly of inorganic nanoparticles (NPs) and nanostructured polymer matrix represents an intricate interplay of enthalpic or entropic forces. Particle size largely affects the phase behavior of the nanocomposite. Theoretical studies indicate that new morphologies would emerge when the particles become comparable to the soft matrix's size, but this has rarely been supported experimentally. By designing a multicomponent blend composed of NPs, block copolymer-based supramolecules, and small molecules, a 3-D ordered lattice beyond the native BCP's morphology was recently reported when the particle is larger than the microdomain of BCP. The blend can accommodate various formulation variables. In this work, when the particle size equals the microdomain size, a symmetry-broken phase appears in a narrow range of particle sizes and compositions, which we named the "train track" structure. In this phase, the NPs aligned into a 3-D hexagonal lattice and packed asymmetrically along the c axis, making the projection of the ac and the bc plane resemble train tracks. Computational studies show that the broken symmetry reduces the polymer chain deformation and stabilizes the metastable hexagonally perforated lamellar morphology. Given the mobility of the multicomponent blend, the system shows a self-sorting behavior: segregating into two macroscopic phases with different nanostructures based on only a few nanometers NP size differences. Smaller NPs form "train track" morphology, while larger NPs form a "simple hexagon" structure, where the NPs take a symmetric hexagonal arrangement. Detailed structural evolution and simulation studies confirm the systematic-wide cooperativity across different components, indicating the strong self-regulation of the multicomponent system.

"train track" morphology↗
DOE OSTI2024

Manipulating Spin–Lattice Coupling in Layered Magnetic Topological Insulator Heterostructure via Interface Engineering

Induced magnetic order in a topological insulator (TI) can be realized either by depositing magnetic adatoms on the surface of a TI or engineering the interface with epitaxial thin film or stacked assembly of 2D van der Waals (vdW) materials. Herein, the observation of spin-phonon coupling in the otherwise non-magnetic TI Bi 2 Te 3 is reported, due to the proximity of FePS 3 (an antiferromagnet (AFM), T N ≈ 120 K), in a vdW heterostructure framework. Temperature-dependent Raman spectroscopic studies reveal deviation from the usual phonon anharmonicity originated from spin-lattice coupling at the Bi 2 Te 3 /FePS 3 interface at/below 60 K in the peak position (self-energy) and linewidth (lifetime) of the characteristic phonon modes of Bi 2 Te 3 (106 and 138 cm –1 ) in the stacked heterostructure. The Ginzburg-Landau (GL) formalism, where the respective phonon frequencies of Bi 2 Te 3 couple to phonons of similar frequencies of FePS 3 in the AFM phase, is adopted to understand the origin of the hybrid magneto-elastic modes. At the same time, the reduction of characteristic TN of FePS 3 from 120 K in isolated flakes to 65 K in the heterostructure, possibly due to the interfacial strain, which leads to smaller Fe-S-Fe bond angles as corroborated by computational studies using density functional theory (DFT). Besides, inserting hexagonal boron nitride within Bi 2 Te 3 /FePS 3 stacking regains the anharmonicity in Bi 2 Te 3 . As a result, controlling interfacial spin-phonon coupling in stacked heterostructure can have potential application in surface code spin logic devices.

36 MATERIALS SCIENCE↗