Search NASA⌕ Search

SEARCH · Search NASA

Results for “TENSOR ANALYSIS”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗
DOE OSTI2026

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗
DOE OSTI2021

A Type–Curve Approach for Evaluating Aquifer Properties by Interpreting Shallow Strain Measured During Well Tests

Strains occur at shallow depths in response to pressure changes during well tests in an underlying aquifer, and recent developments in instrumentation have made it feasible to measure essentially the full strain tensor. Simulations using poroelastic analyses indicate that shallow normal strains are approximately proportional to the logarithm of time when a well is injecting into or pumping from a deep aquifer or reservoir. The drawdown is also a linear function of log time, as shown by the classic Cooper-Jacob type-curve analysis. The time when the semilog straight line intercepts the zero-strain axis is similar to the time determined from the Cooper-Jacob pressure analysis, and it can be used to estimate hydraulic diffusivity, suggesting that horizontal strain data can be used directly to estimate aquifer properties. This approach was validated using measurements from shallow (30-m deep) borehole strainmeters during an injection test at a 530-m-deep sandstone aquifer/reservoir in Oklahoma. The results show intercept times for the shallow normal strain data are essentially the same as for deep pressure data from an equivalent radial distance. The slopes of the semilog plots of the pressure and the strain increase at the same time, suggesting that they both respond to a lateral aquifer boundary. Significantly, though, strain was measured at shallow depths while the pressure data were measured at 530-m depth. This suggests that strain data from shallow depths could be an effective way to improve the characterization of an underlying aquifer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Constraint on Primordial B-modes from the First Flight of the Spider Balloon-borne Telescope

We present the first linear polarization measurements from the 2015 long-duration balloon flight of Spider, which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in the Spider bands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis of Spider and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and r parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, finding r < 0.11 and r < 0.19, respectively. Roughly half the uncertainty in r derives from noise associated with the template subtraction. New data at 280 GHz from Spider's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

FOS: Physical sciences↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

MIONet: Learning Multiple-Input Operators via Tensor Product

As an emerging paradigm in scientific machine learning, neural operators aim to learn operators, via neural networks, that map between infinite-dimensional function spaces. Several neural operators have been recently developed. However, all the existing neural operators are only designed to learn operators defined on a single Banach space; i.e., the input of the operator is a single function. Here, for the first time, we study the operator regression via neural networks for multiple-input operators defined on the product of Banach spaces. We first prove a universal approximation theorem of continuous multiple-input operators. We also provide a detailed theoretical analysis including the approximation error, which provides guidance for the design of the network architecture. Based on our theory and a low-rank approximation, we propose a novel neural operator, MIONet, to learn multiple-input operators. MIONet consists of several branch nets for encoding the input functions and a trunk net for encoding the domain of the output function. Here, we demonstrate that MIONet can learn solution operators involving systems governed by ordinary and partial differential equations. In our computational examples, we also show that we can endow MIONet with prior knowledge of the underlying system, such as linearity and periodicity, to further improve accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Morphological analysis of the polarized synchrotron emission with WMAP and Planck

The bright polarized synchrotron emission, away from the Galactic plane, originates mostly from filamentary structures. We implement a filament finder algorithm which allows the detection of bright elongated structures in polarized intensity maps. We analyse the sky at 23 and 30 GHz as observed respectively by WMAP and Planck. We identify 19 filaments, 13 of which have been previously observed. For each filament, we study the polarization fraction, finding values typically larger than for the areas outside the filaments, excluding the Galactic plane, and a fraction of about 30% is reached in two filaments. We study the polarization spectral indices of the filaments, and find a spectral index consistent with the values found in previous analysis (about -3.1) for more diffuse regions. Decomposing the polarization signals into the E and B families, we find that most of the filaments are detected in P E , but not in P B . We then focus on understanding the statistical properties of the diffuse regions of the synchrotron emission at 23 GHz. Using Minkowski functionals and tensors, we analyse the non-Gaussianity and statistical isotropy of the polarized intensity maps. For a sky coverage corresponding to 80% of the fainter emission, and on scales smaller than 6 degrees (ℓ > 30), the deviations from Gaussianity and isotropy are significantly higher than 3σ. The level of deviation decreases for smaller scales, however, it remains significantly high for the lowest analised scale (~ 1.5°). When 60% sky coverage is analysed, we find that the deviations never exceed 3σ. Finally, we present a simple data-driven model to generate non-Gaussian and anisotropic simulations of the synchrotron polarized emission. The simulations are fitted in order to match the spectral and statistical properties of the faintest 80% sky coverage of the data maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Search for charged-lepton flavor violation in top quark production and decay in pp collisions at $ \sqrt{s} $ = 13 TeV

Results are presented from a search for charged-lepton flavor violating (CLFV) interactions in top quark production and decay in pp collisions at a center-of-mass energy of 13 TeV. The events are required to contain one oppositely charged electron-muon pair in the final state, along with at least one jet identified as originating from a bottom quark. The data correspond to an integrated luminosity of 138 fb$^{-1}$, collected by the CMS experiment at the LHC. This analysis includes both the production (q → eμt) and decay (t → eμq) modes of the top quark through CLFV interactions, with q referring to a u or c quark. These interactions are parametrized using an effective field theory approach. With no significant excess over the standard model expectation, the results are interpreted in terms of vector-, scalar-, and tensor-like CLFV four-fermion effective interactions. Finally, observed exclusion limits are set at 95% confidence levels on the respective branching fractions of a top quark to an eμ pair and an up (charm) quark of 0.13 × 10$^{-6}$ (1.31 × 10$^{-6}$), 0.07 × 10$^{-6}$ (0.89 × 10$^{-6}$), and 0.25 × 10$^{-6}$ (2.59 × 10$^{-6}$) for vector, scalar, and tensor CLFV interactions, respectively.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments [plus supplemental information]

Here, we present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time-correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective vs. single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

14 SOLAR ENERGY↗
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The use of digital thread for reconstruction of local fiber orientation in a compression molded pin bracket via deep learning

A deep convolutional neural network (DCNN) was used for microstructure reconstruction using artificial intelligence (MR-AI) by predicting local average fiber orientation distributions (FOD) in a 3D prepreg platelet molded composite (PPMC) pin bracket. To train the MR-AI model, surface strain fields from residual stresses simulated in PPMC plates were used as the input to the DCNN. A training dataset included PPMC plates with various degrees of global fiber alignment, based on the information obtained from high-fidelity flow simulation of a pin bracket. Further, the MR-AI model was then deployed to analyze FOD in the 3D pin bracket by conducting thermo-elastic residual stress analysis. Initially, the MR-AI model was established entirely on the synthetic simulation data. Then, a μCT scan of a physically molded pin bracket was used to create a finite element model that provided data for additional validation of the DCNN model. For the μCT scan finite element pin bracket the MR-AI model predicted the distribution of fiber orientation tensor components with MAE of 0.10 indicating a global prediction error of 10%. For the flow simulated pin bracket, the MR-AI model predicted the distribution of fiber orientation tensor components with a global prediction error of 11%. The MR-AI model showed the ability to predict regions of varying alignment in the base and flange of the pin bracket. The proposed MR-AI methodology allows for rapid prediction of FOD in geometrically complex parts and offers a promising path to detecting unique fiber orientation states in molded components.

42 ENGINEERING↗
DOE OSTI2023

Third-order Hall effect in the surface states of a topological insulator

Time reversal and inversion symmetric materials fail to yield linear and nonlinear responses since they possess net zero Berry curvature. However, higher-order Hall response can be generated in these systems upon constraining the crystalline symmetries. Motivated by the recently discovered third-order Hall (TOH) response mediated by Berry connection polarizability, namely, the variation of the Berry connection with respect to an applied electric field, here, we investigate the existence of such a Hall effect in the surface states of hexagonal warped topological insulators (e.g., Bi 2 Te 3 ) under the application of only electric field. Using the semiclassical Boltzmann formalism, we investigate the effect of tilt and hexagonal warping on the Berry connection polarizability tensor and, consequently, the TOH effect, provided the Dirac cone remains gapless. We find that the magnitude of the response increases significantly with increasing tilt strength and warping, and therefore, they can provide the tunability of this effect. In addition, we also explore the effect of chemical doping on the TOH response in this system. Interestingly, we show, based on the symmetry analysis, that the TOH can be the leading-order response in this system, which can directly be verified in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

TuckerMPI: A Parallel C++/MPI Software Package for Large-scale Data Compression via the Tucker Tensor Decomposition

With this study, our goal is compression of massive-scale grid-structured data, such as the multi-terabyte output of a high-fidelity computational simulation. For such data sets, we have developed a new software package called TuckerMPI, a parallel C++/MPI software package for compressing distributed data. The approach is based on treating the data as a tensor, i.e., a multidimensional array, and computing its truncated Tucker decomposition, a higher-order analogue to the truncated singular value decomposition of a matrix. The result is a low-rank approximation of the original tensor-structured data. Compression efficiency is achieved by detecting latent global structure within the data, which we contrast to most compression methods that are focused on local structure. In this work, we describe TuckerMPI, our implementation of the truncated Tucker decomposition, including details of the data distribution and in-memory layouts, the parallel and serial implementations of the key kernels, and analysis of the storage, communication, and computational costs. We test the software on 4.5 and 6.7 terabyte data sets distributed across 100 s of nodes (1,000 s of MPI processes), achieving compression ratios between 100 and 200,000×, which equates to 99--99.999% compression (depending on the desired accuracy) in substantially less time than it would take to even read the same dataset from a parallel file system. Moreover, we show that our method also allows for reconstruction of partial or down-sampled data on a single node, without a parallel computer so long as the reconstructed portion is small enough to fit on a single machine, e.g., in the instance of reconstructing/visualizing a single down-sampled time step or computing summary statistics. The code is available at https://gitlab.com/tensors/TuckerMPI.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The structure of quantum singularities on a Cauchy horizon

Spacetime singularities pose a long-standing puzzle in quantum gravity. Unlike Schwarzschild, a generic family of black holes gives rise to a Cauchy horizon on which, even in the Hartle-Hawking state, quantum observables such as ⟨T μν ⟩ — the expectation value of the stress-energy tensor — can diverge, causing a breakdown of semiclassical gravity. Because they are diagnosed within quantum field theory (QFT) on a smooth background, these singularities may provide a better-controlled version of the spacetime singularity problem, and merit further study. Here, I highlight a mildness puzzle of Cauchy horizon singularities: the ⟨T μν ⟩ singularity is significantly milder than expected from symmetry and dimensional analysis. I address the puzzle in a simple spacetime ${\mathcal{W}}_P$, which arises universally near all black hole Cauchy horizons: the past of a codimension-two spacelike plane in flat spacetime. Specifically, I propose an extremely broad QFT construction in which, roughly speaking, Cauchy horizon singularities originate from operator insertions in the causal complement of the spacetime. The construction reproduces well-known outer horizon singularities (e.g., in the Boulware state), and remarkably, when applied to ${\mathcal{W}}_P$, gives rise to a universal mild singularity structure for robust singularities, ones whose leading singular behavior is state-independent. I make non-trivial predictions for all black hole Cauchy horizon singularities using this, and discuss extending the results beyond robust singularities and the strict near Cauchy horizon limit.

black holes↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

G2Aero (Separable shape tensors for aerodynamic design) [SWR-22-44]

G2Aero is a flexible and practical tool for design and deformation of 2D airfoils and 3D blades using data-driven approaches. G2Aero utilizes the geometry of matrix manifolds—specifically the Grassmannian—to build a novel framework for representing physics-based separable deformations of shapes. G2Aero offers the flexibility to generate perturbations in a customizable way over any portion of the blade. The G2Aero framework utilizes data-driven methods based on a curated database of physically relevant airfoils. Specific tools include: - principal geodesic analysis over normal coordinate neighborhoods of matrix manifolds; - a variety of data-regularized deformations to nominal 2D airfoil shapes; - Riemannian interpolation connecting a sequence of airfoil cross-sections to build 3D blades from 2D data; - consistent perturbations over the span of interpolated 3D blades based on dominant modes from the data-driven analysis.

Doronina, Olga↗