Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency estimation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Galaxy Clusters Discovered via the Thermal Sunyaev-Zel’dovich Effect in the 500-square-degree SPTpol Survey

We present a catalog of 689 galaxy cluster candidates detected at significance ξ > 4 via their thermal Sunyaev-Zel’dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and Spitzer satellites, to confirm 544 of these candidates as clusters with ~ 94% purity. The sample has an approximately redshift-independent mass threshold at redshift z > 0.25. The confirmed sample spans 1.5 × 10 14 < M 500 c < 9 × 10 14 M ⊙ /h 70 and 0.03 < z ≲ 1.6 in mass and redshift, respectively, with a median mass of 2.5×10 14 M ⊙ /h 70 and median redshift z = 0.7; 21% of the confirmed clusters are at z > 1. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered ξ by a median value of 0.032, or ~ 0.8% of the ξ = 4 threshold value, and ~ 7% of candidates have a predicted contamination greater than Δξ = 1. With the exception of a small number of systems (< 1%), an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample, representing the deepest SZ-selected cluster sample to date, will be a key component in upcoming astrophysical and cosmological analyses of clusters. In addition to the cluster catalog, we also release the millimeter-wave maps and associated data products used to produce this sample. These maps have depths of 5.3 (11.7) µK CMB -arcmin at 150 (95) GHz and an effective angular resolution of 1. '2 (1. '7). The SPTpol products are available at https:// pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

KIC 12602250: A Low-amplitude Double-mode Delta Scuti Star with Amplitude Modulation

We report a detailed light-curve analysis of the Kepler target Kepler Input Catalog (KIC) 12602250. The results show that KIC 12602250 is a low-amplitude radial double-mode δ-Scuti star with amplitude modulation. The Fourier analysis of the long cadence data (i.e., Q0—Q17, spanning 1471 days) reveals that the variations of the light curve are dominated by the strongest mode with frequency F0 = 11.6141 d{sup -1}, suggesting that KIC 12602250 is a δ-Scuti star. The other independent mode F1 = 14.9741 d{sup -1} is newly detected. The amplitude of the light variations of KIC 12602250 is ~0.06 mag, which indicates that this is a low-amplitude δ-Scuti star; but the ratio of F0/F1 is estimated as 0.7756, which is typical of HADS, and a slow amplitude growth is detected in F1 and f {sub 3}, which could be due to stellar evolution, suggesting that KIC 12602250 could be a post-main-sequence δ Scuti that is crossing the instability strip for the first time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Efficient Precoding for Single Carrier Modulation in Multi-User Massive MIMO Networks

Abstract—Frequency domain (FD) multi-user detection (MUD) and precoding has been shown to be an effective means of approaching the theoretical per-user capacity for single carrier modulation (SCM) schemes in massive MIMO scenarios with highly dispersive channels. When a cyclic prefix is added to the SCM waveform, the circulant structures of the resulting convolutional channel matrix allows for relatively simple expressions for the FD detection and precoding. In this paper, we further refine the computationally efficient minimum mean squared error (MMSE) FD-MUD technique in a time-division duplexing (TDD) massive MIMO setup by pre-computing the scale factor needed to produce an unbiased estimate. The resulting scale factor and the matrix inverse that was calculated for the uplink is then reused for the downlink to form a computationally efficient precoding scheme with superior performance compared to zero-forcing (ZF) precoding. Performance analysis is provided and confirmed through simulation results.

5G and Beyond Communications↗
DOE OSTI2021

Thermal-striping analysis methodology for sodium-cooled reactor design

Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

94 GHz microplasma sustained by a photonic crystal

A photonic crystal (PhC) with a single vacancy sustains an argon microplasma at 94 GHz. The PhC consists of a 7 × 13 array of small alumina rods, each 250 μm in radius. The PhC has a forbidden bandgap from 70–105 GHz. The vacancy creates an internal electromagnetic resonance at 94 GHz that intensifies the electric field of an incident wave to 5 × 10 5 V m –1 . Breakdown of low-pressure argon gas (80 mbar) occurs and forms a stable, ellipsoidal plasma inside the crystal. The central core of the microplasma is probed using a 200 μm optical fiber inserted into the PhC. The rotational temperature of the CH molecule is estimated from a numerical fit of the plasma emission and is found to increase from 450–650 K with both pressure (200–1000 mbar) and incident wave power (100 mW–1100 mW). Under these conditions, the electron density increases with both power and pressure from 0.5–1.5 × 10 20 m –3 . Furthermore, this plasma density range brackets the electron plasma frequency at 94 GHz. The results are compared with a similar microplasma operating at 44 GHz.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Robustness of cosmic birefringence measurement against Galactic foreground emission and instrumental systematics

The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates E modes into E modes in the observed CMB signal through an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles, the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. Here, in this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the EB correlation of polarized foreground emission. Here we propose to correct the bias induced by dust EB by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Cobalt-Group 13 Complexes Catalyze CO 2 Hydrogenation via a Co(-I)/Co(I) Redox Cycle

The Co(-I) dihydrogen complexes, [(η 2 -H 2 )CoML]-, where ML is the group 13 metalloligand, N(o-(NCH 2 PiPr 2 )C 6 H 4 ) 3 M, and M is Al, Ga, or In, were previously reported ( J. Am. Chem. Soc. 2017, 139, 6570-6573). In this work, the related Co(-I) end-on dinitrogen adducts, [(N 2 )CoML]-, were isolated and investigated as precatalysts for CO 2 hydrogenation. The Co–Ga catalyst was highly active, achieving 19,200 formate turnovers with an initial turnover frequency of 27,000 h –1 under 34 atm of 1:1 CO 2 /H 2 and using Verkade’s proazaphosphatrane as a base at ambient temperature. The Co–Al catalyst was moderately active, while the Co–In complex was inactive. Hence, tuning the group 13 ion greatly influences the catalytic activity at the Co site. To elucidate the role of the group 13 support, experimental and theoretical mechanistic studies of the Co–Ga and Co–Al catalysts were conducted. The Co(-I) H 2 species are potent hydride donors with estimated thermodynamic hydricities (ΔG° H– ) of 32.0(1) and 37.4(1) kcal/mol in CH 3 CN for M = Al and Ga, respectively. By acting as masked Co(I) dihydrides, the Co(-I) H 2 species operate via an unusual Co(-I)/Co(I) redox cycle. After hydride transfer to CO 2 , the resulting intermediate is the Co(I) hydride complex, HCoML, which was independently synthesized and structurally characterized for M = Al and Ga. The Gibbs free energy for H 2 binding, ΔG° bind (1 atm), to generate (η 2 -H 2 )HCoML was slightly more favorable for HCoGaL (-4.2(1) kcal/mol) than for HCoAlL (-2.7(1) kcal/mol). In the subsequent step, the deprotonation reaction to regenerate the initial catalyst was much more favorable for (η 2 -H 2 )HCoGaL (pK a of 31.4, CH 3 CN) than for (η 2 -H 2 )HCoAlL (pK a of 34.3). Finally, the straightforward substitution of Al with Ga perturbs the energy profile of the catalytic reaction (|ΔΔG° H– | = 5.4 kcal/mol, |ΔΔG° bind | = 1.5 kcal/mol, and |ΔΔG° K a | = 4.0 kcal/mol) and thus provides a thermodynamic rationale for the higher catalytic efficiency of Co–Ga over Co–Al.

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

Effect of Strain on the Resistivity and Thermal Conductivity of High Purity Niobium

High purity niobium (Nb) is a technologically important material for large-scale accelerator and nano-scale quantum computing applications in microwave frequency range. The high thermal conductivity and low resistivity of Nb are critical to the high performance at temperature range of 0.01–4.0 K. The presence of interstitials such as O, N, H, and C act as scattering centers and alter the mean free path, reducing resistivity and thermal conductivity, and contributing significantly to Nb's thermal performance for temperatures of 2.0 K and above. The residual resistivity ratio (RRR), defined as the ratio of the normal state resistivity at 300 K to that at 4.2 K (Nb, T c = 9.2 K), is an accepted direct estimate of the impurity content of fully recrystallized Nb. Complete re-crystallization of Nb is challenging unless very high temperatures are employed, which is often impractical, hence, in practice, dislocation and dislocation structures impact the thermal performance of Nb due to strong phonon scattering contributions. Furthermore, this paper reports on the degradation of thermal conductivity and RRR of high purity Nb large grain, single crystal with fixed impurity, varying strain, and dislocation content levels. Experimental thermal conductivity data fits the Boltzmann transport equation incorporating dislocation density.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Extreme sea levels at different global warming levels

Abstract The Paris agreement focused global climate mitigation policy on limiting global warming to 1.5 or 2 °C above pre-industrial levels. Consequently, projections of hazards and risk are increasingly framed in terms of global warming levels rather than emission scenarios. Here, we use a multimethod approach to describe changes in extreme sea levels driven by changes in mean sea level associated with a wide range of global warming levels, from 1.5 to 5 °C, and for a large number of locations, providing uniform coverage over most of the world’s coastlines. We estimate that by 2100 ~50% of the 7,000+ locations considered will experience the present-day 100-yr extreme-sea-level event at least once a year, even under 1.5 °C of warming, and often well before the end of the century. The tropics appear more sensitive than the Northern high latitudes, where some locations do not see this frequency change even for the highest global warming levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Virtual Synchronous Generator Control of Multi-port Autonomous Reconfigurable Solar Plants (MARS)

Multi-port autonomous reconfigurable solar power plant (MARS) is an integrated concept for integration of solar photovoltaic (PV) and energy storage systems (ESS) to transmission alternating current (ac) grid and high-voltage direct current (HVdc) links. The integrated development incorporates advanced control methods to provide enhanced grid ancillary services such as fast frequency responses and dynamic voltage support. In this paper, the virtual synchronous generator (VSG) control algorithm of MARS is discussed. The proposed VSG-based control enables enhanced synthetic inertial response and flexible frequency response characteristics of hybrid PV and ESS power plants in accordance with grid ancillary service requirements. Multi-port power electronics interface allows controlled emergency power support from MARS to local transmission ac grid and remote grid through the HVdc link. The performance of VSG control is validated using a reduced-order model of MARS in Simulink. Methods for estimating grid ancillary service capabilities of MARS are also discussed.

Pan, Jiuping↗
DOE OSTI2025

Logarithmic Resilience Risk Metrics That Address the Huge Variations in Blackout Cost

Resilience risk metrics must address the customer cost of the largest blackouts of greatest impact. However, there are huge variations in blackout cost in observed distribution utility data that make it impractical to properly estimate the mean large blackout cost and the corresponding risk. These problems are caused by the heavy tail observed in the distribution of customer costs. To solve these problems, we propose resilience metrics that describe large blackout risk using the mean of the logarithm of the cost of large-cost blackouts, the slope index of the heavy tail, and the frequency of large-cost blackouts.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Automating Bayesian inference and design to quantify acoustic particle levitation

Self-propulsion of micro- and nanoparticles powered by ultrasound provides an attractive strategy for the remote manipulation of colloidal matter using biocompatible energy inputs. Quantitative understanding of particle motion and its dependence on size, shape, and composition requires accurate characterization of the acoustic field, which depends sensitively on the experimental setup. Here, we show how automated experiments based on Bayesian inference and design can accurately and efficiently characterize the acoustic field within resonant chambers used to propel acoustic nanomotors. Repeated cycles of observation, inference, and design (OID) are guided by a physical model that describes the rate at which levitating particles approach the nodal plane. Using video microscopy, we observe the relaxation of tracer particles to this plane following the application of the acoustic field. We use sequential Monte Carlo methods to infer model parameters such as the amplitude and frequency of the resonant chamber while accounting for particle-level measurement noise and population-level heterogeneity in the field. Guided by simulated outcomes, we select the optimal design for the next experiment as to maximize the information gain in the relevant parameters. We show how this iterative process serves to discriminate between competing hypotheses and efficiently converges to accurate parameter estimates using only few automated experiments. We discuss the need for model criticism to ensure the validity of the guiding model throughout automated cycles of observation, inference, and design. Furthermore, this work demonstrates how Bayesian methods can learn the parameters of nonlinear, hierarchical models used to describe video microscopy data of active colloids.

42 ENGINEERING↗
DOE OSTI2026

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗
DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Spatiotemporal Pattern Recognition in the PMU Signals in the WECC system

Phasor measurement unit (PMU) data has been used by multiple power system applications, including state estimation, post event analysis, oscillation detection, model validation, and many others. Still, due to its big data nature and availability to general research institutions, comprehensive understanding of the spatiotemporal patterns and underlying mechanisms are incomplete. This study applies a set of signal processing and machine learning approaches aiming at deciphering the characteristic behaviors of multiple phasor measurement units (PMUs) attributes (e.g., voltage, frequency, rate of change of frequency, phase angle), including their auto-correlation, cross-dependence, similarities and discrepancies across units and temporal scales, and distributions of anomalies and their linkages to potential external factors such as weather events. Data analytics are applied to PMUs from the U.S. Western Electricity Coordinating Council (WECC) system. The PMU measurements, recorded events, outages, and weather extremes are all from real world datasets. The findings from the study and mechanistic understanding of the PMU dynamics help provide guidance on system control or preventing blackouts. The derived metrics can be directly used for adjusting or filtering simulated PMU data used for advanced algorithm development.

Hou, Zhangshuan↗
DOE OSTI2020

Iridium-Doped Nanosized Zn–Al Layered Double Hydroxides as Efficient Water Oxidation Catalysts

Layered double hydroxides (LDHs) are an ideal platform to host catalytic metal centers for water oxidation (WO) owing to the high accessibility of water to the interlayer region, which makes all centers potentially reachable and activated. In this paper, we report the syntheses of three iridium-doped zinc–aluminum LDHs (Ir-LDHs) nanomaterials (1–3, with about 80 nm of planar size and a thickness of 8 nm as derived by field emission scanning electron microscopy and powder X-ray diffraction studies, respectively), carried out in the confined aqueous environment of reverse micelles, through a very simple and versatile procedure. These materials exhibit excellent catalytic performances in WO driven by NaIO 4 at neutral pH and 25 °C, with an iridium content as low as 0.5 mol % (~0.8 wt %), leading to quantitative oxygen yields (based on utilized NaIO 4 , turnover number up to ~10,000). Nanomaterials 1–3 display the highest ever reported turnover frequency values (up to 402 min –1 ) for any heterogeneous and heterogenized catalyst, comparable only to those of the most efficient molecular iridium catalysts, tested under similar reaction conditions. The boost in activity can be traced to the increased surface area and pore volume (>5 times and 1 order of magnitude, respectively, higher than those of micrometric materials of size 0.3–1 μm) estimated for the nanosized particles, which guarantee higher noble metal accessibility. X-ray absorption spectroscopy (XAS) studies suggest that 1–3 nanomaterials, as-prepared and after catalysis, contain a mixture of isolated, single octahedral Ir(III) sites, with no evidence of Ir–Ir scattering from second-nearest neighbors, excluding the presence of IrO 2 nanoparticles. The combination of the results obtained from XAS, elemental analysis, and ionic chromatography strongly suggests that iridium is embedded in the brucite-like structure of LDHs, having four hydroxyls and two chlorides as first neighbors. These results demonstrate that nanometric LDHs can be successfully exploited to engineer efficient WOCs, minimizing the amount of iridium used, consistent with the principle of the noble-metal atom economy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phonon modes and Raman signatures of MnBi 2n Te 3n+1 (n=1,2,3,4) magnetic topological heterostructures

An intrinsic antiferromagnetic topological insulator MnBi 2 Te 4 arises when intercalating a Mn-Te bilayer chain in a topological insulator, Bi 2 Te 3 . We present observations on the inter- and intralayer phonon modes of the generalized MnBi 2n Te 3n+1 (n=1,2,3,4) family using cryogenic low-frequency Raman spectroscopy with various polarization configurations. Two peaks at 66 and 112 cm –1 show abnormal perturbation in Raman linewidths below magnetic transition temperature due to spin-phonon coupling. In MnBi 4 Te 7 , B i2 Te 3 layers induce Davydov splitting of the A1g mode around 137 cm –1 at 5 K. The out-of-plane interlayer force constant estimated using the linear chain model was (3.98±0.14)×10 19 N/m 3 , three times weaker than that of Bi 2 Te 3 . Adding more Bi 2 Te 3 layers, such as MnBi 6 Te 10 and MnBi 8 Te 13 , makes Bi 2 Te 3 properties more dominant than magnetic properties. Our work experimentally and theoretically discovers the dynamics of phonon modes of MnBi 2n Te 3n+ 1 family, facilitating utilization of magnetic topological heterostructures.

74 ATOMIC AND MOLECULAR PHYSICS↗