Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maximum Likelihood”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Deducing neutron star equation of state from telescope spectra with machine-learning-derived likelihoods

The interiors of neutron stars reach densities and temperatures beyond the limits of terrestrial experiments, providing vital laboratories for probing nuclear physics. While the star's interior is not directly observable, its pressure and density determine the star's macroscopic structure which affects the spectra observed in telescopes. The relationship between the observations and the internal state is complex and partially intractable, presenting difficulties for inference. Previous work has focused on the regression from stellar spectra of parameters describing the internal state. We demonstrate a calculation of the full likelihood of the internal state parameters given observations, accomplished by replacing intractable elements with machine learning models trained on samples of simulated stars. Our machine-learning-derived likelihood allows us to perform maximum a posteriori estimation of the parameters of interest, as well as full scans. We demonstrate the technique by inferring stellar mass and radius from an individual stellar spectrum, as well as equation of state parameters from a set of spectra. Our results are more precise than pure regression models, reducing the width of the parameter residuals by 11.8% in the most realistic scenario. The neural networks will be released as a tool for fast simulation of neutron star properties and observed spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗
DOE OSTI2025

A strategy for automated core design to increase economic viability and minimize fuel fragmentation, relocation, and dispersal susceptibility in high-burnup cores

The nuclear industry aims to increase the cycle length of pressurized water reactors from 18 to 24 months to increase power plant capacity factors and economic viability. These cycle length extensions will inherently require fuel rods to exceed the current peak rod average burnup limit of 62 GWd/MTU. A chief concern of operating beyond the current burnup limit is the fuel fragmentation, relocation, and dispersal (FFRD) phenomenon in which pulverized fuel fragments can axially relocate and escape through a burst in the cladding formed during a loss-of-coolant accident. In this work, we demonstrate an approach for automating core design employing an optimization tool based on a penalty-free, parallel simulated annealing algorithm to produce pressurized water reactor core designs with two different optimization objectives. The two objectives were to produce core designs with (1) mitigated FFRD susceptibility while achieving 24-month cycle lengths (2) maximum cycle length with no regard for the likelihood of FFRD. Batch size was considered in tandem with both cases to maximize economic viability. The PARCS nodal model was the primary reactor physics tool used in the optimizations and used nuclear cross sections calculated with 2D Polaris lattice physics models. Reactor performance and safety characteristics of the optimized cores were verified using high-fidelity Virtual Environment for Reactor Applications models. The core designs produced by the optimization tool are compared with each other and to a high-burnup core design produced and analyzed in previous works to highlight the fuel management strategies that may enhance high-burnup reactor safety and economic viability. The optimized cores satisfied their respective objective functions, producing a maximum cycle length of 720 effective full-power days in one core design and one that may reduce FFRD susceptibility by up to 50% based on the first-order approximation to FFRD risk formulated in this work. The optimized cores met most constraints but exceeded the hot channel factor limit, especially in FFRD cases where fresh fuel carried more power. Furthermore, this highlights the need for future lattice-level optimizations and broader assembly options.

Cycle length↗
DOE OSTI2020

A Climatology and Extreme Value Analysis of Large Hail in China

In this paper, hail size records are analyzed at 2254 stations in China and a hail size climatology is developed based on gridded hail observations for the period 1960–2015. It is found that the annual percentiles of hail size records changed sharply and national-wide after 1980, therefore two periods, 1960–79 and 1980–2015, are studied. There are some similarities between the two periods in terms of the characteristics of hail size such as the spatial distribution patterns of mean annual maximum hail size and occurrence week of annual maximum hail size. The 1980–2015 period had higher observation density than the 1960–79 period, but showed smaller mean annual maximum hail size, especially in northern China. In the majority of grid boxes, the annual maximum hail size experienced a decreasing trend during the 1980–2015 period. A Gumbel extreme value model is fitted to each grid box to estimate the return periods of maximum hail size. The scale and location parameter of the fitted Gumbel distributions are higher in eastern China than in western China, thereby reflecting a greater likelihood of large hail in eastern China. In southern China, the maximum hail size exceeds 127 mm for a 10-yr return period, whereas in northern China maximum hail size exceeds this threshold for a 50-yr return period. Furthermore, the Gumbel model is found to potentially underestimate the maximum hail size for certain return periods, but provides a more informed picture of the spatial distribution of extreme hail size and the regional differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Enhancing Gaussian Process Surrogates for Optimization and Posterior Approximation via Random Exploration

This paper proposes novel noise-free Bayesian optimization strategies that rely on a random exploration step to enhance the accuracy of Gaussian process surrogate models. The new algorithms retain the ease of implementation of the classical GP-UCB algorithm, but the additional random exploration step accelerates their convergence, nearly achieving the optimal convergence rate. Furthermore, to facilitate Bayesian inference with intractable likelihoods, we propose to utilize optimization iterates for maximum a posteriori estimation to build a Gaussian process surrogate model for the unnormalized log-posterior density. We provide bounds for the Hellinger distance between the true and the approximate posterior distributions in terms of the number of design points. We demonstrate the effectiveness of our Bayesian optimization algorithms in nonconvex benchmark objective functions, in a machine learning hyperparameter tuning problem, and in a black-box engineering design problem. The effectiveness of our posterior approximation approach is demonstrated in two Bayesian inference problems for parameters of dynamical systems.

Bayesian inference↗
DOE OSTI2025

An iterative CMB lensing estimator minimizing instrumental noise bias

Noise maps from cosmic microwave background (CMB) experiments are generally statistically anisotropic, due to scanning strategies, atmospheric conditions, or instrumental effects. Any mismodeling of this complex noise can bias the reconstruction of the lensing potential and the measurement of the lensing power spectrum from the observed CMB maps. We introduce a new CMB lensing estimator based on the maximum (MAP) reconstruction that is minimally sensitive to these instrumental noise biases. By modifying the likelihood to rely exclusively on correlations between CMB map splits with independent noise realizations, we minimize autocorrelations that contribute to biases. In the regime of many independent splits, this maximum closely approximates the optimal MAP reconstruction of the lensing potential. In simulations, we demonstrate that this method is able to determine lensing observables that are immune to any noise mismodeling with a negligible cost in signal-to-noise ratio. Our estimator enables unbiased and nearly optimal lensing reconstruction for next-generation CMB surveys.

Legrand, Louis [University of Cambridge (United Ki↗
DOE OSTI2022

BayoTech Risk and Modeling Support of NM Gas Site

The BayoTech hydrogen generation system has been evaluated in terms of safety considerations at the NM Gas site. The consequence of a leak in different components in the system was evaluated in terms of plume dispersion and overpressure. Additionally, the likelihood of a leak scenario for different hydrogen components was identified. The worst-case plume dispersion cases, full-bore leaks, resulted in relatively large plumes. However, these cases were noted to be far less likely than the partial break cases that were evaluated. The partial break cases resulted in nearly negligible plume lengths. Similarly, the overpressure analysis of the full-bore break scenarios resulted in much larger overpressures than the partial break cases (which resulted in negligible overpressure at the lot line). There were several cases evaluated in the analysis that represented leak scenarios from both hydrogen and natural gas sources. Generally, the natural gas leak scenarios resulted in a smaller horizontal impact than that of hydrogen leaks. The worst-case consequence from a hydrogen leak resulted from the compressors, storage pods, or dispensing system. To consider the safety features that may isolate the leak, the consequence was evaluated at different times after the leak event to show the reduction of pressure. After 2 seconds, the plume dispersion from this event is contained within the perimeter of the site. The worst-case consequences show that the plume may disperse to adjacent facilities and to the street. When considering both likelihood and consequence, the risk may be considered low because the maximum frequency of a full-bore leak from any component within the hydrogen compound is 8.2 E-5/yr. This means that a full-bore leak is expected to occur less than once every 10,000 years. The risk can be further reduced by implementing mitigative countermeasures, such as CMU walls along the sides of the equipment compound. This would reduce the overall consequence of the worst-case dispersion scenarios (horizontal impact of plume). In terms of siting and safety analysis, the NFPA 2 code was used to provide a high-level evaluation of the current site plan. The most limiting equipment in terms of set-back distance are the compressors/storage units because of the high-pressure hydrogen. The site layout was evaluated for an acceptable location for the compression/storage unit based on NFPA 2 set-back distances. It is important to note that the NFPA 2 set-back distances consider both likelihood and consequence. This is important because the worst-case results evaluated herein also represent the least likely leak scenario. Other site-specific considerations were evaluated, including the parking shade structure with photovoltaic cells and refueling vehicles. These issues were dispositioned and determined to not present a safety risk.

08 HYDROGEN↗
DOE OSTI2022

Pervasive alterations to snow-dominated ecosystem functions under climate change

Climate change projections consistently demonstrate that warming temperatures and dwindling seasonal snowpack will elicit cascading effects on ecosystem function and water resource availability. Despite this consensus, little is known about potential changes in the variability of ecohydrological conditions, which is also required to inform climate change adaptation and mitigation strategies. Considering potential changes in ecohydrological variability is critical to evaluating the emergence of trends, assessing the likelihood of extreme events such as floods and droughts, and identifying when tipping points may be reached that fundamentally alter ecohydrological function. Using a single-model Large Ensemble with sophisticated terrestrial ecosystem representation, we characterize projected changes in the mean state and variability of ecohydrological processes in historically snow-dominated regions of the Northern Hemisphere. Widespread snowpack reductions, earlier snowmelt timing, longer growing seasons, drier soils, and increased fire risk are projected for this century under a high-emissions scenario. In addition to these changes in the mean state, increased variability in winter snowmelt will increase growing-season water deficits and increase the stochasticity of runoff. Thus, with warming, declining snowpack loses its dependable buffering capacity so that runoff quantity and timing more closely reflect the episodic characteristics of precipitation. This results in a declining predictability of annual runoff from maximum snow water equivalent, which has critical implications for ecosystem stress and water resource management. So our results suggest that there is a strong likelihood of pervasive alterations to ecohydrological function that may be expected with climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temporal persistence of postfire flood hazards under present and future climate conditions in southern Arizona, USA

Changes to soil hydraulic properties that reduce infiltration capacity following fire can increase flash flood risks. These risks are exacerbated by rainfall intensification associated with a warming climate. However, the potential effects of climate-change-driven rainfall intensification on postfire floods remain largely unexplored. Using rainfall and runoff observations from a 49.4 km 2 watershed in southern Arizona, USA, and a hydrologic model (KINEROS2), we examined the temporal evolution following a historic fire of three crucial hydrologic parameters: soil saturated hydraulic conductivity (K sp ), net capillary drive (G p ), and hydraulic roughness (n c ). We explored how the effect of fire on these parameters may influence peak flow under future climate scenarios derived from CMIP6, specifically the medium emissions scenario (SSP2-4.5) and high emissions scenario (SSP5-8.5). Results demonstrate an increase in K sp from 11 mm h −1 in the first postfire year to 60 mm h −1 in postfire year 3. G p similarly increased from 19 mm in the first postfire year to 30 mm in the third, while nc was relatively constant. The highest simulated Q p occurred in the first postfire year. Under the SSP2-4.5 scenario, the likelihood of a 100-year flood is projected to be twice as large by the middle of the century relative to its historical magnitude. Simulations further indicate that the maximum expected discharge associated with a postfire flood, as derived from historical data, could be triggered by a 10-year rainstorm under the SSP5-8.5 scenario by the late century. Simulations also demonstrate that rainfall intensification will lead to greater persistence of elevated flood hazards following fire by the late century under both the SSP2-4.5 and the SSP5-8.5 scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fragmentation of photovoltaic backsheets after accelerated weathering exposure

The backsheet of a photovoltaic (PV) module functions primarily to provide safe operation through electrical insulation and protect the interior from environmental exposure. For most applications, backsheets consist of opaque, but highly reflective, films that protect inner layers from degradation from ultraviolet (UV) light. Backsheets can fail if they delaminate or form cracks under strain. In a PV module, the maximum strains that can be applied are only on the order of a few percent. The greatest strains are likely to be highly localized around topographical features. Therefore, to assess the likelihood of a crack being initiated in a backsheet, we apply small amplitude strains to backsheet films under microscopic examination using a laser scanning confocal microscope. Furthermore, this allows us to monitor both the depth and density of cracks as well as the onset strain. Here we perform these measurements after exposure to UV light.

14 SOLAR ENERGY↗
DOE OSTI2023

Simulation Based Inference with Domain Adaptation for Strong Gravitational Lensing

Simulation based inference leverages machine learning to carry out Bayesian inference for systems with intractable likelihoods. However, transitioning a network trained on simulated data to real data runs the risk of encountering domain shift, leading to performance losses. We attempt to implement domain adaptation into the sbi neural posterior estimation framework using the Maximum Mean Discrepancy as an additional network loss. We test two network architectures and use masked autoregressive flow for density estimation. We test the network on a set of 400,000 simulated strong gravitational lensing images generated using deeplenstronomy. The source domain is defined as low noise whereas the target domain has a noise profile sampled from experimentally derived DES survey conditions. We find that, while DA does lead to performance improvements, they are marginal at ~6% less inference error. We also find a similar marginal improvement in uncertainty calibration at around 8%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Simulation based inference with domain adaptation for strong gravitational lensing

Simulation based inference leverages machine learning to carry out Bayesian inference in systems with intractable likelihoods. However, transitioning a network trained on simulated data to real data runs the risk of encountering domain shift, leading to performance losses. We attempt to implement domain adaptation into the sbi neural posterior estimation framework using the Maximum Mean Discrepancy as an additional network loss, using masked autoregressive fow (MAF) as our density estimator. We test the network on a set of 400,000 simulated strong gravitational lensing images generated using deeplenstronomy. The source domain is defined as low noise whereas the target domain has a noise profile sampled from experimentally derived DES survey conditions. We find that SBI appears robust against small changes in the data with similar performance on source and target. Moreover, while DA does lead to performance improvements, they are marginal at 6% less inference error.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Characterizing electrical demand and load diversity of low-power water and space heating appliances in US homes

Home renovation and remodeling projects can involve costly and time consuming electrical infrastructure upgrades at the household level. From the grid perspective they also lead to costly replacement of local infrastructure, such as transformers, and can add stress to the grid at peak times. The emergence of innovative, power-efficient household appliances offers a way to minimize these problems. These appliances are designed for lower power consumption, simplifying installation through standard plug-in connections, reducing the need for new electric circuits/panels/service, and minimizing the peak power demand for the home. Key examples include low-power heat pump water heaters (HPWHs) and cold climate window heat pumps that operate on standard 120V outlets. To assess the real-world impact of these solutions, we compiled and analyzed power metering data from several US field studies. This data provides insights into the effects on peak power demand of selecting lower-power appliances. Our analysis focuses on several key metrics, including the maximum power demand of individual appliances, their operational runtime, continuous operation and load diversity. While individual low-power 120V space and water heating appliances offer significant peak demand reductions compared to 240V heat pump or resistance alternatives, their longer runtimes might increase the likelihood of operation during whole-dwelling peak events, albeit at lower power levels.

Less, Brennan↗
DOE OSTI2020

Planning green infrastructure placement based on projected precipitation data

Continued urbanization has led to tremendous changes on the landscape. These changes have exacerbated the effects of extreme climatic events such as flooding because of constrained water infiltration and increased surface flow. Typical runoff control measures involve sophisticated gray infrastructure that guide excess surface flow into storage and disposal sites. In a dynamic climate system, these measures are not sustainable since they cannot be easily modified to accommodate large volumes of runoff. Green Infrastructure (GI) is an adaptable technique that can be used to minimize runoff, in addition to offering an array of additional benefits (urban heat regulation, aesthetics, improved air quality etc.). Strategic placement of GI is key to achieving maximum utility. While physical site characteristics play a major role in determining suitable GI placement sites, knowledge of future precipitation patterns is crucial to ensure successful flood mitigation. Here, suitable GI sites within the city of Knoxville, Tennessee, were determined based on potential impact of an extreme flood event as indicated by site characteristics. Then, the relative potential likelihood of a flood event was determined based on projected precipitation data and knowledge of existing flood zones. By combining potential impact with likelihood information, low, medium, and high priority GI implementation sites were established. Results indicate that high priority sites are in the central parts of the city with priority decreasing outward. The GI prioritization scheme presented here, offers valuable guidance to city planners and policy makers who wish to exploit the GI approach for flood mitigation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Investigation on Tuning Power–Frequency Droop for Improved Grid–Forming Inverter and Synchronous Generator Transient Load Sharing

Poor transient load sharing has been observed during the parallel operation of synchronous generators (SGs) and droop-controlled grid-forming (GFM) inverter-based resources (IBRs) in islanded microgrids. This can result in overcurrent protections, causing the inverters to trip, which leads to the collapse of voltage and frequency. In this article, small-signal analysis supported by electromagnetic transients (EMTs) simulation is performed on a detailed model of a microgrid containing IBRs and an SG to identify underdamped modes and their sensitivities to the power–frequency (P-F) droop parameter. Further EMT simulation of a more complex microgrid model with multiple GFM IBRs and SGs is performed to determine whether trends identified from the modal analysis extend to cases where the number of SGs in operation may vary. The objective is to investigate the effectiveness of tuning this parameter to reduce the likelihood of overcurrents in response to load and generation disturbances. The work uncovers both fast and slow modes of concern. Primarily, the findings show that increasing the P-F droop value improves the damping of the slower mode while degrading the damping of the faster mode, leading to an upper limit on the damping improvement of the GFM IBR output current. In EMT simulations of the more complex microgrid, the greatest damping improvement still exhibits significant maximum overshoot. Furthermore, the optimal value of the P-F droop parameter is sensitive to the number of SGs in operation. The conclusion is that tuning the P-F droop parameter is not an effective means to avoid the overcurrent tripping of GFM IBRs during response to large load or generation disturbances when operating in parallel to SGs in islanded microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Alleviating prior dependencies for DESI DR1 clustering fits through reparameterization

Bayesian analyses of the full-shape clustering of Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) exhibit prior-volume projection effects, whereby weakly constrained nuisance parameters of the Effective Field Theory of Large Scale Structure (EFTofLSS) shift marginalized cosmological posteriors away from the posterior maximum. We reanalyze DESI DR1 power spectrum multipoles using two complementary mitigation strategies: (i) nonlinear orthogonalization to decorrelate nuisance and cosmological parameter priors, and (ii) a fully reparameterization-invariant Jeffreys prior over all EFTofLSS coefficients, evaluated on-the-fly via closed-form Jacobians. Including data from DESI, Big-Bang Nuclesynthesis and a constraint on $n_{\mathrm{s}}$, baseline priors lead to multi-$σ$ projection in the Hubble parameter $H_{0}$ and dark energy equation of state parameters $w_{0}$ and $w_{a}$; the Jeffreys prior successfully recenters these posteriors to enclose the maximum a posteriori estimate within the 68% credible regions, demonstrating clear mitigation of projection effects for these late-time expansion parameters. A hybrid Jeffreys+baseline-Gaussian configuration controls residual over-broad tails in the physical cold dark matter density $ω_{\mathrm{c}}$ while preserving the volume correction, and is our favoured approach. We compare the credible intervals derived using our methodology to those obtained using Halo Occupation Distribution (HOD)-informed priors and to confidence intervals derived using frequentist profile likelihood analyses, finding agreement in both central values and degeneracy directions in the $w_{0}$--$w_{a}$ plane. This demonstrates that, once projection effects are properly controlled, we can make robust inferences about the late-time cosmological expansion independent of the statistical framework adopted.

Bonici, M. [Waterloo U.; Perimeter Inst. Theor. Ph↗
DOE OSTI2025

Compactly‐Supported Nonstationary Kernels for Computing Exact Gaussian Processes on Big Data

The Gaussian process (GP) is a widely used method for analyzing large-scale data sets, including spatio-temporal measurements of nonlinear processes that are now commonplace in the environmental sciences. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility, and exact methods for inference that prevent application to data sets with more than about 10,000 points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space–time prediction based on more than 1 million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Gaussian processes↗
DOE OSTI2021

Dark Energy Survey year 3 results: covariance modelling and its impact on parameter estimation and quality of fit

ABSTRACT We describe and test the fiducial covariance matrix model for the combined two-point function analysis of the Dark Energy Survey Year 3 (DES-Y3) data set. Using a variety of new ansatzes for covariance modelling and testing, we validate the assumptions and approximations of this model. These include the assumption of Gaussian likelihood, the trispectrum contribution to the covariance, the impact of evaluating the model at a wrong set of parameters, the impact of masking and survey geometry, deviations from Poissonian shot noise, galaxy weighting schemes, and other sub-dominant effects. We find that our covariance model is robust and that its approximations have little impact on goodness of fit and parameter estimation. The largest impact on best-fitting figure-of-merit arises from the so-called fsky approximation for dealing with finite survey area, which on average increases the χ2 between maximum posterior model and measurement by $3.7{{\ \rm per\ cent}}$ (Δχ2 ≈ 18.9). Standard methods to go beyond this approximation fail for DES-Y3, but we derive an approximate scheme to deal with these features. For parameter estimation, our ignorance of the exact parameters at which to evaluate our covariance model causes the dominant effect. We find that it increases the scatter of maximum posterior values for Ωm and σ8 by about $3{{\ \rm per\ cent}}$ and for the dark energy equation-of-state parameter by about $5{{\ \rm per\ cent}}$.

79 ASTRONOMY AND ASTROPHYSICS↗