SigSpec

Last updated August 29, 2023

SigSpec (acronym of SIGnificance SPECtrum) is a statistical technique to provide the reliability of periodicities in a measured (noisy and not necessarily equidistant) time series.^[1] It relies on the amplitude spectrum obtained by the Discrete Fourier transform (DFT) and assigns a quantity called the spectral significance (frequently abbreviated by “sig”) to each amplitude. This quantity is a logarithmic measure of the probability that the given amplitude level would be seen in white noise, in the sense of a type I error. It represents the answer to the question, “What would be the chance to obtain an amplitude like the measured one or higher, if the analysed time series were random?”

SigSpec may be considered a formal extension to the Lomb-Scargle periodogram,^[2]^[3] appropriately incorporating a time series to be averaged to zero before applying the DFT, which is done in many practical applications. When a zero-mean corrected dataset has to be statistically compared to a random sample, the sample mean (rather than the population mean only) has to be zero.

Probability density function (pdf) of white noise in Fourier space

Considering a time series to be represented by a set of $K$ pairs $(t_{k},x_{k})$ , the amplitude pdf of white noise in Fourier space, depending on frequency and phase angle may be described in terms of three parameters, $\alpha _{0}$ , $\beta _{0}$ , $\theta _{0}$ , defining the “sampling profile”, according to

\tan 2\theta _{0}={\frac {K\sum _{k=0}^{K-1}\sin 2\omega t_{k}-2\left(\sum _{k=0}^{K-1}\cos \omega t_{k}\right)\left(\sum _{k=0}^{K-1}\sin \omega t_{k}\right)}{K\sum _{k=0}^{K-1}\cos 2\omega t_{k}-{\big (}\sum _{k=0}^{K-1}\cos \omega t_{k}{\big )}^{2}+{\big (}\sum _{k=0}^{K-1}\sin \omega t_{k}{\big )}^{2}}},

\alpha _{0}={\sqrt {{\frac {2}{K^{2}}}\left(K\sum _{k=0}^{K-1}\cos ^{2}\left(\omega t_{k}-\theta _{0}\right)-\left[\sum _{l=0}^{K-1}\cos \left(\omega t_{k}-\theta _{0}\right)\right]^{2}\right)}},

\beta _{0}={\sqrt {{\frac {2}{K^{2}}}\left(K\sum _{k=0}^{K-1}\sin ^{2}\left(\omega t_{k}-\theta _{0}\right)-\left[\sum _{l=0}^{K-1}\sin \left(\omega t_{k}-\theta _{0}\right)\right]^{2}\right)}}.

In terms of the phase angle in Fourier space, $\theta$ , with

\tan \theta ={\frac {\sum _{k=0}^{K-1}\sin \omega t_{k}}{\sum _{k=0}^{K-1}\cos \omega t_{k}}},

the probability density of amplitudes is given by

\phi (A)={\frac {KA\cdot \operatorname {sock} }{2<x^{2}>}}\exp \left(-{\frac {KA^{2}}{4<x^{2}>}}\cdot \operatorname {sock} \right),

where the sock function is defined by

\operatorname {sock} (\omega ,\theta )=\left[{\frac {\cos ^{2}\left(\theta -\theta _{0}\right)}{\alpha _{0}^{2}}}+{\frac {\sin ^{2}\left(\theta -\theta _{0}\right)}{\beta _{0}^{2}}}\right]

and $<x^{2}>$ denotes the variance of the dependent variable $x_{k}$ .

False-alarm probability and spectral significance

Integration of the pdf yields the false-alarm probability that white noise in the time domain produces an amplitude of at least $A$ ,

\Phi _{\operatorname {FA} }(A)=\exp \left(-{\frac {KA^{2}}{4<x^{2}>}}\cdot \operatorname {sock} \right).

The sig is defined as the negative logarithm of the false-alarm probability and evaluates to

\operatorname {sig} (A)={\frac {KA^{2}\log e}{4<x^{2}>}}\cdot \operatorname {sock} .

It returns the number of random time series one would have to examine to obtain one amplitude exceeding $A$ at the given frequency and phase.

Applications

SigSpec is primarily used in asteroseismology to identify variable stars and to classify stellar pulsation (see references below). The fact that this method incorporates the properties of the time-domain sampling appropriately makes it a valuable tool for typical astronomical measurements containing data gaps.

Related Research Articles

In classical mechanics, a harmonic oscillator is a system that, when displaced from its equilibrium position, experiences a restoring force F proportional to the displacement x:

<span class="mw-page-title-main">Rutherford scattering</span> Elastic scattering of charged particles by the Coulomb force

In particle physics, Rutherford scattering is the elastic scattering of charged particles by the Coulomb interaction. It is a physical phenomenon explained by Ernest Rutherford in 1911 that led to the development of the planetary Rutherford model of the atom and eventually the Bohr model. Rutherford scattering was first referred to as Coulomb scattering because it relies only upon the static electric (Coulomb) potential, and the minimum distance between particles is set entirely by this potential. The classical Rutherford scattering process of alpha particles against gold nuclei is an example of "elastic scattering" because neither the alpha particles nor the gold nuclei are internally excited. The Rutherford formula further neglects the recoil kinetic energy of the massive target nucleus.

In mechanics and physics, simple harmonic motion is a special type of periodic motion an object experiences due to a restoring force whose magnitude is directly proportional to the distance of the object from an equilibrium position and acts towards the equilibrium position. It results in an oscillation that is described by a sinusoid which continues indefinitely.

Kinematics is a subfield of physics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of mathematics. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

In geometry, a solid angle is a measure of the amount of the field of view from some particular point that a given object covers. That is, it is a measure of how large the object appears to an observer looking from that point. The point from which the object is viewed is called the apex of the solid angle, and the object is said to subtend its solid angle at that point.

Chebyshev filters are analog or digital filters that have a steeper roll-off than Butterworth filters, and have either passband ripple or stopband ripple. Chebyshev filters have the property that they minimize the error between the idealized and the actual filter characteristic over the operating frequency range of the filter, but they achieve this with ripples in the passband. This type of filter is named after Pafnuty Chebyshev because its mathematical characteristics are derived from Chebyshev polynomials. Type I Chebyshev filters are usually referred to as "Chebyshev filters", while type II filters are usually called "inverse Chebyshev filters". Because of the passband ripple inherent in Chebyshev filters, filters with a smoother response in the passband but a more irregular response in the stopband are preferred for certain applications.

<span class="mw-page-title-main">Rotating reference frame</span> Concept in classical mechanics

A rotating frame of reference is a special case of a non-inertial reference frame that is rotating relative to an inertial reference frame. An everyday example of a rotating reference frame is the surface of the Earth.

Limb darkening is an optical effect seen in stars and planets, where the central part of the disk appears brighter than the edge, or limb. Its understanding offered early solar astronomers an opportunity to construct models with such gradients. This encouraged the development of the theory of radiative transfer.

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

In physics and engineering, a phasor is a complex number representing a sinusoidal function whose amplitude, angular frequency, and initial phase are time-invariant. It is related to a more general concept called analytic representation, which decomposes a sinusoid into the product of a complex constant and a factor depending on time and frequency. The complex constant, which depends on amplitude and phase, is known as a phasor, or complex amplitude, and sinor or even complexor.

In fluid dynamics, dispersion of water waves generally refers to frequency dispersion, which means that waves of different wavelengths travel at different phase speeds. Water waves, in this context, are waves propagating on the water surface, with gravity and surface tension as the restoring forces. As a result, water with a free surface is generally considered to be a dispersive medium.

<span class="mw-page-title-main">Hopf bifurcation</span> Critical point where a periodic solution arises

In the mathematical theory of bifurcations, a Hopfbifurcation is a critical point where, as a parameter changes, a system's stability switches and a periodic solution arises. More accurately, it is a local bifurcation in which a fixed point of a dynamical system loses stability, as a pair of complex conjugate eigenvalues—of the linearization around the fixed point—crosses the complex plane imaginary axis as a parameter crosses a threshold value. Under reasonably generic assumptions about the dynamical system, the fixed point becomes a small-amplitude limit cycle as the parameter changes.

A pendulum is a body suspended from a fixed support so that it swings freely back and forth under the influence of gravity. When a pendulum is displaced sideways from its resting, equilibrium position, it is subject to a restoring force due to gravity that will accelerate it back toward the equilibrium position. When released, the restoring force acting on the pendulum's mass causes it to oscillate about the equilibrium position, swinging it back and forth. The mathematics of pendulums are in general quite complicated. Simplifying assumptions can be made, which in the case of a simple pendulum allow the equations of motion to be solved analytically for small-angle oscillations.

In mathematics, the axis–angle representation parameterizes a rotation in a three-dimensional Euclidean space by two quantities: a unit vector $e$ indicating the direction (geometry) of an axis of rotation, and an angle of rotation $θ$ describing the magnitude and sense of the rotation about the axis. Only two numbers, not three, are needed to define the direction of a unit vector $e$ rooted at the origin because the magnitude of $e$ is constrained. For example, the elevation and azimuth angles of $e$ suffice to locate it in any particular Cartesian coordinate frame.

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

The Kapitza–Dirac effect is a quantum mechanical effect consisting of the diffraction of matter by a standing wave of light. The effect was first predicted as the diffraction of electrons from a standing wave of light by Paul Dirac and Pyotr Kapitsa in 1933. The effect relies on the wave–particle duality of matter as stated by the de Broglie hypothesis in 1924.

Amplitude amplification is a technique in quantum computing which generalizes the idea behind Grover's search algorithm, and gives rise to a family of quantum algorithms. It was discovered by Gilles Brassard and Peter Høyer in 1997, and independently rediscovered by Lov Grover in 1998.

In optics, the Fraunhofer diffraction equation is used to model the diffraction of waves when the diffraction pattern is viewed at a long distance from the diffracting object, and also when it is viewed at the focal plane of an imaging lens.

In physics, and especially scattering theory, the momentum-transfer cross section is an effective scattering cross section useful for describing the average momentum transferred from a particle when it collides with a target. Essentially, it contains all the information about a scattering process necessary for calculating average momentum transfers but ignores other details about the scattering angle.

<span class="mw-page-title-main">Averaged Lagrangian</span>

In continuum mechanics, Whitham's averaged Lagrangian method – or in short Whitham's method – is used to study the Lagrangian dynamics of slowly-varying wave trains in an inhomogeneous (moving) medium. The method is applicable to both linear and non-linear systems. As a direct consequence of the averaging used in the method, wave action is a conserved property of the wave motion. In contrast, the wave energy is not necessarily conserved, due to the exchange of energy with the mean motion. However the total energy, the sum of the energies in the wave motion and the mean motion, will be conserved for a time-invariant Lagrangian. Further, the averaged Lagrangian has a strong relation to the dispersion relation of the system.

References

↑ P. Reegen (2007). "SigSpec - I. Frequency- and phase-resolved significance in Fourier space". Astronomy and Astrophysics . 467 (3): 1353–1371. arXiv: physics/0703160 . Bibcode:2007A&A...467.1353R. doi:10.1051/0004-6361:20066597. S2CID 15076973.
↑ N. R. Lomb (1976). "Least-squares frequency analysis of unequally spaced data". Astrophysics and Space Science. 39 (2): 447–462. Bibcode:1976Ap&SS..39..447L. doi:10.1007/BF00648343. S2CID 2671466.
↑ J. D. Scargle (1982). "Studies in astronomical time series analysis. II. Statistical aspects of spectral analysis of unevenly spaced data". The Astrophysical Journal. 263: 835–853. Bibcode:1982ApJ...263..835S. doi:10.1086/160554.

M. Breger; S. M. Rucinski; P. Reegen (2007). "The Pulsation of EE Camelopardalis". The Astronomical Journal. 134 (5): 1994–1998. arXiv: 0709.3393 . Bibcode:2007AJ....134.1994B. doi:10.1086/522795. S2CID 120843648.
M. Gruberbauer; K. Kolenberg; J. F. Rowe; D. Huber; J. M. Matthews; P. Reegen; R. Kuschnig; C. Cameron; T. Kallinger; W. W. Weiss; D. B. Guenther; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker (2007). "MOST photometry of the RRdLyrae variable AQLeo: two radial modes, 32 combination frequencies and beyond". Monthly Notices of the Royal Astronomical Society . 379 (4): 1498–1506. arXiv: 0705.4603 . Bibcode:2007MNRAS.379.1498G. doi:10.1111/j.1365-2966.2007.12042.x. S2CID 55678660.
M. Gruberbauer; H. Saio; D. Huber; T. Kallinger; W. W. Weiss; D. B. Guenther; R. Kuschnig; J. M. Matthews; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker (2008). "MOST photometry and modeling of the rapidly oscillating (roAp) star γ Equulei". Astronomy and Astrophysics. 480 (1): 223–232. arXiv: 0801.0863 . Bibcode:2008A&A...480..223G. doi:10.1051/0004-6361:20078830. S2CID 54726017.
D. B. Guenther; T. Kallinger; P. Reegen; W. W. Weiss; J. M. Matthews; R. Kuschnig; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker (2007). "Searching for p-modes in η Bootis & Procyon using MOST satellite data". Communications in Asteroseismology. 151: 5–25. Bibcode:2007CoAst.151....5G. doi: 10.1553/cia151s5 .
D. B. Guenther; T. Kallinger; K. Zwintz; W. W. Weiss; J. Tanner (2007). "Seismology of Pre-Main-Sequence Stars in NGC 6530" (PDF). The Astrophysical Journal. 671 (1): 581–591. Bibcode:2007ApJ...671..581G. doi:10.1086/522880. S2CID 54866017.
D. Huber; H. Saio; M. Gruberbauer; W. W. Weiss; J. F. Rowe; M. Hareter; T. Kallinger; P. Reegen; J. M. Matthews; R. Kuschnig; D. B. Guenther; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker (2008). "MOST photometry of the roAp star 10 Aquilae". Astronomy and Astrophysics. 483 (1): 239–248. arXiv: 0803.1721 . Bibcode:2008A&A...483..239H. doi:10.1051/0004-6361:20079220. S2CID 3032930.
T. Kallinger; D. B. Guenther; J. M. Matthews; W. W. Weiss; D. Huber; R. Kuschnig; A. F. J. Moffat; S. M. Rucinski; D. Sasselov (2008). "Nonradial p-modes in the G9.5 giant ε Ophiuchi? Pulsation model fits to MOST photometry". Astronomy and Astrophysics. 478 (2): 497–505. arXiv: 0711.0837 . Bibcode:2008A&A...478..497K. doi:10.1051/0004-6361:20078171. S2CID 18201762.
T. Kallinger; P. Reegen; W. W. Weiss (2008). "A heuristic derivation of the uncertainty for frequency determination in time series data". Astronomy and Astrophysics. 481 (2): 571–574. arXiv: 0801.0683 . Bibcode:2008A&A...481..571K. doi:10.1051/0004-6361:20077559. S2CID 18481860.
P. Reegen (2005). ""SigSpec - reliable computation of significance in Fourier space", in The A-Star Puzzle, Proceedings IAU Symp. 224, eds. J. Zverko, J. Ziznovsky, S.J. Adelman, W.W. Weiss". The A-Star Puzzle, Proceedings of IAU Symposium 224. Cambridge, UK: Cambridge University Press. pp. 791–798. ISBN 0-521-85018-5.
P. Reegen; M. Gruberbauer; L. Schneider; W. W. Weiss (2008). "Cinderella - Comparison of INDEpendent RELative Least-squares Amplitudes". Astronomy and Astrophysics. 484 (2): 601–608. arXiv: 0710.2963 . Bibcode:2008A&A...484..601R. doi:10.1051/0004-6361:20078855. S2CID 11390524.
C. Schoenaers; A. E. Lynas-Gray (2007). "A new slowly pulsating subdwarf-B star: HD 4539". Communications in Asteroseismology. 151: 67–76. Bibcode:2007CoAst.151...67S. doi: 10.1553/cia151s67 .
M. Zechmeister; M. Kuerster (2009). "The gemeralised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms". Astronomy and Astrophysics. 496 (2): 577–584. arXiv: 0901.2573 . Bibcode:2009A&A...496..577Z. doi:10.1051/0004-6361:200811296. S2CID 10408194.
K. Zwintz; T. Kallinger; D. B. Guenther; M. Gruberbauer; D. Huber; J. Rowe; R. Kuschnig; W. W. Weiss; J. M. Matthews; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker; M. P. Casey (2009). "MOST photometry of the enigmatic PMS pulsator HD 142666". Astronomy and Astrophysics. 494 (3): 1031–1040. arXiv: 0812.1960 . Bibcode:2009A&A...494.1031Z. doi:10.1051/0004-6361:200811116. S2CID 54503935.
K. Zwintz; M. Hareter; R. Kuschnig; P. J. Amado; N. Nesvacil; E. Rodriguez; D. Diaz-Fraile; W. W. Weiss; T. Pribulla; D. B. Guenther; J. M. Matthews; A. F. J. Moffat; S. M. Rucinski; D. Sasselov; G. A. H. Walker (2009). "MOST observations of the young open cluster NGC 2264". Astronomy and Astrophysics. 502: 1239–252. doi:10.1051/0004-6361/200911863. S2CID 123505620.

External links

Website with further information on SigSpec calculation, etc.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] P. Reegen (2007). "SigSpec - I. Frequency- and phase-resolved significance in Fourier space". Astronomy and Astrophysics . 467 (3): 1353–1371. arXiv: physics/0703160 . Bibcode:2007A&A...467.1353R. doi:10.1051/0004-6361:20066597. S2CID 15076973.

[2] N. R. Lomb (1976). "Least-squares frequency analysis of unequally spaced data". Astrophysics and Space Science. 39 (2): 447–462. Bibcode:1976Ap&SS..39..447L. doi:10.1007/BF00648343. S2CID 2671466.

[3] J. D. Scargle (1982). "Studies in astronomical time series analysis. II. Statistical aspects of spectral analysis of unevenly spaced data". The Astrophysical Journal. 263: 835–853. Bibcode:1982ApJ...263..835S. doi:10.1086/160554.

[1]

[2]

[3]