diff --git a/.gitignore b/.gitignore index 0df680f..2c66ccd 100644 --- a/.gitignore +++ b/.gitignore @@ -1,3 +1,8 @@ .venv *.pyc -*.csv \ No newline at end of file +*.csv + +# Quarto generated output +*.quarto/ +*_book/ +*_site/ \ No newline at end of file diff --git a/CHANGELOG.md b/CHANGELOG.md new file mode 100644 index 0000000..13106b9 --- /dev/null +++ b/CHANGELOG.md @@ -0,0 +1,8 @@ +# Changelog + +## 0.1.1 2022-10-14 + +### Added +- Used [Quarto](https://quarto.org) for authoring + * Converted proposal written in docs to quarto markdown. [#2](https://github.com/Thaza-Kun/sarjana/issues/2) +- Data fetching from [FRBSTATS](https://www.herta-experiment.org/frbstats/) in [`notebooks/01-get-data.ipynb`](notebooks/01-get-data.ipynb) \ No newline at end of file diff --git a/pyproject.toml b/pyproject.toml index 01d3301..850f984 100644 --- a/pyproject.toml +++ b/pyproject.toml @@ -1,6 +1,6 @@ [tool.poetry] name = "sarjana-muda" -version = "0.1.0" +version = "0.1.1" description = "" authors = ["Thaza_Kun <61819672+Thaza-Kun@users.noreply.github.com>"] diff --git a/thesis/.gitignore b/thesis/.gitignore new file mode 100644 index 0000000..075b254 --- /dev/null +++ b/thesis/.gitignore @@ -0,0 +1 @@ +/.quarto/ diff --git a/thesis/_quarto.yml b/thesis/_quarto.yml new file mode 100644 index 0000000..60d5aad --- /dev/null +++ b/thesis/_quarto.yml @@ -0,0 +1,23 @@ +project: + type: book + preview: + port: 4000 + +book: + title: "Characteristic Study Of Selected Fast Radio Bursts (FRB) Transients" + author: "Murthadza Aznam" + # TODO To be changed when published + date: '2022-07-14' + chapters: + - index.qmd + - literature.qmd + - methodology.qmd + - references.qmd + +bibliography: references.bib + +format: + html: + theme: cosmo + pdf: + documentclass: scrreprt \ No newline at end of file diff --git a/thesis/index.qmd b/thesis/index.qmd new file mode 100644 index 0000000..8fb95ca --- /dev/null +++ b/thesis/index.qmd @@ -0,0 +1,37 @@ +# Introduction + +Fast radio bursts (FRBs) are bright (50 mJy – 100 Jy), milliseconds long pulses of radio signal that is extra-terrestrial in nature. +The frequency range of FRBs reported is between 400 MHz to 8 GHz [@petroff_fast_2019]. +The underlying process of its emission is still poorly understood. +However, its detection has rapidly increased ever since its identification in Lorimer et al., (2007). +By June 2022, @spanakis-misirlis_frbstats_2021 has reported 806 FRB events including 25 repeaters indexed from the Transient Name Server, the CHIME/FRB Collaboration and several other cited reports. + +## Instruments +One reason for the abundance of the detection of FRBs is that many telescope dedicated to FRB has been built within the last decade. +Listed here are several telescopes that have been known to provide open data for FRB transients: + +### CHIME/FRB +The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a stationary radio telescope located in the southern parts of British Columbia, Canada (49^∘^ 19^'^ 14^''^.52 N, 119^∘^ 37^'^ 25^''^.25 W). +Unlike conventional dish shaped radio telescope, CHIME consists of four 20m x 100m semi-cylindrical reflectors pointed skywards to achieve 8000m2¬ with a sensitivity between 400 to 800 MHz [@the_chimefrb_collaboration_chime_2018]. +This allows the telescope to map ≳200 square degrees of the sky continuously with three different output pipelines: CHIME/Cosmology, CHIME/FRB, and CHIME/Pulsar. +The data relevant to this study is from the CHIME/FRB pipeline. + +### Parkes +The Parkes radio telescope is a 64-meter steerable radio telescope located in the central-west region of New South Wales, Australia. +It was the telescope used by Lorimer for the serendipitous discovery of the FRB in @lorimer_bright_2007. +It has been given the name Murriyang in honour of Australian aboriginals in 2020. +However, the name “Parkes radio telescope” or “Parkes” remains in use within the literature. + +## Problem Statement +This proposal aims to answer questions such as: + +1. What relationship of its properties can be extracted from the population of known FRBs now that we have more than 800 observed transients? +2. Which FRBs exhibit interesting properties compared to the rest of the population? +3. What might be the reason behind the interesting properties of the FRBs in question (2)? + +## Objective +Now that there is an abundant number of FRB detected, and many more are expected to be detected, this proposal aims to: + +1. Study the distribution of FRBs to understand the relationship between its properties. +2. Study the properties of selected samples of FRB to gain insight on its source or emission mechanism. +3. Implement programming methods for the distribution study of FRBs and for the feature study of selected FRBs. \ No newline at end of file diff --git a/thesis/literature.qmd b/thesis/literature.qmd new file mode 100644 index 0000000..4383afd --- /dev/null +++ b/thesis/literature.qmd @@ -0,0 +1,66 @@ +# Literature Review +## Observed Properties +### Dispersion Measure +As radio signals propagate through space, it interacts with mediums it encounters and thus is dispersed. +As a consequence, its arrival time is directly dependent on its frequency with lower frequencies arriving later [@day_pinpointing_2022; @kulkarni_dispersion_2020]. +The measure of this dispersion – its dispersion measure (DM) – is a key observable quantity in studies of FRBs [@keane_fast_2016]. +It is generally accepted that the DM of a FRB transient is due to free electrons encountered, $n_e$, along the path, dl, with the relation, +$$ +\text{DM} = \int_0^d n_e \text{d}l. +$$ {#eq-DM} + +However, it is worth noting that ionized particles, plasma temperature, magnetic fields and its relative motion contribute to this measure even though these other contributions are very small and can be neglected [@kulkarni_dispersion_2020]. +The DM of a given transient is calculated using two observables: its arrival time, $t$, and its frequency, $\nu$; where the DM is related by the slope relation +$$ +\Delta t = a\frac{\text{DM}}{\Delta\nu^2}, +$$ {#eq-DM-slope-a} +or more directly +$$ +\text{DM} = K\Delta\nu^2\Delta t, +$$ {#eq-DM-slope-k} +where $a=K^{-1}$ and $\Delta\nu^2$ is the difference of the highest and lowest frequency. +Users of calculation software should be made aware which conventions are used. +For the sake of precision, both values with its uncertainty are provided here [@kulkarni_dispersion_2020]: + + +$a=4.148 806 4239(11) \text{GH}^2\text{pc}^{-1}\text{cm}^3\text{ms}$ + +$K=241.033 1786(66) \text{GH}^{-2}\text{pc}\text{cm}^{-3}\text{s}$ + +### Fluence And Peak Flux Density +The fluence, $\mathcal{F}$, of a transient is the total energy received by the antenna throughout the duration of the signal typically measured in Jansky seconds, $\text{J}\cdot\text{s}$. +It then is characterized by the expression [@petroff_fast_2019] +$$ +\mathcal{F} = \int_\text{pulse} S(t) \text{d}t, +$$ {#eq-fluence-by-integral} +where $S(t)$ is the flux density of the signal. We can then characterize the peak flux density, $S_\text{peak}$ and its pulse width, $W_\text{pulse}$ like so +$$ +\mathcal{F} = S_\text{peak} W_\text{peak}. +$$ {#eq-fluence-by-peaks} + +### Rotation Measure +Analogous to dispersion measure, where the frequency is dispersed by the interacting medium, the radio wave may also be rotated as electrons interact with its magnetic field component. +This rotation is measurable via its polarization, $\Theta$, which is proportional to its wavelength squared, $\lambda^2$, +$$ +\Theta = \text{RM} \lambda^2, +$$ {#eq-RM} +where RM is the Faraday rotation measure. +This quantity, RM, is analogous to @eq-DM for DM as it is the total interaction along the line of sight [@brentjens_faraday_2005; @feng_frequency-dependent_2022], +$$ +\text{RM} = -0.81\int_0^d B(l)_\parallel n_e(l) \text{d}l. +$$ {#eq-RM-intergral} + +## Derived Properties + +### Distance Constraints +The DM of the FRBs can be used to infer the distance of the source by estimating contributions along the line of sight. +FRBs typically have DM more than the contributions from Milky Way, $\text{DM}_\text{MW}$, suggesting it is an extragalactic phenomenon. +The Milky Way contribution can be calculated using an electron density model associated with @eq-DM. +The YMW16 model in @yao_new_2017 uses known distances of pulsars using independent methods and matching it with their respective DMs. +Its value highly depends on its galactic latitude as there is little material higher up in the latitude so its DM is expected to be small [@thornton_population_2013]. +Non-Milky Way contributions to the observed DM (dubbed dispersion measure excess, $\text{DM}_\text{E}$) come from the intergalactic medium (IGM) and the host galaxy [@deng_cosmological_2014; @petroff_fast_2019], +$$ +\text{DM}_\text{E} = \text{DM} =\text{DM}_\text{MW} + \left(\frac{\text{DM}_\text{Host}}{1+z}\right), +$$ {#eq-DM-excess} +where $z$ is the redshift associated with the host galaxy. +The $\text{DM}_\text{Host}$ of the transient depends its distance from the host galactic core and viewing angle [@thornton_population_2013] because those factor determines how much material interacts with the signal throughout its propagation. diff --git a/thesis/methodology.qmd b/thesis/methodology.qmd new file mode 100644 index 0000000..e476dd0 --- /dev/null +++ b/thesis/methodology.qmd @@ -0,0 +1,49 @@ +# Methodology +## Data Source + +Listed below are open data on FRBs released by their respective telescopes that might be helpful for the research: +Telescope Server Data + +- CHIME/FRB + * Server: CANFAR + - [doi:10.11570/19.0004](doi:10.11570/19.0004) + - [doi:10.11570/19.0005](doi:10.11570/19.0005) + - [doi:10.11570/20.0002](doi:10.11570/20.0002) + - [doi:10.11570/20.0006](doi:10.11570/20.0006) + - [doi:10.11570/20.0006](doi:10.11570/20.0006) + +- Parkes + * Sever: gDCMP + - [https://data-portal.hpc.swin.edu.au/dataset/parkes-frbs-archival-data](https://data-portal.hpc.swin.edu.au/dataset/parkes-frbs-archival-data) + - [https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-high-time-resolution-universe-survey-high-latitude](https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-high-time-resolution-universe-survey-high-latitude) + - [https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-frb-140514](https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-frb-140514) + - [https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-frb-150215](https://data-portal.hpc.swin.edu.au/dataset/fast-radio-burst-data-frb-150215) + +- Lovell + * Server: Zenodo + - [https://zenodo.org/record/3974768#.YGWzqK8zap0](https://zenodo.org/record/3974768#.YGWzqK8zap0) + +- STARE2 + * Server: CalTechDATA + - [https://data.caltech.edu/records/1647](https://data.caltech.edu/records/1647) + + + +Other than that, the FRBSTATS^[@spanakis-misirlis_frbstats_2021 [https://www.herta-experiment.org/frbstats/](https://www.herta-experiment.org/frbstats/)] website provides an open-source live update of known FRBs, its properties, and its references. This will be helpful for statistical analysis for the entire population. It can be accessed programmatically via an Application Programming Interface (API) or downloaded as a comma separated value (.csv) or excel (.xlsx) files. + +## Data Analysis +The study will be using Python3.8 or higher to conduct the analysis of the data. It has good support for `.hd5`, `.msgpack`, `.fits` and other common file formats used in astronomy and statistics. In addition to that, the CHIME/FRB Open Data package that aims to assist the analysis of CHIME/FRB datasets, and the ‘astropy’ package (a reliable package for astronomy) is also written in Python. + +## Expected Outcome +The study expects to output the distribution of selected properties mentioned throughout the proposal and pick out important characteristics from the distribution. + +## Timeline +This study can be divided roughly into four phases and is expected to be 2 years long: + +1. Distribution study. Phase 1 will consist of a statistical study to the whole population of known FRBs. This is done to get a big picture of the landscape of FRBs. +2. Selection of FRBs. Phase 2 is the selection phase where several FRBs with interesting properties will be chosen. The selection phase includes the gathering of observation data on the selected FRB. +3. Feature study. Phase 3 will study the chosen FRBs for features that contribute to the interesting properties mentioned in phase 2. This is expected to be longest of the four phases. +4. Thesis writing and defence. + + + \ No newline at end of file diff --git a/thesis/references.bib b/thesis/references.bib new file mode 100644 index 0000000..f26c42c --- /dev/null +++ b/thesis/references.bib @@ -0,0 +1,174 @@ +@article{brentjens_faraday_2005, + title = {Faraday rotation measure synthesis}, + author = {Brentjens, M. A. and de Bruyn, A. G.}, + year = 2005, + month = oct, + journal = {Astronomy \& Astrophysics}, + volume = 441, + number = 3, + pages = {1217--1228}, + doi = {10.1051/0004-6361:20052990}, + issn = {0004-6361, 1432-0746}, + url = {http://www.aanda.org/10.1051/0004-6361:20052990}, + urldate = {2022-10-14}, + language = {en} +} +@article{day_pinpointing_2022, + title = {Pinpointing the {Origins} of {Fast} {Radios} {Bursts}}, + author = {Day, Cherie}, + year = 2022, + month = jan, + pages = 221, + language = {en}, +} + +@article{deng_cosmological_2014, + title = {Cosmological implications of {Fast} {Radio} {Burst} / {Gamma}-{Ray} {Burst} {Associations}}, + author = {Deng, Wei and Zhang, Bing}, + year = 2014, + month = feb, + journal = {The Astrophysical Journal}, + volume = 783, + number = 2, + pages = {L35}, + doi = {10.1088/2041-8205/783/2/L35}, + issn = {2041-8205, 2041-8213}, + url = {http://arxiv.org/abs/1401.0059}, + urldate = {2022-10-14}, + note = {arXiv:1401.0059 [astro-ph]}, + language = {en} +} +@article{feng_frequency-dependent_2022, + title = {Frequency-dependent polarization of repeating fast radio bursts-implications for their origin}, + author = {Feng, Yi and Li, Di and Yang, Yuan-Pei and Zhang, Yongkun and Zhu, Weiwei and Zhang, Bing and Lu, Wenbin and Wang, Pei and Dai, Shi and Lynch, Ryan S. and Yao, Jumei and Jiang, Jinchen and Niu, Jiarui and Zhou, Dejiang and Xu, Heng and Miao, Chenchen and Niu, Chenhui and Meng, Lingqi and Qian, Lei and Tsai, Chao-Wei and Wang, Bojun and Xue, Mengyao and Yue, Youling and Yuan, Mao and Zhang, Songbo and Zhang, Lei}, + year = 2022, + month = mar, + journal = {Science}, + volume = 375, + number = 6586, + pages = {1266--1270}, + doi = {10.1126/science.abl7759}, + issn = {0036-8075, 1095-9203}, + url = {http://arxiv.org/abs/2202.09601}, + urldate = {2022-10-14}, + note = {arXiv:2202.09601 [astro-ph]}, + language = {en} +} +@article{keane_fast_2016, + title = {A {Fast} {Radio} {Burst} {Host} {Galaxy}}, + author = {Keane, E. F. and Johnston, S. and Bhandari, S. and Barr, E. and Bhat, N. D. R. and Burgay, M. and Caleb, M. and Flynn, C. and Jameson, A. and Kramer, M. and Petroff, E. and Possenti, A. and van Straten, W. and Bailes, M. and Burke-Spolaor, S. and Eatough, R. P. and Stappers, B. W. and Totani, T. and Honma, M. and Furusawa, H. and Hattori, T. and Morokuma, T. and Niino, Y. and Sugai, H. and Terai, T. and Tominaga, N. and Yamasaki, S. and Yasuda, N. and Allen, R. and Cooke, J. and Jencson, J. and Kasliwal, M. M. and Kaplan, D. L. and Tingay, S. J. and Williams, A. and Wayth, R. and Chandra, P. and Perrodin, D. and Berezina, M. and Mickaliger, M. and Bassa, C.}, + year = 2016, + month = feb, + journal = {Nature}, + volume = 530, + number = 7591, + pages = {453--456}, + doi = {10.1038/nature17140}, + issn = {0028-0836, 1476-4687}, + url = {http://arxiv.org/abs/1602.07477}, + urldate = {2022-10-14}, + note = {arXiv:1602.07477 [astro-ph]}, + language = {en} +} +@misc{kulkarni_dispersion_2020, + title = {Dispersion measure: {Confusion}, {Constants} \& {Clarity}}, + shorttitle = {Dispersion measure}, + author = {Kulkarni, S. R.}, + year = 2020, + month = jul, + publisher = {arXiv}, + url = {http://arxiv.org/abs/2007.02886}, + urldate = {2022-10-14}, + note = {arXiv:2007.02886 [astro-ph]}, + language = {en} +} +@article{lorimer_bright_2007, + title = {A bright millisecond radio burst of extragalactic origin}, + author = {Lorimer, D. R. and Bailes, M. and McLaughlin, M. A. and Narkevic, D. J. and Crawford, F.}, + year = 2007, + month = nov, + journal = {Science}, + volume = 318, + number = 5851, + pages = {777--780}, + doi = {10.1126/science.1147532}, + issn = {0036-8075, 1095-9203}, + url = {http://arxiv.org/abs/0709.4301}, + urldate = {2022-10-14}, + note = {arXiv:0709.4301 [astro-ph]}, + language = {en} +} +@article{petroff_fast_2019, + title = {Fast radio bursts}, + author = {Petroff, E. and Hessels, J. W. T. and Lorimer, D. R.}, + year = 2019, + month = dec, + journal = {The Astronomy and Astrophysics Review}, + volume = 27, + number = 1, + pages = 4, + doi = {10.1007/s00159-019-0116-6}, + issn = {0935-4956, 1432-0754}, + url = {http://link.springer.com/10.1007/s00159-019-0116-6}, + urldate = {2022-10-14}, + language = {en} +} +@article{spanakis-misirlis_frbstats_2021, + title = {{FRBSTATS}: {A} web-based platform for visualization of fast radio burst properties}, + shorttitle = {{FRBSTATS}}, + author = {Spanakis-Misirlis, Apostolos}, + year = 2021, + month = jun, + journal = {Astrophysics Source Code Library}, + pages = {ascl:2106.028}, + url = {https://ui.adsabs.harvard.edu/abs/2021ascl.soft06028S}, + urldate = {2022-10-14}, + note = {ADS Bibcode: 2021ascl.soft06028S}, +} +@article{the_chimefrb_collaboration_chime_2018, + title = {The {CHIME} {Fast} {Radio} {Burst} {Project}: {System} {Overview}}, + shorttitle = {The {CHIME} {Fast} {Radio} {Burst} {Project}}, + author = {{The CHIME/FRB Collaboration} and Amiri, M. and Bandura, K. and Berger, P. and Bhardwaj, M. and Boyce, M. M. and Boyle, P. J. and Brar, C. and Burhanpurkar, M. and Chawla, P. and Chowdhury, J. and Cliche, J. F. and Cranmer, M. D. and Cubranic, D. and Deng, M. and Denman, N. and Dobbs, M. and Fandino, M. and Fonseca, E. and Gaensler, B. M. and Giri, U. and Gilbert, A. J. and Good, D. C. and Guliani, S. and Halpern, M. and Hinshaw, G. and Hofer, C. and Josephy, A. and Kaspi, V. M. and Landecker, T. L. and Lang, D. and Liao, H. and Masui, K. W. and Mena-Parra, J. and Naidu, A. and Newburgh, L. B. and Ng, C. and Patel, C. and Pen, U.-L. and Pinsonneault-Marotte, T. and Pleunis, Z. and Ravandi, M. Rafiei and Ransom, S. M. and Renard, A. and Scholz, P. and Sigurdson, K. and Siegel, S. R. and Smith, K. M. and Stairs, I. H. and Tendulkar, S. P. and Vanderlinde, K. and Wiebe, D. V.}, + year = 2018, + month = aug, + journal = {The Astrophysical Journal}, + volume = 863, + number = 1, + pages = 48, + doi = {10.3847/1538-4357/aad188}, + issn = {1538-4357}, + url = {http://arxiv.org/abs/1803.11235}, + urldate = {2022-10-10}, + note = {arXiv:1803.11235 [astro-ph]}, + language = {en} +} +@article{thornton_population_2013, + title = {A {Population} of {Fast} {Radio} {Bursts} at {Cosmological} {Distances}}, + author = {Thornton, D. and Stappers, B. and Bailes, M. and Barsdell, B. and Bates, S. and Bhat, N. D. R. and Burgay, M. and Burke-Spolaor, S. and Champion, D. J. and Coster, P. and D'Amico, N. and Jameson, A. and Johnston, S. and Keith, M. and Kramer, M. and Levin, L. and Milia, S. and Ng, C. and Possenti, A. and van Straten, W.}, + year = 2013, + month = jul, + journal = {Science}, + volume = 341, + number = 6141, + pages = {53--56}, + doi = {10.1126/science.1236789}, + issn = {0036-8075, 1095-9203}, + url = {https://www.science.org/doi/10.1126/science.1236789}, + urldate = {2022-10-14}, + language = {en} +} +@article{yao_new_2017, + title = {A {NEW} {ELECTRON}-{DENSITY} {MODEL} {FOR} {ESTIMATION} {OF} {PULSAR} {AND} {FRB} {DISTANCES}}, + author = {Yao, J. M. and Manchester, R. N. and Wang, N.}, + year = 2017, + month = jan, + journal = {The Astrophysical Journal}, + volume = 835, + number = 1, + pages = 29, + doi = {10.3847/1538-4357/835/1/29}, + issn = {1538-4357}, + url = {https://iopscience.iop.org/article/10.3847/1538-4357/835/1/29}, + urldate = {2022-10-14}, + language = {en} +} diff --git a/thesis/references.qmd b/thesis/references.qmd new file mode 100644 index 0000000..fd6da42 --- /dev/null +++ b/thesis/references.qmd @@ -0,0 +1,4 @@ +# References {.unnumbered} + +:::{#refs} +::: \ No newline at end of file