Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave prediction”

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

Detectability of QCD phase transitions in binary neutron star mergers: Bayesian inference with the next generation gravitational wave detectors

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the nrpmw waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f$$^{peak}_{2}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f$$^{peak}_{2}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f$$^{peak}_{2}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed 1.6⁢σ. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. Furthermore, unless the phase transition is “strong,” disentangling it from other hadronic physics uncertainties will require significant theory improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Universal CMB 𝐵-mode spectrum from early causal tensor sources

Many early Universe scenarios predict postinflationary tensor perturbations from causality-limited, subhorizon sources. While the microphysical details may be different, as long as these sources are bounded in duration and correlation length, their tensor power spectra exhibit a universal scaling behavior at small wave number: 𝒫 ℎ ⁡(𝑘)∝𝑘 3 , corresponding to white noise on superhorizon scales at the time of production. If these early causal tensor sources (ECTs) exclusively produce gravitational waves before redshift 𝑧 ∼10 5 , this scaling is realized on all of the scales observed in the CMB, and thus yields a universal multipole distribution for the 𝐵-mode angular power spectrum. Unlike the scale-invariant distributions of inflationary 𝐵 modes, ECTs generically predict enhanced power on small scales and suppressed power on large scales, which allows these source classes to be distinguished given measurements over a sufficient range of angular scales. In this paper, we introduce a unified framework for characterizing ECTs and demonstrate how their universal infrared scaling manifests in low-frequency observables, including CMB 𝐵 modes and stochastic gravitational wave spectral densities. We illustrate this mapping with representative case studies of this universality class involving first-order phase transitions, topological defects, and enhanced scalar perturbations, which source tensor modes at second order in perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2020

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim (Wave Energy Converter SIMulator) is an open-source code for simulating wave energy converters, which has been actively developed and applied to simulate a wide variety of device archetypes and has become a popular tool since its initial release in 2014. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) within the MATLAB/SIMULINK environment. Figure 1 illustrates a general wave-to-wire model which begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of a single WEC (or array), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modelling the entire wave-to-wire is difficult and encompass multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM) based-frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Table 1 lists additional WEC-Sim functionalities, which are created using prebuilt Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems.

hydrodynamics modeling↗
DOE OSTI2023

Model Predictive Control-Based Trajectory Shaper for Safe and Efficient Adaptive Cruise Control

Recent studies show that commercially-available adaptive cruise control (ACC) systems are string-unstable, indicating that ACC-driven vehicles amplify speed fluctuations from downstream traffic and induce stop-and-go waves. Moreover, it is challenging to revise the original control algorithm of an ACC system to achieve string stability due to its internal complexity and powertrain uncertainties. To achieve desired control performance given a string-unstable ACC system and circumvent revising the original control algorithm, this study proposes a model predictive control-based trajectory shaper (MPC-TS), which only modifies the sensor-measured trajectory information (i.e., position and speed) of the preceding vehicle. The proposed MPC-TS leverages the input shaping technique to generate reference trajectory to improve string stability, while incorporating tracking errors and vehicle acceleration/deceleration magnitude in the MPC cost function and constraining fluctuations of vehicle speed and spacing to ensure desired car-following performance. Numerical experiments validate the control performance of ACC with the proposed MPC-TS in terms of string stability, safety, traffic efficiency, and comfort.

Zhou, Anye↗
DOE OSTI2023

BICEP Array: 150 GHz Detector Module Development

The Background Imaging of Cosmic Extragalactic Polarization (BICEP)/Keck (BK) collaboration is currently leading the quest for the highest-sensitivity measurements of the polarized cosmic microwave background (CMB) anisotropies on a degree scale with a series of cryogenic telescopes, of which BICEP Array (BA) is the latest Stage-3 upgrade with a total of ~ 32,000 detectors. The instrument comprises 4 receivers spanning 30–270 GHz, with the low-frequency 30/40 GHz deployed to the South Pole Station in late 2019. The full complement of receivers is forecast to set the most stringent constraints on the tensor-to-scalar ratio r. Building on these advances, the overarching small-aperture telescope concept is already being used as the reference for further Stage-4 experiment design. This paper describes the development of the BICEP Array 150 GHz detector module and its fabrication requirements, with highlights on the high-density time division multiplexing (TDM) design of the cryogenic circuit boards. The low-impedance wiring required between the detectors and the frst stage of superconducting quantum interference device amplifers is crucial to maintaining a stable bias current on the detectors. Here, a novel multilayer FR4 Printed Circuit Board with superconducting traces, capable of reading out up to 648 detectors, is detailed along with its validation tests. An ultra-high-density TDM detector module concept we developed for a CMB-S4-like experiment that allows up to 1920 detectors to be read out is also presented. TDM has been chosen as the detector readout technology for the Cosmic Microwave Background Stage-4 (CMB-S4) experiment based on its proven low-noise performance, predictable costs, and overall maturity of the architecture. The heritage for TDM is rooted in mm- and sub-mm-wave experiments dating back 20 years and has since evolved to support a multiplexing factor of 64x in Stage-3 experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Enabling Multireference Calculations on Multimetallic Systems with Graphic Processing Units

Modeling multimetallic systems efficiently enables faster prediction of desirable chemical properties and the design of new materials. This work describes an initial implementation for performing multireference wave function method localized active-space self-consistent field (LASSCF) calculations through the use of multiple graphics processing units (GPUs) to accelerate time-to-solution. Density fitting is leveraged to reduce memory requirements, and we demonstrate the ability to fully utilize multi-GPU compute nodes. Performance improvements of 5–10x in total application runtime were observed in LASSCF calculations for multimetallic catalyst systems up to 1200 AOs and an active space of (22e,40o) using up to four NVIDIA A100 GPUs. Furthermore, written with performance portability in mind, a comparable performance is also observed in early runs on the Aurora exascale system using Intel Max Series GPUs.

Algorithms↗
DOE OSTI2020

OC6 Phase I: Investigating the underprediction of low-frequency hydrodynamic loads and responses of a floating wind turbine

Phase I of the OC6 project is focused on examining why offshore wind design tools underpredict the response (loads/motion) of the OC5-DeepCwind semisubmersible at its surge and pitch natural frequencies. Previous investigations showed that the underprediction was primarily related to nonlinear hydrodynamic loading, so two new validation campaigns were performed to separately examine the different hydrodynamic load components. In this paper, we validate a variety of tools against this new test data, focusing on the ability to accurately model the low-frequency loads on a semisubmersible floater when held fixed under wave excitation and when forced to oscillate in the surge direction. However, it is observed that models providing better load predictions in these two scenarios do not necessarily produce a more accurate motion response in a moored configuration.

17 WIND ENERGY↗
DOE OSTI2021

Effective point-charge analysis of crystal fields: Application to rare-earth pyrochlores and tripod kagome magnets 𝑅 ⁢3⁢ Mg⁢ 2 ⁢Sb⁢ 3⁢ O⁢ 14

An indispensable step to understand collective magnetic phenomena in rare-earth compounds is the determination of spatially anisotropic single-ion properties resulting from spin-orbit coupling and crystal field (CF). The CF Hamiltonian has a discrete energy spectrum—accessible to spectroscopic probes such as neutron scattering—controlled by a number of independent parameters reflecting the point symmetry of the magnetic sites. Determining these parameters in low-symmetry systems is often challenging. Here, we describe a general method to analyze CF excitation spectra using adjustable effective point-charges. We benchmark our method to existing neutron-scattering measurements on pyrochlore rare-earth oxides and obtain a universal point-charge model that describes a large family of related materials. We adapt this model to the newly discovered tripod kagome magnets (𝑅 3 ⁢Mg 2 ⁢Sb 3 ⁢O 14 , 𝑅 = Tb, Ho, Er, Yb) for which we report broadband inelastic neutron-scattering spectra. Analysis of these data using adjustable point-charges yields the CF wave functions for each compound. From this, we calculate thermomagnetic properties that accurately reflect our measurements on powder samples and predict the effective gyromagnetic tensor for pseudospin degrees of freedom—a crucial step to understand the exotic collective properties of these kagome magnets at low temperature. We present further applications of our method to other tripod kagome materials and triangular rare-earth compounds 𝑅⁢MgGaO 4 (𝑅 =Yb, Tm). Overall, this study establishes a widely applicable methodology to predict CF and single-ion properties of rare-earth compounds based on interpretable and adjustable models of effective point charges.

Crystal-field theory↗
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim

WEC-Sim is an open-source software for simulating wave energy converters and has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory and Sandia National Laboratories within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with the latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1, which spans a 2-year period from June 2020 to September 2022. New features described here include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim: Preprint

WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
DOE OSTI2020

Charge transfer excitations, pair density waves, and superconductivity in moiré materials

Transition-metal dichalcogenide (TMD) bilayers are a new class of tunable moiré systems attracting interest as quantum simulators of strongly interacting electrons in two dimensions. In particular, recent theory predicts that the correlated insulator observed in WSe 2 /WS 2 at half filling is a charge-transfer insulator similar to cuprates and, upon further hole doping, exhibits a transfer of charge from anionlike to cationlike orbitals at different locations in the moiré unit cell. In this work, we demonstrate that in this doped charge-transfer insulator, tightly bound charge-2$\textit{e}$ excitations can form to lower the total electrostatic repulsion. This composite excitation, which we dub a trimer, consists of a pair of holes bound to a charge-transfer exciton. When the bandwidth of doped holes is small, trimers crystallize into insulating pair density waves at a sequence of commensurate doping levels. Further, when the bandwidth becomes comparable to the pair binding energy, itinerant holes and charge-2$\textit{e}$ trimers interact resonantly, leading to unconventional superconductivity similar to superfluidity in an ultracold Fermi gas near Feshbach resonance. Our theory is broadly applicable to strongly interacting charge-transfer insulators, such as WSe 2 /WS 2 or TMD homobilayers under an applied electric field.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Millimeter-Wave Imaging Technology Development for Real-Time 2D and 3D Fusion Plasma Diagnostics

Including a one year No-Cost Extension, The University of California, Davis has been funded since 2020 to advance innovative millimeter-wave diagnostics for fusion plasma research. This work includes: (1) Developing and applying cutting-edge System-on-Chip millimeter wave Integrated Circuit (IC) technology for Electron Cyclotron Emission Imaging (ECEI) and Microwave Imaging Reflectometry (MIR) on the DIII-D tokamak and the NSTX-U spherical tokamak. (2) Development of radiation-hardened millimeter-wave chip for diagnostics in harsh environments, utilizing wide bandgap Gallium Nitride (GaN) semiconductors. (3) Integrating forward-modeling synthetic diagnostics to enable advanced, real-time plasma prediction and feedback control. Originally planned as a three-year program (2020–2023), the project was granted a one-year no-cost extension, extending its duration through the end of 2024.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Hot-lines topology and the fate of the spin resonance mode in three-dimensional unconventional superconductors

In the quasi-two-dimensional (quasi-2D) copper- and iron-based superconductors, the onset of superconductivity is accompanied by a prominent peak in the magnetic spectrum at momenta close to the wave-vector of the nearby antiferromagnetic state. Such a peak is well described in terms of a spin resonance mode, i.e., a spin-1 exciton theoretically predicted for quasi-2D superconductors with a sign-changing gap. The same theories, however, indicate that such a resonance mode should be absent in a three-dimensional (3D) system with a spherical Fermi surface. This raises the question of the fate of the spin resonance mode in layered unconventional superconductors that are not strongly anisotropic, such as certain heavy-fermion compounds and potentially the newly discovered nickelate superconductor NdNiO 2 . Here, we use the random phase approximation to calculate the dynamical spin susceptibility of 3D superconductors with a d x 2 –y 2 -wave gap symmetry and corrugated cylindrical-like Fermi surfaces. By varying the out-of-plane hopping anisotropy t z /t, we demonstrate that the appearance of a spin resonance mode is determined by the topology of the hot lines, i.e., lines on the Fermi surface that are connected by the magnetic wave vector. For an in-plane antiferromagnetic wave vector, the hot lines undergo a topological transition from open lines to closed loops at a critical t z /t value. The closed hot lines cross the nodal superconducting lines, making the spin resonance mode overdamped and incoherent. In contrast, for an out-of-plane antiferromagnetic wave vector, the hot lines remain open and the spin resonance mode remains sharp. We discuss the experimental implications of our results for the out-of-plane dispersion of the spin resonance mode and, more generally, for inelastic neutron scattering experiments on unconventional superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗
DOE OSTI2024

Observation of Significant Flavor-SU(3) Breaking in the Kaon Wave Function at 12 < Q 2 < 25 GeV 2 and Discovery of the Charmless Decay ψ ( 3770 ) → K S 0 K L 0

We present cross sections for the reaction e + e − → K S 0 K L 0 at center-of-mass energies ranging from 3.51 to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5 fb − 1 . The ratio of neutral-to-charged kaon form factors at large momentum transfers ( 12 < Q 2 < 25 GeV 2 ) is determined to be 0.21 ± 0.01 , which indicates a small but significant effect of flavor-SU(3) breaking in the kaon wave function, and, consequently, excludes the possibility that flavor-SU(3) breaking is the primary reason for the strong experimental violation of the pQCD prediction | F ( π ± ) | / | F ( K ± ) | = f π 2 / f K 2 , where F ( π ± ) and F ( K ± ) are the form factors, and f π and f K are the decay constants of charged pions and kaons, respectively. We also observe a significant signal for the charmless decay ψ ( 3770 ) → K S 0 K L 0 for the first time. Within a 1 σ contour of the likelihood value, the branching fraction for ψ ( 3770 ) → K S 0 K L 0 is determined to be B = ( 2.6 3 − 1.59 + 1.40 ) × 10 − 5 , and the relative phase between the continuum and ψ ( 3770 ) amplitudes is ϕ = ( − 0.3 9 − 0.10 + 0.05 ) π . The branching fraction is in good agreement with the S - and D -wave charmonia mixing scheme proposed in the interpretation of the “ ρ π puzzle” between J / ψ and ψ ( 3686 ) decays. Published by the American Physical Society 2024

Physics↗