SigSpec

Last updated October 04, 2024

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: $where k is a positive constant.$

In optics, polarized light can be described using the Jones calculus, invented by R. C. Jones in 1941. Polarized light is represented by a Jones vector, and linear optical elements are represented by Jones matrices. When light crosses an optical element the resulting polarization of the emerging light is found by taking the product of the Jones matrix of the optical element and the Jones vector of the incident light. Note that Jones calculus is only applicable to light that is already fully polarized. Light which is randomly polarized, partially polarized, or incoherent must be treated using Mueller calculus.

In mechanics and physics, simple harmonic motion is a special type of periodic motion an object experiences by means of 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.

<span class="mw-page-title-main">Bremsstrahlung</span> Electromagnetic radiation due to deceleration of charged particles

In particle physics, bremsstrahlung is electromagnetic radiation produced by the deceleration of a charged particle when deflected by another charged particle, typically an electron by an atomic nucleus. The moving particle loses kinetic energy, which is converted into radiation, thus satisfying the law of conservation of energy. The term is also used to refer to the process of producing the radiation. Bremsstrahlung has a continuous spectrum, which becomes more intense and whose peak intensity shifts toward higher frequencies as the change of the energy of the decelerated particles increases.

In control theory and signal processing, a linear, time-invariant system is said to be minimum-phase if the system and its inverse are causal and stable.

In physics, the Rabi cycle is the cyclic behaviour of a two-level quantum system in the presence of an oscillatory driving field. A great variety of physical processes belonging to the areas of quantum computing, condensed matter, atomic and molecular physics, and nuclear and particle physics can be conveniently studied in terms of two-level quantum mechanical systems, and exhibit Rabi flopping when coupled to an optical driving field. The effect is important in quantum optics, magnetic resonance and quantum computing, and is named after Isidor Isaac Rabi.

In mathematics and signal processing, an analytic signal is a complex-valued function that has no negative frequency components. The real and imaginary parts of an analytic signal are real-valued functions related to each other by the Hilbert transform.

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, and initial phase are time-invariant and whose angular frequency is fixed. 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.

The Kerr–Newman metric describes the spacetime geometry around a mass which is electrically charged and rotating. It is a vacuum solution which generalizes the Kerr metric by additionally taking into account the energy of an electromagnetic field, making it the most general asymptotically flat and stationary solution of the Einstein–Maxwell equations in general relativity. As an electrovacuum solution, it only includes those charges associated with the magnetic field; it does not include any free electric charges.

The Kuramoto model, first proposed by Yoshiki Kuramoto, is a mathematical model used in describing synchronization. More specifically, it is a model for the behavior of a large set of coupled oscillators. Its formulation was motivated by the behavior of systems of chemical and biological oscillators, and it has found widespread applications in areas such as neuroscience and oscillating flame dynamics. Kuramoto was quite surprised when the behavior of some physical systems, namely coupled arrays of Josephson junctions, followed his model.

<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 such that it freely swings 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 towards 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.

Sinusoidal plane-wave solutions are particular solutions to the wave equation.

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.

The direct-quadrature-zerotransformation or zero-direct-quadraturetransformation is a tensor that rotates the reference frame of a three-element vector or a three-by-three element matrix in an effort to simplify analysis. The DQZ transform is the product of the Clarke transform and the Park transform, first proposed in 1929 by Robert H. Park.

Quantum mechanics was first applied to optics, and interference in particular, by Paul Dirac. Richard Feynman, in his Lectures on Physics, uses Dirac's notation to describe thought experiments on double-slit interference of electrons. Feynman's approach was extended to $N$ -slit interferometers for either single-photon illumination, or narrow-linewidth laser illumination, that is, illumination by indistinguishable photons, by Frank Duarte. The $N$ -slit interferometer was first applied in the generation and measurement of complex interference patterns.

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 (1): 1239–252. Bibcode:2009A&A...502..239Z. 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]