Search NASA⌕ Search

SEARCH · Search NASA

Results for “Generalized Linear Model”

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 235 records · Page 13

Stability of Schwarzschild-like solutions in f(R,G) gravity models

We study linear metric perturbations around a spherically symmetric static spacetime for general f(R,G) theories, where R is the Ricci scalar and G is the Gauss-Bonnet term. We find that, unless the determinant of the Hessian of f(R,G) is zero, even-type perturbations have a ghost for any multipole mode. In order for these theories to be plausible alternatives to general relativity, the theory should satisfy the condition that the ghost is massive enough to effectively decouple from the other fields. We study the requirement on the form of f(R,G) which satisfies this condition. We also classify the number of propagating modes both for the odd-type and the even-type perturbations and derive the propagation speeds for each mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗

A closer look at turbulence spreading: How bistability admits intermittent, propagating turbulence fronts

In magnetic fusion plasmas, mounting evidence indicates the possibility of sustained turbulence below the linear stability threshold or more generally global turbulence bistability. The usual reduced models for turbulence spreading are unistable/supercritical and incompatible with this result. The older models further cannot realistically support fronts connecting laminar and turbulent domains. In this work, a minimal model for “subcritical” turbulence spreading is introduced and analyzed. The model may be viewed as phenomenological or derived directly by considering the effect of profile corrugations in an E×B staircase. The model, which is related to the FitzHugh–Nagumo system, supports the robust coexistence of multiple turbulence levels via bistability. We show that this model predicts stronger penetration of turbulence into a linearly stable region as well as the formation of intermittent turbulence fronts that resemble avalanches. We derive the critical size that a localized slug of turbulence must exceed in order to spread. Lastly, we make a prediction of global hysteretic behavior associated with the bistability, which should be testable via experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computation of generalized and exact dispersion relations for longitudinal plasma waves in nonextensive statistics and the effects of the nonextensivity on the oscillation modes and damps

We have derived generalized dispersion relations for longitudinal waves in collisionless thermal plasma using linear Vlasov-Poisson kinetic model and nonextensive distributions for electrons. The Maxwellian limit of the dispersion relations, where the q-nonextensive parameter tends to one, is calculated. The generalized dispersion relations are reduced to polynomials for some specific values of q. The well-known modes of oscillations such as the Langmuir and electron acoustic waves have been obtained by solving the dispersion relations. Some new modes of oscillation are also found. Finally, the dependence of the oscillation modes and damps on q is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MITgcm-AD v2: Open source tangent linear and adjoint modeling framework for the oceans and atmosphere enabled by the Automatic Differentiation tool Tapenade

The Massachusetts Institute of Technology General Circulation Model (MITgcm) is widely used by the climate science community to simulate planetary atmosphere and ocean circulations. A defining feature of the MITgcm is that it has been developed to be compatible with an algorithmic differentiation (AD) tool, TAF, enabling the generation of tangent-linear and adjoint models. These provide gradient information which enables dynamics-based sensitivity and attribution studies, state and parameter estimation, and rigorous uncertainty quantification. Importantly, gradient information is essential for computing comprehensive sensitivities and performing efficient large-scale data assimilation, ensuring that observations collected from satellites and in-situ measuring instruments can be effectively used to optimize a large uncertain control space. As a result, the MITgcm forms the dynamical core of a key data assimilation product employed by the physical oceanography research community: Estimating the Circulation and Climate of the Ocean (ECCO) state estimate. Although MITgcm and ECCO are used extensively within the research community, the AD tool TAF is proprietary and hence inaccessible to a large proportion of these users. The new version 2 (MITgcm-AD v2) framework introduced here is based on the source-to-source AD tool Tapenade, which has recently been open-sourced. Another feature of Tapenade is that it stores required variables by default (instead of recomputing them) which simplifies the implementation of efficient, AD-compatible code. The framework has been integrated with the MITgcm model’s main branch and is now freely available.

Adjoints↗

Efficient Automated Driving Strategies Leveraging Anticipation and Optimal Control

Automated vehicles and advanced driver assistance systems bring computation, sensing, and communication technologies that exceed human abilities in some ways. For example, automated vehicles may sense a panorama all at once, do not suffer from human impairments and distractions, and could wirelessly communicate precise data with neighboring vehicles. Prototype and commercial deployments have demonstrated the capability to relieve human operators of some driving tasks up to and including fully autonomous taxi rides in some areas. The ultimate impact of this technology’s large-scale market penetration on energy efficiency remains unclear, with potential negative factors like road use by empty vehicles competing with positive ones like automatic eco-driving. Fundamentally enabled by historic and look-ahead data, this dissertation addresses the use of automated driving and driver assistance to optimize vehicle motion for energy efficiency. Facets of this problem include car following, co-optimized acceleration and lane change planning, and collaborative multi-agent guidance. Optimal control, especially model predictive control, is used extensively to improve energy efficiency while maintaining safe and timely driving via constraints. Techniques including chance constraints and mixed integer programming help overcome uncertainty and non-convexity challenges. Extensions of these techniques to tractor trailers on sloping roads are provided by making use of linear parameter-varying models. To approach the wheel-input energy eco-driving problem over generally shaped sloping roads with the computational potential for closed-loop implementation, a linear programming formulation is constructed. Distributed and collaborative techniques that enable connected and automated vehicles to accommodate their neighbors in traffic are also explored and compared to centralized control. Using simulations and vehicle-in-the-loop car following experiments, the proposed algorithms are benchmarked against others that do not make use of look-ahead information.

Dollar, Robert Austin↗

Ultraviolet ISM diagnostics for star-forming galaxies. I. Tracers of metallicity and extinction

We have observed a sample of 14 nearby (z∼0.03) star-forming blue compact galaxies (BCGs) in the rest-frame far-UV (∼1150–2200 Å) using the Cosmic Origins Spectrograph on the Hubble Space Telescope. We have also generated a grid of stellar population synthesis models using the Starburst99 evolutionary synthesis code, allowing us to compare observations and theoretical predictions for the Si iv-1400 and C iv-1550 UV indices; both are comprised of a blend of stellar wind and interstellar lines and have been proposed as metallicity diagnostics in the UV. Our models and observations both demonstrate that there is a positive linear correlation with metallicity for both indices, and we find generally good agreement between our observations and the predictions of the Starburst99 models (with the models slightly under-estimating the value of the indices due to contributions from interstellar lines not simulated by a stellar population synthesis code). By combining the rest-frame UV observations with pre-existing rest-frame optical spectrophotometry of our BCG sample, we also directly compare the predictions of metallicity and extinction diagnostics across both wavelength regimes. This comparison reveals a correlation between the UV absorption and optical strong-line diagnostics, offering the first means of directly comparing interstellar medium (ISM) properties determined across different rest-frame regimes. Finally, using our Starburst99 model grid, we determine theoretical values for the short-wavelength UV continuum slope, β{sub 18}, which can be used for determining extinction in rest-frame UV spectra of star-forming galaxies. We consider the implications of these results and discuss future work aimed at parameterizing these and other environmental diagnostics in the UV (a suite of diagnostics that could offer particular utility in the study of star-forming galaxies at high redshift) as well as the development of robust comparisons between ISM diagnostics across a broad wavelength baseline.

79 ASTRONOMY AND ASTROPHYSICS↗

Perturbative Stability and Error-Correction Thresholds of Quantum Codes

Topologically ordered phases are stable to local perturbations, and topological quantum error-correcting codes enjoy thresholds to local errors. We connect the two notions of stability by constructing classical statistical mechanics models for decoding general Calderbank-Shor-Steane codes and classical linear codes. Our construction encodes correction success probabilities under uncorrelated bit-flip and phase-flip errors, and simultaneously describes a generalized ℤ 2 lattice-gauge theory with quenched disorder. We observe that the clean limit of the latter is precisely the discretized imaginary-time path integral of the corresponding quantum code Hamiltonian when the errors are turned into a perturbative 𝑋 or 𝑍 magnetic field. Motivated by error-correction considerations, we define general order parameters for all such generalized ℤ 2 lattice-gauge theories, and show that they are generally lower bounded by success probabilities of error correction. For CSS codes satisfying the low-density parity-check condition and with a sufficiently large code distance, we prove the existence of a low-temperature ordered phase of the corresponding lattice-gauge theories, particularly for those lacking Euclidean spatial locality and/or when there is a nonzero code rate. We further argue that these results provide evidence for stable phases in the corresponding perturbed quantum Hamiltonians, obtained in the limit of continuous imaginary time. To do so, we distinguish space- and timelike defects in the lattice-gauge theory. A high free-energy cost of spacelike defects corresponds to a successful “memory experiment” and suppresses the energy splitting among the ground states, while a high free-energy cost of timelike defects corresponds to a successful “stability experiment” and points to a nonzero gap to local excitations.

quantum error correction↗

Learning genetic perturbation effects with variational causal inference

Advances in sequencing technologies have enhanced the understanding of gene regulation in cells. In particular, Perturb-seq has enabled high-resolution profiling of the transcriptomic response to genetic perturbations at the single-cell level. This understanding has implications in functional genomics and potentially for identifying therapeutic targets. Various computational models have been developed to predict perturbational effects. While deep learning models excel at interpolating observed perturbational data, they tend to overfit in the lack of enough data and may not generalize well to unseen perturbations. In contrast, mechanistic models, such as linear causal models based on gene regulatory networks, hold greater potential for extrapolation, as they encapsulate regulatory information that can predict responses to unseen perturbations. However, their application has been limited to small studies due to overly simplistic assumptions, making them less effective in handling noisy, large-scale single-cell data. We propose a hybrid approach that combines a mechanistic causal model with variational deep learning, termed Single Cell Causal Variational Autoencoder (SCCVAE). The mechanistic model employs a learned regulatory network to represent perturbational changes as shift interventions that propagate through the learned network. SCCVAE integrates this mechanistic causal model into a variational autoencoder, generating rich, comprehensive transcriptomic responses. Our results indicate that SCCVAE exhibits superior performance over current state-of-the-art baselines for extrapolating to predict unseen perturbational responses. Additionally, for the observed perturbations, the latent space learned by SCCVAE allows for the identification of functional perturbation modules and simulation of single-gene knockdown experiments of varying penetrance, presenting a robust tool for interpreting and interpolating perturbational responses at the single-cell level.

59 BASIC BIOLOGICAL SCIENCES↗

A conformal approach for the analysis of the non-linear stability of radiation cosmologies

The conformal Einstein equations for a trace-free (radiation) perfect fluid are derived in terms of the Levi-Civita connection of a conformally rescaled metric. These equations are used to provide a non-linear stability result for de Sitter-like trace-free (radiation) perfect fluid Friedman-Lemaitre-Robertson-Walker cosmological models. The solutions thus obtained exist globally towards the future and are future geodesically complete. - Highlights: Black-Right-Pointing-Pointer We study the Einstein-Euler system in General Relativity using conformal methods. Black-Right-Pointing-Pointer We analyze the structural properties of the associated evolution equations. Black-Right-Pointing-Pointer We establish the non-linear stability of pure radiation cosmological models.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-Scale Inference of Multivariate Regression for Heavy-Tailed and Asymmetric Data

Large-scale multivariate regression is a fundamental statistical tool with a wide range of applications. Here, this study considers the problem of simultaneously testing a large number of general linear hypotheses, encompassing covariate-effect analysis, analysis of variance, and model comparisons. The challenge that accompanies a large number of tests is the ubiquitous presence of heavy-tailed and/or highly skewed measurement noise, which is the main reason for the failure of conventional least squares-based methods. For large-scale multivariate regression, we develop a set of robust inference methods to explore data features such as heavy tailedness and skewness, which are not visible to least squares methods. The new testing procedure is based on the data-adaptive Huber regression and a new covariance estimator of regression estimates. Under mild conditions, we show that our methods produce consistent estimates of the false discovery proportion. Extensive numerical experiments and an empirical study on quantitative linguistics demonstrate the advantage of the proposed method over many state-of-the-art methods when the data are generated from heavy-tailed and/or skewed distributions.

97 MATHEMATICS AND COMPUTING↗

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

Quasi-Stochastic Approximation and Off-Policy Reinforcement Learning: Preprint

The Robbins-Monro stochastic approximation algorithm is a foundation of many algorithmic frameworks for reinforcement learning and often an efficient approach to solving (or approximating the solution to) complex optimal control problems. In many cases, however, practitioners are unable to apply these techniques because of an inherent high variance. This paper aims to provide a general foundation for 'quasi-stochastic approximation,' in which all of the processes under consideration are deterministic, much like quasi-Monte Carlo for variance reduction in simulation. The variance reduction can be substantial, subject to tuning of pertinent parameters in the algorithm. This paper introduces a new coupling argument to establish the optimal rate of convergence provided the gain is sufficiently large. These results are established for linear models and tested also in nonideal settings. A major application of these general results is a new class of reinforcement learning algorithms for deterministic state space models. In this setting, the main contribution is a class of algorithms for approximating the value function for a given policy, using a different policy designed to introduce exploration.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Decoherence of spin-deformed bosonic model

The decoherence rate and some parameters affecting it are investigated for the generalized spin-boson model. We consider the spin-bosonic model when the bosonic environment is modeled by the deformed harmonic oscillators. We show that the state of the environment approaches a non-linear coherent state. Then, we obtain the decoherence rate of a two-level system which is in contact with a deformed bosonic environment which is either in thermal equilibrium or in the ground state. By using some recent realization of f-deformed oscillators, we show that some physical parameters strongly affect the decoherence rate of a two-level system. -- Highlights: •Decoherence of the generalized spin-boson model is considered. •In this model the environment consists of f-oscillators. •Via the interaction, the state of the environment approaches non-linear coherent states. •Effective parameters on decoherence are considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Photoelectron angular distributions for states of any mixed character: An experiment-friendly model for atomic, molecular, and cluster anions

We present a model for laboratory-frame photoelectron angular distributions in direct photodetachment from (in principle) any molecular orbital using linearly polarized light. A transparent mathematical approach is used to generalize the Cooper-Zare central-potential model to anionic states of any mixed character. In the limit of atomic-anion photodetachment, the model reproduces the Cooper-Zare formula. In the case of an initial orbital described as a superposition of s and p-type functions, the model yields the previously obtained s-p mixing formula. The formalism is further advanced using the Hanstorp approximation, whereas the relative scaling of the partial-wave cross-sections is assumed to follow the Wigner threshold law. The resulting model describes the energy dependence of photoelectron anisotropy for any atomic, molecular, or cluster anions, usually without requiring a direct calculation of the transition dipole matrix elements. As a benchmark case, we apply the p-d variant of the model to the experimental results for NO{sup −} photodetachment and show that the observed anisotropy trend is described well using physically meaningful values of the model parameters. Overall, the presented formalism delivers insight into the photodetachment process and affords a new quantitative strategy for analyzing the photoelectron angular distributions and characterizing mixed-character molecular orbitals using photoelectron imaging spectroscopy of negative ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solutions of nonlocal cosmological equations

A general class of gravitational models driven by a nonlocal scalar field with a linear or quadratic potential is considered. The Ostrogradski representation for nonlocal gravitational models and the way of its localization are proposed. We study the action with an arbitrary analytic function F({open_square}{sub g}), which has both simple and double roots. The way of localization allows to find solutions of nonlocal Einstein equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparisons of characteristic timescales and approximate models for Brownian magnetic nanoparticle rotations

Magnetic nanoparticles are promising tools for a host of therapeutic and diagnostic medical applications. The dynamics of rotating magnetic nanoparticles in applied magnetic fields depend strongly on the type and strength of the field applied. There are two possible rotation mechanisms and the decision for the dominant mechanism is often made by comparing the equilibrium relaxation times. This is a problem when particles are driven with high-amplitude fields because they are not necessarily at equilibrium at all. Instead, it is more appropriate to consider the “characteristic timescales” that arise in various applied fields. Approximate forms for the characteristic time of Brownian particle rotations do exist and we show agreement between several analytical and phenomenological-fit models to simulated data from a stochastic Langevin equation approach. We also compare several approximate models with solutions of the Fokker-Planck equation to determine their range of validity for general fields and relaxation times. The effective field model is an excellent approximation, while the linear response solution is only useful for very low fields and frequencies for realistic Brownian particle rotations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

D-Factor: A Quantitative Model of Application Slow-Down in Multi-Resource Shared Systems

Scheduling multiple jobs onto a platform enhances system utilization by sharing resources. The benefits from higher resource utilization include reduced cost to construct, operate, and maintain a system, which often include energy consumption. Maximizing these benefits comes at a price - resource contention among jobs increases job completion time. In this paper, we analyze slow-downs of jobs due to contention for multiple resources in a system; referred to as dilation factor. We observe that multiple-resource contention creates non-linear dilation factors of jobs. From this observation, we establish a general quantitative model for dilation factors of jobs in multi-resource systems. A job is characterized by a vector-valued loading statistics and dilation factors of a job set are given by a quadratic function of their loading vectors. We demonstrate how to systematically characterize a job, maintain the data structure to calculate the dilation factor (loading matrix), and calculate the dilation factor of each job. We validate the accuracy of the model with multiple processes running on a native Linux server, virtualized servers, and with multiple MapReduce workloads co-scheduled in a cluster. Evaluation with measured data shows that the D-factor model has an error margin of less than 16%. We also show that the model can be integrated with an existing on-line scheduler to minimize the makespan of workloads.

99 GENERAL AND MISCELLANEOUS↗