Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Alternating least squares”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Confocal Raman imaging of skin sections containing hair follicles using classical least squares regression and multivariate curve resolution – alternating least squares

Confocal Raman microscopy (CRM) is applied ex vivo for imaging of the spatial distribution of different skin components in skin sections containing hair follicles. For multivariate data analysis, different methods are used in order to spectrally decompose the reference spectra of the skin components (dermis, viable epidermis, stratum corneum and hair). Classical least squares regression (CLS) and multivariate curve resolution–alternating least squares (MCR-ALS) are chosen as suitable methods. In comparison to other optical methods, the advantage of CRM is molecular specificity and dispensability of labelling dyes, which is e.g. necessary in fluorescence microscopy. Therefore, a useful future application of CRM in combination with multivariate data analysis lies in the analysis of penetration routes of topically applied substances, such as cosmetic formulations or drugs into the skin, which is particularly interesting in and around hair follicles. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Efficient CP Rounding Using Alternating Least Squares with QR Decomposition

The CANDECOMP/PARAFAC (CP) decomposition is widely used for analyzing multidimensional data, and the alternating least squares (CP-ALS) algorithm is a common method for its computation. CP rounding is the problem of computing a lower-rank CP decomposition of an input already in a higher-rank CP format. While the normal equations (NE) approach in CP-ALS is efficient for the CP rounding problem and frequently used, it becomes unstable in the presence of ill-conditioned subproblems. This paper presents a new QR-based CP-ALS method for CP rounding that preserves both numerical stability and computational efficiency. Here, our experiments show that the proposed method offers significant speedup over a previous QR-based approach and the Tensor Toolbox's NE-based implementation, particularly for higher-order tensors. Furthermore, our approach demonstrates a marked reduction in error for ill-conditioned problems, with error reductions several orders of magnitude smaller compared to the NE-based method, while achieving faster convergence and more accurate solutions. By using a more numerically stable approach, we can solve more problems in reduced working precision, which enables further reduction in time to solution.

CANDECOMP/PARAFAC↗
DOE OSTI2023

CP decomposition for tensors via alternating least squares with QR decomposition

The CP tensor decomposition is used in applications such as machine learning and signal processing to discover latent low-rank structure in multidimensional data. Computing a CP decomposition via an alternating least squares (ALS) method reduces the problem to several linear least squares problems. The standard way to solve these linear least squares subproblems is to use the normal equations, which inherit special tensor structure that can be exploited for computational efficiency. However, the normal equations are sensitive to numerical ill-conditioning, which can compromise the results of the decomposition. In this paper, we develop versions of the CP-ALS algorithm using the QR decomposition and the singular value decomposition, which are more numerically stable than the normal equations, to solve the linear least squares problems. Our algorithms utilize the tensor structure of the CP-ALS subproblems efficiently, have the same complexity as the standard CP-ALS algorithm when the input is dense and the rank is small, and are shown via examples to produce more stable results when ill-conditioning is present. Our MATLAB implementation achieves the same running time as the standard algorithm for small ranks, and we show that the new methods can obtain lower approximation error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Application of Multivariate Curve Resolution–Alternate Least Square Technique on Extracting Pure Spectral Components from Multiple Emitting Systems: a Case Study

Multiple emitting components in a fluorophoric system often produce complicated emission spectra. Extracting the individual spectral information from the composite spectra is important in order to comprehend the photophysical processes occurring in the multifluorophoric systems. Although the combination of Principal Component Analysis and Multivariate Curve Resolution-Alternate Least Square (PCA/MCR-ALS) technique is a well-known approach for curve deconvolution, its applicability in the spectral deconvolution of vibronically and electronically mixed up emitting systems as well as systems merged up with multiple electronic bands without a priori knowledge of the individual emitting species, is not properly studied. The present work highlights the strength of PCA/MCR-ALS in retrieving pure spectral information from the set of complex spectra arising out of the regular variation of causative factors that result in the variation of spectral composition. The retrieval of the emission bands utilizing the PCA/MCR-ALS technique has been made without having a priori information of the emitting species present in the multifluorophoric systems and the resolved spectra correspond well with the fluorescence spectra of the individual chemical species. The common curve fitting methods such as Gaussian and Lorentzian techniques have been found to be unsuccessful in providing meaningful photophysical information through the retrieved spectra. A comparative study of the curve fitting techniques MCR-ALS, Gaussian and Lorentzian in a set of complicated emission spectra of (i) pyrene and its excimer, (ii) pyrene and its excimer in presence of benzo[a]pyrene, and (iii) fisetin in bile salt medium is presented herein in details.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Use of Fisher's Ratio assisted multivariate curve resolution- alternating least squares for discovery-based analysis using ultrahigh pressure liquid chromatography-high resolution mass spectrometry

Non-targeted analysis of complex chemical mixtures can be difficult considering the convoluted nature of the matrix and the potential unknown chemical differences between samples or classes of samples. Ultrahigh pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOF) is an ideal technique to probe chemical differences for a wide variety of samples. While UHPLC-QTOF can discover minute chemical differences down to low part per billion (ppb) concentrations with a high degree of confidence, the application of high-resolution mass spectrometry can yield massive amounts of information (∼ 10 gb per sample) that cannot be analyzed manually. Therefore, the application of chemometric techniques is mandatory for the interrogation of complex samples. Fisher's ratio (FR) assisted multivariate curve resolution-alternating least squares (MCR-ALS) was used to the discover and identify the chemical differences between two classes of materials: 1) a pond water matrix and 2) the matrix spiked with a pharmaceutical standard mix containing 17 compounds. Thirteen of the seventeen spiked compounds were discovered using FR analysis, and then five were successfully deconvoluted using MCR-ALS wherein the number of curves chosen were automatically determined using singular value decomposition (SVD). In conclusion, the use of an automated FR assisted MCR-ALS will aid in discovering trace levels of chemical components without the need for the researcher to provide potentially biased input which will aid in non-targeted workflow.

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

The controlling role of atmosphere in dawsonite versus gibbsite precipitation from tetrahedral aluminate species

In highly alkaline solution, aluminum speciates as the tetrahedrally coordinated aluminate monomer, Al(OH)4- and/or dimer Al2O(OH)62-, yet precipitates as octahedrally coordinated gibbsite (Al(OH)3). This tetrahedral to octahedral transformation governs Al precipitation, which is crucial to worldwide Al production, and to the processing of caustic high-level radioactive wastes. Despite its significance, the transformation pathway remains unknown. Here we explore the roles of atmospheric water and carbon dioxide in mediating the transformation of the tetrahedrally coordinated potassium aluminate dimer salt (K2Al2O(OH)6) to gibbsite versus potassium dawsonite (KAl(CO3)(OH)2). A combination of in-situ attenuated total reflection infrared spectroscopy, ex-situ micro X-ray diffraction, and multivariate curve resolution-alternating least squares chemometrics analysis reveals that humidity plays a key role in the transformation by limiting the amount of alkalinity neutralization by dissolved CO2. Lower humidity favors higher alkalinity and incorporation of carbonate species in the final Al product to form KAl(CO3)(OH)2. Higher humidity enables more acid generation that destabilizes dawsonite and favors gibbsite as the solubility limiting phase. Because the transformation was restricted to occur in thin water films, the results suggest that transition from tetra- to octahedrally coordinated Al does not have to occur in bulk solution, as has often been hypothesized, but may instead appear on the source mineral surface.

Dembowski, Mateusz↗
DOE OSTI2019

Interpreting mixing relationships in energetic melts to estimate vapor contribution and composition

Heterogeneous, glassy debris can be produced by energetic events such as near surface nuclear explosions and asteroid impacts. In these events, the resultant debris comprises a rapidly-cooled mixture of the projectile or bomb material with proximate materials such as soil and anthropogenic structures. Glasses formed in these events preserve inter and/or intra sample compositional variation reflecting complex mixtures of the precursor materials or their chemically fractionated products. Understanding formation mechanisms and chemical processes relating to energetic glass formation based on study of the end products requires an objective method for calculating the precursor components that mixed together to produce the final product. Here, we compare a set of multivariate linear-least squares approaches to address the unmixing problem of heterogeneous glasses including classical least squares (CLS), principal component analysis (PCA) and two multivariate curve resolution – alternating least squares (MCR-ALS) approaches. We generate a synthetic data set representing spatially resolved, compositional data in a set of heterogeneously mixed glasses to understand model sensitivity to precursor abundance and composition. We explore the impact of incomplete knowledge of the system prior to the formation of the melts on the resulting confidence and accuracy of the solution (i.e., the number, composition, and abundances of the precursor materials). We conclude that a closure and equality constrained MCR-ALS approach succeeds in identifying unrecognized precursor components and is a suitable alternative to CLS approaches when the investigator has incomplete knowledge of the starting precursor chemical compositions. Finally, we demonstrate application of this MCR-ALS approach to measurements of Trinitite.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2017

Gradient-based Optimization for Regression in the Functional Tensor-Train Format

We consider the task of low-multilinear-rank functional regression, i.e., learning a low-rank parametric representation of functions from scattered real-valued data. Our first contribution is the development and analysis of gradient-based optimization procedures, including stochastic gradient descent and quasi-Newton methods, for learning the parameters of a functional tensor-train (FT) that yields improved accuracy over standard alternating least squares methods. The functional tensor-train uses the tensor-train (TT) representation of low-rank arrays as an ansatz for a class of low-multilinear-rank functions. The FT is represented by a set of matrix-valued functions that contain a set of univariate functions, and the regression task is to learn the parameters of these univariate functions. Our second contribution demonstrates that using nonlinearly parameterized univariate functions, e.g., symmetric kernels with moving centers, within each core can outperform the standard approach of using a linear expansion of basis functions. Our final contributions are new rank adaptation and group-sparsity regularization procedures to minimize overfitting. We use several benchmark problems to demonstrate at least an order of magnitude lower accuracy with gradient-based optimization methods than ALS in the low-sample number regime. We also demonstrate an order of magnitude reduction in accuracy on a test problem resulting from nonlinear parameterizations over linear parameterizations. Finally we compare regression performance with 22 other nonparametric and parametric regression methods on 10 real-world data sets. We achieve top-five accuracy for seven of the data sets and best accuracy for two of the data sets. These rankings are the best amongst parametric models and competetive with the best non-parametric methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

The Electrochemical Sodiation of FeSb 2 : New Insights from Operando 57 Fe Synchrotron Mössbauer and X-Ray Absorption Spectroscopy

Operando Sb K–edge X–ray absorption spectroscopy and 57 Fe Synchrotron Mössbauer Spectroscopy, used for the first time in the field of operando energy storage materials, assisted by operando magnetic measurements, were combined to clarify the role of iron and antimony in the electrochemical reaction mechanism of FeSb 2 as negative electrode material for sodium–ion batteries (SIB). Both datasets were analyzed using an innovative chemometric approach involving principal component analysis (PCA) and multivariate curve resolution – alternating least square analysis (MCR–ALS) yielding new insights on the sodiation reaction. Our findings show that the reaction of Na with FeSb 2 during discharge leads to the formation of Na 3 Sb along with superparamagnetic Fe amorphous nanoparticles which contain small amounts of Sb dissolved in their lattice. During the following desodiation, the pristine material FeSb 2 is not recovered while iron nanoparticles grow in size, and continue growing also along the following discharge. Here, even though such iron nanoparticles remain electrochemically inactive, they play a key role in the reduction and stabilization of the polarization as well as in the reversibility of the electrochemical sodiation of antimony.

25 ENERGY STORAGE↗
DOE OSTI2020

The Fe 4+ / 3+ Redox Mechanism in NaFeO 2 : A Simultaneous Operando Nuclear Resonance and X-ray Scattering Study

Simultaneous operando Nuclear Forward Scattering and transmission X-ray diffraction and 57 Fe Mössbauer spectroscopy measurements were carried out in order to investigate the electrochemical mechanism of NaFeO 2 vs. Na metal using a specifically designed in situ cell. The obtained data were analysed using an alternative and innovative data analysis approach based on chemometric tools such as Principal Component Analysis (PCA) and Multivariate Curve Resolution - Alternating Least Squares (MCR-ALS). This approach, which allows the unbiased extraction of all possible information from the operando data, enabled the stepwise reconstruction of the independent “real” components permitting the description of the desodiation mechanism of NaFeO 2 . This wealth of information allows a clear description of the electrochemical reaction at the redox-active iron centres, and thus an improved comprehension of the cycling mechanisms of this material vs. sodium.

25 ENERGY STORAGE↗
DOE OSTI2017

HPC formulations of optimization algorithms for tensor completion

Tensor completion is a powerful tool used to estimate or recover missing values in multi-way data. It has seen great success in domains such as product recommendation and healthcare. Tensor completion is most often accomplished via low-rank sparse tensor factorization, a computationally expensive non-convex optimization problem which has only recently been studied in the context of parallel computing. In this work, we study three optimization algorithms that have been successfully applied to tensor completion: alternating least squares (ALS), stochastic gradient descent (SGD), and coordinate descent (CCD++). We explore opportunities for parallelism on shared- and distributed-memory systems and address challenges such as memory- and operation-efficiency, load balance, cache locality, and communication. Among our advancements are a communication-efficient CCD++ algorithm, an ALS algorithm rich in level-3 BLAS routines, and an SGD algorithm which combines stratification with asynchronous communication. Furthermore, we show that introducing randomization during ALS and CCD++ can accelerate convergence. We evaluate our parallel formulations on a variety of real datasets on a modern supercomputer and demonstrate speedups through 16384 cores. Further, these improvements reduce time-to-solution from hours to seconds on real-world datasets. We show that after our optimizations, ALS is advantageous on parallel systems of small-to-moderate scale, while both ALS and CCD++ provide the lowest time-to-solution on large-scale distributed systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

New methods for trace analysis of gamma-irradiated pentaerythritol tetranitrate

High explosives (HEs) are used in a diverse range of applications in which they could be exposed to various radiation levels that may cause potential chemical changes. This study further evaluated pentaerythritol tetranitrate (PETN) that was previously aged with a low-level 2 kGy dose of gamma irradiation in order to understand chemical changes caused by irradiation. Both unirradiated PETN and gamma-irradiated PETN were analyzed using ultra high-pressure liquid chromatography coupled to quadrupole time of flight mass spectrometry (UHPLC-QTOF). The resulting data were processed in a non-targeted manner using Fisher's ratio analysis and multivariate curve resolution-alternating least squares (MCR-ALS) to aid in discovery and identification of the chemical changes brought about by irradiation without a priori knowledge. The application of using UHPLC-QTOF in combination with chemometric techniques for the analysis of irradiated samples has not previously been performed. Here, in this work, we show how to use this method to provide chemical information that would otherwise not be discernible, such as the discovery of degradation of the various homologues of PETN. Major differences identified with radiolytic aging of the PETN sample included decomposition products that resulted from the degradation of the trigger linkage – the O–NO 2 bonds – resulting in the formation of alcohol and aldehyde groups. Similar degradation was also observed in the PETN homologues as well as interconversion from one homologue to another.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Temporal Ni K-Edge X-ray Absorption Spectroscopy Study Reveals the Kinetics of the Ni Redox Behavior of the Iron-Nickel Oxide Bimetallic OER Catalyst

Operando X-ray absorption spectroscopy (XAS) can be utilized to probe the phase and structural changes of FeNiOx and similar transition metal oxide electrocatalysts during electrocatalytic reactions. However, capturing the temporal changes occurring in the operando chemistry of electrocatalysts has been little studied. Here, in this work, we successfully capture the time-resolved changes at the Ni K-edge for both the X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) regions of an aqueous phase synthesized FeNiO x bimetallic nanoparticle electrocatalyst. During a stepped voltage experiment, the temporal change in the Ni K-edge was observed as the applied voltage was stepped from 0.7 to 0.8 V versus the silver/silver chloride reference electrode, and the change observed is associated with the Ni redox transition from the 2+ to 3+/4+ oxidation state. Individual XAS spectra were obtained in 90 s, and a total of 9 scans post voltage step showed a unique transition from a hydroxide to an oxyhydroxide phase. The shift in absorption edge energy position and the changes in spectral shape for individual scans explain the typically observed broadening of the Ni K-edge XANES spectrum during time-averaged operando XAS, and a kinetic analysis revealed a first-order observed rate constant, |k obs |, of 0.00426 s –1 and a half-life (t 1/2 ) of 163 s. Multivariate curve resolution-alternating least squares analysis followed by linear combination fitting analysis for time-averaged versus time-resolved Ni K-edge changes upon voltage step show distinct differences in estimated contributions from Ni 2+ oxide/hydroxide phases. Detailed EXAFS modeling shows the phase transition from hydroxide to oxyhydroxide on top of the unchanged metallic core.

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

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↗
DOE OSTI2019

Hydration of ferric chloride and nitrate in aqueous solutions: water-mediated ion pairing revealed by Raman spectroscopy

Iron is the most abundant transition metal in the earth's crust and is important for the proper functioning of many technological and natural processes. Despite the importance, a complete microscopic understanding of the hydration of ferric ions and water mediated ion pairing has not been realized. Hydrated Fe(III) is difficult to study due to the process of complexation to the anion and hydrolysis of the hydrating water molecules leading to a heterogeneous solution with diverse speciation. Here, ferric chloride and nitrate aqueous solutions are studied using polarized Raman spectroscopy as a function of concentration and referenced to their respective sodium salt or mineral acid. Perturbed water spectra (PWS) were generated using multivariate curve resolution-alternating least squares (MCR-ALS) to show the residual spectral response uniquely attributable to the hydration of ferric speciation. The hydrogen bonding network associated with the hydrating water molecules in ferric chloride solutions are found to be more similar to hydrochloric acid solutions, whereas in ferric nitrate solutions, the network behaves more similar to sodium nitrate, despite increased acidity. Thus, in the FeNO 3 and FeCl 3 solutions, ion pairing and coordination, respectively, are significantly influencing the hydration spectra signature. These results further reveal concentration dependent changes to the hydrogen bonding network, hydrating water symmetry, and changes to the relative abundance of solvent shared ion pairs that are governed primarily by the ferric salt identity.

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

Operando structural and chemical evolutions of TiS 2 in Na-ion batteries

Titanium disulfide (TiS 2 ) with high electric conductivity, fast rate capability, and good cycling performance is a promising candidate for electrode material in sodium (Na)-ion batteries. Despite the well-studied electrochemical behaviors of TiS 2 in Li-ion batteries, the detailed reaction mechanism of TiS 2 in Na-ion batteries is not yet fully understood due to a more complex multi-phase conversion process. In this work, reactions of TiS 2 in Na-ion batteries are investigated via a multi-modal synchrotron approach: operando X-ray Absorption Spectroscopy (XAS) – including X-ray Absorption Near-Edge Structure (XANES) and Extended X-ray Absorption Fine Structure (EXAFS) – and ex situ X-ray Powder Diffraction (XPD), coupled with computational modeling. Operando XANES spectra indicate that the redox reactions occur in both Ti and S during the electrochemically driven phase transformation. Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) analysis of XAS suggests that different numbers of components are involved in the lithiation and sodiation of TiS 2 , with the sodiation including at least one intermediate phase in addition to the starting material and the final sodiation product. Ex situ XPD and Rietveld refinement further determined and quantified the unknown phases, showing that three phases, TiS 2 , Na 0.55 TiS 2 , and NaTiS 2 , participate in the sodiation of TiS 2 . Operando EXAFS results show the changes in the Ti–Ti coordination number and interatomic distance. This explains the coulombic efficiency decay due to the incomplete recovery of the coordination number of Ti after cycling. Overall, this work reveals the reaction mechanism occurring in Na–TiS 2 batteries with a greater quantitative understanding of the structural evolution. By combining the multi-modal synchrotron approach and computational work, this study provides a framework for studying a broader range of electrochemically driven phase-transformation systems towards advanced energy storage and conversion applications.

25 ENERGY STORAGE↗
DOE OSTI2017

Low-rank canonical-tensor decomposition of potential energy surfaces: application to grid-based diagrammatic vibrational Green's function theory

Here, a new method is proposed for a fast evaluation of high-dimensional integrals of potential energy surfaces (PES) that arise in many areas of quantum dynamics. It decomposes a PES into a canonical low-rank tensor format, reducing its integral into a relatively short sum of products of low-dimensional integrals. The decomposition is achieved by the alternating least squares (ALS) algorithm, requiring only a small number of single-point energy evaluations. Therefore, it eradicates a force-constant evaluation as the hotspot of many quantum dynamics simulations and also possibly lifts the curse of dimensionality. This general method is applied to the anharmonic vibrational zero-point and transition energy calculations of molecules using the second-order diagrammatic vibrational many-body Green's function (XVH2) theory with a harmonic-approximation reference. In this application, high dimensional PES and Green's functions are both subjected to a low-rank decomposition. Evaluating the molecular integrals over a low-rank PES and Green's functions as sums of low-dimensional integrals using the Gauss–Hermite quadrature, this canonical-tensor-decomposition-based XVH2 (CT-XVH2) achieves an accuracy of 0.1 cm -1 or higher and nearly an order of magnitude speedup as compared with the original algorithm using force constants for water and formaldehyde.

74 ATOMIC AND MOLECULAR PHYSICS↗