Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor product”

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.

338 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Determining topological order from infinite projected entangled pair states

In this work, we present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators onto states with well-defined anyon flux are used to compute topological S and T matrices (encoding mutual and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient simulation of low-temperature physics in one-dimensional gapless systems

Here, we discuss the computational efficiency of the finite-temperature simulation with minimally entangled typical thermal states (METTS). To argue that METTS can be efficiently represented as matrix product states, we present an analytic upper bound for the average entanglement Rényi entropy of METTS for a Rényi index 0 < q ≤ 1. In particular, for one-dimensional (1D) gapless systems described by conformal field theories, the upper bound scales as O⁡(cN 0 ⁢log⁡β) where c is the central charge and N is the system size. Furthermore, we numerically find that the average Rényi entropy exhibits a universal behavior characterized by the central charge and is roughly given by half of the analytic upper bound. Based on these results, we show that METTS can provide a speedup compared to employing the purification method to analyze thermal equilibrium states at low temperatures in 1D gapless systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

A co-kurtosis PCA based dimensionality reduction with nonlinear reconstruction using neural networks

For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Finally, results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.

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

Efficient simulation of moiré materials using the density matrix renormalization group

Here we present an infinite density-matrix renormalization group (DMRG) study of an interacting continuum model of twisted bilayer graphene (tBLG) near the magic angle. Because of the long-range Coulomb interaction and the large number of orbital degrees of freedom, tBLG is difficult to study with standard DMRG techniques - even constructing and storing the Hamiltonian already poses a major challenge. To overcome these difficulties, we use a recently developed compression procedure to obtain a matrix product operator representation of the interacting tBLG Hamiltonian which we show is both efficient and accurate even when including the spin, valley, and orbital degrees of freedom. To benchmark our approach, we focus mainly on the spinless, single-valley version of the problem where, at half filling, we find that the ground state is a nematic semimetal. Remarkably, we find that the ground state is essentially a k-space Slater determinant, so that Hartree-Fock and DMRG give virtually identical results for this problem. Our results show that the effects of long-range interactions in magic angle graphene can be efficiently simulated with DMRG and open up a new route for numerically studying strong correlation physics in spinful, two-valley tBLG, and other moiré materials in future work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗
DOE OSTI2022

Quantum-inspired method for solving the Vlasov-Poisson equations

Kinetic simulations of collisionless (or weakly collisional) plasmas using the Vlasov equation are often infeasible due to high-resolution requirements and the exponential scaling of computational cost with respect to dimension. Recently, it has been proposed that matrix product state (MPS) methods, a quantum-inspired but classical algorithm, can be used to solve partial differential equations with exponential speed-up, provided that the solution can be compressed and efficiently represented as a MPS within some tolerable error threshold. Here, in this work, we explore the practicality of MPS methods for solving the Vlasov-Poisson equations for systems with one coordinate in space and one coordinate in velocity, and find that important features of linear and nonlinear dynamics, such as damping or growth rates and saturation amplitudes, can be captured while compressing the solution significantly. Furthermore, by comparing the performance of different mappings of the distribution functions onto the MPS, we develop an intuition of the MPS representation and its behavior in the context of solving the Vlasov-Poisson equations, which will be useful for extending these methods to higher-dimensional problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Drell-Yan q T resummation of fiducial power corrections at N 3 LL

We consider Drell-Yan production pp → V*X → LX at small q T << Q, where q T and Q are the total transverse momentum and invariant mass of the leptonic final state L. Experimental measurements require fiducial cuts on L, which in general introduce enhanced, linear power corrections in q T /Q. We show that they can be unambiguously predicted from factorization, and resummed to the same order as the leading-power contribution. For the fiducial q T spectrum, they constitute the complete linear power corrections. We thus obtain predictions for the fiducial q T spectrum to N 3 LL and next-to-leading-power in q T /Q. Matching to full NNLO ($α$$^{2}_{s}$), we find that the linear power corrections are indeed the dominant ones, and once included by factorization, the remaining fixed-order corrections become almost negligible below q T ≲ 40 GeV. We also discuss the implications for more complicated observables, and provide predictions for the fiducial Φ* spectrum at N 3 LL+NNLO. We find excellent agreement with ATLAS and CMS measurements of q T and Φ*. We also consider the $p$$^{ℓ}_{T}$ spectrum. We show that it develops leptonic power corrections in q T /(Q – 2$p$$^{ℓ}_{T}$), which diverge near the Jacobian peak $p$$^{ℓ}_{T}$ ~ Q/2 and must be kept to all powers to obtain a meaningful result there. Doing so, we obtain for the first time an analytically resummed result for the $p$$^{ℓ}_{T}$ spectrum around the Jacobian peak at N 3 LL+NNLO. Our method is based on performing a complete tensor decomposition for hadronic and leptonic tensors. We show that in practice this is equivalent to often-used recoil prescriptions, for which our results now provide rigorous, formal justification. Our tensor decomposition yields nine Lorentz-scalar hadronic structure functions, which for Z/γ* → ℓℓ or W → ℓν directly map onto the commonly used angular coefficients, but also holds for arbitrary leptonic final states. In particular, for suitably defined Born-projected leptons it still yields a LO-like angular decomposition even when including QED final-state radiation. Finally, we also discuss the application to q T subtractions. Including the unambiguously predicted fiducial power corrections significantly improves their performance, and in particular makes them applicable near kinematic edges where they otherwise break down due to large leptonic power corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
NASA NTRS2011

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of l.2xl2 Earth radii in the first phase and l.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial -wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be· used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations; the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗
NASA NTRS2010

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of 1.2 x 12 Earth radii in the first phase and 1.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations, the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗
NASA NTRS1998

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗
DOE OSTI2023

Study of η(1405)/η(1475) in $J/\psi \to \gamma {K}_S^0{K}_S^0{\pi}^0$ decay

Using a sample of (10 . 09 ± 0 . 04) × 10 9 J/ψ decays collected with the BESIII detector, partial wave analyses of the decay $ J/\psi \to \gamma {K}_S^0{K}_S^0{\pi}^0$ are performed within the ${K}_S^0{K}_S^0{\pi}^0$ invariant mass region below 1.6 GeV/ c 2 . The covariant tensor amplitude method is used in both mass independent and mass dependent approaches. Both analysis approaches exhibit dominant pseudoscalar and axial vector components, and show good consistency for the other individual components. Furthermore, the mass dependent analysis reveals that the ${K}_S^0{K}_S^0{\pi}^0$ invariant mass spectrum for the pseudoscalar component can be well described with two isoscalar resonant states using relativistic Breit-Wigner model, i.e., the η (1405) with a mass of $1391.7\pm {0.7}_{-0.3}^{+11.3}$ MeV/ c 2 and a width of $60.8\pm {1.2}_{-12.0}^{+5.5}$ MeV, and the η (1475) with a mass of $1507.6\pm {1.6}_{-32.2}^{+15.5}$ MeV/ c 2 and a width of $115.8\pm {2.4}_{-10.9}^{+14.8}$ MeV. The first and second uncertainties are statistical and systematic, respectively. Alternate models for the pseudoscalar component are also tested, but the description of the ${K}_S^0{K}_S^0{\pi}^0$ invariant mass spectrum deteriorates significantly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗