Advanced options

Below are examples of different commands, including their before and after plots to demonstrate the desired effects.

Changing the fractional width of the power excess

via –ew & –exwidth

Fractional amount to scale the width of the oscillations envelope by – which is normally calculated w.r.t. solar values.

Before

After

pysyd run --star 9542776 --numax 900

pysyd run --star 9542776 --numax 900 --ew 1.5

_static/examples/9542776_before.png
_static/examples/9542776_after.png

Mitigating known Kepler artefacts

via -k, –kc & –kepcorr

Remove the well-known Kepler short-cadence artefact that occurs at/near the long-cadence nyquist frequency (\(\sim 270 \mu \mathrm{Hz}\)) by simulating white noise

Before

After

pysyd run --star 8045442 --numax 550

pysyd run --star 8045442 --numax 550 --kc

_static/examples/8045442_before.png
_static/examples/8045442_after.png

Hard-wiring the lower/upper limits of the power excess

via –lp & –lowerp

Manually set the lower frequency bound (or limit) of the power excess, which is helpful in the following scenarios:

  1. the width of the power excess is wildly different from that estimated by the solar scaling relation

  2. artefact or strange (typically not astrophysical) feature is close to the power excess and cannot be removed otherwise

  3. power excess is near the nyquist frequency

Before

After

pysyd run --star 10731424 --numax 750

pysyd run --star 10731424 --numax 750 --lp 490

_static/examples/10731424_before.png
_static/examples/10731424_after.png

I’m not sure how I feel about this one

via –npeaks & –peaks

Change the number of peaks chosen in the autocorrelation function (ACF) - this is especially helpful for low S/N cases, where the spectrum is noisy and the ACF has many peaks close the expected spacing (FIX THIS)

Before

After

pysyd run --star 9455860

pysyd run --star 9455860 --npeaks 10

_static/examples/9455860_before.png
_static/examples/9455860_after.png

Provide estimate for numax and save some time

via –numax

Turns out that a majority of the scaling relations used in this software can be written in terms of numax and therefore with the single estimate, we can guess the rest of the parameters (and fairly well, at that!)

If the value of \(\rm \nu_{max}\) is known, this can be provided to bypass the first module and save some time. There are also other ways to go about doing this, please see our notebook tutorial that goes through these different ways.

Before

After

pysyd run --star 5791521

pysyd run --star 5791521 --numax 670

_static/examples/5791521_before.png
_static/examples/5791521_after.png

Setting different frequency limits for the

via –ux & –upperx

Set the upper frequency limit in the power spectrum when estimating \(\rm \nu_{max}\) before the main fitting routine. This is helpful if there are high frequency artefacts that the software latches on to.

Before

After

pysyd run --star 11769801

pysyd run --star 11769801 --ux 3500

_static/examples/11769801_before.png
_static/examples/11769801_after.png

Smooth the echelle diagram by using matplotlib’s built-in interpolator

via -i, –ie & –interpech

Smooth the echelle diagram output by turning on the (bilinear) interpolation, which is helpful for identifying ridges in low S/N cases

Before

After

pysyd run 3112889 --numax 871.52

pysyd run --star 3112889 --numax 871.52 --ie

_static/examples/3112889_before.png
_static/examples/3112889_after.png