Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Active galaxies”

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 163 records · Page 9

The spectra of IceCube neutrino (SIN) candidate sources – II. Source characterization

ABSTRACT Eight years after the first detection of high-energy astrophysical neutrinos by IceCube, we are still almost clueless as regards to their origin, although the case for blazars being neutrino sources is getting stronger. After the first significant association at the $3\!-\!3.5\, \sigma$ level in time and space with IceCube neutrinos, i.e. the blazar TXS 0506+056 at z = 0.3365, some of us have in fact selected a unique sample of 47 blazars, out of which ∼16 could be associated with individual neutrino track events detected by IceCube. Building upon our recent spectroscopy work on these objects, here we characterize them to determine their real nature and check if they are different from the rest of the blazar population. For the first time we also present a systematic study of the frequency of masquerading BL Lacs, i.e. flat-spectrum radio quasars with their broad lines swamped by non-thermal jet emission, in a γ-ray- and IceCube-selected sample, finding a fraction >24 per cent and possibly as high as 80 per cent. In terms of their broad-band properties, our sources appear to be indistinguishable from the rest of the blazar population. We also discuss two theoretical scenarios for neutrino emission, one in which neutrinos are produced in interactions of protons with jet photons and one in which the target photons are from the broad-line region. Both scenarios can equally account for the neutrino–blazar correlation observed by some of us. Future observations with neutrino telescopes and X-ray satellites will test them out.

79 ASTRONOMY AND ASTROPHYSICS↗

Redshifts of radio sources in the Million Quasars Catalogue from machine learning

ABSTRACT With the aim of using machine learning techniques to obtain photometric redshifts based upon a source’s radio spectrum alone, we have extracted the radio sources from the Million Quasars Catalogue. Of these, 44 119 have a spectroscopic redshift, required for model validation, and for which photometry could be obtained. Using the radio spectral properties as features, we fail to find a model which can reliably predict the redshifts, although there is the suggestion that the models improve with the size of the training sample. Using the near-infrared–optical–ultraviolet bands magnitudes, we obtain reliable predictions based on the 12 503 radio sources which have all of the required photometry. From the 80:20 training–validation split, this gives only 2501 validation sources, although training the sample upon our previous SDSS model gives comparable results for all 12 503 sources. This makes us confident that SkyMapper, which will survey southern sky in the u, v, g, r, i, z bands, can be used to predict the redshifts of radio sources detected with the Square Kilometre Array. By using machine learning to impute the magnitudes missing from much of the sample, we can predict the redshifts for 32 698 sources, an increase from 28 to 74 per cent of the sample, at the cost of increasing the outlier fraction by a factor of 1.4. While the ‘optical’ band data prove successful, at this stage we cannot rule out the possibility of a radio photometric redshift, given sufficient data which may be necessary to overcome the relatively featureless radio spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Possible jet reorientation in low-frequency radio structures of blazars

ABSTRACT Blazars are active galactic nuclei (AGNs) observed at small viewing angles to the observer’s line of sight; hence most of them are compact sources in the radio regime. At low radio frequencies, we may observe additional extended structures, considered to be lobes propagating directly towards Earth. However, there are rare cases of blazars whose morphologies are interpreted in the framework of an episodic jet activity with a reorientation of the jet axis. In this paper, we aim to find other candidates for sources of this kind and present Low Frequency Array (LOFAR) observations of seven blazars with rare morphologies composed of kiloparsec-scale jets and elongated multihundred kiloparsec-scale features. Based on their radio images and physical parameters, we propose that jets could change the direction of propagation out of the sky plane during the lifetime of at least some of these sources and with such we may observe morphologically different AGN types simultaneously, similar to the previously discovered blazar SBS B1646+499. Several plausible origin scenarios are presented, including fast realignment and precession. We also derive the 144 MHz fluxes of the compact and extended structures, the total radio luminosities, the spectral indices for the diffuse haloes, and the core-dominance parameters. We also discuss the possibility that these objects are actually giant blazars with rare Fanaroff–Riley I (FR I) morphology. The previously discovered blazar SBS B1646+499 should be classified as a giant independently of its inclination with regard to the observer.

79 ASTRONOMY AND ASTROPHYSICS↗

Coordinated time variability of multi-phase ultra-fast outflows in J132216.25 + 052446.3

ABSTRACT We present a time variability analysis of broad absorption lines (BAL; spread over the velocity range of 5800–29 000 km s−1) seen in the spectrum of J132216.25 + 052446.3 (zem = 2.04806) at ten different epochs spanning over 19 yr. The strongest absorption component (BAL-A; spread over 5800–9900 km s−1) is made up of several narrow components having velocity separations close to C iv doublet splitting. The C iv, N v, and Si iv absorption from BAL-A show correlated optical depth variability without major changes in the velocity structure. A very broad and shallow absorption (BAL-C; spread over the velocity range 15 000–29 000 km s−1) emerged during our monitoring period coinciding with a dimming episode of J1322 + 0524. All the identified absorption lines show correlated variability with the equivalent widths increasing with decreasing flux. This together with the C iv emission line variability is consistent with ionization being the main driver of the correlated variability. The observed UV-continuum variations are weaker than what is required by the photoionization models. This together with a scatter in the C iv equivalent width at a given continuum flux can be understood if variations of the C iv ionizing photons are much larger than that of the UV continuum, the variations in the ionizing photon and UV fluxes are not correlated and/or the covering factor of the flow varies continuously. We suggest BAL-A is produced by a stable clumpy outflow located beyond the broad emission line region and BAL-C is a newly formed wind component located near the accretion disc and both respond to changes in the ionizing continuum.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiple giant eruptions and X-ray emission in the recoiling AGN/LBV candidate SDSS1133

We present a comprehensive analysis of 20 yr worth of multicolour photometric light curves, multiepoch optical spectra, and X-ray data of an off-nuclear variable object SDSS1133 in Mrk 177 at z = 0.0079. The UV-optical light curves reveal that SDSS1133 experienced four outbursts in 2001, 2014, 2019, and 2021. The persistent UV-optical luminosity in the non-outbursting state is ~1041 erg s -1 with small-scale flux variations, and peak luminosities during the outbursts reach ~10 42 erg s -1 . The optical spectra exhibit enduring broad hydrogen Balmer P-Cygni profiles with the absorption minimum at ~-2000 km s -1 , indicating the presence of fast-moving ejecta. Chandra detected weak X-ray emission at a 0.3-10-keV luminosity of L X = 4 × 10 38 erg s -1 after the 2019 outburst. These lines of evidence suggests that SDSS1133 is an extreme luminous blue variable (LBV) star experiencing multiple giant eruptions with interactions of the ejected shell with different shells and/or circumstellar medium (CSM), and disfavours the recoiling active galactic nuclei scenario suggested in the literature. We suggest that pulsational pair-instability may provide a viable explanation for the multiple energetic eruptions in SDSS1133. If the current activity of SDSS1133 is a precursor of a supernova explosion, we may be able to observe a few additional giant eruptions and then the terminal supernova explosion or collapse to a massive black hole in future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

The relation between optical and γ-ray emission in BL Lac sources

ABSTRACT The relativistic jets produced by some Active Galactic Nuclei (AGNs) are among the most efficient persistent sources of non-thermal radiation and represent an ideal laboratory for studying high-energy interactions. In particular, when the relativistic jet propagates along the observer’s line of sight, the beaming effect produces dominant signatures in the observed spectral energy distribution (SED), from the radio domain up to the highest energies, with the further possibility of resulting in radiation-particle multimessenger associations. In this work, we investigate the relationships between the emission of γ rays and the optical spectra of a sample of AGN, selected from BL Lac sources detected by the Fermi Large Area Telescope (Fermi-LAT). We find that there is a close relationship between the optical and γ-ray spectral indices. Despite all the limitations due to the non-simultaneity of the data, this observation strongly supports a substantial role of Synchrotron-Self Compton (SSC) radiation in a single zone leptonic scenario for most sources. This result simplifies the application of theoretical models to explore the physical parameters of the jets in this type of sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficient analysis routines for single and double peaked Type 2 AGN spectra

Driven by the imminent need to rapidly process and classify millions of AGN spectra drawn from next generation astronomical facilities, we present a spectral fitting routine for Type 2 AGN spectra optimized for high volume processing, using the quasar spectral fitting library (qsfit). We analyse an optically selected sample of 813 luminous Type 2 AGN spectra at z < 0.83 from the Sloan Digital Sky Survey (SDSS) to qualify its performance. We report a median narrow line H α/H β Balmer decrement of 4.5 ± 0.8, alluding to the presence of dust in the narrow line region (NLR). We publish a specialized QSFIT fitting routine for high signal-to-noise ratio spectra and general fitting routine for double peaked Type 2 AGN spectra applied on a subsample of 45 spectra from our parent sample. We report a median red and blue peak velocity separation of 390 ± 60 kms −1 . No trend is found for red or blue peaks to exhibit systematically different luminosity or ionization properties. Emission line diagnostics show that the double peaks in all sources are illuminated by an AGN-powered ionizing continuum. Finally, we examine the morphology of host galaxies of our double peaked sample. We find double peaked Type 2 AGN reside in merging systems at a comparable frequency to single peaked AGN. This suggests that the double peaked AGN phenomenon is likely to have a bi-conical outflow origin in the majority of cases. We publicly release the code used for spectral analysis and produced catalogues used in this work.

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray emission of the radio-loud quasar SDSS J121426.52+140258.9: independent variations between optical/UV and X-ray emission

To understand the X-ray emission of active galactic nuclei (AGNs), we explored the optical-to-X-ray variation correlation of a radio-loud quasar (RLQ) SDSS J121426.52+140258.9 (hereafter J1214+1402) with multi-epoch observations by the Swift and XMM–Newton telescopes. With the historical multiband data, we found that the infrared-to-X-ray flux of RLQ J1214+1402 should not be dominated by the beamed-jet emission. The Swift optical/UV and X-ray light curves showed that J1214+1402 has two optical states, with low flux before 2014 April 8 and high flux after 2014 June 11, but has no significant X-ray variations during the time range between 2007 March 9 and 2014 August 4. This result was supported by the XMM–Newton observations in the time overlapping with Swift. Interestingly, the early XMM–Newton data prior to the Swift time present two unusual emission epochs when J1214+1402 has relatively low optical fluxes but has the brightest X-ray fluxes. The overall independence of optical-to-X-ray variation seems hard to describe by the disc–corona model. With the X-ray spectral fitting, we find that the soft X-ray excess in J1214+1402 appears only during the high optical state when the X-ray emission is at a low state. The soft X-ray excess in J1214+1402 is difficult to explain by an ionized accretion disc; instead, it may be related to the warm corona.

79 ASTRONOMY AND ASTROPHYSICS↗

The H I content of red geyser galaxies

Red geysers are a specific type of quiescent galaxy, denoted by twin jets emerging from their galactic centres. These bisymmetric jets possibly inject energy and heat into the surrounding material, effectively suppressing star formation by stabilizing cool gas. In order to confirm the presence and evolutionary consequences of these jets, this paper discusses the scaling, stacking, and conversion of 21-cm H I flux data sourced from the H I -MaNGA survey into H I gas-to-stellar mass (G/S) spectra. Our samples were dominated by non-detections, or galaxies with weak H I signals, and consequently by H I upper limits. The stacking technique discussed successfully resolved emission features in both the red geyser G/S spectrum and the control sample G/S spectrum. From these stacked spectra, we find that on average, red geyser galaxies have G/S of 0.086 ± 0.011 (random) + 0.029 (systematic), while non-red geyser galaxies of similar stellar mass have a G/S ratio of 0.039 ± 0.018 (random) + 0.013 (systematic). Therefore, we find no statistically significant evidence that the H I content of red geysers is different from the general quiescent population.

79 ASTRONOMY AND ASTROPHYSICS↗

XMM–Newton study of the Sculptor dwarf spheroidal galaxy

ABSTRACT We analysed observations with XMM–Newton in the field of the Sculptor dwarf spheroidal galaxy (dSph). The aim of the study was the classification of X-ray binaries and accreting white dwarfs belonging to the Sculptor dSph. Using different methods of X-ray timing and spectral analyses, together with an extensive multiwavelength study of the optical and infrared counterparts of the X-ray sources, we classified the sources detected with XMM–Newton in the field of the Sculptor dSph. The long-term variability of the sources has been studied over two XMM–Newton observations. None of the members of the Sculptor dSph show significant long-term variability over these two observations. We also searched for periodicity and pulsation using the Lomb–Scargle and Rayleigh Z$^{2}_{n}$ techniques. No signals of pulsation or periodicity have been found for the X-ray sources. The results show the presence of a noticeable number of background X-ray sources in the field of this galaxy. We classified 43 sources as active galactic nuclei, galaxies, and galaxy candidates. Three Galactic foreground stars have been identified in the field of the Sculptor dSph, and one of them is an M-dwarf candidate. Moreover, we classified four symbiotic-star candidates and three quiescent low-mass X-ray binary candidates in the Sculptor dSph. The luminosity of these X-ray sources is ∼1033−35 erg s−1.

79 ASTRONOMY AND ASTROPHYSICS↗

Production efficiencies of sound waves in the intracluster medium driven by AGN jets

Feedback from active galactic nuclei (AGNs) is believed to be the most promising solution to the cooling flow problem in cool-core clusters. Dissipation of sound waves is considered as one of the possible heating mechanisms; however, its relative contribution to heating remains unclear. To estimate the energy budget for heating, we perform 3D hydrodynamic simulations of AGN jet injections in a Perseus-like cluster and quantify the amount of energy stored in the forms of weak shocks and waves. For this work, e find that, for a single jet injection with typical parameters in cool-core clusters, ${\sim}9{{\ \rm per\ cent}}$ of the total jet energy is stored in compressional waves (including both shocks and waves). However, due to the destructive effects among randomly phased waves as well as the dissipation of shock energies, in our simulations including self-regulated AGN feedback, no more than $3{{\ \rm per\ cent}}$ of the total injected energy goes into compressional waves. We further separate the energy contribution from shocks and waves and find that, for a single outburst, the shocks can only contribute to ${\sim}20{-}30{{\ \rm per\ cent}}$ of the total compressional energy in the inner radii and quickly dissipate away. In the self-regulated case where shocks are repeatedly generated, shocks completely dominate over sound waves in the inner region and can still provide ${\sim}40{-}50{{\ \rm per\ cent}}$ of the total compressional energy at outer radii. Our results suggest that the production of sound waves is not as efficient as what was previously found, and thus sound wave dissipation may be a subdominant source of heating in cool-core clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

A practicable estimation of opening angle of dust torus in Type-1.9 AGN with double-peaked broad Hα

In this manuscript, an independent method is proposed to estimate opening angle of dust torus in AGN, through unique properties of Type-1.9 AGN with double-peaked broad Hα (Type-1.9 DPAGN) coming from central accretion disc. Type-1.9 AGN without broad Hβ can be expected by the commonly accepted unified model of AGN, considering central BLRs seriously obscured by dust torus with its upper boundary in the line of sight. For the unique Type-1.9 DPAGN, accretion disc originations of double-peaked broad Hα can be applied to determine the inclination angle of the central accretion disc, which is well accepted as substitute of the half opening angle of the central dust torus. Then, among low redshift Type-1.9 DPAGN in SDSS, SDSS J1607+3319 at redshift 0.063 is collected, and the half opening angle of the central dust torus is determined to be around 46 ± 4°, after considering the disfavoured BBH system to explain the double-peaked broad Hα through long-term none variabilities and disfavoured local physical conditions to explain the disappearance of broad Hβ through virial BH mass properties. The results indicate that more detailed studying on dust torus of AGN can be appropriately done through Type-1.9 DPAGN in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗

Time variability of ultra fast BAL outflows using SALT: C IV equivalent width analysis

We study the time variability (over ≤7.3 yr) of ultra fast outflows (UFOs) detected in a sample of 64 C IV broad absorption line (BAL) quasars (with 80 distinct BAL components) monitored using the Southern African Large Telescope. By comparing the properties of the quasar in our sample with those of a control sample of non-BAL quasars, we show that the distributions of black hole mass are different and the bolometric luminosities and optical photometric variations of UFO BAL quasars are slightly smaller compared to that of non-BAL quasars. The detection fraction of C IV equivalent width (W) variability (∼95 per cent), the fractional variability amplitude $\left(\frac{\Delta W}{W}\right)$ and the fraction of ‘highly variable’ BAL (i.e. $\big|$$\frac{\Delta W}{W}$$\big|$ > 0.67) components (∼33 per cent) are higher in our sample compared to the general BAL population. The scatter in $\frac{\Delta W}{W}$ and the fraction of ‘highly variable’ BALs increase with the time-scale probed. The $\frac{\Delta W}{W}$ distribution is asymmetric at large time scales. We attribute this to the BAL strengthening time-scales being shorter than the weakening time-scales. The BAL variability amplitude correlates strongly with the BAL properties compared to the quasar properties. BALs with low W, high-velocity, shallow profiles, and low-velocity width tend to show more variability. When multiple BAL components are present, a correlated variability is seen between low- and high-velocity components with the latter showing a larger amplitude variations. We find an anticorrelation between the fractional variations in the continuum flux and W. While this suggests photoionization induced variability, the scatter in continuum flux is much smaller than that of W.

79 ASTRONOMY AND ASTROPHYSICS↗

OzDES Reverberation Mapping Program: Hβ lags from the 6-yr survey

ABSTRACT Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This R–L relation is used to estimate single-epoch virial black hole masses, and has been proposed to use to standardize AGN to determine cosmological distances. We present reverberation measurements made with Hβ from the 6-yr Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at 0.12 < z < 0.71, probing higher redshifts than the bulk of Hβ measurements made to date. Our fit to the R–L relation has a slope of α = 0.41 ± 0.03 and an intrinsic scatter of σ = 0.23 ± 0.02 dex. The results from our multi-object spectroscopic survey are consistent with previous measurements made by dedicated source-by-source campaigns, and with the observed dependence on accretion rate. Future surveys, including LSST, TiDES, and SDSS-V, which will be revisiting some of our observed fields, will be able to build on the results of our first-generation multi-object reverberation mapping survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing quasar viewing angle with the variability structure function

Given the anisotropic emission from quasar accretion discs, their viewing angle affects estimates of the quasar luminosity, black hole mass, and Eddington ratio. Discs appear overluminous when viewed pole-on and underluminous when viewed at high inclination. In radio-quiet quasars, the viewing angle is usually unknown, although spectroscopic indicators have been proposed. Here, we use a recently discovered universality in the variability structure function (SF) of quasar light curves (LCs), where all quasars show the same SF when clocks run in units of orbital time-scale. As an offset from the mean relation can be caused by incorrect orbital time-scales and thus incorrect luminosities, we correlate these offsets with suggested inclination indicators. We derive SFs from NASA/Asteroid Terrestrial-impact Last Alert System (ATLAS) LCs spanning ~6 yr of observation, using a sample of 183 luminous quasars with measured H β lines as well as 753 quasars with C $\tiny{IV}$ and Mg $\tiny{II}$ lines. Starting from the proposed orientation indicators, we expect quasars with narrower H β lines and with more blueshifted C $\tiny{IV}$ lines to be viewed more pole-on and thus appear overluminous. In contrast, our SF analysis finds that presumed pole-on discs appear underluminous, consistently for both line indicators. We discuss possible explanations for the behaviour of quasars with highly blueshifted C $\tiny{IV}$ lines irrespective of inclination angle, including dusty outflows that might render the accretion disc underluminous and flatter disc temperature profiles with longer orbital time-scales than in thin-disc models but reach no satisfying conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolution and feedback of AGN jets of different cosmic ray composition

Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic ray proton (CRp), cosmic ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics, and heating, and model their expected cavity-power versus radio-luminosity relation (P cav –L R ). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within ∼30 Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the P cav –L R plane, but the values are broadly consistent with observed ranges for FR-I sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗

Follow-up observations for IceCube-170922A: Detection of rapid near-infrared variability and intensive monitoring of TXS 0506+056

Abstract We present our follow-up observations to search for an electromagnetic counterpart of the IceCube high-energy neutrino IceCube-170922A. Monitoring observations of a likely counterpart, TXS 0506+056, are also described. First, we quickly took optical and near-infrared images of seven flat-spectrum radio sources within the IceCube error region right after the neutrino detection and found a rapid flux decline of TXS 0506+056 in Kanata/HONIR J-band data. Motivated by this discovery, intensive follow-up observations of TXS 0506+056 were continuously performed, including our monitoring imaging observations, spectroscopic observations, and polarimetric observations in optical and near-infrared wavelengths. TXS 0506+056 showed a large-amplitude (∼1.0 mag) variability in a time scale of several days or longer, although no significant variability was detected in a time scale of a day or shorter. TXS 0506+056 also showed a bluer-when-brighter trend in optical and near-infrared wavelengths. Structure functions of the variabilities were examined and indicate that TXS 0506+056 is not a special blazar in terms of optical variability. Polarization measurement results of TXS 0506+056 are also discussed.

(galaxies:) BL Lacertae objects: general↗

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗