Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Slitless spectrophotometry with forward modelling: Principles and application to measuring atmospheric transmission

Context.In the next decade, many optical surveys will aim to answer the question of the nature of dark energy by measuring its equation-of-state parameter at the per mill level. This requires trusting the photometric calibration of the survey with a precision never reached so far on many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or for each season or for the full survey on average, can help reach the per mill precision for the magnitudes. Aims.This work aims at proving the ability to use slitless spectroscopy for standard-star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by theVera C. RubinObservatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in theRubinAuxiliary Telescope. Methods.The theoretical basis of slitless spectrophotometry is at the heart of our forward-model approach to extract spectroscopic information from slitless data. We developed a publicly available software calledSpectractor, which implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to obtain the spectrum. Results.We show through simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data by solving specific issues and illustrating that our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from the data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure in simulations and test it on the limited available data set.

Astronomy & Astrophysics↗
DOE OSTI2025

Quantifying dispersity in size and shape of nanoparticles from small-angle scattering data using machine learning based CREASE

Here, we use machine learning (ML) enhanced computational reverse engineering analysis of scattering experiments (CREASE) to interpret small-angle X-ray scattering (SAXS) data obtained from a system of nanoparticles without a priori knowledge of their exact shapes (e.g. spheres or ellipsoids), sizes (0.5–50 nm) and distributions. The SAXS measurements yielded three categories of scattering profiles exhibiting 'strong', 'weak' and 'no' features. Diminishing features (e.g. broadening or disappearing peaks) in scattering profiles have always been attributed to the presence of significant dispersity in the system. Such featureless SAXS data are not suitable for traditional analysis using analytical models. If one were to fit a relevant analytical model (e.g. the lmfit analytical model for polydisperse spheres) to these 'weak' and 'no' SAXS profiles from our nanoparticle systems, one would obtain non-unique interpretations of the data. Relying on electron microscopy to identify the distributions of nanoparticle shapes and sizes is also unfeasible, especially in high-throughput synthesis and characterization loops. In such situations, to identify the distributions of particle sizes and shapes that could be present in the sample, one must rely on methods like ML-CREASE to interpret the data quickly and output all relevant interpretations about the structure present in the system. The ML-CREASE optimization loop takes the experimental scattering profile as input and outputs multiple candidate solutions whose computed scattering profiles match the SAXS profile input. The ML-CREASE method outputs distributions of relevant structural features, such as the volume fraction of the nanoparticles in the system and the mean and standard deviation of the particle size and aspect ratio, assuming a type of distribution (e.g. normal, log-normal) for size and aspect ratio. We find that, for the SAXS profiles analyzed here, accounting for the shape dispersity along with size dispersity of the nanoparticles using ML-CREASE improved the match between the computed scattering profiles and input experimental profiles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Buzzard to Cardinal: Improved Mock Catalogs for Large Galaxy Surveys

Abstract We present the Cardinal mock galaxy catalogs, a new version of the Buzzard simulation that has been updated to support ongoing and future cosmological surveys, including the Dark Energy Survey (DES), DESI, and LSST. These catalogs are based on a one-quarter sky simulation populated with galaxies out to a redshift of z = 2.35 to a depth of m r = 27. Compared to the Buzzard mocks, the Cardinal mocks include an updated subhalo abundance matching model that considers orphan galaxies and includes mass-dependent scatter between galaxy luminosity and halo properties. This model can simultaneously fit galaxy clustering and group–galaxy cross-correlations measured in three different luminosity threshold samples. The Cardinal mocks also feature a new color assignment model that can simultaneously fit color-dependent galaxy clustering in three different luminosity bins. We have developed an algorithm that uses photometric data to further improve the color assignment model and have also developed a novel method to improve small-scale lensing below the ray-tracing resolution. These improvements enable the Cardinal mocks to accurately reproduce the abundance of galaxy clusters and the properties of lens galaxies in the DES data. As such, these simulations will be a valuable tool for future cosmological analyses based on large sky surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Deuterium Permeation Through a Self-Supported Palladium-Silver Membrane in Helium Gas Mixtures

Palladium membranes and membrane reactors can separate and purify tritium from impurities in the plasma exhaust processing system for the fusion energy fuel cycle. Membranes can also act as a continuous separation method to remove tritium from helium streams in the breeder blanket tritium extraction system, such as from the purge gas of solid breeders. To investigate the potential of these membranes for these applications, we performed a deuterium permeation campaign with a self-supported palladium-silver (Pd-25Ag wt%) tube of 0.15 m length, 76 µ m wall thickness, and 3.0 mm inner diameter. A gas mixture of 3.95% D 2 and a balance of He was delivered to the inside of the tube and D permeated radially outwards through the membrane into a vacuum chamber. Permeation experiments were conducted at 300 ° C, 350 ° C, 400° C, and 450°C, from 100 to 1000 sccm total flow rate, and with total pressures of 90, 150, 190, and 250 kPa. Further, deuterium permeation was consistently lower than predicted from diffusion-limited permeation models, thus we developed a transport model that included gas-phase mass transfer and surface reactions to model experimental results. The dissociation constant was optimized to fit the developed model to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evaluating Cell Temperature Models and the Effect of Wind Speed in PV System Capacity Testing: Preprint

Capacity testing is a routine procedure for assessing a photovoltaic system's performance relative to expectations. The most common test method involves fitting a regression model that predicts system output power using operating weather conditions including wind speed. Structural modifications to the regression model to incorporate wind in different ways improved the model's ability to fit measured system performance, but the observed improvements were small and unlikely to change the result of a capacity test. However, the results showed that the choice of reporting wind speed and inclusion or exclusion of wind speed in the performance model used as the test benchmark can significantly change the test result.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Searches for New Physics With Muon Conversion at Fermilab and Triboson Production at the LHC

We report on several efforts to search for physics beyond the standard model of particle physics at broad energy scales. The Mu2e experiment at Fermilab will search for charged lepton flavor violation via the muon to electron conversion process, which is suppressed in the Standard Model. Mu2e will be operated at a low energy, yet can probe New Physics at very high mass scales (O(1e3 - 1e4 ) TeV). At high energies, the CMS experiment at the CERN LHC continues to deliver an impressive suite of Standard Model measurements and limits on a variety of New Physics signatures. Mu2e is under construction and slated to collect its first physics data in the coming years. This thesis describes work done during the construction phase of Mu2e and focuses on two critical areas: magnetic field modeling and statistical analysis. We describe a novel method for field modeling which we validate using a simulated dataset representing the expected magnetic field in the Detector Solenoid. This method blends a standard least-squares fitting technique that utilizes physically motivated analytical model functions with a novel physics informed network that is constructed to obey Maxwell’s equations. We show the technique can model the field with an accuracy of 10−7 despite the presence of injected noise in the pseudo-measurements at the 10−5 level. We then present preliminary results of the calibration of 3D Hall probes at the sub-10−4 level. These probes will be used to directly measure the Mu2e Detector Solenoid magnetic field on a sparse grid; these measurements serve as the input to the field model fitting. Finally, we describe the first implementation of both an unbinned shape analysis and a Bayesian interpretation applied to Mu2e pseudo-data. Up to 20% tighter limits can be set by the shape analysis compared to a standard cut & count analysis. The AlCap experiment collected data at PSI in 2015 to measure several important quantities related to nuclear muon capture on an aluminum target, which is a significant background process for Mu2e. The neutron emission from muon capture can introduce background hits in the Mu2e detectors and can increase radiation damage in various elements of the apparatus. We present measurements of the neutron group fluence and mean neutron multiplicity for muon capture on aluminum nuclei. Finally, we discuss an analysis of triboson production at CMS using an Effective Field Theory framework. Standard Model triboson production, which was first observed at CMS in 2020, has a relatively small cross section and provides direct access to both anomalous triple gauge couplings and quartic gauge couplings. These couplings, interpreted in the Standard Model Effective Field Theory, are studied in the present work. We target the boosted regime where the background rate is low and yields are enhanced when dimension-6 and dimension-8 Wilson coefficients are non-zero. We do not observe an excess in the data and therefore set bounds on the Wilson coefficients. For dimension-6 coefficients the tightest observed (expected) bounds are set on cW /Λ2 where Λ is the mass scale of new physics; the bounds are [−0.13, 0.12] TeV−2 ([−0.12, 0.12] TeV−2 ) at 95% CL. The tightest bounds in dimension-8 are set on fT,0 / Λ4 ; the observed (expected) bounds at 95% CL are [−0.63, 0.69] TeV−4 ([−0.54, 0.62] TeV−4 ). Additional results are presented which include scenarios where multiple Wilson coefficients are non-zero, the application of signal model clipping to address unitarity violation in Effective Field Theories, and a novel template fit developed for easier reinterpretation of our results.

Kampa, Cole Erik [Northwestern U. (main)] (ORCID:0↗
DOE OSTI2023

Single-Point Extrapolation to the Complete Basis Set Limit through Deep Learning

Machine learning (ML) offers an attractive method for making predictions about molecular systems while circumventing the need to run expensive electronic structure calculations. Once trained on ab initio data, the promise of ML is to deliver accurate predictions of molecular properties that were previously computationally infeasible. In this work, we develop and train a graph neural network model to correct the basis set incompleteness error (BSIE) between a small and large basis set at the RHF and B3LYP levels of theory. Our results show that, when compared to fitting to the total potential, an ML model fitted to correct the BSIE is better at generalizing to systems not seen during training. We test this ability by training on single molecules while evaluating on molecular complexes. We also show that ensemble models yield better behaved potentials in situations where the training data is insufficient. However, even when only fitting to the BSIE, acceptable performance is only achieved when the training data sufficiently resemble the systems one wants to make predictions on. As a result, the test error of the final model trained to predict the difference between the cc-pVDZ and cc-pV5Z potential is 0.184 kcal/mol for the B3LYP density functional, and the ensemble model accurately reproduces the large basis set interaction energy curves on the S66x8 dataset.

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

An Initial Microstructurally Informed Model of High Burnup Structure Formation in UO 2 Fuel

The microstructure of a UO 2 fuel pellet changes as burnup increases, impacting fuel performance. Predicting and characterizing high burnup structure (HBS) and dark zone formation is a key part of supporting burnup limit extensions for light water reactors. This paper describes a model developed through fitting radially resolved pellet data obtained from recently published microstructural characterization data. The model predicts grain size and grain character, in addition to pore density and size, with fitting dependencies on power history variables. Separately fitting power history variables to microstructural parameters allows for insight into the underlying physical phenomena for future model development. Additionally, experimental data have been correlated to an HBS fraction to facilitate the development of a model capable of predicting a total fuel restructured fraction at the engineering scale. In conclusion, this two-step approach provides a coupling from reactor power history to microstructural data to fractional HBS and creates a basis to model HBS-dependent parameters in a fuel performance code.

High burnup structure↗
DOE OSTI2020

Temporal covariance structure of multi-spectral phenotypes and their predictive ability for end-of-season traits in maize

Vegetation indices (VIs) derived from multi-spectral imaging (MSI) platforms can be used to study properties of crop canopy, providing non-destructive phenotypes that could be used to better understand growth curves throughout the growing season. To investigate the amount of variation present in several VIs and their relationship with important end-of-season traits, genetic and residual (co)variances for VIs, grain yield and moisture were estimated using data collected from maize hybrid trials. The VIs considered were Normalized Difference Vegetation Index (NDVI), Green NDVI, Red Edge NDVI, Soil-Adjusted Vegetation Index, Enhanced Vegetation Index and simple Ratio of Near Infrared to Red (Red) reflectance. Genetic correlations of VIs with grain yield and moisture were used to fit multi-trait models for prediction of end-of-season traits and evaluated using within site/year cross-validation. To explore alternatives to fitting multiple phenotypes from MSI, random regression models with linear splines were fit using data collected in 2016 and 2017. Heritability estimates ranging from (0.10 to 0.82) were observed, indicating that there exists considerable amount of genetic variation in these VIs. Furthermore, strong genetic and residual correlations of the VIs, NDVI and NDRE, with grain yield and moisture were found. Considerable increases in prediction accuracy were observed from the multi-trait model when using NDVI and NDRE as a secondary trait. Finally, random regression with a linear spline function shows potential to be used as an alternative to mixed models to fit VIs from multiple time points.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Hanford low-activity waste glass composition-temperature-melt viscosity relationships

This study developed a model for predicting viscosity of alkali-alumino-borosilicate glass melts as functions of composition and temperature. The model is based on a total of 3935 viscosity-temperature data from 574 glasses with viscosity values ranging from 2.53 to 7260 Poise (P) in the temperature range of 900–1260°C. Several different model forms were surveyed, including those based on Arrhenius, Vogel-Fulcher-Tammann, Avramov-Milchev, and Mauro-Yue-Ellison-Gupta-Allen. For each of these models, combinations of the temperature-independent parameters were fitted to composition. It was found that generally fitting more than one temperature-independent parameter as functions of composition resulted in overfitting. The Avramov-Milchev-based model was found to best represent the Hanford low-activity waste glass melt viscosity data based on model fit and validation statistics. A 21-term partial quadratic mixture model was recommended for use. This model predicts melt viscosity with a root-mean square error of .1736 ln(P), which is similar to the error in viscosity measurements from replicate glass analyses of .1383 ln(P). Viscosity was found to be most increased by SiO 2 > Al 2 O 3 > ZrO 2 > SnO 2 and most decreased by Li 2 O > Na 2 O > B 2 O 3 > CaO > K 2 O > MgO, at temperatures from 900 to 1260°C.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Direct measurement of storage and loss behavior in AFM force–distance experiments using the modified Fourier transformation

Force–distance curve experiments are commonly performed in atomic force microscopy (AFM) to obtain the viscoelastic characteristics of materials, such as the storage and loss moduli or compliances. The classic methods used to obtain these characteristics consist of fitting a viscoelastic material model to the experimentally obtained AFM data. Here, we demonstrate a new method that utilizes the modified discrete Fourier transform to approximate the storage and loss behavior of a material directly from the data, without the need for a fit. Additionally, one may still fit a model to the resulting storage and loss behavior if a parameterized description of the material is desired. In contrast to fitting the data to a model chosen a priori, departing from a model-free description of the material's frequency behavior guides the selection of the model, such that the user may choose the one that is most appropriate for the particular material under study. To this end, we also include modified Fourier domain descriptions of commonly used viscoelastic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

eDNAjoint: An R package for interpreting paired or semi‐paired environmental DNA and traditional survey data in a Bayesian framework

Abstract Environmental DNA (eDNA) sampling is increasingly used in surveys of species distribution as a potentially sensitive and efficient monitoring method. Yet access to modelling tools designed specifically for interpreting this new data type lags behind its ubiquity. While occupancy modelling software has dominated the analytical landscape for eDNA data analysis of single species, this type of model may not always be the most appropriate. The rate of eDNA detection often corresponds to species density, rather than just occupancy, and researchers often have access to observations from non‐genetic sampling methods at the same sites. To provide users access to a modelling framework designed to maximize the use of all available data, we developed an R package, eDNAjoint . The package provides an easy‐to‐use interface for fitting a ‘joint’ model that integrates data from paired or semi‐paired eDNA and traditional surveys in a Bayesian framework. The model can be used to estimate parameters like the probability of a false positive eDNA detection and mean catch rate at a site, and the package allows access to multiple model variations and Bayesian prior customization. Additional functionality can be used for model selection, summarising posteriors and comparing the relative sensitivities of the two survey methods. We demonstrate the use of eDNAjoint by fitting a variation of the model with site‐level covariates that scale the sensitivity of eDNA sampling relative to traditional sampling. The example workflow uses binary eDNA and seine count data for the endangered tidewater goby ( Eucyclogobius newberryi ) from a study by Schmelzle and Kinziger (2016). This use case includes a prior sensitivity analysis and an evaluation of the relationship between detection rates and environmental variables. eDNAjoint has the potential to greatly increase the range of users who will be able to rigorously analyse eDNA and traditional survey data in a Bayesian framework, understand if and how eDNA can improve monitoring practices, and gain confidence in the interpretability of eDNA data.

Keller, Abigail G. [Department of Environment Scie↗
DOE OSTI2024

Steady-state fuel performance analyses for the preliminary fuel concept of general atomics fast modular reactor

Here this manuscript presents the fuel performance analysis results of the General Atomics Fast Modular Reactor (FMR) based on an axi-symmetric (2D-RZ) geometry. Three fuel performance model sets that fit the FMR fuel specifications best, i.e., a BISON baseline model set, a BISON diffusion enhancement model set, and a BISON-FASTGRASS model set, were identified and evaluated against a series of relevant experimental cases featuring high burnup and low irradiation temperature conditions. The three BISON-based model sets were then utilized to conduct a comprehensive fuel performance analysis of the FMR fuel under normal operation including the shutdown/restarting periods for refueling. The evaluation of the fuel performance parameters, represented by temperature, internal pressure, stress, and strain, shows that the FMR fuel maintains its thermal and mechanical integrity during normal operation. Technology gaps and limitations are also discussed to guide future efforts for extending the performance analysis to transient scenarios as well as improving the fuel performance evaluation through experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Vibrational spectroscopy of shock compressed condensed phase nitrous oxide: Frequency shifts in the ν 1 mode

Over many years, we have published vibrational spectra of molecular materials at extreme pressure and temperature obtained using shock compression and coherent Raman spectroscopy methods. For diatomic molecules, we were able to extract vibrational temperatures from the intensities of spectrally resolved vibrational hot bands. Larger molecules in the condensed phase suffer band broadening effects that obscure the vibrational hot bands as the anharmonicities are typically smaller than the widths of the bands. This inability to resolve the hot bands inhibits the extraction of vibrational temperature and the measurement of the vibrational frequency of the fundamental. Here, we use a hot band model based on gas-phase anharmonic coupling coefficients to fit coherent anti-Stokes Raman spectra of the ν 1 vibrational mode of shock compressed condensed phase N 2 O with shock pressures and temperatures estimated from literature equations of state and compare to fits from a model using a single Gaussian peak. We report the resulting vibrational frequency shifts with shock pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Highly selective Si 3 N 4 /SiO 2 etching using an NF 3 /N 2 /O 2 /H 2 remote plasma. II. Surface reaction mechanism

Developing processes for highly selective etching of silicon nitride (Si 3 N 4 ) with respect to silicon dioxide (SiO 2 ) is a major priority for semiconductor fabrication processing. In this paper and in Paper I [Volynets et al., J. Vac. Sci. Technol. A 38, 023007 (2020)], mechanisms are discussed for highly selective Si3N4 etching in a remote plasma based on experimental and theoretical investigations. The Si 3 N 4 /SiO 2 etch selectivity of up to 380 was experimentally produced using a remote plasma sustained in NF 3 /N 2 /O 2 /H 2 mixtures. A selectivity strongly depends on the flow rate of H 2 , an effect attributed to the formation of HF molecules in vibrationally excited states that accelerate etching reactions. Based on experimental measurements and zero-dimensional plasma simulations, an analytical etching model was developed for etch rates as a function of process parameters. Reaction rates and sticking coefficients were provided by quantum chemistry models and also fitted to the experimental results. Etch rates from the analytical model show good agreement with the experimental results and demonstrate why certain etchants accelerate or inhibit the etch process. In particular, the modeling shows the important role of HF molecules in the first vibrationally excited state [HF( v = 1)] in achieving high Si 3 N 4 /SiO 2 selectivity.

42 ENGINEERING↗
DOE OSTI2025

Data and code from: Multivariate bayesian regression model for predicting disposed ash composition at U.S. coal fired power stations

This dataset contains the code and data files needed for implementation of a Multivariate Bayesian Regression model, described in Jin et al. (2025), for the historical prediction of the chemical composition of disposed coal ash at U.S. coal fired power plants as a function of annualized coal purchase data. The integrated coal supply data file (CoalSupplyDataset.csv) represents a compilation of monthly fuel purchase records for the period 1973-2022 at major U.S. power stations. These records were obtained from the U.S. Energy Information Administration. The CSV file also contains, for each coal purchase record, the coal region of the mine as defined by the U.S. Geological Survey. Data entry errors and data gaps in the EIA records were corrected as described in Jin et al. This CSV file represents the integrated coal supply data after corrections were made. The model structure and fitting parameters are encoded in pickle file format (Bayesian.pkl). The model was developed with the coal supply data and coal ash composition data, apportioned according to the Stratified Shuffle Split for training and testing subsets. The model was built using Python and the PyMC library. Reference Publication: Jin, Z.; Huang, J.; Hower, J.C.; Hsu-Kim, H.(2025). Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites. Environmental Science & Technology.

Coal ash composition↗
DOE OSTI2020

Modelling the Milky Way – I. Method and first results fitting the thick disc and halo with DES-Y3 data

ABSTRACT We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from trilegal models to real data. We apply this technique, which we call mwfitting, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ∼2300 deg2. We further explore a smaller subset (∼1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of mwfitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗