Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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 955 records · Page 53

8760-Based Method for Representing Variable Generation Capacity Value in Capacity Expansion Models

Capacity expansion models (CEMs) are widely used to evaluate the least-cost portfolio of electricity generators, transmission, and storage needed to reliably serve load over many years or decades. CEMs can be computationally complex and are often forced to estimate key parameters using simplified methods to achieve acceptable solve times or for other reasons. In this paper, we discuss one of these parameters -- capacity value (CV). We first provide a high-level motivation for and overview of CV. We next describe existing modeling simplifications and an alternate approach for estimating CV that utilizes hourly '8760' data of load and VG resources. We then apply this 8760 method to an established CEM, the National Renewable Energy Laboratory's (NREL's) Regional Energy Deployment System (ReEDS) model (Eurek et al. 2016). While this alternative approach for CV is not itself novel, it contributes to the broader CEM community by (1) demonstrating how a simplified 8760 hourly method, which can be easily implemented in other power sector models when data is available, more accurately captures CV trends than a statistical method within the ReEDS CEM, and (2) providing a flexible modeling framework from which other 8760-based system elements (e.g., demand response, storage, and transmission) can be added to further capture important dynamic interactions, such as curtailment.

14 SOLAR ENERGY↗

A steady-state alignment front in an accretion disk subjected to Lense–Thirring torques

Using only physical mechanisms, i.e., 3D magnetohydrodynamics (MHD) with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian gravitomagnetism. In this simulation, we have shown how, in the presence of MHD turbulence, a time-steady transition can be achieved between an inner disk region aligned with the equatorial plane of the central mass’s spin and an outer region orbiting in a different plane. The position of the equilibrium orientation transition is determined by a balance between gravitomagnetic torque and warp-induced inward mixing of misaligned angular momentum from the outer disk. If the mixing is interpreted in terms of diffusive transport, the implied diffusion coefficient is ≃(0.6–0.8)c{sub s}{sup 2}/Ω for sound speed c{sub s} and orbital frequency Ω. This calibration permits estimation of the orientation transition’s equilibrium location given the central mass, its spin parameter, and the disk’s surface density and scaleheight profiles. However, the alignment front overshoots before settling into an equilibrium, signaling that a diffusive model does not fully represent the time-dependent properties of alignment fronts under these conditions. Because the precessional torque on the disk at the alignment front is always comparable to the rate at which misaligned angular momentum is brought inward to the front by warp-driven radial motions, no break forms between the inner and outer portions of the disk in our simulation. Our results also raise questions about the applicability to MHD warped disks of the traditional distinction between “bending wave” and “diffusive” regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Combined Experimental and Computational Approach to Predict the Glass-Water Reaction

The use of mineral and glass dissolution rates measured in laboratory experiments to predict the weathering of primary minerals and volcanic and nuclear waste glasses in field studies requires the construction of rate models that accurately describe the weathering process over geologic timescales. Additionally, the need to model the long-term behavior of nuclear waste glass for the purpose of estimating radionuclide release rates requires that rate models be validated with long-term experiments. Several long-term test methods have been developed to accelerate the glass-water reaction [drip test, vapor hydration test, product consistency test B, and pressurized unsaturated flow (PUF)], thereby reducing the duration required to evaluate long-term performance. Currently, the PUF test is the only method that mimics the unsaturated hydraulic properties expected in a subsurface disposal facility and simultaneously monitors the glass-water reaction. PUF tests are being conducted to accelerate the weathering of glass and validate the model parameters being used to predict long-term glass behavior. A one-dimensional reactive chemical transport simulation of glass dissolution and secondary phase formation during a 1.5-year-long PUF experiment was conducted with the Subsurface Transport Over Reactive Multiphases (STORM) code. Results show that parameterization of the computer model by combining direct bench scale laboratory measurements and thermodynamic data provides an integrated approach to predicting glass behavior over the length of the experiment. Over the 1.5-year-long test duration, the rate decreased from 0.2 to 0.01 g/(m2 day) based on B release for low-activity waste glass LAWA44. The observed decrease is approximately two orders of magnitude higher than the decrease observed under static conditions with the SON68 glass (estimated to be a decrease by four orders of magnitude) and suggests that the gel-layer properties are less protective under these dynamic conditions.

aluminoborosilicate nuclear waste glass↗

Dynamic Model of Basic Oxygen Steelmaking Process Based on Multizone Reaction Kinetics: Modeling of Decarburization

In a previous study by the authors (Rout et al. in Metall Mater Trans B 49:537–557, 2018), a dynamic model for the BOF, employing the concept of multizone kinetics was developed. In the current study, the kinetics of decarburization reaction is investigated. The jet impact and slag–metal emulsion zones were identified to be primary zones for carbon oxidation. The dynamic parameters in the rate equation of decarburization such as residence time of metal drops in the emulsion, interfacial area evolution, initial size, and the effects of surface-active oxides have been included in the kinetic rate equation of the metal droplet. A modified mass-transfer coefficient based on the ideal Langmuir adsorption equilibrium has been proposed to take into account the surface blockage effects of SiO{sub 2} and P{sub 2}O{sub 5} in slag on the decarburization kinetics of a metal droplet in the emulsion. Further, a size distribution function has been included in the rate equation to evaluate the effect of droplet size on reaction kinetics. The mathematical simulation indicates that decarburization of the droplet in the emulsion is a strong function of the initial size and residence time. A modified droplet generation rate proposed previously by the authors has been used to estimate the total decarburization rate by slag–metal emulsion. The model’s prediction shows that about 76 pct of total carbon is removed by reactions in the emulsion, and the remaining is removed by reactions at the jet impact zone. The predicted bath carbon by the model has been found to be in good agreement with the industrially measured data.

36 MATERIALS SCIENCE↗

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗

The XXL Survey: XVII. X-ray and Sunyaev-Zel’dovich properties of the redshift 2.0 galaxy cluster XLSSC 122

We present results from a 100 ks XMM-Newton observation of galaxy cluster XLSSC 122, the first massive cluster discovered through its X-ray emission at z ≈ 2. The data provide the first precise constraints on the bulk thermodynamic properties of such a distant cluster, as well as an X-ray spectroscopic confirmation of its redshift. We measure an average temperature of kT = 5.0 ± 0.7 keV; a metallicity with respect to solar of Z/Z ⊙ = 0:33 -0.17 +0.19 , consistent with lower-redshift clusters; and a redshift of z = 1.99 -0.07 -0.06 , consistent with the earlier photo-z estimate. The measured gas density profile leads to a mass estimate at r 500 of M 500 = (6.3 ± 1.5) x 10 13 M ⊙ . From CARMA 30 GHz data, we measure the spherically integrated Compton parameter within r 500 to be Y 500 = (3.6 ± 0.4) x 10 -12 . We compare the measured properties of XLSSC 122 to lower-redshift cluster samples, and find good agreement when assuming the simplest (self-similar) form for the evolution of cluster scaling relations. While a single cluster provides limited information, this result suggests that the evolution of the intracluster medium in the most massive, well developed clusters is remarkably simple, even out to the highest redshifts where they have been found. At the same time, our data reaffirm the previously reported spatial offset between the centers of the X-ray and SZ signals for XLSSC 122, suggesting a disturbed configuration. Higher spatial resolution data could thus provide greater insights into the internal dynamics of this system.

79 ASTRONOMY AND ASTROPHYSICS↗

The electro-mechanical effect from charge dynamics on polymeric insulation lifetime

For polymeric material used as electrical insulation, the presence of space charges could be the consequence of material degradations that are thermally activated but increased by the application of an electric field. The dynamics of space charge, therefore, can be potentially used to characterize the material. In this direction, a new aging model in which parameters have clear physical meanings has been developed and applied to the material to extrapolate the lifetime. The kinetic equation has been established based on charge trapping and detrapping of the injected charge from the electrodes. The local electromechanical energy stored in the region surrounding the trap is able to reduce the trap-depth with a value related to the electric field. At a level where the internal electric field exceeds the detrapping field in the material, an electron can be efficiently detrapped and the released energy from detrapping process can cause a weak bond or chain scission i.e. material degradation. The model has been applied to the electro-thermally aged low density polyethylene film samples, showing well fitted result, as well as interesting relationships between parameter estimates and insulation morphology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reaction rates for reaction-diffusion kinetics on unstructured meshes

The reaction-diffusion master equation is a stochastic model often utilized in the study of biochemical reaction networks in living cells. It is applied when the spatial distribution of molecules is important to the dynamics of the system. A viable approach to resolve the complex geometry of cells accurately is to discretize space with an unstructured mesh. Diffusion is modeled as discrete jumps between nodes on the mesh, and the diffusion jump rates can be obtained through a discretization of the diffusion equation on the mesh. Reactions can occur when molecules occupy the same voxel. Here in this paper, we develop a method for computing accurate reaction rates between molecules occupying the same voxel in an unstructured mesh. For large voxels, these rates are known to be well approximated by the reaction rates derived by Collins and Kimball, but as the mesh is refined, no analytical expression for the rates exists. We reduce the problem of computing accurate reaction rates to a pure preprocessing step, depending only on the mesh and not on the model parameters, and we devise an efficient numerical scheme to estimate them to high accuracy. We show in several numerical examples that as we refine the mesh, the results obtained with the reaction-diffusion master equation approach those of a more fine-grained Smoluchowski particle-tracking model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global Data-Driven Determination of Baryon Transition Form Factors

Hadronic resonances emerge from strong interactions encoding the dynamics of quarks and gluons. The structure of these resonances can be probed by virtual photons parametrized in transition form factors. Here, in this study, twelve N* and Δ transition form factors at the pole are extracted from data with the center-of-mass energy from πN threshold to 1.8 GeV, and the photon virtuality 0 ≤ Q 2 /GeV 2 ≤ 8. For the first time, these results are determined from a simultaneous analysis of more than one state, i.e., ~10 5 π⁢N, η⁢N, and K⁢Λ electroproduction data. In addition, about 5×10 4 data in the hadronic sector as well as photoproduction serve as boundary conditions. For the Δ⁡(1232) and N⁡(1440) states our results are in qualitative agreement with previous studies, while the transition form factors at the poles of some higher excited states are estimated for the first time. Realistic uncertainties are determined by further exploring the parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Static and Dynamic Contributions to the Splitting of the Eu{sup 2+} Ground State in SrMoO{sub 4}

The EPR spectrum of an Eu{sup 2+} impurity center in a SrMoO{sub 4} single crystal in the temperature range T = 1.8, 111–300 K has been studied, and the temperature changes in the spin Hamiltonian parameters describing the EPR spectrum of odd europium isotopes have been determined. It is shown that small temperature changes in the diagonal parameters of the spin Hamiltonian (for odd Eu{sup 2+} isotopes) b{sub 2}{sup 0}(T) = b{sub 2}(F) + b{sub 2}(L) and P{sub 2}{sup 0}(T) = P{sub 2}(F) + P{sub 2}(L) are explained by the compensation of spin–phonon contributions b{sub 2}(F) and P{sub 2}(F) by the contributions of the lattice thermal expansion b{sub 2}(L) and P{sub 2}(L). The quantities b{sub 2}(L) and P{sub 2}(L) that are dependent on the static lattice parameters at a given temperature, are estimated in terms of the superposition Newman model. Then, the spin–phonon b{sub 2}(F) and P{sub 2}(F) contributions determined by the lattice ion vibrations are separated. An analysis shows that b{sub 2}{sup 0}(F) and P{sub 2}{sup 0}(F) > 0, b{sub 2}(L) and P{sub 2}(L) < 0, and the temperature behavior of the spin–phonon contribution is well described by G. Pfister’s model of local vibrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cutting out the middleman: calibrating and validating a dynamic vegetation model (ED2-PROSPECT5) using remotely sensed surface reflectance

Canopy radiative transfer is the primary mechanism by which models relate vegetation composition and state to the surface energy balance, which is important to light- and temperature-sensitive plant processes as well as understanding land–atmosphere feedbacks. In addition, certain parameters (e.g., specific leaf area, SLA) that have an outsized influence on vegetation model behavior can be constrained by observations of shortwave reflectance, thus reducing model predictive uncertainty. Importantly, calibrating against radiative transfer outputs allows models to directly use remote sensing reflectance products without relying on highly derived products (such as MODIS leaf area index) whose assumptions may be incompatible with the target vegetation model and whose uncertainties are usually not well quantified. Here, we created the EDR model by coupling the two-stream representation of canopy radiative transfer in the Ecosystem Demography model version 2 (ED2) with a leaf radiative transfer model (PROSPECT-5) and a simple soil reflectance model to predict full-range, high-spectral-resolution surface reflectance that is dependent on the underlying ED2 model state. We then calibrated this model against estimates of hemispherical reflectance (corrected for directional effects) from the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and survey data from 54 temperate forest plots in the northeastern United States. The calibration significantly reduced uncertainty in model parameters related to leaf biochemistry and morphology and canopy structure for five plant functional types.

54 ENVIRONMENTAL SCIENCES↗

Optimization of hot workability and microstructure control in a 12Cr-10Ni precipitation hardenable stainless steel: An approach using processing maps

Highlights: • Deformation processing map was developed for a 12Cr10Ni precipitation hardenable martensitic stainless steel and optimized hot working parameters were suggested to obtain defect free products. • The ‘stable’ domain of the processing maps were correlated to the dynamic recrystallization in the austenite phase and further confirmed by detailed microstructural investigations. • The DRx domain showed a refinement of the martensite packet size and of the prior transformation austenite grain size, essential to obtain the desired mechanical properties. - Abstract: Optimization of hot working parameters and corresponding microstructure control were explored in a 12Cr-10Ni precipitation hardenable martensitic stainless steel. From hot compression tests, 1073–1373 K (800°–1100 °C) and five (0.001–10 s{sup −1}) different strain rates, processing maps were developed, and various domains of microstructure evolution (in the austenite phase) were identified. The stable domain of the processing maps appeared to indicate dynamic recrystallization (DRx) in the austenite phase: as evident by the appearance of inflection point in the strain hardening rate versus flow stress plots. The suspected DRx domain was further confirmed by detailed microstructural investigations of the respective water-quenched specimens. The DRx domain showed a refinement of the martensite packet size and of the prior transformation austenite grain size. The later was shown with optical microstructure and reconstructed (modelled) microtexture data of the prior austenite grains. Additionally, the DRx domain showed a significant difference in activation energy as estimated from appropriate constitutive modelling.

36 MATERIALS SCIENCE↗

TRANSVERSE OSCILLATIONS OF LOOPS WITH CORONAL RAIN OBSERVED BY HINODE/SOLAR OPTICAL TELESCOPE

The condensations composing coronal rain, falling down along loop-like structures observed in cool chromospheric lines such as H{alpha} and Ca II H, have long been a spectacular phenomenon of the solar corona. However, considered a peculiar sporadic phenomenon, it has not received much attention. This picture is rapidly changing due to recent high-resolution observations with instruments such as the Hinode/Solar Optical Telescope (SOT), CRISP of the Swedish 1-m Solar Telescope, and the Solar Dynamics Observatory. Furthermore, numerical simulations have shown that coronal rain is the loss of thermal equilibrium of loops linked to footpoint heating. This result has highlighted the importance that coronal rain can play in the field of coronal heating. In this work, we further stress the importance of coronal rain by showing the role it can play in the understanding of the coronal magnetic field topology. We analyze Hinode/SOT observations in the Ca II H line of a loop in which coronal rain puts in evidence in-phase transverse oscillations of multiple strand-like structures. The periods, amplitudes, transverse velocities, and phase velocities are calculated, allowing an estimation of the energy flux of the wave and the coronal magnetic field inside the loop through means of coronal seismology. We discuss the possible interpretations of the wave as either standing or propagating torsional Alfven or fast kink waves. An estimate of the plasma beta parameter of the condensations indicates a condition that may allow the often observed separation and elongation processes of the condensations. We also show that the wave pressure from the transverse wave can be responsible for the observed low downward acceleration of coronal rain.

79 ASTRONOMY AND ASTROPHYSICS↗

CONDITIONS OF PASSAGE AND ENTRAPMENT OF TERRESTRIAL PLANETS IN SPIN-ORBIT RESONANCES

The dynamical evolution of terrestrial planets resembling Mercury in the vicinity of spin-orbit resonances is investigated using comprehensive harmonic expansions of the tidal torque taking into account the frequency-dependent quality factors and Love numbers. The torque equations are integrated numerically with a small step in time, including the oscillating triaxial torque components but neglecting the layered structure of the planet and assuming a zero obliquity. We find that a Mercury-like planet with a current value of orbital eccentricity (0.2056) is always captured in 3:2 resonance. The probability of capture in the higher 2:1 resonance is approximately 0.23. These results are confirmed by a semi-analytical estimation of capture probabilities as functions of eccentricity for both prograde and retrograde evolutions of spin rate. As follows from analysis of equilibrium torques, entrapment in 3:2 resonance is inevitable at eccentricities between 0.2 and 0.41. Considering the phase space parameters at the times of periastron, the range of spin rates and phase angles for which an immediate resonance passage is triggered is very narrow, and yet a planet like Mercury rarely fails to align itself into this state of unstable equilibrium before it traverses 2:1 resonance.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamic and Kinetic Modeling of Co-utilization of Glucose and Xylose for 2,3-BDO Production by Zymomonas mobilis

Prior engineering of the ethanologen Zymomonas mobilis has enabled it to metabolize xylose and to produce 2,3-butanediol (2,3-BDO) as a dominant fermentation product. When co-fermenting with xylose, glucose is preferentially utilized, even though xylose metabolism generates ATP more efficiently during 2,3-BDO production on a BDO-mol basis. To gain a deeper understanding of Z. mobilis metabolism, we first estimated the kinetic parameters of the glucose facilitator protein of Z. mobilis by fitting a kinetic uptake model, which shows that the maximum transport capacity of glucose is seven times higher than that of xylose, and glucose is six times more affinitive to the transporter than xylose. With these estimated kinetic parameters, we further compared the thermodynamic driving force and enzyme protein cost of glucose and xylose metabolism. It is found that, although 20% more ATP can be yielded stoichiometrically during xylose utilization, glucose metabolism is thermodynamically more favorable with 6% greater cumulative Gibbs free energy change, more economical with 37% less enzyme cost required at the initial stage and sustains the advantage of the thermodynamic driving force and protein cost through the fermentation process until glucose is exhausted. Glucose-6-phosphate dehydrogenase (g6pdh), glyceraldehyde-3-phosphate dehydrogenase (gapdh) and phosphoglycerate mutase (pgm) are identified as thermodynamic bottlenecks in glucose utilization pathway, as well as two more enzymes of xylose isomerase and ribulose-5-phosphate epimerase in xylose metabolism. Acetolactate synthase is found as potential engineering target for optimized protein cost supporting unit metabolic flux. Pathway analysis was then extended to the core stoichiometric matrix of Z. mobilis metabolism. Growth was simulated by dynamic flux balance analysis and the model was validated showing good agreement with experimental data. Dynamic FBA simulations suggest that a high agitation is preferable to increase 2,3-BDO productivity while a moderate agitation will benefit the 2,3-BDO titer. Taken together, this work provides thermodynamic and kinetic insights of Z. mobilis metabolism on dual substrates, and guidance of bioengineering efforts to increase hydrocarbon fuel production.

23-butanediol↗

Dynamical constraints on the dark matter distribution in the Milky Way

An accurate knowledge of the dark matter distribution in the Milky Way is of crucial importance for galaxy formation studies and current searches for particle dark matter. In this paper we set new dynamical constraints on the Galactic dark matter profile by comparing the observed rotation curve, updated with a comprehensive compilation of kinematic tracers, with that inferred from a wide range of observation-based morphologies of the bulge, disc and gas. The generalised Navarro-Frenk-White (NFW) and Einasto dark matter profiles are fitted to the data in order to determine the favoured ranges of local density, slope and scale radius. For a representative baryonic model, a typical local circular velocity v{sub 0}=230 km/s and a distance of the Sun to the Galactic centre R{sub 0}=8 kpc, we find a local dark matter density ρ{sub 0} = 0.420{sup +0.021}{sub −0.018} (2σ) ± 0.025 GeV/cm{sup 3} (ρ{sub 0} = 0.420{sup +0.019}{sub −0.021} (2σ) ± 0.026 GeV/cm{sup 3}) for NFW (Einasto), where the second error is an estimate of the systematic due to baryonic modelling. Apart from the Galactic parameters, the main sources of uncertainty inside and outside the solar circle are baryonic modelling and rotation curve measurements, respectively. Upcoming astronomical observations are expected to reduce all these uncertainties substantially over the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗