-
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.
-
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.
-
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.
-
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.
-
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.
-
Phase offset method of ptychographic contrast reversal correction
Authors:
Christoph Hofer,
Chuang Gao,
Tamazouzt Chennit,
Biao Yuan,
Timothy J. Pennycook
Abstract:
The contrast transfer function of direct ptychography methods such as the single side band (SSB) method are single signed, yet these methods still sometimes exhibit contrast reversals, most often where the projected potentials are strong. In thicker samples central focusing often provides the best ptychographic contrast as this leads to defocus variations within the sample canceling out. However f…
▽ More
The contrast transfer function of direct ptychography methods such as the single side band (SSB) method are single signed, yet these methods still sometimes exhibit contrast reversals, most often where the projected potentials are strong. In thicker samples central focusing often provides the best ptychographic contrast as this leads to defocus variations within the sample canceling out. However focusing away from the entrance surface is often undesirable as this degrades the annular dark field (ADF) signal. Here we discuss how phase wrap asymptotes in the frequency response of SSB ptychography give rise to contrast reversals, without the need for dynamical scattering, and how these can be counteracted by manipulating the phases such that the asymptotes are either shifted to higher frequencies or damped via amplitude modulation. This is what enables post collection defocus correction of contrast reversals. However, the phase offset method of counteracting contrast reversals we introduce here is generally found to be superior to post collection application of defocus, with greater reliability and generally stronger contrast. Importantly, the phase offset method also works for thin and thick samples where central focusing does not. Finally, the independence of the method from focus is useful for optical sectioning involving ptychography, improving interpretability by better disentangling the effects of strong potentials and focus.
△ Less
Submitted 21 December, 2023; v1 submitted 12 October, 2023;
originally announced October 2023.
-
A Detection of Cosmological 21 cm Emission from CHIME in Cross-correlation with eBOSS Measurements of the Lyman-$α$ Forest
Authors:
CHIME Collaboration,
Mandana Amiri,
Kevin Bandura,
Arnab Chakraborty,
Matt Dobbs,
Mateus Fandino,
Simon Foreman,
Hyoyin Gan,
Mark Halpern,
Alex S. Hill,
Gary Hinshaw,
Carolin Höfer,
T. L. Landecker,
Zack Li,
Joshua MacEachern,
Kiyoshi Masui,
Juan Mena-Parra,
Nikola Milutinovic,
Arash Mirhosseini,
Laura Newburgh,
Anna Ordog,
Sourabh Paul,
Ue-Li Pen,
Tristan Pinsonneault-Marotte,
Alex Reda
, et al. (6 additional authors not shown)
Abstract:
We report the detection of 21 cm emission at an average redshift $\bar{z} = 2.3$ in the cross-correlation of data from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) with measurements of the Lyman-$α$ forest from eBOSS. Data collected by CHIME over 88 days in the $400-500$~MHz frequency band ($1.8 < z < 2.5$) are formed into maps of the sky and high-pass delay filtered to suppress the…
▽ More
We report the detection of 21 cm emission at an average redshift $\bar{z} = 2.3$ in the cross-correlation of data from the Canadian Hydrogen Intensity Mapping Experiment (CHIME) with measurements of the Lyman-$α$ forest from eBOSS. Data collected by CHIME over 88 days in the $400-500$~MHz frequency band ($1.8 < z < 2.5$) are formed into maps of the sky and high-pass delay filtered to suppress the foreground power, corresponding to removing cosmological scales with $k_\parallel \lesssim 0.13\ \text{Mpc}^{-1}$ at the average redshift. Line-of-sight spectra to the eBOSS background quasar locations are extracted from the CHIME maps and combined with the Lyman-$α$ forest flux transmission spectra to estimate the 21 cm-Lyman-$α$ cross-correlation function. Fitting a simulation-derived template function to this measurement results in a $9σ$ detection significance. The coherent accumulation of the signal through cross-correlation is sufficient to enable a detection despite excess variance from foreground residuals $\sim6-10$ times brighter than the expected thermal noise level in the correlation function. These results are the highest-redshift measurement of \tcm emission to date, and set the stage for future 21 cm intensity mapping analyses at $z>1.8$.
△ Less
Submitted 8 September, 2023;
originally announced September 2023.
-
Reliable phase quantification in focused probe electron ptychography of thin materials
Authors:
Christoph Hofer,
Timothy J. Pennycook
Abstract:
Electron ptychography provides highly sensitive, dose efficient phase images which can be corrected for aberrations after the data has been acquired. This is crucial when very precise quantification is required, such as with sensitivity to charge transfer due to bonding. Drift can now be essentially eliminated as a major impediment to focused probe ptychography, which benefits from the availabilit…
▽ More
Electron ptychography provides highly sensitive, dose efficient phase images which can be corrected for aberrations after the data has been acquired. This is crucial when very precise quantification is required, such as with sensitivity to charge transfer due to bonding. Drift can now be essentially eliminated as a major impediment to focused probe ptychography, which benefits from the availability of easily interpretable simultaneous Z-contrast imaging. However challenges have remained when quantifying the ptychographic phases of atomic sites. The phase response of a single atom has a negative halo which can cause atoms to reduce in phase when brought closer together. When unaccounted for, as in integrating methods of quantification, this effect can completely obscure the effects of charge transfer. Here we provide a new method of quantification that overcomes this challenge, at least for 2D materials, and is robust to experimental parameters such as noise, sample tilt.
△ Less
Submitted 26 July, 2023;
originally announced July 2023.
-
On central focusing for contrast optimization in direct electron ptychography of thick samples
Authors:
C. Gao,
C. Hofer,
T. J. Pennycook
Abstract:
Ptychography provides high dose efficiency images that can reveal light elements next to heavy atoms. However, despite ptychography having an otherwise single signed contrast transfer function, contrast reversals can occur when the projected potential becomes strong for both direct and iterative inversion ptychography methods. It has recently been shown that these reversals can often be counteract…
▽ More
Ptychography provides high dose efficiency images that can reveal light elements next to heavy atoms. However, despite ptychography having an otherwise single signed contrast transfer function, contrast reversals can occur when the projected potential becomes strong for both direct and iterative inversion ptychography methods. It has recently been shown that these reversals can often be counteracted in direct ptychography methods by adapting the focus. Here we provide an explanation of why the best contrast is often found with the probe focused to the middle of the sample. The phase contribution due to defocus at each sample slice above and below the central plane in this configuration effectively cancels out, which can prevent contrast reversals when dynamical scattering effects are not overly strong. In addition we show that the convergence angle can be an important consideration for removal of contrast reversals in relatively thin samples.
△ Less
Submitted 14 June, 2023;
originally announced June 2023.
-
Detecting charge transfer at defects in 2D materials with electron ptychography
Authors:
Christoph Hofer,
Jacob Madsen,
Toma Susi,
Timothy J. Pennycook
Abstract:
Charge transfer between atoms is fundamental to chemical bonding but has remained very challenging to detect directly in real space. Atomic-resolution imaging of charge density is not sufficient by itself, as the change in the density due to bonding is very subtle compared to the total local charge density. Sufficiently high sensitivity, precision and accuracy are required, which we demonstrate he…
▽ More
Charge transfer between atoms is fundamental to chemical bonding but has remained very challenging to detect directly in real space. Atomic-resolution imaging of charge density is not sufficient by itself, as the change in the density due to bonding is very subtle compared to the total local charge density. Sufficiently high sensitivity, precision and accuracy are required, which we demonstrate here for the detection of charge transfer at defects in two-dimensional WS\textsubscript{2} via high-speed electron ptychography and its ability to correct errors due to residual lens aberrations.
△ Less
Submitted 19 January, 2024; v1 submitted 11 January, 2023;
originally announced January 2023.
-
Characterization of the John A. Galt telescope for radio holography with CHIME
Authors:
Alex Reda,
Tristan Pinsonneault-Marotte,
Meiling Deng,
Mandana Amiri,
Kevin Bandura,
Arnab Chakraborty,
Simon Foreman,
Mark Halpern,
Alex S. Hill,
Carolin Höfer,
Joseph Kania,
T. L. Landecker,
Joshua MacEachern,
Kiyoshi Masui,
Juan Mena-Parra,
Nikola Milutinovic,
Laura Newburgh,
Anna Ordog,
Sourabh Paul,
J. Richard Shaw,
Seth R. Siegel,
Rick Smegal,
Haochen Wang,
Dallas Wulf
Abstract:
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) will measure the 21 cm emission of astrophysical neutral hydrogen to probe large scale structure at redshifts z=0.8-2.5. However, detecting the 21 cm signal beneath substantially brighter foregrounds remains a key challenge. Due to the high dynamic range between 21 cm and foreground emission, an exquisite calibration of instrument systemat…
▽ More
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) will measure the 21 cm emission of astrophysical neutral hydrogen to probe large scale structure at redshifts z=0.8-2.5. However, detecting the 21 cm signal beneath substantially brighter foregrounds remains a key challenge. Due to the high dynamic range between 21 cm and foreground emission, an exquisite calibration of instrument systematics, notably the telescope beam, is required to successfully filter out the foregrounds. One technique being used to achieve a high fidelity measurement of the CHIME beam is radio holography, wherein signals from each of CHIME's analog inputs are correlated with the signal from a co-located reference antenna, the 26 m John A. Galt telescope, as the 26 m Galt telescope tracks a bright point source transiting over CHIME. In this work we present an analysis of several of the Galt telescope's properties. We employ driftscan measurements of several bright sources, along with background estimates derived from the 408 MHz Haslam map, to estimate the Galt system temperature. To determine the Galt telescope's beam shape, we perform and analyze a raster scan of the bright radio source Cassiopeia A. Finally, we use early holographic measurements to measure the Galt telescope's geometry with respect to CHIME for the holographic analysis of the CHIME and Galt interferometric data set.
△ Less
Submitted 30 September, 2022; v1 submitted 28 July, 2022;
originally announced July 2022.
-
Overcoming contrast reversals in focused probe ptychography of thick materials: an optimal pipeline for efficiently determining local atomic structure in materials science
Authors:
C. Gao,
C. Hofer,
D. Jannis,
A. Béché,
J. Verbeeck,
T. J. Pennycook
Abstract:
Ptychography provides highly efficient imaging in scanning transmission electron microscopy (STEM), but questions have remained over its applicability to strongly scattering samples such as those most commonly seen in materialsscience. Although contrast reversals can appear in ptychographic phase images as the projected potentials of the sample increase, we show here how these can be easily overco…
▽ More
Ptychography provides highly efficient imaging in scanning transmission electron microscopy (STEM), but questions have remained over its applicability to strongly scattering samples such as those most commonly seen in materialsscience. Although contrast reversals can appear in ptychographic phase images as the projected potentials of the sample increase, we show here how these can be easily overcome by a small amount of defocus. The amount of defocus is small enough that it can exist naturally when focusing using the annular dark field (ADF) signal, but can also be adjusted post acquisition. The ptychographic images of strongly scattering materials are clearer at finite doses than other STEM techniques, and can better reveal light atomic columns within heavy lattices. In addition data for ptychography can now be collected simultaneously with the fastest of ADF scans. This combination of sensitivity and interpretability presents an ideal workflow for materials science.
△ Less
Submitted 4 August, 2022; v1 submitted 26 May, 2022;
originally announced May 2022.
-
Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
Authors:
CHIME Collaboration,
Mandana Amiri,
Kevin Bandura,
Tianyue Chen,
Meiling Deng,
Matt Dobbs,
Mateus Fandino,
Simon Foreman,
Mark Halpern,
Alex S. Hill,
Gary Hinshaw,
Carolin Höfer,
Joseph Kania,
T. L. Landecker,
Joshua MacEachern,
Kiyoshi Masui,
Juan Mena-Parra,
Nikola Milutinovic,
Arash Mirhosseini,
Laura Newburgh,
Anna Ordog,
Ue-Li Pen,
Tristan Pinsonneault-Marotte,
Ava Polzin,
Alex Reda
, et al. (8 additional authors not shown)
Abstract:
We present a detection of 21-cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). Radio observations acquired over 102 nights are used to construct maps which are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRG), emission line galaxies (ELG), and quasars…
▽ More
We present a detection of 21-cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). Radio observations acquired over 102 nights are used to construct maps which are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRG), emission line galaxies (ELG), and quasars (QSO) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes Factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood-ratio test, yields a detection significance of $7.1σ$ (LRG), $5.7σ$ (ELG), and $11.1σ$ (QSO). These are the first 21-cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (HI), defined as $\mathcal{A}_{\rm HI}\equiv 10^{3}\,Ω_\mathrm{HI}\left(b_\mathrm{HI}+\langle\,fμ^{2}\rangle\right)$, where $Ω_\mathrm{HI}$ is the cosmic abundance of HI, $b_\mathrm{HI}$ is the linear bias of HI, and $\langle\,fμ^{2}\rangle=0.552$ encodes the effect of redshift-space distortions at linear order. We find $\mathcal{A}_\mathrm{HI}=1.51^{+3.60}_{-0.97}$ for LRGs $(z=0.84)$, $\mathcal{A}_\mathrm{HI}=6.76^{+9.04}_{-3.79}$ for ELGs $(z=0.96)$, and $\mathcal{A}_\mathrm{HI}=1.68^{+1.10}_{-0.67}$ for QSOs $(z=1.20)$, with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and find a non-zero bias $Δ\,v= -66 \pm 20 \mathrm{km/s}$ for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at $z=1.30$ producing the highest redshift 21-cm intensity mapping measurement thus far.
△ Less
Submitted 2 February, 2022;
originally announced February 2022.
-
Using the Sun to Measure the Primary Beam Response of the Canadian Hydrogen Intensity Mapping Experiment
Authors:
CHIME Collaboration,
Mandana Amiri,
Kevin Bandura,
Anja Boskovic,
Jean-François Cliche,
Meiling Deng,
Matt Dobbs,
Mateus Fandino,
Simon Foreman,
Mark Halpern,
Alex S. Hill,
Gary Hinshaw,
Carolin Höfer,
Joseph Kania,
T. L. Landecker,
Joshua MacEachern,
Kiyoshi Masui,
Juan Mena-Parra,
Laura Newburgh,
Anna Ordog,
Tristan Pinsonneault-Marotte,
Ava Polzin,
Alex Reda,
J. Richard Shaw,
Seth R. Siegel
, et al. (5 additional authors not shown)
Abstract:
We present a beam pattern measurement of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) made using the Sun as a calibration source. As CHIME is a pure drift scan instrument, we rely on the seasonal North-South motion of the Sun to probe the beam at different elevations. This semiannual range in elevation, combined with the radio brightness of the Sun, enables a beam measurement which s…
▽ More
We present a beam pattern measurement of the Canadian Hydrogen Intensity Mapping Experiment (CHIME) made using the Sun as a calibration source. As CHIME is a pure drift scan instrument, we rely on the seasonal North-South motion of the Sun to probe the beam at different elevations. This semiannual range in elevation, combined with the radio brightness of the Sun, enables a beam measurement which spans ~7,200 square degrees on the sky without the need to move the telescope. We take advantage of observations made near solar minimum to minimize the impact of solar variability, which is observed to be <10% in intensity over the observation period. The resulting data set is highly complementary to other CHIME beam measurements -- both in terms of angular coverage and systematics -- and plays an important role in the ongoing program to characterize the CHIME primary beam.
△ Less
Submitted 3 May, 2022; v1 submitted 27 January, 2022;
originally announced January 2022.
-
An Overview of CHIME, the Canadian Hydrogen Intensity Mapping Experiment
Authors:
The CHIME Collaboration,
Mandana Amiri,
Kevin Bandura,
Anja Boskovic,
Tianyue Chen,
Jean-François Cliche,
Meiling Deng,
Nolan Denman,
Matt Dobbs,
Mateus Fandino,
Simon Foreman,
Mark Halpern,
David Hanna,
Alex S. Hill,
Gary Hinshaw,
Carolin Höfer,
Joseph Kania,
Peter Klages,
T. L. Landecker,
Joshua MacEachern,
Kiyoshi Masui,
Juan Mena-Parra,
Nikola Milutinovic,
Arash Mirhosseini,
Laura Newburgh
, et al. (18 additional authors not shown)
Abstract:
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a drift scan radio telescope operating across the 400-800 MHz band. CHIME is located at the Dominion Radio Astrophysical Observatory near Penticton, BC Canada. The instrument is designed to map neutral hydrogen over the redshift range 0.8 to 2.5 to constrain the expansion history of the Universe. This goal drives the design features of…
▽ More
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a drift scan radio telescope operating across the 400-800 MHz band. CHIME is located at the Dominion Radio Astrophysical Observatory near Penticton, BC Canada. The instrument is designed to map neutral hydrogen over the redshift range 0.8 to 2.5 to constrain the expansion history of the Universe. This goal drives the design features of the instrument. CHIME consists of four parallel cylindrical reflectors, oriented north-south, each 100 m $\times$ 20 m and outfitted with a 256 element dual-polarization linear feed array. CHIME observes a two degree wide stripe covering the entire meridian at any given moment, observing 3/4 of the sky every day due to Earth rotation. An FX correlator utilizes FPGAs and GPUs to digitize and correlate the signals, with different correlation products generated for cosmological, fast radio burst, pulsar, VLBI, and 21 cm absorber backends. For the cosmology backend, the $N_\mathrm{feed}^2$ correlation matrix is formed for 1024 frequency channels across the band every 31 ms. A data receiver system applies calibration and flagging and, for our primary cosmological data product, stacks redundant baselines and integrates for 10 s. We present an overview of the instrument, its performance metrics based on the first three years of science data, and we describe the current progress in characterizing CHIME's primary beam response. We also present maps of the sky derived from CHIME data; we are using versions of these maps for a cosmological stacking analysis as well as for investigation of Galactic foregrounds.
△ Less
Submitted 23 May, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
Aligned Stacking of Nanopatterned 2D materials -- Towards 3D printing at atomic resolution
Authors:
Jonas Haas,
Finn Ulrich,
Christoph Hofer,
Xiao Wang,
Kai Braun,
Jannik C. Meyer
Abstract:
Two-dimensional materials can be combined by placing individual layers on top of each other, so that they are bound only by their van der Waals interaction. The sequence of layers can be chosen arbitrarily, enabling an essentially atomic-level control of the material and thereby a wide choice of properties along one dimension. However, simultaneous control over the structure in the in-plane direct…
▽ More
Two-dimensional materials can be combined by placing individual layers on top of each other, so that they are bound only by their van der Waals interaction. The sequence of layers can be chosen arbitrarily, enabling an essentially atomic-level control of the material and thereby a wide choice of properties along one dimension. However, simultaneous control over the structure in the in-plane directions is so far still rather limited. Here, we combine spatially controlled modifications of 2D materials, using focused electron irradiation or electron beam induced etching, with the layer-by-layer assembly of van der Waals heterostructures. A novel assembly process makes it possible to structure each layer with an arbitrary pattern prior to the assembly into the heterostructure. Moreover, it enables a stacking of the layers with accurate lateral alignment, with an accuracy of currently 10nm, under observation in an electron microscope. Together, this enables the fabrication of almost arbitrary 3D structures with highest spatial resolution.
△ Less
Submitted 18 October, 2021; v1 submitted 15 October, 2021;
originally announced October 2021.
-
Towards Exotic Layered Materials: 2D Cuprous Iodide
Authors:
Kimmo Mustonen,
Christoph Hofer,
Peter Kotrusz,
Alexander Markevich,
Martin Hulman,
Clemens Mangler,
Toma Susi,
Timothy J. Pennycook,
Karol Hricovini,
Christine M. Richter,
Jannik C. Meyer,
Jani Kotakoski,
Viera Skakalova
Abstract:
Heterostructures composed of two-dimensional (2D) materials are already opening many new possibilities in such fields of technology as electronics and magnonics, but far more could be achieved if the number and diversity of 2D materials is increased. So far, only a few dozen 2D crystals have been extracted from materials that exhibit a layered phase in ambient conditions, omitting entirely the lar…
▽ More
Heterostructures composed of two-dimensional (2D) materials are already opening many new possibilities in such fields of technology as electronics and magnonics, but far more could be achieved if the number and diversity of 2D materials is increased. So far, only a few dozen 2D crystals have been extracted from materials that exhibit a layered phase in ambient conditions, omitting entirely the large number of layered materials that may exist in other temperatures and pressures. Here, we demonstrate how these structures can be stabilized in 2D van der Waals stacks under room temperature via growing them directly in graphene encapsulation by using graphene oxide as the template material. Specifically, we produce an ambient stable 2D structure of copper and iodine, a material that normally only occurs in layered form at elevated temperatures between 645 and 675 K. Our results establish a route to the production of more exotic phases of materials that would otherwise be difficult or impossible to stabilize for experiments in ambient.
△ Less
Submitted 3 December, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
Picometer-precision few-tilt ptychotomography of 2D materials
Authors:
Christoph Hofer,
Kimmo Mustonen,
Viera Skakalova,
Timothy J. Pennycook
Abstract:
From ripples to defects, edges and grain boundaries, the 3D atomic structure of 2D materials is critical to their properties. However the damage inflicted by conventional 3D analysis precludes its use with fragile 2D materials, particularly for the analysis of local defects. Here we dramatically increase the potential for precise local 3D atomic structure analysis of 2D materials, with both greatl…
▽ More
From ripples to defects, edges and grain boundaries, the 3D atomic structure of 2D materials is critical to their properties. However the damage inflicted by conventional 3D analysis precludes its use with fragile 2D materials, particularly for the analysis of local defects. Here we dramatically increase the potential for precise local 3D atomic structure analysis of 2D materials, with both greatly improved dose efficiency and sensitivity to light elements. We demonstrate light atoms can now be located in complex 2D materials with picometer precision at doses 30 times lower than previously possible. Moreover we demonstrate this using a material, WS$_2$, in which the light S atoms are in fact practically invisible to conventional methods. The key advance is combining the concept of few tilt tomography with highly dose efficient ptychography in scanning transmission electron microscopy. We further demonstrate the method experimentally with the even more challenging and newly discovered 2D CuI, leveraging a new extremely high temporal resolution camera.
△ Less
Submitted 21 June, 2022; v1 submitted 10 August, 2021;
originally announced August 2021.
-
Topological Attention for Time Series Forecasting
Authors:
Sebastian Zeng,
Florian Graf,
Christoph Hofer,
Roland Kwitt
Abstract:
The problem of (point) forecasting $ \textit{univariate} $ time series is considered. Most approaches, ranging from traditional statistical methods to recent learning-based techniques with neural networks, directly operate on raw time series observations. As an extension, we study whether $\textit{local topological properties}$, as captured via persistent homology, can serve as a reliable signal t…
▽ More
The problem of (point) forecasting $ \textit{univariate} $ time series is considered. Most approaches, ranging from traditional statistical methods to recent learning-based techniques with neural networks, directly operate on raw time series observations. As an extension, we study whether $\textit{local topological properties}$, as captured via persistent homology, can serve as a reliable signal that provides complementary information for learning to forecast. To this end, we propose $\textit{topological attention}$, which allows attending to local topological features within a time horizon of historical data. Our approach easily integrates into existing end-to-end trainable forecasting models, such as $\texttt{N-BEATS}$, and in combination with the latter exhibits state-of-the-art performance on the large-scale M4 benchmark dataset of 100,000 diverse time series from different domains. Ablation experiments, as well as a comparison to a broad range of forecasting methods in a setting where only a single time series is available for training, corroborate the beneficial nature of including local topological information through an attention mechanism.
△ Less
Submitted 19 July, 2021;
originally announced July 2021.
-
Event driven 4D STEM acquisition with a Timepix3 detector: microsecond dwell time and faster scans for high precision and low dose applications
Authors:
Daen Jannis,
Christoph Hofer,
Chuang Gao,
Xiaobin Xie,
Armand Béché,
Timothy J. Pennycook,
Jo Verbeeck
Abstract:
Four dimensional scanning transmission electron microscopy (4D STEM) records the scattering of electrons in a material in great detail. The benefits offered by 4D STEM are substantial, with the wealth of data it provides facilitating for instance high precision, high electron dose efficiency phase imaging via center of mass or ptychography based analysis. However the requirement for a 2D image of…
▽ More
Four dimensional scanning transmission electron microscopy (4D STEM) records the scattering of electrons in a material in great detail. The benefits offered by 4D STEM are substantial, with the wealth of data it provides facilitating for instance high precision, high electron dose efficiency phase imaging via center of mass or ptychography based analysis. However the requirement for a 2D image of the scattering to be recorded at each probe position has long placed a severe bottleneck on the speed at which 4D STEM can be performed. Recent advances in camera technology have greatly reduced this bottleneck, with the detection efficiency of direct electron detectors being especially well suited to the technique. However even the fastest frame driven pixelated detectors still significantly limit the scan speed which can be used in 4D STEM, making the resulting data susceptible to drift and hampering its use for low dose beam sensitive applications. Here we report the development of the use of an event driven Timepix3 direct electron camera that allows us to overcome this bottleneck and achieve 4D STEM dwell times down to 100~ns; orders of magnitude faster than what has been possible with frame based readout. We characterise the detector for different acceleration voltages and show that the method is especially well suited for low dose imaging and promises rich datasets without compromising dwell time when compared to conventional STEM imaging.
△ Less
Submitted 8 December, 2021; v1 submitted 6 July, 2021;
originally announced July 2021.
-
Dissecting Supervised Contrastive Learning
Authors:
Florian Graf,
Christoph D. Hofer,
Marc Niethammer,
Roland Kwitt
Abstract:
Minimizing cross-entropy over the softmax scores of a linear map composed with a high-capacity encoder is arguably the most popular choice for training neural networks on supervised learning tasks. However, recent works show that one can directly optimize the encoder instead, to obtain equally (or even more) discriminative representations via a supervised variant of a contrastive objective. In thi…
▽ More
Minimizing cross-entropy over the softmax scores of a linear map composed with a high-capacity encoder is arguably the most popular choice for training neural networks on supervised learning tasks. However, recent works show that one can directly optimize the encoder instead, to obtain equally (or even more) discriminative representations via a supervised variant of a contrastive objective. In this work, we address the question whether there are fundamental differences in the sought-for representation geometry in the output space of the encoder at minimal loss. Specifically, we prove, under mild assumptions, that both losses attain their minimum once the representations of each class collapse to the vertices of a regular simplex, inscribed in a hypersphere. We provide empirical evidence that this configuration is attained in practice and that reaching a close-to-optimal state typically indicates good generalization performance. Yet, the two losses show remarkably different optimization behavior. The number of iterations required to perfectly fit to data scales superlinearly with the amount of randomly flipped labels for the supervised contrastive loss. This is in contrast to the approximately linear scaling previously reported for networks trained with cross-entropy.
△ Less
Submitted 2 March, 2023; v1 submitted 17 February, 2021;
originally announced February 2021.
-
The CHIME Pulsar Project: System Overview
Authors:
CHIME/Pulsar Collaboration,
M. Amiri,
K. M. Bandura,
P. J. Boyle,
C. Brar,
J. F. Cliche,
K. Crowter,
D. Cubranic,
P. B. Demorest,
N. T. Denman,
M. Dobbs,
F. Q. Dong,
M. Fandino,
E. Fonseca,
D. C. Good,
M. Halpern,
A. S. Hill,
C. Höfer,
V. M. Kaspi,
T. L. Landecker,
C. Leung,
H. -H. Lin,
J. Luo,
K. W. Masui,
J. W. McKee
, et al. (20 additional authors not shown)
Abstract:
We present the design, implementation and performance of a digital backend constructed for the Canadian Hydrogen Intensity Mapping Experiment (CHIME) that uses accelerated computing to observe radio pulsars and transient radio sources. When operating, the CHIME correlator outputs 10 independent streams of beamformed data for the CHIME/Pulsar backend that digitally track specified celestial positio…
▽ More
We present the design, implementation and performance of a digital backend constructed for the Canadian Hydrogen Intensity Mapping Experiment (CHIME) that uses accelerated computing to observe radio pulsars and transient radio sources. When operating, the CHIME correlator outputs 10 independent streams of beamformed data for the CHIME/Pulsar backend that digitally track specified celestial positions. Each of these independent streams are processed by the CHIME/Pulsar backend system which can coherently dedisperse, in real-time, up to dispersion measure values of 2500 pc/cm$^{-3}$ . The tracking beams and real-time analysis system are autonomously controlled by a priority-based algorithm that schedules both known sources and positions of interest for observation with observing cadences as small as one day. Given the distribution of known pulsars and radio-transient sources, the CHIME/Pulsar system can monitor up to 900 positions once per sidereal day and observe all sources with declinations greater than $-20^\circ$ once every $\sim$2 weeks. We also discuss the science program enabled through the current modes of data acquisition for CHIME/Pulsar that centers on timing and searching experiments.
△ Less
Submitted 10 June, 2021; v1 submitted 13 August, 2020;
originally announced August 2020.
-
Parameterized Complexity of Min-Power Asymmetric Connectivity
Authors:
Matthias Bentert,
Roman Haag,
Christian Hofer,
Tomohiro Koana,
André Nichterlein
Abstract:
We investigate parameterized algorithms for the NP-hard problem Min-Power Asymmetric Connectivity (MinPAC) that has applications in wireless sensor networks. Given a directed arc-weighted graph, MinPAC asks for a strongly connected spanning subgraph minimizing the summed vertex costs. Here, the cost of each vertex is the weight of its heaviest outgoing arc in the chosen subgraph. We present linear…
▽ More
We investigate parameterized algorithms for the NP-hard problem Min-Power Asymmetric Connectivity (MinPAC) that has applications in wireless sensor networks. Given a directed arc-weighted graph, MinPAC asks for a strongly connected spanning subgraph minimizing the summed vertex costs. Here, the cost of each vertex is the weight of its heaviest outgoing arc in the chosen subgraph. We present linear-time algorithms for the cases where the number of strongly connected components in a so-called obligatory subgraph or the feedback edge number in the underlying undirected graph is constant. Complementing these results, we prove that the problem is W[2]-hard with respect to the solution cost, even on restricted graphs with one feedback arc and binary arc weights.
△ Less
Submitted 29 May, 2020;
originally announced May 2020.
-
A bright millisecond-duration radio burst from a Galactic magnetar
Authors:
The CHIME/FRB Collaboration,
:,
B. C. Andersen,
K. M. Bandura,
M. Bhardwaj,
A. Bij,
M. M. Boyce,
P. J. Boyle,
C. Brar,
T. Cassanelli,
P. Chawla,
T. Chen,
J. -F. Cliche,
A. Cook,
D. Cubranic,
A. P. Curtin,
N. T. Denman,
M. Dobbs,
F. Q. Dong,
M. Fandino,
E. Fonseca,
B. M. Gaensler,
U. Giri,
D. C. Good,
M. Halpern
, et al. (47 additional authors not shown)
Abstract:
Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen…
▽ More
Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.
△ Less
Submitted 15 June, 2020; v1 submitted 20 May, 2020;
originally announced May 2020.
-
Topologically Densified Distributions
Authors:
Christoph D. Hofer,
Florian Graf,
Marc Niethammer,
Roland Kwitt
Abstract:
We study regularization in the context of small sample-size learning with over-parameterized neural networks. Specifically, we shift focus from architectural properties, such as norms on the network weights, to properties of the internal representations before a linear classifier. Specifically, we impose a topological constraint on samples drawn from the probability measure induced in that space.…
▽ More
We study regularization in the context of small sample-size learning with over-parameterized neural networks. Specifically, we shift focus from architectural properties, such as norms on the network weights, to properties of the internal representations before a linear classifier. Specifically, we impose a topological constraint on samples drawn from the probability measure induced in that space. This provably leads to mass concentration effects around the representations of training instances, i.e., a property beneficial for generalization. By leveraging previous work to impose topological constraints in a neural network setting, we provide empirical evidence (across various vision benchmarks) to support our claim for better generalization.
△ Less
Submitted 17 May, 2021; v1 submitted 12 February, 2020;
originally announced February 2020.
-
Periodic activity from a fast radio burst source
Authors:
The CHIME/FRB Collaboration,
M. Amiri,
B. C. Andersen,
K. M. Bandura,
M. Bhardwaj,
P. J. Boyle,
C. Brar,
P. Chawla,
T. Chen,
J. F. Cliche,
D. Cubranic,
M. Deng,
N. T. Denman,
M. Dobbs,
F. Q. Dong,
M. Fandino,
E. Fonseca,
B. M. Gaensler,
U. Giri,
D. C. Good,
M. Halpern,
J. W. T. Hessels,
A. S. Hill,
C. Höfer,
A. Josephy
, et al. (48 additional authors not shown)
Abstract:
Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadicall…
▽ More
Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.
△ Less
Submitted 18 June, 2020; v1 submitted 28 January, 2020;
originally announced January 2020.
-
CHIME/FRB Detection of Eight New Repeating Fast Radio Burst Sources
Authors:
The CHIME/FRB Collaboration,
:,
B. C. Andersen,
K. Bandura,
M. Bhardwaj,
P. Boubel,
M. M. Boyce,
P. J. Boyle,
C. Brar,
T. Cassanelli,
P. Chawla,
D. Cubranic,
M. Deng,
M. Dobbs,
M. Fandino,
E. Fonseca,
B. M. Gaensler,
A. J. Gilbert,
U. Giri,
D. C. Good,
M. Halpern,
A. S. Hill,
G. Hinshaw,
C. Höfer,
A. Josephy
, et al. (33 additional authors not shown)
Abstract:
We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent…
▽ More
We report on the discovery of eight repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 103.5 to 1281 pc cm$^{-3}$. They display varying degrees of activity: six sources were detected twice, another three times, and one ten times. These eight repeating FRBs likely represent the bright and/or high-rate end of a distribution of infrequently repeating sources. For all sources, we determine sky coordinates with uncertainties of $\sim$10$^\prime$. FRB 180916.J0158+65 has a burst-averaged DM = $349.2 \pm 0.3$ pc cm$^{-3}$ and a low DM excess over the modelled Galactic maximum (as low as $\sim$20 pc cm$^{-3}$); this source also has a Faraday rotation measure (RM) of $-114.6 \pm 0.6$ rad m$^{-2}$, much lower than the RM measured for FRB 121102. FRB 181030.J1054+73 has the lowest DM for a repeater, $103.5 \pm 0.3$ pc cm$^{-3}$, with a DM excess of $\sim$ 70 pc cm$^{-3}$. Both sources are interesting targets for multi-wavelength follow-up due to their apparent proximity. The DM distribution of our repeater sample is statistically indistinguishable from that of the first 12 CHIME/FRB sources that have not repeated. We find, with 4$σ$ significance, that repeater bursts are generally wider than those of CHIME/FRB bursts that have not repeated, suggesting different emission mechanisms. Our repeater events show complex morphologies that are reminiscent of the first two discovered repeating FRBs. The repetitive behavior of these sources will enable interferometric localizations and subsequent host galaxy identifications.
△ Less
Submitted 21 October, 2019; v1 submitted 9 August, 2019;
originally announced August 2019.
-
Connectivity-Optimized Representation Learning via Persistent Homology
Authors:
Christoph Hofer,
Roland Kwitt,
Mandar Dixit,
Marc Niethammer
Abstract:
We study the problem of learning representations with controllable connectivity properties. This is beneficial in situations when the imposed structure can be leveraged upstream. In particular, we control the connectivity of an autoencoder's latent space via a novel type of loss, operating on information from persistent homology. Under mild conditions, this loss is differentiable and we present a…
▽ More
We study the problem of learning representations with controllable connectivity properties. This is beneficial in situations when the imposed structure can be leveraged upstream. In particular, we control the connectivity of an autoencoder's latent space via a novel type of loss, operating on information from persistent homology. Under mild conditions, this loss is differentiable and we present a theoretical analysis of the properties induced by the loss. We choose one-class learning as our upstream task and demonstrate that the imposed structure enables informed parameter selection for modeling the in-class distribution via kernel density estimators. Evaluated on computer vision data, these one-class models exhibit competitive performance and, in a low sample size regime, outperform other methods by a large margin. Notably, our results indicate that a single autoencoder, trained on auxiliary (unlabeled) data, yields a mapping into latent space that can be reused across datasets for one-class learning.
△ Less
Submitted 21 June, 2019;
originally announced June 2019.
-
Graph Filtration Learning
Authors:
Christoph D. Hofer,
Florian Graf,
Bastian Rieck,
Marc Niethammer,
Roland Kwitt
Abstract:
We propose an approach to learning with graph-structured data in the problem domain of graph classification. In particular, we present a novel type of readout operation to aggregate node features into a graph-level representation. To this end, we leverage persistent homology computed via a real-valued, learnable, filter function. We establish the theoretical foundation for differentiating through…
▽ More
We propose an approach to learning with graph-structured data in the problem domain of graph classification. In particular, we present a novel type of readout operation to aggregate node features into a graph-level representation. To this end, we leverage persistent homology computed via a real-valued, learnable, filter function. We establish the theoretical foundation for differentiating through the persistent homology computation. Empirically, we show that this type of readout operation compares favorably to previous techniques, especially when the graph connectivity structure is informative for the learning problem.
△ Less
Submitted 17 May, 2021; v1 submitted 27 May, 2019;
originally announced May 2019.
-
High-resolution remote thermography using luminescent low-dimensional tin-halide perovskites
Authors:
Sergii Yakunin,
Bogdan M. Benin,
Yevhen Shynkarenko,
Olga Nazarenko,
Maryna I. Bodnarchuk,
Dmitry N. Dirin,
Christoph Hofer,
Stefano Cattaneo,
Maksym V. Kovalenko
Abstract:
While metal-halide perovskites have recently revolutionized research in optoelectronics through a unique combination of performance and synthetic simplicity, their low-dimensional counterparts can further expand the field with hitherto unknown and practically useful optical functionalities. In this context, we present the strong temperature dependence of the photoluminescence (PL) lifetime of low-…
▽ More
While metal-halide perovskites have recently revolutionized research in optoelectronics through a unique combination of performance and synthetic simplicity, their low-dimensional counterparts can further expand the field with hitherto unknown and practically useful optical functionalities. In this context, we present the strong temperature dependence of the photoluminescence (PL) lifetime of low-dimensional, perovskite-like tin-halides, and apply this property to thermal imaging with a high precision of 0.05 °C. The PL lifetimes are governed by the heat-assisted de-trapping of self-trapped excitons, and their values can be varied over several orders of magnitude by adjusting the temperature (up to 20 ns °C-1). Typically, this sensitive range spans up to one hundred centigrade, and it is both compound-specific and shown to be compositionally and structurally tunable from -100 to 110 ° C going from [C(NH2)3]2SnBr4 to Cs4SnBr6 and (C4N2H14I)4SnI6. Finally, through the innovative implementation of cost-effective hardware for fluorescence lifetime imaging (FLI), based on time-of-flight (ToF) technology, these novel thermoluminophores have been used to record thermographic videos with high spatial and thermal resolution.
△ Less
Submitted 21 May, 2019;
originally announced May 2019.
-
The Metaphysics of Machian Frame-Dragging
Authors:
Antonio Vassallo,
Carl Hoefer
Abstract:
The paper investigates the kind of dependence relation that best portrays Machian frame-dragging in general relativity. The question is tricky because frame-dragging relates local inertial frames to distant distributions of matter in a time-independent way, thus establishing some sort of non-local link between the two. For this reason, a plain causal interpretation of frame-dragging faces huge cha…
▽ More
The paper investigates the kind of dependence relation that best portrays Machian frame-dragging in general relativity. The question is tricky because frame-dragging relates local inertial frames to distant distributions of matter in a time-independent way, thus establishing some sort of non-local link between the two. For this reason, a plain causal interpretation of frame-dragging faces huge challenges. The paper will shed light on the issue by using a generalized structural equation model analysis in terms of manipulationist counterfactuals recently applied in the context of metaphysical enquiry by Schaffer (2016) and Wilson (2017). The verdict of the analysis will be that frame-dragging is best understood in terms of a novel type of dependence relation that is half-way between causation and grounding.
△ Less
Submitted 30 January, 2019;
originally announced January 2019.
-
A Second Source of Repeating Fast Radio Bursts
Authors:
The CHIME/FRB Collaboration,
:,
M. Amiri,
K. Bandura,
M. Bhardwaj,
P. Boubel,
M. M. Boyce,
P. J. Boyle,
C. Brar,
M. Burhanpurkar,
T. Cassanelli,
P. Chawla,
J. F. Cliche,
D. Cubranic,
M. Deng,
N. Denman,
M. Dobbs,
M. Fandino,
E. Fonseca,
B. M. Gaensler,
A. J. Gilbert,
A. Gill,
U. Giri,
D. C. Good,
M. Halpern
, et al. (36 additional authors not shown)
Abstract:
The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by…
▽ More
The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
Observations of Fast Radio Bursts at Frequencies down to 400 Megahertz
Authors:
CHIME/FRB Collaboration,
:,
Mandana Amiri,
Kevin Bandura,
Mohit Bhardwaj,
Paula Boubel,
Michelle M. Boyce,
Patrick J. Boyle,
Charanjot Brar,
Maya Burhanpurkar,
Pragya Chawla,
Jean F. Cliche,
Davor Cubranic,
Meiling Deng,
Nolan Denman,
Matthew Dobbs,
M. Fandino,
Emmanuel Fonseca,
Bryan M. Gaensler,
Adam J. Gilbert,
Utkarsh Giri,
Deborah C. Good,
Mark Halpern,
David Hanna,
Alexander S. Hill
, et al. (31 additional authors not shown)
Abstract:
Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low…
▽ More
Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.
△ Less
Submitted 14 January, 2019;
originally announced January 2019.
-
Direct visualization of the 3D structure of silicon impurities in graphene
Authors:
Christoph Hofer,
Viera Skakalova,
Mohammad Reza Ahmadpour Monazam,
Clemens Mangler,
Jani Kotakoski,
Toma Susi,
Jannik C. Meyer
Abstract:
We directly visualize the three-dimensional (3D) geometry and dynamics of silicon impurities in graphene as well as their dynamics by aberration-corrected scanning transmission electron microscopy. By acquiring images when the sample is tilted, we show that an asymmetry of the atomic position of the heteroatom in the projection reveals the non-planarity of the structure. From a sequence of images,…
▽ More
We directly visualize the three-dimensional (3D) geometry and dynamics of silicon impurities in graphene as well as their dynamics by aberration-corrected scanning transmission electron microscopy. By acquiring images when the sample is tilted, we show that an asymmetry of the atomic position of the heteroatom in the projection reveals the non-planarity of the structure. From a sequence of images, we further demonstrate that the Si atom switches between up- and down- configurations with respect to the graphene plane, with an asymmetric cross-section. We further analyze the 3D structure and dynamics of a silicon tetramer in graphene. Our results clarify the out-of-plane structure of impurities in graphene by direct experimental observation and open a new route to study their dynamics in three dimensions.
△ Less
Submitted 24 September, 2018;
originally announced September 2018.
-
Revealing the 3D Structure of Graphene Defects
Authors:
Christoph Hofer,
Christian Kramberger,
Mohammad Reza Ahmadpour Monazam,
Clemens Mangler,
Andreas Mittelberger,
Giacomo Argentero,
Jani Kotakoski,
Jannik C. Meyer
Abstract:
We demonstrate insights into the three-dimensional structure of defects in graphene, in particular grain boundaries, obtained via a new approach from two transmission electron microscopy images recorded at different angles. The structure is obtained through an optimization process where both the atomic positions as well as the simulated imaging parameters are iteratively changed until the best pos…
▽ More
We demonstrate insights into the three-dimensional structure of defects in graphene, in particular grain boundaries, obtained via a new approach from two transmission electron microscopy images recorded at different angles. The structure is obtained through an optimization process where both the atomic positions as well as the simulated imaging parameters are iteratively changed until the best possible match to the experimental images is found. We first demonstrate that this method works using an embedded defect in graphene that allows direct comparison to the computationally predicted three-dimensional shape. We then applied the method to a set of grain boundary structures with misorientation angles nearly spanning the whole available range (2.6-29.8°). The measured height variations at the boundaries reveal a strong correlation with the misorientation angle with lower angles resulting in stronger corrugation and larger kink angles. Our results allow for the first time a direct comparison with theoretical predictions for the corrugation at grain boundaries and we show that the measured kink angles are significantly smaller than the largest predicted ones.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Analysis of point defects in graphene using low dose scanning transmission electron microscopy imaging and maximum likelihood reconstruction
Authors:
Christian Kramberger,
Andreas Mittelberger,
Christoph Hofer,
Jannik C. Meyer
Abstract:
Freestanding graphene displays an outstanding resilience to electron irradiation at low electron energies. Point defects in graphene are, however, subject to beam driven dynamics. This means that high resolution micrographs of point defects, which usually require a high electron irradiation dose might not represent the intrinsic defect population. Here, we capture the inital defects formed by ejec…
▽ More
Freestanding graphene displays an outstanding resilience to electron irradiation at low electron energies. Point defects in graphene are, however, subject to beam driven dynamics. This means that high resolution micrographs of point defects, which usually require a high electron irradiation dose might not represent the intrinsic defect population. Here, we capture the inital defects formed by ejecting carbon atoms under electron irradiation, by imaging with very low doses and subsequent reconstruction of the frequently occuring defects via a maximum likelihood algorithm.
△ Less
Submitted 4 May, 2018;
originally announced May 2018.
-
A parallel multigrid solver for multi-patch Isogeometric Analysis
Authors:
Christoph Hofer,
Stefan Takacs
Abstract:
Isogeometric Analysis (IgA) is a framework for setting up spline-based discretizations of partial differential equations, which has been introduced around a decade ago and has gained much attention since then. If large spline degrees are considered, one obtains the approximation power of a high-order method, but the number of degrees of freedom behaves like for a low-order method. One important in…
▽ More
Isogeometric Analysis (IgA) is a framework for setting up spline-based discretizations of partial differential equations, which has been introduced around a decade ago and has gained much attention since then. If large spline degrees are considered, one obtains the approximation power of a high-order method, but the number of degrees of freedom behaves like for a low-order method. One important ingredient to use a discretization with large spline degree, is a robust and preferably parallelizable solver. While numerical evidence shows that multigrid solvers with standard smoothers (like Gauss Seidel) does not perform well if the spline degree is increased, the multigrid solvers proposed by the authors and their co-workers proved to behave optimal both in the grid size and the spline degree. In the present paper, the authors want to show that those solvers are parallelizable and that they scale well in a parallel environment.
△ Less
Submitted 7 April, 2018;
originally announced April 2018.
-
The CHIME Fast Radio Burst Project: System Overview
Authors:
The CHIME/FRB Collaboration,
:,
M. Amiri,
K. Bandura,
P. Berger,
M. Bhardwaj,
M. M. Boyce,
P. J. Boyle,
C. Brar,
M. Burhanpurkar,
P. Chawla,
J. Chowdhury,
J. F. Cliche,
M. D. Cranmer,
D. Cubranic,
M. Deng,
N. Denman,
M. Dobbs,
M. Fandino,
E. Fonseca,
B. M. Gaensler,
U. Giri,
A. J. Gilbert,
D. C. Good,
S. Guliani
, et al. (28 additional authors not shown)
Abstract:
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful…
▽ More
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a novel transit radio telescope operating across the 400-800-MHz band. CHIME is comprised of four 20-m x 100-m semi-cylindrical paraboloid reflectors, each of which has 256 dual-polarization feeds suspended along its axis, giving it a >200 square degree field-of-view. This, combined with wide bandwidth, high sensitivity, and a powerful correlator makes CHIME an excellent instrument for the detection of Fast Radio Bursts (FRBs). The CHIME Fast Radio Burst Project (CHIME/FRB) will search beam-formed, high time-and frequency-resolution data in real time for FRBs in the CHIME field-of-view. Here we describe the CHIME/FRB backend, including the real-time FRB search and detection software pipeline as well as the planned offline analyses. We estimate a CHIME/FRB detection rate of 2-42 FRBs/sky/day normalizing to the rate estimated at 1.4-GHz by Vander Wiel et al. (2016). Likely science outcomes of CHIME/FRB are also discussed. CHIME/FRB is currently operational in a commissioning phase, with science operations expected to commence in the latter half of 2018.
△ Less
Submitted 29 March, 2018;
originally announced March 2018.
-
Discontinuous Galerkin Isogeometric Analysis for segmentations generating overlapping regions
Authors:
Christoph Hofer,
Ioannis Toulopoulos
Abstract:
In the Isogeometric Analysis (IGA) framework, the computational domain has very often a multipatch representation. The multipatch domain can be obtained by a volume segmentation of a boundary represented domain, e.g., provided by a Computer Aided Design (CAD) model. Typically, small gap and overlapping regions can appear at the patch interfaces of such multipatch representations. In the current wo…
▽ More
In the Isogeometric Analysis (IGA) framework, the computational domain has very often a multipatch representation. The multipatch domain can be obtained by a volume segmentation of a boundary represented domain, e.g., provided by a Computer Aided Design (CAD) model. Typically, small gap and overlapping regions can appear at the patch interfaces of such multipatch representations. In the current work we consider multipatch representations having only small overlapping regions between the patches.
We develop a Discontinuous Galerkin (DG)- IGA method which can be immediately applied to these representations. Our method appropriately connects the fluxes of the one face of the overlapping region with the flux of the opposite face. We provide a theoretical justification of our approach by splitting the whole error into two components: the first is related to the incorrect representation of the patches (consistency error) and the second to the approximation properties of the IGA space. We show bounds for both components of the error. We verify the theoretical error estimates in a series of numerical examples.
△ Less
Submitted 26 March, 2018;
originally announced March 2018.
-
Automated image acquisition for low-dose STEM at atomic resolution
Authors:
Andreas Mittelberger,
Christian Kramberger,
Christoph Hofer,
Clemens Mangler,
Jannik C. Meyer
Abstract:
Beam damage is a major limitation in electron microscopy that becomes increasingly severe at higher resolution. One possible route to circumvent radiation damage, which forms the basis for single-particle electron microscopy and related techniques, is to distribute the dose over many identical copies of an object. For the acquisition of low-dose data, ideally no dose should be applied to the regio…
▽ More
Beam damage is a major limitation in electron microscopy that becomes increasingly severe at higher resolution. One possible route to circumvent radiation damage, which forms the basis for single-particle electron microscopy and related techniques, is to distribute the dose over many identical copies of an object. For the acquisition of low-dose data, ideally no dose should be applied to the region of interest prior to the acquisition of data. We present an automated approach that can collect large amounts of data efficiently by acquiring images in an user-defined area-of-interest with atomic resolution. We demonstrate that the stage mechanics of the Nion UltraSTEM, combined with an intelligent algorithm to move the sample, allows the automated acquisition of atomically resolved images from micron-sized areas of a graphene substrate. Moving the sample stage automatically in a regular pattern over the area-of-interest enables the collection of data from pristine sample regions without exposing them to the electron beam before recording an image. Therefore, it is possible to obtain data with minimal dose (no prior exposure from focusing), which is only limited by the minimum signal needed for data processing. This enables us to prevent beam induced damage in the sample and to acquire large datasets within a reasonable amount of time.
△ Less
Submitted 16 March, 2018; v1 submitted 15 March, 2018;
originally announced March 2018.
-
Robust Preconditioning for Space-Time Isogeometric Analysis of Parabolic Evolution Problems
Authors:
Christoph Hofer,
Ulrich Langer,
Martin Neumüller
Abstract:
We propose and investigate new robust preconditioners for space-time Isogeometric Analysis of parabolic evolution problems. These preconditioners are based on a time parallel multigrid method. We consider a decomposition of the space-time cylinder into time-slabs which are coupled via a discontinuous Galerkin technique. The time-slabs provide the structure for the time-parallel multigrid solver. T…
▽ More
We propose and investigate new robust preconditioners for space-time Isogeometric Analysis of parabolic evolution problems. These preconditioners are based on a time parallel multigrid method. We consider a decomposition of the space-time cylinder into time-slabs which are coupled via a discontinuous Galerkin technique. The time-slabs provide the structure for the time-parallel multigrid solver. The most important part of the multigrid method is the smoother. We utilize the special structure of the involved operator to decouple its application into several spatial problems by means of generalized eigenvalue or Schur decompositions. Some of these problems have a symmetric saddle point structure, for which we present robust preconditions. Finally, we present numerical experiments confirming the robustness of our space-time IgA solver.
△ Less
Submitted 26 February, 2018;
originally announced February 2018.
-
Deep Learning with Topological Signatures
Authors:
Christoph Hofer,
Roland Kwitt,
Marc Niethammer,
Andreas Uhl
Abstract:
Inferring topological and geometrical information from data can offer an alternative perspective on machine learning problems. Methods from topological data analysis, e.g., persistent homology, enable us to obtain such information, typically in the form of summary representations of topological features. However, such topological signatures often come with an unusual structure (e.g., multisets of…
▽ More
Inferring topological and geometrical information from data can offer an alternative perspective on machine learning problems. Methods from topological data analysis, e.g., persistent homology, enable us to obtain such information, typically in the form of summary representations of topological features. However, such topological signatures often come with an unusual structure (e.g., multisets of intervals) that is highly impractical for most machine learning techniques. While many strategies have been proposed to map these topological signatures into machine learning compatible representations, they suffer from being agnostic to the target learning task. In contrast, we propose a technique that enables us to input topological signatures to deep neural networks and learn a task-optimal representation during training. Our approach is realized as a novel input layer with favorable theoretical properties. Classification experiments on 2D object shapes and social network graphs demonstrate the versatility of the approach and, in case of the latter, we even outperform the state-of-the-art by a large margin.
△ Less
Submitted 16 February, 2018; v1 submitted 13 July, 2017;
originally announced July 2017.
-
Inexact Dual-Primal Isogeometric Tearing and Interconnecting Methods
Authors:
Christoph Hofer,
Ulrich Langer,
Stefan Takacs
Abstract:
In this paper, we investigate inexact variants of dual-primal isogeometric tearing and interconnecting methods for solving large-scale systems of linear equations arising from Galerkin isogeometric discretizations of elliptic boundary value problems. The considered methods are extensions of standard finite element tearing and interconnecting methods to isogeometric analysis. The algorithms are imp…
▽ More
In this paper, we investigate inexact variants of dual-primal isogeometric tearing and interconnecting methods for solving large-scale systems of linear equations arising from Galerkin isogeometric discretizations of elliptic boundary value problems. The considered methods are extensions of standard finite element tearing and interconnecting methods to isogeometric analysis. The algorithms are implemented by means of energy minimizing primal subspaces. We discuss the replacement of local sparse direct solvers by iterative methods, particularly, multigrid solvers. We investigate the incorporation of these iterative solvers into different formulations of the algorithm. Finally, we present numerical examples comparing the performance of these inexact versions.
△ Less
Submitted 12 May, 2017;
originally announced May 2017.
-
Limits on the ultra-bright Fast Radio Burst population from the CHIME Pathfinder
Authors:
CHIME Scientific Collaboration,
Mandana Amiri,
Kevin Bandura,
Philippe Berger,
J. Richard Bond,
Jean-François Cliche,
Liam Connor,
Meiling Deng,
Nolan Denman,
Matt Dobbs,
Rachel Simone Domagalski,
Mateus Fandino,
Adam J Gilbert,
Deborah C. Good,
Mark Halpern,
David Hanna,
Adam D. Hincks,
Gary Hinshaw,
Carolin Höfer,
Gilbert Hsyu,
Peter Klages,
T. L. Landecker,
Kiyoshi Masui,
Juan Mena-Parra,
Laura Newburgh
, et al. (13 additional authors not shown)
Abstract:
We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope,…
▽ More
We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope, $α\equiv-\frac{\partial \log N}{\partial \log S}$, is quite small. A 256-input incoherent beamformer was deployed on the CHIME Pathfinder for this purpose. If the FRB distribution were described by a single power-law with $α=0.7$, we would expect an FRB detection every few days, making this the fastest survey on sky at present. We collected 1268 hours of data, amounting to one of the largest exposures of any FRB survey, with over 2.4\,$\times$\,10$^5$\,deg$^2$\,hrs. Having seen no bursts, we have constrained the rate of extremely bright events to $<\!13$\,sky$^{-1}$\,day$^{-1}$ above $\sim$\,220$\sqrt{(τ/\rm ms)}$ Jy\,ms for $τ$ between 1.3 and 100\,ms, at 400--800\,MHz. The non-detection also allows us to rule out $α\lesssim0.9$ with 95$\%$ confidence, after marginalizing over uncertainties in the GBT rate at 700--900\,MHz, though we show that for a cosmological population and a large dynamic range in flux density, $α$ is brightness-dependent. Since FRBs now extend to large enough distances that non-Euclidean effects are significant, there is still expected to be a dearth of faint events and relative excess of bright events. Nevertheless we have constrained the allowed number of ultra-intense FRBs. While this does not have significant implications for deeper, large-FoV surveys like full CHIME and APERTIF, it does have important consequences for other wide-field, small dish experiments.
△ Less
Submitted 20 April, 2017; v1 submitted 26 February, 2017;
originally announced February 2017.
-
Analysis of discontinuous Galerkin dual-primal isogeometric tearing and interconnecting methods
Authors:
Christoph Hofer
Abstract:
In this paper, we present the analysis of the discontinuous Galerkin dual-primal isogeometric tearing and interconnecting method (dG-IETI-DP) for the two-dimensional case where we only consider vertex primal variables. The dG-IETI-DP method is a combination of the dual-primal isogeometric tearing and interconnecting method (IETI-DP) with the discontinuous Galerkin (dG) method. We use the dG method…
▽ More
In this paper, we present the analysis of the discontinuous Galerkin dual-primal isogeometric tearing and interconnecting method (dG-IETI-DP) for the two-dimensional case where we only consider vertex primal variables. The dG-IETI-DP method is a combination of the dual-primal isogeometric tearing and interconnecting method (IETI-DP) with the discontinuous Galerkin (dG) method. We use the dG method only on the interfaces to couple different patches. This enables us to handle non-matching grids on patch interfaces as well as segmentation crimes (gaps and overlaps) between the patches. The purpose of this paper is to derive quasi-optimal bounds for the condition number of the preconditioned system with respect to the maximal ratio $H/h:=\max_k(H_k/h_k)$ of subdomain diameter and meshsize. We show that the constant is independent of $h_k$ and $H_k$, but depends on the ratio of meshsizes of neighbouring patches $h_\ell/h_k$.
△ Less
Submitted 30 November, 2016;
originally announced November 2016.
-
Parallelization of continuous and discontinuous Galerkin dual-primal Isogeometric tearing and interconnecting methods
Authors:
Christoph Hofer
Abstract:
In this paper we investigate the parallelization of dual-primal isogeometric tearing and interconnecting (IETI-DP) type methods for solving large-scale continuous and discontinuous Galerkin systems of equations arising from Isogeometric analysis of elliptic boundary value problems. These methods are extensions of the finite element tearing and interconnecting methods to isogeometric analysis. The…
▽ More
In this paper we investigate the parallelization of dual-primal isogeometric tearing and interconnecting (IETI-DP) type methods for solving large-scale continuous and discontinuous Galerkin systems of equations arising from Isogeometric analysis of elliptic boundary value problems. These methods are extensions of the finite element tearing and interconnecting methods to isogeometric analysis. The algorithms are implemented by means of energy minimizing primal subspaces. We discuss how these methods can efficiently be parallelized in a distributed memory setting. Weak and strong scaling studies presented for two and three dimensional problems show an excellent parallel efficiency.
△ Less
Submitted 24 November, 2016;
originally announced November 2016.
-
Discontinuous Galerkin Isogeometric Analysis on Non-matching Segmentation: Error Estimates and Efficient Solvers
Authors:
Christoph Hofer,
Ulrich Langer,
Ioannis Toulopoulos
Abstract:
The Isogeometric Analysis (IgA) of boundary value problems in complex domains often requires a decomposition of the computational domain into patches such that each of which can be parametrized by the so-called geometrical mapping. In this paper, we develop discontinuous Galerkin (dG) IgA techniques for solving elliptic diffusion problems on decompositions that can include non-matching parametriza…
▽ More
The Isogeometric Analysis (IgA) of boundary value problems in complex domains often requires a decomposition of the computational domain into patches such that each of which can be parametrized by the so-called geometrical mapping. In this paper, we develop discontinuous Galerkin (dG) IgA techniques for solving elliptic diffusion problems on decompositions that can include non-matching parametrizations of the interfaces, i.e., the interfaces of the adjacent patches may be not identical. The lack of the exact parametrization of the patches leads to the creation of gap and overlapping regions between the patches. This does not allow the immediate use of the classical numerical fluxes that are known in the literature. The unknown normal fluxes of the solution on the non-matching interfaces are approximated by Taylor expansions using the values of the solution computed on the boundary of the patches These approximations are used in order to build up the numerical fluxes of the final dG IgA scheme and to couple the local patch-wise discrete problems. The resulting linear systems are solved by using efficient domainecomposition methods based on the tearing and interconnecting technology. We present numerical results of a series of test problems that validate the theoretical estimates presented.
△ Less
Submitted 12 October, 2016;
originally announced October 2016.
-
Isotope analysis in the transmission electron microscope
Authors:
Toma Susi,
Christoph Hofer,
Giacomo Argentero,
Gregor T. Leuthner,
Timothy J. Pennycook,
Clemens Mangler,
Jannik C. Meyer,
Jani Kotakoski
Abstract:
The Ångström-sized probe of the scanning transmission electron microscope can visualize and collect spectra from single atoms. This can unambiguously resolve the chemical structure of materials, but not their isotopic composition. Here we differentiate between two isotopes of the same element by quantifying how likely the energetic imaging electrons are to eject atoms. First, we measure the displa…
▽ More
The Ångström-sized probe of the scanning transmission electron microscope can visualize and collect spectra from single atoms. This can unambiguously resolve the chemical structure of materials, but not their isotopic composition. Here we differentiate between two isotopes of the same element by quantifying how likely the energetic imaging electrons are to eject atoms. First, we measure the displacement probability in graphene grown from either $^{12}$C or $^{13}$C and describe the process using a quantum mechanical model of lattice vibrations coupled with density functional theory simulations. We then test our spatial resolution in a mixed sample by ejecting individual atoms from nanoscale areas spanning an interface region that is far from atomically sharp, mapping the isotope concentration with a precision better than 20%. Although we use a scanning instrument, our method should be applicable to any atomic resolution transmission electron microscope and to other low-dimensional materials.
△ Less
Submitted 16 August, 2016;
originally announced August 2016.
-
Holographic Beam Mapping of the CHIME Pathfinder Array
Authors:
Philippe Berger,
Laura B. Newburgh,
Mandana Amiri,
Kevin Bandura,
Jean-Francois Cliche,
Liam Connor,
Meiling Deng,
Nolan Denman,
Matt Dobbs,
Mateus Fandino,
Adam J. Gilbert,
Deborah Good,
Mark Halpern,
David Hanna,
Adam D. Hincks,
Gary Hinshaw,
Carolin Hofer,
Andre M. Johnson,
Tom L. Landecker,
Kiyoshi W. Masui,
Juan Mena Parra,
Niels Oppermann,
Ue-Li Pen,
Jeffrey B. Peterson,
Andre Recnik
, et al. (10 additional authors not shown)
Abstract:
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from th…
▽ More
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.
△ Less
Submitted 5 July, 2016;
originally announced July 2016.
-
Dual-Primal Isogeometric Tearing and Interconnecting solvers for multipatch dG-IgA equations
Authors:
Christoph Hofer,
Ulrich Langer
Abstract:
In this paper we consider a new version of the dual-primal isogeometric tearing and interconnecting (IETI-DP) method for solving large-scale linear systems of algebraic equations arising from discontinuous Galerkin (dG) isogeometric analysis of diffusion problems on multipatch domains with non-matching meshes. The dG formulation is used to couple the local problems across patch interfaces. The pur…
▽ More
In this paper we consider a new version of the dual-primal isogeometric tearing and interconnecting (IETI-DP) method for solving large-scale linear systems of algebraic equations arising from discontinuous Galerkin (dG) isogeometric analysis of diffusion problems on multipatch domains with non-matching meshes. The dG formulation is used to couple the local problems across patch interfaces. The purpose of this paper is to present this new method and provide numerical examples indicating a polylogarithmic condition number bound for the preconditioned system and showing an incredible robustness with respect to large jumps in the diffusion coefficient across the interfaces.
△ Less
Submitted 8 January, 2016;
originally announced January 2016.