Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal modes”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Parametric amplification of spin-motion coupling in three-dimensional trapped-ion crystals

Three-dimensional (3D) crystals offer a route to scaling up trapped-ion systems for quantum sensing and quantum simulation applications; however, engineering coherent spin-motion couplings and effective spin-spin interactions in large crystals poses technical challenges associated with decoherence and prolonged timescales to generate appreciable entanglement. Here, we explore the possibility of speeding up these interactions in 3D crystals via parametric amplification. For this purpose, we derive a general Hamiltonian for the parametric amplification of spin-motion coupling that is broadly applicable to normal modes with motion transverse to or along the spatial extent of the crystal. Unlike in lower-dimensional crystals, we find that the ability to faithfully (uniformly) amplify the spin-spin interactions in 3D crystals depends on the physical implementation of the spin-motion coupling. We consider the light-shift gate, and the so-called phase-insensitive and phase-sensitive Mølmer-Sørensen (MS) gates, and we find that only the phase-sensitive MS gate can be faithfully amplified in general 3D crystals. We discuss a situation where nonuniform amplification can be advantageous. We also reconsider the effect of counter-rotating terms on parametric amplification and find that they are not as detrimental as previous studies suggest.

Optical parametric oscillators & amplifiers↗
DOE OSTI2020

Thermal conductance enhanced via inelastic phonon transport by atomic vacancies at Cu/Si interfaces

Understanding and controlling heat transfer across interfaces has become an important issue for the performance of micro- and nanoscale electronics, as well as achieving a high figure of merit for thermoelectrics. Intrinsic and extrinsic defects can have a significant impact on thermal transport in bulk materials and across interfaces, but the mechanism is not well understood. Here, nonequilibrium molecular dynamics simulations are used to determine the impact of interfacial atomic vacancies on thermal transport across a Cu/Si junction. In contrast to the reduction in thermal transport typically seen with bulk defects, we find that by introducing atomic vacancies at a concentration of 6.3% near the interface in either or both materials, the interfacial thermal conductance can be increased by up to 76%. By controlling the initial positions of the vacancies and keeping track of their movements and population, we find that interfacial thermal transport is dependent on temperature, vacancy concentration, and distribution, and a positive correlation between the conductance and point defect activities (extent of vacancy migration, rate of Frenkel defect creation and annihilation) is observed. Further calculations based on the phonon density of states and normal mode decomposition reveal that the increase in interfacial thermal conductance originates primarily from high-frequency phonons, supported by enhanced inelastic phonon transport which contributes to more than 60% of the increase. Our findings suggest a practical way to manipulate inelastic phonon conversion through the presence of defects, which provides an alternative perspective on improving thermal transport between materials with a large lattice mismatch.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optically and thermally driven huge lattice orbital and spin angular momenta from spinning fullerenes

Lattice vibration in solids may carry angular momentum. But unlike the intrinsic spin of electrons, the lattice vibration is rarely rotational. To induce angular momentum, one needs to find a material that can accommodate a twisted normal mode, two orthogonal modes, or excitation of magnons. If excitation is too strong, one may exceed the Lindemann limit, so the material melts. Therefore these methods are not ideal. Here, we theoretically propose a route to phonon angular momentum in a molecular crystal C 60 . We find that a single laser pulse is able to inject a significant amount of angular momentum to C 60 , and the momentum transfer is helicity dependent. Changing from right-circularly polarized light to left-circularly polarized light switches the direction of phonon angular momentum. On the ultrafast timescale, the orbital angular momentum change closely resembles the displacive excitation of coherent phonons, with a cosine-function dependence on time, different from the spin counterpart. Atomic displacements, even under strong laser excitation, remain far below the Lindemann criterion. Under thermal excitation, spinning C 60 even at room temperature generates a huge angular momentum close to several hundred h¯. Our finding opens the door to a large group of fullerenes, from C 60 , C 70 to their endohedral derivatives, where angular momentum can be generated through light or temperature. Furthermore, this paves the way to the phononic control electronic spin and harvesting thermal energy through phonon angular momentum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic and vibrational properties of bulk Cr 2 Ge 2 Te 6 from first-principles calculations

Here we use density functional theory based ab initio calculations to investigate the structural, vibrational, magnetic, and electronic properties of the layered ferromagnet Cr 2 Ge 2 Te 6 (CGT) that has attracted attention for potential spintronic applications. We optimize the structure for a fixed c/a ratio. Our results are in excellent agreement with experimental data on structure, phonons, and electronic properties. The use of van der Waals interactions and relativistic spin-orbit coupling yields accurate lattice constants and interlayer distances. Simultaneously, the adopted theoretical methods lead to a rigorous description of the vibrational normal modes, as well as the valence bands, that are in excellent agreement with Raman spectroscopy and angle-resolved photoemission spectroscopy spectra, respectively. The magnetic moment is slightly overestimated, and the magnetic anisotropy has the correct sign but is greater in magnitude than that in experiment. We find that inclusion of the on-site Coulomb repulsion on Cr−3d orbitals worsens agreement with experiment, in contrast to previous studies. Our results provide a path toward ab initio analysis of magnetic heterostructures based on CGT.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics and transport in the boundary-driven dissipative Klein-Gordon chain

Motivated by experiments on chains of superconducting qubits, we consider the dynamics of a classical Klein-Gordon chain coupled to coherent driving and subject to dissipation solely at its boundaries. As the strength of the boundary driving is increased, this minimal classical model recovers the main features of the "dissipative phase transition" seen experimentally. Between the transmitting and non-transmitting regimes on either side of this transition (which support ballistic and diffusive energy transport respectively), we observe additional dynamical regimes of interest. These include a regime of superdiffusive energy transport at weaker driving strengths, together with a "dissipative breather crystal" regime at stronger driving strengths, which is characterized by emergent translation symmetry, ballistic energy transport, and coherent oscillations of a nonlinear normal mode. Here, we propose a non-local Lyapunov exponent as an experimentally measurable diagnostic of many-body chaos in this system, and more generally in open systems that are only coupled to an environment at their boundaries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Coupled fission fragment angular momenta

Here, nuclear fission produces fragments whose spins are coupled to the relative angular motion via angular momentum conservation. It is shown how ensembles of such spins can readily be obtained by either direct microcanonical sampling or by sampling of the associated normal modes of rotation. The resulting distribution of the spin-spin opening angle is illustrated in various three- and two-dimensional scenarios and it is demonstrated how recent mutually conflicting model calculations can be well reproduced with different assumptions about the scission geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Search for the semileptonic decays Ξ c 0 → Ξ 0 ℓ + ℓ − at Belle

Using the full data sample of 980 fb − 1 collected with the Belle detector at the KEKB asymmetric energy electron-positron collider, we report the results of the first search for the rare semileptonic decays Ξ c 0 → Ξ 0 ℓ + ℓ − ( ℓ = e or μ ). No significant signals are observed in the Ξ 0 ℓ + ℓ − invariant-mass distributions. Taking the decay Ξ c 0 → Ξ − π + as the normalization mode, we report 90% credibility upper limits on the branching fraction ratios B ( Ξ c 0 → Ξ 0 e + e − ) / B ( Ξ c 0 → Ξ − π + ) < 6.7 × 10 − 3 and B ( Ξ c 0 → Ξ 0 μ + μ − ) / B ( Ξ c 0 → Ξ − π + ) < 4.3 × 10 − 3 based on the phase-space assumption for signal decays. The 90% credibility upper limits on the absolute branching fractions of B ( Ξ c 0 → Ξ 0 e + e − ) and B ( Ξ c 0 → Ξ 0 μ + μ − ) are found to be 9.9 × 10 − 5 and 6.5 × 10 − 5 , respectively. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2023

Locating the Source of Forced Oscillations in Transmission Power Grids

A forced oscillation event in power grids refers to a state where malfunctioning or abnormally operating equipment causes persisting periodic disturbances in the system. While power grids are designed to damp most perturbations during standard operations, some of them can excite normal modes of the system and cause significant energy transfers across the system, creating large oscillations thousands of miles away from the source. Localization of the source of such disturbances remains an outstanding challenge due to a limited knowledge of the system parameters outside of the zone of responsibility of system operators. Here, we propose a new method for locating the source of forced oscillations that addresses this challenge by performing a simultaneous dynamic model identification using a principled maximum likelihood approach. We illustrate the validity of the algorithm on a variety of examples where forcing leads to resonance conditions in the system dynamics. Our results establish that an accurate knowledge of system parameters is not required for a successful inference of the source and frequency of a forced oscillation. We anticipate that our method will find a broader application in general dynamical systems that can be well described by their linearized dynamics over short periods of time.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Nitrogen in black phosphorus structure

Group V elements in crystal structure isostructural to black phosphorus with unique puckered two-dimensional layers exhibit exciting physical and chemical phenomena. However, as the first element of group V, nitrogen has never been found in the black phosphorus structure. Here, we report the synthesis of the black phosphorus–structured nitrogen at 146 GPa and 2200 K. Metastable black phosphorus–structured nitrogen was retained after quenching it to room temperature under compression and characterized in situ during decompression to 48 GPa, using synchrotron x-ray diffraction and Raman spectroscopy. We show that the original molecular nitrogen is transformed into extended single-bonded structure through gauche and trans conformations. Raman spectroscopy shows that black phosphorus–structured nitrogen is strongly anisotropic and exhibits high Raman intensities in two A g normal modes. Synthesis of black phosphorus–structured nitrogen provides a firm base for exploring new type of high-energy-density nitrogen and a new direction of two-dimensional nitrogen.

58 GEOSCIENCES↗
DOE OSTI2021

Vibronic and Environmental Effects in Simulations of Optical Spectroscopy

Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.

Chemistry↗
DOE OSTI2021

Role of Pr-Vacancies and O-Interstitials on the Activity and Stability of (Pr 1−x Ln x ) 2 NiO 4 (Ln = La, Nd, Pm, Sm, Gd, Tb, Dy, and Ho) towards Oxygen Reduction Reactions: A DFT Study

Praseodymium nickelate, Pr 2 NiO 4 (PNO), is a promising electrode to promote oxygen reduction reaction (ORR) in a solid oxide fuel cell, but it exhibits phase transformation during electrochemical operation. The origin of the simultaneous phase transformation and high electrochemical performance remains obscure. We carried out a systematic density functional theory study to elucidate the mechanism for this conjugated phenomenon. Charge, electronic structure, and normal-mode analysis suggest the presence of peroxide. Our study shows that the formation of peroxide (O 2 2– ) is attributed to both oxygen interstitials and Pr vacancies. The peroxide species limits the oxygen ion migration due to the additional energy required to break its O–O bond, which leads to a decrease in ORR activity. Subsequently, we investigate the diffusion paths of Pr-ions while comparing them with those of other Ln 3+ ions (La, Nd, Pm, Sm, Gd, Tb, Dy, and Ho) in PNO. The formation energies for various Ln 3+ cation occupancies are calculated, as well as segregation energies in CeO 2 (111) surfaces. Lastly, criteria for effective Ln 3+ dopants are developed. La, Nd, and Pm are proposed as potential substituents in PNO to obtain a stable structure.

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

sta-imex (scripts to investigate stability of IMEX methods)

Set of Matlab scripts to investigate properties of IMEX time stepping methods. Scripts include discretization of a system of linearized equations. One script plots numerical dispersion against analytical dispersion of the system. Another set of scripts uses the system to break it into stiff and non stiff part to investigate stability of Implicit-Explicit Runge-Kutta methods. Scripts are based on Thuburn, J. and Woollings, T. J., Vertical Discretizations for Compressible Euler Equation Atmospheric Models Giving Optimal Representation of Normal Modes , 2005 and Lock, S.-J. and Wood, N. and Weller, H., Numerical analyses of Runge-Kutta implicit-explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models, 2014 expanding ideas from these two papers. Scripts are used to develop IMEX schemes. SAND2020-12594 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Guba, Oksana↗
DOE OSTI2023

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗
DOE OSTI2021

Application of Shift Ex-Core Calculations to a Detailed MSR Model

This report documents the collaborative research conducted between Argonne National Laboratory (ANL) and Oak Ridge National Laboratory (ORNL) on applying the Monte Carlo neutronics code Shift to calculate the radiation dose rates for a detailed molten salt reactor (MSR) model based on the Molten Salt Breeder Reactor (MSBR). In this work scope, the ANL and ORNL teams worked together and applied Shift to calculate the radiological environments of this detailed MSR facility modeled with explicit geometries. The radiological conditions within the MSR facility were modeled for when the reactor is at two different operation modes: normal full power and drained state. The FW-CADIS hybrid method in Shift was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Dose rate maps showed that in the current MSBR numerical model, potential pathways exist for neutrons and gammas to stream through the 8-foot concrete shield to reach the top of the MSBR reactor cell. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the Shift simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates, which suggests that a steel liner in the tank model must be added for calculating the dose rates in the drain cell.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Unveiling the Significance of Correlations in K-Space and Configuration Space for Drift Wave Turbulence in Tokamaks

Turbulence and transport phenomena play a crucial role in the confinement and stability of tokamak plasmas. Turbulent fluctuations in certain physical quantities, such as density or temperature fluctuations, can have a wide range of spatial scales, and understanding their correlation length is important for predicting and controlling the behavior of the plasma. The correlation length in the radial direction is identified as the critical length in real space. The dynamics in real space are of significant interest because transport in configuration space is primarily focused on them. When investigating transport caused by the E×B drift, the correlation length in real space represents the size of E×B whirls. It was numerically discovered that in drift wave turbulence, this length is inversely proportional to the normalized mode number of the fastest growing mode relative to the drift frequency. Considerable time was required before a proper analytical derivation of this condition was accomplished. Therefore, a connection has been established between phenomena occurring in real space and those occurring in k-space. Although accompanied by a turbulent spectrum in k-space with a substantial width, transport in real space is uniquely determined by the correlation length, allowing for accurate transport calculations through the dynamics of a single mode. Naturally, the dynamics are subject to nonlinear effects, with resonance broadening in frequency being the most significant nonlinear effect. Thus, mode number space is once again involved. Resonance broadening leads to the detuning of waves from particles, permitting a fluid treatment. It should be emphasized that the consideration here involves the total electric field, including the induction part, which becomes particularly important at higher beta plasmas.

Weiland, Jan (ORCID:0000000253496025)↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗