Showing 1–50 of 117 results for author: Mertens, F

Ionospheric contributions to the excess power in high-redshift 21-cm power-spectrum observations with LOFAR

Authors: S. A. Brackenhoff, M. Mevius, L. V. E. Koopmans, A. Offringa, E. Ceccotti, J. K. Chege, B. K. Gehlot, S. Ghosh, C. Höfer, F. G. Mertens, S. Munshi, S. Zaroubi

Abstract: The turbulent ionosphere causes phase shifts to incoming radio waves on a broad range of temporal and spatial scales. When an interferometer is not sufficiently calibrated for the direction-dependent ionospheric effects, the time-varying phase shifts can cause the signal to decorrelate. The ionosphere's influence over various spatiotemporal scales introduces a baseline-dependent effect on the inte… ▽ More The turbulent ionosphere causes phase shifts to incoming radio waves on a broad range of temporal and spatial scales. When an interferometer is not sufficiently calibrated for the direction-dependent ionospheric effects, the time-varying phase shifts can cause the signal to decorrelate. The ionosphere's influence over various spatiotemporal scales introduces a baseline-dependent effect on the interferometric array. We study the impact of baseline-dependent decorrelation on high-redshift observations with the Low Frequency Array (LOFAR). Datasets with a range of ionospheric corruptions are simulated using a thin-screen ionosphere model, and calibrated using the state-of-the-art LOFAR Epoch of Reionisation pipeline. For the first time ever, we show the ionospheric impact on various stages of the calibration process including an analysis of the transfer of gain errors from longer to shorter baselines using realistic end-to-end simulations. We find that direction-dependent calibration for source subtraction leaves excess power of up to two orders of magnitude above the thermal noise at the largest spectral scales in the cylindrically averaged auto-power spectrum under normal ionospheric conditions. However, we demonstrate that this excess power can be removed through Gaussian process regression, leaving no excess power above the ten per cent level for a $5~$km diffractive scale. We conclude that ionospheric errors, in the absence of interactions with other aggravating effects, do not constitute a dominant component in the excess power observed in LOFAR Epoch of Reionisation observations of the North Celestial Pole. Future work should therefore focus on less spectrally smooth effects, such as beam modelling errors. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 22 pages, 20 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The impact of lossy data compression on the power spectrum of the high redshift 21-cm signal with LOFAR

Authors: J. K. Chege, L. V. E. Koopmans, A. R. Offringa, B. K. Gehlot, S. A. Brackenhoff, E. Ceccotti, S. Ghosh, C. Höfer, F. G. Mertens, M. Mevius, S. Munshi

Abstract: Current radio interferometers output multi-petabyte-scale volumes of data per year making the storage, transfer, and processing of this data a sizeable challenge. This challenge is expected to grow with the next-generation telescopes such as the Square Kilometre Array. Lossy compression of interferometric data post-correlation can abate this challenge. However, since high-redshift 21-cm studies im… ▽ More Current radio interferometers output multi-petabyte-scale volumes of data per year making the storage, transfer, and processing of this data a sizeable challenge. This challenge is expected to grow with the next-generation telescopes such as the Square Kilometre Array. Lossy compression of interferometric data post-correlation can abate this challenge. However, since high-redshift 21-cm studies impose strict precision requirements, the impact of such lossy data compression on the 21-cm signal power spectrum statistic should be understood. We apply Dysco visibility compression, a technique to normalize and quantize specifically designed for radio interferometric data. We establish the level of the compression noise in the power spectrum in comparison to the thermal noise as well as its coherency behavior. Finally, for optimal compression results, we compare the compression noise obtained from different compression settings to a nominal 21-cm signal power. From a single night of observation, we find that the noise introduced due to the compression is more than five orders of magnitude lower than the thermal noise level in the power spectrum. The noise does not affect calibration. The compression noise shows no correlation with the sky signal and has no measurable coherent component. The level of compression error in the power spectrum ultimately depends on the compression settings. Dysco visibility compression is found to be of insignificant concern for 21-cm power spectrum studies. Hence, data volumes can be safely reduced by factors of $\sim 4$ and with insignificant bias to the final power spectrum. Data from SKA-low will likely be compressible by the same factor as LOFAR, owing to the similarities of the two instruments. The same technique can be used to compress data from other telescopes, but a small adjustment of the compression parameters might be required. △ Less

Submitted 16 July, 2024; originally announced July 2024.

Comments: 12 pages, 9 figures, and 3 tables; submitted to Astronomy and Astrophysics (A&A)

Beyond the horizon: Quantifying the full sky foreground wedge in the cylindrical power spectrum

Authors: S. Munshi, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, E. Ceccotti, S. A. Brackenhoff, J. K. Chege, B. K. Gehlot, S. Ghosh, C. Höfer, M. Mevius

Abstract: One of the main obstacles preventing the detection of the redshifted 21-cm signal from neutral hydrogen in the early Universe is the astrophysical foreground emission, which is several orders of magnitude brighter than the signal. The foregrounds, due to their smooth spectra, are expected to predominantly occupy a region in the cylindrical power spectrum known as the foreground wedge. However, the… ▽ More One of the main obstacles preventing the detection of the redshifted 21-cm signal from neutral hydrogen in the early Universe is the astrophysical foreground emission, which is several orders of magnitude brighter than the signal. The foregrounds, due to their smooth spectra, are expected to predominantly occupy a region in the cylindrical power spectrum known as the foreground wedge. However, the conventional equations describing the extent of the foreground wedge are derived under a flat-sky approximation. This assumption breaks down for tracking wide-field instruments, thus rendering these equations inapplicable in these situations. In this paper, we derive equations for the full sky foreground wedge and show that the foregrounds can potentially extend far beyond what the conventional equations suggest. We also derive the equations that describe a specific bright source in the cylindrical power spectrum space. The validity of both sets of equations is tested against numerical simulations. Many current and upcoming interferometers (e.g., LOFAR, NenuFAR, MWA, SKA) are wide-field phase-tracking instruments. These equations give us new insights into the nature of foreground contamination in the cylindrical power spectra estimated using wide-field instruments. Additionally, they allow us to accurately associate features in the power spectrum to foregrounds or instrumental effects. The equations are also important for correctly selecting the "EoR window" for foreground avoidance analyses, and for planning 21-cm observations. In future analyses, it is recommended to use these updated horizon lines to indicate the foreground wedge in the cylindrical power spectrum accurately. The new equations for generating the updated wedge lines are made available in a Python library, pslines. △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 21 pages, 13 figures, and 2 tables; submitted to Astronomy and Astrophysics (A&A)

Probing the intergalactic medium during the Epoch of Reionization using 21-cm signal power spectra

Authors: Raghunath Ghara, Abinash Kumar Shaw, Saleem Zaroubi, Benedetta Ciardi, Garrelt Mellema, Léon V. E. Koopmans, Anshuman Acharya, Madhurima Choudhury, Sambit K. Giri, Ilian T. Iliev, Qing-Bo Ma, Florent Mertens

Abstract: The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent fe… ▽ More The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21-cm signal power spectrum measurements. Under the assumption of high spin temperature, we modelled the redshift evolution of the ratio of EoR 21-cm brightness temperature power spectrum and the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. We have tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation $\lesssim 0.1$. Our framework allows us to interpret 21-cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework can efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 16 pages, 13 figures, 2 tables, Accepted for publication in Astronomy and Astrophysics

Report number: NORDITA 2024-009

Journal ref: A&A 687, A252 (2024)

YOLO-CIANNA: Galaxy detection with deep learning in radio data. I. A new YOLO-inspired source detection method applied to the SKAO SDC1

Authors: D. Cornu, P. Salomé, B. Semelin, A. Marchal, J. Freundlich, S. Aicardi, X. Lu, G. Sainton, F. Mertens, F. Combes, C. Tasse

Abstract: The upcoming Square Kilometer Array (SKA) will set a new standard regarding data volume generated by an astronomical instrument, which is likely to challenge widely adopted data analysis tools that scale inadequately with the data size. This study aims to develop a new source detection and characterization method for massive radio astronomical datasets by adapting modern deep-learning object detec… ▽ More The upcoming Square Kilometer Array (SKA) will set a new standard regarding data volume generated by an astronomical instrument, which is likely to challenge widely adopted data analysis tools that scale inadequately with the data size. This study aims to develop a new source detection and characterization method for massive radio astronomical datasets by adapting modern deep-learning object detection techniques. These approaches have proved their efficiency on complex computer vision tasks, and we seek to identify their specific strengths and weaknesses when applied to astronomical data. We introduce YOLO-CIANNA, a highly customized deep-learning object detector designed specifically for astronomical datasets. This paper presents the method and describes all the low-level adaptations required to address the specific challenges of radio-astronomical images. We demonstrate this method's capabilities using simulated 2D continuum images from the SKAO SDC1 dataset. Our method outperforms every other published result on the specific SDC1 dataset. Using the SDC1 metric, we improve the challenge-winning score by +139\% and the score of the only other post-challenge participation by +61\%. Our catalog has a detection purity of 94\% while detecting 40 to 60 \% more sources than previous top-score results. The trained model can also be forced to reach 99\% purity in post-process and still detect 10 to 30\% more sources than the other top-score methods. It is also capable of real-time detection, with a peak prediction speed of 500 images of 512x512 pixels per second on a single GPU. YOLO-CIANNA achieves state-of-the-art detection and characterization results on the simulated SDC1 dataset. The method is open source and included in the wider CIANNA framework. We provide scripts to train and apply this method to the SDC1 dataset in the CIANNA repository. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 40 pages

21-cm Signal from the Epoch of Reionization: A Machine Learning upgrade to Foreground Removal with Gaussian Process Regression

Authors: Anshuman Acharya, Florent Mertens, Benedetta Ciardi, Raghunath Ghara, Léon V. E. Koopmans, Sambit K. Giri, Ian Hothi, Qing-Bo Ma, Garrelt Mellema, Satyapan Munshi

Abstract: In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is… ▽ More In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is misestimated. To obtain better covariance models, we propose to use a kernel trained on the {\tt GRIZZLY} simulations using a Variational Auto-Encoder (VAE) based algorithm. In this work, we explore the abilities of this Machine Learning based kernel (VAE kernel) used with GPR, by testing it on mock signals from a variety of simulations, exploring noise levels corresponding to $\approx$10 nights ($\approx$141 hours) and $\approx$100 nights ($\approx$1410 hours) of observations with LOFAR. Our work suggests the possibility of successful extraction of the 21-cm signal within 2$σしぐま$ uncertainty in most cases using the VAE kernel, with better recovery of both shape and power than with previously used covariance models. We also explore the role of the excess noise component identified in past applications of GPR and additionally analyse the possibility of redshift dependence on the performance of the VAE kernel. The latter allows us to prepare for future LOFAR observations at a range of redshifts, as well as compare with results from other telescopes. △ Less

Submitted 29 November, 2023; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 13 pages, 7 figures, 3 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Report number: NORDITA 2023-074

Laser-induced Demagnetization in van der Waals $XY$- and Ising-like Antiferromagnets NiPS$_3$ and FePS$_3$

Authors: D. V. Kuntu, E. A. Arkhipova, L. A. Shelukhin, F. Mertens, M. A. Prosnikov, I. A. Eliseyev, A. N. Smirnov, V. Yu. Davydov, S. Mañas-Valero, E. Coronado, M. Cinchetti, A. M. Kalashnikova

Abstract: The critical behaviour of laser-induced changes in magnetic ordering is studied experimentally in two-dimensional zigzag antiferromagnets $XY$-like NiPS$_3$ and Ising-like FePS$_3$. To examine laser-induced dynamics in flakes of these compounds, we employ time-resolved exchange linear dichroism effect sensitive to zigzag magnetic ordering and independent of the orientation of the antiferromagnetic… ▽ More The critical behaviour of laser-induced changes in magnetic ordering is studied experimentally in two-dimensional zigzag antiferromagnets $XY$-like NiPS$_3$ and Ising-like FePS$_3$. To examine laser-induced dynamics in flakes of these compounds, we employ time-resolved exchange linear dichroism effect sensitive to zigzag magnetic ordering and independent of the orientation of the antiferromagnetic vector. In both compounds laser excitation in the vicinity of the absorption edge induces partial quenching of the antiferromagnetic ordering manifested by exchange linear dichroism reduction. The amplitude of the effect varies with temperature as the derivative of the antiferromagnetic vector and exhibits a critical behaviour with the exponents corresponding to $XY$- and Ising-models for NiPS$_3$ and FePS$_3$, respectively. Critical slowing down of the demagnetization in the vicinity of Néel temperature is found, however, only in FePS$_3$. In contrast, the increase of the demagnetization time near the ordering temperature in NiPS$_3$ is minor. We show that the difference in the demagnetization times correlates well with the spin specific heat in both compounds. Beyond the range of slowing down, the demagnetization times in NiPS$_3$ and FePS$_3$ are comparable, about 5 - 10 ps, and are longer than those reported earlier for CoPS$_3$ and considerably shorter than for MnPS$_3$. This points to the importance of the unquenched angular momentum of transition-metal ions in laser-induced demagnetization process. △ Less

Submitted 31 January, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. Materials 8, 014408 (2024)

First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR

Authors: S. Munshi, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, B. Semelin, D. Aubert, R. Barkana, A. Bracco, S. A. Brackenhoff, B. Cecconi, E. Ceccotti, S. Corbel, A. Fialkov, B. K. Gehlot, R. Ghara, J. N. Girard, J. M. Grießmeier, C. Höfer, I. Hothi, R. Mériot, M. Mevius, P. Ocvirk, A. K. Shaw, G. Theureau, S. Yatawatta , et al. (2 additional authors not shown)

Abstract: The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 2… ▽ More The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$σしぐま$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future. △ Less

Submitted 30 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Journal ref: A&A 681, A62 (2024)

Transient RFI environment of LOFAR-LBA at 72-75 MHz: Impact on ultra-widefield AARTFAAC Cosmic Explorer observations of the redshifted 21-cm signal

Authors: B. K. Gehlot, L. V. E. Koopmans, S. A. Brackenhoff, E. Ceccotti, S. Ghosh, C. Höfer, F. G. Mertens, M. Mevius, S. Munshi, A. R. Offringa, V. N. Pandey, A. Rowlinson, A. Shulevski, R. A. M. J. Wijers, S. Yatawatta, S. Zaroubi

Abstract: Measurement of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn (CD) and Epoch of Reionisation (EoR) promises to unveil a wealth of information about the astrophysical processes during the first billion years of evolution of the universe. The AARTFAAC Cosmic Explorer (ACE) utilises the AARTFAAC wide-field imager of LOFAR to measure the power spectrum of the intensity fluctuatio… ▽ More Measurement of the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn (CD) and Epoch of Reionisation (EoR) promises to unveil a wealth of information about the astrophysical processes during the first billion years of evolution of the universe. The AARTFAAC Cosmic Explorer (ACE) utilises the AARTFAAC wide-field imager of LOFAR to measure the power spectrum of the intensity fluctuations of the redshifted 21-cm signal from the CD at z~18. The RFI from various sources contaminates the observed data and it is crucial to exclude the RFI-affected data in the analysis for reliable detection. In this work, we investigate the impact of non-ground-based transient RFI using cross-power spectra and cross-coherence metrics to assess the correlation of RFI over time and investigate the level of impact of transient RFI on the ACE 21-cm power spectrum estimation. We detected moving sky-based transient RFI sources that cross the field of view within a few minutes and appear to be mainly from aeroplane communication beacons at the location of the LOFAR core in the 72-75 MHz band, by inspecting filtered images. This transient RFI is mostly uncorrelated over time and is only expected to dominate over the thermal noise for an extremely deep integration time of 3000 hours or more with a hypothetical instrument that is sky temperature dominated at 75 MHz. We find no visible correlation over different k-modes in Fourier space in the presence of noise for realistic thermal noise scenarios. We conclude that the sky-based transient RFI from aeroplanes, satellites and meteorites at present does not pose a significant concern for the ACE analyses at the current level of sensitivity and after integrating over the available 500 hours of observed data. However, it is crucial to mitigate or filter such transient RFI for more sensitive experiments aiming for significantly deeper integration. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 13 pages, 9 figures, and 3 tables; accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 681, A71 (2024)

A novel radio imaging method for physical spectral index modelling

Authors: E. Ceccotti, A. R. Offringa, L. V. E. Koopmans, R. Timmerman, S. A. Brackenhoff, B. K. Gehlot, F. G. Mertens, S. Munshi, V. N. Pandey, R. J. van Weeren, S. Yatawatta, S. Zaroubi

Abstract: We present a new method, called "forced-spectrum fitting", for physically-based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multi-frequency deconvolution of WS… ▽ More We present a new method, called "forced-spectrum fitting", for physically-based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multi-frequency deconvolution of WSClean, where the pre-determined spectrum is fitted. The component magnitude is evaluated by performing a modified weighted linear least-squares fit. We test this method on a simulated LOFAR-HBA observation of the 3C196 QSO and a real LOFAR-HBA observation of the 4C+55.16 FRI galaxy. We compare the results from the forced-spectrum fitting with traditional joined-channel deconvolution using polynomial fitting. Because no prior spectral information was available for 4C+55.16, we demonstrate a method for extracting spectral indices in the observed frequency band using "clustering". The models generated by the forced-spectrum fitting are used to improve the calibration of the datasets. The final residuals are comparable to existing multi-frequency deconvolution methods, but the output model agrees with the provided spectral index map, embedding correct spectral information. While forced-spectrum fitting does not solve the determination of the spectral information itself, it enables the construction of accurate multi-frequency models that can be used for wide-band calibration and subtraction. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 17 pages, 9 figures, 5 tables. Accepted for publication in MNRAS

Retrieving the 21-cm signal from the Epoch of Reionization with learnt Gaussian process kernels

Authors: Florent G. Mertens, Jérôme Bobin, Isabella P. Carucci

Abstract: Direct detection of the Cosmic Dawn and Epoch of Reionization via the redshifted 21-cm line of neutral Hydrogen will have unprecedented implications for studying structure formation in the early Universe. This exciting goal is challenged by the difficulty of extracting the faint 21-cm signal buried beneath bright astrophysical foregrounds and contaminated by numerous systematics. Here, we focus on… ▽ More Direct detection of the Cosmic Dawn and Epoch of Reionization via the redshifted 21-cm line of neutral Hydrogen will have unprecedented implications for studying structure formation in the early Universe. This exciting goal is challenged by the difficulty of extracting the faint 21-cm signal buried beneath bright astrophysical foregrounds and contaminated by numerous systematics. Here, we focus on improving the Gaussian Process Regression (GPR) signal separation method originally developed for LOFAR observations. We address a key limitation of the current approach by incorporating covariance prior models learnt from 21-cm signal simulations using Variational Autoencoder (VAE) and Interpolatory Autoencoder (IAE). Extensive tests are conducted to evaluate GPR, VAE-GPR, and IAE-GPR in different scenarios. Our findings reveal that the new method outperforms standard GPR in component separation tasks. Moreover, the improved method demonstrates robustness when applied to signals not represented in the training set. It also presents a certain degree of resilience to data systematics, highlighting its ability to effectively mitigate their impact on the signal recovery process. However, our findings also underscore the importance of accurately characterizing and understanding these systematics to achieve successful detection. Our generative approaches provide good results even with limited training data, offering a valuable advantage when a large training set is not feasible. Comparing the two algorithms, IAE-GPR shows slightly higher fidelity in recovering power spectra compared to VAE-GPR. These advancements highlight the strength of generative approaches and optimise the analysis techniques for future 21-cm signal detection at high redshifts. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 14 pages, 9 figues, submitted to MNRAS, any comments welcome

Deep learning approach for identification of HII regions during reionization in 21-cm observations -- II. foreground contamination

Authors: Michele Bianco, Sambit. K. Giri, David Prelogović, Tianyue Chen, Florent G. Mertens, Emma Tolley, Andrei Mesinger, Jean-Paul Kneib

Abstract: The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 fra… ▽ More The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 framework, an enhanced version of our convolutional neural network, built to identify neutral and ionized regions in the 21-cm signal contaminated with foreground emission. We trained our neural network on 21-cm image data processed by a foreground removal method based on Principal Component Analysis achieving an average classification accuracy of 71 per cent between redshift $z=7$ to $11$. We tested SegU-Net v2 against various foreground removal methods, including Gaussian Process Regression, Polynomial Fitting, and Foreground-Wedge Removal. Results show comparable performance, highlighting SegU-Net v2's independence on these pre-processing methods. Statistical analysis shows that a perfect classification score with $AUC=95\%$ is possible for $8<z<10$. While the network prediction lacks the ability to correctly identify ionized regions at higher redshift and differentiate well the few remaining neutral regions at lower redshift due to low contrast between 21-cm signal, noise and foreground residual in images. Moreover, as the photon sources driving reionization are expected to be located inside ionised regions, we show that SegU-Net v2 can be used to correctly identify and measure the volume of isolated bubbles with $V_{\rm ion}>(10\, {\rm cMpc})^3$ at $z>9$, for follow-up studies with infrared/optical telescopes to detect these sources. △ Less

Submitted 28 February, 2024; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 19 pages, 10 figures, 3 tables, 2 appendixes

Report number: NORDITA 2023-013

Studying the Multi-frequency Angular Power Spectrum of the Cosmic Dawn 21-cm Signal

Authors: Abinash Kumar Shaw, Raghunath Ghara, Saleem Zaroubi, Rajesh Mondal, Garrelt Mellema, Florent Mertens, Léon V. E. Koopmans, Benoît Semelin

Abstract: The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along… ▽ More The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along every direction in the signal volume. Unlike the 3D power spectrum which mixes the cosmic evolution of the 21-cm signal along the LoS $k$ modes, MAPS keeps the evolution information disentangled. Here we first study the impact of different underlying physical processes during cosmic dawn on the behaviour of the 21-cm MAPS using simulations of various different scenarios and models. We also make error predictions in 21-cm MAPS measurements considering only the system noise and cosmic variance for mock observations of HERA, NenuFAR and SKA-Low. We find that $100~{\rm h}$ of HERA observations will be able to measure 21-cm MAPS at $\geq 3σしぐま$ for $\ell \lesssim 1000$ with $0.1\,{\rm MHz}$ channel-width. The better sensitivity of SKA-Low allows reaching this sensitivity up to $\ell \lesssim 3000$. Note that due to the difference in the frequency coverage of the various experiements, the CD-EoR model considered for NenuFAR is different than those used for the HERA and SKA-Low predictions. Considering NenuFAR with the new model, measurements $\geq 2σしぐま$ are possible only for $\ell \lesssim 600$ with $0.2\,{\rm MHz}$ channel-width and for a ten times longer observation time of $t_{\rm obs} = 1000~{\rm h}$. However, for the range $300 \lesssim \ell \lesssim 600$ and $t_{\rm obs}=1000~{\rm h}$ more than $3σしぐま$ measurements are still possible for NenuFAR when combining consecutive frequency channels within a $5 ~{\rm MHz}$ band. △ Less

Submitted 13 April, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 21 pages, 12 figures, 3 tables, accepted for publication in MNRAS

Accurate modelling of the Lyman-$αあるふぁ$ coupling for the 21-cm signal, observability with NenuFAR and SKA

Authors: Benoit Semelin, Romain Mériot, Florent Mertens, Léon V. E. Koopmans, Dominique Aubert, Rennan Barkana, Anastasia Fialkov, Satyapan Munshi, Pierre Ocvirk

Abstract: The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluct… ▽ More The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$αあるふぁ$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$αあるふぁ$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$αあるふぁ$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 672, A162 (2023)

Assessing the impact of two independent direction-dependent calibration algorithms on the LOFAR 21-cm signal power spectrum

Authors: H. Gan, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, M. Mevius, V. N. Pandey, S. A. Brackenhoff, E. Ceccotti, B. Ciardi, B. K. Gehlot, R. Ghara, S. K. Giri, I. T. Iliev, S. Munshi

Abstract: Detecting the 21-cm signal from the Epoch of Reionisation (EoR) is challenging due to the strong astrophysical foregrounds, ionospheric effects, radio frequency interference and instrumental effects. Understanding and calibrating these effects are crucial for the detection. In this work, we introduce a newly developed direction-dependent (DD) calibration algorithm DDECAL and compare its performanc… ▽ More Detecting the 21-cm signal from the Epoch of Reionisation (EoR) is challenging due to the strong astrophysical foregrounds, ionospheric effects, radio frequency interference and instrumental effects. Understanding and calibrating these effects are crucial for the detection. In this work, we introduce a newly developed direction-dependent (DD) calibration algorithm DDECAL and compare its performance with an existing algorithm, SAGECAL, in the context of the LOFAR-EoR 21-cm power spectrum experiment. In our data set, the North Celestial Pole (NCP) and its flanking fields were observed simultaneously. We analyse the NCP and one of its flanking fields. The NCP field is calibrated by the standard pipeline, using SAGECAL with an extensive sky model and 122 directions, and the flanking field is calibrated by DDECAL and SAGECAL with a simpler sky model and 22 directions. Additionally, two strategies are used for subtracting Cassiopeia A and Cygnus A. The results show that DDECAL performs better at subtracting sources in the primary beam region due to the application of a beam model, while SAGECAL performs better at subtracting Cassiopeia A and Cygnus A. This indicates that including a beam model during DD calibration significantly improves the performance. The benefit is obvious in the primary beam region. We also compare the 21-cm power spectra on two different fields. The results show that the flanking field produces better upper limits compared to the NCP in this particular observation. Despite the minor differences between DDECAL and SAGECAL due to the beam application, we find that the two algorithms yield comparable 21-cm power spectra on the LOFAR-EoR data after foreground removal. Hence, the current LOFAR-EoR 21-cm power spectrum limits are not likely to depend on the DD calibration method. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 28 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 669, A20 (2023)

Statistical analysis of the causes of excess variance in the 21 cm signal power spectra obtained with the Low-Frequency Array

Authors: H. Gan, L. V. E Koopmans, F. G. Mertens, M. Mevius, A. R. Offringa, B. Ciardi, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, G. Mellema, V. N. Pandey, S. Zaroubi

Abstract: The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still abov… ▽ More The detection of the 21 cm signal of neutral hydrogen from the Epoch of Reionization (EoR) is challenging due to bright foreground sources, radio frequency interference (RFI), the ionosphere, and instrumental effects. Even after correcting for these effects in the calibration step and applying foreground removal techniques, the remaining residuals in the observed 21 cm power spectra are still above the thermal noise, which is referred to as the "excess variance." We study potential causes of this excess variance based on 13 nights of data obtained with the Low-Frequency Array (LOFAR). We focused on the impact of gain errors, the sky model, and ionospheric effects on the excess variance by correlating the relevant parameters such as the gain variance over time or frequency, local sidereal time (LST), diffractive scale, and phase structure-function slope with the level of excess variance. Our analysis shows that excess variance has an LST dependence, which is related to the power from the sky. And the simulated Stokes I power spectra from bright sources and the excess variance show a similar progression over LST with the minimum power appearing at LST bin 6h to 9h. This LST dependence is also present in sky images of the residual Stokes I of the observations. In very-wide sky images, we demonstrate that the extra power comes exactly from the direction of bright and distant sources Cassiopeia A and Cygnus A with the array beam patterns. These results suggest that the level of excess variance in the 21 cm signal power spectra is related to sky effects and, hence, it depends on LST. In particular, very bright and distant sources such as Cassiopeia A and Cygnus A can dominate the effect. This is in line with earlier studies and offers a path forward toward a solution since the correlation between the sky-related effects and the excess variance is non-negligible. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 29 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 663, A9 (2022)

Degree-Scale Galactic Radio Emission at 122 MHz around the North Celestial Pole with LOFAR-AARTFAAC

Authors: B. K. Gehlot, L. V. E. Koopmans, A. R. Offringa, H. Gan, R. Ghara, S. K. Giri, M. Kuiack, F. G. Mertens, M. Mevius, R. Mondal, V. N. Pandey, A. Shulevski, R. A. M. J. Wijers, S. Yatawatta

Abstract: Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galacti… ▽ More Aims: Contamination from bright diffuse Galactic thermal and non-thermal radio emission poses crucial challenges in experiments aiming to measure the 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization. If not included in calibration, this diffuse emission can severely impact the analysis and signal extraction in 21-cm experiments. We examine large-scale diffuse Galactic emission at 122~MHz, around the North Celestial Pole, using the AARTFAAC-HBA system. Methods: In this pilot project, we present the first-ever wide-field image produced with a single sub-band of the data recorded with the AARTFAAC-HBA. We demonstrate two methods: multiscale CLEAN and shapelet decomposition, to model the diffuse emission revealed in the image. We use angular power spectrum metrics to quantify different components of the emission and compare the performance of the two diffuse structure modelling approaches. Results: We observe that the point sources dominate the angular power spectrum ($\ell(\ell+1)C_{\ell}/2πぱい\equiv Δでるた^2(\ell)$) of the emission in the field on scales $\ell\gtrsim 60$ ($\lesssim 3$~degree). The angular power spectrum after subtraction of compact sources is flat within $20\lesssim \ell \lesssim200$ range, suggesting that the residual power is dominated by the diffuse emission on scales $\ell\lesssim200$. The residual diffuse emission has a brightness temperature variance of $Δでるた^2_{\ell=180} = (145.64 \pm 13.61)~{\rm K}^2$ at 122~MHz on angular scales of 1~degree, and is consistent with a power-law following $C_{\ell}\propto \ell^{-2.0}$ in $20\lesssim \ell \lesssim200$ range. We also find that, in the current setup, the multiscale CLEAN is suitable to model the compact and diffuse structures on a wide range of angular scales, whereas the shapelet decomposition method better models the large scales, which are of the order of a few degrees and wider. △ Less

Submitted 7 April, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 11 pages, 9 figures, accepted (A&A)

Journal ref: A&A 662, A97 (2022)

A numerical study of 21-cm signal suppression and noise increase in direction-dependent calibration of LOFAR data

Authors: M. Mevius, F. Mertens, L. V. E. Koopmans, A. R. Offringa, S. Yatawatta, M. A. Brentjens, E. Chapman, B. Ciardi, H. Gan, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, G. Mellema, V. N. Pandey, S. Zaroubi

Abstract: We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 λらむだ$. We show that suppression of faint signals can effectively be avoided by calibrating these short base… ▽ More We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 λらむだ$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

Journal ref: MNRAS 2021

Ultrafast coherent THz lattice dynamics coupled to spins in a van der Waals antiferromagnetic flake

Authors: Fabian Mertens, David Mönkebüscher, Eugenio Coronado, Samuel Mañas-Valero, Carla Boix-Constant, Alberta Bonanni, Margherita Matzer, Rajdeep Adhikari, Alexandra M. Kalashnikova, Davide Bossini, Mirko Cinchetti

Abstract: A coherent THz optical lattice mode is triggered by femtosecond laser pulses in the antiferromagnetic van der Waals semiconductor FePS$_3$. The 380 nm thick exfoliated flake was placed on a substrate and laser-driven lattice and spin dynamics were investigated as a function of the excitation photon energy and sample temperature. The pump-probe spectroscopic measurements reveal that the photo-induc… ▽ More A coherent THz optical lattice mode is triggered by femtosecond laser pulses in the antiferromagnetic van der Waals semiconductor FePS$_3$. The 380 nm thick exfoliated flake was placed on a substrate and laser-driven lattice and spin dynamics were investigated as a function of the excitation photon energy and sample temperature. The pump-probe spectroscopic measurements reveal that the photo-induced phonon is generated by a displacive mechanism. The amplitude of the phononic signal decreases as the sample is heated up to the Néel temperature and vanishes as the phase transition to the paramagnetic phase occurs. This evidence confirms that the excited lattice mode is intimately connected to the long-range magnetic order. Therefore our work discloses a pathway towards a femtosecond coherent manipulation of the magneto-crystalline anisotropy in a van der Waals antiferromagnet. In fact, it is reported that by applying a magnetic field the induced phonon mode hybridizes via the Kittel-mechanism with zone-centre magnons. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: The following article has been submitted to Applied Physics Letters. After it is published, it will be found at https://aip.scitation.org/journal/apl

SKA-Low Intensity Mapping Pathfinder Updates: Deeper 21 cm Power Spectrum Limits from Improved Analysis Frameworks

Authors: Nichole Barry, Gianni Bernardi, Bradley Greig, Nicholas Kern, Florent Mertens

Abstract: The Square Kilometre Array (SKA) is a planned radio interferometer of unprecedented scale that will revolutionize low-frequency radio astronomy when completed. In particular, one of its core science drivers is the systematic mapping of the Cosmic Dawn and Epoch of Reionization, which mark the birth of the first stars and galaxies in the Universe and their subsequent ionization of primordial interg… ▽ More The Square Kilometre Array (SKA) is a planned radio interferometer of unprecedented scale that will revolutionize low-frequency radio astronomy when completed. In particular, one of its core science drivers is the systematic mapping of the Cosmic Dawn and Epoch of Reionization, which mark the birth of the first stars and galaxies in the Universe and their subsequent ionization of primordial intergalactic hydrogen, respectively. The SKA will offer the most sensitive view of these poorly understood epochs using the redshifted 21 cm hyperfine signal from intergalactic hydrogen. However, significant technical challenges stand in the way of realizing this scientific promise. These mainly involve the mitigation of systematics coming from astrophysical foregrounds, terrestrial radio interference, and the instrumental response. The Low Frequency Array, the Murchison Widefield Array and the Hydrogen Epoch of Reionization Array are SKA pathfinder experiments that have developed a variety of strategies for addressing these challenges, each with unique characteristics that stem largely from their different instrumental designs. We outline these various directions, highlighting key differences and synergies, and discuss how these relate to the future of low-frequency intensity mapping with the SKA. We also briefly summarize the challenges associated with modeling the 21 cm signal and discuss the methodologies being proposed for inferring constraints on astrophysical models. △ Less

Submitted 9 January, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: Version accepted to JATIS

Exploring the Cosmic Dawn with NenuFAR

Authors: F. G. Mertens, B. Semelin, L. V. E. Koopmans

Abstract: The exploration of the Cosmic Dawn, the period of the Universe during which the first stars and galaxies were formed, is one of the last frontiers of modern astronomy and cosmology. The redshifted 21-cm line emission from neutral hydrogen is a unique probe and can open this era for astrophysical and cosmological studies. The tentative detection of the 21-cm global signal by the EDGES team at… ▽ More The exploration of the Cosmic Dawn, the period of the Universe during which the first stars and galaxies were formed, is one of the last frontiers of modern astronomy and cosmology. The redshifted 21-cm line emission from neutral hydrogen is a unique probe and can open this era for astrophysical and cosmological studies. The tentative detection of the 21-cm global signal by the EDGES team at $z \sim 17$ underlines the need for an interferometric detection of this signal to discriminate between the numerous models trying to explain this unexpected discovery. The NenuFAR Cosmic-Dawn Key-Science Program aims to perform this detection in the redshift range $z \sim 15 - 31$ with a novel SKA precursor, the NenuFAR radio telescope located at the Station de Radioastronomie de Nançay and that started operating in 2019. Due to its compactness it is particularly sensitive to the large scale of the 21-cm signal. Only 100 hours of observation are needed to reach the level of the most extreme models, while 1000 hours are needed for the more standard models. Observations have already started, accumulating to almost 500 hours on the North Celestial Pole field. In this contribution, we introduce the project, our first results and the developments in calibration and RFI mitigation specific to this new instrument. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 4 pages, SF2A 2021 proceeding

A detailed kinematic study of 3C 84 and its connection to Gamma-rays

Authors: Jeffrey A. Hodgson, Bindu Rani, Junghwan Oh, Alan Marscher, Svetlana Jorstad, Yosuke Mizuno, Jongho Park, Sang-Sung Lee, Sascha Trippe, Florent Mertens

Abstract: 3C 84 (NGC 1275) is the bright radio core of the Perseus Cluster. Even in the absence of strong relativistic effects, the source has been detected at Gamma-rays up to TeV energies. Despite its intensive study, the physical processes responsible for the high-energy emission in the source remain unanswered. We present a detailed kinematics study of the source and its connection to Gamma-ray emission… ▽ More 3C 84 (NGC 1275) is the bright radio core of the Perseus Cluster. Even in the absence of strong relativistic effects, the source has been detected at Gamma-rays up to TeV energies. Despite its intensive study, the physical processes responsible for the high-energy emission in the source remain unanswered. We present a detailed kinematics study of the source and its connection to Gamma-ray emission. The sub-parsec scale radio structure is dominated by slow-moving features in both the eastern and western lanes of the jet. The jet appears to have accelerated to its maximum speed within less than 125 000 gravitational radii. The fastest reliably detected speed in the jet was ~0.9 c. This leads to a minimum Lorentz factor of ~1.35. Our analysis suggests the presence of multiple high-energy sites in the source. If Gamma-rays are associated with kinematic changes in the jet, they are being produced in both eastern and western lanes in the jet. Three Gamma-ray flares are contemporaneous with epochs where the slowly moving emission region splits into two sub-regions. We estimate the significance of these events being associated as ~2-3 sigma. We tested our results against theoretical predictions for magnetic reconnection-induced mini-jets and turbulence and find them compatible. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 15 pages, 21 figures, 13 tables, Accepted to ApJ

Ultrafast amplification and non-linear magneto-elastic coupling of coherent magnon modes in an antiferromagnet

Authors: D. Bossini, M. Pancaldi, L. Soumah, M. Basini, F. Mertens, M. Cinchetti, T. Satoh, O. Gomonay, S. Bonetti

Abstract: We investigate the role of domain walls in the ultrafast magnon dynamics of an antiferromagnetic NiO single crystal in a pump-probe experiment with variable pump photon energy. Analysing the amplitude of the energy-dependent photo-induced ultrafast spin dynamics, we detect a yet unreported coupling between the material's characteristic THz- and a GHz-magnon modes. We explain this unexpected coupli… ▽ More We investigate the role of domain walls in the ultrafast magnon dynamics of an antiferromagnetic NiO single crystal in a pump-probe experiment with variable pump photon energy. Analysing the amplitude of the energy-dependent photo-induced ultrafast spin dynamics, we detect a yet unreported coupling between the material's characteristic THz- and a GHz-magnon modes. We explain this unexpected coupling between two orthogonal eigenstates of the corresponding Hamiltonian by modelling the magneto-elastic interaction between spins in different domains. We find that such interaction, in the non-linear regime, couples the two different magnon modes via the domain walls and it can be optically exploited via the exciton-magnon resonance. △ Less

Submitted 12 August, 2021; v1 submitted 20 March, 2021; originally announced March 2021.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 127, 077202 (2021)

Temperature dependence of the interface spin Seebeck effect

Authors: F. N. Kholid, D. Hamara, M. Terschanski, F. Mertens, D. Bossini, M. Cinchetti, L. McKenzie-Sell, J. Patchett, D. Petit, R. Cowburn, J. Robinson, J. Barker, C. Ciccarelli

Abstract: We performed temperature-dependent optical pump - THz emission measurements in Y3Fe5O12 (YIG)|Pt from 5 K to room temperature in the presence of an externally applied magnetic field. We study the temperature dependence of the spin Seebeck effect and observe a continuous increase as temperature is decreased, opposite to what is observed in electrical measurements where the spin Seebeck effect is su… ▽ More We performed temperature-dependent optical pump - THz emission measurements in Y3Fe5O12 (YIG)|Pt from 5 K to room temperature in the presence of an externally applied magnetic field. We study the temperature dependence of the spin Seebeck effect and observe a continuous increase as temperature is decreased, opposite to what is observed in electrical measurements where the spin Seebeck effect is suppressed as 0K is approached. By quantitatively analysing the different contributions we isolate the temperature dependence of the spin-mixing conductance and observe features that are correlated to the bands of magnon spectrum in YIG. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 11 pages, 9 figures

21cm signal predictions at Cosmic Dawn and Reionization with coupled radiative-hydrodynamics

Authors: Nicolas J. F. Gillet, Dominique Aubert, Florent G. Mertens, Pierre Ocvirk

Abstract: The process of heating and reionization of the Universe at high redshift links small scale structures/galaxy formation and large scale inter-galactic medium properties. Even if the first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21cm brightness temperature g… ▽ More The process of heating and reionization of the Universe at high redshift links small scale structures/galaxy formation and large scale inter-galactic medium properties. Even if the first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21cm brightness temperature global signal and fluctuations. The need of large scale simulations is therefore strong to understand the properties of the IGM that will be observed. But at the same time the urge to resolve the structures responsible of those process is important. We introduce in this study, a coupled hydro-radiative transfer simulations of the Cosmic Dawn and Reionization with a simple sub-grid star formation process developed and calibrated on the state of the art simulation CoDaII. This scheme permits to follow consistently dark matter, hydrodynamics and radiative transfer evolution's on large scales, while the sub-grid models bridges to the galaxy formation scale. We process the simulation to produce 21cm signal as close as possible to the observations. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Comparing Foreground Removal Techniques for Recovery of the LOFAR-EoR 21cm Power Spectrum

Authors: Ian Hothi, Emma Chapman, Jonathan R. Pritchard, F. G. Mertens, L. V. E Koopmans, B. Ciardi, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, S. Zaroubi

Abstract: We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limi… ▽ More We compare various foreground removal techniques that are being utilised to remove bright foregrounds in various experiments aiming to detect the redshifted 21cm signal of neutral hydrogen from the Epoch of Reionization. In this work, we test the performance of removal techniques (FastICA, GMCA, and GPR) on 10 nights of LOFAR data and investigate the possibility of recovering the latest upper limit on the 21cm signal. Interestingly, we find that GMCA and FastICA reproduce the most recent 2$σしぐま$ upper limit of $Δでるた^2_{21} <$ (73)$^2$ mK$^2$ at $k=0.075~ h \mathrm{cMpc}^{-1}$, which resulted from the application of GPR. We also find that FastICA and GMCA begin to deviate from the noise-limit at \textit{k}-scales larger than $\sim 0.1 ~h \mathrm{cMpc}^{-1}$. We then replicate the data via simulations to see the source of FastICA and GMCA's limitations, by testing them against various instrumental effects. We find that no single instrumental effect, such as primary beam effects or mode-mixing, can explain the poorer recovery by FastICA and GMCA at larger \textit{k}-scales. We then test scale-independence of FastICA and GMCA, and find that lower \textit{k}-scales can be modelled by a smaller number of independent components. For larger scales ($k \gtrsim 0.1~h \mathrm{cMpc}^{-1}$), more independent components are needed to fit the foregrounds. We conclude that, the current usage of GPR by the LOFAR collaboration is the appropriate removal technique. It is both robust and less prone to overfitting, with future improvements to GPR's fitting optimisation to yield deeper limits. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 15 pages, 11 figures, accepted by MNRAS on 02 November 2020

The AARTFAAC Cosmic Explorer: observations of the 21-cm power spectrum in the EDGES absorption trough

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, A. Shulevski, M. Mevius, M. A. Brentjens, M. Kuiack, V. N. Pandey, A. Rowlinson, A. M. Sardarabadi, H. K. Vedantham, R. A. M. J. Wijers, S. Yatawatta, S. Zaroubi

Abstract: The 21-cm absorption feature reported by the EDGES collaboration is several times stronger than that predicted by traditional astrophysical models. If genuine, a deeper absorption may lead to stronger fluctuations on the 21-cm signal on degree scales (up to 1~Kelvin in rms), allowing these fluctuations to be detectable in nearly 50~times shorter integration times compared to previous predictions.… ▽ More The 21-cm absorption feature reported by the EDGES collaboration is several times stronger than that predicted by traditional astrophysical models. If genuine, a deeper absorption may lead to stronger fluctuations on the 21-cm signal on degree scales (up to 1~Kelvin in rms), allowing these fluctuations to be detectable in nearly 50~times shorter integration times compared to previous predictions. We commenced the "AARTFAAC Cosmic Explorer" (ACE) program, that employs the AARTFAAC wide-field imager, to measure or set limits on the power spectrum of the 21-cm fluctuations in the redshift range $z = 17.9-18.6$ ($Δでるたνにゅー= 72.36-75.09$~MHz) corresponding to the deep part of the EDGES absorption feature. Here, we present first results from two LST bins: 23.5-23.75h and 23.5-23.75h, each with 2~h of data, recorded in `semi drift-scan' mode. We demonstrate the application of the new ACE data-processing pipeline (adapted from the LOFAR-EoR pipeline) on the AARTFAAC data. We observe that noise estimates from the channel and time-differenced Stokes~$V$ visibilities agree with each other. After 2~h of integration and subtraction of bright foregrounds, we obtain $2σしぐま$ upper limits on the 21-cm power spectrum of $Δでるた_{21}^2 < (8139~\textrm{mK})^2$ and $Δでるた_{21}^2 < (8549~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$ for the two LST bins. Incoherently averaging the noise bias-corrected power spectra for the two LST bins yields an upper limit of $Δでるた_{21}^2 < (7388~\textrm{mK})^2$ at $k = 0.144~h\,\textrm{cMpc}^{-1}$. These are the deepest upper limits thus far at these redshifts. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Interpreting LOFAR 21-cm signal upper limits at z~9.1 in the context of high-z galaxy and reionisation observations

Authors: Bradley Greig, Andrei Mesinger, Léon V. E. Koopmans, Benedetta Ciardi, Garrelt Mellema, Saleem Zaroubi, Sambit K. Giri, Raghunath Ghara, Abhik Ghosh, Ilian T. Iliev, Florent G. Mertens, Rajesh Mondal, André R. Offringa, Vishambhar N. Pandey

Abstract: Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, whic… ▽ More Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2σしぐま$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 12 pages, 4 pages, 1 table. Submitted to MNRAS, comments welcome

Foreground modelling via Gaussian process regression: an application to HERA data

Authors: Abhik Ghosh, Florent Mertens, Gianni Bernardi, Mário G. Santos, Nicholas S. Kern, Christopher L. Carilli, Trienko L. Grobler, Léon V. E. Koopmans, Daniel C. Jacobs, Adrian Liu, Aaron R. Parsons, Miguel F. Morales, James E. Aguirre, Joshua S. Dillon, Bryna J. Hazelton, Oleg M. Smirnov, Bharat K. Gehlot, Siyanda Matika, Paul Alexander, Zaki S. Ali, Adam P. Beardsley, Roshan K. Benefo, Tashalee S. Billings, Judd D. Bowman, Richard F. Bradley , et al. (48 additional authors not shown)

Abstract: The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize… ▽ More The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz. △ Less

Submitted 12 May, 2020; v1 submitted 13 April, 2020; originally announced April 2020.

Comments: 15 pages, 15 figures, 1 table, Accepted to MNRAS

Tight Constraints on the Excess Radio Background at $z = 9.1$ from LOFAR

Authors: R. Mondal, A. Fialkov, C. Fling, I. T. Iliev, R. Barkana, B. Ciardi, G. Mellema, S. Zaroubi, L. V. E Koopmans, F. G. Mertens, B. K. Gehlot, R. Ghara, A. Ghosh, S. K. Giri, A. Offringa, V. N. Pandey

Abstract: The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at 1.42 GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit… ▽ More The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at 1.42 GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen. We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on 141 hours of observations with LOFAR to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 95% C.L.) that the cosmological radio background can be at most 9.6% of the CMB at 1.42 GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-9.6$% over the CMB (at 1.42 GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$ and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data constrain scenarios with inefficient X-ray sources and cases where the Universe was ionized by stars in massive halos only. △ Less

Submitted 13 August, 2020; v1 submitted 1 April, 2020; originally announced April 2020.

Comments: Published in Monthly Notices of the Royal Astronomical Society (MNRAS). Available at https://doi.org/10.1093/mnras/staa2422

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 498, Issue 3, November 2020, Pages 4178-4191

Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\boldsymbol{z \approx 9.1}$ from LOFAR

Authors: F. G. Mertens, M. Mevius, L. V. E Koopmans, A. R. Offringa, G. Mellema, S. Zaroubi, M. A. Brentjens, H. Gan, B. K. Gehlot, V. N. Pandey, A. M. Sardarabadi, H. K. Vedantham, S. Yatawatta, K. M. B. Asad, B. Ciardi, E. Chapman, S. Gazagnes, R. Ghara, A. Ghosh, S. K. Giri, I. T. Iliev, V. Jelić, R. Kooistra, R. Mondal, J. Schaye , et al. (1 additional authors not shown)

Abstract: A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due… ▽ More A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$σしぐま$ upper limit of $Δでるた^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power. △ Less

Submitted 19 February, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 27 pages, 14 figues, accepted in MNRAS (updated with reference to accompanying paper)

Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations

Authors: R. Ghara, S. K. Giri, G. Mellema, B. Ciardi, S. Zaroubi, I. T. Iliev, L. V. E. Koopmans, E. Chapman, S. Gazagnes, B. K. Gehlot, A. Ghosh, V. Jelic, F. G. Mertens, R. Mondal, J. Schaye, M. B. Silva, K. M. B. Asad, R. Kooistra, M. Mevius, A. R. Offringa, V. N. Pandey, S. Yatawatta

Abstract: We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY an… ▽ More We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 22 pages, 15 Figures, 5 tables, Accepted for publication in MNRAS

Reliable computational prediction of supramolecular ordering of complex molecules under electrochemical conditions

Authors: Benedikt Hartl, Shubham Sharma, Oliver Brügner, Stijn F. L. Mertens, Michael Walter, Gerhard Kahl

Abstract: We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on a metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide reference data for the energies (evaluated for archetypical configurations) to fit the parameters of a conceptually much simpler and computationally less expensive mo… ▽ More We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on a metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide reference data for the energies (evaluated for archetypical configurations) to fit the parameters of a conceptually much simpler and computationally less expensive model of the molecules: classical, spherical particles, representing the respective atomic entities, a soft but perfectly conductive wall potential represents the metallic surface. Stage two feeds the energies that emerge from this classical model into highly efficient and reliable optimization techniques to identify via energy minimization the ordered ground state configurations of the molecules. We demonstrate the power of our approach by successfully reproducing, on a semi-quantitative level, the intricate supramolecular ordering observed experimentally for PQP$^+$ and ClO$_4^-$ molecules at an Au(111)-electrolyte interface, including the formation of open-porous, self-hosts--guest, and stratified bilayer phases as a function of the electric field at the solid--liquid interface. We also discuss the role of the perchlorate ions in the self-assembly process, whose positions could not be identified in the related experimental investigations. △ Less

Submitted 30 April, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

Authors: Fred Cooper, Avinash Khare, Niurka R. Quintero, Bernardo Sánchez-Rey, Franz G. Mertens, Avadh Saxena

Abstract: The damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $κかっぱ$ is analyzed, when the external force is periodic in space and given by $f(x) =r\cos(K x)$, both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order m… ▽ More The damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $κかっぱ$ is analyzed, when the external force is periodic in space and given by $f(x) =r\cos(K x)$, both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force $λらむだ=2 πぱい/K$, and the soliton width $l_s$, which is a function of the nonlinearity $κかっぱ$ as well as the initial frequency $ωおめが_0$ of the solitary wave, there is a transition (at fixed $ωおめが_0$) from trapped to unbound behavior of the soliton, which depends on the parameters $r$ and $K$ of the external force and the nonlinearity parameter $κかっぱ$. We previously studied this phenomena when $κかっぱ=1$ (2019 J. Phys. A: Math. Theor. {\bf 52} 285201) where we showed that for $λらむだ\gg l_s$ the soliton oscillates in an effective potential, while for $λらむだ\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $κかっぱ$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $r$ and $K$ when $κかっぱ=1/2,1,2$ at fixed value of the frequency $ωおめが_0$. Since the solitary wave gets narrower for fixed $ωおめが_0$ as a function of $κかっぱ$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $κかっぱ$. △ Less

Submitted 29 November, 2019; originally announced December 2019.

Comments: 13 pages, 4 figures

Precision requirements for interferometric gridding in 21-cm power spectrum analysis

Authors: A. R. Offringa, F. Mertens, S. van der Tol, B. Veenboer, B. K. Gehlot, L. V. E. Koopmans, M. Mevius

Abstract: We analyse the accuracy of radio interferometric gridding of visibilities with the aim to quantify the Epoch of Reionization (EoR) 21-cm power spectrum bias caused by gridding, ultimately to determine the suitability of different imaging algorithms and gridding settings for 21-cm power spectrum analysis. We simulate realistic LOFAR data, and construct power spectra with convolutional gridding and… ▽ More We analyse the accuracy of radio interferometric gridding of visibilities with the aim to quantify the Epoch of Reionization (EoR) 21-cm power spectrum bias caused by gridding, ultimately to determine the suitability of different imaging algorithms and gridding settings for 21-cm power spectrum analysis. We simulate realistic LOFAR data, and construct power spectra with convolutional gridding and w-stacking, w-projection, image domain gridding and without w-correction. These are compared against directly Fourier transformed data. The influence of oversampling, kernel size, w-quantization, kernel windowing function and image padding are quantified. The gridding excess power is measured with a foreground subtraction strategy, for which foregrounds have been subtracted using Gaussian progress regression, as well as with a foreground avoidance strategy. Constructing a power spectrum that has a bias significantly lower compared to the expected EoR signals is possible with the tested methods, but requires a kernel oversampling factor > 4000 and, when using w-correction, > 500 w-quantization levels. These values are higher than typical values used for imaging, but are computationally feasible. The kernel size and padding factor parameters are less crucial. Among the tested methods, image domain gridding shows the highest accuracy with the lowest imaging time. LOFAR 21-cm power spectrum results are not affected by gridding. Image domain gridding is overall the most suitable algorithm for 21-cm EoR experiments, including for future SKA EoR analyses. Nevertheless, convolutional gridding with tuned parameters results in sufficient accuracy. This holds also for w-stacking for wide-field imaging. The w-projection algorithm is less suitable because of the kernel oversampling requirements, and a faceting approach is unsuitable due to the resulting spatial discontinuities. △ Less

Submitted 29 August, 2019; originally announced August 2019.

Comments: Accepted for publication in A&A ; abstract abridged

Journal ref: A&A 631, A12 (2019)

Peering into the Dark (Ages) with Low-Frequency Space Interferometers

Authors: Leon Koopmans, Rennan Barkana, Mark Bentum, Gianni Bernardi, Albert-Jan Boonstra, Judd Bowman, Jack Burns, Xuelei Chen, Abhirup Datta, Heino Falcke, Anastasia Fialkov, Bharat Gehlot, Leonid Gurvits, Vibor Jelić, Marc Klein-Wolt, Léon Koopmans, Joseph Lazio, Daan Meerburg, Garrelt Mellema, Florent Mertens, Andrei Mesinger, André Offringa, Jonathan Pritchard, Benoit Semelin, Ravi Subrahmanyan , et al. (6 additional authors not shown)

Abstract: Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at… ▽ More Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged] △ Less

Submitted 12 August, 2019; originally announced August 2019.

Comments: 24 pages, 8 figures, 3 tables; White paper submitted to ESA Voyage 2050

The TeV-emitting radio galaxy 3C 264. VLBI kinematics and SED modeling

Authors: B. Boccardi, G. Migliori, P. Grandi, E. Torresi, F. Mertens, V. Karamanavis, R. Angioni, C. Vignali

Abstract: In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE; > 100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferr… ▽ More In March 2018, Mukherjee (2018) reported the detection by VERITAS of very-high-energy emission (VHE; > 100 GeV) from 3C 264. This is the sixth, and second most distant, radio galaxy ever detected in the TeV regime. In this article we present a radio and X-ray analysis of the jet in 3C 264. We determine the main physical parameters of the parsec-scale flow and explore the implications of the inferred kinematic structure for radiative models of this gamma-ray emitting jet. The radio data set is comprised of VLBI observations at 15 GHz from the MOJAVE program, and cover a time period of ~2 years. Through a segmented wavelet decomposition method we estimate the apparent displacement of individual plasma features; we then perform a pixel-based analysis of the stacked image to determine the jet shape. The X-ray data set includes all available observations from the Chandra, XMM, and Swift satellites, and is used, together with archival data in the other bands, to build the SED. Proper motion is mostly detected along the edges of the flow, which appears strongly limb-brightened. The apparent speeds increase as a function of distance from the core up to a maximum of ~$11.5$ c. This constrains the jet viewing angle to assume relatively small values ($θしーた\lesssim10^{\circ}$). In the acceleration region, extending up to a de-projected distance of ~$4.8\times10^4$ Schwarzschild radii (~$11$ pc), the jet is collimating ($r\propto z^{0.40\pm 0.04}$), as predicted for a magnetically-driven plasma flow. By assuming that the core region is indeed magnetically dominated ($U_B/U_e>1$), the SED and the jet power can be well reproduced in the framework of leptonic models, provided that the high-energy component is associated to a second emitting region. The possibility that this region is located at the end of the acceleration zone, either in the jet layer or in the spine, is explored in the modeling. △ Less

Submitted 16 May, 2019; originally announced May 2019.

Comments: 12 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 627, A89 (2019)

The impact of interference excision on 21-cm Epoch of Reionization power spectrum analyses

Authors: A. R. Offringa, F. Mertens, L. V. E. Koopmans

Abstract: We investigate the implications of interference detection for experiments that are pursuing a detection of the redshifted 21-cm signals from the Epoch of Reionization. Interference detection causes samples to be sporadically flagged and rejected. As a necessity to reduce the data volume, flagged samples are typically (implicitly) interpolated during time or frequency averaging or uv-gridding. This… ▽ More We investigate the implications of interference detection for experiments that are pursuing a detection of the redshifted 21-cm signals from the Epoch of Reionization. Interference detection causes samples to be sporadically flagged and rejected. As a necessity to reduce the data volume, flagged samples are typically (implicitly) interpolated during time or frequency averaging or uv-gridding. This so-far unexplored systematic biases the 21-cm power spectrum, and it is important to understand this bias for current 21-cm experiments as well as the upcoming SKA Epoch of Reionization experiment. We analyse simulated data using power spectrum analysis and Gaussian process regression. We find that the combination of flagging and averaging causes tiny spectral fluctuations, resulting in `flagging excess power'. This excess power does not substantially average down over time and, without extra mitigation techniques, can exceed the power of realistic models of the 21-cm reionization signals in LOFAR observations. We mitigate the bias by i) implementing a novel way to average data using a Gaussian-weighted interpolation scheme; ii) using unitary instead of inverse-variance weighting of visibilities; and iii) using low-resolution forward modelling of the data. After these modifications, which have been integrated in the LOFAR EoR processing pipeline, the excess power reduces by approximately three orders of magnitude, and is no longer preventing a detection of the 21-cm signals. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: 11 pages, 8 figures, accepted for publication in MNRAS

The first power spectrum limit on the 21-cm signal of neutral hydrogen during the Cosmic Dawn at z=20-25 from LOFAR

Authors: B. K. Gehlot, F. G. Mertens, L. V. E. Koopmans, M. A. Brentjens, S. Zaroubi, B. Ciardi, A. Ghosh, M. Hatef, I. T. Iliev, V. Jelić, R. Kooistra, F. Krause, G. Mellema, M. Mevius, M. Mitra, A. R. Offringa, V. N. Pandey, A. M. Sardarabadi, J. Schaye, M. B. Silva, H. K. Vedantham, S. Yatawatta

Abstract: Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz fr… ▽ More Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2σしぐま$ upper limit on the 21-cm power spectrum of $Δでるた_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $Δでるた_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes. △ Less

Submitted 20 July, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 19 pages, 16 figures, accepted for publication in MNRAS

Global millimeter VLBI array survey of ultracompact extragalactic radio sources at 86 GHz

Authors: Dhanya G. Nair, Andrei P. Lobanov, Thomas P. Krichbaum, Eduardo Ros, J. Anton Zensus, Yuri Y. Kovalev, Sang-Sung Lee, Florent Mertens, Yoshiaki Hagiwara, Michael Bremer, Michael Lindqvist, Pablo de Vicente

Abstract: (abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $λらむだ= 3\,$mm) reach a resolution of about 50 $μみゅー$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the G… ▽ More (abridged) Very long baseline interferometry (VLBI) observations at 86$\,$GHz (wavelength, $λらむだ= 3\,$mm) reach a resolution of about 50 $μみゅー$as, probing the collimation and acceleration regions of relativistic outflows in active galactic nuclei. To extend the statistical studies of compact extragalactic jets, a large global 86 GHz VLBI survey of 162 radio sources was conducted in 2010-2011 using the Global Millimeter VLBI Array. The survey data attained a typical baseline sensitivity of 0.1 Jy and a typical image sensitivity of 5 mJy/beam, providing successful detections and images for all of the survey targets. For 138 objects, the survey provides the first ever VLBI images made at 86 GHz. Gaussian model fitting of the visibility data was applied to represent the structure of the sources. The Gaussian model-fit-based estimates of brightness temperature ($T_\mathrm{b}$) at the jet base (core) and in moving regions (jet components) downstream from the core were compared to the estimates of $T_\mathrm{b}$ limits made directly from the visibility data, demonstrating a good agreement between the two methods. The apparent brightness temperature estimates for the jet cores in our sample range from $2.5 \times 10^{9}\,$K to $ 1.3\times 10^{12}\,$K. A population model with a single intrinsic value of brightness temperature, $T_\mathrm{0}$, is applied to reproduce the observed $T_\mathrm{b}$ distribution. It yields $T_\mathrm{0} = (3.77^{+0.10}_{-0.14}) \times 10^{11}\,$K for the jet cores, implying that the inverse Compton losses dominate the emission. In the jet components, $T_\mathrm{0} =(1.42^{+0.16}_{-0.19})\times 10^{11}\,$K is found, slightly higher than the equipartition limit of $\sim5\times 10^{10}\,$K expected for these jet regions. For objects with sufficient structural detail detected, the adiabatic energy losses dominate the observed changes of $T_\mathrm{b}$ along the jet. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 25 pages, 12 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 622, A92 (2019)

Nuclear spin dynamics influenced and detected by electron spin polarization in CdTe/CdMgTe quantum wells

Authors: E. Evers, T. Kazimierczuk, F. Mertens, D. R. Yakovlev, G. Karczewski, T. Wojtowicz, J. Kossut, M. Bayer, A. Greilich

Abstract: Nuclear spin coherence and relaxation dynamics of all constituent isotopes of an n-doped CdTe/(Cd,Mg)Te quantum well structure are studied employing optically detected nuclear magnetic resonance. Using time-resolved pump-probe Faraday ellipticity, we generate and detect the coherent spin dynamics of the resident electrons. The photogenerated electron spin polarization is transferred into the nucle… ▽ More Nuclear spin coherence and relaxation dynamics of all constituent isotopes of an n-doped CdTe/(Cd,Mg)Te quantum well structure are studied employing optically detected nuclear magnetic resonance. Using time-resolved pump-probe Faraday ellipticity, we generate and detect the coherent spin dynamics of the resident electrons. The photogenerated electron spin polarization is transferred into the nuclear spin system, which becomes polarized and acts back on the electron spins as the Overhauser field. Under the influence of resonant radio frequency pulses, we trace the coherent spin dynamics of the nuclear isotopes $^{111}$Cd, $^{113}$Cd, and $^{125}$Te. We measure nuclear Rabi oscillations, the inhomogeneous dephasing time $T_2^*$, the spin coherence time $T_2$, and the longitudinal relaxation time $T_1$. Furthermore, we investigate the influence of the laser excitation and the corresponding electron spin polarization on the nuclear spin relaxation time and find a weak extension of this time induced by interaction with the electron spins. △ Less

Submitted 3 April, 2019; v1 submitted 25 July, 2018; originally announced July 2018.

Comments: 5 pages, 2 figures

Journal ref: Phys. Rev. B 99, 045303 (2019)

Copper underpotential deposition on boron nitride nanomesh

Authors: Stijn F. L. Mertens

Abstract: The boron nitride nanomesh is a corrugated monolayer of hexagonal boron nitride (h-BN) on Rh(111), which so far has been studied mostly under ultrahigh vacuum conditions. Here, we investigate how copper underpotential deposition (upd) can be used to quantify defects in the boron nitride monolayer and to assess the potential window of the nanomesh, which is important to explore its functionality un… ▽ More The boron nitride nanomesh is a corrugated monolayer of hexagonal boron nitride (h-BN) on Rh(111), which so far has been studied mostly under ultrahigh vacuum conditions. Here, we investigate how copper underpotential deposition (upd) can be used to quantify defects in the boron nitride monolayer and to assess the potential window of the nanomesh, which is important to explore its functionality under ambient and electrochemical conditions. In dilute sulfuric acid, the potential window of h-BN/Rh(111) is close to 1 volt, i.e. larger than that of the Rh substrate, and is limited by molecular hydrogen evolution on the negative and by oxidative removal on the positive side. From copper upd on pristine h-BN/Rh(111) wafer samples, we estimate a collective defect fraction on the order of 0.08-0.7% of the geometric area, which may arise from line and point defects in the h-BN layer that are created during its chemical vapour deposition. Overpotential deposition (opd) is demonstrated to have significant consequences on the defect area. We hypothesise that this non-innocent Cu electrodeposition involves intercalation originating at initial defects, causing irreversible delamination of the h-BN layer; this effect may be used for 2D material nanoengineering. On the relevant timescale, upd itself does not alter the defect area on repeated cycling; therefore, metal upd may find use as a general tool to determine the collective defect area in hybrids between 2D materials and various substrate metals. △ Less

Submitted 24 June, 2018; originally announced June 2018.

Comments: 19 pages, 6 figures

Journal ref: Electrochimica Acta 246, 2017, 730-736

Speed-of-light pulses in a nonlinear Weyl equation

Authors: J. Cuevas-Maraver, P. G. Kevrekidis, F. G. Mertens, A. Saxena

Abstract: We introduce a prototypical nonlinear Weyl equation, motivated by recent developments in massless Dirac fermions, topological semimetals and photonics. We study the dynamics of its pulse solutions and find that a localized one-hump initial condition splits into a localized two-hump pulse, while an associated phase structure emerges in suitable components of the spinor field. For times larger than… ▽ More We introduce a prototypical nonlinear Weyl equation, motivated by recent developments in massless Dirac fermions, topological semimetals and photonics. We study the dynamics of its pulse solutions and find that a localized one-hump initial condition splits into a localized two-hump pulse, while an associated phase structure emerges in suitable components of the spinor field. For times larger than a transient time $t_s$ this pulse moves with the speed of light (or Fermi velocity in Weyl semimetals), effectively featuring linear wave dynamics and maintaining its shape (both in two and three dimensions). We show that for the considered nonlinearity, this pulse represents an exact solution of the nonlinear Weyl (NLW) equation. Finally, we comment on the generalization of the results to a broader class of nonlinearities and on their emerging potential for observation in different areas of application. △ Less

Submitted 6 December, 2018; v1 submitted 8 May, 2018; originally announced May 2018.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. E 100, 022210 (2019)

Localization of magnon modes in a curved magnetic nanowire

Authors: Yuri Gaididei, Volodymyr P. Kravchuk, Franz G. Mertens, Oleksandr V. Pylypovskyi, Avadh Saxena, Denis D. Sheka, Oleksii M. Volkov

Abstract: Spin waves in magnetic nanowires can be bound by a local bending of the wire. The eigenfrequency of a truly local magnon mode is determined by the curvature: a general analytical expression is established for any infinitesimally weak localized curvature of the wire. The interaction of the local mode with spin waves, propagating through the bend, results in scattering features, which is well confir… ▽ More Spin waves in magnetic nanowires can be bound by a local bending of the wire. The eigenfrequency of a truly local magnon mode is determined by the curvature: a general analytical expression is established for any infinitesimally weak localized curvature of the wire. The interaction of the local mode with spin waves, propagating through the bend, results in scattering features, which is well confirmed by spin-lattice simulations. △ Less

Submitted 4 July, 2018; v1 submitted 21 January, 2018; originally announced January 2018.

Comments: ReVTeX4-1: 11 pages, 4 figures

Journal ref: Low Temperature Physics, 2018, v. 44, No. 7, pp. 634-643

Robust Foregrounds Removal for 21-cm Experiments

Authors: F. Mertens, A. Ghosh, L. V. E. Koopmans

Abstract: Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emissi… ▽ More Direct detection of the Epoch of Reionization via the redshifted 21-cm line will have unprecedented implications on the study of structure formation in the early Universe. To fulfill this promise current and future 21-cm experiments will need to detect the weak 21-cm signal over foregrounds several order of magnitude greater. This requires accurate modeling of the galactic and extragalactic emission and of its contaminants due to instrument chromaticity, ionosphere and imperfect calibration. To solve for this complex modeling, we propose a new method based on Gaussian Process Regression (GPR) which is able to cleanly separate the cosmological signal from most of the foregrounds contaminants. We also propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere. Using this method, chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: Subbmited to the Proceedings of the IAUS333, Peering Towards Cosmic Dawn, 4 pages, 2 figures

Statistical 21-cm Signal Separation via Gaussian Process Regression Analysis

Authors: F. G. Mertens, A. Ghosh, L. V. E. Koopmans

Abstract: Detecting and characterizing the Epoch of Reionization and Cosmic Dawn via the redshifted 21-cm hyperfine line of neutral hydrogen will revolutionize the study of the formation of the first stars, galaxies, black holes and intergalactic gas in the infant Universe. The wealth of information encoded in this signal is, however, buried under foregrounds that are many orders of magnitude brighter. Thes… ▽ More Detecting and characterizing the Epoch of Reionization and Cosmic Dawn via the redshifted 21-cm hyperfine line of neutral hydrogen will revolutionize the study of the formation of the first stars, galaxies, black holes and intergalactic gas in the infant Universe. The wealth of information encoded in this signal is, however, buried under foregrounds that are many orders of magnitude brighter. These must be removed accurately and precisely in order to reveal the feeble 21-cm signal. This requires not only the modeling of the Galactic and extra-galactic emission, but also of the often stochastic residuals due to imperfect calibration of the data caused by ionospheric and instrumental distortions. To stochastically model these effects, we introduce a new method based on `Gaussian Process Regression' (GPR) which is able to statistically separate the 21-cm signal from most of the foregrounds and other contaminants. Using simulated LOFAR-EoR data that include strong instrumental mode-mixing, we show that this method is capable of recovering the 21-cm signal power spectrum across the entire range $k = 0.07 - 0.3 \ \rm{h\, cMpc^{-1}}$. The GPR method is most optimal, having minimal and controllable impact on the 21-cm signal, when the foregrounds are correlated on frequency scales $\gtrsim 3$\,MHz and the rms of the signal has $σしぐま_{\mathrm{21cm}} \gtrsim 0.1\,σしぐま_{\mathrm{noise}}$. This signal separation improves the 21-cm power-spectrum sensitivity by a factor $\gtrsim 3$ compared to foreground avoidance strategies and enables the sensitivity of current and future 21-cm instruments such as the {\sl Square Kilometre Array} to be fully exploited. △ Less

Submitted 7 May, 2018; v1 submitted 29 November, 2017; originally announced November 2017.

Comments: 15 pages, accepted to MNRAS, augmented with new subsection and appendix, no changes to conclusion

Deconvolving the Wedge: Maximum-Likelihood Power Spectra via Spherical-Wave Visibility Modeling

Authors: Abhik Ghosh, Florent Mertens, Leon V. E. Koopmans

Abstract: Direct detection of the Epoch of Reionization (EoR) via the red-shifted 21-cm line will have unprecedented implications on the study of structure formation in the infant Universe. To fulfill this promise, current and future 21-cm experiments need to detect this weak EoR signal in the presence of foregrounds that are several orders of magnitude larger. This requires extreme noise control and improv… ▽ More Direct detection of the Epoch of Reionization (EoR) via the red-shifted 21-cm line will have unprecedented implications on the study of structure formation in the infant Universe. To fulfill this promise, current and future 21-cm experiments need to detect this weak EoR signal in the presence of foregrounds that are several orders of magnitude larger. This requires extreme noise control and improved wide-field high dynamic-range imaging techniques. We propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere, or equivalently in the $uvw$-domain. The method uses the one-to-one relation between spherical waves and spherical harmonics (SpH). It consistently handles signals from the entire sky, and does not require a $w$-term correction. The spherical-harmonics coefficients represent the sky-brightness distribution and the visibilities in the $uvw$-domain, and provide a direct estimate of the spatial power spectrum. Using these spectrally-smooth SpH coefficients, bright foregrounds can be removed from the signal, including their side-lobe noise, which is one of the limiting factors in high dynamics range wide-field imaging. Chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum, compared to a power spectrum computed directly from the images of the foreground visibilities where the wedge is clearly present. We illustrate our method using simulated LOFAR observations, finding an excellent reconstruction of the input EoR signal with minimal bias. △ Less

Submitted 13 November, 2017; v1 submitted 20 September, 2017; originally announced September 2017.

Comments: 13 pages, 8 figures. Replaced to match accepted MNRAS version; few typos corrected & textual clarification added (no changes to results)

Response of exact solutions of the nonlinear Schrodinger equation to small perturbations in a class of complex external potentials having supersymmetry and parity-time symmetry

Authors: Fred Cooper, John F. Dawson, Franz G. Mertens, Edward Arevalo, Niurka R. Quintero, Bogdan Mihaila, Avinash Khare, Avadh Saxena

Abstract: We discuss the effect of small perturbation on nodeless solutions of the nonlinear \Schrodinger\ equation in 1+1 dimensions in an external complex potential derivable from a parity-time symmetric superpotential that was considered earlier [Phys.~Rev.~E 92, 042901 (2015)]. In particular we consider the nonlinear partial differential equation… ▽ More We discuss the effect of small perturbation on nodeless solutions of the nonlinear \Schrodinger\ equation in 1+1 dimensions in an external complex potential derivable from a parity-time symmetric superpotential that was considered earlier [Phys.~Rev.~E 92, 042901 (2015)]. In particular we consider the nonlinear partial differential equation $\{ \, \rmi \, \partial_t + \partial_x^2 + g |ψぷさい(x,t)|^2 - V^{+}(x) \, \} \, ψぷさい(x,t) = 0$, where $V^{+}(x) = \qty( -b^2 - m^2 + 1/4 ) \, \sech^2(x) - 2 i \, m \, b \, \sech(x) \, \tanh(x)$ represents the complex potential. Here we study the perturbations as a function of $b$ and $m$ using a variational approximation based on a dissipation functional formalism. We compare the result of this variational approach with direct numerical simulation of the equations. We find that the variational approximation works quite well at small and moderate values of the parameter $b m$ which controls the strength of the imaginary part of the potential. We also show that the dissipation functional formalism is equivalent to the generalized traveling wave method for this type of dissipation. △ Less

Submitted 5 July, 2017; originally announced July 2017.

Comments: 18 pages, 6 figures

Speed-of-light pulses in the massless nonlinear Dirac equation with a potential

Authors: Niurka R. Quintero, Franz G. Mertens, Fred Cooper, Avadh Saxena, A. R. Bishop

Abstract: We consider the massless nonlinear Dirac (NLD) equation in $1+1$ dimension with scalar-scalar self-interaction $\frac{g^2}{2} (\barΨぷさい Ψぷさい)^2$ in the presence of three external electromagnetic potentials $V(x)$, a potential barrier, a constant potential, and a potential well. By solving numerically the NLD equation, we find that, for all three cases, after a short transit time, the initial pulse break… ▽ More We consider the massless nonlinear Dirac (NLD) equation in $1+1$ dimension with scalar-scalar self-interaction $\frac{g^2}{2} (\barΨぷさい Ψぷさい)^2$ in the presence of three external electromagnetic potentials $V(x)$, a potential barrier, a constant potential, and a potential well. By solving numerically the NLD equation, we find that, for all three cases, after a short transit time, the initial pulse breaks into two pulses which are solutions of the massless linear Dirac equation traveling in opposite directions with the speed of light. During this splitting the charge and the energy are conserved, whereas the momentum is conserved when the solutions possess specific symmetries. For the case of the constant potential, we derive exact analytical solutions of the massless NLD equation that are also solutions of the massless linearized Dirac equation. △ Less

Submitted 31 May, 2017; originally announced May 2017.

Comments: 11 pages, 7 figures

Journal ref: Phys. Rev. E 96, 052219 (2017)

Non-autonomous bright solitons and their stability in Rabi coupled binary Bose-Einstein Condensates

Authors: T. Kanna, R. Babu Mareeswaran, Franz G. Mertens

Abstract: The dynamics of non-autonomous bright matter-wave solitons in Rabi coupled binary Bose-Einstein condensates is explored. By performing a unitary and similarity/lens-type transformation, we reduce the non-autonomous Gross-Pitaevskii equation into the celebrated Manakov model. Then, we construct the exact bright solitons of the non-autonomous Gross-Pitaevskii system in the presence of time dependent… ▽ More The dynamics of non-autonomous bright matter-wave solitons in Rabi coupled binary Bose-Einstein condensates is explored. By performing a unitary and similarity/lens-type transformation, we reduce the non-autonomous Gross-Pitaevskii equation into the celebrated Manakov model. Then, we construct the exact bright solitons of the non-autonomous Gross-Pitaevskii system in the presence of time dependent nonlinearities for two specific forms, namely hyperbolic nonlinearities, which are of physical interest. The experimental possibilities of realizing the forms of temporally modulated potentials corresponding to these time dependent nonlinearities and supporting such localized structures are explored. Our study on the propagation of one soliton shows that the amplitude, velocity and shape of the bright soliton are altered by the time-dependent scattering length. We also analyse the non-trivial energy sharing collision of Manakov solitons in the presence of Rabi coupling and aforementioned nonlinearities. We find that breathers arise in two-soliton collisions and the nature of energy sharing collisions is altered from that of Manakov system due to Rabi coupling only. Further in the presence of time-dependent nonlinearities the collision scenario is again altered significantly. Finally, the stability of these localized structures is examined using a recently developed powerful analytic method by Quintero et. al., [Phys. Rev. E 91, 012905 (2015)] and it is shown that the non- autonomous bright solitons are indeed stable. The evolution of position and velocity is also studied. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Comments: 20 pages, 8 figures