Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 307 records · Page 17

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗

The Atacama Cosmology Telescope: DR4 maps and cosmological parameters

In this paper, we present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013–2016 at 98 and 150 GHz. The maps cover more than 17,000 deg 2 , the deepest 600 deg 2 with noise levels below 10µK-arcmin. We use the power spectrum derived from almost 6,000 deg 2 of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, H 0 . By combining ACT data with large-scale information from WMAP we measure H 0 = 67.6±1.1 km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find H 0 = 67.9 ± 1.5 km/s/Mpc). The ΛCDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1σ; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with ΛCDM predictions to within 1.5–2.2σ. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis.

60 APPLIED LIFE SCIENCES↗

Noble gas constraints on spent fuel irradiation histories

Fission gas isotopic compositions are sensitive to a variety of reactor operating parameters that include the neutron flux, neutron fluence at different neutron energies, and operating temperature. Measurements of fission gas isotopic compositions thus have potential to constrain reactor simulations for nuclear forensics and safeguards applications. In this paper, we present Kr and Xe isotope measurements from a suite of samples obtained from locations that span the axial length a fuel pin with a well-characterized irradiation history and compare these data to spatially resolved reactor simulations. Here, we observed positive correlations between fluence sensitive isotopic ratios and burnup and between a flux sensitive ratio and power, although some discrepancies are observed between the measured data and model predictions. These differences may be due to simplifications in the model and/or inaccuracies in the cross sections. A much broader measurement to model comparison is required to better understand the discrepancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mechanistic Insight into Tunable Spin Relaxation in Two-Dimensional Type-II Ligand-Perovskite Heterostructures

Two-dimensional (2D) metal-halide perovskites with spin-dependent optical properties hold great promise for spintronic and quantum applications. However, their spin lifetimes, especially for n = 1 2D perovskites, are typically limited to subpicosecond time scales due to rapid spin relaxation driven by strong spin−orbit coupling (SOC), electron−hole exchange interactions, and phonon-mediated scattering. Here, we demonstrate that type-II ligand-perovskite heterostructures overcome these constraints by reducing electron−hole wave function overlap and exciton binding energy. Compared to the type-I 2D perovskite (PEA) 2 PbI 4 with a spin lifetime of 0.29 ps at room temperature, our engineered type-II systems achieve substantially extended spin lifetimes, ∼6.37 ps for (4Tm) 2 PbI 4 and ∼18.47 ps for (4TCNm) 2 PbI 4 . In both materials, spatial charge separation across the perovskite−ligand interface mitigates the Bir−Aronov− Pikus (BAP) mechanism. Temperature- and fluence-dependent measurements reveal Elliott−Yafet (EY)-dominated spin relaxation in (4Tm) 2 PbI 4 , consistent with the observation of coherent phonon oscillation, whereas (4TCNm) 2 PbI 4 exhibits D’yakonov−Perel (DP)-dominated spin relaxation, with weaker phonon coupling further suppressing the EY relaxation, enabling spin lifetimes up to ∼126.81 ps at 5 K. Our findings establish a structural design framework for tailoring spin dynamics in 2D perovskites, offering a promising strategy to engineering spin and optoelectronic properties via rational ligand engineering.

Excitons↗

Quasar clustering at redshift 6

Large-scale surveys over the last years have revealed about 300 quasi-stellar objects (QSOs) at redshifts above 6. Follow-up observations have identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, and high CIV-MgII Velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some BHs likely have high masses at birth and grow efficiently. Our aim is to quantify quasar pairs at high redshift for a large sample of objects. This provides a new key constraint on a combination of parameters related to the origin and assembly for the most massive BHs: formation efficiency and clustering, growth efficiency, and the relative contribution of BH mergers. We observed 116 spectroscopically confirmed QSOs around redshift 6 with the simultaneous seven-channel imager Gamma-ray Burst Optical/Near-infrared Detector in order to search for companions. Applying colour-colour cuts identical to those which led to the spectroscopically confirmed QSOs, we performed Le PHARE fits to the 26 best QSO pair candidates, and obtained spectroscopic observations for 11 of them. We do not find any QSO pair with a companion brighter than M 1450 (AB) < –26 mag within our 0.1–3.3 h –1 cMpc search radius, in contrast to the serendipitous findings in the redshift range 4–5. However, a small fraction of such pairs at this luminosity and redshift is consistent with indications from present-day cosmological-scale galaxy evolution models. In turn, the incidence of L- and T-type brown dwarfs, which occupy a similar colour space to z ~ 6 QSOs, is higher than expected, by a factor of 5 and 20, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

Deciphering phase evolution in complex metal oxide thin films via high-throughput materials synthesis and characterization

Discovery of structure-property relationships in thin film alloys of complex metal oxides enabled by high-throughput materials synthesis and characterization facilities is demonstrated here with a case-study. In this study, thin films of binary transition metal oxides (Ti–Zn) are prepared by pulsed laser deposition with continuously varying Ti:Zn ratio, creating combinatorial samples for exploration of the properties of this material family. The atomic structure and electronic properties are probed by spatially resolved techniques including x-ray absorption near edge structures (XANES) and x-ray fluorescence (XRF) at the Ti and Zn K-edge, x-ray diffraction, and spectroscopic ellipsometry. The observed properties as a function of Ti:Zn ratio are resolved into mixtures of five distinguishable phases by deploying multivariate curve resolution analysis on the XANES spectral series, under constraints set by results from the other characterization techniques. First-principles computations based on density function theory connect the observed properties of each distinct phase with structural and spectral characteristics of crystalline polymorphs of Ti–Zn oxide. Continuous tuning of the optical absorption edge as a function of Ti:Zn ratio, including the unusual observation of negative optical bowing, exemplifies a functional property of the film correlated to the phase evolution.

36 MATERIALS SCIENCE↗

Multi-objective Optimization Paradigm for Toroidal Inductors with Spatially Tuned Permeability

Spatially tuning core permeability of an electromagnetic device enables superior performance. A permeability profile can be heuristically selected to improve the flux distribution in a device with a given geometry, but in order to fully leverage the capacity of spatial dependent permeability engineering, the geometry and the permeability should be optimized simultaneously. The work herein presented sets forth a multi-physics design optimization paradigm that includes the permeability profile tuning in the context of an inductor design. This approach enables the determination of Pareto optimal fronts consisting of a set of optimal solutions against competing objectives (e.g. mass and loss) under imposed constraints. To this end, analytical solutions of the heat transfer and electromagnetic formulations are derived for toroidal inductors. The software implemented in Matlab 2018b is available online as an attachment to this paper.

42 ENGINEERING↗

T -matrix analysis of static Wilson line correlators from lattice QCD at finite temperature

Here, we utilize a previously constructed thermodynamic T -matrix approach to the quark-gluon plasma (QGP) to derive constaints on its input by using results on Wilson line correlators (WLCs) of a static quark-antiquark pair from 2 + 1-flavor lattice-QCD (lQCD) computations with realistic pion mass. The self-consistent T -matrix results, which include previous constraints from the lQCD equation of state in the light-parton sector, can describe the lQCD data for WLCs fairly well once refinements of its driving kernel are applied. In particular, the input potential requires less screening than what has been inferred from previous T-matrix analyses. Pertinent predictions for the spectral and transport properties of the QGP are discussed, including the spatial diffusion coefficient for heavy quarks; the latter turns out to have a rather weak temperature dependence, in approximate agreement with recent unquenched lQCD results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LensWatch. I. Resolved HST Observations and Constraints on the Strongly Lensed Type Ia Supernova 2022qmx (“SN Zwicky”)

Supernovae (SNe) that have been multiply imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multicycle "LensWatch (www.lenswatch.org)" program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, and F814W) but unresolved with WFC3/IR F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens-model-predicted time delays (≲1 day), and delays estimated with the single epoch of Hubble Space Telescope colors (≲3.5 days), including the uncertainty from chromatic microlensing (~1–1.5 days). Our lens models converge to an Einstein radius of ${\theta }_{{\rm{E}}}=({0.168}_{-0.005}^{+0.009})^{\prime\prime} $, the smallest yet seen in a lensed SN system. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by $\sim {1.7}_{-0.6}^{+0.8}$ mag and $\sim {0.9}_{-0.6}^{+0.8}$ mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models.

79 ASTRONOMY AND ASTROPHYSICS↗

Joint Survey Processing. I. Compact Oddballs in the COSMOS Field—Low-luminosity Quasars at z > 6?

The faint-end slope of the quasar luminosity function at z ~ 6 and its implication on the role of quasars in reionizing the intergalactic medium at early times has been an outstanding problem for some time. The identification of faint high-redshift quasars with luminosities of <10 44.5 erg s –1 is challenging. They are rare (few per square degree), and the separation of these unresolved quasars from late-type stars and compact star-forming galaxies is difficult from ground-based observations alone. In addition, source confusion becomes significant at >25 mag, with ~30% of sources having their flux contaminated by foreground objects when the seeing resolution is ~0''.7. We mitigate these issues by performing a pixel-level joint processing of ground and space-based data from Subaru/Hyper-SuprimeCam (HSC) and Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). We create a deconfused catalog over the 1.64 deg 2 of the COSMOS field, after accounting for spatial varying point-spread functions and astrometric differences between the two data sets. We identify twelve low-luminosity (M UV ~ –21 mag) z > 6 quasar candidates through (i) their red color measured between ACS/F814W and HSC/i band and (ii) their compactness in the space-based data. Nondetections of our candidates in Hubble DASH data argues against contamination from late-type stars. Our constraints on the faint end of the quasar luminosity function at z ~ 6.4 suggest a negligibly small contribution to reionization compared to the star-forming galaxy population. The confirmation of our candidates and the evolution of number density with redshift could provide better insights into how supermassive galaxies grew in the first billion years of cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗

Testing cosmology with double source lensing

Double source lensing provides a dimensionless ratio of distance ratios, a "remote viewing" of cosmology through distances relative to the gravitational lens, beyond the observer. We use this to test the cosmological framework, particularly with respect to spatial curvature and the distance duality relation. We derive a consistency equation for constant spatial curvature, allowing not only the investigation of flat vs curved but of the Friedmann-Lemaître-Robertson-Walker framework itself. For distance duality, we demonstrate that the evolution of the lens mass profile slope must be controlled to ≳ 5 times tighter fractional precision than a claimed distance duality violation. Furthermore, using LensPop forecasts of double source lensing systems in Euclid and LSST surveys we also explore constraints on dark energy equation of state parameters and any evolution of the lens mass profile slope.

79 ASTRONOMY AND ASTROPHYSICS↗

An empirical calibration of the serpentine-water oxygen isotope fractionation at T = 25–100°C

The ability to precisely constrain the physical–chemical conditions of serpentinization, such as temperature and fluid source, is limited by the accuracy of calibrations for oxygen isotope fractionation between serpentine and water – i.e., 1000 lnα( serpentine-water ) – which disagree by up to 20‰ when extrapolated to T<200°C. In this study, we present a new empirical calibration of 1000 lnα( serpentine-water ) aiming to improve applications of oxygen isotope geochemistry to very low-T serpentinization (T<100°C). We used the high-spatial resolution capabilities of Secondary Ion Mass Spectrometry (SIMS) to analyze oxygen isotope ratios in mineral pairs of calcite + serpentine, quartz + serpentine and talc + serpentine co-crystallized at scales≤50μm in six serpentinite samples from the Samail ophiolite (Oman). SIMS analysis shows that the mineral pairs are relatively homogenous in oxygen isotope ratios with variability in δ 18 O values ≤2‰ (2s). Clumped isotope thermometry and petrological constraints indicate crystallization temperatures from 25 to 100°C for the investigated samples. These independent constraints on temperature allowed us to derive 1000 lnα( serpentine-water ) by combining Δ 18 O( quartz-serpentine ), Δ 18 O( talc-serpentine ) and Δ 18 O( calcite-serpentine ) measured by SIMS in the investigated samples with Δ 18 O( quartz-water ), Δ 18 O( talc-water ) and Δ 18 O( calcite-water ) calibrations available in the literature. Our empirical calibration of 1000 lnα( serpentine-water ) = 1.04±0.20 (2SE) × 10 6 /T 2 (T in K), from T = 25 to 100°C, is within uncertainty of former high-temperature empirical calibrations extrapolated to T<100°C and matches experimental calibrations when extrapolated to T=250°C. In conclusion, the new serpentine-water calibration enables more accurate reconstructions of fluid-rock interactions occurring during low-temperature serpentinization processes in various tectonic settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optically multiplexed neutron time-of-flight technique for inertial confinement fusion

Neutron time-of-flight (nTOF) detectors are crucial in diagnosing the performance of inertial confinement fusion (ICF) experiments, which implode targets of deuterium–tritium fuel to achieve thermonuclear conditions. These detectors utilize the fusion neutron energy spectrum to extract key measurements, including the hotspot ion temperature and fuel areal density. Previous work [Danly et al., Rev. Sci. Instrum. 94, 043502 (2023)] has demonstrated adding 1D spatial resolution to an nTOF-like detector using a neutron aperture and streak camera to measure the ion temperature profile of an ICF implosion. By contrast, the study presented herein explores modifying the 1D detector to use a fast photomultiplier tube (PMT) to validate the design of a 2D spatially resolved instrument based on reconstruction from 1D profiles. The modification would collect time-of-flight traces from separate scintillators in an imaging array with one PMT using optical fibers of varying lengths to time-multiplex the signals. This technique has been demonstrated in ride-along experiments on the OMEGA laser with 20 fiber-coupled scintillator channels connected to a Photek PMT210. Results provide constraints on the fiber lengths and PMT gating requirements to promote pulse fidelity throughout all channels. Calibration of the detector to fixed nTOFs can provide a preliminary estimate of the instrument response function (IRF), although measurement of the IRF is currently under way. These results suggest that nTOF signals can potentially be time-multiplexed with fibers so long as the design is strategic to mitigate signal-to-noise reduction, modal dispersion, and charge build-up in the PMT, which has implications beyond ion temperature imaging.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

11-th order of accuracy for numerical solution of 3-D Poisson equation with irregular interfaces on unfitted Cartesian meshes

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

97 MATHEMATICS AND COMPUTING↗

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗

CO2ROpt (Feedstock Allocation Optimization Model) [SWR-22-13]

Reducing greenhouse gas emissions remains a critical challenge for the transportation and chemical sectors. The electrochemical reduction of carbon dioxide is a potential pathway for production of fuels and chemicals that uses atmospheric carbon dioxide as a feedstock. Here we present an analysis of the potential for carbon dioxide captured from point sources and via direct air capture to be utilized in electrochemical reduction under different market scenarios. Using a numerical optimization framework and national data sets, scenarios for nationwide production of ethylene, formate, and carbon monoxide are evaluated at scales comparable to current production of these molecules in the United States. We show that developing a network for production of these products at scale requires capture and utilization of significant portions of the carbon dioxide that is currently emitted from large stationary point sources. Because carbon dioxide point sources are spatially and compositionally variable, their use for carbon dioxide reduction depends on electricity prices, capture cost, and location. If the power sector in the United States is decarbonized, carbon dioxide supply decreases significantly, increasing the importance of utilizing other carbon dioxide streams, and increasing the likelihood that direct air capture plays a role in supplying carbon dioxide feedstocks. These findings provide a first insight into the feedstock and energy constraints for carbon dioxide reduction.

Badgett, Alex↗

Active Galactic Nuclei Feedback in SDSS-IV MaNGA: AGNs Have Suppressed Central Star Formation Rates

Despite the importance of feedback from active galactic nuclei (AGNs) in models of galaxy evolution, observational constraints on the influence of AGN feedback on star formation remain weak. To this end, we have compared the star formation trends of 279 low-redshift AGN galaxies with 558 inactive control galaxies using integral field unit spectroscopy from the Sloan Digital Sky Survey-IV Mapping Nearby Galaxies at Apache Point Observatory survey. With a Gaussian-process-based methodology, we reconstruct nonparametric star formation histories in spatially resolved spaxels covering the face of each galaxy. Based on the galaxy-wide star formation rates (SFRs) alone, we find no obvious signatures of AGN feedback. However, the AGN galaxies have significantly suppressed central (kiloparsec-scale) SFRs, lying up to a factor of 2 below those of the control galaxies, providing direct observational evidence of AGN feedback suppressing star formation. The suppression of central SFRs in the AGN galaxies began in the central regions ~6 Gyr ago (redshift z ~ 0.7), taking place over a few gigayears. A small subset of the AGN galaxies were rapidly driven to quiescence shortly before being observed (in the last 500 Myr), potentially indicating instances of AGN-driven feedback. More frequently, however, star formation continues in the AGN galaxies, with suppression primarily in the central regions. This is suggestive of a picture in which integrated (gigayear-timescale) AGN feedback can significantly affect central star formation, but may be inefficient in driving galaxy-wide quenching in low-redshift galaxies, instead leaving them in the green valley.

79 ASTRONOMY AND ASTROPHYSICS↗