Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ACCRETION DISKS”

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 577 records · Page 32

iPTF Discovery of the Rapid "Turn-On" of a Luminous Quasar

We present a radio-quiet quasar at z = 0.237 discovered ''turning on'' by the intermediate Palomar Transient Factory (iPTF). The transient, iPTF 16bco, was detected by iPTF in the nucleus of a galaxy with an archival Sloan Digital Sky Survey spectrum with weak narrow-line emission characteristic of a low-ionization nuclear emission-line region (LINER). Our follow-up spectra show the dramatic appearance of broad Balmer lines and a power-law continuum characteristic of a luminous (L(sub bol) approx. equal 10(exp 45) erg/s) type 1 quasar 12 yr. later. Our photometric monitoring with PTF from 2009-2012 and serendipitous X-ray observations from the XMM-Newton Slew Survey in 2011 and 2015 constrain the change of state to have occurred less than 500 days before the iPTF detection. An enhanced broad H(alpha)/[O III] (lambda)5007 line ratio in the type 1 state relative to other changing-look quasars also is suggestive of the most rapid change of state yet observed in a quasar. We argue that the greater than 10 increase in Eddington ratio inferred from the brightening in UV and X-ray continuum flux is more likely due to an intrinsic change in the accretion rate of a preexisting accretion disk than an external mechanism such as variable obscuration, microlensing, or the tidal disruption of a star. However, further monitoring will be helpful in better constraining the mechanism driving this change of state. The rapid ''turn-on'' of the quasar is much shorter than the viscous infall timescale of an accretion disk and requires a disk instability that can develop around a approx. equal 10(exp 8)solar mass black hole on timescales less than 1 yr.

Gezari, S.↗

The standard model and some new directions

A 'standard' model of Active Galactic Nuclei (AGN), based upon a massive black hole surrounded by a thin accretion disk, is defined. It is argued that, although there is good evidence for the presence of black holes and orbiting gas, most of the details of this model are either inadequate or controversial. Magnetic field may be responsible for the confinement of continuum and line-emitting gas, for the dynamical evolution of accretion disks and for the formation of jets. It is further argued that gaseous fuel is supplied in molecular form and that this is responsible for thermal re-radiation, equatorial obscuration and, perhaps, the broad line gas clouds. Stars may also supply gas close to the black hole, especially in low power AGN and they may be observable in discrete orbits as probes of the gravitational field. Recent observations suggest that magnetic field, stars, dusty molecular gas and orientation effects must be essential components of a complete description of AGN. The discovery of quasars with redshifts approaching 5 is an important clue to the mechanism of galaxy formation.

Blandford, R. D.↗

Automatic Search of Cataclysmic Variables Based on LightGBM in LAMOST-DR7

The search for special and rare celestial objects has always played an important role in astronomy. Cataclysmic Variables (CVs) are special and rare binary systems with accretion disks. Most CVs are in the quiescent period, and their spectra have the emission lines of Balmer series, HeI, and HeII. A few CVs in the outburst period have the absorption lines of Balmer series. Owing to the scarcity of numbers, expanding the spectral data of CVs is of positive significance for studying the formation of accretion disks and the evolution of binary star system models. At present, the research for astronomical spectra has entered the era of Big Data. The Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST) has produced more than tens of millions of spectral data. the latest released LAMOST-DR7 includes 10.6 million low-resolution spectral data in 4926 sky regions, providing ideal data support for searching CV candidates. To process and analyze the massive amounts of spectral data, this study employed the Light Gradient Boosting Machine (LightGBM) algorithm, which is based on the ensemble tree model to automatically conduct the search in LAMOST-DR7. Finally, 225 CV candidates were found and four new CV candidates were verified by SIMBAD and published catalogs. This study also built the Gradient Boosting Decision Tree (GBDT), Adaptive Boosting (AdaBoost), and eXtreme Gradient Boosting (XGBoost) models and used Accuracy, Precision, Recall, the F1-score, and the ROC curve to compare the four models comprehensively. Experimental results showed that LightGBM is more efficient. The search for CVs based on LightGBM not only enriches the existing CV spectral library, but also provides a reference for the data mining of other rare celestial objects in massive spectral data.

79 ASTRONOMY AND ASTROPHYSICS↗

Temporal behavior of Hercules X-1 - The long-term variability of its pulse profile and the 35 day X-ray intensity modulation

Data collected by the low-energy detectors of the Hard X-Ray and Low Energy Gamma-Ray Experiment on HEAO 1 were analyzed to study the variability in the temporal behavior of Hercules X-1. Data obtained during thirteen pointed observations distributed in 35-day phase provided details on the time and the energy dependence of the pulse shape and the spectra. It was found that a structural change at the inner edge of the accretion disk, but not at the large radii, was likely to cause the 35-day as well as the observed short-term rapid variations of the pulse profile, which had a considerably more pronounced time dependence than energy dependence. It is suggested that the variable mass-channeling scheme and an intervening effect at the inner accretion disk are directly responsible for the 35-day X-ray modulation.

Soong, Y.↗

Turbulization of shear flows in astrophysics

The problem of instability in shear flows is treated in a number of works which present a global analysis of irrotational (potential) modes in a 2-D compressible shear flow, implying that, in the presence of reflecting boundaries, those modes can increase and that the characteristic time of the increase is several orders as high as the dynamic time in the medium. Due to the slowness of the increase of irrotational modes and to the problematic character of the theory's applicability for thin accretion disks, an attempt was made to develop an alternative theory of turbulization shear flows. It is quite contrary to the above theory. Instead of the global analysis, a local one is made far from the flow boundaries. Instead of irrotational perturbations, vortical ones are examined. Lastly, the medium is considered to be incompressible. The flow in the accretion disk can be modeled by a plane Couette flow, without regard for boundary conditions: it can be modeled by a plane shear flow in infinite incompressible fluid. According to this scenario, there is a critical perturbation level in the free shear flow. If the level of initial perturbations exceeds the critical one, the flow is turbulized.

Chagelishvili, G. D.↗

Multi-Epoch XMM Observations of NGC4258

The goal of this project is to better understand whether the well known, warped accretion disk around the supermassive black hole in NGC4258 is responsible for the observed X-ray absorbing column of gas. The active galactic nuclei (AGN) paradigm suggests that doughnut-like tori obscure black holes, not warped disks. Using radio interferometry, we have mapped the disk and its orientation. Through time monitoring radio observations we are following the rotation and variability in the disk structure. The disk crosses the line of sight to the black hole, if it is clumpy, then we will observe variations in the X-ray spectrum of NGC4258 on month to year time scales. We have reduced 14 epochs of very long base interferometry (VLBI) observations that we use to characterize the structure of the maser-laden accretion disk and its evolution. Five epochs of data remain to be reduced in the next year. We have received distributions of data for two of our four intended epochs of X-ray Multimirror Mission (XMM) observation. Because the first two epochs are relatively close in time, and because community understanding of how to calibrate XMM data has until recently been relatively immature, we are waiting until we receive the third epoch of data before we reduce the several datasets. Delivery of the third epoch should take place late this year. The last epoch of XMM observation will take place in mid-2002. In the next year, we will begin to combine constraints on disk clumpiness with estimates of X-ray variability observed by XMM.

Greenhill, Lincoln J.↗

Accretion Flows in Magnetic White Dwarf Systems

We received Type A and B funding under the NASA Astrophysics Data Program for the analysis and interpretation of hard x-ray data obtained by the Rossi X-ray Timing Explorer and other NASA sponsored missions for Intermediate Polars (IPS) and Polars. For some targets, optical data was available. We reduced and analyzed the X-ray spectra and the X-ray and optical (obtained at the Cerro Tololo Inter-American Observatory) timing data using detailed shock models (which we constructed) to place constraints on the properties of the accreting white dwarfs, the high energy emission mechanisms of white dwarfs, and the large-scale accretion flows of Polars and IPS. IPS and Polars are white dwarf mass-transfer binaries, members of the larger class of cata,clysmic variables. They differ from the bulk of the cataclysmic variables in that they contain strongly magnetic white dwarfs; the white dwarfs in Polars have B, = 7 to 230 MG and those in IPS have B, less than 10 MG. The IPS and Polars are both examples of funneled accretion flows in strong magnetic field systems. The IPS are similar to x-ray pulsars in that accretion disks form in the systems which are disrupted by the strong stellar magnetic fields of the white dwarfs near the stellar surface from where the plasma is funneled to the surface of the white dwarf. The localized hot spots formed at the footpoints of the funnels coupled with the rotation of the white dwarf leads to coherent pulsed x-ray emission. The Polars offer an example of a different accretion topology; the magnetic field of the white dwarf controls the accretion flow from near the inner Lagrangian point of the system directly to the stellar surface. Accretion disks do not form. The strong magnetic coupling generally leads to synchronous orbital/rotational motion in the Polars. The physical system in this sense resembles the Io/Jupiter system. In both IPS and Polars, pulsed emission from the infrared to x-rays is produced as the funneled flows merge onto the white dwarfs through the formation of strong radiating shock waves. A comparative study of the IPS and Polars can elucidate the primary effects of the magnetic fields on the dynamics and thermodynamics in accreting white dwarf systems.

Imamura, James N.↗

Poynting Robertson Battery and the Chiral Magnetic Fields of AGN Jets

We propose that the magnetic fields in the accretion disks of active galactic nuclei (AGNs) are generated by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with the AGN photons (the Poynting Robertson battery). This process provides a unique relation between the polarity of the poloidal B field to the angular velocity Omega of the accretion disk (B is parallel to Omega), a relation absent in the more popular dynamo B-field generation. This then leads to a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry comes about by chance being approx.0.06 %. This lends support to the hypothesis that the universe is seeded by B fields that are generated in AGNs via this mechanism and subsequently injected into intergalactic space by the jet outflows.

Kazanas, Demosthenes↗

Observations of the bright novalike variable IX Velorum with the Hopkins Ultraviolet Telescope

The Hopkins Ultraviolet Telescope, an experiment flown on the Space Shuttle as part of the Astro-1 mission, was used to obtain a spectrum of the novalike variable IX Vel (= CPD -48 deg 1577) in the wavelength range 830-1860 A. The observation revealed a rich absorption-line and continuum spectrum that peaks near 1050 A at a flux of 1.6 x 10(exp -11) ergs/sq cm/s/A. In the sub-Lyman-alpha region, some of the more prominent absorption lines are S VI lambda lambda-933, 945, C III lambda-977, Lyman-beta, O VI lambda lambda-1032, 1038, P V lambda lambda-1118, 1128, and C III lambda-1176. No emission was detected below the Lyman limit. The overall continuum shape of IX Vel in the FUV can be approximated using models of an optically thick accretion disk in which the integrated spectrum has been constructed by summing model stellar atmospheres or proper disk model spectra. However, if the distance to IX Vel is approximately 95 pc, standard disk models without reddening cannot simultaneously reproduce the color and flux in the UV. While interstellar reddening can reconcile this difference, the amount of reddening appears inconsistent with the absence of a 2200 A bump in the spectrum and the very low H I column density measured along the line of sight. Improved fits to the data can be obtained by modifying the accretion disk stucture within three white dwarf radii. None of the models reproduces the profiles of the Li- and Na-like ions, which are observed as strong but relatively narrow absorption lines, and which are almost surely due to a wind above the disk.

Long, Knox S.↗

Simultaneous NuSTAR/Chandra Observations of the Bursting Pulsar GRO J1744-28 During its Third Reactivation

We report on a 10 ks simultaneous Chandra/High Energy Transmission Grating (HETG)-Nuclear Spectroscopic Telescope Array (NuSTAR) observation of the Bursting Pulsar, GRO J1744-28, during its third detected outburst since discovery and after nearly 18 yr of quiescence. The source is detected up to 60 keV with an Eddington persistent flux level. Seven bursts, followed by dips, are seen with Chandra, three of which are also detected with NuSTAR. Timing analysis reveals a slight increase in the persistent emission pulsed fraction with energy (from 10% to 15%) up to 10 keV, above which it remains constant. The 0.5-70 keV spectra of the persistent and dip emission are the same within errors and well described by a blackbody (BB), a power-law (PL) with an exponential rolloff, a 10 keV feature, and a 6.7 keV emission feature, all modified by neutral absorption. Assuming that the BB emission originates in an accretion disk, we estimate its inner (magnetospheric) radius to be about 4 x 10(exp 7) cm, which translates to a surface dipole field B approximately 9 x 10(exp 10) G. The Chandra/HETG spectrum resolves the 6.7 keV feature into (quasi-)neutral and highly ionized Fe XXV and Fe XXVI emission lines. XSTAR modeling shows these lines to also emanate from a truncated accretion disk. The burst spectra, with a peak flux more than an order of magnitude higher than Eddington, are well fit with a PL with an exponential rolloff and a 10 keV feature, with similar fit values compared to the persistent and dip spectra. The burst spectra lack a thermal component and any Fe features. Anisotropic (beamed) burst emission would explain both the lack of the BB and any Fe components.

pulsars: general aEuro" stars: individual (GRO J17↗

Spectral energy distributions of the brightest Palomar-Green quasars at intermediate redshifts

We have combined low-dispersion International Ultraviolet Explorer (IUE) spectra with the optical/near-IR spectrophotometry of Neugebauer et al. (1987) in order to study the spectral energy distributions of seven of the brightest Palomar-Green (PG) quasars at intermediate redshifts (Z(sub em) greater than or equal to 0.9 and less than or equal to 1.5). Some of these PG quasars are barely detectable in long IUE exposures, so we have used the Gaussian Extraction (GEX) technique to maximize the signal-to-noise of the IUE data, and we have co-added all spectra available from the IUE archive for each QSO unless the ultraviolet spectra varied significantly from one exposure to the next. We have corrected the spectral energy distributions for Milky Way reddening using the observed neutral hydrogen column densities on each sight line and the gas-to-dust relation recently derived by Diplas & Savage. Six of the seven quasars are detected down to lambda much less than 700 A in the rest frame, and consequently continuum reddening due to dust in the immediate vicinity of the quasar can have a dramatic effect on the spectral energy distributions. In order to explore the possible importance of intrinsic continuum reddening, we have assembled a heuristic extinction curve which extends to lambda much less than 912 A. Using this heuristic extinction curve, we derive reasonable upper limits on the intrinsic E(B-V) for each quasar. We briefly discuss some of the implications of the derived intrinsic continuum reddening limits. We use geometrically thin accretion disk models to derive the black hole masses and accretion rates implied by the spectral energy distributions. Even if we neglect intrinsic reddening, we find that a large fraction of the quasars require super-Eddington accretion rates (which is not consistent with the thin disk assumption). Comparison of the data in this paper to a large body of data from the literature on the accretion disk M(sub BH) - M dot grid calculated by Wandel & Petrosian reveals that our quasars are among the brightest in the sky at 1450 A, and ostensibly suggests that the fraction of quasars which require super-Eddington accretion rates is much smaller than the fraction that we derive from our data alone. However, intrinsic continuum reddening has been ignored in this comparison, and a small amount of intrinsic reddening will push more of the quasars into the super-Eddington regime. We also plot the recent reverberation monitoring results on NGC 5548 and NGC 3783 on the Wandel & Petrosian grid, and we find that these Seyfert galaxies appear to vary along lines of constant M(sub BH). Continuum flux from two of the quasars in our main sample, PG 1338+416 and PG 1630+377, is detected at lambda(sub rest) less than 584 A. These quasars can in principle be used for the He I Gunn-Peterson test, but the S/N of IUE spectra of individual objects is usually too low to place interesting limits on the Gunn-Peterson optical depth. In order to improve the S/N, we have formed a composite spectrum from the spectra of five quasars detected with IUE at lambda(sub rest) less than 584 A, and we have used this composite spectrum to place a tighter limit on tau(sub GP, He I). We briefly discuss intermediate-redshift Lyman limit systems (Z(sub LL) greater than or equal to 0.5 and less than or equal to 1.6) detected in the IUE spectra of five quasars, including lower limits on N(H I) in each Lyman limit system.

Tripp, Todd M.↗

The Invariant Twist of Magnetic Fields in the Relativistic Jets of Active Galactic Nuclei

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in Active Galactic Nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1 %. This lends support to the hypothesis that the Universe is seeded by B fields that are generated in AGN via this mechanism

Contopoulos, Ioannis↗

Iron K lines from low-mass X-ray binaries

Models are presented for the 6-7 keV iron line emission from low-mass X-ray binaries. A simplified model for an accretion disk corona is used to examine the dependence of the observable line properties, line width and mean energy, on the radial distance of the emission region from the X-ray source, and on the fraction of the X-rays from the source which reach the disk surface. The effects of blending of multiple line components and of Comptonization of the line profile are included in numerical calculations of the emitted profile shape. The results of these calculations, when compared with the line properties observed from several low-mass X-ray binaries, suggest that the broadening is dominated either by rotation or by Compton scattering through a greater optical depth than is expected from an accretion disk corona.

Kallman, T.↗

A reprocessing model for the ultraviolet and optical light from 4U 1820-30

We show that the recently discovered optical and ultraviolet light from the X-ray burst source 4U 1820-30 in the globular cluster NGC 6624 is due to reprocessing of the X-rays in the outer regions of an optically thick, geometrically thin accretion disk. We suggest that observation of orbital modulation of the reprocessed light, due to the variable contribution made as the heated face of the companion turns toward and away from the observer, would provide constraints on the inclination of the binary orbit, and we suggest that detection in the reprocessed flux of the 'red noise' already observed in the X-rays would provide useful constraints on the geometry and physics of the accretion disk.

Arons, Jonathan↗

Quasar Astrophysics with the Space Interferometry Mission

Optical astrometry of quasars and active galaxies can provide key information on the spatial distribution and variability of emission in compact nuclei. The Space Interferometry Mission (SIM PlanetQuest) will have the sensitivity to measure a significant number of quasar positions at the microarcsecond level. SIM will be very sensitive to astrometric shifts for objects as faint as V = 19. A variety of AGN phenomena are expected to be visible to SIM on these scales, including time and spectral dependence in position offsets between accretion disk and jet emission. These represent unique data on the spatial distribution and time dependence of quasar emission. It will also probe the use of quasar nuclei as fundamental astrometric references. Comparisons between the time-dependent optical photocenter position and VLBI radio images will provide further insight into the jet emission mechanism. Observations will be tailored to each specific target and science question. SIM will be able to distinguish spatially between jet and accretion disk emission; and it can observe the cores of galaxies potentially harboring binary supermassive black holes resulting from mergers.

instrumentation↗

Physical Structure of Four Symbiotic Binaries

Disk accretion powers many astronomical objects, including pre-main sequence stars, interacting binary systems, and active galactic nuclei. Unfortunately, models developed to explain the behavior of disks and their surroundings - boundary layers, jets, and winds - lack much predictive power, because the physical mechanism driving disk evolution - the viscosity - is not understood. Observations of many types of accreting systems are needed to constrain the basic physics of disks and provide input for improved models. Symbiotic stars are an attractive laboratory for studying physical phenomena associated with disk accretion. These long period binaries (P(sub orb) approx. 2-3 yr) contain an evolved red giant star, a hot companion, and an ionized nebula. The secondary star usually is a white dwarf accreting material from the wind of its red giant companion. A good example of this type of symbiotic is BF Cygni: our analysis shows that disk accretion powers the nuclear burning shell of the hot white dwarf and also manages to eject material perpendicular to the orbital plane (Mikolajewska, Kenyon, and Mikolajewski 1989). The hot components in other symbiotic binaries appear powered by tidal overflow from a very evolved red giant companion. We recently completed a study of CI Cygni and demonstrated that the accreting secondary is a solar-type main sequence star, rather than a white dwarf (Kenyon et aL 1991). This project continued our study of symbiotic binary systems. Our general plan was to combine archival ultraviolet and optical spectrophotometry with high quality optical radial velocity observations to determine the variation of line and continuum sources as functions of orbital phase. We were very successful in generating orbital solutions and phasing UV+optical spectra for five systems: AG Dra, V443 Her, RW Hya, AG Peg, and AX Per. Summaries of our main results for these systems appear below. A second goal of our project was to consider general models for the outbursts of symbiotic stars, with an emphasis on understanding the differences between disk-driven and nuclear-powered eruptions.

Kenyon, Scott J.↗

The average X-ray/gamma-ray spectra of Seyfert galaxies from Ginga and OSSE and the origin of the cosmic X-ray background

We have obtained the first average 2-500 keV spectra of Seyfert galaxies, using the data from Ginga and Compton Gamma-Ray Observatory's (CGRO) Oriented Scintillation Spectrometer Experiment (OSSE). Our sample contains three classes of objects with markedly different spectra: radio-quiet Seyfert 1's and 2's, and radio-loud Seyfert 1's. The average radio-quiet Seyfert 1 spectrum is well-fitted by a power law continuum with the energy spectral index alpha approximately equals 0.9, a Compton reflection component corresponding to a approximately 2 pi covering solid angle, and ionized absorption. There is a high-energy cutoff in the incident power law continuum: the e-folding energy is E(sub c) approximately equals 0.6(sup +0.8 sub -0.3) MeV. The simplest model that describes this spectrum is Comptonization in a relativistic optically-thin thermal corona above the surface of an accretion disk. Radio-quiet Seyfert 2's show strong netural absorption, and there is an indication that their X-ray power laws are intrinsically harder. Finally, the radio-loud Seyfert spectrum has alpha approximately equals 0.7, moderate neutral absorption E(sub C) = 0.4(sup +0.7 sub -0.2) MeV, and no or little Compton reflection. This is incompatible with the radio-quiet Seyfert 1 spectrum, and probably indicating that the X-rays are beamed away from the accretion disk in these objects. The average spectra of Seyferts integrated over redshift with a power-law evolution can explain the hard X-ray spectrum of the cosmic background.

Zdziarski, Andrzej A.↗

Extended X-Ray Emission around Quasars at Intermediate Redshift

We compare the optical to soft X-ray spectral energy distribution (SED) of a sample of bright low-redshift (0.048 less than z less than 0.155), radio-quiet quasars, with a range of thermal models which have been proposed to explain the optical/UV/soft X-ray quasar emission: (a) optically thin emission from an ionized plasma, (b) optically thick emission from the innermost regions of an accretion disk in Schwarzschild and Kerr geometries. We presented ROSAT PSPC observations of these quasars in an earlier paper. Here our goals are to search for the signature of thermal emission in the quasar SED, and to investigate whether a single component is dominating at different frequencies. We find that isothermal optically thin plasma models can explain the observed soft X-ray color and the mean OUV color. However, they predict an ultraviolet (1325 Angstrom) luminosity a factor of 3 to 10 times lower than observed. Pure disk models, even in a Kerr geometry, do not have the necessary flexibility to account for the observed OUV and soft X-ray luminosities. Additional components are needed both in the optical and in the soft X-rays (e.g. a hot corona can explain the soft X-ray color). The most constrained modification of pure disk models, is the assumption of an underlying power law component extending from the infrared (3 micrometers) to the X-ray. This can explain both the OUV and soft X-ray colors and luminosities and does not exceed the 3 micrometers luminosity, where a contribution from hot dust is likely to be important. We also discuss the possibility that the observed soft X-ray color and luminosity are dominated by reflection from the ionized surface of the accretion disk. While modifications of both optically thin plasma models and pure disk models might account for the observed SED, we do not find any strong evidence that the OUV bump and soft X-ray emission are one and the same component. Likewise, we do not find any strong argument which definitely argues in favor of thermal models.

Fiore, Fabrizio↗