Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized canonical correlation 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.

Cell-Edge Detection via Selective Cooperation and Generalized Canonical Correlation

Improving the uplink quality of service for users located around the boundaries between cells is a key challenge in cellular systems. Existing approaches relying on power control throttle the rates of cell-center users, while multi-user detection requires accurate channel estimates for the cell-edge users, which is another challenge due to their low received signal-to-noise ratio (SNR). Utilizing the fact that cell-edge user signals are weak but common (received at roughly equal power) at different base stations (BSs), this paper establishes a connection between cell-edge user detection and generalized canonical correlation analysis (GCCA). It puts forth a GCCA-based method that leverages selective BS cooperation to recover the cell-edge user signal subspace even at low SNR. The cell-edge user signals can then be extracted from the resulting mixture via algebraic signal processing techniques. The paper includes theoretical analysis showing why GCCA recovers the correct subspace containing the cell-edge user signals under mild conditions. The proposed method can also identify the number of cell-edge users in the system, i.e., the common subspace dimension. Simulations reveal significant performance improvement relative to various multiuser detection techniques. Cell-edge detection performance is further studied as a function of how many / which BSs are selected, and it is shown that using the closest three BS is always the best choice.

base station cooperation↗

Geodesic-based relaxation for deep canonical correlation analysis

Deep canonical correlation analysis (DCCA) is often applied to paired data samples from diverse sources to extract meaningful common information. However, when the data sources are heterogeneous, some of the useful information may be complementary but not exactly common. In spite of this fact, existing techniques learn maximally correlated representations from multiple views and are formulated so that they aim to yield identical latent subspaces for each view. This approach is sub-optimal in estimating the true signal subspaces for heterogeneous data sources. We propose a residual relaxation for deep canonical correlation analysis (RDCCA) based on a subspace distance metric, which generalizes the existing problem formulation and extracts representations that are better estimates of the actual, non-identical subspaces. We demonstrate that when using such a relaxation, the learned representations are closer to the true ones and that RDCCA outperforms CCA and DCCA in scenarios with heterogeneous data.

Kuschel, Maurice↗

Troposphere-Stratosphere Dynamic Coupling Under Strong and Weak Polar Vortex Conditions

The relationship between Northern Hemisphere (NH) tropospheric and stratospheric wave-like anomalies of spherical zonal wave number (ZWN) 1 is studied by applying Canonical Correlation Analysis (CCA). A lag-correlation technique is used with 10-day lowpass filtered daily time series of 50- and 500-hPa geopotential heights. Generally stratospheric circulation is determined by ultralong tropospheric planetary waves. During winter seasons characterized either by any anomalously strong or weak polar winter vortex different propagation characteristics for waves of ZWN 1 are observed. The non-linear perspective of the results have implications for medium range weather forecast and climate sensitivity experiments.

Perlwitz, Judith↗

Harvest and nitrogen effects on bioenergy feedstock quality of grass-legume mixtures on Conservation Reserve Program grasslands

Perennial grass mixtures established on Conservation Reserve Program (CRP) lands can be an important source of feedstock for bioenergy production. This study aimed to evaluate management practices for optimizing the quality of bioenergy feedstock and stand persistence of grass-legume mixtures under diverse environments. A 5-year field study (2008–2012) was conducted to assess the effects of two harvest timings (at anthesis vs after complete senescence) and three nitrogen (N) rates (0, 56, 112 kg N ha -1 ) on biomass chemical compositions (i.e., cell wall components, ash, volatiles, total carbon, and N contents) and the feedstock energy potential, examined by the theoretical ethanol yield (TEY) and the total TEY (i.e., the product of biomass yield and TEY, L ha -1 ), of cool-season mixtures in Georgia and Missouri and a warm-season mixture in Kansas. The canonical correlation analysis (CCA) was used to investigate the effect of vegetative species transitions on feedstock quality. Although environmental variations (mainly precipitation) greatly influenced the management effect on chemical compositions, the delayed harvest after senescence generally improved feedstock quality. In particular, the overall cell wall concentrations and TEY of the warm-season mixtures increased by approximately 7%. Additional N supplies improved the total TEY of both mixtures by ~1.6–4.2 L ha -1 per 1.0 kg N ha -1 input but likely lowered the feedstock quality, particularly for the cool-season mixture. The cell wall concentrations of cool-season mixture reduced by approximately 3%–6%. The CCA results indicated that the increased legume compositions (under low N input) likely enhanced lignin but reduced ash concentrations. This field research demonstrated that with proper management, grass-legume mixtures on CRP lands can provide high-quality feedstock for bioenergy productions.

09 BIOMASS FUELS↗

Machine Learning for Mapping Multipactor Susceptibility in RF Systems: Capabilities and Generalization Constraints

Multipactor is a surface-driven electron avalanche phenomenon that degrades the performance and reliability of radio-frequency (RF) systems in particle accelerator and vacuum electronics applications. Multipactor behavior in a given device structure is conventionally assessed through susceptibility charts, which provide a parameter-space characterization of the instability. In this work, we assess the capabilities of machine-learning (ML) models to learn and predict such susceptibility charts and analyze the constraints governing their generalization across materials. Using a simulation-derived dataset spanning six distinct secondary-electron-yield material profiles in a canonical two-surface planar geometry, we train supervised regression models and artificial neural networks to predict the time-averaged electron growth rate, δavg, across the relevant parameter space. Model performance is evaluated using metrics that explicitly probe the structure of susceptibility charts, including Intersection over Union, Structural Similarity Index, and correlation analysis. Tree-based ensemble models outperform neural-network models in reconstructing susceptibility regions and in generalizing across material domains. Principal-component analysis reveals disjoint material feature distributions, indicating that the piecewise mode structure of multipactor susceptibility is difficult to represent with a single global model and that generalization is constrained by data coverage rather than by model complexity. An exhaustive reduced-coverage study further shows that sparse material-space coverage can yield mean performance in the same general range but producing large variability in the susceptibility-region overlap. These results clarify the capabilities of ML-based surrogate models for parameter-space characterization of multipactor discharge. They also provide guidance for their appropriate use in RF system design.

43 PARTICLE ACCELERATORS↗

Geochemistry and metagenomics analyses of bacterial community structure in selected waste dumpsites in Lagos Metropolis, Nigeria

Dumpsites are reservoirs of persistent organic pollutants (POPs) and heavy metals (HMs), constituting environmental hazards to humanity. Autochthonous microorganisms in dumpsites exhibit various degrees of responses to contaminants. Unfortunately, there is a dearth of information on the types and concentration of pollutants and the array of microorganisms in these dumpsites which may play important roles in the metabolism of such pollutants or other community processes. Therefore, determining the microbial community structure in such contaminated sites across a municipality is essential for profiling the taxa that would serve as consensus degraders of the pollutants. In this study, soil samples from three dumpsites (Cele, CS; Solous, SS; and Computer Village, CVS) were characterized for geochemical properties using GC-MS, MP-AES, and other analytical protocols, while the dynamics of bacterial communities were evaluated based on their 16S rRNA gene barcodes. A significant difference in the bacterial communities was observed among the dumpsites in relation to the extent of pollution caused by POPs and HMs. CVS, with the highest HM contamination, was rich in Actinobacteria (41.7%) and Acidobacteria (10.2%), in contrast to CS and SS. Proteobacteria (34.1%) and Firmicutes (20%) were the dominant phyla in CS (highest POP contamination), while Bacteroidetes (45.5%) and Proteobacteria (39.9%) were dominant in SS soil.Bacilluswas the dominant genus in the most polluted dumpsite. Canonical correspondence analysis revealed that polycyclic aromatic hydrocarbons (PAHs) and HMs shaped the structure of the bacterial operational taxonomic units (OTUs) in the most polluted dumpsite. Out of a total of 706 OTUs, 628 OTUs exhibited a significant correlation (>50%) with benzo(b)fluoranthene, azobenzene, dibenzofurans, pyrene, dibenzo(a,l)pyrene, Cu, and Zn. In particular, Proteobacteria (Achromobactersp. andSerratiasp.), Bacteroidetes (Zhouiasp.), and Firmicutes (Bacillussp.) were suggested to be pivotal to the ecophysiology of dumpsite soils contaminated with POPs and HMs. The results generally underscored the importance of metagenomic and physicochemical analyses of polluted systems in enabling correlations for useful prediction of drivers of such ecosystems. This will further improve our understanding of the metabolic potential and adaptation of organisms in such systems.

Environmental Sciences & Ecology↗

Validation of a Global Geospace Model With a Systems Science Approach Based on Canonical Correlation Analysis

A systems science approach based on canonical correlation analysis (CCA) is applied as a new, behavioral way to validate global geospace models. The biggest novelty of the technique is that it validates models at a system level, whereby a side‐by‐side comparison is performed of CCA applied to a 30‐day observational and the corresponding simulation data sets comprising quiet, moderate and active times. The simulation used the Multiscale Atmosphere‐Geospace Environment (MAGE) model. It is shown that (a) CCA must be combined with sensitivity analysis to be effective, (b) the MAGE model generally reproduces the observed behavior (more so for quieter time intervals), quantified by the intercorrelations between different variables and (c) the technique identifies the SuperMAG SML index as a quantity for which refinements of the model are needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Path-Based Dictionary Augmentation: A Framework for Improving $k$ -Sparse Image Processing

In this study, we have previously shown that augmenting orthogonal matching pursuit (OMP) with an additional step in the identification stage of each pursuit iteration yields improved $k$ -sparse reconstruction and denoising performance relative to baseline OMP. At each iteration a “path” or geodesic, is generated between the two dictionary atoms that are most correlated with the residual and from this path a new atom that has a greater correlation to the residual than either of the two bracketing atoms is selected. Here, we provide new computational results illustrating improvements in sparse coding and denoising on canonical datasets using both learned and structured dictionaries. The two methods of constructing a path are investigated for each dictionary type: the Euclidean geodesic formed by a linear combination of the two atoms and the 2-Wasserstein geodesic corresponding to the optimal transport map between the atoms. We prove here the existence of a higher-correlation atom in the Euclidean case under assumptions on the two bracketing atoms and introduce algorithmic modifications to improve the likelihood that the bracketing atoms meet those conditions. Although, we demonstrate our augmentation on OMP alone, in general it may be applied to any reconstruction algorithm that relies on the selection and sorting of high-similarity atoms during an analysis or identification phase.

97 MATHEMATICS AND COMPUTING↗

Fundamental microscopic properties as predictors of large-scale quantities of interest: Validation through grain boundary energy trends

Correlations between fundamental microscopic properties computable from first principles, which we term canonical properties, and complex large-scale quantities of interest (QoIs) provide an avenue to predictive materials discovery. Here, we propose that such correlations can be efficiently discovered through simulations utilizing approximate interatomic potentials (IPs), which serve as an ensemble of “synthetic materials”. As a proof of principle we build a regression model relating canonical properties to the symmetric tilt grain boundary (GB) energy curves in face-centered cubic crystals, characterized by the scaling factor in the universal lattice matching model of Runnels et al. (2016), which we take to be our QoI. Our analysis recovers known correlations of GB energy to other properties and discovers new ones. We also demonstrate, using available density functional theory (DFT) GB energy data, that the regression model constructed from IP data is consistent with DFT results, confirming the assumption that the IPs and DFT belong to same statistical pool and thereby validating the approach. Regression models constructed in this fashion can be used to predict large-scale QoIs based on first-principles data and provide a general method for training IPs for QoIs beyond the scope of first-principles calculations.

36 MATERIALS SCIENCE↗

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Disentangling the Hard X-ray Background ROSAT HRI

The goal of our investigation was to improve our understanding of the properties of the discrete X-ray sources that produce the X-ray background (XRB). Many surveys have shown that the XRB at energies of 0.5-3.0 keV is dominated by emission from extragalactic point sources and that a significant fraction of the XRB at higher energies also is produced by discrete sources. In spite of the fact that the bulk of the 0.5-10 keV XRB was demonstrated to arise from extragalactic point sources, the spectral shape of the background presented a difficulty, referred to as the "spectral paradox". Spectra for classes of individual sources generally have been found to be incompatible with the observed energy index, alpha = 0.4 of the XRB over the 2-10 keV energy range. For the 0.3-3.0 keV Einstein band, Macca'caro et al. (1988) derived a mean energy spectral index of a approx. 0.95 for 599 extragalactic sources and for a subset of X-ray selected AGN, found alpha = 1.03(sup +0-05, sub -0.06) . Wilkes and Elvis showed that radio "loudness" was strongly correlated with the source spectrum, such that radio-load quasars exhibited flatter spectra (alpha approx. 0.5), while radio-quiet quasars had steeper spectra (alpha approx. 1). Studies of moderately faint sources in the 0.1-2.0 keV ROSAT band also found rather steep spectra (alpha = 0.96 +/- 0.11 for sources with an average flux of 1.5 x 10(exp -14) ergs/sq cm sec). At higher energies and much higher fluxes, energy spectra of individual AGN suggested a "canonical" alpha = 0.7 energy spectrum. Thus, the best evidence suggested that known classes of AGN could not readily explain the observed X-ray background spectrum. In our ROSAT PSPC analysis, we studied not only the traditional log N - log S, but also the spectral properties of the sources. We computed hardness ratios for individual sources and performed spectral fits to the summed source spectra, averaged in flux bins from 10(exp -15) to 10(exp -12) ergs/sq cm sec. We found that the hardness ratios and source spectra consistently show a hardening in the source spectra as the source fluxes decrease, confirming the suggestion by Hasinger et al. (1993).

Murray, Stephen S.↗

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnon damping and mode softening in quantum double-exchange ferromagnets

We present a comprehensive analysis of the magnetic excitations and electronic properties of fully quantum double-exchange ferromagnets, i.e. systems where ferromagnetic (FM) ordering emerges from the competition between spin, charge, and orbital degrees of freedom, but without the canonical approximation of using classical localized spins. Specifically, we investigate spin excitations within the Kondo lattice-like model, as well as a two-orbital Hubbard Hamiltonian in proximity to the orbital-selective Mott phase. Computational analysis of the magnon dispersion, damping, and spectral weight within these models reveals unexpected phenomena, such as magnon mode softening and the anomalous decoherence of magnetic excitations as observed in earlier experimental efforts, but explained here without the use of the phononic degrees of freedom. We show that these effects are intrinsically linked to incoherent spectral features near the Fermi level, which arise due to the quantum nature of the local (on-site) triplets. This incoherent spectrum leads to a Stoner-like continuum on which spin excitations scatter, governing magnon lifetime and strongly influencing the dynamical spin structure factor. Our study explores the transition from coherent to incoherent magnon spectra by varying the electron density. Furthermore, we demonstrate that the magnitude of the localized spin mitigates decoherence by suppressing the incoherent spectral contributions near the Fermi level. We also discuss the effective J 1 –J 2 spin Hamiltonian, which can accurately describe the large doping region characterized by the magnon-mode softening. Finally, we show that this behavior is also present in multiorbital models with partially filled orbitals, namely, in systems without localized spin moments, provided that the model is in a strong coupling regime. Our results potentially have far-reaching implications for understanding FM ordering in various multi-band systems. These findings establish a previously unknown direct connection between the electronic correlations of those materials and spin excitations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗