Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional Principal Component 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 163 records · Page 9

Principal Component-Based Radiative Transfer Model (PCRTM) for Hyperspectral Sensors: Theoretical Concept - Part I

Modern infrared satellite sensors such as Atmospheric Infrared Sounder (AIRS), Cosmic Ray Isotope Spectrometer (CrIS), Thermal Emission Spectrometer (TES), Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) and Infrared Atmospheric Sounding Interferometer (IASI) are capable of providing high spatial and spectral resolution infrared spectra. To fully exploit the vast amount of spectral information from these instruments, super fast radiative transfer models are needed. This paper presents a novel radiative transfer model based on principal component analysis. Instead of predicting channel radiance or transmittance spectra directly, the Principal Component-based Radiative Transfer Model (PCRTM) predicts the Principal Component (PC) scores of these quantities. This prediction ability leads to significant savings in computational time. The parameterization of the PCRTM model is derived from properties of PC scores and instrument line shape functions. The PCRTM is very accurate and flexible. Due to its high speed and compressed spectral information format, it has great potential for super fast one-dimensional physical retrievals and for Numerical Weather Prediction (NWP) large volume radiance data assimilation applications. The model has been successfully developed for the National Polar-orbiting Operational Environmental Satellite System Airborne Sounder Testbed - Interferometer (NAST-I) and AIRS instruments. The PCRTM model performs monochromatic radiative transfer calculations and is able to include multiple scattering calculations to account for clouds and aerosols.

Liu, Xu↗

Probing Molecular Packing of Amorphous Pharmaceutical Solids Using X-ray Atomic Pair Distribution Function and Solid-State NMR

The structural investigation of amorphous pharmaceuticals is of paramount importance in comprehending their physicochemical stability. However, it has remained a relatively underexplored realm primarily due to the limited availability of high-resolution analytical tools. Here in this study, we utilized the combined power of X-ray pair distribution functions (PDFs) and solid-state nuclear magnetic resonance (ssNMR) techniques to probe the molecular packing of amorphous posaconazole and its amorphous solid dispersion at the molecular level. Leveraging synchrotron X-ray PDF data and employing the empirical potential structure refinement (EPSR) methodology, we unraveled the existence of a rigid conformation and discerned short-range intermolecular C–F contacts within amorphous posaconazole. Encouragingly, our ssNMR 19 F– 13 C distance measurements offered corroborative evidence supporting these findings. Furthermore, employing principal component analysis on the X-ray PDF and ssNMR data sets enabled us to gain invaluable insights into the chemical nature of the intermolecular interactions governing the drug–polymer interplay. These outcomes not only furnish crucial structural insights facilitating the comprehension of the underlying mechanisms governing the physicochemical stability but also underscore the efficacy of synergistically harnessing X-ray PDF and ssNMR techniques, complemented by robust modeling strategies, to achieve a high-resolution exploration of amorphous structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microbeam X-Ray Absorption Near-Edge Spectroscopy of Alloying Elements in the Oxide Layers of Irradiated Zircaloy-2

Hydrogen pickup of zirconium-based fuel cladding and structural materials during in-reactor corrosion can degrade fuel components because the ingress of hydrogen can lead to the formation of brittle hydrides. In the boiling water reactor (BWR) environment, Zircaloy-2 fuel cladding and structural components such as water rods and channels can experience accelerated hydrogen pickup, whereas Zircaloy-4 components exposed to similar conditions do not. Because the principal difference between the two alloys is that Zircaloy-2 contains nickel, accelerated hydrogen pickup has been hypothesized to result from the presence of nickel. However, an understanding of the mechanism by which this acceleration occurs is still lacking. We investigated the link between hydrogen pickup and the oxidation behavior of alloying elements when incorporated into the oxide layers formed on zirconium alloys when corroded in the reactor. Synchrotron radiation microbeam X-ray absorption near-edge spectroscopy (XANES) at the Advanced Photon Source was performed on carefully selected BWR-corroded Zircaloy-2 water rods at an assembly-averaged burnup ranging from 32.8 to 74.6 GWd/MTU to determine the oxidation states of alloying elements, such as iron and nickel, within the oxide layers as a function of distance from the oxide-metal interface at high burnup. Samples were chosen for comparison based on having similar oxide thicknesses, processing, elevation, reactors, and fluences but different hydrogen pickup fractions. Examinations of the oxide layers formed on these samples showed that (1) the oxidation states of these alloying elements changed with distance from the oxide-metal interface, (2) these elements exhibited delayed oxidation relative to the host zirconium, and (3) nickel in Zircaloy-2 remained metallic in the oxide layer at a longer distance from the oxide-metal interface than iron. An analysis of these results showed an apparent correlation between the delayed oxidation of nickel and higher hydrogen pickup of Zircaloy-2 at high burnup.

alloying elements↗

WE-AB-207B-07: Dose Cloud: Generating “Big Data” for Radiation Therapy Treatment Plan Optimization Research

Purpose: To provide a tool to generate large sets of realistic virtual patient geometries and beamlet doses for treatment optimization research. This tool enables countless studies exploring the fundamental interplay between patient geometry, objective functions, weight selections, and achievable dose distributions for various algorithms and modalities. Methods: Generating realistic virtual patient geometries requires a small set of real patient data. We developed a normalized patient shape model (PSM) which captures organ and target contours in a correspondence-preserving manner. Using PSM-processed data, we perform principal component analysis (PCA) to extract major modes of variation from the population. These PCA modes can be shared without exposing patient information. The modes are re-combined with different weights to produce sets of realistic virtual patient contours. Because virtual patients lack imaging information, we developed a shape-based dose calculation (SBD) relying on the assumption that the region inside the body contour is water. SBD utilizes a 2D fluence-convolved scatter kernel, derived from Monte Carlo simulations, and can compute both full dose for a given set of fluence maps, or produce a dose matrix (dose per fluence pixel) for many modalities. Combining the shape model with SBD provides the data needed for treatment plan optimization research. Results: We used PSM to capture organ and target contours for 96 prostate cases, extracted the first 20 PCA modes, and generated 2048 virtual patient shapes by randomly sampling mode scores. Nearly half of the shapes were thrown out for failing anatomical checks, the remaining 1124 were used in computing dose matrices via SBD and a standard 7-beam protocol. As a proof of concept, and to generate data for later study, we performed fluence map optimization emphasizing PTV coverage. Conclusions: We successfully developed and tested a tool for creating customizable sets of virtual patients suitable for large-scale radiation therapy optimization research.

60 APPLIED LIFE SCIENCES↗

An Empirical Cumulus Parameterization Scheme for a Global Spectral Model

Realistic vertical heating and drying profiles in a cumulus scheme is important for obtaining accurate weather forecasts. A new empirical cumulus parameterization scheme based on a procedure to improve the vertical distribution of heating and moistening over the tropics is developed. The empirical cumulus parameterization scheme (ECPS) utilizes profiles of Tropical Rainfall Measuring Mission (TRMM) based heating and moistening derived from the European Centre for Medium- Range Weather Forecasts (ECMWF) analysis. A dimension reduction technique through rotated principal component analysis (RPCA) is performed on the vertical profiles of heating (Q1) and drying (Q2) over the convective regions of the tropics, to obtain the dominant modes of variability. Analysis suggests that most of the variance associated with the observed profiles can be explained by retaining the first three modes. The ECPS then applies a statistical approach in which Q1 and Q2 are expressed as a linear combination of the first three dominant principal components which distinctly explain variance in the troposphere as a function of the prevalent large-scale dynamics. The principal component (PC) score which quantifies the contribution of each PC to the corresponding loading profile is estimated through a multiple screening regression method which yields the PC score as a function of the large-scale variables. The profiles of Q1 and Q2 thus obtained are found to match well with the observed profiles. The impact of the ECPS is investigated in a series of short range (1-3 day) prediction experiments using the Florida State University global spectral model (FSUGSM, T126L14). Comparisons between short range ECPS forecasts and those with the modified Kuo scheme show a very marked improvement in the skill in ECPS forecasts. This improvement in the forecast skill with ECPS emphasizes the importance of incorporating realistic vertical distributions of heating and drying in the model cumulus scheme. This also suggests that in the absence of explicit models for convection, the proposed statistical scheme improves the modeling of the vertical distribution of heating and moistening in areas of deep convection.

Rajendran, K.↗

Passive microwave observations of the Wedell Sea during austral winter and early spring

The results of multispectral passive microwave observations (6.7 to 90-GHz) are presented from the cruises of the FS Polarstern in the Weddell Sea from July to December 1986. This paper includes primarily the analysis of radiometric observations taken at ice station sites. Averaged emissivity spectra for first-year (FY) ice were relatively constant throughout the experiment and were not statistically different from FY ice signatures in the Arctic. Detailed ice characterization was carried out at each site to compare the microwave signatures of the ice with the physical properties. Absorption optical depths of FY ice were found to be sufficiently high that only the structure in the upper portions of the ice contributed significantly to interstation emissivity variations. The emissivities at 90-GHz, e(90), had the greatest variance. Both e(90) at vertical polarization and GR(sub e)(90, 18.7)(defined as (e(sub V)(90)-e(sub V)(18.7))/e(sub V)(90 + e(sub V)(18.7)) depended on the scattering optical depth which is a function of the snow grain diameter and layer thickness. The variance showed a latitude dependence and is probably due to an increase in the strength of snow metamorphism nearer the northern edge of the ice pack. The contribution of variations of near-surface brine volume to the emissivity was not significant over the range of values encountered at the station sites. Emissivity spectra are presented for a range of thin ice types. Unsupervised principal component analysis produced three significant eigenvectors and showed a separation among four different surface types: open water, thin ice, FY ice, and FY ice with a thick snow cover. A comparison with SMMR satellite data showed that average ice concentrations derived from the ship's ice watch log were consistent with the satellite concentrations. The surface based emissivities for FY ice were also compared with emissivities calculated from scanning multichannel microwave radiometer (SMMR) satellite radiances. Best agreement was found at 6.7 and 10-GHz, while at 18 and 37-GHz, SMMR emissivities were slightly lower than surface based results. For the three lower frequencies agreement was found within a confidence limit of 95% and for 37-GHz within about 90%.

Grenfell, T. C.↗

Seasonal Variability of Saturn's Atmosphere

The seasonal variability of Saturn's clouds and weather layer, currently displaying a variety of phenomena (convective storms, planetary waves, giant storms and lightning-induced events, etc.) is not yet fully understood. Variations of Saturn's radiance at 5.2 microns, a spectral region dominated by thermal emission in an atmospheric window containing weak gaseous absorption, contain a strong axisymmetric component as well as large discrete features at low and mid-latitudes that are several degrees colder than the planetary average and uncorrelated with features at shorter wavelengths that are dominated by reflected sunlight (Yanamandra-Fisher et al., 2001. Icarus, Vol. 150). The characterization of several fundamental atmospheric properties and processes, however, remains incomplete, namely: How do seasons affect (a) the global distribution of gaseous constituents and aerosols; and (b) temperatures and the stability against convection and large scale-atmospheric transport? Do 5-micron clouds have counterparts at other altitude levels? What changes occur during the emergence of Great White Storms? Data acquired at the NASA/IRTF and NAOJ/Subaru from 1995 - 2011; since 2004, high-resolution multi-spectral and high-spatial imaging data acquired by the NASA/ESA Cassini mission, represents half a Saturnian year or two seasons. With the addition of detailed multi-spectral data sets acquired by amateur observers, we study these dramatic phenomena to better understand the timeline of the evolution of these events. Seasonal (or temporal) trends in the observables such as albedo of the clouds, thermal fields of the atmosphere as function of altitude, development of clouds, hazes and global abundances of various hydrocarbons in the atmosphere can now be modeled. We will present results of our ongoing investigation for the search and characterization of periodicities over half a Saturnian year, based on a non-biased a priori approach and time series techniques (such as Principal Component Analysis, PCA and Lomb-Scargle periodograms, LSP).

Yanamandra-Fisher, Padma A.↗

Probing the Hydrogen Enhanced Near-Field Emission of ITO without a Vacuum-Gap

In-situ monitoring of the multi-component gas streams in high temperature energy conversion devices offer the promises to higher efficiency via improved understanding of the chemical environments during device operation. While conventional resistive based metal oxide semiconductor gas sensors suffer from strong cross-sensitivity, optical sensing approaches offer intrinsic advantages to achieve gas selectivity based on wavelength specific interactions. This manuscript describes a novel method to achieve multicomponent gas sensing during gas exposure of H2, CO2, CH4and CO in humid high temperature environments. A single sensor element comprised of a perovskite La0.3Sr0.7TiO3(LSTO) oxide thin film layer coated on silica optical fiber was used. The sensing responses consisted of two wavelength-specific near infrared (NIR) mechanisms, namely broadband absorption associated with the metal oxide layer, and wavelength localized thermal emission responses associated with the hydroxyl defects within the silica fiber. Principal component analysis (PCA) was applied to couple the two mechanisms to achieve selective gas identification. Successful discrimination of H2and CO2on a single fiber sensor was achieved, where the results are both stable and reversible. This design demonstrates that by coupling multiple optical mechanisms on a single oxide coated fiber sensor, simple platforms can also achieve multi-component sensing functionality without the added complexity of a sensor array. Thus, it suggests a new approach to construct simple, robust and functional sensor designs capable of gas discrimination and quantification in multi-component gas streams.

01 COAL, LIGNITE, AND PEAT↗

MO-DE-207B-03: Improved Cancer Classification Using Patient-Specific Biological Pathway Information Via Gene Expression Data

Purpose: To develop an efficient, pathway-based classification system using network biology statistics to assist in patient-specific response predictions to radiation and drug therapies across multiple cancer types. Methods: We developed PICS (Pathway Informed Classification System), a novel two-step cancer classification algorithm. In PICS, a matrix m of mRNA expression values for a patient cohort is collapsed into a matrix p of biological pathways. The entries of p, which we term pathway scores, are obtained from either principal component analysis (PCA), normal tissue centroid (NTC), or gene expression deviation (GED). The pathway score matrix is clustered using both k-means and hierarchical clustering, and a clustering is judged by how well it groups patients into distinct survival classes. The most effective pathway scoring/clustering combination, per clustering p-value, thus generates various ‘signatures’ for conventional and functional cancer classification. Results: PICS successfully regularized large dimension gene data, separated normal and cancerous tissues, and clustered a large patient cohort spanning six cancer types. Furthermore, PICS clustered patient cohorts into distinct, statistically-significant survival groups. For a suboptimally-debulked ovarian cancer set, the pathway-classified Kaplan-Meier survival curve (p = .00127) showed significant improvement over that of a prior gene expression-classified study (p = .0179). For a pancreatic cancer set, the pathway-classified Kaplan-Meier survival curve (p = .00141) showed significant improvement over that of a prior gene expression-classified study (p = .04). Pathway-based classification confirmed biomarkers for the pyrimidine, WNT-signaling, glycerophosphoglycerol, beta-alanine, and panthothenic acid pathways for ovarian cancer. Despite its robust nature, PICS requires significantly less run time than current pathway scoring methods. Conclusion: This work validates the PICS method to improve cancer classification using biological pathways. Patients are classified with greater specificity and physiological relevance as compared to current gene-specific approaches. Focus now moves to utilizing PICS for pan-cancer patient-specific treatment response prediction.

60 APPLIED LIFE SCIENCES↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Machine learning prediction of the mechanical properties of refractory multicomponent alloys based on a dataset of phase and first principles simulation

In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.

36 MATERIALS SCIENCE↗

Benefits of adaptive radiation therapy in lung cancer as a function of replanning frequency

Purpose: To quantify the potential benefit associated with daily replanning in lung cancer in terms of normal tissue dose sparing and to characterize the tradeoff between adaptive benefit and replanning frequency. Methods: A set of synthetic images and contours, derived from weekly active breathing control images of 12 patients who underwent radiation therapy treatment for nonsmall cell lung cancer, is generated for each fraction of treatment using principal component analysis in a way that preserves temporal anatomical trends (e.g., tumor regression). Daily synthetic images and contours are used to simulate four different treatment scenarios: (1) a “no-adapt” scenario that simulates delivery of an initial plan throughout treatment, (2) a “midadapt” scenario that implements a single replan for fraction 18, (3) a “weekly adapt” scenario that simulates weekly adaptations, and (4) a “full-adapt” scenario that simulates daily replanning. An initial intensity modulated radiation therapy plan is created for each patient and replanning is carried out in an automated fashion by reoptimizing beam apertures and weights. Dose is calculated on each image and accumulated to the first in the series using deformable mappings utilized in synthetic image creation for comparison between simulated treatments. Results: Target coverage was maintained and cord tolerance was not exceeded for any of the adaptive simulations. Average reductions in mean lung dose (MLD) and volume of lung receiving 20 Gy or more (V20{sub lung}) were 65 ± 49 cGy (p = 0.000 01) and 1.1% ± 1.2% (p = 0.0006), respectively, for all patients. The largest reduction in MLD for a single patient was 162 cGy, which allowed an isotoxic escalation of the target dose of 1668 cGy. Average reductions in cord max dose, mean esophageal dose (MED), dose received by 66% of the heart (D66{sub heart}), and dose received by 33% of the heart (D33{sub heart}), were 158 ± 280, 117 ± 121, 37 ± 77, and 99 ± 120 cGy, respectively. Average incremental reductions in MLD for the midadapt, weekly adapt, and full-adapt treatments were 38, 18, and 8 cGy, respectively. Incremental reductions in MED for the same treatments were 57, 37, and 23 cGy. Reductions in MLD and MED for the full-adapt treatment were correlated with the absolute decrease in the planning target volume (r = 0.34 and r = 0.26). Conclusions: Adaptive radiation therapy for lung cancer yields clinically relevant reductions in normal tissue doses for frequencies of adaptation ranging from a single replan up to daily replanning. Increased frequencies of adaptation result in additional benefit while magnitude of benefit decreases.

60 APPLIED LIFE SCIENCES↗

Tools for Performing SBG hyperspectral Observing System Simulation Experiment

One of NASA’s Decadal Survey mission, Surface Biology and Geology (SBG), will include a hyperspectral remote sensing imager, which has a very high spatial resolution and a wide spectral coverage (from UV to Near IR). Unprecedented large data volumes will be generated by the SBG hyperspectral instrument. Before the launch of the new satellite, an Observing System Simulation Experiment (OSSE) can be used to study different designs of the new satellite system. One of the key components in an OSSE study is a radiative transfer model (RTM) or forward model. In this presentation, we will describe a Principal Component-based Radiative Transfer Model (PCRTM) which is capable of simulating atmospheric (TOA) radiance or reflectance spectra from far IR to visible and UV spectral regions (50 wavenumber to 30000 wavenumber) quickly and accurately. Multiple scattering from multiple layers of clouds/aerosols are included in the model. The PCRTM has a very good accuracy relative to reference line-by-line radiative transfer models (LBLRTM), and it saves 3-4 orders of magnitude computational time relative to LBLRTM or MODTRAN. The PCRTM model has been successfully used to analyze large volumes of data from hyperspectral sensors such as AIRS, CrIS, and IASI. It has also been used to perform OSSE studies for the Climate Absolute Radiance and Refractivity Observatory (CLARREO) mission. Another useful tool for the OSSE is surface Bidirectional Reflectance Distribution Function (BRDF) database. It is very crucial for the SBG OSSE to include realistic BRDF spectra. Currently, most of the surface reflectance spectra such as those in the ECOSIS and ECOSTRESS are measured at specific observation geometries. We have developed a hyperspectral bidirectional reflectance (HSBR) model which combines Ross-Li BRDF model with the existing reflectance spectral libraries using a principal component analysis. This HSBR model can provide realistic BRDF spectra under various observation conditions. It can also be used to generate realistic BRDF spectra using measurements from multi-band imagers or spectrometers such as MODIS or VIIRS.

Xu Liu↗

Identification of phases, symmetries and defects through local crystallography

Here we report that advances in electron and probe microscopies allow 10 pm or higher precision in measurements of atomic positions. This level of fidelity is sufficient to correlate the length (and hence energy) of bonds, as well as bond angles to functional properties of materials. Traditionally, this relied on mapping locally measured parameters to macroscopic variables, for example, average unit cell. This description effectively ignores the information contained in the microscopic degrees of freedom available in a high-resolution image. Here we introduce an approach for local analysis of material structure based on statistical analysis of individual atomic neighbourhoods. Clustering and multivariate algorithms such as principal component analysis explore the connectivity of lattice and bond structure, as well as identify minute structural distortions, thus allowing for chemical description and identification of phases. This analysis lays the framework for building image genomes and structure–property libraries, based on conjoining structural and spectral realms through local atomic behaviour.

36 MATERIALS SCIENCE↗

Variational quantum state eigensolver

Extracting eigenvalues and eigenvectors of exponentially large matrices will be an important application of near-term quantum computers. The variational quantum eigensolver (VQE) treats the case when the matrix is a Hamiltonian. Here, we address the case when the matrix is a density matrix ρ. We introduce the variational quantum state eigensolver (VQSE), which is analogous to VQE in that it variationally learns the largest eigenvalues of ρ as well as a gate sequence V that prepares the corresponding eigenvectors. VQSE exploits the connection between diagonalization and majorization to define a cost function C=Tr(ρ~H) where H is a non-degenerate Hamiltonian. Due to Schur-concavity, C is minimized when ρ~=VρV† is diagonal in the eigenbasis of H. VQSE only requires a single copy of ρ (only n qubits) per iteration of the VQSE algorithm, making it amenable for near-term implementation. We heuristically demonstrate two applications of VQSE: (1) Principal component analysis, and (2) Error mitigation.

97 MATHEMATICS AND COMPUTING↗

Model orthogonalization and Bayesian forecast mixing via principal component analysis

One can improve predictability in the unknown domain by combining forecasts of imperfect complex computational models using a Bayesian statistical machine learning framework. In many cases, however, the models used in the mixing process are similar. In addition to contaminating the model space, the existence of such similar, or even redundant, models during the multimodeling process can result in misinterpretation of results and deterioration of predictive performance. In this paper we describe a method based on the principal component analysis that eliminates model redundancy. We show that by adding model orthogonalization to the proposed Bayesian model combination framework, one can arrive at better prediction accuracy and reach excellent uncertainty quantification performance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tectonic tremor characterized by principal-component analysis in the vicinity of central Chile and Argentina

No conclusive evidence has been presented to date for tectonic tremor (TT) in the vicinity of central Chile, where the Nazca Plate is subducting beneath the South American Plate. Subduction in our experimental location (roughly 35.5° S, 70.5° W) is steep and fairly unobstructed compared to the flattened and more seismogenic behavior to the north. We seek to identify TT in our experimental area, whose geodynamics are comparable to tremor-rich subduction zones such as Cascadia and the Nankai Trough. Our method combines time-series visual inspection, frequency-spectrum analysis, waveform cross-correlation, and 3-component (3C) signal covariance to explore the presence of TT in this region. We have identified TT using stations in central Chile and the Malargüe region, Argentina. The TT exhibits similar features to other TT observations worldwide. Waveform characteristics for the TT in our study, particularly dimension of the 3C signal covariance, vary as a function of apparent source location. The duration of one episode of identified TT was about 10 h, which may indicate that the plate interface where tremor generates is strongly coupled. Here, we conclude that our observations reflect features of the local propagation, rather than the tremor source itself.

58 GEOSCIENCES↗

Cardiovascular group

As a starting point, the group defined a primary goal of maintaining in flight a level of systemic oxygen transport capacity comparable to each individual's preflight upright baseline. The goal of maintaining capacity at preflight levels would seem to be a reasonable objective for several different reasons, including the maintenance of good health in general and the preservation of sufficient cardiovascular reserve capacity to meet operational demands. It is also important not to introduce confounding variables in whatever other physiological studies are being performed. A change in the level of fitness is likely to be a significant confounding variable in the study of many organ systems. The principal component of the in-flight cardiovascular exercise program should be large-muscle activity such as treadmill exercise. It is desirable that at least one session per week be monitored to assure maintenance of proper functional levels and to provide guidance for any adjustments of the exercise prescription. Appropriate measurements include evaluation of the heart-rate/workload or the heart-rate/oxygen-uptake relationship. Respiratory gas analysis is helpful by providing better opportunities to document relative workload levels from analysis of the interrelationships among VO2, VCO2, and ventilation. The committee felt that there is no clear evidence that any particular in-flight exercise regimen is protective against orthostatic hypotension during the early readaptation phase. Some group members suggested that maintenance of the lower body muscle mass and muscle tone may be helpful. There is also evidence that late in-flight interventions to reexpand blood volume to preflight levels are helpful in preventing or minimizing postflight orthostatic hypotension.

Blomqvist, Gunnar↗