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 433 records · Page 24

A precessing relativistic jet model for 3C 449

It is shown that the radio structure of 3C 449 can be matched with a model in which the jets are precessing and have relativistic (beta greater-than or equal to 0.4) velocities. The best-fit model implies a precession period of about 100,000 yr and a cone angle which increases with time. A similar model may be relevant for the radio structure of 3C 31. A brief discussion of the implications for 3C 449 is given.

Gower, A. C.↗

Seeing Double: ASASSN-18bt Exhibits a Two-component Rise in the Early-time K2 Light Curve

On 2018 February 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B(max)=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and prediscovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System. The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN Ia light curve, and it unambiguously shows a ∼4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a nondegenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe Ia light curves. Finally, we use Swift X-ray nondetections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant-density CSM, these nondetections constrain ρ<4.5×10(exp 5)per cu.cm at a radius of 4×10(exp 15) cm from the progenitor star. Assuming a wind-like environment, we place mass loss limits of M˙ < 8 x 10(exp -6) M(ʘ)per yr for v(w)=100 km/s, ruling out some symbiotic progenitor systems. This work highlights the power of well-sampled early-time data and the need for immediate multiband, high-cadence follow-up for progress in understanding SNe Ia.

B. J. Shappee↗

A Misperception in Reliability Growth Modelling

The Duane reliability growth model is n(t)/t = k t^-alpha (1) The reliability growth rate is alpha, the downward slope of n(t)/t versus t. It usually varies from 0.2 to 0.6. k is a constant. Crow used a 56-failure data set to illustrate reliability growth.1 A graphical Duane model fit to this data gives n(t)/t = 0.640 t^-0.283 (2) A problem in using the Duane-Crow reliability growth model is that it assumes that reliability growth continues and the failure rate decreases throughout the test period. It is more usual that reliability growth stops when the failure rated is low enough. Growth testing is often followed by testing with a low constant failure rate due to rare or uncorrectable failure modes. As more and more low constant rate acceptable failures accumulate after the period of reliability growth, the reliability growth time exponent alpha decreases toward zero. This occurs if constant rate failures are treated as occurring during the reliability growth period. It is more accurate to model a period of initial reliability growth followed by testing without repair to more accurately determine the final constant failure rate. This is done in the abcd model. n(t)/t = a t^-b + c from t = 0 to td (3) = c + d after td, where d = a td^-b (4) The term a t^-b describes the continuous reliability growth that continues out to time td and c is the constant uncorrected failure rate. The parameter d represents an additional constant failure rate due to correctable but uncorrected failure modes. After the reliability growth process is terminated, the failure rate n(t)/t = c + d.

Harry W Jones↗

Spectral and Temporal Properties of Galactic Black Hole Systems

Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).

Wheeler, J. Craig↗

Q -score as a reliability measure for protein, nucleic acid and small-molecule atomic coordinate models derived from 3DEM maps

Atomic coordinate models are important for the interpretation of 3D maps produced with cryoEM and cryoET (3D electron microscopy; 3DEM). In addition to visual inspection of such maps and models, quantitative metrics can inform about the reliability of the atomic coordinates, in particular how well the model is supported by the experimentally determined 3DEM map. A recently introduced metric, Q-score, was shown to correlate well with the reported resolution of the map for well fitted models. Here, we present new statistical analyses of Q-score based on its application to ∼10 000 maps and models archived in the EMDB (Electron Microscopy Data Bank) and PDB (Protein Data Bank). Further, we introduce two new metrics based on Q-score to represent each map and model relative to all entries in the EMDB and those with similar resolution. We explore through illustrative examples of proteins, nucleic acids and small molecules how Q-scores can indicate whether the atomic coordinates are well fitted to 3DEM maps and also whether some parts of a map may be poorly resolved due to factors such as molecular flexibility, radiation damage and/or conformational heterogeneity. These examples and statistical analyses provide a basis for how Q-scores can be interpreted effectively in order to evaluate 3DEM maps and atomic coordinate models prior to publication and archiving.

B factors↗

Partitioning controls on Amazon forest photosynthesis between environmental and biotic factors at hourly to interannual timescales

Gross ecosystem productivity (GEP) in tropical forests varies both with the environment and with biotic changes in photosynthetic infrastructure, but our understanding of the relative effects of these factors across timescales is limited. Here, we used a statistical model to partition the variability of seven years of eddy covariance-derived GEP in a central Amazon evergreen forest into two main causes: variation in environmental drivers (solar radiation, diffuse light fraction, and vapor pressure deficit) that interact with model parameters that govern photosynthesis and biotic variation in canopy photosynthetic light-use efficiency associated with changes in the parameters themselves. Our fitted model was able to explain most of the variability in GEP at hourly (R 2 = 0.77) to interannual (R 2 = 0.80) timescales. At hourly timescales, we found that 75% of observed GEP variability could be attributed to environmental variability. When aggregating GEP to the longer timescales (daily, monthly, and yearly), however, environmental variation explained progressively less GEP variability: At monthly timescales, it explained only 3%, much less than biotic variation in canopy photosynthetic light-use efficiency, which accounted for 63%. These results challenge modeling approaches that assume GEP is primarily controlled by the environment at both short and long timescales. Our approach distinguishing biotic from environmental variability can help to resolve debates about environmental limitations to tropical forest photosynthesis. For example, we found that biotically regulated canopy photosynthetic light-use efficiency (associated with leaf phenology) increased with sunlight during dry seasons (consistent with light but not water limitation of canopy development) but that realized GEP was nonetheless lower relative to its potential efficiency during dry than wet seasons (consistent with water limitation of photosynthesis in given assemblages of leaves). Lastly, this work highlights the importance of accounting for differential regulation of GEP at different timescales and of identifying the underlying feedbacks and adaptive mechanisms.

54 ENVIRONMENTAL SCIENCES↗

The Y-type Brown Dwarfs: Estimates of Mass and Age from New Astrometry, Homogenized Photometry, and Near-infrared Spectroscopy

The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer ( WISE ) led to the discovery of extremely cold, low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST , ESO, Spitzer , and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put onto the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that nonequilibrium atmospheric chemistry is important for these objects. Nonequilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 ≤ T {sub eff} (K) ≤ 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 ≤ T {sub eff} (K) ≤ 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the T {sub eff} = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5–8 Jupiter masses and an age of 0.3–6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6–8.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

High resolution, low temperature photoabsorption cross-section of C2H2 with application to Saturn's atmosphere

New laboratory observations of the VUV absorption cross-section of C2H2, obtained under physical conditions approximating stratospheres of the giant planets, were combined with IUE observations of the albedo of Saturn, for which improved data reduction techniques have been used, to produce new models for that atmosphere. When the effects of C2H2 absorption are accounted for, additional absorption by other molecules is required. The best-fitting model also includes absorption by PH3, H2O, C2H6 and CH4. A small residual disagreement near 1600 A suggests that an additional trace species may be required to complete the model.

Caldwell, John↗

The Distribution of Cosmic-Ray Sources in the Galaxy, Gamma-Rays and the Gradient in the CO-to-H2 Relation

We present a solution to the apparent discrepancy between the radial gradient in the diffuse Galactic gamma-ray emissivity and the distribution of supernova remnants, believed to be the sources of cosmic rays. Recent determinations of the pulsar distribution have made the discrepancy even more apparent. The problem is shown to be plausibly solved by a variation in the Wco-to-N(H2) scaling factor. If this factor increases by a factor of 5-10 from the inner to the outer Galaxy, as expected from the Galactic metallicity gradient and supported by other evidence, we show that the source distribution required to match the radial gradient of gamma-rays can be reconciled with the distribution of supernova remnants as traced by current studies of pulsars. The resulting model fits the EGRET gamma-ray profiles extremely well in longitude, and reproduces the mid-latitude inner Galaxy intensities better than previous models.

Strong, A. W.↗

Integrated Nondestructive Assay Systems to Estimate Plutonium in Spent Fuel Assemblies

An integrated nondestructive assay (NDA) system combining active (neutron generator) and passive neutron detection and passive gamma (PG) detection is being analyzed in order to estimate the amount of plutonium, verify initial enrichment, burnup, and cooling time, and detect partial defects in a spent fuel assembly (SFA). Active signals are measured using the differential die-away (DDA), delayed neutron (DN), and delayed gamma (DG) techniques. Passive signals are measured using total neutron (TN) counts and both gross and spectral resolved gamma counts. To quantify how a system of several NDA techniques is expected to perform, all of the relevant NDA techniques listed above were simulated as a function of various reactor conditions such as initial enrichment, burnup, cooling time, assembly shuffling pattern, reactor operating conditions (including temperature, pressure, and the presence of burnable poisons) by simulating the NDA response for five sets of light water reactor assemblies. In this study, we compare the performance of several exploratory model-fitting options (including neural networks, adaptive regression with splines, iterative bias reduction smoothing, projection pursuit regression, and regression with quadratic terms and interaction terms) to relate data simulated with measurement and model error effects from various subsets of the NDA techniques to the total Pu mass. Isotope masses for SFAs and expected detector responses (DRs) for several NDA techniques are simulated using MCNP, and the DRs become inputs to the fitting process. Such responses include eight signals from DDA, one from DN, one from TN, and up to seven from PG; the DG signal will be examined separately. Results are summarized using the root-mean-squared estimation error for plutonium mass in held-out subsets of the data for a range of model and measurement error variances. Different simulation assumptions lead to different spent fuel libraries relating DRs to Pu mass. Some results for training with one library and testing with another library are also given.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Modeling the Dynamics and Migratory Pathways of Virus-Specific Antibody-Secreting Cell Populations in Primary Influenza Infection

The B cell response to influenza infection of the respiratory tract contributes to viral clearance and establishes profound resistance to reinfection by related viruses. Numerous studies have measured virus-specific antibody-secreting cell (ASC) frequencies in different anatomical compartments after influenza infection and provided a general picture of the kinetics of ASC formation and dispersion. However, the dynamics of ASC populations are difficult to determine experimentally and have received little attention. Here, we applied mathematical modeling to investigate the dynamics of ASC growth, death, and migration over the 2-week period following primary influenza infection in mice. Experimental data for model fitting came from high frequency measurements of virus-specific IgM, IgG, and IgA ASCs in the mediastinal lymph node (MLN), spleen, and lung. Model construction was based on a set of assumptions about ASC gain and loss from the sampled sites, and also on the directionality of ASC trafficking pathways. Most notably, modeling results suggest that differences in ASC fate and trafficking patterns reflect the site of formation and the expressed antibody class. Essentially all early IgA ASCs in the MLN migrated to spleen or lung, whereas cell death was likely the major reason for IgM and IgG ASC loss from the MLN. In contrast, the spleen contributed most of the IgM and IgG ASCs that migrated to the lung, but essentially none of the IgA ASCs. This finding points to a critical role for regional lymph nodes such as the MLN in the rapid generation of IgA ASCs that seed the lung. Results for the MLN also suggest that ASC death is a significant early feature of the B cell response. Overall, our analysis is consistent with accepted concepts in many regards, but it also indicates novel features of the B cell response to influenza that warrant further investigation.

59 BASIC BIOLOGICAL SCIENCES↗

PANCHROMATIC IMAGING OF A TRANSITIONAL DISK: THE DISK OF GM AUR IN OPTICAL AND FUV SCATTERED LIGHT

We have imaged GM Aurigae with the Hubble Space Telescope , detected its disk in scattered light at 1400 and 1650 Å, and compared these with observations at 3300 Å, 5550 Å, 1.1 μ m, and 1.6 μ m. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 Å shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that the surface of the entire disk is populated with submicron grains. We have compiled a spectral energy distribution from 0.1 μ m to 1 mm and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that could be due to a planet of mass <9 M {sub J} interior to 24 au. A unique cylindrical structure is detected in the far-UV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb brightened and extends 190 ± 35 au above the disk midplane. The inner radius of the limb brightening is 40 ± 10 au, just beyond the submillimeter cavity wall.

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting biomass of hyperdiverse and structurally complex central Amazonian forests - a virtual approach using extensive field data

Notice on corrigendum: This paper has a corresponding corrigendum published. Please read the corrigendum first. Old-growth forests are subject to substantial changes in structure and species composition due to the intensification of human activities, gradual climate change and extreme weather events. Trees store ca. 90 % of the total aboveground biomass (AGB) in tropical forests and precise tree biomass estimation models are crucial for management and conservation. In the central Amazon, predicting AGB at large spatial scales is a challenging task due to the heterogeneity of successional stages, high tree species diversity and inherent variations in tree allometry and architecture. We parameterized generic AGB estimation models applicable across species and a wide range of structural and compositional variation related to species sorting into height layers as well as frequent natural disturbances. We used 727 trees (diameter at breast height ≥ 5 cm) from 101 genera and at least 135 species harvested in a contiguous forest near Manaus, Brazil. Sampling from this data set we assembled six scenarios designed to span existing gradients in floristic composition and size distribution in order to select models that best predict AGB at the landscape level across successional gradients. We found that good individual tree model fits do not necessarily translate into reliable predictions of AGB at the landscape level. When predicting AGB (dry mass) over scenarios using our different models and an available pantropical model, we observed systematic biases ranging from -31 % (pantropical) to +39 %, with root-mean-square error (RMSE) values of up to 130 Mg ha -1 (pantropical). Our first and second best models had both low mean biases (0.8 and 3.9 %, respectively) and RMSE (9.4 and 18.6 Mg ha -1 ) when applied over scenarios. Predicting biomass correctly at the landscape level in hyperdiverse and structurally complex tropical forests, especially allowing good performance at the margins of data availability for model construction/calibration, requires the inclusion of predictors that express inherent variations in species architecture. The model of interest should comprise the floristic composition and size-distribution variability of the target forest, implying that even generic global or pantropical biomass estimation models can lead to strong biases. Reliable biomass assessments for the Amazon basin (i.e., secondary forests) still depend on the collection of allometric data at the local/regional scale and forest inventories including species-specific attributes, which are often unavailable or estimated imprecisely in most regions.

54 ENVIRONMENTAL SCIENCES↗

Panchromatic Imaging of a Transitional Disk: The Disk of GM AUR in Optical and FUV Scattered Light

We have imaged GM Aurigae with the Hubble Space Telescope, detected its disk in scattered light at 1400 and 1650 A, and compared these with observations at 3300 A, 5550 A, 1.1 micron, and 1.6 micron. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 A shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that the surface of the entire disk is populated with submicron grains. We have compiled a spectral energy distribution from 0.1 micron to 1 mm and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that could be due to a planet of mass less than 9M(sub J) interior to 24 au. A unique cylindrical structure is detected in the far- UV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb brightened and extends 190 +/- 35 au above the disk midplane. The inner radius of the limb brightening is 40 +/- 10 au, just beyond the submillimeter cavity wall.

circumstellar matter aEuro" protoplanetary disks a↗

Shot model parameters for Cygnus X-1 through phase portrait fitting

Shot models for systems having about 1/f power density spectrum are developed by utilizing a distribution of shot durations. Parameters of the distribution are determined by fitting the power spectrum either with analytic forms for the spectrum of a shot model with a given shot profile, or with the spectrum derived from numerical realizations of trial shot models. The shot fraction is specified by fitting the phase portrait, which is a plot of intensity at a given time versus intensity at a delayed time and in principle is sensitive to different shot profiles. These techniques have been extensively applied to the X-ray variability of Cygnus X-1, using HEAO 1 A-2 and an Exosat ME observation. The power spectra suggest models having characteristic shot durations lasting from milliseconds to a few seconds, while the phase portrait fits give shot fractions of about 50 percent. Best fits to the portraits are obtained if the amplitude of the shot is a power-law function of the duration of the shot. These fits prefer shots having a symmetric exponential rise and decay. Results are interpreted in terms of a distribution of magnetic flares in the accretion disk.

Lochner, James C.↗

Supersymmetric B – L inflation near the conformal coupling

We investigate a novel scenario of cosmological inflation in a gauged B-L extended minimal supersymmetric Standard Model with R-symmetry. We use a noncanonical Kähler potential and a superpotential, both preserving the R-symmetry to construct a model of slow-roll inflation. The model is controlled by two real parameters: the nonminimal coupling ξ that originates from the Kähler potential, and the breaking scale v of the U (1)B-L symmetry. We compute the spectrum of the cosmological microwave background radiation and show that the prediction of the model fits well the recent Planck satellite observation for a wide range of the parameter space. We also find that the typical reheating temperature of the model is low enough to avoid the gravitino problem but nevertheless allows sufficient production of the baryon asymmetry if we take into account the effect of resonance enhancement. The model is free from cosmic strings that impose stringent constraints on generic U (1)B-L based scenarios, as in our scenario the U (1)B-L symmetry is broken from the onset.

Arai, Masato↗

Empirical modeling of dopability in diamond-like semiconductors

Carrier concentration optimization has been an enduring challenge when developing newly discovered semiconductors for applications (e.g., thermoelectrics, transparent conductors, photovoltaics). This barrier has been particularly pernicious in the realm of high-throughput property prediction, where the carrier concentration is often assumed to be a free parameter and the limits are not predicted due to the high computational cost. In this work, we explore the application of machine learning for high-throughput carrier concentration range prediction. Bounding the model within diamond-like semiconductors, the learning set was developed from experimental carrier concentration data on 127 compounds ranging from unary to quaternary. The data were analyzed using various statistical and machine learning methods. Accurate predictions of carrier concentration ranges in diamond-like semiconductors are made within approximately one order of magnitude on average across both p- and n-type dopability. The model fit to empirical data is analyzed to understand what drives trends in carrier concentration and compared with previous computational efforts. Lastly, dopability predictions from this model are combined with high-throughput quality factor predictions to identify promising thermoelectric materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗