Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data fitting”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A HYBRID MODE MODEL OF THE BLAZHKO EFFECT, SHOWN TO ACCURATELY FIT KEPLER DATA FOR RR Lyr

The waveform for Blazhko stars can be substantially different during the ascending and descending parts of the Blazhko cycle. A hybrid model, consisting of two component oscillators of the same frequency, is proposed as a means to fit the data over the entire cycle. One component exhibits a sawtooth-like velocity waveform while the other is nearly sinusoidal. One method of generating such a hybrid is presented: a nonlinear model is developed for the first overtone mode, which, if excited to large amplitude, is found to drop strongly in frequency and become highly non-sinusoidal. If the frequency drops sufficiently to become equal to the fundamental frequency, the two can become phase locked and form the desired hybrid. A relationship is assumed between the hybrid mode velocity and the observed light curve, which is approximated as a power series. An accurate fit of the hybrid model is made to actual Kepler data for RR Lyr. The sinusoidal component may tend to stabilize the period of the hybrid which is found in real Blazhko data to be extremely stable. It is proposed that the variations in amplitude and phase might result from a nonlinear interaction with a third mode, possibly a nonradial mode at 3/2 the fundamental frequency. The hybrid model also applies to non-Blazhko RRab stars and provides an explanation for the light curve bump. A method to estimate the surface gravity is also proposed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

MO-FG-303-02: BEST IN PHYSICS (THERAPY): Cherenkov Emission Dosimetry: Feasibility for Electron Radiotherapy

Purpose: To investigate from first principles, corroborated by Monte Carlo simulations and experimental measurements, the feasibility of developing a relative Cherenkov emission (CE) dosimetry protocol for electron beam radiotherapy. Methods: Monte Carlo (MC) simulations of mono-energetic electrons incident on water were carried out in Geant4. Percent depth Cherenkov emission (PDCE) and dose (PDD) distributions were scored for incidence energies of 4, 6, 9, 12, 15, and 18 MeV. PDCE-to-PDD analytical conversion models were developed from least-squares data fits generated for PDD as a function of PDCE at the same depth and at different depths. Experimental techniques for validation of these models are examined. Results: Same-depth PDD versus PDCE data fits indicate that although the relationship is linear to first order (correlation r > 0.9 for all energies), it is much more accurately approximated by separate linear and quadratic models for the build-up and drop-off regions, respectively (r > 0.999), which is theoretically underpinned. To understand the source of this relationship and its basis for developing robust conversion models, an approximate quadratic first-principles model was derived and found in agreement with MC/measured data (20% deviation at worst). Conversely, data fits of PDD versus different-depth PDCE unveiled a depth-invariant effective point of measurement of 1.5–2.1 mm downstream with 4–18 MeV incidence, respectively (r > 0.999 in the drop-off region). We present an analytical first-principles justification for this shift. This method led to errors of <1% in drop-off region PDD (<2% for PDD<20% with 4 MeV incidence) and <0.2 mm in practical range prediction. Conclusion: We present robust quantitative prediction models, derived from first-principles and supported by simulation and measurement, for relative dose from Cherenkov emission by high-energy electrons. This constitutes a major step towards development of protocols for routine clinical quality assurance as well as real-time in vivo Cherenkov dosimetry in radiotherapy. The authors acknowledge partial support by Fonds de recherche du Quebec - Nature et technologies (FRQNT), CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council of Canada (NSERC), CREATE Integrated Sensor Systems grant of NSERC, the Canadian Institutes of Health Research (CIHR), and NSERC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-685: Evaluation of Tumor Hypoxic Fraction Using Serial Volumetric Imaging During Radiation Therapy

Purpose: To develop a tumor response model which could be uses to compute tumor hypoxic fraction using serial volumetric tumor imaging. This algorithm may be used for treatment response assessment and also for guidance of more expensive PET imaging of hypoxia. Methods: Previously developed two-level cell population tumor response model was modified to include a third cell level describing hypoxic and necrotic cells. This third level was considered constant value during radiotherapy treatment; therefore, inclusion additional parameter did not compromise stability of model fitting to imaging data. Fitting the model to serial volumetric imaging data was performed using a least squares objective function and simulated annealing algorithm. The problem of reconstruction of radiobiological parameters from serial imaging data was considered as inverse ill-posed problem described by the Fredholm integral equation of the first kind. Variational regularization was used to stabilize solutions. Results: To evaluate performance of the algorithm, we used a set of serial CT imaging data on tumor-volume for 14 head and neck cancer patients. The hypoxic fractions were reconstructed for each patient and the distribution of hypoxic fractions was compared to the distribution of initial hypoxic fractions previously measured using histograph. The measured and reconstructed from imaging data distributions of hypoxic fractions are in good agreement. The reconstructed distribution of cell surviving fraction was also in better agreement with in vitro data than previously obtained using the two-level cell population model. Conclusion: Our results indicate that it is possible to evaluate the initial hypoxic tumor fraction using serial volumetric imaging and a tumor response model. This algorithm can be used for treatment response assessment and guidance of more expensive PET imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

A Parallelised ROOT for Future HEP Data Processing

In the coming years, HEP data processing will need to exploit parallelism on present and future hardware resources to sustain the bandwidth requirements. As one of the cornerstones of the HEP software ecosystem, ROOT embraced an ambitious parallelisation plan which delivered compelling results. In this contribution the strategy is characterised as well as its evolution in the medium term. The units of the ROOT framework are discussed where task and data parallelism have been introduced, with runtime and scaling measurements. We will give an overview of concurrent operations in ROOT, for instance in the areas of I/O (reading and writing of data), fitting / minimization, and data analysis. This paper introduces the programming model and use cases for explicit and implicit parallelism, where the former is explicit in user code and the latter is implicitly managed by ROOT internally.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Dynamics of neutrals and ions in an ultrafast laser produced Zn plasma

Optical time of flight dynamics of neutrals and ions in an ultrafast laser produced zinc plasma generated by irradiating a solid zinc target using 100 fs laser pulses is investigated. An acceleration of ions is observed which arises from internal Coulomb forces acting between charged species in the plasma. Some of the fast ions recombine with electrons in the plasma and generate fast neutrals. Plasma plume imaging performed at various ambient pressures indicates adiabatic expansion at lower pressures and plume front deceleration at higher pressures: at lower pressures the plume front-time (R-t) plot displays a linear expansion, shock wave model fits to the data at 5 Torr and at higher pressures the data fits better to the drag model. Furthermore, around an intermediate pressure of 10 Torr, the R-t plot fits to the shock wave model at earlier stages of plasma expansion, while it fits to the drag model at the later stages. These investigations provide relevant information on the acceleration of ions and neutrals in an expanding zinc plasma plume produced by ultrafast laser pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Bottom-up design of actinide materials from molecular clusters: Demonstration of a general-purpose simulation capability leveraging machine-learned atomic potentials

Actinide thin-film coatings such as uranium dioxide (UO 2 ) play an important role in nuclear reactors and other mission-relevant applications, but realization of their potential requires a deep fundamental understanding of the chemical vapor deposition (CVD) processes used for their growth. The slow experimental progress can be attributed, in part, to the standard safety guidelines associated with handling uranium byproducts, which are often corrosive, toxic, and radioactive. Accurate simulation techniques, when used in concert with experiment, can improve laboratory safety, material durability, and deliverable timeframes. However, state-of-the-art computational methods are either insufficiently accurate or intractably expensive. To remedy this situation, in this project we suggested a machine-learning (ML) accelerated workflow for simulating molecular clustering toward deposition. As a benchmark test case, we considered molecular clustering in steam and assessed independent components of our workflow by comparing with measured thermodynamic properties of water. After analyzing each component individually and finding no fundamental barrier to realization of the workflow, we attempted to integrate the ML component, a Sandia-developed tool called FitSNAP. As this was the first application of FitSNAP to atoms and molecules in the gas phase at Sandia, the method required more fitting data than was originally anticipated. Systematic improvements were made by including in the fit data diatomic potentials, molecular single-bond-breaking curves, and symmetry-constrained intermolecular potentials. We concluded that our strategy provides a feasible pathway toward modeling CVD and related processes, but that extensive training data must be generated before it can be of practical use.

36 MATERIALS SCIENCE↗
DOE OSTI2013

The eigenmode perspective of NMR spin relaxation in proteins

We developed in recent years the two-body (protein and probe) coupled-rotator slowly relaxing local structure (SRLS) approach for elucidating protein dynamics from NMR spin relaxation. So far we used as descriptors the set of physical parameters that enter the SRLS model. They include the global (protein-related) diffusion tensor, D{sub 1}, the local (probe-related) diffusion tensor, D{sub 2}, and the local coupling/ordering potential, u. As common in analyzes based on mesoscopic dynamic models, these parameters have been determined with data-fitting techniques. In this study, we describe structural dynamics in terms of the eigenmodes comprising the SRLS time correlation functions (TCFs) generated by using the best-fit parameters as input to the Smoluchowski equation. An eigenmode is a weighted exponential with decay constant given by an eigenvalue of the Smoluchowski operator, and weighting factor determined by the corresponding eigenvector. Obviously, both quantities depend on the SRLS parameters as determined by the SRLS model. Unlike the set of best-fit parameters, the eigenmodes represent patterns of motion of the probe-protein system. The following new information is obtained for the typical probe, the {sup 15}N−{sup 1}H bond. Two eigenmodes, associated with the protein and the probe, dominate when the time scale separation is large (i.e., D{sub 2} ≫ D{sub 1}), the tensorial properties are simple, and the local potential is either very strong or very weak. When the potential exceeds these limits while the remaining conditions are preserved, new eigenmodes arise. The multi-exponentiality of the TCFs is associated in this case with the restricted nature of the local motion. When the time scale separation is no longer large, the rotational degrees of freedom of the protein and the probe become statistically dependent (coupled dynamically). The multi-exponentiality of the TCFs is associated in this case with the restricted nature of both the local and the global motion. The effects of local diffusion axiality, potential strength, and extent of mode-coupling on the eigenmode setup are investigated. We detect largely global motional or largely local motional eigenmodes. In addition, we detect mixed eigenmodes associated with correlated/prograde or anti-correlated/retrograde rotations of the global (D{sub 1}) and local (D{sub 2}) motional modes. The eigenmode paradigm is applied to N−H bond dynamics in the β-sheet residue K19, and the α-helix residue A34, of the third immunoglobulin-binding domain of streptococcal protein G. The largest contribution to the SRLS TCFs is made by mixed anti-correlated D{sub 1} and D{sub 2} eigenmodes. The next largest contribution is made by D{sub 1}-dominated eigenmodes. Eigenmodes dominated by the local motion contribute appreciably to A34 and marginally to K19. Correlated D{sub 1} and D{sub 2} eigenmodes contribute exclusively to K19 and do not contribute above 1% to A34. The differences between K19 and A34 are delineated and rationalized in terms of the best-fit SRLS parameters and mode-mixing. It may be concluded that eigenmode analysis is complementary and supplementary to data-fitting-based analysis.

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

A semblance measure for model comparison

Algorithmic and computational advances have made it possible that geophysical survey and earth model design can be aided by many systematic trial inverse-modelling runs with synthetic data. Such may, for example, come up in machine-learning approaches. Automated image appraisal pertaining to such applications will involve common statistical tests for goodness-of-data fit as a primary evaluation method. However, solution non-uniqueness may render multiple images equivalent in terms of their data fit, requiring secondary categorizers. A logical choice for classifying synthetic-imaging results quantifies the goodness of model fit where a known reference model replaces the observational input. The task of model intercomparison in terms of measuring the resemblance to the reference model poses challenges to common distance-based metrics like root mean square error and mean absolute error. First, distance-based metrics can introduce spurious contributions when smooth models with fuzzy target contours are to be compared against a sharp reference. Second, large differences due to parameter-estimation overshoots can dominate distance metrics. Here, we propose a remedy that is referred to as semblance and is based on the idea of logistic functions, where a binary-dependent variable adds non-zero or zero accumulation terms for the, respectively, passing or failing of preset target thresholds. This classifying approach is amenable to an objective where model feature recognition is primary. Numerical comparisons to distance-based metrics provide evidence for the advantages of the semblance in view of this objective. Geophysical imaging in conjunction with machine-learning is seen as a benefitting upcoming application area.

58 GEOSCIENCES↗
DOE OSTI2021

Multi-Omics Driven Metabolic Network Reconstruction and Analysis of Lignocellulosic Carbon Utilization in Rhodosporidium toruloides

An oleaginous yeast Rhodosporidium toruloides is a promising host for converting lignocellulosic biomass to bioproducts and biofuels. In this work, we performed multi-omics analysis of lignocellulosic carbon utilization in R. toruloides and reconstructed the genome-scale metabolic network of R. toruloides . High-quality metabolic network models for model organisms and orthologous protein mapping were used to build a draft metabolic network reconstruction. The reconstruction was manually curated to build a metabolic model using functional annotation and multi-omics data including transcriptomics, proteomics, metabolomics, and RB-TDNA sequencing. The multi-omics data and metabolic model were used to investigate R. toruloides metabolism including lipid accumulation and lignocellulosic carbon utilization. The developed metabolic model was validated against high-throughput growth phenotyping and gene fitness data, and further refined to resolve the inconsistencies between prediction and data. We believe that this is the most complete and accurate metabolic network model available for R. toruloides to date.

09 BIOMASS FUELS↗
DOE OSTI2016

Bayesian evidence and predictivity of the inflationary paradigm

In this paper we consider the issue of paradigm evaluation by applying Bayes' theorem along the following nested hierarchy of progressively more complex structures: i) parameter estimation (within a model), ii) model selection and comparison (within a paradigm), iii) paradigm evaluation. In such a hierarchy the Bayesian evidence works both as the posterior's normalization at a given level and as the likelihood function at the next level up. Whilst raising no objections to the standard application of the procedure at the two lowest levels, we argue that it should receive a considerable modification when evaluating paradigms, when testability and fitting data are equally important. By considering toy models we illustrate how models and paradigms that are difficult to falsify are always favoured by the Bayes factor. We argue that the evidence for a paradigm should not only be high for a given dataset, but exceptional with respect to what it would have been, had the data been different. With this motivation we propose a measure which we term predictivity , as well as a prior to be incorporated into the Bayesian framework, penalising unpredictivity as much as not fitting data. We apply this measure to inflation seen as a whole, and to a scenario where a specific inflationary model is hypothetically deemed as the only one viable as a result of information alien to cosmology (e.g. Solar System gravity experiments, or particle physics input). We conclude that cosmic inflation is currently hard to falsify, but that this could change were external/additional information to cosmology to select one of its many models. We also compare this state of affairs to bimetric varying speed of light cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fatty Acid and Alcohol Metabolism in Pseudomonas putida: Functional Analysis Using Random Barcode Transposon Sequencing.

With its ability to catabolize a wide variety of carbon sources and a growing engineering toolkit, Pseudomonas putida KT2440 is emerging as an important chassis organism for metabolic engineering. Despite advances in our understanding of the organism, many gaps remain in our knowledge of the genetic basis of its metabolic capabilities. The gaps are particularly noticeable in our understanding of both fatty acid and alcohol catabolism, where many paralogs putatively coding for similar enzymes coexist, making biochemical assignment via sequence homology difficult. To rapidly assign function to the enzymes responsible for these metabolisms, we leveraged random barcode transposon sequencing (RB-Tn-Seq). Global fitness analyses of transposon libraries grown on 13 fatty acids and 10 alcohols produced strong phenotypes for hundreds of genes. Fitness data from mutant pools grown on fatty acids of varying chain lengths indicated specific enzyme substrate preferences and enabled us to hypothesize that DUF1302/DUF1329 family proteins potentially function as esterases. From the data, we also postulate catabolic routes for the two biogasoline molecules isoprenol and isopentanol, which are catabolized via leucine metabolism after initial oxidation and activation with coenzyme A (CoA). Because fatty acids and alcohols may serve as both feedstocks and final products of metabolic-engineering efforts, the fitness data presented here will help guide future genomic modifications toward higher titers, rates, and yields. IMPORTANCE To engineer novel metabolic pathways into P. putida , a comprehensive understanding of the genetic basis of its versatile metabolism is essential. Here, we provide functional evidence for the putative roles of hundreds of genes involved in the fatty acid and alcohol metabolism of the bacterium. These data provide a framework facilitating precise genetic changes to prevent product degradation and to channel the flux of specific pathway intermediates as desired.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Polarimetric Geometric Modeling for mm-VLBI Observations of Black Holes

The Event Horizon Telescope (EHT) is a millimeter very long baseline interferometry (VLBI) array that has imaged the apparent shadows of the supermassive black holes M87* and Sagittarius A*. Polarimetric data from these observations contain a wealth of information on the black hole and accretion flow properties. In this work, we develop polarimetric geometric modeling methods for mm-VLBI data, focusing on approaches that fit data products with differing degrees of invariance to broad classes of calibration errors. We establish a fitting procedure using a polarimetric “m-ring” model to approximate the image structure near a black hole. By fitting this model to synthetic EHT data from general relativistic magnetohydrodynamic models, we show that the linear and circular polarization structure can be successfully approximated with relatively few model parameters. We then fit this model to EHT observations of M87* taken in 2017. In total intensity and linear polarization, the m-ring fits are consistent with previous results from imaging methods. In circular polarization, the m-ring fits indicate the presence of event-horizon-scale circular polarization structure, with a persistent dipolar asymmetry and orientation across several days. The same structure was recovered independently of observing band, used data products, and model assumptions. Despite this broad agreement, imaging methods do not produce similarly consistent results. Our circular polarization results, which imposed additional assumptions on the source structure, should thus be interpreted with some caution. Polarimetric geometric modeling provides a useful and powerful method to constrain the properties of horizon-scale polarized emission, particularly for sparse arrays like the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Detection of the pairwise kinematic Sunyaev-Zel'dovich effect with BOSS DR11 and the Atacama Cosmology Telescope

Here, we present a new measurement of the kinematic Sunyaev-Zel'dovich effect using data from the Atacama Cosmology Telescope (ACT) and the Baryon Oscillation Spectroscopic Survey (BOSS). Using 600 square degrees of overlapping sky area, we evaluate the mean pairwise baryon momentum associated with the positions of 50,000 bright galaxies in the BOSS DR11 Large Scale Structure catalog. A non-zero signal arises from the large-scale motions of halos containing the sample galaxies. The data fits an analytical signal model well, with the optical depth to microwave photon scattering as a free parameter determining the overall signal amplitude. We estimate the covariance matrix of the mean pairwise momentum as a function of galaxy separation, using microwave sky simulations, jackknife evaluation, and bootstrap estimates. The most conservative simulation-based errors give signal-to-noise estimates between 3.6 and 4.1 for varying galaxy luminosity cuts. We discuss how the other error determinations can lead to higher signal-to-noise values, and consider the impact of several possible systematic errors. Estimates of the optical depth from the average thermal Sunyaev-Zel'dovich signal at the sample galaxy positions are broadly consistent with those obtained from the mean pairwise momentum signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Projected Multi-Agent Consensus Equilibrium (PMACE) With Application to Ptychography

Multi-Agent Consensus Equilibrium (MACE) formulates an inverse imaging problem as a balance among multiple update agents such as data-fitting terms and denoisers. However, each such agent operates on a separate copy of the full image, leading to redundant memory use and slow convergence when each agent affects only a small subset of the full image. In this article, we extend MACE to Projected Multi-Agent Consensus Equilibrium (PMACE), in which each agent updates only a projected component of the full image, thus greatly reducing memory use for some applications. We describe PMACE in terms of an equilibrium problem and an equivalent fixed-point problem and show that in most cases the PMACE equilibrium is not the solution of an optimization problem. To demonstrate the value of PMACE, we apply it to the problem of ptychography, in which a sample is reconstructed from the diffraction patterns resulting from coherent X-ray illumination at multiple overlapping spots. In our PMACE formulation, each spot corresponds to a separate data-fitting agent, with the final solution found as an equilibrium among all the agents. In conclusion, our results demonstrate that the PMACE reconstruction algorithm generates more accurate reconstructions at a lower computational cost than existing ptychography algorithms when the spots are sparsely sampled.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Iron oxide hydroxide nanoflower assisted removal of arsenic from water

Graphical abstract: Non-magnetic polycrystalline iron oxide hydroxide nanoparticle with flower like morphology is found to play as an effective adsorbent media to remove As(III) from 300 μg L{sup −1} to less than 10 μg L{sup −1} from drinking water over wide range of pH. TEM image clearly reveals that the nanoparticle looks flower like morphology with average particle size less than 20 nm. The maximum sorption capacity of the sorbent is found to be 475 μg g{sup −1} for arsenic at room temperature and the data fitted to different isotherm models indicate the heterogeneity of the adsorbent surface. The material can be regenerated up to 70% using dilute hydrochloric acid and it would be utilized for de-arsenification purposes. - Highlights: • The work includes synthesis of iron oxide hydroxide nanoflower and its applicability for the removal of arsenic from water. • The nanoparticle was characterized using modern instrumental methods like FESEM, TEM, BET, XRD, etc. • The maximum sorption capacity of the sorbent is found to be 475 μg g{sup −1} for arsenic at room temperature. • The sorption is multilayered on the heterogeneous surface of the nano adsorbent. • The mechanism of arsenic removal of IOH nanoflower follows both adsorption and ion-exchange. - Abstract: Non-magnetic polycrystalline iron oxide hydroxide nanoparticle with flower like morphology is found to play as an effective adsorbent media to remove As(III) from 300 μg L{sup −1} to less than 10 μg L{sup −1} from drinking water over wide range of pH. The nanoparticle was characterized by X-ray powder diffraction analysis (XRD), BET surface area, FTIR, FESEM and TEM images. TEM image clearly reveals flower like morphology with average particle size less than 20 nm. The nanoflower morphology is also supported by FESEM images. The maximum sorption capacity of the sorbent is found to be 475 μg g{sup −1} for arsenic and the data fitted to different isotherm models indicate the heterogeneity of the adsorbent surface. Study on adsorption kinetics shows that adsorption of arsenic onto iron oxide hydroxide nanoflower follows pseudo-second order kinetic. The material can be regenerated up to 70% using dilute hydrochloric acid and it would be utilized for de-arsenification purposes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lithium-ion battery physics and statistics-based state of health model

A pseudo-2d model using COMSOL Multiphysics® software is developed to simulate performance and performance degradation of Li-ion batteries consisting of layered and olivine cathodes with graphite anode when subjected to peak shaving grid service. Multiple degradation pathways are considered, including solid electrolyte interphase (SEI) formation and breakdown at the anode, cathode dissolution and its synergistic effect on SEI formation at the anode. The model is validated by simulating commercial cylindrical cell performance. A global model is developed to simulate performance across all chemistries, along with individual chemistry models using global model parameters as initial values. There is good agreement between these models for various optimization parameters such as SEI equilibrium potential, cathode dissolution exchange current density, solvent diffusivity in the SEI and SEI ionic conductivity. To circumvent time constraints related to the COMSOL model, a 0d global model is developed which fits data well and provides more clarity on differences in cathode dissolution exchange current density. Again, good agreement for various optimization parameters is obtained among the COMSOL global & individual chemistry models and the 0-d model. The lessons learned from the physics-based model is used to develop a top down statistics-based model using current, voltage and anode volumetric change per mole lithium intercalated, along with their interactions as degradation predictors. This model predicts out of sample degradation for multiple grid services and electric vehicle drive cycle with high accuracy and provides the pathway to develop an efficient battery management system combining machine learning and findings from physics-based computationally intensive algorithms.

Crawford, Aladsair J.↗
DOE OSTI2015

New axion and hidden photon constraints from a solar data global fit

We present a new statistical analysis that combines helioseismology (sound speed, surface helium and convective radius) and solar neutrino observations (the {sup 8}B and {sup 7}Be fluxes) to place upper limits to the properties of non standard weakly interacting particles. Our analysis includes theoretical and observational errors, accounts for tensions between input parameters of solar models and can be easily extended to include other observational constraints. We present two applications to test the method: the well studied case of axions and axion-like particles and the more novel case of low mass hidden photons. For axions we obtain an upper limit at 3σ for the axion-photon coupling constant of g{sub aγ} < 4.1 · 10{sup −10} GeV{sup −1}. For hidden photons we obtain the most restrictive upper limit available accross a wide range of masses for the product of the kinetic mixing and mass of χ m < 1.8 ⋅ 10{sup −12} eV at 3σ. Both cases improve the previous solar constraints based on the Standard Solar Models showing the power of using a global statistical approach.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling the partitioning of amphiphilic molecules and co-solvents in biomembranes

We report amphiphilic co-solvents can have a significant impact on the structure, organization and physical properties of lipid bilayers. Describing the mutual impact of partitioning and induced structure changes is therefore a crucial consideration for a range of topics such as anesthesia and other pharmacokinetic effects, as well as microbial solvent tolerance in the production of biofuels and other fermentation products, where molecules such as ethanol, butanol or acetic acid might be generated. Small-angle neutron scattering (SANS) is a key method for studying lipid and polymer bilayer structures, with many models for extracting bilayer structure (thickness, area per lipid etc.) from scattering data in use today. However, the molecular details of co-solvent partitioning are conflated with induced changes to bilayer structure, making interpretation and modeling of the scattering curves a challenge with the existing set of models. To address this, a model of a bilayer structure is presented which invokes a two-term partition constant accounting for the localization of the co-solvent within the bilayer. This model was validated using a series of SANS measurements of lipid vesicles in the presence of the co-solvent tetrahydrofuran (THF), showing several strategies of how to deploy the two-parameter partition constant model to describe scattering data and extract both structure and partitioning information from the data. Molecular dynamics simulations are then used to evaluate assumptions of the model, provide additional molecular scale details and illustrate its complementary nature to the data fitting procedure. This approach results in estimates of the partition coefficient for THF in 1,2-dimyristoyl-sn-glycero-3-phosphocholine at 35°C, along with an estimate of the fraction of THF residing in the hydrophobic core of the membrane. The authors envision that this model will be applicable to a wide range of other bilayer/amphiphile interactions and provide the associated code needed to implement this model as a fitting algorithm for scattering data in the SasView suite.

59 BASIC BIOLOGICAL SCIENCES↗