Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Universality Classes of Interaction Structures for NK Fitness Landscapes

Kauffman’s NK-model is a paradigmatic example of a class of stochastic models of genotypic fitness landscapes that aim to capture generic features of epistatic interactions in multilocus systems. Genotypes are represented as sequences of L binary loci. The fitness assigned to a genotype is a sum of contributions, each of which is a random function defined on a subset of k≤L loci. These subsets or neighborhoods determine the genetic interactions of the model. Whereas earlier work on the NK model suggested that most of its properties are robust with regard to the choice of neighborhoods, recent work has revealed an important and sometimes counter-intuitive influence of the interaction structure on the properties of NK fitness landscapes. Here we review these developments and present new results concerning the number of local fitness maxima and the statistics of selectively accessible (that is, fitness-monotonic) mutational pathways. In particular, we develop a unified framework for computing the exponential growth rate of the expected number of local fitness maxima as a function of L, and identify two different universality classes of interaction structures that display different asymptotics of this quantity for large k. Moreover, we show that the probability that the fitness landscape can be traversed along an accessible path decreases exponentially in L for a large class of interaction structures that we characterize as locally bounded. Finally, we discuss the impact of the NK interaction structures on the dynamics of evolution using adaptive walk models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Making Steppingstones out of Stumbling Blocks: A Bayesian Model Evidence Estimator with Application to Groundwater Transport Model Selection

Bayesian model evidence (BME) is a measure of the average fit of a model to observation data given all the parameter values that the model can assume. By accounting for the trade-off between goodness-of-fit and model complexity, BME is used for model selection and model averaging purposes. For strict Bayesian computation, the theoretically unbiased Monte Carlo based numerical estimators are preferred over semi-analytical solutions. This study examines five BME numerical estimators and asks how accurate estimation of the BME is important for penalizing model complexity. The limiting cases for numerical BME estimators are the prior sampling arithmetic mean estimator (AM) and the posterior sampling harmonic mean (HM) estimator, which are straightforward to implement, yet they result in underestimation and overestimation, respectively. We also consider the path sampling methods of thermodynamic integration (TI) and steppingstone sampling (SS) that sample multiple intermediate distributions that link the prior and the posterior. Although TI and SS are theoretically unbiased estimators, they could have a bias in practice arising from numerical implementation. For example, sampling errors of some intermediate distributions can introduce bias. We propose a variant of SS, namely the multiple one-steppingstone sampling (MOSS) that is less sensitive to sampling errors. We evaluate these five estimators using a groundwater transport model selection problem. SS and MOSS give the least biased BME estimation at an efficient computational cost. If the estimated BME has a bias that covariates with the true BME, this would not be a problem because we are interested in BME ratios and not their absolute values. On the contrary, the results show that BME estimation bias can be a function of model complexity. Thus, biased BME estimation results in inaccurate penalization of more complex models, which changes the model ranking. This was less observed with SS and MOSS as with the three other methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

HZ Her: Stellar radius from X-ray eclipse observations, evolutionary state, and a new distance

Observations of HZ Her/Her X-1 by the Rossi X-Ray Timing Explorer (RXTE) covering high state eclipses of the neutron star are analyzed here. Models of the eclipse are used to measure the radius and atmospheric scale height of HZ Her, the stellar companion to the neutron star. The radius is 2.58-3.01 × 10{sup 11} cm, depending on system inclination and mass ratio (q), with an accuracy of ∼1 part in 1000 for given inclination and q. We fit Kurucz model stellar atmosphere models to archival optical observations. The resulting effective temperature (T {sub eff}) of the unheated face of HZ Her is determined to be in the 2σ range of 7720 K-7865 K, and metallicity (log (Z/Z {sub ☉})) in the range of –0.27 to +.03. The model atmosphere surface flux and new radius yield a new distance to HZ Her/Her X-1, depending on system inclination and q: a best-fit value of 6.1 kpc with upper and lower limits of 5.7 kpc and 7.0 kpc. We calculate stellar evolution models for the range of allowed masses (from orbital parameters) and allowed metallicities (from optical spectrum fits). The stellar models agree with T {sub eff} and the radius of HZ Her for two narrow ranges of mass: 2.15-2.20 M {sub ☉} and 2.35-2.45 M {sub ☉}. This lower mass range implies a low neutron star mass (1.3 M {sub ☉}), whereas the higher mass range implies a high neutron star mass (1.5-1.7 M {sub ☉}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Forecasting climate change impacts on plant populations over large spatial extents

Plant population models are powerful tools for predicting climate change impacts in one location, but are difficult to apply at landscape scales. Here, we overcome this limitation by taking advantage of two recent advances: remotely sensed, species-specific estimates of plant cover and statistical models developed for spatiotemporal dynamics of animal populations. Using computationally efficient model reparameterizations, we fit a spatiotemporal population model to a 28-year time series of sagebrush (Artemisia spp.) percent cover over a 2.5 × 5 km landscape in southwestern Wyoming while formally accounting for spatial autocorrelation. We include interannual variation in precipitation and temperature as covariates in the model to investigate how climate affects the cover of sagebrush. We then use the model to forecast the future abundance of sagebrush at the landscape scale under projected climate change, generating spatially explicit estimates of sagebrush population trajectories that have, until now, been impossible to produce at this scale. Our broadscale and long-term predictions are rooted in small-scale and short-term population dynamics and provide an alternative to predictions offered by species distribution models that do not include population dynamics. Finally, our approach, which combines several existing techniques in a novel way, demonstrates the use of remote sensing data to model population responses to environmental change that play out at spatial scales far greater than the traditional field study plot.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Adsorption Studies of Dimetridazole and Metronidazole onto Biochar Derived from Sugarcane Bagasse: Kinetic, Equilibrium, and Mechanisms

Biochar derived from sugarcane bagasse was prepared at different pyrolysis temperatures and holding times for adsorption of dimetridazole and metronidazole from aqueous solution. The optimal pyrolysis temperature and holding time were 500 °C and 120 min, respectively. The removal efficiencies of dimetridazole and metronidazole were 98.0% and 88.8%, respectively, at biochar dosage of 4.0 g L{sup −1} and solution temperature of 30 °C. Kinetics data were fitted by four kinetics models, and results indicated that adsorption was governed by pseudo-second-order kinetics. Adsorption is a multi-step process that involves film diffusion and pore filling. In addition, five isotherm models were employed to describe the adsorption equilibrium. Brunauer–Emmett–Teller model presented better fitting for the adsorption equilibrium, and the maximum adsorption capacities were 72.17 and 23.61 mg g{sup − 1} for dimetridazole and metronidazole, respectively. The calculated values of ΔG{sup 0} and ΔH{sup 0} indicated the non-spontaneous and exothermic nature of the adsorption process at the range of temperature studied. Thermodynamic studies also revealed that physical and chemical adsorptions were co-action. The adsorption mechanisms of dimetridazole and metronidazole onto biochar were mainly hydrogen bonding and π–π interaction. All the results revealed that sugarcane bagasse biochar can be used as alternative to costly adsorbents for the removal of dimetridazole and metronidazole from aqueous solution.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Kinetic model for quartz and spinel dissolution during melting of high-level-waste glass batch

The dissolution of quartz particles and the growth and dissolution of crystalline phases during the conversion of batch to glass potentially affects both the glass melting process and product quality. Crystals of spinel exiting the cold cap to molten glass below can be troublesome during the vitrification of iron-containing high-level wastes. To estimate the distribution of quartz and spinel fractions within the cold cap, we used kinetic models that relate fractions of these phases to temperature and heating rate. Fitting the model equations to data showed that the heating rate, apart from affecting quartz and spinel behavior directly, also affects them indirectly via concurrent processes, such as the formation and motion of bubbles. Because of these indirect effects, it was necessary to allow one kinetic parameter (the pre-exponential factor) to vary with the heating rate. The resulting kinetic equations are sufficiently simple for the detailed modeling of batch-to-glass conversion as it occurs in glass melters. The estimated fractions and sizes of quartz and spinel particles as they leave the cold cap, determined in this study, will provide the source terms needed for modeling the behavior of these solid particles within the flow of molten glass in the melter.

Batch melting, Cold cap, Quartz, Spinel, Waste gla↗
DOE OSTI2022

The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can therefore help to understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on the following galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies (log M*/M ⊙ ≳ 11.4) from the Sloan Digital Sky Survey spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sérsic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame r-band luminosity, half-light radius, and rest-frame u - g colour, using KiDS weak lensing measurements and a Bayesian hierarchical approach. For the sake of robustness with respect to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing the total luminosity and half-light radius with the luminosity within a 10 kpc aperture, L r, 10 , and the light-weighted surface brightness slope, Γ 10 . Results. We did not detect any correlation between the halo mass and either the half-light radius or colour at fixed redshift and luminosity. Using the robust surface brightness parameterisation, we found that the halo mass correlates weakly with L r,10 and anti-correlates with Γ 10 . At fixed redshift, L r, 10 and Γ 10 , the difference in the average halo mass between galaxies at the 84th percentile and 16th percentile of the colour distribution is 0.00 ± 0.11 dex. Conclusion. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than the halo mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

High-frequency characterization of thermionic charge transport in silicon-on-insulator nanowire transistors

We report on DC and microwave electrical transport measurements in silicon-on-insulator nano-transistors at low and room temperature. At low source-drain voltage, the DC current and radio frequency response show signs of conductance quantization. We attribute this to Coulomb blockade resulting from barriers formed at the spacer-gate interfaces. We show that at high bias transport occurs thermionically over the highest barrier: Transconductance traces obtained from microwave scattering-parameter measurements at liquid helium and room temperature are accurately fitted by a thermionic model. From the fits we deduce the ratio of gate capacitance and quantum capacitance, as well as the electron temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Radial velocity variations of photometrically quiet, chromospherically inactive Kepler stars: A link between RV jitter and photometric flicker

We compare stellar photometric variability, as measured from Kepler light curves by Basri et al., with measurements of radial velocity (RV) rms variations of all California Planet Search overlap stars. We newly derive rotation periods from the Kepler light curves for all of the stars in our study sample. The RV variations reported herein range from less than 4 to 135 m s{sup –1}, yet the stars all have amplitudes of photometric variability less than 3 mmag, reflecting the preference of the RV program for chromospherically 'quiet' stars. Despite the small size of our sample, we find with high statistical significance that the RV rms manifests strongly in the Fourier power spectrum of the light curve: stars that are noisier in RV have a greater number of frequency components in the light curve. We also find that spot models of the observed light curves systematically underpredict the observed RV variations by factors of ∼2-1000, likely because the low-level photometric variations in our sample are driven by processes not included in simple spot models. The stars best fit by these models tend to have simpler light curves, dominated by a single relatively high-amplitude component of variability. Finally, we demonstrate that the RV rms behavior of our sample can be explained in the context of the photometric variability evolutionary diagram introduced by Bastien et al. We use this diagram to derive the surface gravities of the stars in our sample, revealing many of them to have moved off the main sequence. More generally, we find that the stars with the largest RV rms are those that have evolved onto the 'flicker floor' sequence in that diagram, characterized by relatively low amplitude but highly complex photometric variations which grow as the stars evolve to become subgiants.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Global study of the simplest scalar phantom dark matter model

We present a global study of the simplest scalar phantom dark matter model. The best fit parameters of the model are determined by simultaneously imposing (i) relic density constraint from WMAP, (ii) 225 live days data from direct experiment XENON100, (iii) upper limit of gamma-ray flux from Fermi-LAT indirect detection based on dwarf spheroidal satellite galaxies, and (iv) the Higgs boson candidate with a mass about 125 GeV and its invisible branching ratio no larger than 40% if the decay of the Higgs boson into a pair of dark matter is kinematically allowed. The allowed parameter space is then used to predict annihilation cross sections for gamma-ray lines, event rates for three processes mono-b jet, single charged lepton and two charged leptons plus missing energies at the Large Hadron Collider, as well as to evaluate the muon anomalous magnetic dipole moment for the model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

ALMA MEASUREMENTS OF THE HNC AND HC{sub 3}N DISTRIBUTIONS IN TITAN'S ATMOSPHERE

We present spectrally and spatially resolved maps of HNC and HC{sub 3}N emission from Titan's atmosphere, obtained using the Atacama Large Millimeter/submillimeter Array on 2013 November 17. These maps show anisotropic spatial distributions for both molecules, with resolved emission peaks in Titan's northern and southern hemispheres. The HC{sub 3}N maps indicate enhanced concentrations of this molecule over the poles, consistent with previous studies of Titan's photochemistry and atmospheric circulation. Differences between the spectrally integrated flux distributions of HNC and HC{sub 3}N show that these species are not co-spatial. The observed spectral line shapes are consistent with HNC being concentrated predominantly in the mesosphere and above (at altitudes z ≳ 400 km), whereas HC{sub 3}N is abundant at a broader range of altitudes (z ≈ 70-600 km). From spatial variations in the HC{sub 3}N line profile, the locations of the HC{sub 3}N emission peaks are shown to be variable as a function of altitude. The peaks in the integrated emission from HNC and the line core (upper atmosphere) component of HC{sub 3}N (at z ≳ 300 km) are found to be asymmetric with respect to Titan's polar axis, indicating that the mesosphere may be more longitudinally variable than previously thought. The spatially integrated HNC and HC{sub 3}N spectra are modeled using the NEMESIS planetary atmosphere code and the resulting best-fitting disk-averaged vertical mixing ratio profiles are found to be in reasonable agreement with previous measurements for these species. Vertical column densities of the best-fitting gradient models for HNC and HC{sub 3}N are 1.9 × 10{sup 13} cm{sup –2} and 2.3 × 10{sup 14} cm{sup –2}, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The GMRT archive atomic gas survey – II. Mass modelling and dark matter halo properties across late-type spirals

Studying the kinematics and mass modelling of galaxies from H i 21 cm data provides valuable insights into the properties of both the baryonic components and the dark matter halo in nearby galaxies. Despite many observational studies, mass modelling of galaxies remains challenging due to different limitations. For example, most of the previous studies involving mass modelling are based on rotation curves derived from 2D velocity fields from H i or H α spectroscopic observation which are often affected by beam smearing and projection effect. However, kinematic modelling done by fitting the ‘Tilted ring model’ to 3D data cube is not affected by these issues. In this study, we present and compare 3D kinematic modelling of a pilot sample of 11 galaxies from the GMRT archive atomic gas survey (GARCIA) using two different publicly available pipelines. We model the observed H i rotation curve using 3.6-μm infrared data and SDSS r-band data for stellar contribution, H i surface density profile for gas, and Navarro–Frenk–White profile for dark matter halo; and employ the Markov chain Monte Carlo optimization method for parameter estimation. Further, to validate our analysis, we revisit important scaling relations, e.g. the M gas –M star relation, M star –M halo relation, M gas –M halo relation and Baryonic Tully–Fisher relation. Here, the scaling relations from our analysis are broadly consistent with that reported in the literature. A larger sample of galaxies from GARCIA in the near future will allow studying these scaling relations in greater details.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Real-space refinement in PHENIX for cryo-EM and crystallography

This article describes the implementation of real-space refinement in the phenix.real_space_refine program from the PHENIX suite. The use of a simplified refinement target function enables very fast calculation, which in turn makes it possible to identify optimal data-restraint weights as part of routine refinements with little runtime cost. Refinement of atomic models against low-resolution data benefits from the inclusion of as much additional information as is available. In addition to standard restraints on covalent geometry, phenix.real_space_refine makes use of extra information such as secondary-structure and rotamer-specific restraints, as well as restraints or constraints on internal molecular symmetry. The re-refinement of 385 cryo-EM-derived models available in the Protein Data Bank at resolutions of 6 Å or better shows significant improvement of the models and of the fit of these models to the target maps.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Antibody Responses during Hepatitis B Viral Infection

Hepatitis B is a DNA virus that infects liver cells and can cause both acute and chronic disease. It is believed that both viral and host factors are responsible for determining whether the infection is cleared or becomes chronic. Here we investigate the mechanism of protection by developing a mathematical model of the antibody response following hepatitis B virus (HBV) infection. We fitted the model to data from seven infected adults identified during acute infection and determined the ability of the virus to escape neutralization through overproduction of non-infectious subviral particles, which have HBs proteins on their surface, but do not contain nucleocapsid protein and viral nucleic acids. We showed that viral clearance can be achieved for high anti-HBV antibody levels, as in vaccinated individuals, when: (1) the rate of synthesis of hepatitis B subviral particles is slow; (2) the rate of synthesis of hepatitis B subviral particles is high but either anti-HBV antibody production is fast, the antibody affinity is high, or the levels of pre-existent HBV-specific antibody at the time of infection are high, as could be attained by vaccination. We further showed that viral clearance can be achieved for low equilibrium antiHBV antibody levels, as in unvaccinated individuals, when a strong cellular immune response controls early infection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Stellar Mass and Stellar Mass-to-light Ratio–Color Relations for Low Surface Brightness Galaxies

We estimate the stellar mass for a sample of low surface brightness galaxies (LSBGs) by fitting their multiband spectral energy distributions (SEDs) to the stellar population synthesis model. The derived stellar masses (log M {sub *}/M {sub ⊙}) span from 7.1 to 11.1, with a mean of log M {sub *}/M {sub ⊙} = 8.5, which is lower than that for normal galaxies. The stellar mass-to-light ratio (γ*) in each band varies little with the absolute magnitude but increases with higher M {sub *}. This trend of γ* with M {sub *} is even stronger in bluer bands. In addition, the γ* for our LSBGs slightly declines from the r band to the longer-wavelength bands. The log γ{sub ∗}{sup j} (j = g, r, i, and z) have relatively tight relations with optical colors of g − r and g − i. Compared with several representative γ*–color relations (MLCRs) from the literature, our MLCRs based on LSBG data are consistently among those literature MLCRs previously defined on diverse galaxy samples, and the existing minor differences between the MLCRs are caused by the differences in the SED model ingredients (including initial mass function, star formation history, and stellar population model), line fitting techniques, galaxy samples, and photometric zero-point, rather than the galaxy surface brightness itself, which distinguishes LSBGs from high surface brightness galaxies. Our LSBGs would be very likely to follow those representative MLCRs previously defined in diverse galaxy populations, if those main ingredients were taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗