Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 415 records · Page 23

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Design Basis Document / Owner’s Technical Specification for Nitrate Salt Systems in CSP Projects: Volume 1 Specifications for Parabolic Trough Projects, Volume 2 - Specifications for Central Receiver Projects, Volume 3 - Narrative

Design Basis Documents / Owner’s Technical Specifications are developed for parabolic trough and central receiver power plants using nitrate salt as the heat transport fluid and the thermal storage medium. The goals were to 1) distill the successful experience with nitrate salt systems from as many commercial projects as possible, 2) provide technical bases for equipment design/selection that an owner can impose on an EPC contractor, 3) compile this information in one location to provide guidance on salt systems that is as broadly applicable to as many projects as possible, and 4) work toward an industry consensus on a Design Basis Document that reflects the lessons learned from earlier commercial projects. The document allows an owner to provide design guidance, and to impose a minimum set of requirements, above and beyond those in the normal Codes and Standards, on an EPC contractor in an effort to avoid a repeat of past mistakes. If successful, this would allow salt systems to achieve the same levels of reliability and availability as commercial parabolic trough plants using diphenyl oxide / biphenyl (i.e., Dowtherm A, Therminol VP-1) as the heat transfer fluid. The report consists of 3 volumes: • Volume 1 - Specifications for Parabolic Trough Projects • Volume 2 - Specifications for Central Receiver Projects • Volume 3 - Narrative Volumes 1 and 2 include discussions of the following topics: • Plant functional requirements for the salt systems • Plant operating states, and transitions between states, for the salt systems • Risk analysis of the principal salt components • Plant requirements for the salt systems to meet the functional and the operating requirements • Type of specification for the principal salt components: functional; or prescriptive • Current state of the art for salt systems. Volume 3 discusses a number of cases in which the salt equipment at commercial parabolic trough and central receiver projects has not met the projected levels of reliability and availability. Possible reasons for the sources of the problems are discussed, as is a range of possible alternate designs that could avoid the known problems.

14 SOLAR ENERGY↗

Quantifying the Impact of Metal Population Distribution in MFI-Supported Mo Catalysts for Methane Dehydroaromatization

Precise evaluation of intrinsic kinetic behavior in Mo/MFI catalysts for methane dehydroaromatization (MDA) is confounded by variations in Mo dispersion and speciation, which are influenced by metal loading and zeolite acidity. This work advances the utility of H 2 -temperature programmed reduction (H 2 -TPR) for characterizing Mo/MFI catalysts by enabling quantitative comparison of MoO x populations distinguished by reduction behavior and linked to initial catalytic performance. A comprehensive H 2 reduction pathway is established through systematic H 2 - TPR studies varying catalyst composition (1–10 wt % Mo, Si/Al = 15, 40, ∞), supplemented by UV-Raman spectroscopy, X-ray powder diffraction, N 2 physisorption, NH 3 -TPD, and advanced spectroscopic analysis (in situ XAS with principal component analysis/multivariate curve resolution—alternating least squares). Two low-temperature H 2 -TPR regions capture distinct Mo populations undergoing initial Mo(VI)→Mo(IV) reduction: a lower-temperature population (Mo-RI) associated primarily with highly dispersed, anchored MoO x species expected to predominantly reside within MFI channels, and a higher-temperature population (Mo-RII) corresponding to a broader set of MoO x species that becomes increasingly bulk-like/extrazeolitic at higher Mo loading. Quantification of these populations provides practical, kinetically relevant descriptors for comparing initial MDA rates across catalysts with varying Mo loading, Si/Al ratio, and MoO x heterogeneity. Application of lower-temperature Mo-RI estimates to kinetic measurements reveals a minimum threshold of ∼0.12 × 10 –3 mol Mo-RI/g cat , above which initial forward benzene rates normalized to this population converge despite differences in metal loading and zeolite Brønsted acidity. This threshold coincides with a transition toward a common C 2 -mediated benzene-forming regime, as indicated by approach-to-equilibrium analysis of methane-to-ethane, ethane-to-ethylene, and ethylene-to-benzene reaction steps. Above this threshold, normalized initial benzene rates are nearly invariant with increasing Mo-RII/Mo-RI population ratio, indicating that excess Mo-RII populations, including bulk-like/extrazeolitic MoO x domains present at higher loading, do not measurably suppress benzene formation associated with anchored and mostly channel-confined Mo population under the initial-rate conditions examined.

Mo/MFI↗

Carbon fiber classification using raman spectroscopy

Carbon fiber characterization processes are described that include multi-condition Raman spectroscopy-based examination combined with multivariate data analyses. Methods are a nondestructive material characterization approach that can provide predictions as to carbon fiber bulk physical properties, as well as identification of unknown carbon fiber materials for quality control purposes. The framework of the multivariate analysis methods includes a principal component-based identification protocol including comparison of Raman spectral data from an unknown carbon fiber with a data library of multiple principal component spaces.

Houk, Amanda L.↗

Data Analytics Applied to Coal Fired Boilers for Detecting Leaks

Data analytics were used to detect boiler leaks from five different coal-fired boilers including both subcritical and supercritical systems. Discriminant functions were developed that detected leaks up to two weeks prior to forced plant shutdowns for repairs. The leaks were identified to occur at different sections of the boiler for each plant, including waterwalls, economizer and superheater using conventional process measurement data. Leaking conditions were detected with a high degree of confidence (≪ 1% misclassified observations) and were able to distinguish normal operations from those time periods with steam leaks even while operating the power plants in power cycling mode.Multivariable statistical analyses, including Principal Component (PCA), cluster, and Fischer Discriminant Analysis (FDA) were used to characterize the leak occurrence. Normal and operational states with steam leaks were provided in the original process datasets. These datasets were split into two different groups for training and validation purposes. The data were sorted chronologically, and every third observation was assigned to training the Discriminant Function Model (DFM) while the rest were reserved for validation. PCA was used to reduce dimensionality of the original datasets. Canonical and FDA analyses were used to investigate the relationship between process variables. The outcome of the analyses revealed that nearly 35,000 observations were classified correctly; less than 0.05% of total observations were misclassified to be leaking, i.e. both false positives and false negatives.

Indrawan, Natarianto↗

Quantifying market volume sensitivity to material property modifications in polyhydroxybutyrate: A parametric analysis approach

Polyhydroxybutyrate (PHB), a biodegradable biopolymer, represents a promising alternative to petroleum-based thermoplastics. However, despite consistent market growth, PHB faces persistent commercialization challenges that limit widespread adoption. Existing research has focused predominantly on optimizing PHB production processes, leaving a critical gap in understanding which material property modifications would most effectively enhance market competitiveness. This study addresses this gap by systematically analyzing the relationship between polymer material properties and market performance using U.S. market data from 2008 to 2021 for 21 thermoplastic polymers across 19 material properties. We employed principal component regression to identify property modifications that could maximize market volume while reducing CO 2 emissions. Our parametric analysis revealed that two specific material properties – Hardness Shore A and Sheet Extrusion Temperature – significantly influence PHB marketability across different price points. Market simulations demonstrated that a 10% increase in Hardness Shore A could increase PHB market volume by 431.5 million kg while reducing emissions by 188.7 kg CO 2 . A similar 10% increase to Sheet Extrusion Temperature could yield a 297.5 million kg volume increase and a 99.2 kg CO 2 reduction in emissions. Critically, this approach is agnostic to the specific methods required to achieve these property changes, instead providing material scientists with quantitative, data-driven targets for R&D prioritization. Here, this framework offers a novel methodology for evaluating biopolymer competitiveness and supporting strategic decisions to accelerate PHB market adoption and contribute to decarbonization of the plastics industry.

09 BIOMASS FUELS↗

Practical Guide to Chemometric Analysis of Optical Spectroscopic Data

The methodology and mathematical treatment of several classic multivariate methods for the analysis of spectroscopic data is demonstrated in a straightforward way that can be used as a basis for teaching an undergraduate introductory course on chemometric analysis. The multivariate techniques of classical least squares (CLS), principal component regression (PCR), and partial least squares (PLS), as well as the univariate Beer’s law method have been described and compared, building students’ understanding by starting with the univariate method and progressing step by step into the multivariate methods. Equations for the production of regression vectors from training set spectral data is described and their use demonstrated for the prediction of constituent concentrations on a separate validation set of spectra. Extreme care is taken to ensure consistency in variable formatting of data matrices. This provides a key foundation to understanding how spectral data are manipulated using these different mathematical approaches for building quantitative regression models. Each method is applied to a real-world data set, and the results are discussed to show students the types of information that can be gleaned from each method. A training set comprised of 20 infrared absorbance spectra containing 3 constituents (benzene, polystyrene, and gasoline) of known composition are used to demonstrate the matrix operations for each regression method. A separate set of 12 real-world napalm samples (containing benzene, polystyrene and gasoline) are used as a validation set to demonstrate the ability to utilize the regression models on an unknown dataset. A toolbox (PNNL Chemometric Toolbox) written in MATLAB language is supplied in the Supplemental Information file and can be used as a companion for understanding the development and deployment of the chemometric algorithms described in this paper. The datasets of the infrared spectra are also supplied, allowing users to build and inspect the chemometric models on their own. Finally, the Toolbox includes scripts to assist users in loading their own datasets into MATLAB and performing CLS, PCR, and PLS on their data.

Upper-Division Undergraduate, Analytical Chemistry↗

Probing the Electronic Structure of a Thorium Nitride Complex by Solid-State 15N NMR Spectroscopy

The solid-state 15N NMR powder spectra of the thorium nitride complex, [K(18-crown-6)(THF)2][(R2N)3Th(m-15N)(Th(NR2)3] ([K][1-15N], R = SiMe3), and the thorium amide complex, [Th(NR2)3(15NH2)] (2-15N) were recorded. The spectrum for [K][1-15N] represents the first reported solid-state 15N NMR data for an actinide nitride complex. The experimentally measured tensor spans are 847 ppm for [K][1-15N] and 237 ppm for 2-15N. Both shielding tensors exhibit a near-zero asymmetry parameter, which for [K][1-15N] is consistent with a local rotational symmetry of its 15N-labelled nitride ligand. For 2-15N, the lack of asymmetry can be rationalized by a quasi-free Th-NH2 bond rotation in the solid-state. DFT calculations overestimate the tensor span somewhat for [K][1-15N], but provide isotropic shifts in good agreement with both the solid-state and solution values for both complexes. Natural localized molecular orbital (NLMO) analyses of the nuclear shielding reveal that the larger tensor span in [K][1-15N] vs. 2-15N is primarily a consequence of more pronounced covalency of the N-Th bonds, and large spin-orbit coupling due to significant Th 5f orbital contribution to those bonds, impacting the principal components of the shielding tensor perpendicular to the Th-N-Th axis. Overall, our analysis confirms the involvement of the 5f orbitals in Th-N multiple bonds, and further demonstrates the value of solid-state NMR spectroscopy for interrogating actinide-ligand bonding.

actinides, covalent bonding, Nuclear Magnetic Reso↗

Spectral Properties and Energy Injection in Mercury's Magnetotail Current Sheet

Mercury's magnetotail hosts a thin and highly dynamic current sheet (CS), where magnetic reconnection and strong fluctuations frequently occur. Here, we statistically analyze magnetic field power spectra across 370 magnetotail CSs observed by MESSENGER. About 20% of the events are quasi-laminar, showing single power-law spectra, whereas ∼80% are turbulent, exhibiting a spectral break separating inertial and kinetic ranges. A dawn–dusk asymmetry is identified: inertial-range slopes are systematically shallower on the dawnside, whereas kinetic-range slopes are steeper, indicating more developed turbulence there, consistent with the higher occurrence of reconnection-related processes on the dawnside. Component analysis shows that the transverse components, orthogonal to the tail-aligned principal field ( B X ), display shallow slopes near −1 in the inertial range, suggesting energy injection at ion scales rather than a classical inertial range. These results demonstrate that Mercury's unique plasma environment fundamentally reshapes the initiation of turbulence and the redistribution of energy in the magnetotail.

Mercury↗

Leveraging visible and near-infrared spectroelectrochemistry to calibrate a robust model for Vanadium(IV/V) in varying nitric acid and temperature levels

Spectroelectrochemistry and optimal design of experiments can be used to rapidly build accurate models for species quantification and enable a greater level of process awareness. Optical spectroscopy can provide vital elemental and molecular information, but several hurdles must be overcome before it can become a widely adopted analytical method for remote analysis in the nuclear field. Analytes with varying oxidation state, acid concentration, and fluctuating temperature must be efficiently accounted for to minimize time and resources in restrictive hot cell environments. The classic one-factor-at-a-time approach is not suitable for frequent calibration/maintenance operations in this setting. Therefore, a novel alternative was developed to characterize a system containing vanadium(IV/V) (0.01–0.1 M), nitric acid (0.1–4 M), and varying temperatures (20–45 °C). Here, spectroelectrochemistry methods were used to acquire a sample set selected by optimal design of experiments. This new approach allows for the accurate analysis of vanadium and HNO 3 concentration by leveraging UV–Vis–NIR absorption spectroscopy with robust and accurate chemometric models. The top model's root mean squared error of prediction percent values were 3.47%, 4.06%, 3.40%, and 10.9% for V(IV), V(V), HNO 3 , and temperature, respectively. These models, efficiently developed using the designed approach, exhibited strong predictive accuracy for vanadium and acid with varying oxidation states and temperature using only spectrophotometry, which advances current technology for real-world hot cell applications. Additionally, Nernstian analysis of the V(IV/V) standard potential was performed using traditional absorbance methods and multivariate curve resolution (MCR). The successful tests demonstrated that MCR Nernst tests may be valuable in highly convoluted spectral systems to better understand the redox processes' behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray absorption spectroscopy study of Mn reference compounds for Mn speciation in terrestrial surface environments

Abstract X-ray absorption spectroscopy (XAS) offers great potential to identify and quantify Mn species in surface environments by means of linear combination fit (LCF), fingerprint, and shell-fit analyses of bulk Mn XAS spectra. However, these approaches are complicated by the lack of a comprehensive and accessible spectrum library. Additionally, molecular-level information on Mn coordination in some potentially important Mn species occurring in soils and sediments is missing. Therefore, we investigated a suite of 32 natural and synthetic Mn reference compounds, including Mn oxide, oxyhydroxide, carbonate, phosphate, and silicate minerals, as well as organic and adsorbed Mn species, by Mn K-edge X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy. The ability of XAS to infer the average oxidation state (AOS) of Mn was assessed by comparing XANES-derived AOS with the AOS obtained from redox titrations. All reference compounds were studied for their local (<5 Å) Mn coordination environment using EXAFS shell-fit analysis. Statistical analyses were employed to clarify how well and to what extent individual Mn species (groups) can be distinguished by XAS based on spectral uniqueness. Our results show that LCF analysis of normalized XANES spectra can reliably quantify the Mn AOS within ~0.1 v.u. in the range +2 to +4. These spectra are diagnostic for most Mn species investigated, but unsuitable to identify and quantify members of the manganate and Mn(III)-oxyhydroxide groups. First-derivative XANES fingerprinting allows the unique identification of pyrolusite, ramsdellite, and potentially lithiophorite within the manganate group. However, XANES spectra of individual Mn compounds can vary significantly depending on chemical composition and/or crystallinity, which limits the accuracy of XANES-based speciation analyses. In contrast, EXAFS spectra provide a much better discriminatory power to identify and quantify Mn species. Principal component and cluster analyses of k2-weighted EXAFS spectra of Mn reference compounds implied that EXAFS LCF analysis of environmental samples can identify and quantify at least the following primary Mn species groups: (1) Phyllo- and tectomanganates with large tunnel sizes (2 × 2 and larger; hollandite sensu stricto, romanèchite, todorokite); (2) tectomanganates with small tunnel sizes (2 × 2 and smaller; cryptomelane, pyrolusite, ramsdellite); (3) Mn(III)-dominated species (nesosilicates, oxyhydroxides, organic compounds, spinels); (4) Mn(II) species (carbonate, phosphate, and phyllosilicate minerals, adsorbed and organic species); and (5) manganosite. All Mn compounds, except for members of the manganate group (excluding pyrolusite) and adsorbed Mn(II) species, exhibit unique EXAFS spectra that would allow their identification and quantification in mixtures. Therefore, our results highlight the potential of Mn K-edge EXAFS spectroscopy to assess bulk Mn speciation in soils and sediments. A complete XAS-based speciation analysis of bulk Mn in environmental samples should preferably include the determination of Mn valences following the “Combo” method of Manceau et al. (2012), EXAFS LCF analyses based on principal component and target transformation results, as well as EXAFS shell-fit analyses for the validation of LCF results. For this purpose, all 32 XAS reference spectra are provided in the Online Materials1 for further use by the scientific community.

Geochemistry & Geophysics↗

Global teleconnections influencing large-scale drought in the United States using SVDI

Understanding recent large-scale drought patterns and the mechanisms producing extreme drought events is vital for future drought forecasts and understanding future drought risks. Increasingly, vapor pressure deficit (VPD) has been used as an important measure of evaporative demand and proxy for drought detection. In this study, VPD is used to calculate the new Standardized VPD Drought Index (SVDI) with NASA North American Land Data Assimilation System (NLDAS) data. Previous studies have shown that SVDI accurately identifies the timing and magnitude short-term droughts in the United States (U.S). In the present study, SVDI is now used to identify large-scale drought patterns between 1980 and 2021 and drought variability driven by selected global teleconnections originating in the Pacific and Atlantic Oceans. Spatial drought characteristics were extracted from SVDI using empirical orthogonal function (EOF) analysis. Then a k-means clustering algorithm was applied to both EOF principal components and primary teleconnections, including the El Nino-Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO) to identify drought events driven by the Pacific Ocean. Results show that the SVDI is useful in evaluating large-scale drought variability in the U.S. related to global teleconnections, and that mechanisms influencing summer drought patterns in the Western and Southwestern U.S. are driven by a tropical-extratropical interactions originating in the equatorial Pacific Ocean related to ENSO dynamics with interdecadal variability modulated by PDO. The large-scale droughts in the Central and Southern U.S., like those in 2011 and 2012, on the other hand, are driven by the North Pacific Ocean warm pool during a strong negative PDO, which subsequently influenced variability in the Bermuda-Azores High in the Atlantic Ocean. In summer 2011, the Bermuda-Azores High weakened, reducing the onshore winds and moisture transport along the eastern Gulf of Mexico and contributing to ongoing drought in the region. The Northern Pacific and Atlantic Ocean sea surface temperatures (SSTs) have increased between 1980 and 2021. In conclusion, as SSTs continue to rise in the Northern Pacific Ocean, one consequence of the coupled North Pacific warm pool and atmospheric dynamics, is to increase summer drought variability over a large region in the southern and midwestern U.S. under global warming.

54 ENVIRONMENTAL SCIENCES↗

pnnl/Chemometric_Toolbox

The PNNL Chemometric Toolbox is a software collection of common MATLAB scripts that implement core chemometric algorithms for regression analysis. The three core regression techniques within this toolbox are classical least squares (CLS), principal component regression (PCR), and partial least squares (PLS). This is demonstrated using a supplied dataset of infrared (FTIR) spectral data with their corresponding concentrations

Smith, Ian↗

GOOML - Finding Optimization Opportunities for Geothermal Operations: Preprint

Geothermal Operational Optimization with Machine Learning (GOOML) is a transferable and extensible component-based geothermal asset modeling framework that considers complex steamfield relationships and identifies optimization prospects using a data-driven approach. We have used this framework to develop digital twins that provide steamfield operators with an operational environment to analyze and understand historical and forecasted power production, explore new steamfield configuration possibilities, and seek optimal asset management for real world applications. The GOOML modeling software is built on a generic component-based systems framework that allows for both historical and forecast analysis. A GOOML model can perform historical data-assimilation using first-principal thermodynamics to create a meaningful data model. Historical production data can then be coupled with a forecast framework to train machine-learning models of steamfield components to predict future outputs. This modeling environment enables digital exploration of steamfield design configurations and operational scenarios. GOOML digital twins have been developed for steamfields in New Zealand and the United States representing differing power generation and field conditions. These digital twins have been validated by comparing hindcast predictions against historical production data. Reinforcement learning experiments were conducted to demonstrate the ability to programmatically explore the operations space using machine learning agents. Our initial results are compelling; two to five percent increases in annual energy production were demonstrated by the GOOML models with no additional infrastructure build required. GOOML offers a new approach to geothermal operations by applying state-of-the-art machine learning algorithms, comprehensive data analytics, and interaction with digital twins. Through application of these tools, operators will realize greater availability and higher net generation which will increase the cost effectiveness of geothermal energy projects.

access↗

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗

Isotopic Signatures of Lithium Carbonate and Lithium Hydroxide Monohydrate Measured Using Raman Spectroscopy

Lithium isotopic ratios have wide ranging applications as chemical signatures, including improved understanding of geochemical processes and battery development. Measurement of isotope ratios using optical spectroscopies would provide an alternative to traditional mass spectrometric methods, which are expensive and often limited to a chemical laboratory. In this work, Raman spectra of 7 Li 2 CO 3 , 6 Li 2 CO 3 , 7 LiOH*H 2 O and 6 LiOH*H 2 O have been measured to determine the effect of lithium isotope substitution on the Raman molecular vibrations. Thirteen peaks were observed in the spectrum of lithium carbonate, with discernable isotopic shifts occurring in eleven of the thirteen vibrations, two of which have not been previously reported in the literature. The spectrum of lithium hydroxide monohydrate contained nine peaks, with discernable isotopic shifts occurring in eight of the nine vibrations, four of which have not been previously reported in the literature. The Raman spectral data reported here for lithium carbonate and lithium hydroxide monohydrate are in agreement with the previously reported works in the literature, in which the Raman active modes of these molecules were first identified and assigned. However, due to the stability and resolution of the detection system used in this work, isotopic shifts with a magnitude less than one wavenumber have been identified. Principal Component Regression was used to evaluate the sensitivity to isotopic content of small Raman peak shifts in Li 2 CO 3 and indicates differences greater than 2 atom-% could be reliably determined. These measurements add to the body of work on lithium isotope Raman spectroscopy for these two compounds and increases the number of Raman bands which could be used for lithium isotope content analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗

Unsupervised Machine Learning for Exploratory Data Analysis of Exoplanet Transmission Spectra

Abstract Transit spectroscopy is a powerful tool for decoding the chemical compositions of the atmospheres of extrasolar planets. In this paper, we focus on unsupervised techniques for analyzing spectral data from transiting exoplanets. After cleaning and validating the data, we demonstrate methods for: (i) initial exploratory data analysis, based on summary statistics (estimates of location and variability); (ii) exploring and quantifying the existing correlations in the data; (iii) preprocessing and linearly transforming the data to its principal components; (iv) dimensionality reduction and manifold learning; (v) clustering and anomaly detection; and (vi) visualization and interpretation of the data. To illustrate the proposed unsupervised methodology, we use a well-known public benchmark data set of synthetic transit spectra. We show that there is a high degree of correlation in the spectral data, which calls for appropriate low-dimensional representations. We explore a number of different techniques for such dimensionality reduction and identify several suitable options in terms of summary statistics, principal components, etc. We uncover interesting structures in the principal component basis, namely well-defined branches corresponding to different chemical regimes of the underlying atmospheres. We demonstrate that those branches can be successfully recovered with a K-means clustering algorithm in a fully unsupervised fashion. We advocate for lower-dimensional representations of the spectroscopic data in terms of the main principal components, in order to reveal the existing structure in the data and quickly characterize the chemical class of a planet.

Matchev, Konstantin T. (ORCID:0000000341829096)↗