Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formation and Evolution of Compact Objects”

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.

48 records · Page 3

On the Anticorrelation Between High Accretion Luminosity and Radio Jet Ejection in Gro J1655-40 and Other Objects

A model is described in which radio jet formation in accreting objects is suppressed by processes which occur when the accretion rate approaches the Eddington limit. This was motivated by GRO J1655-40 and other objects which show an anticorrelation between high luminosity and the onset of a radio jet. The jet production mechanism employed is the Blandford-Payne MHD acceleration process, seeded by an e'e- pair wind. Observations indicate that all key ingredients of this mechanism are, or should be, present in these sources. Observed jet velocities and total powers are consistent with theoretical and numerical predictions of this model. The primary jet suppression mechanism proposed is the Papaloizou-Pringle compressible shear instability, which should disrupt the jet-producing region of the disk when the accretion rate approaches about one-third Eddington. The turn-on of the jet in GRO J 1 655 - 40 is consistent with this estimate. When super-Eddington, the disk should also drive an optically thick, subrelativistic wind, which may be a secondary jet suppression mechanism. The possible presence of such a wind is seen in the early spectral evolution of GRO J1655-40 and in the broad absorption lines of certain quasi-stellar objects (QSOs). Important tests of the model would be independent measurements of the compact object masses in these sources, a comparison of normal and broad absorption line (BAL) QSO X-ray spectra to see if the latter objects are significantly cooler, and a low-frequency search around BAL QSOs to see if at least some have fossil radio sources.

Meier, D.↗

High-mass X-ray binary populations. 1: Galactic modeling

Modern stellar evolutionary tracks are used to calculate the evolution of a very large number of massive binary star systems (M(sub tot) greater than or = 15 solar mass) which cover a wide range of total masses, mass ratios, and starting separations. Each binary is evolved accounting for mass and angular momentum loss through the supernova of the primary to the X-ray binary phase. Using the observed rate of star formation in our Galaxy and the properties of massive binaries, we calculate the expected high-mass X-ray binary (HMXRB) population in the Galaxy. We test various massive binary evolutionary scenarios by comparing the resulting HMXRB predictions with the X-ray observations. A major goal of this study is the determination of the fraction of matter lost from the system during the Roche lobe overflow phase. Curiously, we find that the total numbers of observable HMXRBs are nearly independent of this assumed mass-loss fraction, with any of the values tested here giving acceptable agreement between predicted and observed numbers. However, comparison of the period distribution of our HMXRB models with the observed period distribution does reveal a distinction among the various models. As a result of this comparison, we conclude that approximately 70% of the overflow matter is lost from a massive binary system during mass transfer in the Roche lobe overflow phase. We compare models constructed assuming that all X-ray emission is due to accretion onto the compact object from the donor star's wind with models that incorporate a simplified disk accretion scheme. By comparing the results of these models with observations, we conclude that the formation of disks in HMXRBs must be relatively common. We also calculate the rate of formation of double degenerate binaries, high velocity detached compact objects, and Thorne-Zytkow objects.

Dalton, William W.↗

Search for Subsolar Mass Ultracompact Binaries in Advanced LIGO’s Second Observing Run

Here, we present a search for subsolar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational-wave candidates consistent with subsolar mass ultracompact binaries with at least one component between $0.2 M_{⊙}-1.0 M_{⊙}$. We use the null result to constrain the binary merger rate of ($0.2 M_{⊙}, 0.2 M_{⊙}$) binaries to be less than $3.7×10^{5} Gpc^{-3} yr^{-1}$ and the binary merger rate of ($1.0 M_{⊙}, 1.0 M_{⊙}$) binaries to be less than $5.2×10^{3} Gpc^{-3} yr^{-1}$. Subsolar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with subsolar masses. Assuming a particular primordial black hole (PBH) formation model, we constrain a population of merging $0.2 M_{⊙}$ black holes to account for less than 16% of the dark matter density and a population of merging $1.0 M_{⊙}$ black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict subsolar mass binaries.

79 ASTRONOMY AND ASTROPHYSICS↗

Optical observations of NGC 2915: A nearby blue compact dwarf galaxy

This paper presents B and R band Charge Coupled Device (CCD) images and medium resolution spectroscopy of NGC 2915, a relatively isolated BCD (blue compact dwarf) galaxy at a distance of approximately 5 Mpc. NGC 2915 contains two stellar populations: a high surface brightness blue core population and a red diffuse population. The core population contains all of the H II, and numerous embedded objects. It is the locus of current high mass star formation. The brightest embedded objects are likely to be young ionizing clusters, while many of the fainter objects are likely to be individual supergiant stars with masses up to approximately 25 solar mass, or blends of a few such stars. Curious aligned structures on the SE side of the galaxy are seen and their nature discussed. The spectrum of the core is dominated by bright narrow emission lines like that of a high excitation and low metallicity (less than half solar) H II region. The continuum is flat, with Balmer and Ca II features seen in absorption. The velocity of the Ca II features suggest contamination by galactic interstellar absorption. There is a significant velocity gradient in the spectra, probably indicative of rotation. Outside of its core, NGC 2915 resembles a dE (dwarf elliptical) galaxy, in that it has an exponential surface brightness profile, is red ((B-R)(sub 0) = 1.65), and has a low extrapolated central surface brightness (B(0)(sub c) = 22.44). NGC 2915's properties are compared with other BCDs, concentrating on two morphologically similar BCDs that are near enough to resolve into stars: NGC 1705 and NGC 5253. It is noted that the presence of winds in BCDs invalidates closed box chemical evolution models and the remaining constraints on star formation duration are relatively weak. Some BCDs, including NGC 2915, may be able to maintain their present star formation rate for Gyr time scales. This suggests that the overall evolution of these BCDs may be much slower than the approximately 10 Myr burst time scales commonly quoted. However, shortly after the formation of a massive (10(exp 6) solar mass) cluster a BCD will have all the properties of strong starburst galaxy).

Meurer, G. R.↗

Formation and Evolution of Binary Systems Containing Collapsed Stars

This research includes theoretical studies of the formation and evolution of five types of interacting binary systems. Our main focus has been on developing a number of comprehensive population synthesis codes to study the following types of binary systems: (i) cataclysmic variables (#3, #8, #12, #15), (ii) low- and intermediate-mass X-ray binaries (#13, #20, #21), (iii) high-mass X-ray binaries (#14, #17, #22), (iv) recycled binary millisecond pulsars in globular clusters (#5, #10, #ll), and (v) planetary nebulae which form in interacting binaries (#6, #9). The numbers in parentheses refer to papers published or in preparation that are listed in this paper. These codes take a new unified approach to population synthesis studies. The first step involves a Monte Carlo selection of the primordial binaries, including the constituent masses, and orbital separations and eccentricities. Next, a variety of analytic methods are used to evolve the primary star to the point where either a dynamical episode of mass transfer to the secondary occurs (the common envelope phase), or the system evolves down an alternate path. If the residual core of the primary is greater than 2.5 solar mass, it will evolve to Fe core collapse and the production of a neutron star and a supernova explosion. In the case of systems involving neutron stars, a kick velocity is chosen randomly from an appropriate distribution and added to the orbital dynamics which determine the state of the binary system after the supernova explosion. In the third step, all binaries which commence stable mass transfer from the donor star (the original secondary in the binary system) to the compact object, are followed with a detailed binary evolution code. Finally, we include all the relevant dynamics of the binary system. For example, in the case of LMXBs, the binary system, with its recoil velocity from the supernova explosion, is followed in time through its path in the Galactic potential. For our globular cluster study, detailed three-body calculations are used to follow the interactions of the binary with passing field stars.

Rappaport, Saul↗

The electrospray: Fundamentals and combustion applications

Liquid fuel dispersion in practical systems is typically achieved by spraying the fuel into a polydisperse distribution of droplets evaporating and burning in a turbulent gaseous environment. In view of the nearly unsurmountable difficulties of this two-phase flow, it would be useful to use an experimental arrangement that allow a systematic study of spray evolution and burning in configurations of gradually increasing levels of complexity, starting from laminar sprays to fully turbulent ones. An Electrostatic Spray (ES) of charged droplets lends itself to this type of combustion experiments under well-defined conditions and can be used to synthesize gradually more complex spray environments. In its simplest configuration, a liquid is fed into a small metal tube maintained at several kilovolts relative to a ground electrode few centimeters away. Under the action of the electric field, the liquid meniscus at the outlet of the capillary takes a conical shape, with a thin jet emerging from the cone tip. This jet breaks up farther downstream into a fine spray of charged droplets. Several advantages distinguish the electrospray from alternative atomization techniques: the self-dispersion property of the spray due to coulombic repulsion; the absence of droplet coalescence; the potential control of the trajectories of charged droplets by suitable disposition of electrostatic fields; and the decoupling of atomization, which is strictly electrostatic, from gas flow processes. Furthermore, as recently shown in our laboratory, the electrospray can produce quasi-monodisperse droplets over a very broad size range (1-100 microns). The ultimate objective of this research project is to study the formation and burning of electrosprays of liquid fuels first in laminar regimes and then in turbulent ones. Combustion will eventually be investigated in conditions of three-dimensional droplet-droplet interaction, for which experimental studies have been limited to either qualitative observations in sprays or more quantitative observations on simplified systems consisting of a small number of droplets or droplet arrays. The compactness and potential controllability of this spray generaiton system makes it appealing for studies to be undertaken in the next two years on electrospray combustion in reduced-gravity environments such as those achievable at NASA microgravity test facilities.

Gomez, Alessandro↗

Observations of Fast-Moving Features in the Debris Disk of AU Mic on a Three-Year Timescale: Confirmation and New Discoveries

Context. The nearby and young M star AU Mic is surrounded by a debris disk in which we previously identified a series of large-scale arch-like structures that have never been seen before in any other debris disk and that move outward at high velocities. Aims. We initiated a monitoring program with the following objectives: (1) track the location of the structures and better constrain their projected speeds, (2) search for new features emerging closer in, and ultimately (3) understand the mechanism responsible for the motion and production of the disk features. Methods. AU Mic was observed at 11 different epochs between August 2014 and October 2017 with the IR camera and spectrograph of SPHERE. These high-contrast imaging data were processed with a variety of angular, spectral, and polarimetric differential imaging techniques to reveal the faintest structures in the disk. We measured the projected separations of the features in a systematic way for all epochs. We also applied the very same measurements to older observations from the Hubble Space Telescope (HST) with the visible cameras STIS and ACS. Results. The main outcomes of this work are (1) the recovery of the five southeastern broad arch-like structures we identified in our first study, and confirmation of their fast motion (projected speed in the range 4–12 km/s); (2) the confirmation that the very first structures observed in 2004 with ACS are indeed connected to those observed later with STIS and now SPHERE; (3) the discovery of two new very compact structures at the northwest side of the disk (at 0.40′′ and 0.55′′ in May 2015) that move to the southeast at low speed; and (4) the identification of a new arch-like structure that might be emerging at the southeast side at about 0.4′′ from the star (as of May 2016). Conclusions. Although the exquisite sensitivity of SPHERE allows one to follow the evolution not only of the projected separation, but also of the specific morphology of each individual feature, it remains difficult to distinguish between possible dynamical scenarios that may explain the observations. Understanding the exact origin of these features, the way they are generated, and their evolution over time is certainly a significant challenge in the context of planetary system formation around M stars.

Boccaletti, A.↗

Advection algorithms for quantum neutrino moment transport

Neutrino transport in compact objects is an inherently challenging multidimensional problem. Here, this difficulty is compounded if one includes flavor transformation—an intrinsically quantum phenomenon requiring one to follow the coherence between flavors and thus necessitating the introduction of complex numbers. To reduce the computational burden, simulations of compact objects that include neutrino transport often make use of momentum-angle-integrated moments (the lowest order ones being commonly referred to as the energy density and flux) and these quantities can be generalized to include neutrino flavor, i.e., they become quantum moments. Numerous finite-volume approaches to solving the moment evolution equations for classical neutrino transport have been developed based on solving a Riemann problem at cell interfaces. In this paper we describe our generalization of a Riemann solver for quantum moments, specifically decomposing complex numbers in terms of a (signed) magnitude and phase instead of real and imaginary parts. We then test our new algorithm in numerous cases showing a neutrino fast flavor instability, varying from toy models with analytic solutions to snapshots from neutron star merger simulations. Compared to previous algorithms for neutrino transport with flavor mixing, we find uniformly smaller growth rates of the flavor transformation along with concomitantly larger length-scales, and that the results are a better match with the growth rates seen from multiangle codes.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring the Overabundance of ULXs in Metal- and Dust-Poor Local Lyman Break Analogs

We have studied high-mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z greater than 2) Lyman break galaxies and, within the larger sample of Lyman break analogs (LBAs), they are sufficiently nearby (less than 87 Mpc) to be spatially resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12 + log[O/H] = 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. In this paper, we have performed an in-depth study of the only two LBAs that have spatially resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L(sub X) approximately greater than 10(exp 39) erg s(exp −1); ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on eight detected ULXs. Compared with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al., Haro 11 and VV 114 host approximately equal to 4 times more L(sub X) greater than 10(exp 40) erg s(exp −1) sources than expected given their SFRs. We simulate the effects of source blending from crowded lower-luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF normalizations and bright-end slopes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright-end slope (gamma(sub 2) = 1.90) than the standard XLF (gamma(sub 2) 2.73). If we conservatively assume that the brightest X-ray source from each galaxy is powered by an accreting supermassive black hole rather than an HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF.

galaxies: evolution aEuro" galaxies: individual↗

Cassini atmospheric chemistry mapper. Volume 1. Investigation and technical plan

The Cassini Atmospheric Chemistry Mapper (ACM) enables a broad range of atmospheric science investigations for Saturn and Titan by providing high spectral and spatial resolution mapping and occultation capabilities at 3 and 5 microns. ACM can directly address the major atmospheric science objectives for Saturn and for Titan, as defined by the Announcement of Opportunity, with pivotal diagnostic measurements not accessible to any other proposed Cassini instrument. ACM determines mixing ratios for atmospheric molecules from spectral line profiles for an important and extensive volume of the atmosphere of Saturn (and Jupiter). Spatial and vertical profiles of disequilibrium species abundances define Saturn's deep atmosphere, its chemistry, and its vertical transport phenomena. ACM spectral maps provide a unique means to interpret atmospheric conditions in the deep (approximately 1000 bar) atmosphere of Saturn. Deep chemistry and vertical transport is inferred from the vertical and horizontal distribution of a series of disequilibrium species. Solar occultations provide a method to bridge the altitude range in Saturn's (and Titan's) atmosphere that is not accessible to radio science, thermal infrared, and UV spectroscopy with temperature measurements to plus or minus 2K from the analysis of molecular line ratios and to attain an high sensitivity for low-abundance chemical species in the very large column densities that may be achieved during occultations for Saturn. For Titan, ACM solar occultations yield very well resolved (1/6 scale height) vertical mixing ratios column abundances for atmospheric molecular constituents. Occultations also provide for detecting abundant species very high in the upper atmosphere, while at greater depths, detecting the isotopes of C and O, constraining the production mechanisms, and/or sources for the above species. ACM measures the vertical and horizontal distribution of aerosols via their opacity at 3 microns and, particularly, at 5 microns. ACM recovers spatially-resolved atmospheric temperatures in Titan's troposphere via 3- and 5-microns spectral transitions. Together, the mixing ratio profiles and the aerosol distributions are utilized to investigate the photochemistry of the stratosphere and consequent formation processes for aerosols. Finally, ring opacities, observed during solar occultations and in reflected sunlight, provide a measurement of the particle size and distribution of ring material. ACM will be the first high spectral resolution mapping spectrometer on an outer planet mission for atmospheric studies while retaining a high resolution spatial mapping capability. ACM, thus, opens an entirely new range of orbital scientific studies of the origin, physio-chemical evolution and structure of the Saturn and Titan atmospheres. ACM provides high angular resolution spectral maps, viewing nadir and near-limb thermal radiation and reflected sunlight; sounds planetary limbs, spatially resolving vertical profiles to several atmospheric scale heights; and measures solar occultations, mapping both atmospheres and rings. ACM's high spectral and spatial resolution mapping capability is achieved with a simplified Fourier Transform spectrometer with a no-moving parts, physically compact design. ACM's simplicity guarantees an inherent stability essential for reliable performance throughout the lengthy Cassini Orbiter mission.

Smith, William Hayden↗

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Large Primtive Asteroids: Thermal and Dynamical Context

Primitiveasteroids, most of which are located in the outer belt and Jupiter’s Trojan clouds, provide information related to the origin and evolution of the solar system and the conditions in which the solar nebula was formed. These asteroids are widely thought to be the origin of the least-altered car-bonaceous chondrite that allow us to put crucial con-straints on the current dynamical and thermal theories of the formation and evolution of the early solar sys-tem. The nature of surface composition of large and low-albedo asteroids, like (1) Ceres, (10) Hygiea and (52) Europa, is still under intense debate and different interpretations have beenput forth to explain the ab-sorption features in these objects (e.g.,[1, 2, 3,4,5,6]). Laboratory and spectroscopic experiments on me-teorites that represent all nine carbonaceous chondrite types also found no spectral matches for these large asteroids (e.g.,[7]). Previous studies of asteroid Ceres (the largest primitive asteroid in the solar system) have been conducted to constrain and estimate its surface composition (e.g., [2, 5]). Using linear mixing,[2]found hydroxide brucite, serpentines, and carbonates, to be consistent with Ceres’ ground-based spectra.[5]estimated the surface composition of Ceres and found evidence of widespread NH3-phyllosilicates across its surface using best-fit solutions to Dawn’s NIR spectra. The presence of NH3-phyllosilicates implies that mate-rial from the outer solar system was incorporated into large primitive asteroids, either during their formation at great heliocentric distance or by incorporation of material transported into the Main Belt region. Here we presentnew largeand primitive asteroidsthat share the same spectral similaritieswith the largest asteroid in the solar system, Ceres.We also present the context of these new observations in terms oftheirthermal and dynamicalevolution.Large Primitive Asteroids: Dynamical and Thermal Context:Orbital Distribution of LargePrimitiveAsteroids: Constraining the mineralogy and surface composition of large primitive asteroids will place crucial con-straints on existing dynamical and thermal theories of the formation of the early solar system. We identified several additional asteroids in the Europa-like group in addition to asteroids Europa, Euphrosyne, and Patien-tia, that wereidentified in [8]using NASA IRTFtele-scope (e.g.,Figure 1). These new asteroids along with the already-observed Ceres-and Europa-like group members are localized in the 2.8 < a < 3.4 AU region and characterized by larger sizes, showing an interest-ing orbital distribution (Figure 2). Ceres-likeand Eu-ropa-like groups, which include the largest asteroids in the solar system, show an interesting orbital distribu-tion. These groups are located in the 2.6 < a < 3.6 AU region that contains the snow-line. The snow-line’s location may have been driftedinwards due to the disk’s cooling and evolution[9,10]. Recent dynamical models[11,12]suggested that a substantial fraction of primitive asteroids originated between or beyond the giant planets (a > 5 AU), where water ice would have been stable, and then implanted in the outer Main Belt region because of the giant planets’ growth.Figure 1. Two asteroids, 94 Aurora and 423 Diotima, showing spectra similar to asteroids (1)Ceres and (52)Europa.Thermal modeling and evolution of primitive aster-oids:Primitive water-rich asteroids are thought to be originally composed of mixtures of anhydrous materi-als and water ice that waslater melted by heating sources such as the decay of 26Al, reacting with anhy-drous materials to form H2O/OH-rich minerals. Calcu-lations of the evolution of the temperature and struc-ture of icy planetesimals were performedusing a 1D finite differences thermal evolution model[13, 14] for 26Al-heated planetesimals. In particular, thermally activated compaction due to hot pressing of bodies with an initially unconsolidated porous structure is included. An ice-rich initial composition that leads to a material dominated by phyllosilicates upon aqueous alteration (with 25 vol% H2O and a rock fraction that contains 85 vol% phyllosilicates and 15 vol% olivine upon aque-ous alteration, similar to CIand CM chondrites) was assumed. Atypical initial porosity of 40%[15] is re-duced following the change of the strain rate that is calculated as Voigt approximation from the strain rates of components[16]. Material properties (thermal con-ductivity, density, heat capacity, etc.) correspond to the composition assumed and are adjusted with tempera-ture and porosity. Melting of the water ice as well as water-rock separation are included[14]. Both short-and long-lived radionuclides are consideredas heat sources. Figure 3 shows the maximum temperature calculated as a function of radius and accretion time. A variety of internal structures is obtained, ranging from primordial (no melting of water ice) over partially melted or partially differentiated (melting of water ice, hydration, formation of a rocky core and water ocean below an undifferentiated layer) to completely differ-entiated ones (rocky core, water mantle, Enceladus-like case). The heating and differentiation of planetesi-mals is determined by the availability of 26Al, i.e., by the accretion time t0 relative to the formation of the calcium-aluminum-rich inclusions (CAIs), such that maximum temperatures and structures vary stronglyfor t0< 6 Myr rel. to CAIs. However, for a later accre-tion only the size of the body determines its maximum tempera-ture and structure due to the nearly constant heating by long-lived radionuclides. Average densities of Ceres-and Europa-like group membersimply high-ly porous interiors and, consequently, relatively late accretion at t0> 3 Myr rel. to CAIs with a maximum temperature of < 600 K (Figure3)

D Takir↗