Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “eclipse”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

2M17091769+3127589: A Mass-transfer Binary with an Extreme Mass Ratio

We present the orbital solution of a peculiar double-lined spectroscopic and eclipsing binary system, 2M17091769+3127589. This solution was obtained by a simultaneous fit of both APOGEE radial velocities and TESS and ASAS-SN light curves to determine masses and radii. This system consists of an M=0.256{sub -0.006}{sup +0.010} M {sub ⊙}, R=3.961{sub -0.032}{sup +0.049} R {sub ⊙} red giant and a hotter M=1.518{sub -0.031}{sup +0.057} M {sub ⊙}, R=2.608{sub -0.321}{sup +0.034} R {sub ⊙} subgiant. Modeling with the MESA evolutionary codes indicates that the system likely formed 5.26 Gyrs ago, with a M = 1.2 M {sub ⊙} primary that is now the system’s red giant and a M = 1.11 M {sub ⊙} secondary that is now a more massive subgiant. Due to Roche-lobe overflow as the primary ascends the red giant branch, the more evolved “primary” (i.e., originally the more massive star of the pair) is now only one sixth as massive as the “secondary.” Such a difference between the initial and the current mass ratio is one of the most extreme detected so far. Evolutionary modeling suggests the system is still engaged in mass transfer, at a rate of M-dot ~10{sup -9} M {sub ⊙}/yr, and it provides an example of a less evolved precursor to some of the systems that consist of white dwarfs and blue stragglers.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Timing of five PALFA-discovered millisecond pulsars

We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the five pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 years. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ∼0.1 M{sub ⊙} companion and possible eclipses that last ∼10% of the orbit. The position of PSR J1957+2516 is also coincident with a near-infrared source. All five MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of γ-ray pulsations in a fold of Fermi Gamma-Ray Space Telescope data. These five MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

APPARENT NON-COEVALITY AMONG THE STARS IN UPPER SCORPIO: RESOLVING THE PROBLEM USING A MODEL OF MAGNETIC INHIBITION OF CONVECTION

Two eclipsing binaries in the USco association have recently yielded precise values of masses and radii for four low-mass members of the association. Standard evolution models would require these dM4.5–dM5 stars to have ages which are younger than those of more massive stars in the association by factors which appear (in extreme cases) to be as large as ∼3. Are the stars in the association therefore non-coeval? We suggest that the answer is no: by incorporating the effects of magnetic inhibition of convective onset, we show that the stars in USco can be restored to coevality provided the four low-mass member stars have vertical surface fields in the range 200–700 G. Fields of such magnitude have already been measured on the surface of certain solar-type stars in other young clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

THE AGE OF THE KIC 7177553 SYSTEM

KIC 7177553 is a quadruple system containing two binaries of orbital periods 16.5 and 18 days. All components have comparable masses and are slowly rotating with spectral types of ∼G2V. The longer period binary is eclipsing with component masses and radii M {sub 1} = 1.043 ± 0.014 M {sub ⊙}, R {sub 1} = 0.940 ± 0.005 R {sub ⊙} and M {sub 2} = 0.986 ± 0.015 M {sub ⊙}, R {sub 2} = 0.941 ± 0.005 R {sub ⊙}. The essentially equal radii measurements are inconsistent with the two stars being on the man sequence at the same age using standard nonmagnetic stellar evolution models. Instead a consistent scenario is found if the stars are in their pre-main-sequence phase of evolution and have an age of 32–36 Myr. We have also computed evolutionary models of magnetic stars, but we find that our nonmagnetic models fit the empirical radii and effective temperatures better than the magnetic models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Transient Transit Signature Associated with the Young Star RIK-210

We find transient transit-like dimming events within the K2 time series photometry of the young star RIK-210 in the Upper Scorpius OB association. These dimming events are variable in depth, duration, and morphology. High spatial resolution imaging revealed that the star is single and radial velocity monitoring indicated that the dimming events cannot be due to an eclipsing stellar or brown dwarf companion. Archival and follow-up photometry suggest the dimming events are transient in nature. The variable morphology of the dimming events suggests they are not due to a single spherical body. The ingress of each dimming event is always shallower than egress, as one would expect for an orbiting body with a leading tail. The dimming events are periodic and synchronous with the stellar rotation. However, we argue it is unlikely the dimming events could be attributed to anything on the stellar surface based on the observed depths and durations. Variable obscuration by a protoplanetary disk is unlikely on the basis that the star is not actively accreting and lacks the infrared excess associated with an inner disk. Rather, we explore the possibilities that the dimming events are due to magnetospheric clouds, a transiting protoplanet surrounded by circumplanetary dust and debris, eccentric orbiting bodies undergoing periodic tidal disruption, or an extended field of dust or debris near the corotation radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Pruning The ELM Survey: Characterizing Candidate Low-mass White Dwarfs through Photometric Variability

We assess the photometric variability of nine stars with spectroscopic T {sub eff} and log g values from the ELM Survey that locates them near the empirical extremely low-mass (ELM) white dwarf instability strip. We discover three new pulsating stars: SDSS J135512.34+195645.4, SDSS J173521.69+213440.6, and SDSS J213907.42+222708.9. However, these are among the few ELM Survey objects that do not show radial velocity (RV) variations that confirm the binary nature expected of helium-core white dwarfs. The dominant 4.31 hr pulsation in SDSS J135512.34+195645.4 far exceeds the theoretical cut-off for surface reflection in a white dwarf, and this target is likely a high-amplitude δ Scuti pulsator with an overestimated surface gravity. We estimate the probability to be less than 0.0008 that the lack of measured RV variations in four of eight other pulsating candidate ELM white dwarfs could be due to low orbital inclination. Two other targets exhibit variability as photometric binaries. Partial coverage of the 19.342 hr orbit of WD J030818.19+514011.5 reveals deep eclipses that imply a primary radius >0.4 R {sub ⊙}—too large to be consistent with an ELM white dwarf. The only object for which our time series photometry adds support to ELM white dwarf classification is SDSS J105435.78−212155.9, which has consistent signatures of Doppler beaming and ellipsoidal variations. We conclude that the ELM Survey contains multiple false positives from another stellar population at T {sub eff}≲9000 K, possibly related to the sdA stars recently reported from SDSS spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Activity and Rotation of Kepler-17

Magnetic activity on stars manifests itself in the form of dark spots on the stellar surface, which cause modulations of a few percent in the light curve of the star as it rotates. When a planet eclipses its host star, it might cross in front of one of these spots, creating a “bump” in the transit light curve. By modeling these spot signatures, it is possible to determine the physical properties of the spots such as size, temperature, and location. In turn, monitoring of the spots’ longitude provides estimates of the stellar rotation and differential rotation. This technique was applied to the star Kepler-17, a solar–type star orbited by a hot Jupiter. The model yields the following spot characteristics: average radius of 49 ± 10 Mm, temperatures of 5100 ± 300 K, and surface area coverage of 6 ± 4%. The rotation period at the transit latitude, −5{sup ∘}, occulted by the planet was found to be 11.92 ± 0.05 day, slightly smaller than the out-of-transit average period of 12.4 ± 0.1 day. Adopting a solar-like differential rotation, we estimated the differential rotation of Kepler-17 to be ΔΩ=0.041±0.005 rd day{sup −1}, which is close to the solar value of 0.050 rd day{sup −1}, and a relative differential rotation of ΔΩ/Ω=8.0±0.9%. Because Kepler-17 is much more active than our Sun, it appears that, for this star, larger rotation rate is more effective in the generation of magnetic fields than shear.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Spectrum of SS 433 in the H and K Bands

SS 433 is an X-ray binary and the source of sub-relativistic, precessing, baryonic jets. We present high-resolution spectrograms of SS 433 in the infrared H and K bands. The spectrum is dominated by hydrogen and helium emission lines. The precession phase of the emission lines from the jet continues to be described by a constant period, P{sub jet}=162.375 days. The limit on any secularly changing period is | P-dot |≲10{sup −5}. The He i λ2.0587 μm line has complex and variable P-Cygni absorption features produced by an inhomogeneous wind with a maximum outflow velocity near 900 km s{sup −1}. The He ii emission lines in the spectrum also arise in this wind. The higher members of the hydrogen Brackett lines show a double-peaked profile with symmetric wings extending more than ±1500 km s{sup −1} from the line center. The lines display radial velocity variations in phase with the radial velocity variation expected of the compact star, and they show a distortion during disk eclipse that we interpret as a rotational distortion. We fit the line profiles with a model in which the emission comes from the surface of a symmetric, Keplerian accretion disk around the compact object. The outer edge of the disk has velocities that vary from 110 to 190 km s{sup −1}. These comparatively low velocities place an important constraint on the mass of the compact star: its mass must be less than 2.2 M{sub ⊙} and is probably less than 1.6 M{sub ⊙}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Mass and p-factor of the Type II Cepheid OGLE-LMC-T2CEP-098 in a Binary System

We present the results of a study of the type II Cepheid (P {sub puls} = 4.974 days) in the eclipsing binary system OGLE-LMC-T2CEP-098 (P {sub orb} = 397.2 days). The Cepheid belongs to the peculiar W Vir group, for which the evolutionary status is virtually unknown. It is the first single-lined system with a pulsating component analyzed using the method developed by Pilecki et al. We show that the presence of a pulsator makes it possible to derive accurate physical parameters of the stars even if radial velocities can be measured for only one of the components. We have used four different methods to limit and estimate the physical parameters, eventually obtaining precise results by combining pulsation theory with the spectroscopic and photometric solutions. The Cepheid radius, mass, and temperature are 25.3±0.2 R{sub ⊙}, 1.51±0.09 M{sub ⊙}, and 5300±100 K, respectively, while its companion has a similar size (26.3 R{sub ⊙}), but is more massive (6.8 M{sub ⊙}) and hotter (9500 K). Our best estimate for the p-factor of the Cepheid is 1.30±0.03. The mass, position on the period–luminosity diagram, and pulsation amplitude indicate that the pulsating component is very similar to the Anomalous Cepheids, although it has a much longer period and is redder in color. The very unusual combination of the components suggest that the system has passed through a mass-transfer phase in its evolution. More complicated internal structure would then explain its peculiarity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Precision Determination of the Effect of Metallicity on Cepheid Absolute Magnitudes in VIJHK Bands from Magellanic Cloud Cepheids

Using high-quality observed period–luminosity ( P – L ) relations in both Magellanic Clouds in the VIJHK s bands and optical and near-infrared Wesenheit indices, we determine the effect of metallicity on Cepheid P – L relations by comparing the relative distance between the LMC and SMC as determined from the Cepheids to the difference in distance between the Clouds that has been derived with very high accuracy from late-type eclipsing binary systems. Within an uncertainty of 3%, which is dominated by the uncertainty on the mean difference in metallicity between the Cepheid populations in the LMC and SMC, we find metallicity effects smaller than 2% in all bands and in the Wesenheit indices, consistent with a zero metallicity effect. This result is valid for the metallicity range from −0.35 dex to −0.75 dex corresponding to the mean [Fe/H] values for classical Cepheids in the LMC and SMC, respectively. Yet most Cepheids in galaxies beyond the Local Group and located in the less crowded outer regions of these galaxies do fall into this metallicity regime, making our result important for applications to determine the distances to spiral galaxies well beyond the Local Group. Our result supports previous findings that indicated a very small metallicity effect on the near-infrared absolute magnitudes of classical Cepheids, and resolves the dispute about the size and sign of the metallicity effect in the optical spectral range. It also resolves one of the most pressing problems in the quest toward a measurement of the Hubble constant with an accuracy of 1% from the Cepheid–supernova Ia method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Thermal Phase Curve Offset on Tidally and Nontidally Locked Exoplanets: A Shallow Water Model

Using a shallow water model with time-dependent forcing, we show that the peak of an exoplanet thermal phase curve is, in general, offset from the secondary eclipse when the planet is rotating. That is, the planetary hot spot is offset from the point of maximal heating (the substellar point) and may lead or lag the forcing; the extent and sign of the offset are functions of both the rotation rate and orbital period of the planet. We also find that the system reaches a steady state in the reference frame of the moving forcing. The model is an extension of the well-studied Matsuno–Gill model into a full spherical geometry and with a planetary-scale translating forcing representing the insolation received on an exoplanet from a host star. The speed of the gravity waves in the model is shown to be a key metric in evaluating the phase curve offset. If the velocity of the substellar point (relative to the planet’s surface) exceeds that of the gravity waves, then the hot spot will lag the substellar point, as might be expected by consideration of forced gravity wave dynamics. However, when the substellar point is moving slower than the internal wave speed of the system, the hottest point may lead the passage of the forcing. We provide an interpretation of this result by consideration of the Rossby and Kelvin wave dynamics, as well as, in the very slowly rotating case, a one-dimensional model that yields an analytic solution. Finally, we consider the inverse problem of constraining planetary rotation rate from an observed phase curve.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Asteroseismology and Gaia: Testing Scaling Relations Using 2200 Kepler Stars with TGAS Parallaxes

We present a comparison of parallaxes and radii from asteroseismology and Gaia DR1 (TGAS) for 2200 Kepler stars spanning from the main sequence to the red-giant branch. We show that previously identified offsets between TGAS parallaxes and distances derived from asteroseismology and eclipsing binaries have likely been overestimated for parallaxes ≲5--10 mas (≈90%–98% of the TGAS sample). The observed differences in our sample can furthermore be partially compensated by adopting a hotter T{sub eff} scale (such as the infrared flux method) instead of spectroscopic temperatures for dwarfs and subgiants. Residual systematic differences are at the ≈2% level in parallax across three orders of magnitude. We use TGAS parallaxes to empirically demonstrate that asteroseismic radii are accurate to ≈5% or better for stars between ≈0.8--8 R{sub ⊙}. We find no significant offset for main-sequence (≲1.5 R{sub ⊙}) and low-luminosity RGB stars (≈3–8 R{sub ⊙}), but seismic radii appear to be systematically underestimated by ≈5% for subgiants (≈1.5–3 R{sub ⊙}). We find no systematic errors as a function of metallicity between [Fe/H]≈−0.8 to +0.4 dex, and show tentative evidence that corrections to the scaling relation for the large frequency separation (Δν) improve the agreement with TGAS for RGB stars. Finally, we demonstrate that beyond ≈3 kpc asteroseismology will provide more precise distances than end-of-mission Gaia data, highlighting the synergy and complementary nature of Gaia and asteroseismology for studying galactic stellar populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evolution of Spin, Orbital, and Superorbital Modulations of 4U 0114+650

We report a systematic analysis of the spin, orbital, and superorbital modulations of 4U 0114+650, a high-mass X-ray binary that consists of one of the slowest spinning neutron stars. Using the dynamic power spectrum, we found that the spin period varied dramatically and is anticorrelated with the long-term X-ray flux variation that can be observed using the Rossi X-ray Timing Explorer ASM, Swift BAT, and the Monitor of All-sky X-ray Image. The spin-up rate over the entire data set is consistent with previously reported values; however, the local spin-up rate is considerably higher. The corresponding local spin-up timescale is comparable to the local spin-up rate of OAO 1657−415, indicating that 4U 0114+650 could also have a transient disk. Moreover, the spin period evolution shows two ∼1000-day spin-down/random-walk epochs that appeared together with depressions of the superorbital modulation amplitude. This implies that the superorbital modulation was closely related to the presence of the accretion disk, which is not favored in the spin-down/random-walk epochs because the accretion is dominated by the direct wind accretion. The orbital period is stable during the entire time span; however, the orbital profile significantly changes with time. We found that the depth of the dip near the inferior conjunction of the companion is highly variable, which disfavors the eclipsing scenario. Moreover, the dip was less obvious during the spin-down/random-walk epochs, indicating its correlation with the accretion disk. Further monitoring in both X-ray and optical bands could reveal the establishment of the accretion disk in this system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Toward Space-like Photometric Precision from the Ground with Beam-shaping Diffusers

We demonstrate a path to hitherto unachievable differential photometric precisions from the ground, both in the optical and near-infrared (NIR), using custom-fabricated beam-shaping diffusers produced using specialized nanofabrication techniques. Such diffusers mold the focal plane image of a star into a broad and stable top-hat shape, minimizing photometric errors due to non-uniform pixel response, atmospheric seeing effects, imperfect guiding, and telescope-induced variable aberrations seen in defocusing. This PSF reshaping significantly increases the achievable dynamic range of our observations, increasing our observing efficiency and thus better averages over scintillation. Diffusers work in both collimated and converging beams. We present diffuser-assisted optical observations demonstrating 62{sub −16}{sup +26} ppm precision in 30 minute bins on a nearby bright star 16 Cygni A (V = 5.95) using the ARC 3.5 m telescope—within a factor of ∼2 of Kepler's photometric precision on the same star. We also show a transit of WASP-85-Ab (V = 11.2) and TRES-3b (V = 12.4), where the residuals bin down to 180{sub −41}{sup +66} ppm in 30 minute bins for WASP-85-Ab—a factor of ∼4 of the precision achieved by the K2 mission on this target—and to 101 ppm for TRES-3b. In the NIR, where diffusers may provide even more significant improvements over the current state of the art, our preliminary tests demonstrated 137{sub −36}{sup +64} ppm precision for a K {sub S} = 10.8 star on the 200 inch Hale Telescope. These photometric precisions match or surpass the expected photometric precisions of TESS for the same magnitude range. This technology is inexpensive, scalable, easily adaptable, and can have an important and immediate impact on the observations of transits and secondary eclipses of exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Model of the Pulsating Extremely Low-mass White Dwarf Precursor WASP 0247–25B

We present an analysis of the evolutionary and pulsation properties of the extremely low-mass white dwarf precursor (B) component of the double-lined eclipsing system WASP 0247−25. Given that the fundamental parameters of that star have been obtained previously at a unique level of precision, WASP 0247−25B represents the ideal case for testing evolutionary models of this newly found category of pulsators. Taking into account the known constraints on the mass, orbital period, effective temperature, surface gravity, and atmospheric composition, we present a model that is compatible with these constraints and show pulsation modes that have periods very close to the observed values. Importantly, these modes are predicted to be excited. Although the overall consistency remains perfectible, the observable properties of WASP 0247−25B are closely reproduced. A key ingredient of our binary evolutionary models is represented by rotational mixing as the main competitor against gravitational settling. Depending on assumptions made about the values of the degree index ℓ for the observed pulsation modes, we found three possible seismic solutions. We discuss two tests, rotational splitting and multicolor photometry, that should readily identify the modes and discriminate between these solutions. However, this will require improved temporal resolution and higher S/N observations, which are currently unavailable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Mind Your Ps and Qs: The Interrelation between Period (P) and Mass-ratio (Q) Distributions of Binary Stars

We compile observations of early-type binaries identified via spectroscopy, eclipses, long-baseline interferometry, adaptive optics, common proper motion, etc. Each observational technique is sensitive to companions across a narrow parameter space of orbital periods P and mass ratios q = M{sub comp}/M {sub 1}. After combining the samples from the various surveys and correcting for their respective selection effects, we find that the properties of companions to O-type and B-type main-sequence (MS) stars differ among three regimes. First, at short orbital periods P ≲ 20 days (separations a ≲ 0.4 au), the binaries have small eccentricities e ≲ 0.4, favor modest mass ratios 〈q〉≈0.5, and exhibit a small excess of twins q > 0.95. Second, the companion frequency peaks at intermediate periods log P (days) ≈ 3.5 (a ≈ 10 au), where the binaries have mass ratios weighted toward small values q ≈ 0.2–0.3 and follow a Maxwellian “thermal” eccentricity distribution. Finally, companions with long orbital periods log P (days) ≈ 5.5–7.5 (a ≈ 200–5000 au) are outer tertiary components in hierarchical triples and have a mass ratio distribution across q ≈ 0.1–1.0 that is nearly consistent with random pairings drawn from the initial mass function. We discuss these companion distributions and properties in the context of binary-star formation and evolution. We also reanalyze the binary statistics of solar-type MS primaries, taking into account that 30% ± 10% of single-lined spectroscopic binaries likely contain white dwarf companions instead of low-mass stellar secondaries. The mean frequency of stellar companions with q > 0.1 and log P (days) < 8.0 per primary increases from 0.50 ± 0.04 for solar-type MS primaries to 2.1 ± 0.3 for O-type MS primaries. We fit joint probability density functions f(M{sub 1},q,P,e)≠f(M{sub 1})f(q)f(P)f(e) to the corrected distributions, which can be incorporated into binary population synthesis studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Wayne—A Simulator for HST WFC3 IR Grism Spectroscopy

Wayne is an algorithm that simulates Hubble Space Telescope Wide Field Camera 3 (WFC3) grism spectroscopic frames, including sources of noise and systematics. It can simulate both staring and spatial scan modes, and observations such as the transit and the eclipse of an exoplanet. Unlike many other instrument simulators, the focus of Wayne is on creating frames with realistic systematics in order to test the effectiveness of different data analysis methods in a variety of different scenarios. This approach is critical for method validation and optimizing observing strategies. In this paper we describe the implementation of Wayne for WFC3 in the near-infrared channel with the G102 and G141 grisms. We compare the simulations to real data obtained for the exoplanet HD 209458b, to verify the accuracy of the simulation. The software is now available as open source at https://github.com/ucl-exoplanets/wayne.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Searching for Rapid Orbital Decay of WASP-18b

The WASP-18 system, with its massive and extremely close-in planet, WASP-18b ( M{sub p} = 10.3 M{sub J}, a = 0.02 au, P = 22.6 hr), is one of the best-known exoplanet laboratories to directly measure Q ′, the modified tidal quality factor and proxy for efficiency of tidal dissipation, of the host star. Previous analysis predicted a rapid orbital decay of the planet toward its host star that should be measurable on the timescale of a few years, if the star is as dissipative as is inferred from the circularization of close-in solar-type binary stars. We have compiled published transit and secondary eclipse timing (as observed by WASP, TRAPPIST, and Spitzer ) with more recent unpublished light curves (as observed by TRAPPIST and Hubble Space Telescope ) with coverage spanning nine years. We find no signature of a rapid decay. We conclude that the absence of rapid orbital decay most likely derives from Q ′ being larger than was inferred from solar-type stars and find that Q ′ ≥ 1 × 10{sup 6}, at 95% confidence; this supports previous work suggesting that F stars, with their convective cores and thin convective envelopes, are significantly less tidally dissipative than solar-type stars, with radiative cores and large convective envelopes.

79 ASTRONOMY AND ASTROPHYSICS↗