-
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.
-
Interactive Treatment Planning in High Dose-Rate Brachytherapy for Gynecological Cancer
Authors:
Huan Liu,
Chang M Ma,
Xun Jia,
Chenyang Shen,
Peter Klages,
Kevin Albuquerque
Abstract:
High dose-rate brachytherapy (HDRBT) is widely used for gynecological cancer treatment. Although commercial treatment planning systems (TPSs) have inverse optimization modules, it takes several iterations to adjust planning objectives to achieve a satisfactory plan. Interactive plan-modification modules enable modifying the plan and visualizing results in real time, but they update plans based on…
▽ More
High dose-rate brachytherapy (HDRBT) is widely used for gynecological cancer treatment. Although commercial treatment planning systems (TPSs) have inverse optimization modules, it takes several iterations to adjust planning objectives to achieve a satisfactory plan. Interactive plan-modification modules enable modifying the plan and visualizing results in real time, but they update plans based on simple geometrical or heuristic algorithms, which cannot ensure resulting plan optimality. This project develops an interactive plan optimization module for HDRBT of gynecological cancer. By efficiently solving an optimization problem in real time, it allows a user to visualize a plan and interactively modify it to improve quality. We formulated an optimization problem with an objective function containing a weighted sum of doses to normal organs subject to user-specified target coverage. A user interface was developed that allows a user to adjust organ weights using scroll bars. With a simple mouse click, the optimization problem is solved in seconds with a highly efficient alternating-direction method of multipliers and a warm start optimization strategy. Resulting clinically relevant D2cc of organs are displayed immediately. This allows a user to intuitively adjust plans with satisfactory quality. We tested the effectiveness of our development in cervix cancer cases treated with a tandem-and-ovoid applicator. It took a maximum of 3 seconds to solve the optimization problem in each instance. With interactive optimization capability, a satisfactory plan can be obtained in <1 min. In our clinic, although the time for plan adjustment was typically <10min with simple interactive plan modification tools in TPS, the resulting plans do not ensure optimality. Our plans achieved on average 5% lower D2cc than clinical plans, while maintaining target coverage.
△ Less
Submitted 10 September, 2021;
originally announced September 2021.
-
Comparison of Patch-Based Conditional Generative Adversarial Neural Net Models with Emphasis on Model Robustness for Use in Head and Neck Cases for MR-Only planning
Authors:
Peter Klages,
Ilyes Benslimane,
Sadegh Riyahi,
Jue Jiang,
Margie Hunt,
Joe Deasy,
Harini Veeraraghavan,
Neelam Tyagi
Abstract:
A total of twenty paired CT and MR images were used in this study to investigate two conditional generative adversarial networks, Pix2Pix, and Cycle GAN, for generating synthetic CT images for Headand Neck cancer cases. Ten of the patient cases were used for training and included such common artifacts as dental implants; the remaining ten testing cases were used for testing and included a larger r…
▽ More
A total of twenty paired CT and MR images were used in this study to investigate two conditional generative adversarial networks, Pix2Pix, and Cycle GAN, for generating synthetic CT images for Headand Neck cancer cases. Ten of the patient cases were used for training and included such common artifacts as dental implants; the remaining ten testing cases were used for testing and included a larger range of image features commonly found in clinical head and neck cases. These features included strong metal artifacts from dental implants, one case with a metal implant, and one case with abnormal anatomy. The original CT images were deformably registered to the mDixon FFE MR images to minimize the effects of processing the MR images. The sCT generation accuracy and robustness were evaluated using Mean Absolute Error (MAE) based on the Hounsfield Units (HU) for three regions (whole body, bone, and air within the body), Mean Error (ME) to observe systematic average offset errors in the sCT generation, and dosimetric evaluation of all clinically relevant structures. For the test set the MAE for the Pix2Pix and Cycle GAN models were 92.4 $\pm$ 13.5 HU, and 100.7 $\pm$ 14.6 HU, respectively, for the body region, 166.3 $\pm$ 31.8 HU, and 184 $\pm$ 31.9 HU, respectively, for the bone region, and 183.7 $\pm$ 41.3 HU and 185.4 $\pm$ 37.9 HU for the air regions. The ME for Pix2Pix and Cycle GAN were 21.0 $\pm$ 11.8 HU and 37.5 $\pm$ 14.9 HU, respectively. Absolute Percent Mean/Max Dose Errors were less than 2% for the PTV and all critical structures for both models, and DRRs generated from these models looked qualitatively similar to CT generated DRRs showing these methods are promising for MR-only planning.
△ Less
Submitted 27 February, 2019; v1 submitted 1 February, 2019;
originally announced February 2019.
-
Intelligent Inverse Treatment Planning via Deep Reinforcement Learning, a Proof-of-Principle Study in High Dose-rate Brachytherapy for Cervical Cancer
Authors:
Chenyang Shen,
Yesenia Gonzalez,
Peter Klages,
Nan Qin,
Hyunuk Jung,
Liyuan Chen,
Dan Nguyen,
Steve B. Jiang,
Xun Jia
Abstract:
Inverse treatment planning in radiation therapy is formulated as optimization problems. The objective function and constraints consist of multiple terms designed for different clinical and practical considerations. Weighting factors of these terms are needed to define the optimization problem. While a treatment planning system can solve the optimization problem with given weights, adjusting the we…
▽ More
Inverse treatment planning in radiation therapy is formulated as optimization problems. The objective function and constraints consist of multiple terms designed for different clinical and practical considerations. Weighting factors of these terms are needed to define the optimization problem. While a treatment planning system can solve the optimization problem with given weights, adjusting the weights for high plan quality is performed by human. The weight tuning task is labor intensive, time consuming, and it critically affects the final plan quality. An automatic weight-tuning approach is strongly desired. The weight tuning procedure is essentially a decision making problem. Motivated by the tremendous success in deep learning for decision making with human-level intelligence, we propose a novel framework to tune the weights in a human-like manner. Using treatment planning in high-dose-rate brachytherapy as an example, we develop a weight tuning policy network (WTPN) that observes dose volume histograms of a plan and outputs an action to adjust organ weights, similar to the behaviors of a human planner. We train the WTPN via end-to-end deep reinforcement learning. Experience replay is performed with the epsilon greedy algorithm. Then we apply the trained WTPN to guide treatment planning of testing patient cases. The trained WTPN successfully learns the treatment planning goals to guide the weight tuning process. On average, the quality score of plans generated under the WTPN's guidance is improved by ~8.5% compared to the initial plan with arbitrary weights, and by 10.7% compared to the plans generated by human planners. To our knowledge, this is the first tool to adjust weights for the treatment planning in a human-like fashion based on learnt intelligence. The study demonstrates potential feasibility to develop intelligent treatment planning system via deep reinforcement learning.
△ Less
Submitted 25 November, 2018;
originally announced November 2018.
-
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.
-
CHIME FRB: An application of FFT beamforming for a radio telescope
Authors:
C. Ng,
K. Vanderlinde,
A. Paradise,
P. Klages,
K. Masui,
K. Smith,
K. Bandura,
P. J. Boyle,
M. Dobbs,
V. Kaspi,
A. Renard,
J. R. Shaw,
I. Stairs,
I. Tretyakov
Abstract:
We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four for…
▽ More
We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four formed beams along East-West via exact phasing, tiling a sky area of ~250 square degrees. A zero-padding approximation is employed to improve chromatic beam alignment across the wide bandwidth of 400 to 800 MHz. We up-channelize the data in order to achieve fine spectral resolution of $Δν$=24 kHz and time cadence of 0.983 ms, desirable for detecting transient and dispersed signals such as those from FRBs.
△ Less
Submitted 15 February, 2017;
originally announced February 2017.
-
GPU Kernels for High-Speed 4-Bit Astrophysical Data Processing
Authors:
Peter Klages,
Kevin Bandura,
Nolan Denman,
Andre Recnik,
Jonathan Sievers,
Keith Vanderlinde
Abstract:
Interferometric radio telescopes often rely on computationally expensive O(N^2) correlation calculations; fortunately these computations map well to massively parallel accelerators such as low-cost GPUs. This paper describes the OpenCL kernels developed for the GPU based X-engine of a new hybrid FX correlator. Channelized data from the F-engine is supplied to the GPUs as 4-bit, offset-encoded real…
▽ More
Interferometric radio telescopes often rely on computationally expensive O(N^2) correlation calculations; fortunately these computations map well to massively parallel accelerators such as low-cost GPUs. This paper describes the OpenCL kernels developed for the GPU based X-engine of a new hybrid FX correlator. Channelized data from the F-engine is supplied to the GPUs as 4-bit, offset-encoded real and imaginary integers. Because of the low bit width of the data, two values may be packed into a 32-bit register, allowing multiplication and addition of more than one value with a single fused multiply-add instruction. With this data and calculation packing scheme, as many as 5.6 effective tera-operations per second (TOPS) can be executed on a 4.3 TOPS GPU. The kernel design allows correlations to scale to large numbers of input elements, limited only by maximum buffer sizes on the GPU. This code is currently working on-sky with the CHIME Pathfinder Correlator in BC, Canada.
△ Less
Submitted 20 March, 2015;
originally announced March 2015.
-
A GPU-based Correlator X-engine Implemented on the CHIME Pathfinder
Authors:
Nolan Denman,
Mandana Amiri,
Kevin Bandura,
Jean-François Cliche,
Liam Connor,
Matt Dobbs,
Mateus Fandino,
Mark Halpern,
Adam Hincks,
Gary Hinshaw,
Carolin Höfer,
Peter Klages,
Kiyoshi Masui,
Juan Mena Parra,
Laura Newburgh,
Andre Recnik,
J. Richard Shaw,
Kris Sigurdson,
Kendrick Smith,
Keith Vanderlinde
Abstract:
We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction o…
▽ More
We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today.
△ Less
Submitted 11 June, 2015; v1 submitted 20 March, 2015;
originally announced March 2015.
-
An Efficient Real-time Data Pipeline for the CHIME Pathfinder Radio Telescope X-Engine
Authors:
Andre Recnik,
Kevin Bandura,
Nolan Denman,
Adam D. Hincks,
Gary Hinshaw,
Peter Klages,
Ue-Li Pen,
Keith Vanderlinde
Abstract:
The CHIME Pathfinder is a new interferometric radio telescope that uses a hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs across a 400 MHz radio band. A software framework is presented to manage this real-time data flow, which allows each of 16 processing servers to handle 51.2 Gb/s of astronomi…
▽ More
The CHIME Pathfinder is a new interferometric radio telescope that uses a hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs across a 400 MHz radio band. A software framework is presented to manage this real-time data flow, which allows each of 16 processing servers to handle 51.2 Gb/s of astronomical data, plus 8 Gb/s of ancillary data. Each server receives data in the form of UDP packets from an FPGA F-engine over the eight 10 GbE links, combines data from these packets into large (32MB-256MB) buffered frames, and transfers them to multiple GPU co-processors for correlation. The results from the GPUs are combined and normalized, then transmitted to a collection server, where they are merged into a single file. Aggressive optimizations enable each server to handle this high rate of data; allowing the efficient correlation of 25 MHz of radio bandwidth per server. The solution scales well to larger values of N by adding additional servers.
△ Less
Submitted 15 June, 2015; v1 submitted 20 March, 2015;
originally announced March 2015.
-
Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder
Authors:
Kevin Bandura,
Graeme E. Addison,
Mandana Amiri,
J. Richard Bond,
Duncan Campbell-Wilson,
Liam Connor,
Jean-Francois Cliche,
Greg Davis,
Meiling Deng,
Nolan Denman,
Matt Dobbs,
Mateus Fandino,
Kenneth Gibbs,
Adam Gilbert,
Mark Halpern,
David Hanna,
Adam D. Hincks,
Gary Hinshaw,
Carolin Hofer,
Peter Klages,
Tom L. Landecker,
Kiyoshi Masui,
Juan Mena,
Laura B. Newburgh,
Ue-Li Pen
, et al. (9 additional authors not shown)
Abstract:
A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryo…
▽ More
A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.
△ Less
Submitted 9 June, 2014;
originally announced June 2014.
-
Calibrating CHIME, A New Radio Interferometer to Probe Dark Energy
Authors:
Laura B. Newburgh,
Graeme E. Addison,
Mandana Amiri,
Kevin Bandura,
J. Richard Bond,
Liam Connor,
Jean-François Cliche,
Greg Davis,
Meiling Deng,
Nolan Denman,
Matt Dobbs,
Mateus Fandino,
Heather Fong,
Kenneth Gibbs,
Adam Gilbert,
Elizabeth Griffin,
Mark Halpern,
David Hanna,
Adam D. Hincks,
Gary Hinshaw,
Carolin Höfer,
Peter Klages,
Tom Landecker,
Kiyoshi Masui,
Juan Mena Parra
, et al. (10 additional authors not shown)
Abstract:
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark…
▽ More
The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.
△ Less
Submitted 9 June, 2014;
originally announced June 2014.