Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radii”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

PARET/ANL v7.7 User Guide

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of many reactors both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL v7.7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Continuous-Energy Verification of MCNP Calculations Using One-Group Spherical and Slab Criticality Benchmarks

This work presents a continuous-energy Monte Carlo verification study of one group spherical and slab criticality benchmarks using the MCNP ®1 code. Classical tabulations and newly generated benchmark solutions obtained by direct numerical evaluation by the author are considered. The benchmarks span weakly to strongly multiplying regimes and provide analytically defined critical radii as functions of a single parameter, c .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Anomalous Hall effects in chiral superconductors

We report theoretical results for the electronic contribution to thermal and electrical transport for chiral superconductors belonging to even or odd-parity E 1 and E 2 representations of the tetragonal and hexagonal point groups. Chiral superconductors exhibit novel properties that depend on the topology of the order parameter and Fermi surface, and—as we highlight—the structure of the impurity potential. An anomalous thermal Hall effect is predicted and shown to be sensitive to the winding number, ν, of the chiral order parameter via Andreev scattering that transfers angular momentum from the chiral condensate to excitations that scatter off the random potential. For heat transport in a chiral superconductor with isotropic impurity scattering, i.e., point-like impurities, a transverse heat current is obtained for ν = ± 1, but vanishes for | ν | > 1. This is not a universal result. For finite-size impurities with radii of order or greater than the Fermi wavelength, R ≥ ℏ / p f, the thermal Hall conductivity is finite for chiral order with | ν | ≥ 2, and determined by a specific Fermi-surface average of the differential cross-section for electron-impurity scattering. Our results also provide quantitative formulae for analyzing and interpreting thermal transport measurements for superconductors predicted to exhibit broken time-reversal and mirror symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Proton electromagnetic generalized polarizabilities

Electromagnetic polarizabilities are fundamental properties of the proton that characterize its response to an external electromagnetic (EM) field. The generalization of the EM polarizabilities to non-zero four-momentum transfer opens up a powerful path to study the internal structure of the proton. They map out the spatial distribution of the polarization densities in the proton, provide access to key dynamical mechanisms that contribute to the electric and magnetic polarizability effects, and allow for the determination of fundamental characteristics of the system, such as the electric and magnetic polarizability radii. This article reviews our knowledge about proton EM generalized polarizabilities (GPs). An introduction is given to the basic concepts and the theoretical framework, which is then followed by a discussion that emphasizes the recent developments and findings of the virtual Compton scattering (VCS) experiments and future perspectives on the topic.

hadrons↗
DOE OSTI2023

Effects of the evolving early Moon and Earth magnetospheres

Recently it has been identified that our Moon had an extensive magnetosphere for several hundred million years soon after it was formed when the Moon was within 20 Earth Radii (R E ) from the Earth. Some aspects of the interaction between the early Earth-Moon magnetospheres are investigated by mapping the interconnected field lines between the Earth and the Moon and investigating how the early lunar magnetosphere affects the magnetospheric dynamics within the coupled magnetospheres over time. So long as the magnetosphere of the Moon remains strong as it moves away from the Earth in the antialigned dipole configuration, the extent of the Earth’s open field lines decreases. As a result, at times it significantly changes the structure of the field-aligned current system, pushing the polar cusp significantly northward, and forcing magnetotail reconnection sites into the deeper tail region. In addition, the combined magnetospheres of the Earth and the Moon greatly extend the number of closed field lines enabling a much larger plasmasphere to exist and connecting the lunar polar cap with closed field lines to the Earth. That configuration supports the transfer of plasma between the Earth and the Moon potentially creating a time capsule of the evolution of volatiles with depth. This paper only touches on the evolution of the early Earth and Moon magnetospheres, which has been a largely neglected space physics problem and has great potential for complex follow-on studies using more advanced tools and due to the expected new lunar data coming in the next decade through the Artemis Program.

Astronomy & Astrophysics↗
DOE OSTI2023

Influences of Cloud Microphysics on the Components of Solar Irradiance in the WRF-Solar Model

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Crystal Chemistry and Structural Variations for Zircon Samples from Various Localities

This study investigates the variations of structural parameters and chemistry of a partially metamict and seven detrital zircon samples from different localities using single-crystal X-ray diffraction, synchrotron high-resolution powder X-ray diffraction, and electron-probe micro-analysis techniques. The unit-cell parameters for the eight zircon samples vary linearly with increasing unit-cell volume, V. A zircon sample from the Canadian Arctic Islands has the smallest unit-cell parameters, bond distances, ideal stoichiometric composition, unaffected by α-radiation damage, and is chemically pure. A zircon sample from Jemaa, Nigeria has the largest unit-cell parameters because of the effect of α-radiation doses received over a long time (2384 Ma). All the samples show good correlations between Zr and Si apfu (atom per formula unit) versus unit-cell volume, V. The α-radiation doses in the samples are lower than ~3.5 × 1015 α-decay events/mg. Substitutions of other cations at the Zr and Si sites control the variations of the structural parameters. Relatively large unit-cell parameters and bond distances occur because the Zr site accommodates other cations that have larger ionic radii than the Zr atom. Geological age increases the radiation doses in zircon and it is related to V.

58 GEOSCIENCES↗
DOE OSTI2025

Tuning Solvation Dynamics of Electrolytes at Their Eutectic Point Through Halide Identity

Deep eutectic solvents (DESs) are regarded as highly promising solvent systems for redox flow batteries. DESs, composed of choline halides (ChX, X = F − , Cl − , Br − , I − ) and ethylene glycol (EG), exhibit distinct physicochemical properties at their eutectic points, including halide-dependent phase behavior, viscosity, polarity, conductivity, and solvation dynamics. In this study, we investigate the effects of the halide identity on the solvation properties of ChX:EG mixtures at varying mol % of ChX salt content. The solvatochromic polarity based on ET(30) measurements indicates higher polarity for larger halides (I − > Br − ) than for smaller halides (Cl − > F − ), which exhibit larger compensating solvation shells. The ionic conductivity follows the trend of the solvent fluidity (the inverse of the viscosity), namely ChCl > ChBr > ChI > ChF, influenced by the ion mobility and solvodynamic radii. Measurements of the liquidus temperatures (T L ) reveal that the system with ChCl exhibits the deepest eutectic point (at ~20 mol % ChCl), while ChBr and ChI have shallower minima at ~10 mol % ChBr and ~3 mol % ChI, respectively. ChF does not display a eutectic transition but instead appears to readily supercool at salt concentrations above 30 mol % ChF. Consistent with the phase transition measurements, femtosecond transient absorption spectroscopy shows that in the ChCl system, the solvation dynamics become faster with an increasing salt concentration up to ~16.67 mol %, after which the dynamics slow down with further increases in the salt content. The ChF-based system exhibits similar behavior, though with slower dynamics. In contrast, the solvation dynamics of the systems containing ChBr and ChI monotonously slow down with an increasing salt concentration, in agreement with the phase transition measurements, which show that the eutectic points occur at low salt concentrations. These measurements suggest that the solvent composition and, in particular, the identity of the halide anion play a significant role in the solvation behavior of these ethylene-glycol-based DESs, offering a foundation for tuning the DES properties for specific applications.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Progress in Constraining Nuclear Symmetry Energy Using Neutron Star Observables Since GW170817

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λ l , for several high-order multipoles (l = 5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M ⊙ NS, R 1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, R M . We find that there is an optimal value of M for every M such that the uncertainty in the estimate of R M is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of 0.369{sub −0.021}{sup +0.031} M{sub J}, 0.243{sub −0.030}{sup +0.022} M{sub J}, 0.353{sub −0.027}{sup +0.038} M{sub J}, and 0.1254±0.0039 M{sub J}, respectively, and radii of 1.117±0.014 R{sub J}, 0.800±0.015 R{sub J}, 0.765±0.013 R{sub J}, and 0.7224±0.0032 R{sub J}, respectively. The planets orbit close to their host stars with orbital periods of 3.9228 days, 3.1317 days, 4.1480 days, and 7.3279 days, respectively. The hosts are main-sequence K dwarfs with masses of 0.674{sub −0.012}{sup +0.016} M{sub ⊙}, 0.7279 ±0.0066 M{sub ⊙}, 0.7133±0.0075 M{sub ⊙}, and 0.7311±0.0028, and with V-band magnitudes of V=14.829±0.010, 14.35±0.11, 14.998±0.040 and 12.469±0.010. The super-Neptune HATS-72b (a.k.a. WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP, and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS, and Magellan/PFS. We also incorporate TESS observations of the warm Saturn–hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A, and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved [[Formula]] neighbor in Gaia DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of 0.22 M{sub ⊙} and a current projected physical separation of ∼1400 au.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Probing UV-sensitive Pathways for CN and HCN Formation in Protoplanetary Disks with the Hubble Space Telescope

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and submillimeter wavelengths, we present an analysis of new and archival Hubble Space Telescope (HST), Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure Lyα and far-UV (FUV) continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g., HCN, N{sub 2}). Semi-forbidden C ii] λ2325 and UV-fluorescent H{sub 2} emission were also measured to constrain inner disk populations of C{sup +} and vibrationally excited H{sub 2}. We find a significant positive correlation between 14 μm HCN emission and fluxes from the FUV continuum and C ii] λ2325, consistent with model predictions requiring N{sub 2} photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between submillimeter CN emission and both C ii] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and Lyα or UV-H{sub 2} emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with the James Webb Space Telescope will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Ringed Pole-on Outflow from DO Tauri Revealed by ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations including the CO(2−1) line and 1.3 mm continuum emission from the surroundings of the young stellar object DO Tauri. The ALMA CO molecular data show three different series of rings at different radial velocities. These rings have radii around 220 and 800 au. We make individual fits to the rings and note that their centers are aligned with DO Tauri and its optical high-velocity jet. In addition, we notice that the velocity of these structures increases with the separation from the young star. We discuss the data under the hypothesis that the rings represent velocity cuts through three outflowing shells that are possibly driven by a wide-angle wind, dragging the environment material along a direction close to the line of sight (i = 19°). We estimate the dynamical ages, the mass, the momentum, and the energy of each individual outflow shell and those of the whole outflow. The results are in agreement with those found in outflows from Class II sources. We make a rough estimate for the size of the jet/wind launching region, which needs to be of ≲15 au. We report the physical characteristics of DO Tauri’s disk continuum emission (almost face-on and with a projected major axis in the north–south direction) and its velocity gradient orientation (north–south), indicative of disk rotation for a 1–2 M {sub ⊙} central star. Finally, we show a Hubble Space Telescope [S ii] image of the optical jet and report a measurement of its orientation in the plane of the sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

APOGEE Net: Improving the Derived Spectral Parameters for Young Stars through Deep Learning

Machine learning allows for efficient extraction of physical properties from stellar spectra that have been obtained by large surveys. The viability of machine-learning approaches has been demonstrated for spectra covering a variety of wavelengths and spectral resolutions, but most often for main-sequence (MS) or evolved stars, where reliable synthetic spectra provide labels and data for training. Spectral models of young stellar objects (YSOs) and low-mass MS stars are less well-matched to their empirical counterparts, however, posing barriers to previous approaches to classify spectra of such stars. In this work, we generate labels for YSOs and low-mass MS stars through their photometry. We then use these labels to train a deep convolutional neural network to predict logg, T {sub eff}, and Fe/H for stars with Apache Point Observatory Galactic Evolution Experiment (APOGEE) spectra in the DR14 data set. This “APOGEE Net” has produced reliable predictions of logg for YSOs, with uncertainties of within 0.1 dex and a good agreement with the structure indicated by pre-MS evolutionary tracks, and it correlates well with independently derived stellar radii. These values will be useful for studying pre-MS stellar populations to accurately diagnose membership and ages.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fundamental Parameters of ∼30,000 M dwarfs in LAMOST DR1 Using Data-driven Spectral Modeling

M dwarfs are the most common type of star in the Galaxy, and because of their small size are favored targets for searches of Earth-sized transiting exoplanets. Current and upcoming all-sky spectroscopic surveys, such as the Large Sky Area Multi Fiber Spectroscopic Telescope (LAMOST), offer an opportunity to systematically determine physical properties of many more M dwarfs than has been previously possible. Here, we present new effective temperatures, radii, masses, and luminosities for 29,678 M dwarfs with spectral types M0–M6 in the first data release (DR1) of LAMOST. We derived these parameters from the supervised machine-learning code, The Cannon, trained with 1388 M dwarfs in the Transiting Exoplanet Survey Satellite Cool Dwarf Catalog that were also present in LAMOST with high signal-to-noise ratio (>250) spectra. Our validation tests show that the output parameter uncertainties are strongly correlated with the signal-to-noise of the LAMOST spectra, and we achieve typical uncertainties of 110 K in T{sub eff} (∼3%), 0.065 R{sub ⊙} (∼14%) in radius, 0.054 M{sub ⊙} (∼12%) in mass, and 0.012 L{sub ⊙} (∼20%) in luminosity. The model presented here can be rapidly applied to future LAMOST data releases, significantly extending the samples of well-characterized M dwarfs across the sky using new and exclusively data-based modeling methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Kepler Data Analysis: Non-Gaussian Noise and Fourier Gaussian Process Analysis of Stellar Variability

We develop a statistical analysis model of Kepler stellar flux data in the presence of planet transits, non-Gaussian noise, and stellar variability. We first develop a model for the Kepler noise probability distribution in the presence of outliers, which make the noise probability distribution non-Gaussian. We develop a signal likelihood analysis based on this probability distribution, in which we model the signal as a sum of the star variability and planetary transits. We argue that these components need to be modeled together if optimal signal is to be extracted from the data. For the stellar variability model we develop an optimal Gaussian process analysis using a Fourier-based Wiener filter approach, where the power spectrum is non-parametric and learned from the data. We develop high dimensional optimization of the objective function, where we jointly optimize all the model parameters, including thousands of star variability modes, and planet transit parameters. We apply the method to Kepler-90 data and show that it gives a better match to the stellar variability than the existing methods, and robustly handles noise outliers. As a consequence, the planet radii have a higher value than what the existing methods give, including splines and celerite.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Keck/NIRC2 L’-band Imaging of Jovian-mass Accreting Protoplanets around PDS 70

We present L’-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semimajor axis of 20{sub −4}{sup +3} au and PDS 70 c to have a semimajor axis of 34{sub −6}{sup +12} au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2–3 R {sub Jup}. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 M {sub Jup} and a mean mass accretion rate between 3 × 10{sup −7} and 8 × 10{sup −7} M {sub Jup}/yr. For PDS 70 c, we computed a mass between 1 and 3 M {sub Jup} and mean mass accretion rate between 1 × 10{sup −7} and 5 × 10{sup −7} M {sub Jup}/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets’ SEDs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

KR Persei, a Mid-F Eclipsing Binary with a One-day Period

KR Per is a partially eclipsing binary with an orbital period of 0.9960798 days, very close to one sidereal day, making it difficult to obtain extensive phase coverage in a reasonable amount of time. We used the Wilson–Devinney program to determine its orbital elements and stellar absolute dimensions from recently acquired radial velocities and differential BVRI observations that were supplemented with previous differential UBV measurements and published times of minima. The two components are each F5 V stars with masses of 1.466±0.015 M{sub ⊙} and 1.458±0.015 M{sub ⊙}. The radii are 1.855±0.021 R{sub ⊙} and 1.824±0.022 R{sub ⊙}. The orbital period of the eclipsing system is variable and more times of minima observations are needed. A comparison with evolutionary tracks indicates that the system has an age of 2.1 ± 0.1 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗