Showing 1–50 of 149 results for author: Braun, R

Model-independent searches of new physics in DARWIN with a semi-supervised deep learning pipeline

Authors: J. Aalbers, K. Abe, M. Adrover, S. Ahmed Maouloud, L. Althueser, D. W. P. Amaral, B. Andrieu, E. Angelino, D. Antón Martin, B. Antunovic, E. Aprile, M. Babicz, D. Bajpai, M. Balzer, E. Barberio, L. Baudis, M. Bazyk, N. F. Bell, L. Bellagamba, R. Biondi, Y. Biondi, A. Bismark, C. Boehm, K. Boese, R. Braun , et al. (209 additional authors not shown)

Abstract: We present a novel deep learning pipeline to perform a model-independent, likelihood-free search for anomalous (i.e., non-background) events in the proposed next generation multi-ton scale liquid Xenon-based direct detection experiment, DARWIN. We train an anomaly detector comprising a variational autoencoder and a classifier on extensive, high-dimensional simulated detector response data and cons… ▽ More We present a novel deep learning pipeline to perform a model-independent, likelihood-free search for anomalous (i.e., non-background) events in the proposed next generation multi-ton scale liquid Xenon-based direct detection experiment, DARWIN. We train an anomaly detector comprising a variational autoencoder and a classifier on extensive, high-dimensional simulated detector response data and construct a one-dimensional anomaly score optimised to reject the background only hypothesis in the presence of an excess of non-background-like events. We benchmark the procedure with a sensitivity study that determines its power to reject the background-only hypothesis in the presence of an injected WIMP dark matter signal, outperforming the classical, likelihood-based background rejection test. We show that our neural networks learn relevant energy features of the events from low-level, high-dimensional detector outputs, without the need to compress this data into lower-dimensional observables, thus reducing computational effort and information loss. For the future, our approach lays the foundation for an efficient end-to-end pipeline that eliminates the need for many of the corrections and cuts that are traditionally part of the analysis chain, with the potential of achieving higher accuracy and significant reduction of analysis time. △ Less

Submitted 1 October, 2024; originally announced October 2024.

Comments: 10 Figures, 3 Tables, 23 Pages (incl. references)

Evaporation-driven tear film thinning and breakup in two space dimensions

Authors: Qinying Chen, Tobin A. Driscoll, Richard J. Braun

Abstract: Evaporation profiles have a strong effect on tear film thinning and breakup (TBU), a key factor in dry eye disease (DED). In experiments, TBU is typically seen to occur in patterns that locally can be circular (spot), linear (streak), or intermediate . We investigate a two-dimensional (2D) model of localized TBU using a Fourier spectral collocation method to observe how the evaporation distributio… ▽ More Evaporation profiles have a strong effect on tear film thinning and breakup (TBU), a key factor in dry eye disease (DED). In experiments, TBU is typically seen to occur in patterns that locally can be circular (spot), linear (streak), or intermediate . We investigate a two-dimensional (2D) model of localized TBU using a Fourier spectral collocation method to observe how the evaporation distribution affects the resulting dynamics of tear film thickness and osmolarity, among other variables. We find that the dynamics are not simply an addition of individual 1D solutions of independent TBU events, and we show how the TBU quantities of interest vary continuously from spots to streaks for the shape of the evaporation distribution. We also find a significant speedup by using a proper orthogonal decomposition to reduce the dimension of the numerical system. The speedup will be especially useful for future applications of the model to inverse problems, allowing the clinical observation at scale of quantities that are thought to be important to DED but not directly measurable in vivo within TBU locales. △ Less

Submitted 29 August, 2024; originally announced August 2024.

MSC Class: 76-10; 76Z05

Subsurface Scattering for 3D Gaussian Splatting

Authors: Jan-Niklas Dihlmann, Arjun Majumdar, Andreas Engelhardt, Raphael Braun, Hendrik P. A. Lensch

Abstract: 3D reconstruction and relighting of objects made from scattering materials present a significant challenge due to the complex light transport beneath the surface. 3D Gaussian Splatting introduced high-quality novel view synthesis at real-time speeds. While 3D Gaussians efficiently approximate an object's surface, they fail to capture the volumetric properties of subsurface scattering. We propose a… ▽ More 3D reconstruction and relighting of objects made from scattering materials present a significant challenge due to the complex light transport beneath the surface. 3D Gaussian Splatting introduced high-quality novel view synthesis at real-time speeds. While 3D Gaussians efficiently approximate an object's surface, they fail to capture the volumetric properties of subsurface scattering. We propose a framework for optimizing an object's shape together with the radiance transfer field given multi-view OLAT (one light at a time) data. Our method decomposes the scene into an explicit surface represented as 3D Gaussians, with a spatially varying BRDF, and an implicit volumetric representation of the scattering component. A learned incident light field accounts for shadowing. We optimize all parameters jointly via ray-traced differentiable rendering. Our approach enables material editing, relighting and novel view synthesis at interactive rates. We show successful application on synthetic data and introduce a newly acquired multi-view multi-light dataset of objects in a light-stage setup. Compared to previous work we achieve comparable or better results at a fraction of optimization and rendering time while enabling detailed control over material attributes. Project page https://sss.jdihlmann.com/ △ Less

Submitted 22 August, 2024; originally announced August 2024.

Comments: Project page: https://sss.jdihlmann.com/

Open-Source Conversational AI with SpeechBrain 1.0

Authors: Mirco Ravanelli, Titouan Parcollet, Adel Moumen, Sylvain de Langen, Cem Subakan, Peter Plantinga, Yingzhi Wang, Pooneh Mousavi, Luca Della Libera, Artem Ploujnikov, Francesco Paissan, Davide Borra, Salah Zaiem, Zeyu Zhao, Shucong Zhang, Georgios Karakasidis, Sung-Lin Yeh, Pierre Champion, Aku Rouhe, Rudolf Braun, Florian Mai, Juan Zuluaga-Gomez, Seyed Mahed Mousavi, Andreas Nautsch, Xuechen Liu , et al. (7 additional authors not shown)

Abstract: SpeechBrain is an open-source Conversational AI toolkit based on PyTorch, focused particularly on speech processing tasks such as speech recognition, speech enhancement, speaker recognition, text-to-speech, and much more. It promotes transparency and replicability by releasing both the pre-trained models and the complete "recipes" of code and algorithms required for training them. This paper prese… ▽ More SpeechBrain is an open-source Conversational AI toolkit based on PyTorch, focused particularly on speech processing tasks such as speech recognition, speech enhancement, speaker recognition, text-to-speech, and much more. It promotes transparency and replicability by releasing both the pre-trained models and the complete "recipes" of code and algorithms required for training them. This paper presents SpeechBrain 1.0, a significant milestone in the evolution of the toolkit, which now has over 200 recipes for speech, audio, and language processing tasks, and more than 100 models available on Hugging Face. SpeechBrain 1.0 introduces new technologies to support diverse learning modalities, Large Language Model (LLM) integration, and advanced decoding strategies, along with novel models, tasks, and modalities. It also includes a new benchmark repository, offering researchers a unified platform for evaluating models across diverse tasks. △ Less

Submitted 16 October, 2024; v1 submitted 29 June, 2024; originally announced July 2024.

Comments: Accepted to the Journal of Machine Learning research (JMLR), Machine Learning Open Source Software

High fidelity distribution of triggered polarization-entangled telecom photons via a 36km intra-city fiber network

Authors: Tim Strobel, Stefan Kazmaier, Tobias Bauer, Marlon Schäfer, Ankita Choudhary, Nand Lal Sharma, Raphael Joos, Cornelius Nawrath, Jonas H. Weber, Weijie Nie, Ghata Bhayani, Lukas Wagner, André Bisquerra, Marc Geitz, Ralf-Peter Braun, Caspar Hopfmann, Simone L. Portalupi, Christoph Becher, Peter Michler

Abstract: Fiber-based distribution of triggered, entangled, single-photon pairs is a key requirement for the future development of terrestrial quantum networks. In this context, semiconductor quantum dots (QDs) are promising candidates for deterministic sources of on-demand polarization-entangled photon pairs. So far, the best QD polarization-entangled-pair sources emit in the near-infrared wavelength regim… ▽ More Fiber-based distribution of triggered, entangled, single-photon pairs is a key requirement for the future development of terrestrial quantum networks. In this context, semiconductor quantum dots (QDs) are promising candidates for deterministic sources of on-demand polarization-entangled photon pairs. So far, the best QD polarization-entangled-pair sources emit in the near-infrared wavelength regime, where the transmission distance in deployed fibers is limited. Here, to be compatible with existing fiber network infrastructures, bi-directional polarization-conserving quantum frequency conversion (QFC) is employed to convert the QD emission from \unit[780]{nm} to telecom wavelengths. We show the preservation of polarization entanglement after QFC (fidelity to Bell state $F_{φ^+, conv}=0.972\pm0.003$) of the biexciton transition. As a step towards real-world applicability, high entanglement fidelities ($F_{φ^+, loop}=0.945\pm0.005$) after the propagation of one photon of the entangled pair along a \unit[35.8]{km} field installed standard single mode fiber link are reported. Furthermore, we successfully demonstrate a second polarization-conversing QFC step back to \unit[780]{nm} preserving entanglement ($F_{φ^+, back}=0.903\pm0.005$). This further prepares the way for interfacing quantum light to various quantum memories. △ Less

Submitted 27 May, 2024; v1 submitted 23 May, 2024; originally announced May 2024.

On Modeling Tear Breakup Dynamics with a Nematic Lipid Layer

Authors: M. J. Taranchuk, R. J. Braun

Abstract: One of the main roles of the lipid layer (LL) of the tear film (TF) is to help prevent evaporation of the aqueous layer (AL). The LL thickness, composition, and structure all contribute to its barrier function. It is believed that the lipid layer is primarily nonpolar with a layer of polar lipids at the LL/AL interface. There is evidence that the nonpolar region of the LL may have liquid crystalli… ▽ More One of the main roles of the lipid layer (LL) of the tear film (TF) is to help prevent evaporation of the aqueous layer (AL). The LL thickness, composition, and structure all contribute to its barrier function. It is believed that the lipid layer is primarily nonpolar with a layer of polar lipids at the LL/AL interface. There is evidence that the nonpolar region of the LL may have liquid crystalline characteristics. We investigate the structure and function of the LL via a model of the tear film with two layers, using extensional flow of a nematic liquid crystal for the LL and shear-dominated flow of a Newtonian AL. Evaporation is taken into account, and is affected by the LL thickness, internal arrangement of its rod-like molecules, and external conditions. We conduct a detailed parameter study with a focus on the evaporative resistance parameter, the Marangoni number, and primary liquid crystal parameters including the Leslie viscosities and director angle. This new model responds similarly to previous Newtonian models in some respects; however, incorporating internal structure via the orientation of the liquid crystal molecules affects both evaporation and flow. As a result, we see new effects on TF dynamics and breakup. △ Less

Submitted 25 April, 2024; originally announced April 2024.

On the Effect of Liquid Crystal Orientation in the Lipid Layer on Tear Film Thinning and Breakup

Authors: M. J. Taranchuk, R. J. Braun

Abstract: The human tear film (TF) is thin multilayer fluid film that is critical for clear vision and ocular surface health. Its dynamics are strongly affected by a floating lipid layer and, in health, that layer slows evaporation and helps create a more uniform tear film over the ocular surface. The tear film lipid layer (LL) may have liquid crystalline characteristics and plays important roles in the hea… ▽ More The human tear film (TF) is thin multilayer fluid film that is critical for clear vision and ocular surface health. Its dynamics are strongly affected by a floating lipid layer and, in health, that layer slows evaporation and helps create a more uniform tear film over the ocular surface. The tear film lipid layer (LL) may have liquid crystalline characteristics and plays important roles in the health of the tear film. Previous models have treated the lipid layer as a Newtonian fluid in extensional flow. We extend previous models to include extensional flow of a thin nematic liquid crystal atop a Newtonian aqueous layer with insoluble surfactant between them. We derive the resulting system of nonlinear partial differential equations for thickness of the LL and aqueous layers, surfactant transport and velocity in the LL. Evaporation is taken into account, and is affected by the LL thickness, internal arrangement of its rod-like molecules, and external conditions. Despite the complexity, this system still represents a significant reduction of the full system. We solve the system numerically via collocation with finite difference discretization in space together with implicit time stepping. We analyze solutions for different internal LL structures and show significant effect of the orientation. Orienting the molecules close to the normal direction to the TF surface results in slower evaporation, and other orientations have an effect on flow, showing that this type of model has promise for predicting TF dynamics. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: Copyright (2024) M.J. Taranchuk and R.J. Braun. This article is distributed under a Creative Commons Attribution (CC BY) License

Linking QKD testbeds across Europe

Authors: Max Brauer, Rafael J. Vicente, Jaime S. Buruaga, Ruben B. Mendez, Ralf-Peter Braun, Marc Geitz, Piotr Rydlichkowski, Hans H. Brunner, Fred Fung, Momtchil Peev, Antonio Pastor, Diego Lopez, Vicente Martin, Juan P. Brito

Abstract: Quantum-key-distribution (QKD) networks are gaining importance and it has become necessary to analyze the most appropriate methods for their long-distance interconnection. In this paper, four different methods of interconnecting remote QKD networks are proposed. The methods are used to link three different QKD testbeds in Europe, located in Berlin, Madrid, and Poznan. Although long-distance QKD li… ▽ More Quantum-key-distribution (QKD) networks are gaining importance and it has become necessary to analyze the most appropriate methods for their long-distance interconnection. In this paper, four different methods of interconnecting remote QKD networks are proposed. The methods are used to link three different QKD testbeds in Europe, located in Berlin, Madrid, and Poznan. Although long-distance QKD links are only emulated, the used methods can serve as a blueprint for a secure interconnection of distant QKD networks in the future. Specifically, the presented approaches combine, in a transparent way, different fiber and satellite physical media, as well as common standards of key-delivery interfaces. The testbed interconnections are designed to increase the security by utilizing multipath techniques and multiple hybridizations of QKD and post quantum cryptography (PQC) algorithms. △ Less

Submitted 10 January, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Automatized Self-Supervised Learning for Skin Lesion Screening

Authors: Vullnet Useini, Stephanie Tanadini-Lang, Quentin Lohmeyer, Mirko Meboldt, Nicolaus Andratschke, Ralph P. Braun, Javier Barranco García

Abstract: The incidence rates of melanoma, the deadliest form of skin cancer, have been increasing steadily worldwide, presenting a significant challenge to dermatologists. Early detection of melanoma is crucial for improving patient survival rates, but identifying suspicious lesions through ugly duckling (UD) screening, the current method used for skin cancer screening, can be challenging and often require… ▽ More The incidence rates of melanoma, the deadliest form of skin cancer, have been increasing steadily worldwide, presenting a significant challenge to dermatologists. Early detection of melanoma is crucial for improving patient survival rates, but identifying suspicious lesions through ugly duckling (UD) screening, the current method used for skin cancer screening, can be challenging and often requires expertise in pigmented lesions. To address these challenges and improve patient outcomes, an artificial intelligence (AI) decision support tool was developed to assist dermatologists in identifying UD from wide-field patient images. The tool uses a state-of-the-art object detection algorithm to identify and extract all skin lesions from patient images, which are then sorted by suspiciousness using a self-supervised AI algorithm. A clinical validation study was conducted to evaluate the tool's performance, which demonstrated an average sensitivity of 93% for the top-10 AI-identified UDs on skin lesions selected by the majority of experts in pigmented skin lesions. The study also found that dermatologists confidence increased, and the average majority agreement with the top-10 AI-identified UDs improved to 100% when assisted by AI. The development of this AI decision support tool aims to address the shortage of specialists, enable at-risk patients to receive faster consultations and understand the impact of AI-assisted screening. The tool's automation can assist dermatologists in identifying suspicious lesions and provide a more objective assessment, reducing subjectivity in the screening process. The future steps for this project include expanding the dataset to include histologically confirmed melanoma cases and increasing the number of participants for clinical validation to strengthen the tool's reliability and adapt it for real-world consultation. △ Less

Submitted 11 November, 2023; originally announced November 2023.

Comments: 11 pages, 4 figures

Burnt area extraction from high-resolution satellite images based on anomaly detection

Authors: Oscar David Rafael Narvaez Luces, Minh-Tan Pham, Quentin Poterek, Rémi Braun

Abstract: Wildfire detection using satellite images is a widely studied task in remote sensing with many applications to fire delineation and mapping. Recently, deep learning methods have become a scalable solution to automate this task, especially in the field of unsupervised learning where no training data is available. This is particularly important in the context of emergency risk monitoring where fast… ▽ More Wildfire detection using satellite images is a widely studied task in remote sensing with many applications to fire delineation and mapping. Recently, deep learning methods have become a scalable solution to automate this task, especially in the field of unsupervised learning where no training data is available. This is particularly important in the context of emergency risk monitoring where fast and effective detection is needed, generally based on high-resolution satellite data. Among various approaches, Anomaly Detection (AD) appears to be highly potential thanks to its broad applications in computer vision, medical imaging, as well as remote sensing. In this work, we build upon the framework of Vector Quantized Variational Autoencoder (VQ-VAE), a popular reconstruction-based AD method with discrete latent spaces, to perform unsupervised burnt area extraction. We integrate VQ-VAE into an end-to-end framework with an intensive post-processing step using dedicated vegetation, water and brightness indexes. Our experiments conducted on high-resolution SPOT-6/7 images provide promising results of the proposed technique, showing its high potential in future research on unsupervised burnt area extraction. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 10 pages, accepted to the MACLEAN workshop of ECML/PKDD 2023

CANDID: Correspondence AligNment for Deep-burst Image Denoising

Authors: Arijit Mallick, Raphael Braun, Hendrik PA Lensch

Abstract: With the advent of mobile phone photography and point-and-shoot cameras, deep-burst imaging is widely used for a number of photographic effects such as depth of field, super-resolution, motion deblurring, and image denoising. In this work, we propose to solve the problem of deep-burst image denoising by including an optical flow-based correspondence estimation module which aligns all the input bur… ▽ More With the advent of mobile phone photography and point-and-shoot cameras, deep-burst imaging is widely used for a number of photographic effects such as depth of field, super-resolution, motion deblurring, and image denoising. In this work, we propose to solve the problem of deep-burst image denoising by including an optical flow-based correspondence estimation module which aligns all the input burst images with respect to a reference frame. In order to deal with varying noise levels the individual burst images are pre-filtered with different settings. Exploiting the established correspondences one network block predicts a pixel-wise spatially-varying filter kernel to smooth each image in the original and prefiltered bursts before fusing all images to generate the final denoised output. The resulting pipeline achieves state-of-the-art results by combining all available information provided by the burst. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Paper accepted and presented as a poster on 20th Conference on Robots and Vision 2023

Mildly Relativistic Motion in the Radio Quiet Quasar PG 1351+640

Authors: Ailing Wang, Tao An, Shaoguang Guo, Luis C. Ho, Willem A. Baan, Robert Braun, Sina Chen, Xiaopeng Cheng, Philippa Hartley, Jun Yang, Yingkang Zhang

Abstract: Measuring the proper motion of the emission component in radio-quiet quasars (RQQs) could help to distinguish between the origins of the radio emission and to understand whether the jet production mechanism is the same in radio-loud quasars (RLQs) and RQQs. PG 1351+640 is one of the few RQQs suitable for proper motion studies: it has two compact components on milli-arcsecond scales, a flat-spectru… ▽ More Measuring the proper motion of the emission component in radio-quiet quasars (RQQs) could help to distinguish between the origins of the radio emission and to understand whether the jet production mechanism is the same in radio-loud quasars (RLQs) and RQQs. PG 1351+640 is one of the few RQQs suitable for proper motion studies: it has two compact components on milli-arcsecond scales, a flat-spectrum core and a steep-spectrum jet; both components are >2 mJy at 5 GHz and are well suited for Very Long Baseline Array (VLBA) observations. We compare recent VLBA observations with that made seventeen years ago and find no significant change in the core-jet separation between 2005 and 2015 (a proper motion of 0.003 mas yr-1). However, the core-jet separation increased significantly between 2015 and 2022, inferring a jet proper motion velocity of 0.063 mas yr-1, which corresponds to an apparent transverse velocity of 0.37c. The result suggests that the jet of the RQQ PG 1351+640 is mildly relativistic and oriented at a relatively small viewing angle. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: The article has been published by Oxford University Press: https://academic.oup.com/mnrasl/advance-article-abstract/doi/10.1093/mnrasl/slad051/7146829?utm_source=advanceaccess&utm_campaign=mnrasl&utm_medium=email

Modeling Transient Changes in Circadian Rhythms

Authors: Ziyu Zhao, Dae-Sung Hwangbo, Sumit Saurabh, Clark Rosensweig, Ravi Allada, William L. Kath, Rosemary Braun

Abstract: The circadian clock can adapt itself to external cues, but the molecular mechanisms and regulatory networks governing circadian oscillations' transient adjustments are still largely unknown. Here we consider the specific case of circadian oscillations transiently responding to a temperature change. Using a framework motivated by Floquet theory, we model the mRNA expression level of the fat body fr… ▽ More The circadian clock can adapt itself to external cues, but the molecular mechanisms and regulatory networks governing circadian oscillations' transient adjustments are still largely unknown. Here we consider the specific case of circadian oscillations transiently responding to a temperature change. Using a framework motivated by Floquet theory, we model the mRNA expression level of the fat body from Drosophila melanogaster following a step change from 25C to 18C. Using the method we infer the adaptation rates of individual genes as they adapt to the new temperature. To deal with heteroskedastic noise and outliers present in the expression data we employ quantile regression and wild bootstrap for significance testing. Model selection with finite-sample corrected Akaike Information Criterion (AICc) is performed additionally for robust inference. We identify several genes with fast transition rates as potential sources of temperature-mediated responses in the circadian system of fruit flies, and the constructed network suggests that the proteasome may play important roles in governing these responses. △ Less

Submitted 14 April, 2023; originally announced April 2023.

Extensional Flow of a Free Film of Nematic Liquid Crystal with Moderate Elasticity

Authors: M. J. Taranchuk, L. J. Cummings, T. A. Driscoll, R. J. Braun

Abstract: Motivated by problems arising in tear film dynamics, we present a model for the extensional flow of thin sheets of nematic liquid crystal. The rod-like molecules of these substances impart an elastic contribution to its response. We rescale a weakly elastic model due to Cummings et al. [European Journal of Applied Mathematics 25 (2014): 397-423] to describe a case of moderate elasticity. The resul… ▽ More Motivated by problems arising in tear film dynamics, we present a model for the extensional flow of thin sheets of nematic liquid crystal. The rod-like molecules of these substances impart an elastic contribution to its response. We rescale a weakly elastic model due to Cummings et al. [European Journal of Applied Mathematics 25 (2014): 397-423] to describe a case of moderate elasticity. The resulting system of two nonlinear partial differential equations for sheet thickness and axial velocity is nonlinear and fourth order in space, but still represents a significant reduction of the full system. We analyze solutions arising from several different boundary conditions, motivated by the underlying application, with particular focus on dynamics and underlying mechanisms under stretching. We solve the system numerically, via collocation with either finite difference or Chebyshev spectral discretization in space, together with implicit time stepping. At early times, depending on the initial film shape, pressure either aids or opposes extensional flow, which changes the shape of the sheet and may result in the loss of a minimum or maximum at the moving end. We contrast this finding with the cases of weak elasticity and Newtonian flow, where the sheet retains all extrema from the initial condition throughout time. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: Copyright (2023) M.J. Taranchuk, L.J. Cummings, T.A. Driscoll, and R.J. Braun. This article is distributed under a Creative Commons Attribution (CC BY) License. The following article has been submitted to Physics of Fluids. After it is published, it will be found at https://aip.scitation.org/journal/phf

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

A minimal model of peripheral clocks reveals differential circadian re-entrainment in aging

Authors: Yitong Huang, Yuanzhao Zhang, Rosemary Braun

Abstract: The mammalian circadian system comprises a network of cell-autonomous oscillators, spanning from the central clock in the brain to peripheral clocks in other organs. These clocks are tightly coordinated to orchestrate rhythmic physiological and behavioral functions. Dysregulation of these rhythms is a hallmark of aging, yet it remains unclear how age-related changes lead to more easily disrupted c… ▽ More The mammalian circadian system comprises a network of cell-autonomous oscillators, spanning from the central clock in the brain to peripheral clocks in other organs. These clocks are tightly coordinated to orchestrate rhythmic physiological and behavioral functions. Dysregulation of these rhythms is a hallmark of aging, yet it remains unclear how age-related changes lead to more easily disrupted circadian rhythms. Using a two-population model of coupled oscillators that integrates the central clock and the peripheral clocks, we derive simple mean-field equations that can capture many aspects of the rich behavior found in the mammalian circadian system. We focus on three age-associated effects which have been posited to contribute to circadian misalignment: attenuated input from the sympathetic pathway, reduced responsiveness to light, and a decline in the expression of neurotransmitters. We find that the first two factors can significantly impede re-entrainment of the clocks following a perturbation, while a weaker coupling within the central clock does not affect the recovery rate. Moreover, using our minimal model, we demonstrate the potential of using the feed-fast cycle as an effective intervention to accelerate circadian re-entrainment. These results highlight the importance of peripheral clocks in regulating the circadian rhythm and provide fresh insights into the complex interplay between aging and the resilience of the circadian system. △ Less

Submitted 27 December, 2022; originally announced December 2022.

Fitting ODE models of tear film breakup

Authors: Tobin A. Driscoll, Richard J. Braun, Rayanne A. Luke, Dominick Sinopoli, Aashish Phatak, Julianna Dorsch, Carolyn G. Begley, Deborah Awisi-Gyau

Abstract: The contribution of different physical effects to tear breakup (TBU) in subjects with no self-reported history of dry eye are quantified. An automated system using a convolutional neural network is deployed on fluorescence (FL) imaging videos to identify multiple likely TBU instances in each trial. Once identified, extracted FL intensity data was fit by mathematical models that included tangential… ▽ More The contribution of different physical effects to tear breakup (TBU) in subjects with no self-reported history of dry eye are quantified. An automated system using a convolutional neural network is deployed on fluorescence (FL) imaging videos to identify multiple likely TBU instances in each trial. Once identified, extracted FL intensity data was fit by mathematical models that included tangential flow along the eye, evaporation, osmosis and FL intensity of emission from the tear film. The mathematical models consisted of systems of ordinary differential equations for the aqueous layer thickness, osmolarity, and the FL concentration. Optimizing the fit of the models to the FL intensity data determined the mechanism(s) driving each instance of TBU and produced an estimate of the osmolarity within TBU. Fits were produced for 467 instances of potential TBU from 15 non-DED subjects. The results showed a distribution of causes of TBU in these healthy subjects, as reflected by estimated flow and evaporation rates, which appear to agree well with previously published data. Final osmolarity depended strongly on the TBU mechanism, generally increasing with evaporation rate but complicated by the dependence on flow. The results suggest that it might be possible to classify individual subjects and provide a baseline for comparison and potential classification of dry eye disease subjects. △ Less

Submitted 21 February, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

A Robust, Performance-Portable Discontinuous Galerkin Method for Relativistic Hydrodynamics

Authors: Forrest W. Glines, Kristian R. C. Beckwith, Joshua R. Braun, Eric C. Cyr, Curtis C. Ober, Matthew Bettencourt, Keith L. Cartwright, Sidafa Conde, Sean T. Miller, Nicholas Roberds, Nathan V. Roberts, Matthew S. Swan, Roger Pawlowski

Abstract: In this work, we present a discontinuous-Galerkin method for evolving relativistic hydrodynamics. We include an exploration of analytical and iterative methods to recover the primitive variables from the conserved variables for the ideal equation of state and the Taub-Matthews approximation to the Synge equation of state. We also present a new operator for enforcing a physically permissible conser… ▽ More In this work, we present a discontinuous-Galerkin method for evolving relativistic hydrodynamics. We include an exploration of analytical and iterative methods to recover the primitive variables from the conserved variables for the ideal equation of state and the Taub-Matthews approximation to the Synge equation of state. We also present a new operator for enforcing a physically permissible conserved state at all basis points within an element while preserving the volume average of the conserved state. We implement this method using the Kokkos performance-portability library to enable running at performance on both CPUs and GPUs. We use this method to explore the relativistic Kelvin- Helmholtz instability compared to a finite volume method. Last, we explore the performance of our implementation on CPUs and GPUs. △ Less

Submitted 29 April, 2022; originally announced May 2022.

Demonstration of latency-aware 5G network slicing on optical metro networks

Authors: B. Shariati, L. Velasco, J. -J. Pedreño-Manresa, A. Dochhan, R. Casellas, A. Muqaddas, O. González de Dios, L. Luque Canto, B. Lent, J. E. López de Vergara, S. López-Buedo, F. Moreno, P. Pavón, M. Ruiz, S. K. Patri, A. Giorgetti, F. Cugini, A. Sgambelluri, R. Nejabati, D. Simeonidou, R. -P. Braun, A. Autenrieth, J. -P. Elbers, J. K. Fischer, R. Freund

Abstract: The H2020 METRO-HAUL European project has architected a latency-aware, cost-effective, agile, and programmable optical metro network. This includes the design of semidisaggregated metro nodes with compute and storage capabilities, which interface effectively with both 5G access and multi-Tbit/s elastic optical networks in the core. In this paper, we report the automated deployment of 5G services,… ▽ More The H2020 METRO-HAUL European project has architected a latency-aware, cost-effective, agile, and programmable optical metro network. This includes the design of semidisaggregated metro nodes with compute and storage capabilities, which interface effectively with both 5G access and multi-Tbit/s elastic optical networks in the core. In this paper, we report the automated deployment of 5G services, in particular, a public safety video surveillance use case employing low-latency object detection and tracking using on-camera and on-the-edge analytics. The demonstration features flexible deployment of network slice instances, implemented in terms of European Telecommunications Standards Institute (ETSI) network function virtualization network services. We summarize the key findings in a detailed analysis of end-to-end quality of service, service setup time, and soft-failure detection time. The results show that the round-trip time over an 80 km link is under 800s and the service deployment time is under 180s. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 10 pages

Journal ref: Journal of Optical Communication and Networking, 14, A81-A90 (2022)

Neural-PIL: Neural Pre-Integrated Lighting for Reflectance Decomposition

Authors: Mark Boss, Varun Jampani, Raphael Braun, Ce Liu, Jonathan T. Barron, Hendrik P. A. Lensch

Abstract: Decomposing a scene into its shape, reflectance and illumination is a fundamental problem in computer vision and graphics. Neural approaches such as NeRF have achieved remarkable success in view synthesis, but do not explicitly perform decomposition and instead operate exclusively on radiance (the product of reflectance and illumination). Extensions to NeRF, such as NeRD, can perform decomposition… ▽ More Decomposing a scene into its shape, reflectance and illumination is a fundamental problem in computer vision and graphics. Neural approaches such as NeRF have achieved remarkable success in view synthesis, but do not explicitly perform decomposition and instead operate exclusively on radiance (the product of reflectance and illumination). Extensions to NeRF, such as NeRD, can perform decomposition but struggle to accurately recover detailed illumination, thereby significantly limiting realism. We propose a novel reflectance decomposition network that can estimate shape, BRDF, and per-image illumination given a set of object images captured under varying illumination. Our key technique is a novel illumination integration network called Neural-PIL that replaces a costly illumination integral operation in the rendering with a simple network query. In addition, we also learn deep low-dimensional priors on BRDF and illumination representations using novel smooth manifold auto-encoders. Our decompositions can result in considerably better BRDF and light estimates enabling more accurate novel view-synthesis and relighting compared to prior art. Project page: https://markboss.me/publication/2021-neural-pil/ △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: Project page: https://markboss.me/publication/2021-neural-pil/ Video: https://youtu.be/AsdAR5u3vQ8 - Accepted at NeurIPS 2021

Stroke recovery phenotyping through network trajectory approaches and graph neural networks

Authors: Sanjukta Krishnagopal, Keith Lohse, Robynne Braun

Abstract: Stroke is a leading cause of neurological injury characterized by impairments in multiple neurological domains including cognition, language, sensory and motor functions. Clinical recovery in these domains is tracked using a wide range of measures that may be continuous, ordinal, interval or categorical in nature, which presents challenges for standard multivariate regression approaches. This has… ▽ More Stroke is a leading cause of neurological injury characterized by impairments in multiple neurological domains including cognition, language, sensory and motor functions. Clinical recovery in these domains is tracked using a wide range of measures that may be continuous, ordinal, interval or categorical in nature, which presents challenges for standard multivariate regression approaches. This has hindered stroke researchers' ability to achieve an integrated picture of the complex time-evolving interactions amongst symptoms. Here we use tools from network science and machine learning that are particularly well-suited to extracting underlying patterns in such data, and may assist in prediction of recovery patterns. To demonstrate the utility of this approach, we analyzed data from the NINDS tPA trial using the Trajectory Profile Clustering (TPC) method to identify distinct stroke recovery patterns for 11 different neurological domains at 5 discrete time points. Our analysis identified 3 distinct stroke trajectory profiles that align with clinically relevant stroke syndromes, characterized both by distinct clusters of symptoms, as well as differing degrees of symptom severity. We then validated our approach using graph neural networks to determine how well our model performed predictively for stratifying patients into these trajectory profiles at early vs. later time points post-stroke. We demonstrate that trajectory profile clustering is an effective method for identifying clinically relevant recovery subtypes in multidimensional longitudinal datasets, and for early prediction of symptom progression subtypes in individual patients. This paper is the first work introducing network trajectory approaches for stroke recovery phenotyping, and is aimed at enhancing the translation of such novel computational approaches for practical clinical application. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: 20 pages, 5 figures

Apertif, Phased Array Feeds for the Westerbork Synthesis Radio Telescope

Authors: W. A. van Cappellen, T. A. Oosterloo, M. A. W. Verheijen, E. A. K. Adams, B. Adebahr, R. Braun, K. M. Hess, H. Holties, J. M. van der Hulst, B. Hut, E. Kooistra, J. van Leeuwen, G. M. Loose, R. Morganti, V. A. Moss, E. Orrú, M. Ruiter, A. P. Schoenmakers, N. J. Vermaas, S. J. Wijnholds, A. S. van Amesfoort, M. J. Arts, J. J. Attema, L. Bakker, C. G. Bassa , et al. (65 additional authors not shown)

Abstract: We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the s… ▽ More We describe the APERture Tile In Focus (Apertif) system, a phased array feed (PAF) upgrade of the Westerbork Synthesis Radio Telescope which has transformed this telescope into a high-sensitivity, wide field-of-view L-band imaging and transient survey instrument. Using novel PAF technology, up to 40 partially overlapping beams can be formed on the sky simultaneously, significantly increasing the survey speed of the telescope. With this upgraded instrument, an imaging survey covering an area of 2300 deg2 is being performed which will deliver both continuum and spectral line data sets, of which the first data has been publicly released. In addition, a time domain transient and pulsar survey covering 15,000 deg2 is in progress. An overview of the Apertif science drivers, hardware and software of the upgraded telescope is presented, along with its key performance characteristics. △ Less

Submitted 30 September, 2021; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 29 pages, 42 figures, accepted for publication by A&A

Journal ref: A&A 658, A146 (2022)

Improving callsign recognition with air-surveillance data in air-traffic communication

Authors: Iuliia Nigmatulina, Rudolf Braun, Juan Zuluaga-Gomez, Petr Motlicek

Abstract: Automatic Speech Recognition (ASR) can be used as the assistance of speech communication between pilots and air-traffic controllers. Its application can significantly reduce the complexity of the task and increase the reliability of transmitted information. Evidently, high accuracy predictions are needed to minimize the risk of errors. Especially, high accuracy is required in recognition of key in… ▽ More Automatic Speech Recognition (ASR) can be used as the assistance of speech communication between pilots and air-traffic controllers. Its application can significantly reduce the complexity of the task and increase the reliability of transmitted information. Evidently, high accuracy predictions are needed to minimize the risk of errors. Especially, high accuracy is required in recognition of key information, such as commands and callsigns, used to navigate pilots. Our results prove that the surveillance data containing callsigns can help to considerably improve the recognition of a callsign in an utterance when the weights of probable callsign n-grams are reduced per utterance. In this paper, we investigate two approaches: (1) G-boosting, when callsigns weights are adjusted at language model level (G) and followed by the dynamic decoder with an on-the-fly composition, and (2) lattice rescoring when callsign information is introduced on top of lattices generated using a conventional decoder. Boosting callsign n-grams with the combination of two methods allowed us to gain 28.4% of absolute improvement in callsign recognition accuracy and up to 74.2% of relative improvement in WER of callsign recognition. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: Submitted to Interspeech 2021

1.71 Tb/s Single-Channel and 56.51 Tb/s DWDM Transmission over 96.5 km Field-Deployed SSMF

Authors: Fabio Pittala, Ralf-Peter Braun, Georg Boecherer, Patrick Schulte, Maximilian Schaedler, Stefano Bettelli, Stefano Calabro, Maxim Kuschnerov, Andreas Gladisch, Fritz-Joachim Westphal, Changsong Xie, Rongfu Chen, Qibing Wang, Bofang Zheng

Abstract: We report an industry leading optical dense wavelength division multiplexing (DWDM) field trial with line rates per channel exceeding 1.66 Tb/s using 130 GBaud dual-polarization probabilistic constellation shaping 256-ary quadrature amplitude modulation (DP-PCS256QAM) in a high capacity data center interconnect (DCI) scenario. This research trial was performed on 96.5 km of field-deployed standard… ▽ More We report an industry leading optical dense wavelength division multiplexing (DWDM) field trial with line rates per channel exceeding 1.66 Tb/s using 130 GBaud dual-polarization probabilistic constellation shaping 256-ary quadrature amplitude modulation (DP-PCS256QAM) in a high capacity data center interconnect (DCI) scenario. This research trial was performed on 96.5 km of field-deployed standard single mode G.652 fiber infrastructure of Deutsche Telekom in Germany employing Erbium-doped fiber amplifier (EDFA)-only amplification. A total of 34 channels were transmitted with 150 GHz spacing for a total fiber capacity of 56.51 Tb/s and a spectral efficiency higher than 11bit/s/Hz. In the single-channel transmission scenario 1.71 Tb/s was achieved over the same link. In addition, we successfully demonstrate record net bitrates of 1.88 Tb/s in back-to-back (B2B) using 130 GBaud DP-PCS400QAM. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: This work has been submitted to the IEEE Photonics Technology Letters (PTL) for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

A Comparison of Methods for OOV-word Recognition on a New Public Dataset

Authors: Rudolf A. Braun, Srikanth Madikeri, Petr Motlicek

Abstract: A common problem for automatic speech recognition systems is how to recognize words that they did not see during training. Currently there is no established method of evaluating different techniques for tackling this problem. We propose using the CommonVoice dataset to create test sets for multiple languages which have a high out-of-vocabulary (OOV) ratio relative to a training set and release a n… ▽ More A common problem for automatic speech recognition systems is how to recognize words that they did not see during training. Currently there is no established method of evaluating different techniques for tackling this problem. We propose using the CommonVoice dataset to create test sets for multiple languages which have a high out-of-vocabulary (OOV) ratio relative to a training set and release a new tool for calculating relevant performance metrics. We then evaluate, within the context of a hybrid ASR system, how much better subword models are at recognizing OOVs, and how much benefit one can get from incorporating OOV-word information into an existing system by modifying WFSTs. Additionally, we propose a new method for modifying a subword-based language model so as to better recognize OOV-words. We showcase very large improvements in OOV-word recognition and make both the data and code available. △ Less

Submitted 16 July, 2021; originally announced July 2021.

All-Indoor Optical Customer Premises Equipment for Fixed Wireless Access

Authors: Dominic Schulz, Julian Hohmann, Peter Hellwig, Christoph Kottke, Ronald Freund, Volker Jungnickel, Ralf-Peter Braun, Frank Geilhardt

Abstract: We demonstrate an LED-based optical wireless link for fixed wireless access applications, at data rates of 1.5 Gbit/s over 50 m. Transmission between indoor equipment and outdoor access point is possible through metal-coated insulation glass. We demonstrate an LED-based optical wireless link for fixed wireless access applications, at data rates of 1.5 Gbit/s over 50 m. Transmission between indoor equipment and outdoor access point is possible through metal-coated insulation glass. △ Less

Submitted 26 May, 2021; originally announced July 2021.

Comments: 3 pages, 2 figures, Optical Fiber Communication Conference 2021

ACM Class: B.4.1

fasano.franceschini.test: An Implementation of a Multidimensional KS Test in R

Authors: Connor Puritz, Elan Ness-Cohn, Rosemary Braun

Abstract: The Kolmogorov-Smirnov (KS) test is a nonparametric statistical test used to test for differences between univariate probability distributions. The versatility of the KS test has made it a cornerstone of statistical analysis across many scientific disciplines. However, the test proposed by Kolmogorov and Smirnov does not easily extend to multidimensional distributions. Here we present the fasano.f… ▽ More The Kolmogorov-Smirnov (KS) test is a nonparametric statistical test used to test for differences between univariate probability distributions. The versatility of the KS test has made it a cornerstone of statistical analysis across many scientific disciplines. However, the test proposed by Kolmogorov and Smirnov does not easily extend to multidimensional distributions. Here we present the fasano.franceschini.test package, an R implementation of a multidimensional two-sample KS test described by Fasano and Franceschini (1987). The fasano.franceschini.test package provides a test that is computationally efficient, applicable to data of any dimension and type (continuous, discrete, or mixed), and that performs competitively with similar R packages. △ Less

Submitted 18 November, 2022; v1 submitted 19 June, 2021; originally announced June 2021.

Comments: 12 pages, 5 figures

Mechanistic determination of tear film thinning via fitting simplified models to tear breakup

Authors: Rayanne A. Luke, Richard J. Braun, Carolyn G. Begley

Abstract: Purpose: To determine whether evaporation, tangential flow, or a combination of the two cause tear film breakup in a variety of instances; to estimate related breakup parameters that cannot be measured in breakup during subject trials; and to validate our procedure against previous work. Methods: Five ordinary differential equation models for tear film thinning were designed that model evaporation… ▽ More Purpose: To determine whether evaporation, tangential flow, or a combination of the two cause tear film breakup in a variety of instances; to estimate related breakup parameters that cannot be measured in breakup during subject trials; and to validate our procedure against previous work. Methods: Five ordinary differential equation models for tear film thinning were designed that model evaporation, osmosis, and various types of flow. Eight tear film breakup instances of five healthy subjects that were identified in fluorescence images in previous work were fit with these five models. The fitting procedure used a nonlinear least squares optimization that minimized the difference of the computed theoretical fluorescent intensity from the models and the experimental fluorescent intensity from the images. The optimization was conducted over the evaporation rate and up to three flow rate parameters. The smallest norm of the difference was determined to correspond to the model that best explained the tear film dynamics. Results: All of the breakup instances were best fit by models with time-dependent flow. Our optimal parameter values and thinning rate and fluid flow profiles compare well with previous partial differential equation model results in most instances. Conclusion: Our fitting procedure suggests that the combination of the Marangoni effect and evaporation cause most of the breakup instances. Comparison with results from previous work suggests that the simplified models can capture the essential tear film dynamics in most cases, thereby validating this procedure as one that could be used on many other instances. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 28 pages, 11 figures, 6 tables

MSC Class: 92-10

Journal ref: Modeling and Artificial Intelligence in Ophthalmology (2021)

NeRD: Neural Reflectance Decomposition from Image Collections

Authors: Mark Boss, Raphael Braun, Varun Jampani, Jonathan T. Barron, Ce Liu, Hendrik P. A. Lensch

Abstract: Decomposing a scene into its shape, reflectance, and illumination is a challenging but important problem in computer vision and graphics. This problem is inherently more challenging when the illumination is not a single light source under laboratory conditions but is instead an unconstrained environmental illumination. Though recent work has shown that implicit representations can be used to model… ▽ More Decomposing a scene into its shape, reflectance, and illumination is a challenging but important problem in computer vision and graphics. This problem is inherently more challenging when the illumination is not a single light source under laboratory conditions but is instead an unconstrained environmental illumination. Though recent work has shown that implicit representations can be used to model the radiance field of an object, most of these techniques only enable view synthesis and not relighting. Additionally, evaluating these radiance fields is resource and time-intensive. We propose a neural reflectance decomposition (NeRD) technique that uses physically-based rendering to decompose the scene into spatially varying BRDF material properties. In contrast to existing techniques, our input images can be captured under different illumination conditions. In addition, we also propose techniques to convert the learned reflectance volume into a relightable textured mesh enabling fast real-time rendering with novel illuminations. We demonstrate the potential of the proposed approach with experiments on both synthetic and real datasets, where we are able to obtain high-quality relightable 3D assets from image collections. The datasets and code is available on the project page: https://markboss.me/publication/2021-nerd/ △ Less

Submitted 26 August, 2021; v1 submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted at ICCV 2021

Parameter Estimation for Mixed-Mechanism Tear Film Thinning

Authors: Rayanne A. Luke, Richard J. Braun, Tobin A. Driscoll, Deborah Awisi-Gyau, Carolyn G. Begley

Abstract: Etiologies of tear breakup include evaporation-driven, divergent flow-driven, and a combination of these two. A mathematical model incorporating evaporation and lipid-driven tangential flow is fit to fluorescence imaging data. The lipid-driven motion is hypothesized to be caused by localized excess lipid, or "globs." Tear breakup quantities such as evaporation rates and tangential flow rates canno… ▽ More Etiologies of tear breakup include evaporation-driven, divergent flow-driven, and a combination of these two. A mathematical model incorporating evaporation and lipid-driven tangential flow is fit to fluorescence imaging data. The lipid-driven motion is hypothesized to be caused by localized excess lipid, or "globs." Tear breakup quantities such as evaporation rates and tangential flow rates cannot currently be directly measured during breakup. We determine such variables by fitting mathematical models for tear breakup and the computed fluorescent intensity to experimental intensity data gathered in vivo. Parameter estimation is conducted via least squares minimization of the difference between experimental data and computed answers using either the trust-region-reflective or Levenberg-Marquardt algorithm. Best-fit determination of tear breakup parameters supports the notion that evaporation and divergent tangential flow can cooperate to drive breakup. The resulting tear breakup is typically faster than purely evaporative cases. Many instances of tear breakup may have similar causes, which suggests that interpretation of experimental results may benefit from considering multiple mechanisms. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: 47 pages, 24 figures, submitted to Bulletin of Mathematical Biology

MSC Class: 76M21 (Primary) 76A20 (Secondary)

Journal ref: Bulletin of Mathematical Biology (2021)

Automatic Speech Recognition Benchmark for Air-Traffic Communications

Authors: Juan Zuluaga-Gomez, Petr Motlicek, Qingran Zhan, Karel Vesely, Rudolf Braun

Abstract: Advances in Automatic Speech Recognition (ASR) over the last decade opened new areas of speech-based automation such as in Air-Traffic Control (ATC) environment. Currently, voice communication and data links communications are the only way of contact between pilots and Air-Traffic Controllers (ATCo), where the former is the most widely used and the latter is a non-spoken method mandatory for ocean… ▽ More Advances in Automatic Speech Recognition (ASR) over the last decade opened new areas of speech-based automation such as in Air-Traffic Control (ATC) environment. Currently, voice communication and data links communications are the only way of contact between pilots and Air-Traffic Controllers (ATCo), where the former is the most widely used and the latter is a non-spoken method mandatory for oceanic messages and limited for some domestic issues. ASR systems on ATCo environments inherit increasing complexity due to accents from non-English speakers, cockpit noise, speaker-dependent biases, and small in-domain ATC databases for training. Hereby, we introduce CleanSky EC-H2020 ATCO2, a project that aims to develop an ASR-based platform to collect, organize and automatically pre-process ATCo speech-data from air space. This paper conveys an exploratory benchmark of several state-of-the-art ASR models trained on more than 170 hours of ATCo speech-data. We demonstrate that the cross-accent flaws due to speakers' accents are minimized due to the amount of data, making the system feasible for ATC environments. The developed ASR system achieves an averaged word error rate (WER) of 7.75% across four databases. An additional 35% relative improvement in WER is achieved on one test set when training a TDNNF system with byte-pair encoding. △ Less

Submitted 13 August, 2020; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Accepted to: 21st INTERSPEECH conference (Shanghai, October 25-29)

Anticipated Performance of the Square Kilometre Array -- Phase 1 (SKA1)

Authors: Robert Braun, Anna Bonaldi, Tyler Bourke, Evan Keane, Jeff Wagg

Abstract: The Square Kilometre Array (SKA), currently under design, will be a transformational facility for studying the Universe at centimetre and metre wavelengths in the next decade and beyond. This paper provides the current best estimate of the anticipated performance of SKA Phase 1 (SKA1), using detailed design work, before actual on-sky measurements have been made. It will be updated as new informati… ▽ More The Square Kilometre Array (SKA), currently under design, will be a transformational facility for studying the Universe at centimetre and metre wavelengths in the next decade and beyond. This paper provides the current best estimate of the anticipated performance of SKA Phase 1 (SKA1), using detailed design work, before actual on-sky measurements have been made. It will be updated as new information becomes available. The information contained in this paper takes precedent over any previous documents. △ Less

Submitted 29 December, 2019; originally announced December 2019.

Comments: 27 Pages, 10 Tables, 14 Figures

Long-Term Latency Measurement of Deployed Fiber

Authors: Florian Azendorf, Annika Dochhan, Ralf-Peter Braun, Michael Eiselt

Abstract: Using a Correlation-OTDR we measured the latency of fibers in a deployed cable and calculated the time coefficient of the fiber temperature changes. Annual temperature variations of 25K were estimated for the deployed fiber. Using a Correlation-OTDR we measured the latency of fibers in a deployed cable and calculated the time coefficient of the fiber temperature changes. Annual temperature variations of 25K were estimated for the deployed fiber. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Comments: This project has received funding from the European Union Horizon 2020 research and innovation programme under grant agreement No 762055 (BlueSpace Project)

Square Kilometre Array Science Data Challenge 1

Authors: Anna Bonaldi, Robert Braun

Abstract: The Square Kilometre Array (SKA, https://skatelescope.org) will be the world's largest radio telescope. SKA Science Data Challenges will be regularly issued to the community as part of the science preparatory activities. The purpose of these challenges is to inform the development of the data reduction workflows, to allow the science community to get familiar with the standard products the SKA wil… ▽ More The Square Kilometre Array (SKA, https://skatelescope.org) will be the world's largest radio telescope. SKA Science Data Challenges will be regularly issued to the community as part of the science preparatory activities. The purpose of these challenges is to inform the development of the data reduction workflows, to allow the science community to get familiar with the standard products the SKA will deliver, and optimise their analyses to extract science from them. These challenges may consist of real data from currently operating radio facilities or of simulated SKA data. The purpose of this document is to provide information on how the SKA Science data challenge #1 (SDC1) has been produced and to set the challenge for the community. For more information on how to take part in the challenge and to download the data see https://astronomers.skatelescope.org/ska-science-data-challenge-1/ △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: 7 pages, extracted from the document released by the SKA Office: SKA-TEL-SKO-00001001

Time-lagged Ordered Lasso for network inference

Authors: Phan Nguyen, Rosemary Braun

Abstract: Accurate gene regulatory networks can be used to explain the emergence of different phenotypes, disease mechanisms, and other biological functions. Many methods have been proposed to infer networks from gene expression data but have been hampered by problems such as low sample size, inaccurate constraints, and incomplete characterizations of regulatory dynamics. Since expression regulation is dyna… ▽ More Accurate gene regulatory networks can be used to explain the emergence of different phenotypes, disease mechanisms, and other biological functions. Many methods have been proposed to infer networks from gene expression data but have been hampered by problems such as low sample size, inaccurate constraints, and incomplete characterizations of regulatory dynamics. Since expression regulation is dynamic, time-course data can be used to infer causality, but these datasets tend to be short or sparsely sampled. In addition, temporal methods typically assume that the expression of a gene at a time point depends on the expression of other genes at only the immediately preceding time point, while other methods include additional time points without any constraints to account for their temporal distance. These limitations can contribute to inaccurate networks with many missing and anomalous links. We adapted the time-lagged Ordered Lasso, a regularized regression method with temporal monotonicity constraints, for \textit{de novo} reconstruction. We also developed a semi-supervised method that embeds prior network information into the Ordered Lasso to discover novel regulatory dependencies in existing pathways. We evaluated these approaches on simulated data for a repressilator, time-course data from past DREAM challenges, and a HeLa cell cycle dataset to show that they can produce accurate networks subject to the dynamics and assumptions of the time-lagged Ordered Lasso regression. △ Less

Submitted 10 December, 2018; v1 submitted 30 June, 2018; originally announced July 2018.

Single nucleotide polymorphisms that modulate microRNA regulation of gene expression in tumors

Authors: Gary Wilk, Rosemary Braun

Abstract: Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with trait diversity and disease susceptibility, yet the functional properties of many genetic variants and their molecular interactions remains unclear. It has been hypothesized that SNPs in microRNA binding sites may disrupt gene regulation by microRNAs (miRNAs), short non-coding RNAs that bi… ▽ More Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with trait diversity and disease susceptibility, yet the functional properties of many genetic variants and their molecular interactions remains unclear. It has been hypothesized that SNPs in microRNA binding sites may disrupt gene regulation by microRNAs (miRNAs), short non-coding RNAs that bind to mRNA and downregulate the target gene. While a number of studies have been conducted to predict the location of SNPs in miRNA binding sites, to date there has been no comprehensive analysis of how SNP variants may impact miRNA regulation of genes. Here we investigate the functional properties of genetic variants and their effects on miRNA regulation of gene expression in cancer. Our analysis is motivated by the hypothesis that distinct alleles may cause differential binding (from miRNAs to mRNAs or from transcription factors to DNA) and change the expression of genes. We previously identified pathways--systems of genes conferring specific cell functions--that are dysregulated by miRNAs in cancer, by comparing miRNA-pathway associations between healthy and tumor tissue. We draw on these results as a starting point to assess whether SNPs in genes on dysregulated pathways are responsible for miRNA dysregulation of individual genes in tumors. Using an integrative analysis that incorporates miRNA expression, mRNA expression, and SNP genotype data, we identify SNPs that appear to influence the association between miRNAs and genes, which we term "regulatory QTLs (regQTLs)": loci whose alleles impact the regulation of genes by miRNAs. We describe the method, apply it to analyze four cancer types (breast, liver, lung, prostate) using data from The Cancer Genome Atlas (TCGA), and provide a tool to explore the findings. △ Less

Submitted 20 March, 2018; v1 submitted 8 March, 2018; originally announced March 2018.

Integrative analysis reveals disrupted pathways regulated by microRNAs in cancer

Authors: Gary Wilk, Rosemary Braun

Abstract: MicroRNAs (miRNAs) are small endogenous regulatory molecules that modulate gene expression post-transcriptionally. Although differential expression of miRNAs have been implicated in many diseases (including cancers), the underlying mechanisms of action remain unclear. Because each miRNA can target multiple genes, miRNAs may potentially have functional implications for the overall behavior of entir… ▽ More MicroRNAs (miRNAs) are small endogenous regulatory molecules that modulate gene expression post-transcriptionally. Although differential expression of miRNAs have been implicated in many diseases (including cancers), the underlying mechanisms of action remain unclear. Because each miRNA can target multiple genes, miRNAs may potentially have functional implications for the overall behavior of entire pathways. Here we investigate the functional consequences of miRNA dysregulation through an integrative analysis of miRNA and mRNA expression data using a novel approach that incorporates pathway information a priori. By searching for miRNA-pathway associations that differ between healthy and tumor tissue, we identify specific relationships at the systems-level which are disrupted in cancer. Our approach is motivated by the hypothesis that if a miRNA and pathway are associated, then the expression of the miRNA and the collective behavior of the genes in a pathway will be correlated. As such, we first obtain an expression-based summary of pathway activity using Isomap, a dimension reduction method which can articulate nonlinear structure in high-dimensional data. We then search for miRNAs that exhibit differential correlations with the pathway summary between phenotypes as a means of finding aberrant miRNA-pathway coregulation in tumors. We apply our method to cancer data using gene and miRNA expression datasets from The Cancer Genome Atlas (TCGA) and compare ${\sim}10^5$ miRNA-pathway relationships between healthy and tumor samples from four tissues (breast, prostate, lung, and liver). Many of the flagged pairs we identify have a biological basis for disruption in cancer. △ Less

Submitted 18 August, 2017; originally announced August 2017.

Journal ref: Nucleic Acids Res. 2017

Semi-supervised network inference using simulated gene expression dynamics

Authors: Phan Nguyen, Rosemary Braun

Abstract: Motivation: Inferring the structure of gene regulatory networks from high--throughput datasets remains an important and unsolved problem. Current methods are hampered by problems such as noise, low sample size, and incomplete characterizations of regulatory dynamics, leading to networks with missing and anomalous links. Integration of prior network information (e.g., from pathway databases) has th… ▽ More Motivation: Inferring the structure of gene regulatory networks from high--throughput datasets remains an important and unsolved problem. Current methods are hampered by problems such as noise, low sample size, and incomplete characterizations of regulatory dynamics, leading to networks with missing and anomalous links. Integration of prior network information (e.g., from pathway databases) has the potential to improve reconstructions. Results: We developed a semi--supervised network reconstruction algorithm that enables the synthesis of information from partially known networks with time course gene expression data. We adapted PLS-VIP for time course data and used reference networks to simulate expression data from which null distributions of VIP scores are generated and used to estimate edge probabilities for input expression data. By using simulated dynamics to generate reference distributions, this approach incorporates previously known regulatory relationships and links the network to the dynamics to form a semi-supervised approach that discovers novel and anomalous connections. We applied this approach to data from a sleep deprivation study with KEGG pathways treated as prior networks, as well as to synthetic data from several DREAM challenges, and find that it is able to recover many of the true edges and identify errors in these networks, suggesting its ability to derive posterior networks that accurately reflect gene expression dynamics. △ Less

Submitted 1 December, 2017; v1 submitted 2 June, 2017; originally announced June 2017.

Network-based identification of disease genes in expression data: the GeneSurrounder method

Authors: Sahil D. Shah, Rosemary Braun

Abstract: The advent of high--throughput transcription profiling technologies has enabled identification of genes and pathways associated with disease, providing new avenues for precision medicine. A key challenge is to analyze this data in the context of the regulatory networks and pathways that control cellular processes, while still obtaining insights that can be used to design new diagnostic and therape… ▽ More The advent of high--throughput transcription profiling technologies has enabled identification of genes and pathways associated with disease, providing new avenues for precision medicine. A key challenge is to analyze this data in the context of the regulatory networks and pathways that control cellular processes, while still obtaining insights that can be used to design new diagnostic and therapeutic interventions. While classical differential expression analysis provides specific and hence targetable gene-level insights, it does not include any systems-level information. On the other hand, pathway analyses integrate systems-level information with expression data, but are often limited in their ability to indicate specific molecular targets. We introduce GeneSurrounder, an analysis method that takes into account the complex structure of interaction networks to identify specific genes that disrupt pathway activity in a disease-specific manner. GeneSurrounder integrates transcriptomic data and pathway network information in a novel two-step procedure to detect genes that (i) appear to influence the expression of other genes local to it in the network and (ii) are part of a subnetwork of differentially expressed genes. Combined, this evidence can be used to pinpoint specific genes that have a mechanistic role in the phenotype of interest. Applying GeneSurrounder to three distinct ovarian cancer studies using a global KEGG network, we show that our method is able to identify biologically relevant genes and genes missed by single-gene association tests, integrate pathway and expression data, and yield more consistent results across multiple studies of the same phenotype than competing methods. △ Less

Submitted 9 January, 2019; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: We have extended the application and evaluation of our GeneSurrounder method to a second disease (gene expression data from bladder cancer) and added additional analyses of GeneSurrounder's ability to identify known cancer-associated genes

Simulation of Parabolic Flow on an Eye-Shaped Domain with Moving Boundary

Authors: Tobin A. Driscoll, Richard J. Braun, Joseph K. Brosch

Abstract: During the upstroke of a normal eye blink, the upper lid moves and paints a thin tear film over the exposed corneal and conjunctival surfaces. This thin tear film may be modeled by a nonlinear fourth-order PDE derived from lubrication theory. A challenge in the numerical simulation of this model is to include both the geometry of the eye and the movement of the eyelid. A pair of orthogonal and con… ▽ More During the upstroke of a normal eye blink, the upper lid moves and paints a thin tear film over the exposed corneal and conjunctival surfaces. This thin tear film may be modeled by a nonlinear fourth-order PDE derived from lubrication theory. A challenge in the numerical simulation of this model is to include both the geometry of the eye and the movement of the eyelid. A pair of orthogonal and conformal maps transform a square into an approximate representation of the exposed ocular surface of a human eye. A spectral collocation method on the square produces relatively efficient solutions on the eye-shaped domain via these maps. The method is demonstrated on linear and nonlinear second-order diffusion equations and shown to have excellent accuracy as measured pointwise or by conservation checks. Future work will use the method for thin-film equations on the same type of domain. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: 19 pages, 9 figures

Computed Flow and Fluorescence Over the Ocular Surface

Authors: Longfei Li, R. J. Braun, W. D. Henshaw, P. E. King-Smith

Abstract: Fluorescein is perhaps the most commonly used substance to visualize tear film thickness and dynamics; better understanding of this process aids understanding of dry eye syndrome which afflicts millions of people. We study a mathematical model for tear film flow, evaporation, solutal transport and fluorescence over the exposed ocular surface during the interblink. Transport of the fluorescein ion… ▽ More Fluorescein is perhaps the most commonly used substance to visualize tear film thickness and dynamics; better understanding of this process aids understanding of dry eye syndrome which afflicts millions of people. We study a mathematical model for tear film flow, evaporation, solutal transport and fluorescence over the exposed ocular surface during the interblink. Transport of the fluorescein ion by fluid flow in the tear film affects the intensity of fluorescence via changes in concentration and tear film thickness. Evaporation causes increased osmolarity and potential irritation over the ocular surface; it also alters fluorescein concentration and thus fluorescence. Using thinning rates from in vivo measurements together with thin film equations for flow and transport of multiple solutes, we compute dynamic results for tear film quantities of interest. We compare our computed intensity distributions with in vivo observations. A number of experimental features are recovered by the model. △ Less

Submitted 31 October, 2016; originally announced October 2016.

Handicap principle implies emergence of dimorphic mating displays

Authors: Sara M. Clifton, Rosemary I. Braun, Daniel M. Abrams

Abstract: Species spanning the animal kingdom have evolved extravagant and costly ornaments to attract mating partners. Zahavi's handicap principle offers an elegant explanation for this: ornaments signal individual quality, and must be costly to ensure honest signalling, making mate selection more efficient. Here we incorporate the assumptions of the handicap principle into a mathematical model and show th… ▽ More Species spanning the animal kingdom have evolved extravagant and costly ornaments to attract mating partners. Zahavi's handicap principle offers an elegant explanation for this: ornaments signal individual quality, and must be costly to ensure honest signalling, making mate selection more efficient. Here we incorporate the assumptions of the handicap principle into a mathematical model and show that they are sufficient to explain the heretofore puzzling observation of bimodally distributed ornament sizes in a variety of species. △ Less

Submitted 31 October, 2016; v1 submitted 27 September, 2015; originally announced September 2015.

Comments: 7 pages, 4 figures, supplementary information included

Journal ref: Proceedings of the Royal Society B 283, 1970: November 30 2016

Discovery of HI gas in a young radio galaxy at $z = 0.44$ using the Australian Square Kilometre Array Pathfinder

Authors: J. R. Allison, E. M. Sadler, V. A. Moss, M. T. Whiting, R. W. Hunstead, M. B. Pracy, S. J. Curran, S. M. Croom, M. Glowacki, R. Morganti, S. S. Shabala, M. A. Zwaan, G. Allen, S. W. Amy, P. Axtens, L. Ball, K. W. Bannister, S. Barker, M. E. Bell, D. C. -J. Bock, R. Bolton, M. Bowen, B. Boyle, R. Braun, S. Broadhurst , et al. (78 additional authors not shown)

Abstract: We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The… ▽ More We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$. △ Less

Submitted 16 August, 2015; v1 submitted 4 March, 2015; originally announced March 2015.

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

Filaments of the radio cosmic web: opportunities and challenges for SKA

Authors: Franco Vazza, Chiara Ferrari, Annalisa Bonafede, Marcus Brüggen, Claudio Gheller, Robert Braun, Shea Brown

Abstract: The detection of the diffuse gas component of the cosmic web remains a formidable challenge. In this work we study synchrotron emission from the cosmic web with simulated SKA1 observations, which can represent an fundamental probe of the warm-hot intergalactic medium. We investigate radio emission originated by relativistic electrons accelerated by shocks surrounding cosmic filaments, assuming dif… ▽ More The detection of the diffuse gas component of the cosmic web remains a formidable challenge. In this work we study synchrotron emission from the cosmic web with simulated SKA1 observations, which can represent an fundamental probe of the warm-hot intergalactic medium. We investigate radio emission originated by relativistic electrons accelerated by shocks surrounding cosmic filaments, assuming diffusive shock acceleration and as a function of the (unknown) large-scale magnetic fields. The detection of the brightest parts of large ($>10 \rm Mpc$) filaments of the cosmic web should be within reach of the SKA1-LOW, if the magnetic field is at the level of a $\sim 10$ percent equipartition with the thermal gas, corresponding to $\sim 0.1 μG$ for the most massive filaments in simulations. In the course of a 2-years survey with SKA1-LOW, this will enable a first detection of the "tip of the iceberg" of the radio cosmic web, and allow for the use of the SKA as a powerful tool to study the origin of cosmic magnetism in large-scale structures. On the other hand, the SKA1-MID and SKA1-SUR seem less suited for this science case at low redshift ($z \leq 0.4$), owing to the missing short baselines and the consequent lack of signal from the large-scale brightness fluctuations associated with the filaments. In this case only very long exposures ($\sim 1000$ hr) may enable the detection of $\sim 1-2$ filament for field of view in the SKA1-SUR PAF Band1. △ Less

Submitted 1 January, 2015; originally announced January 2015.

Comments: 10 pages, 4 figures, to appear in Proceedings of 'Advancing Astrophysics with the SKA (AASKA14) - Cosmic Magnetism' Chapters

Network Methods for Pathway Analysis of Genomic Data

Authors: Rosemary Braun, Sahil Shah

Abstract: Rapid advances in high-throughput technologies have led to considerable interest in analyzing genome-scale data in the context of biological pathways, with the goal of identifying functional systems that are involved in a given phenotype. In the most common approaches, biological pathways are modeled as simple sets of genes, neglecting the network of interactions comprising the pathway and treatin… ▽ More Rapid advances in high-throughput technologies have led to considerable interest in analyzing genome-scale data in the context of biological pathways, with the goal of identifying functional systems that are involved in a given phenotype. In the most common approaches, biological pathways are modeled as simple sets of genes, neglecting the network of interactions comprising the pathway and treating all genes as equally important to the pathway's function. Recently, a number of new methods have been proposed to integrate pathway topology in the analyses, harnessing existing knowledge and enabling more nuanced models of complex biological systems. However, there is little guidance available to researches choosing between these methods. In this review, we discuss eight topology-based methods, comparing their methodological approaches and appropriate use cases. In addition, we present the results of the application of these methods to a curated set of ten gene expression profiling studies using a common set of pathway annotations. We report the computational efficiency of the methods and the consistency of the results across methods and studies to help guide users in choosing a method. We also discuss the challenges and future outlook for improved network analysis methodologies. △ Less

Submitted 7 November, 2014; originally announced November 2014.

Comments: Review article

The Australian Square Kilometre Array Pathfinder: System Architecture and Specifications of the Boolardy Engineering Test Array

Authors: A. W. Hotan, J. D. Bunton, L. Harvey-Smith, B. Humphreys, B. D. Jeffs, T. Shimwell, J. Tuthill, M. Voronkov, G. Allen, S. Amy, K. Ardern, P. Axtens, L. Ball, K. Bannister, S. Barker, T. Bateman, R. Beresford, D. Bock, R. Bolton, M. Bowen, B. Boyle, R. Braun, S. Broadhurst, D. Brodrick, K. Brooks , et al. (76 additional authors not shown)

Abstract: This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interfe… ▽ More This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interferometer, fitted with prototype signal processing hardware capable of forming at least 9 dual-polarisation beams simultaneously, allowing several square degrees to be imaged in a single pointed observation. The main purpose of the test array is to develop beamforming and wide-field calibration methods for use with the full telescope, but it will also be capable of limited early science demonstrations. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Comments: 17 pages, 12 figures, accepted for publication in PASA

A Kinematic Measurement of Ram Pressure in the Outer Disk of Regular Galaxies

Authors: Sebastian Haan, Robert Braun

Abstract: While most ram pressure studies have focused on ram pressure stripping in galaxy clusters, we devise a novel approach based on a kinematic measurement of ram pressure perturbations in HI velocity fields for intergalactic material (IGM) densities and relative velocities that are one to two orders of magnitude lower than in galaxies showing ram pressure stripping. Our model evaluates ram pressure in… ▽ More While most ram pressure studies have focused on ram pressure stripping in galaxy clusters, we devise a novel approach based on a kinematic measurement of ram pressure perturbations in HI velocity fields for intergalactic material (IGM) densities and relative velocities that are one to two orders of magnitude lower than in galaxies showing ram pressure stripping. Our model evaluates ram pressure induced kinematic terms in gas disks with constant inclination as well as those with a warped geometry. Ram pressure perturbations are characterized by kinematic modes of even order, m=0 and m=2, corresponding to a ram wind perpendicular and parallel to the gas disk, respectively. Long-term consequences of ram pressure, such as warped disks as well as uncertainties in the disk geometry typically generate uneven modes (m=1 and m=3), that are clearly distinguishable from the kinematic ram pressure terms. We have applied our models to three nearby isolated galaxies, utilizing Markov Chain Monte Carlo fitting routines to determine ram pressure perturbations in the velocity fields of NGC 6946 and NGC 3621 of ~30km s$^{-1}$ (effective line-of-sight velocity change) at HI column densities below (4-10)$\times$10$^{20}$cm$^{-2}$ (at radial scales greater than ~15kpc). In contrast, NGC 628 is dominated by a strongly warped disk. Our model fits reveal the three-dimensional vector of the galaxies' movement with respect to the IGM rest-frame and provide constraints on the product of speed with IGM density, opening a new window for extragalactic velocity measurements and studies of the intergalactic medium. △ Less

Submitted 12 June, 2014; originally announced June 2014.

Comments: 34 pages, 28 figures, accepted for publication in MNRAS

The Radio Continuum-Star Formation Rate Relation in WSRT SINGS Galaxies

Authors: Volker Heesen, Elias Brinks, Adam K. Leroy, George Heald, Robert Braun, Frank Bigiel, Rainer Beck

Abstract: We study the spatially resolved Radio Continuum-Star Formation Rate (RC-SFR) relation using state-of-the-art star-formation (SF) tracers in a sample of 17 THINGS galaxies. We use hybrid Sigma_SFR maps (GALEX FUV plus Spitzer 24 mu), RC maps at 22/18 cm from the WSRT SINGS survey, and H-alpha maps to correct for thermal RC emission. We compare azimuthally averaged radial profiles of the RC and FUV/… ▽ More We study the spatially resolved Radio Continuum-Star Formation Rate (RC-SFR) relation using state-of-the-art star-formation (SF) tracers in a sample of 17 THINGS galaxies. We use hybrid Sigma_SFR maps (GALEX FUV plus Spitzer 24 mu), RC maps at 22/18 cm from the WSRT SINGS survey, and H-alpha maps to correct for thermal RC emission. We compare azimuthally averaged radial profiles of the RC and FUV/MIR-based Sigma_SFR maps and study pixel-by-pixel correlations at fixed linear scales of 1.2 and 0.7 kpc. The ratio of the integrated SFRs from the RC emission to that of the FUV/MIR-based SF tracers is R_int = 0.78 +/- 0.38, consistent with Condon's relation. We find a tight correlation between the radial profiles of the radio and FUV/MIR-based Sigma_SFR for the entire extent of the disk. The ratio R of the azimuthally averaged radio to FUV/MIR-based Sigma_SFR agrees with the integrated ratio with only small quasi-random fluctuations as function of radius. Pixel-by-pixel plots show a tight correlation in log-log diagrams of radio to FUV/MIR-based Sigma_SFR, with a typical standard deviation of a factor of two. Averaged over our sample we find (Sigma_SFR)_RC ~ (Sigma_SFR)_hyb^{0.63+/-0.25} implying that data points with high Sigma_SFR are relatively radio dim, whereas the reverse is true for low Sigma_SFR. We interpret this as a result of spectral ageing of CRe, which is supported by the radio spectral index: data points dominated by young CRe are relatively radio dim, those dominated by old CRe are relatively radio bright. The ratio of radio to FUV/MIR-based integrated SFR is independent of global galaxy parameters, suggesting that we can use RC emission as a universal SF tracer for galaxies, if we restrict ourselves to global or azimuthally averaged measurements. A magnetic field-SFR relation, B ~ SFR_hyb^{0.30+/-0.02}, holding both globally and locally, can explain our results. (abridged) △ Less

Submitted 4 April, 2014; v1 submitted 7 February, 2014; originally announced February 2014.

Comments: 41 pages, 30 figures, 3 tables, AJ accepted, pages 22-41 contain appendices with image atlas, replaced to match published version

Journal ref: AJ, 147, 103 (2014)

On the Formation of Warped Gas Disks in Galaxies

Authors: Sebastian Haan, Robert Braun

Abstract: We consider the most commonly occurring circumstances which apply to galaxies, namely membership in galaxy groups of about $10^{13}h^{-1} M_\odot$ total mass, and estimate the accompanying physical conditions of intergalactic medium (IGM) density and the relative galaxy-IGM space velocity. We then investigate the dynamical consequences of such a typical galaxy-IGM interaction on a rotating gaseous… ▽ More We consider the most commonly occurring circumstances which apply to galaxies, namely membership in galaxy groups of about $10^{13}h^{-1} M_\odot$ total mass, and estimate the accompanying physical conditions of intergalactic medium (IGM) density and the relative galaxy-IGM space velocity. We then investigate the dynamical consequences of such a typical galaxy-IGM interaction on a rotating gaseous disk within the galaxy potential. We find that the rotating outer disk is systematically distorted into a characteristic "warp" morphology, of the type that has been well-documented in the majority of well-studied nearby systems. The distortion is established rapidly, within two rotation periods, and is long-lived, surviving for at least ten. A second consequence of the interaction is the formation of a one arm retrograde spiral wave pattern that propagates in the disk. We suggest that the ubiquity of the warp phenomenon might be used to reconstruct both the IGM density profile and individual member orbits within galaxy groups. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: 7 pages, 7 figures, accepted for publication in MNRAS Letters. For supplementary material, see http://www.atnf.csiro.au/people/Sebastian.Haan/publications.html

The temperature of the diffuse HI in the Milky Way I: High resolution HI 21 cm absorption studies

Authors: Nirupam Roy, Nissim Kanekar, Robert Braun, Jayaram N. Chengalur

Abstract: We have carried out deep, high velocity resolution, interferometric Galactic HI 21 cm absorption spectroscopy towards 32 compact extra-galactic radio sources with the Giant Metrewave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). The optical depth spectra for most sources have root mean square noise values \lesssim 10^{-3} per 1 km/s velocity channel and are thus suffi… ▽ More We have carried out deep, high velocity resolution, interferometric Galactic HI 21 cm absorption spectroscopy towards 32 compact extra-galactic radio sources with the Giant Metrewave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). The optical depth spectra for most sources have root mean square noise values \lesssim 10^{-3} per 1 km/s velocity channel and are thus sufficiently sensitive to detect absorption by warm neutral hydrogen with HI column densities N_HI \gtrsim 10^{20} cm^{-2}, spin temperatures T_s \leq 5000 K, and line widths equal to the thermal width (20 km/s). HI 21 cm absorption was detected against all background sources but one, B0438-436. The spectra of sources observed separately with GMRT and WSRT show excellent agreement, indicating that spectral baseline problems and contamination from HI 21 cm emission are negligible. This paper presents the absorption spectra, the emission spectra along neighbouring sightlines from the Leiden-Argentine-Bonn survey, and the derived spin temperature spectra. On every sightline, the maximum spin temperature detected (at \ge 3 sigma significance) even at a velocity resolution of 1 km/s is \gtrsim 1000 K, indicating that we are detecting the warm neutral medium along most sightlines. This is by far the largest sample of Galactic HI 21 cm absorption spectra of this quality, providing a sensitive probe of physical conditions in the neutral atomic ISM. △ Less

Submitted 16 September, 2013; originally announced September 2013.

Comments: 14 pages, 5 figures. Accepted for publication in MNRAS. The definitive version will be available at http://mnras.oxfordjournals.org/

Journal ref: 2013, MNRAS, 436, 2352