Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damping”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental 3-dimensional tracking of the dynamics of a single electron in the Fermilab Integrable Optics Test Accelerator (IOTA)

We present the results of experimental studies on the transverse and longitudinal dynamics of a single electron in the IOTA storage ring. IOTA is a flexible machine dedicated to beam physics experiments with electrons and protons. A method was developed to reliably inject and circulate a controlled number of electrons in the ring. A key beam diagnostic system is the set of sensitive high-resolution digital cameras for the detection of synchrotron light emitted by the electrons. With 60–130 electrons in the machine, we measured beam lifetime and derived an absolute calibration of the optical system. At exposure times of 0.5 s, the cameras were sensitive to individual electrons. Camera images were used to reconstruct the time evolution of oscillation amplitudes of a single electron in all 3 degrees of freedom. The evolution of amplitudes directly showed the interplay between synchrotron-radiation damping, quantum excitations, and scattering with the residual gas. From the distribution of measured single-electron oscillation amplitudes, we deduced transverse emittances, momentum spread, damping times, and beam energy. Estimates of residual-gas density and composition were calculated from the measured distributions of vertical scattering angles. Combining scattering and lifetime data, we also provide an estimate of the aperture of the ring. To our knowledge, this is the first time that the dynamics of a single electron are tracked in all three dimensions with digital cameras in a storage ring.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Status and future plans for C 3 R&D

C 3 is an opportunity to realize an e + e - collider for the study of the Higgs boson at √s = 250 GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C 3 is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C 3 is the main linac. This paper presents the staged approach towards a facility to demonstrate C 3 technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The primary goal of the C 3 Demonstration R&D Plan is to reduce technical and cost risk by building and operating the key components of C 3 at an adequate scale. This R&D plan starts with the engineering design, and demonstration of one cryomodule and will culminate in the construction of a 3 cryomodule linac with pre-production prototypes. This R&D program would also demonstrate the linac rf fundamentals including achievable gradient and gradient stability over a full electron bunch train and breakdown rates. It will also investigate beam dynamics including energy spread, wakefields, and emittance growth. This work will be critical to confirm the suitability of the C 3 beam parameters for the physics reach and detector performance in preparation for a Conceptual Design Report (CDR), as well as for follow-on technology development and industrialization. The C 3 Demonstration R&D Plan will open up significant new scientific and technical opportunities based on development of high-gradient and high-efficiency accelerator technology. It will push this technology to operate both at the GeV scale and mature the technology to be reliable and provide high-brightness electron beams. The timeline for progressing with C 3 technology development will be governed by practical limitations on both the technical progress and resource availability. It consists of four stages: Stage 0) Ongoing fundamental R&D on structure prototypes, damping and vibrations. Stage 1) Advancing the engineering maturity of the design and developing start-to-end simulations including space-charge and wakefield effects. This stage will include testing of strucutres operating at cryogenic temperatures. Beam tests would be performed with high beam current to test full beam loading. Stage 2) Production and testing of the first cryomodule at cryogenic temperatures. This would provide sufficient experimental data to compile a CDR and it is anticipated for Stage 2 to last 3 years and to culminate with the transport of photo-electrons through the first cryomodule. Stage 3) Updates to the engineering design of the cryomodules, production of the second and third cryomodule and their installation. Lower charge and lower emittance beams will be used to investigate emittance growth. The successful full demonstration of the 3 cryomodules to deliver up to a 3 GeV beam and achieve the C 3 five gradient will allow a comprehensive and robust evaluation of the technical design of C 3 as well as mitigate technical, schedule, and cost risks required to proceed with a Technical Design Report (TDR).

radiation hardened magnets↗
DOE OSTI2021

Effects of resonant magnetic perturbations on radial electric fields in DIII-D tokamak

Here, gyrokinetic simulations of DIII-D tokamak equilibrium find that resonant magnetic perturbation (RMP) drives a neoclassical non-ambipolar electron particle flux, which causes a rapid change of equilibrium radial electric fields consistent with experimental observations during the suppression of the edge localized mode (ELM). The simulation results provide a support for the conjecture that RMP-induced changes of radial electric fields lead to the enhanced turbulent transport at the pedestal top during the ELM suppression (Taimourzadeh et al 2019 Nucl. Fusion 59 046005). Furthermore, gyrokinetic simulations of collisionless damping of zonal flows show that resonant responses to the RMP decrease the residual level of the zonal flows and damp the geodesic acoustic mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Reionization history constraints from neural network based predictions of high-redshift quasar continua

ABSTRACT Observations of the early Universe suggest that reionization was complete by z ∼ 6, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Lyα damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Lyα, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening intergalactic medium (IGM). In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Lyα. Our Quasar Spectra from Artificial Neural Network based predictive Regression Algorithm(QSANNdRA) improves the error in this reconstruction compared to the state-of-the-art principal component analysis (PCA) based model in the literature by 14.2 per cent on average, and provides an improvement of 6.1 per cent on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1 per cent and 16.8 per cent when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at z = 7.0851 and ULAS J1342+0928 at z = 7.5413, respectively. Using our more accurate reconstructions of these two z > 7 quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{H\, \small{I}} = 0.25^{+0.05}_{-0.05}$ at z = 7.0851 and $\bar{x}_\mathrm{H\, \small{I}} = 0.60^{+0.11}_{-0.11}$ at z = 7.5413. Our results are consistent with the literature and favour a rapid end to reionization.

Ďurovčíková, Dominika↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dynamical instability in multi-orbiter systems with gas friction

ABSTRACT Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

Li, Jiaru (ORCID:0000000155507421)↗
DOE OSTI2025

Comparative analysis of plasmon modes in layered Lindhard metals and strange metals

The enigmatic strange metal remains one of the central unsolved problems of 21st century science. Understanding this phase of matter requires knowledge of the momentum- and energy-resolved dynamic charge susceptibility 𝜒⁡(𝑞,𝜔), especially at finite momentum. Inelastic electron scattering (EELS), performed in either transmission or reflection geometry, is a powerful probe of 𝜒⁡(𝑞,𝜔). For the prototypical strange metal Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝑥 , transmission- and reflection EELS, and infrared (IR) spectroscopy agree at 𝑞∼0, all revealing a highly damped plasmon near 1 eV. At larger 𝑞, however, EELS results show unresolved discrepancies. Since IR data are highly reproducible, it is advantageous to use IR data to calculate what the expected EELS response should be at modest 𝑞. Building on prior momentum-resolved reflection geometry M-EELS work [Chen et al., Phys. Rev. B 109, 045108 (2024)], we extend this approach to transmission EELS for finite stacks of metallic layers, comparing a “textbook” Lindhard metal to a strange metal. In the Lindhard case, the low-𝑞 response is dominated by long-lived, standing wave plasmon modes arising from interlayer Coulomb coupling, with in-plane dispersions that resemble the well-known Fetter modes of layered metals. This behavior depends only on the geometry and the long-range nature of the Coulomb interaction and is largely insensitive to layer details. At larger 𝑞, the response reflects the microscopic properties of individual layers. For the strange metal, calculations based on IR data predict a highly damped plasmon with weak dispersion and no distinct surface mode. While our results match IR and M-EELS at low 𝑞, they do not reproduce any published EELS spectra at large 𝑞, highlighting unresolved discrepancies that demand further experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Measuring the growth of Alfvén wave parametric decay instability using counter-propagating waves: Theory and simulations

The parametric decay instability (PDI) of Alfvén waves—where a pump Alfvén wave decays into a backward-propagating child Alfvén wave and a forward ion acoustic wave—is a fundamental nonlinear wave-wave interaction and holds significant implications for space and laboratory plasmas. However, to date there has been no direct experimental measurement of PDI. Here, we propose a novel and experimentally viable scheme to quantify the growth of Alfvén wave PDI on a linear device using a large pump Alfvén wave and a small counter-propagating seed Alfvén wave, with the seed-wave frequency tuned to match the backward Alfvén wave generated by standard PDI. Using hybrid simulations, we show that energy transfer from the pump to the seed reduces the latter's spatial damping. By comparing seed-wave amplitudes with and without the pump wave, this damping reduction can be used as a direct and reliable proxy for PDI growth. The method is validated in our simulations across a range of plasma and wave parameters and agrees well with theoretical predictions. Notably, the scheme exhibits no threshold for PDI excitation and is, in principle, readily implementable under current laboratory conditions. Finally, this scheme is a critical step toward solving the challenge of experimentally accessing Alfvén wave PDI and provides an elegant method that may be used to validate fundamental theories of parametric instabilities in controlled laboratory settings.

heliosphere↗
DOE OSTI2022

Theory of magnon polaritons in quantum Ising materials

Here, we present a theory of magnon-polaritons in quantum Ising materials, and develop a formalism describing the coupling between light and matter as an Ising system is tuned through its quantum critical point. The theory is applied to Ising materials having multilevel single-site Hamiltonians, in which multiple magnon modes are present, such as the insulating Ising magnet LiHoF 4 . We find that the magnon-photon coupling strengths may be tuned by the applied transverse field, with the coupling between the soft mode present in the quantum Ising material and a photonic resonator mode diverging at the quantum critical point of the material. A fixed system of spins will not exhibit the diamagnetic response expected when light is coupled to mobile spins or atoms. Without the diamagnetic response, one expects a divergent magnon-photon coupling strength to lead to a superradiant quantum phase transition. However, this neglects the effects of damping and decoherence present in any real system. We show that damping and decoherence may block the superradiant quantum phase transition, and lead to weak coupling between the soft magnon mode and the resonator mode. The results of the theory are applied to experimental data on the model system LiHoF 4 in a microwave resonator.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Experimental evidence for the two-path description of neutron spin echo

Here, we describe an experiment that strongly supports a two-path interferometric model in which the spin-up and spin-down components of each neutron propagate coherently along spatially separated parallel paths in a typical neutron spin-echo small-angle scattering (SESANS) experiment. Specifically, we show that the usual semi-classical, single-path treatment of Larmor precession of a polarized neutron in an external magnetic field predicts a damping as a function of the spin-echo length of the SESANS signal obtained with a periodic phase grating when the transverse width of the neutron wave packet is finite. However, no such damping is observed experimentally, implying either that the Larmor model is incorrect or that the transverse extent of the wave packet is very large. In contrast, we demonstrate theoretically that a quantum-mechanical interferometric model in which the two mode-entangled (i.e., intraparticle entangled) spin states of a single neutron are separated in space when they interact with the grating accurately predicts the measured SESANS signal, which is independent of the wave packet width.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Enhanced antiferromagnetic resonance linewidth in NiO/Pt and NiO/Pd

In this paper, we investigate the enhancement of the antiferromagnetic damping in the sintered NiO-HM (where HM is heavy metal) granular systems having NiO/HM interfaces, where HM = Pt or Pd. Under the assumption of the spin pumping model, we derive the effective interfacial damping conductance g eff , the parameter which characterizes the upper-bound estimate of the spin pumping effect, to be 12 ± 1 nm ₋2 and 5 ± 1 nm ₋2 for the NiO/Pt and the NiO/Pd interfaces, respectively, at room temperature. We report $g_{\text{eff}}$ experimentally derived in this study are an important milestone in antiferromagnetic spintronics, giving a guideline for various spin current transfer and spin interaction phenomena with antiferromagnets where the spin mixing conductance is involved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Density functionals combined with van der Waals corrections for graphene adsorbed on layered materials

Standard density functionals like Perdew-Burke-Ernzerhof (PBE) or Strongly Constrained and Appropriately Normed (SCAN) need a correction to account for long-range van der Waals (vdW) interaction. The damped Zaremba-Kohn model (dZK) starts from a formula for the vdW interaction of a distant atom with a solid surface, both with known dielectric properties, damps this formula at short range, and then treats an adsorbed molecule or atomic layer as a collection of renormalized atoms. Here, we extend this model to graphene adsorbed on semiconducting layered materials [bulk graphite and hBN, and multilayer transition metal dichalcogenides (TMD)] by including the C 4 asymptotic term and multiple electrostatic image effects due to the two surfaces of the substrate slabs in the vdW calculations. The resulting SCAN-vdW-dZK and PBE-vdW-dZK give approximately the same results for the systems considered here, in agreement with available reference values. The predicted binding energies are roughly 25% lower than those from SCAN+rVV10 (revised Vydrov and Van Voorhis 2010), and ~15% lower than those from PBE+rVV10. Since SCAN+rVV10 usually overbinds, the predicted binding energies by SCAN-vdW-dZK and PBE-vdW-dZK are expected to be closer to the true values. The predicted equilibrium binding distances from SCAN-vdW-dZK and PBE-vdW-dZK are slightly larger (~ 0.1 Å) than those from SCAN+rVV10, and close to those from SCAN. The binding energy depends upon the number of substrate layers more strongly in vdW-dZK than in rVV10. The C 4 -term contributions can be 40% of the total vdW interactions, and the C 5 term contributes about 10%. The effects of images and the back surfaces of slabs can contribute about 4–10%. The vdW interaction energy power laws from the vdW-dZK model for graphene adsorbed on multilayer MoS 2 show slowly varying decay to the pairwise exponent –4 with increasing separation $\textit{D}$, very similar to those obtained from the random-phase approximation and renormalization group approaches by Ambrosetti et al. Both the PBE-vdW-dZK and SCAN-vdW-dZK give a greater increase in interlayer binding energy when the TMD substrate changes from monolayer to four-layer for the graphene/TMD adsorption systems. This is consistent with the relevant electron-energy-loss spectroscopy experimental results showing increased dielectric response from the substrate with increasing substrate layer number.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Independence of the spin current from the Néel vector orientation in antiferromagnet CoO

Spin pumping from ferromagnetic Fe into antiferromagnetic CoO across a Ag spacer layer was studied using ferromagnetic resonance (FMR) in Py/CoO/Ag/Fe/Ag(001). The thin Py film on top of CoO permits an alignment of the CoO Neél vector through field cooling in two otherwise equivalent [110] and [110] crystalline axes which are parallel and perpendicular to the Fe magnetization direction, respectively. Fe FMR linewidth is measured as a function of Ag thickness in 10-20-GHz frequency range and in 180-330 K temperature range. Furthermore, we find that there exists an anisotropy in the Fe FMR damping for parallel and perpendicular alignment of the Fe and CoO spins. However, such anisotropic damping exists only at thin Ag thickness where there exists a magnetic interlayer coupling between Fe and CoO, and vanishes at thick Ag thickness where the interlayer coupling becomes negligible but permitting spin-current transmission into CoO. Our result indicates the absence of anisotropic spin current for parallel and perpendicular alignment of the Fe and CoO spin axes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Sound in a System of Chiral One-Dimensional Fermions

In this paper, we consider a system of one-dimensional fermions moving in one direction, such as electrons at the edge of a quantum Hall system. At sufficiently long timescales the system is brought to equilibrium by weak interactions between the particles, which conserve their total number, energy, and momentum. The time evolution of the system near equilibrium is described by hydrodynamics based on the three conservation laws. We find that the system supports three sound modes. In the low-temperature limit one mode is a pure oscillation of particle density, analogous to the ordinary sound. The other two modes involve oscillations of both particle and entropy densities. In the presence of disorder, the first sound mode is strongly damped at frequencies below the momentum relaxation rate, whereas the other two modes remain weakly damped.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dynamical properties of a driven dissipative dimerized $S$ =1/2 chain

We consider the dynamical properties of a gapped quantum spin system coupled to the electric field of a laser, which drives the resonant excitation of specific phonon modes that modulate the magnetic interactions. Here, we deduce the quantum master equations governing the time-evolution of both the lattice and spin sectors, by developing a Lindblad formalism with bath operators providing an explicit description of their respective phonon-mediated damping terms. We investigate the nonequilibrium steady states (NESS) of the spin system established by a continuous driving, delineating parameter regimes in driving frequency, damping, and spin-phonon coupling for the establishment of physically meaningful NESS and their related nontrivial properties. Focusing on the regime of generic weak spin-phonon coupling, we characterize the NESS by their frequency and wave-vector content, explore their transient and relaxation behavior, and discuss the energy flow, the system temperature, and the critical role of the type of bath adopted. Our study lays a foundation for the quantitative modeling of experiments currently being designed to control coherent many-body spin states in quantum magnetic materials.

1-dimensional spin chains↗
DOE OSTI2021

Unifying femtosecond and picosecond single-pulse magnetic switching in Gd-Fe-Co

Many questions are still open regarding the physical mechanisms behind the magnetic switching in Gd-Fe-Co alloys by single optical pulses. Phenomenological models suggest a femtosecond scale exchange relaxation between sublattice magnetization as the driving mechanism for switching. The recent observation of thermally induced switching in Gd-Fe-Co by using both several picosecond optical laser pulse as well as electric current pulses has questioned this previous understanding. This has raised the question of whether or not the same switching mechanics are acting at the femtosecond and picosecond scales. In this work, we aim at filling this gap in the understanding of the switching mechanisms behind thermal single-pulse switching. To that end, we have studied experimentally thermal single-pulse switching in Gd-Fe-Co alloys, for a wide range of system parameters, such as composition, laser power, and pulse duration. We provide a quantitative description of the switching dynamics using atomistic spin dynamics methods with excellent agreement between the model and our experiments across a wide range of parameters and timescales, ranging from femtoseconds to picoseconds. Furthermore, we find distinct element-specific damping parameters as a key ingredient for switching with long picosecond pulses and argue that switching with pulse durations as long as 15 ps is possible due to a low damping constant of Gd. Our findings can be easily extended to speed up dynamics in other contexts where ferrimagnetic Gd-Fe-Co alloys have been already demonstrated to show fast and energy-efficient processes, e.g., domain-wall motion in a track and spin-orbit torque switching in spintronics devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

van der Waals corrected density functionals for cylindrical surfaces: Ammonia and nitrogen dioxide adsorbed on a single-walled carbon nanotube

In this work, we extend the damped Zaremba-Kohn model (dZK) for long-range dispersion interaction between a molecule and a planar surface to molecules adsorbed on a curved cylindrical surface, and employ this extended model as an additive correction to the semilocal density functionals PBE (Perdew-Burke-Ernzerhof) and SCAN (strongly constrained and appropriately normed). The resulting PBE+vdW (van der Waals)-dZK and SCAN+vdW-dZK are applied to two systems, NH3 and NO2 molecules adsorbed on a single-wall carbon nanotube (CNT), for calculations of binding energies and equilibrium distances. For comparison, the results from vdW nonlocal functionals, such as SCAN+rVV10 and PBE+rVV10, are also presented. The binding energies from PBE+rVV10 (Vydrov and Van Voorhis), SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK are about 70–115 meV for the system of CNT + NH 3 and 300–500 meV for the system of CNT + NO 2 . The results from PBE+vdW-dZK and SCAN+vdW-dZK are closer to each other than those from PBE+rVV10 and SCAN+rVV10 are. The relatively closer results from PBE+vdW-dZK and SCAN+vdW-dZK indicate the consistency of our developed vdW–dZK model for cylindrical surfaces. All methods, including PBE, SCAN, PBE+rVV10, SCAN+rVV10, PBE+vdW-dZK, and SCAN+vdW-dZK, give approximately the same binding energy differences between two adsorption configurations (types I and II) for the two systems. This implies that the two adsorption sites have approximately the same adsorption stability. The exponent of the vdW interaction power law from our vdW-dZK model for the two systems is about 0 at short distance, largely due to the damping factor, and tends slowly to –4 to –4.5 at distances D about 20–50 Å. At even larger distances, the vdW power-law exponent approaches –5. This feature is very similar to the one calculated with random-phase approximation and renormalization group approaches, supporting the applicability of our methods. Our developed vdW-dZK method provides a highly efficient and reliable method for large systems with cylindrical surfaces, such as vdW interactions with nanotubes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Quasi-two-dimensional ferromagnetism and anisotropic interlayer couplings in the magnetic topological insulator MnBi 2 Te 4

MnBi 2 Te 4 (MBT) is a promising van der Waals layered antiferromagnetic (AFM) topological insulator that combines a topologically nontrivial inverted Bi-Te band gap with ferromagnetic (FM) layers of Mn ions. Here, the inelastic neutron scattering on single crystals reported here describes rather complex magnetism in MBT. The magnetic anisotropy that controls the bulk and surface magnetic field response of MBT is found to have contributions from both single-ion and interlayer two-ion terms. A description of the quasi-two-dimensional intralayer FM spin waves requires long-range, competing FM and AFM interactions and anomalous damping. While this might suggest carrier-mediated magnetic coupling, abinitio calculations in insulating MBT also find long-range interactions, and classical spin dynamics simulations suggest that magnetic vacancies are at least partially responsible for observations of anomalous damping near the zone boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗