Showing 1–50 of 70 results for author: Reid, T

Improving tracking algorithms with machine learning: a case for line-segment tracking at the High Luminosity LHC

Authors: Jonathan Guiang, Slava Krutelyov, Manos Vourliotis, Yanxi Gu, Avi Yagil, Balaji Venkat Sathia Narayanan, Matevz Tadel, Philip Chang, Mayra Silva, Gavin Niendorf, Peter Wittich, Tres Reid, Peter Elmer

Abstract: In this work, we present a study on ways that tracking algorithms can be improved with machine learning (ML). We base this study on the line segment tracking (LST) algorithm that we have designed to be naturally parallelized and vectorized in order to efficiently run on modern processors. LST has been developed specifically for the CMS Experiment at the LHC, towards the High Luminosity LHC (HL-LHC… ▽ More In this work, we present a study on ways that tracking algorithms can be improved with machine learning (ML). We base this study on the line segment tracking (LST) algorithm that we have designed to be naturally parallelized and vectorized in order to efficiently run on modern processors. LST has been developed specifically for the CMS Experiment at the LHC, towards the High Luminosity LHC (HL-LHC) upgrade. Moreover, we have already shown excellent efficiency and performance results as we iteratively improve LST, leveraging a full simulation of the CMS detector. At the same time, promising deep-learning-based tracking algorithms, such as Graph Neural Networks (GNNs), are being pioneered on the simplified TrackML dataset. These results suggest that parts of LST could be improved or replaced by ML. Thus, a thorough, step-by-step investigation of exactly how and where ML can be utilized, while still meeting realistic HL-LHC performance and efficiency constraints, is implemented as follows. First, a lightweight neural network is used to replace and improve upon explicitly defined track quality selections. This neural network is shown to be highly efficient and robust to displaced tracks while having little-to-no impact on the runtime of LST. These results clearly establish that ML can be used to improve LST without penalty. Next, exploratory studies of GNN track-building algorithms are described. In particular, low-level track objects from LST are considered as nodes in a graph, where edges represent higher-level objects or even entire track candidates. Then, an edge-classifier GNN is trained, and the efficiency of the resultant edge scores is compared with that of the existing LST track quality selections. These GNN studies provide insights into the practicality and performance of using more ambitious and complex ML algorithms for HL-LHC tracking at the CMS Experiment. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: Connecting the Dots Workshop (CTD 2023)

Graph Configurations and Independent Bondage Numbers of Planar Graphs

Authors: E. G. K. M. Gamlath, Bing Wei, Talmage James Reid

Abstract: The independent domination number of a finite graph G is the minimum cardinality of an independent dominating set of vertices. The independent bondage number of G is the minimum cardinality of a set of edges whose deletion results in a graph with a larger independent domination number than that of G. In this research, we enhance the existing upper bound on the independent bondage number of a plana… ▽ More The independent domination number of a finite graph G is the minimum cardinality of an independent dominating set of vertices. The independent bondage number of G is the minimum cardinality of a set of edges whose deletion results in a graph with a larger independent domination number than that of G. In this research, we enhance the existing upper bound on the independent bondage number of a planar graph with a minimum degree of at least three by identifying specific configurations within such planar graphs. △ Less

Submitted 3 February, 2024; originally announced February 2024.

MSC Class: 05C35 ACM Class: G.2.2

Latest Development of Electropolishing Optimization for 650 MHz Niobium Cavity

Authors: V. Chouhan, D. Bice, D. Burk, S. Chandrasekaran, A. Cravatta, P. Dubiel, G. V. Eremeev, F. Furuta, O. Melnychuk, A. Netepenko, M. K. Ng, J. Ozelis, H. Park, T. Ring, G. Wu, B. Guilfoyle, M. P. Kelly, T. Reid

Abstract: Electropolishing (EP) of 1.3 GHz niobium superconducting RF cavities is conducted to achieve a desired smooth and contaminant-free surface that yields good RF performance. Achieving a smooth surface of a large-sized elliptical cavity with the standard EP conditions was found to be challenging. This work aimed to conduct a systematic parametric EP study to understand the effects of various EP param… ▽ More Electropolishing (EP) of 1.3 GHz niobium superconducting RF cavities is conducted to achieve a desired smooth and contaminant-free surface that yields good RF performance. Achieving a smooth surface of a large-sized elliptical cavity with the standard EP conditions was found to be challenging. This work aimed to conduct a systematic parametric EP study to understand the effects of various EP parameters on the surface of 650 MHz niobium cavities used in the Proton Improvement Plan-II (PIP-II) linear accelerator. Parameters optimized in this study provided a smooth surface of the cavities. The electropolished cavity showed significantly a higher accelerating gradient meeting baseline requirement and qualified for further surface treatment to improve the cavity quality factor. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Comments: SRF2023

Report number: FERMILAB-CONF-23-274-TD

Generalizing mkFit and its Application to HL-LHC

Authors: Giuseppe Cerati, Peter Elmer, Patrick Gartung, Leonardo Giannini, Matti Kortelainen, Vyacheslav Krutelyov, Steven Lantz, Mario Masciovecchio, Tres Reid, Allison Reinsvold Hall, Daniel Riley, Matevz Tadel, Emmanouil Vourliotis, Peter Wittich, Avi Yagil

Abstract: mkFit is an implementation of the Kalman filter-based track reconstruction algorithm that exploits both thread- and data-level parallelism. In the past few years the project transitioned from the R&D phase to deployment in the Run-3 offline workflow of the CMS experiment. The CMS tracking performs a series of iterations, targeting reconstruction of tracks of increasing difficulty after removing hi… ▽ More mkFit is an implementation of the Kalman filter-based track reconstruction algorithm that exploits both thread- and data-level parallelism. In the past few years the project transitioned from the R&D phase to deployment in the Run-3 offline workflow of the CMS experiment. The CMS tracking performs a series of iterations, targeting reconstruction of tracks of increasing difficulty after removing hits associated to tracks found in previous iterations. mkFit has been adopted for several of the tracking iterations, which contribute to the majority of reconstructed tracks. When tested in the standard conditions for production jobs, speedups in track pattern recognition are on average of the order of 3.5x for the iterations where it is used (3-7x depending on the iteration). Multiple factors contribute to the observed speedups, including vectorization and a lightweight geometry description, as well as improved memory management and single precision. Efficient vectorization is achieved with both the icc and the gcc (default in CMSSW) compilers and relies on a dedicated library for small matrix operations, Matriplex, which has recently been released in a public repository. While the mkFit geometry description already featured levels of abstraction from the actual Phase-1 CMS tracker, several components of the implementations were still tied to that specific geometry. We have further generalized the geometry description and the configuration of the run-time parameters, in order to enable support for the Phase-2 upgraded tracker geometry for the HL-LHC and potentially other detector configurations. The implementation strategy and high-level code changes required for the HL-LHC geometry are presented. Speedups in track building from mkFit imply that track fitting becomes a comparably time consuming step of the tracking chain. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Advanced surface treatments for medium-velocity superconducting RF cavities for high accelerating gradient continuous-wave operation

Authors: K. McGee, S. Kim, K. Elliott, A. Ganshyn, W. Hartung, P. Ostroumov, A. Taylor, T. Xu, M. Martinello, G. V. Eremeev, A. Netepenko, F. Furuta, O. Melnychuk, M. P. Kelly, B. Guilfoyle, T. Reid

Abstract: Nitrogen-doping and furnace-baking are advanced high-Q0 recipes developed for 1.3 GHz TESLA-type cavities. These treatments will significantly benefit the high-Q0 linear accelerator community if they can be successfully adapted to different cavity styles and frequencies. Strong frequency- and geometry- dependence of these recipes makes the technology transfer amongst different cavity styles and fr… ▽ More Nitrogen-doping and furnace-baking are advanced high-Q0 recipes developed for 1.3 GHz TESLA-type cavities. These treatments will significantly benefit the high-Q0 linear accelerator community if they can be successfully adapted to different cavity styles and frequencies. Strong frequency- and geometry- dependence of these recipes makes the technology transfer amongst different cavity styles and frequencies far from straightforward, and requires rigorous study. Upcoming high-Q0 continuous-wave linear accelerator projects, such as the proposed Michigan State University Facility for Rare Isotope Beam Energy Upgrade, and the underway Fermilab's Proton Improvement Plan-II, could benefit enormously from adapting these techniques to their beta_opt = 0.6 ~650 MHz 5-cell elliptical superconducting rf cavities, operating at an accelerating gradient of around ~17 MV/m. This is the first investigation of the adaptation of nitrogen doping and medium temperature furnace baking to prototype 644 MHz beta_opt = 0.65 cavities, with the aim of demonstrating the high-Q0 potential of these recipes in these novel cavities for future optimization as part of the FRIB400 project R&D. We find that nitrogen-doping delivers superior Q0, despite the sub-GHz operating frequency of these cavities, but is sensitive to the post-doping electropolishing removal step and experiences elevated residual resistance. Medium temperature furnace baking delivers reasonable performance with decreased residual resistance compared to the nitrogen doped cavity, but may require further recipe refinement. The gradient requirement for the FRIB400 upgrade project is comfortably achieved by both recipes. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 16 pages, 5 figures

Report number: FERMILAB-PUB-23-377-TD

Demonstration of Niobium Tin in 218 MHz Low-beta Quarter Wave Accelerator Cavity

Authors: T. B. Petersen, M. P. Kelly, T. Reid, M. Kedzie, B. Guilfoyle, G. Chen, S. Posen, B. Tennis, G. Eremeev

Abstract: A 218 MHz quarter wave niobium cavity has been fabricated for the purpose of demonstrating Nb3Sn technology on a low-beta accelerator cavity. Niobiumtin has been established as a promising next generation SRF material, but development has focused primarily in high-beta elliptical cell cavities. This material has a significantly higher TC than niobium, allowing for design of higher frequency quarte… ▽ More A 218 MHz quarter wave niobium cavity has been fabricated for the purpose of demonstrating Nb3Sn technology on a low-beta accelerator cavity. Niobiumtin has been established as a promising next generation SRF material, but development has focused primarily in high-beta elliptical cell cavities. This material has a significantly higher TC than niobium, allowing for design of higher frequency quarter wave cavities (that are subsequently smaller) as well as for significantly lowered cooling requirements (possibly leading to cryocooler based designs). The fabrication, initial cold testing, and Nb3Sn coating are discussed as well as test plans and details of future applications. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 21st Intl Conf Radio Frequency Superconductivity (SRF 2023)

Report number: FERMILAB-CONF-23-344-TD

Continuous Ultraviolet to Blue-Green Astrocomb

Authors: Yuk Shan Cheng, Kamalesh Dadi, Toby Mitchell, Samantha Thompson, Nikolai Piskunov, Lewis D. Wright, Corin B. E. Gawith, Richard A. McCracken, Derryck T. Reid

Abstract: The characterization of Earth-like exoplanets and precision tests of cosmological models using next-generation telescopes such as the ELT will demand precise calibration of astrophysical spectrographs in the visible region, where stellar absorption lines are most abundant. Astrocombs--lasers providing a broadband sequence of ultra-narrow, drift-free, regularly spaced optical frequencies on a multi… ▽ More The characterization of Earth-like exoplanets and precision tests of cosmological models using next-generation telescopes such as the ELT will demand precise calibration of astrophysical spectrographs in the visible region, where stellar absorption lines are most abundant. Astrocombs--lasers providing a broadband sequence of ultra-narrow, drift-free, regularly spaced optical frequencies on a multi-GHz grid--promise an atomically-traceable, versatile calibration scale, but their realization is challenging because of the need for ultra-broadband frequency conversion of mode-locked infrared lasers into the blue-green region. Here, we introduce a new concept achieving a broad, continuous spectrum by combining second-harmonic generation and sum-frequency-mixing in an aperiodically-poled MgO:PPLN waveguide to generate gap-free 390-520 nm light from a 1 GHz Ti:sapphire laser frequency comb. We lock a low-dispersion Fabry-Perot etalon to extract a sub-comb of bandwidth from 392-472 nm with a spacing of 30 GHz, visualizing the thousands of resulting comb modes on a high resolution cross-dispersion spectrograph. Complementary experimental data and simulations demonstrate the effectiveness of the approach for eliminating the spectral gaps present in second-harmonic-only conversion, in which weaker fundamental frequencies are suppressed by the quadratic \{chi}^((2)) nonlinearity. Requiring only ~100 pJ pulse energies, our concept establishes a practical new route to broadband UV-visible generation at GHz repetition rates. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 14 pages; 4 figures

Journal ref: NATURE COMMUNICATIONS 15:1466 (2024)

Inducing Quantum Phase Transitions in Non-Topological Insulators Via Atomic Control of Sub-Structural Elements

Authors: Thomas K. Reid, S. Pamir Alpay, Alexander V. Balatsky, Sanjeev K. Nayak

Abstract: Topological insulators (TIs) are an important family of quantum materials that exhibit a Dirac point (DP) in the surface band structure but have a finite band gap in bulk. A large degree of spin-orbit interaction and low bandgap is a prerequisite for stabilizing DPs on selective atomically flat cleavage planes. Tuning of the DP in these materials has been suggested via modifications to the atomic… ▽ More Topological insulators (TIs) are an important family of quantum materials that exhibit a Dirac point (DP) in the surface band structure but have a finite band gap in bulk. A large degree of spin-orbit interaction and low bandgap is a prerequisite for stabilizing DPs on selective atomically flat cleavage planes. Tuning of the DP in these materials has been suggested via modifications to the atomic structure of the entire system. Using the example of As$_2$Te$_3$ and ZnTe$_5$, which are not TIs, we show that a quantum phase transition can be induced in atomically flat and stepped surfaces, for As$_2$Te$_3$ and ZrTe$_5$, respectively. This is achieved by establishing a framework for controlling electronic properties that is focused on local perturbations at key locations that we call sub-structural elements (SSEs). We exemplify this framework through a novel method of isovalent sublayer anion doping and biaxial strain. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 13 pages, 3 figures, 1 table

Journal ref: Phys. Rev. B 108, 235402 (2023)

Trust in Shared Automated Vehicles: Study on Two Mobility Platforms

Authors: Shashank Mehrotra, Jacob G Hunter, Matthew Konishi, Kumar Akash, Zhaobo Zheng, Teruhisa Misu, Anil Kumar, Tahira Reid, Neera Jain

Abstract: The ever-increasing adoption of shared transportation modalities across the United States has the potential to fundamentally change the preferences and usage of different mobilities. It also raises several challenges with respect to the design and development of automated mobilities that can enable a large population to take advantage of this emergent technology. One such challenge is the lack of… ▽ More The ever-increasing adoption of shared transportation modalities across the United States has the potential to fundamentally change the preferences and usage of different mobilities. It also raises several challenges with respect to the design and development of automated mobilities that can enable a large population to take advantage of this emergent technology. One such challenge is the lack of understanding of how trust in one automated mobility may impact trust in another. Without this understanding, it is difficult for researchers to determine whether future mobility solutions will have acceptance within different population groups. This study focuses on identifying the differences in trust across different mobility and how trust evolves across their use for participants who preferred an aggressive driving style. A dual mobility simulator study was designed in which 48 participants experienced two different automated mobilities (car and sidewalk). The results found that participants showed increasing levels of trust when they transitioned from the car to the sidewalk mobility. In comparison, participants showed decreasing levels of trust when they transitioned from the sidewalk to the car mobility. The findings from the study help inform and identify how people can develop trust in future mobility platforms and could inform the design of interventions that may help improve the trust and acceptance of future mobility. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: https://trid.trb.org/view/2117834

Journal ref: Transportation Research Board 102nd Annual Meeting, Washington DC, United States, 1-12 Jan 2023, No. TRBAM-23-04456. 2023

The ESO's Extremely Large Telescope Working Groups

Authors: Paolo Padovani, Michele Cirasuolo, Remco van der Burg, Faustine Cantalloube, Elizabeth George, Markus Kasper, Kieran Leschinski, Carlos Martins, Julien Milli, Sabine Möhler, Mark Neeser, Benoit Neichel, Angel Otarola, Rubén Sánchez-Janssen, Benoit Serra, Alain Smette, Elena Valenti, Christophe Verinaud, Joël Vernet, Olivier Absil, Guido Agapito, Morten Andersen, Carmelo Arcidiacono, Matej Arko, Pierre Baudoz , et al. (60 additional authors not shown)

Abstract: Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations an… ▽ More Since 2005 ESO has been working with its community and industry to develop an extremely large optical/infrared telescope. ESO's Extremely Large Telescope, or ELT for short, is a revolutionary ground-based telescope that will have a 39-metre main mirror and will be the largest visible and infrared light telescope in the world. To address specific topics that are needed for the science operations and calibrations of the telescope, thirteen specific working groups were created to coordinate the effort between ESO, the instrument consortia, and the wider community. We describe here the goals of these working groups as well as their achievements so far. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 8 pages, 4 figures

Journal ref: The Messenger, vol. 189, (December 2022) p. 23-30

Hong-Ou-Mandel interference with a diode-pumped 1-GHz Ti:sapphire laser

Authors: Imogen Morland, Hanna Ostapenko, Feng Zhu, Derryck T. Reid, Jonathan Leach

Abstract: Correlated photon pairs generated through spontaneous parametric down-conversion (SPDC) are a key resource in quantum optics. In many quantum optics applications, such as satellite quantum key distribution (QKD), a compact, high repetition rate pump laser is required. Here we demonstrate the use of a compact, GHz-rate diode-pumped three-element Kerr-lens-modelocked Ti:sapphire laser for the genera… ▽ More Correlated photon pairs generated through spontaneous parametric down-conversion (SPDC) are a key resource in quantum optics. In many quantum optics applications, such as satellite quantum key distribution (QKD), a compact, high repetition rate pump laser is required. Here we demonstrate the use of a compact, GHz-rate diode-pumped three-element Kerr-lens-modelocked Ti:sapphire laser for the generation of correlated photon pairs at 790 nm. We verify the presence of indistinguishable photons produced via SPDC using Hong-Ou-Mandel (HOM) interferometry and observe a dip in coincidence counts with a visibility of 81.8\%. △ Less

Submitted 14 February, 2023; originally announced February 2023.

Comments: 5 pages, 3 figures

Localization & Mapping Requirements for Level 2+ Autonomous Vehicles

Authors: Tyler G. R. Reid, Andrew Neish, Brian Manning

Abstract: Autonomous vehicles are being deployed with a spectrum of capability, extending from driver assistance features for the highway in personal vehicles (SAE Level 2+) to fully autonomous fleet ride sharing services operating in complex city environments (SAE Level 4+). This spectrum of autonomy often operates in different physical environments with different degrees of assumed driver in-the-loop over… ▽ More Autonomous vehicles are being deployed with a spectrum of capability, extending from driver assistance features for the highway in personal vehicles (SAE Level 2+) to fully autonomous fleet ride sharing services operating in complex city environments (SAE Level 4+). This spectrum of autonomy often operates in different physical environments with different degrees of assumed driver in-the-loop oversight and hence have very different system and subsystem requirements. At the heart of SAE Level 2 to 5 systems is localization and mapping, which ranges from road determination for feature geofencing or high-level routing, through lane determination for advanced driver assistance, to where-in-lane positioning for full vehicle control. We assess localization and mapping requirements for different levels of autonomy and supported features. This work provides a framework for system decomposition, including the level of redundancy needed to achieve the target level of safety. We examine several representative autonomous and assistance features and make recommendations on positioning requirements as well map georeferencing and information integrity. △ Less

Submitted 5 January, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: ION ITM 2023

Web Archiving as Entertainment

Authors: Travis Reid, Michael L. Nelson, Michele C. Weigle

Abstract: We want to make web archiving entertaining so that it can be enjoyed like a spectator sport. To this end, we have been working on a proof of concept that involves gamification of the web archiving process and integrating video games and web archiving. Our vision for this proof of concept involves a web archiving live stream and a gaming live stream. We are creating web archiving live streams that… ▽ More We want to make web archiving entertaining so that it can be enjoyed like a spectator sport. To this end, we have been working on a proof of concept that involves gamification of the web archiving process and integrating video games and web archiving. Our vision for this proof of concept involves a web archiving live stream and a gaming live stream. We are creating web archiving live streams that make the web archiving process more transparent to viewers by live streaming the web archiving and replay sessions to video game live streaming platforms like Twitch, Facebook Gaming, and YouTube. We also want to live stream gameplay from games where the gameplay is influenced by web archiving and replay performance. So far we have created web archiving live streams that show the web archiving and replay sessions for two web archive crawlers and gaming live streams that show gameplay influenced by the web archiving performance from the web archiving live stream. We have also applied the gaming concept of speedruns, where a player attempts to complete a game as quickly as possible. This could make a web archiving live stream more entertaining, because we can have a competition between two crawlers to see which crawler is faster at archiving a set of URIs. △ Less

Submitted 3 November, 2022; originally announced November 2022.

Comments: This is an extended version of a paper from ICADL 2022. 20 pages and 10 figures

The Interaction Gap: A Step Toward Understanding Trust in Autonomous Vehicles Between Encounters

Authors: Jacob G. Hunter, Matthew Konishi, Neera Jain, Kumar Akash, Xingwei Wu, Teruhisa Misu, Tahira Reid

Abstract: Shared autonomous vehicles (SAVs) will be introduced in greater numbers over the coming decade. Due to rapid advances in shared mobility and the slower development of fully autonomous vehicles (AVs), SAVs will likely be deployed before privately-owned AVs. Moreover, existing shared mobility services are transitioning their vehicle fleets toward those with increasingly higher levels of driving auto… ▽ More Shared autonomous vehicles (SAVs) will be introduced in greater numbers over the coming decade. Due to rapid advances in shared mobility and the slower development of fully autonomous vehicles (AVs), SAVs will likely be deployed before privately-owned AVs. Moreover, existing shared mobility services are transitioning their vehicle fleets toward those with increasingly higher levels of driving automation. Consequently, people who use shared vehicles on an "as needed" basis will have infrequent interactions with automated driving, thereby experiencing interaction gaps. Using human trust data of 25 participants, we show that interaction gaps can affect human trust in automated driving. Participants engaged in a simulator study consisting of two interactions separated by a one-week interaction gap. A moderate, inverse correlation was found between the change in trust during the initial interaction and the interaction gap, suggesting people "forget" some of their gained trust or distrust in automation during an interaction gap. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 5 pages, 3 figures

Journal ref: Proceedings of the Human Factors and Ergonomics Society Annual Meeting, 2022, 66(1), 147-151

SmOOD: Smoothness-based Out-of-Distribution Detection Approach for Surrogate Neural Networks in Aircraft Design

Authors: Houssem Ben Braiek, Ali Tfaily, Foutse Khomh, Thomas Reid, Ciro Guida

Abstract: Aircraft industry is constantly striving for more efficient design optimization methods in terms of human efforts, computation time, and resource consumption. Hybrid surrogate optimization maintains high results quality while providing rapid design assessments when both the surrogate model and the switch mechanism for eventually transitioning to the HF model are calibrated properly. Feedforward ne… ▽ More Aircraft industry is constantly striving for more efficient design optimization methods in terms of human efforts, computation time, and resource consumption. Hybrid surrogate optimization maintains high results quality while providing rapid design assessments when both the surrogate model and the switch mechanism for eventually transitioning to the HF model are calibrated properly. Feedforward neural networks (FNNs) can capture highly nonlinear input-output mappings, yielding efficient surrogates for aircraft performance factors. However, FNNs often fail to generalize over the out-of-distribution (OOD) samples, which hinders their adoption in critical aircraft design optimization. Through SmOOD, our smoothness-based out-of-distribution detection approach, we propose to codesign a model-dependent OOD indicator with the optimized FNN surrogate, to produce a trustworthy surrogate model with selective but credible predictions. Unlike conventional uncertainty-grounded methods, SmOOD exploits inherent smoothness properties of the HF simulations to effectively expose OODs through revealing their suspicious sensitivities, thereby avoiding over-confident uncertainty estimates on OOD samples. By using SmOOD, only high-risk OOD inputs are forwarded to the HF model for re-evaluation, leading to more accurate results at a low overhead cost. Three aircraft performance models are investigated. Results show that FNN-based surrogates outperform their Gaussian Process counterparts in terms of predictive performance. Moreover, SmOOD does cover averagely 85% of actual OODs on all the study cases. When SmOOD plus FNN surrogates are deployed in hybrid surrogate optimization settings, they result in a decrease error rate of 34.65% and a computational speed up rate of 58.36 times, respectively. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Physics-Guided Adversarial Machine Learning for Aircraft Systems Simulation

Authors: Houssem Ben Braiek, Thomas Reid, Foutse Khomh

Abstract: In the context of aircraft system performance assessment, deep learning technologies allow to quickly infer models from experimental measurements, with less detailed system knowledge than usually required by physics-based modeling. However, this inexpensive model development also comes with new challenges regarding model trustworthiness. This work presents a novel approach, physics-guided adversar… ▽ More In the context of aircraft system performance assessment, deep learning technologies allow to quickly infer models from experimental measurements, with less detailed system knowledge than usually required by physics-based modeling. However, this inexpensive model development also comes with new challenges regarding model trustworthiness. This work presents a novel approach, physics-guided adversarial machine learning (ML), that improves the confidence over the physics consistency of the model. The approach performs, first, a physics-guided adversarial testing phase to search for test inputs revealing behavioral system inconsistencies, while still falling within the range of foreseeable operational conditions. Then, it proceeds with physics-informed adversarial training to teach the model the system-related physics domain foreknowledge through iteratively reducing the unwanted output deviations on the previously-uncovered counterexamples. Empirical evaluation on two aircraft system performance models shows the effectiveness of our adversarial ML approach in exposing physical inconsistencies of both models and in improving their propensity to be consistent with physics domain knowledge. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Study on Electropolishing Conditions for 650 MHz Niobium SRF Cavity

Authors: V. Chouhan, D. Bice, F. Furuta, M. Martinello, M. K. Ng, H. Park, T. Ring, G. Wu, B. Guilfoyle, M. P. Kelly, T. Reid

Abstract: The PIP II linear accelerator includes different types of niobium SRF cavities including 650 MHz elliptical low (0.61) and high (0.92) beta cavities. The elliptical cavity surface is processed with the electropolishing method. The elliptical cavities especially the low-$βべーた$ 650 MHz cavities showed a rough equator surface after the EP was per-formed with the standard EP conditions. This work was foc… ▽ More The PIP II linear accelerator includes different types of niobium SRF cavities including 650 MHz elliptical low (0.61) and high (0.92) beta cavities. The elliptical cavity surface is processed with the electropolishing method. The elliptical cavities especially the low-$βべーた$ 650 MHz cavities showed a rough equator surface after the EP was per-formed with the standard EP conditions. This work was focused to study the effect of different EP parameters, including cathode surface area, temperature and voltage, and optimize them to improve the cavity surface. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Report number: FERMILAB-CONF-22-554-TD

Statistical Properties of the Probabilistic Numeric Linear Solver BayesCG

Authors: Tim W. Reid, Ilse C. F. Ipsen, Jon Cockayne, Chris J. Oates

Abstract: We analyse the calibration of BayesCG under the Krylov prior, a probabilistic numeric extension of the Conjugate Gradient (CG) method for solving systems of linear equations with symmetric positive definite coefficient matrix. Calibration refers to the statistical quality of the posterior covariances produced by a solver. Since BayesCG is not calibrated in the strict existing notion, we propose in… ▽ More We analyse the calibration of BayesCG under the Krylov prior, a probabilistic numeric extension of the Conjugate Gradient (CG) method for solving systems of linear equations with symmetric positive definite coefficient matrix. Calibration refers to the statistical quality of the posterior covariances produced by a solver. Since BayesCG is not calibrated in the strict existing notion, we propose instead two test statistics that are necessary but not sufficient for calibration: the Z-statistic and the new S-statistic. We show analytically and experimentally that under low-rank approximate Krylov posteriors, BayesCG exhibits desirable properties of a calibrated solver, is only slightly optimistic, and is computationally competitive with CG. △ Less

Submitted 7 August, 2022; originally announced August 2022.

Comments: 40 Pages

MSC Class: 65F10; 62F15; 15A06

Line Segment Tracking in the HL-LHC

Authors: Gavin Niendorf, Tres Reid, Peter Wittich, Peter Elmer, Bei Wang, Philip Chang, Yanxi Gu, Vyacheslav Krutelyov, Balaji Venkat Sathia Narayanan, Matevz Tadel, Emmanouil Vourliotis, Avi Yagil

Abstract: The major challenge posed by the high instantaneous luminosity in the High Luminosity LHC (HL-LHC) motivates efficient and fast reconstruction of charged particle tracks in a high pile-up environment. While there have been efforts to use modern techniques like vectorization to improve the existing classic Kalman Filter based reconstruction algorithms, Line Segment Tracking takes a fundamentally di… ▽ More The major challenge posed by the high instantaneous luminosity in the High Luminosity LHC (HL-LHC) motivates efficient and fast reconstruction of charged particle tracks in a high pile-up environment. While there have been efforts to use modern techniques like vectorization to improve the existing classic Kalman Filter based reconstruction algorithms, Line Segment Tracking takes a fundamentally different approach by doing a bottom-up reconstruction of tracks. Small track stubs from adjoining detector regions are constructed, and then these track stubs that are consistent with typical track trajectories are successively linked. Since the production of these track stubs is localized, they can be made in parallel, which lends way into using architectures like GPUs and multi-CPUs to take advantage of the parallelism. The algorithm is implemented in the context of the CMS Phase-2 Tracker and runs on NVIDIA Tesla V100 GPUs. Good physics and timing performance has been obtained, and stepping stones for the future are elaborated. △ Less

Submitted 28 September, 2022; v1 submitted 17 July, 2022; originally announced July 2022.

Comments: Connecting the Dots Workshop (CTD 2022) v2 : Implemented referee comments from internal review

Q-factor optimization for high-beta 650 MHz cavities for PIP-II

Authors: M. Martinello, D. J. Bice, C. Boffo, S. K. Chandrasekeran, G. V. Eremeev, F. Furuta, A. Grassellino, O. Melnychuk, D. A. Sergatskov, G. Wu, T. C. Reid

Abstract: High Q-factors are of utmost importance to minimize losses of superconducting radio-frequency cavities deployed in continuous wave particle accelerators. This study elucidates the surface treatment that can maximize the Q-factors in high-beta 650 MHz elliptical niobium cavities. State-of-the-art surface treatments are applied in many single-cell cavities, and surface resistance studies are perform… ▽ More High Q-factors are of utmost importance to minimize losses of superconducting radio-frequency cavities deployed in continuous wave particle accelerators. This study elucidates the surface treatment that can maximize the Q-factors in high-beta 650 MHz elliptical niobium cavities. State-of-the-art surface treatments are applied in many single-cell cavities, and surface resistance studies are performed to understand the microwave dissipation at this unexplored frequency. The nitrogen doping treatment is confirmed to be necessary to maximize the Q-factors at medium RF fields. We applied this treatment in five-cell high-beta 650 MHz cavities and demonstrated that extremely high Q-factors were obtained at medium RF fields with this treatment. We also demonstrated that adding a cold electropolishing step after N-doping is crucial to push the quench field of multicell cavities to higher gradients. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Journal ref: J. Appl. Phys. 130, 174501 (2021)

Randomized algorithms for rounding in the Tensor-Train format

Authors: Hussam Al Daas, Grey Ballard, Paul Cazeaux, Eric Hallman, Agnieszka Miedlar, Mirjeta Pasha, Tim W. Reid, Arvind K. Saibaba

Abstract: The Tensor-Train (TT) format is a highly compact low-rank representation for high-dimensional tensors. TT is particularly useful when representing approximations to the solutions of certain types of parametrized partial differential equations. For many of these problems, computing the solution explicitly would require an infeasible amount of memory and computational time. While the TT format makes… ▽ More The Tensor-Train (TT) format is a highly compact low-rank representation for high-dimensional tensors. TT is particularly useful when representing approximations to the solutions of certain types of parametrized partial differential equations. For many of these problems, computing the solution explicitly would require an infeasible amount of memory and computational time. While the TT format makes these problems tractable, iterative techniques for solving the PDEs must be adapted to perform arithmetic while maintaining the implicit structure. The fundamental operation used to maintain feasible memory and computational time is called rounding, which truncates the internal ranks of a tensor already in TT format. We propose several randomized algorithms for this task that are generalizations of randomized low-rank matrix approximation algorithms and provide significant reduction in computation compared to deterministic TT-rounding algorithms. Randomization is particularly effective in the case of rounding a sum of TT-tensors (where we observe 20x speedup), which is the bottleneck computation in the adaptation of GMRES to vectors in TT format. We present the randomized algorithms and compare their empirical accuracy and computational time with deterministic alternatives. △ Less

Submitted 8 October, 2021; originally announced October 2021.

GNSS Radio Occultation on Aerial Platforms with Commercial Off-The-Shelf Receivers

Authors: Bryan C. Chan, Ashish Goel, Jonathan Kosh, Tyler G. R. Reid, Corey R. Snyder, Paul M. Tarantino, Saraswati Soedarmadji, Widyadewi Soedarmadji, Kevin Nelson, Feiqin Xie, Michael Vergalla

Abstract: In recent decades, GNSS Radio Occultation soundings have proven an invaluable input to global weather forecasting. The success of government-sponsored programs such as COSMIC is now complemented by commercial low-cost cubesat implementations. The result is access to more than 10,000 soundings per day and improved weather forecasting accuracy. This movement towards commercialization has been suppor… ▽ More In recent decades, GNSS Radio Occultation soundings have proven an invaluable input to global weather forecasting. The success of government-sponsored programs such as COSMIC is now complemented by commercial low-cost cubesat implementations. The result is access to more than 10,000 soundings per day and improved weather forecasting accuracy. This movement towards commercialization has been supported by several agencies, including the National Aeronautics and Space Administration (NASA), National Oceanic and Atmospheric Administration (NOAA) and the U.S. Air Force (USAF) with programs such as the Commercial Weather Data Pilot (CWDP). This has resulted in further interest in commercially deploying GNSS-RO on complementary platforms. Here, we examine a so far underutilized platform: the high-altitude weather balloon. Such meteorological radiosondes are deployed twice daily at over 900 locations globally and form an essential in-situ data source as a long-standing input to weather forecasting models. Adding GNSS-RO capability to existing radiosonde platforms would greatly expand capability, allowing for persistent and local area monitoring, a feature particularly useful for hurricane and other severe weather monitoring. A prohibitive barrier to entry to this inclusion is cost and complexity as GNSS-RO traditionally requires highly specialized and sensitive equipment. This paper describes a multi-year effort to develop a low-cost and scalable approach to balloon GNSS-RO based on Commercial-Off-The-Shelf (COTS) GNSS receivers. We present hardware prototypes and data processing techniques which demonstrate the technical feasibility of the approach through results from several flight testing campaigns. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 16 pages, 25 figures, ION GNSS+ 2021

Journal ref: Proceedings of the 34th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2021)

Two-Photon Dual-Comb LiDAR

Authors: Hollie Wright, Jinghua Sun, David McKendrick, Nick Weston, Derryck T. Reid

Abstract: The interferometric signals produced in conventional dual-comb laser ranging require femtosecond lasers with long-term f_CEO stability, and are limited to an upper sampling rate by radio-frequency aliasing considerations. By using cross-polarized dual combs and two-photon detection we demonstrate carrier-phase-insensitive cross-correlations at sampling rates of up to 12x the conventional dual-comb… ▽ More The interferometric signals produced in conventional dual-comb laser ranging require femtosecond lasers with long-term f_CEO stability, and are limited to an upper sampling rate by radio-frequency aliasing considerations. By using cross-polarized dual combs and two-photon detection we demonstrate carrier-phase-insensitive cross-correlations at sampling rates of up to 12x the conventional dual-comb aliasing limit, recording these in a digitizer-based acquisition system to implement ranging with sub-100-nm precision. We then extend this concept to show how the high data burden of conventional dual-comb acquisition can be eliminated by using a simple microcontroller as a ns-precision stopwatch to record the time intervals separating the two-photon cross-correlation pulses, providing real-time and continuous LiDAR-like distance metrology capable of sub-100 nm precision and dynamic acquisition for unlimited periods. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 9 pages, 6 figures

Study of energy response and resolution of the ATLAS Tile Calorimeter to hadrons of energies from 16 to 30 GeV

Authors: Jalal Abdallah, Stylianos Angelidakis, Giorgi Arabidze, Nikolay Atanov, Johannes Bernhard, Romeo Bonnefoy, Jonathan Bossio, Ryan Bouabid, Fernando Carrio, Tomas Davidek, Michal Dubovsky, Luca Fiorini, Francisco Brandan Garcia Aparisi, Tancredi Carli, Alexander Gerbershagen, Hazal Goksu, Haleh Hadavand, Siarhei Harkusha, Dingane Hlaluku, Michael James Hibbard, Kevin Hildebrand, Juansher Jejelava, Andrey Kamenshchikov, Stergios Kazakos, Tomas Kello , et al. (46 additional authors not shown)

Abstract: Three spare modules of the ATLAS Tile Calorimeter were exposed to test beams from the Super Proton Synchrotron accelerator at CERN in 2017. The measurements of the energy response and resolution of the detector to positive pions and kaons and protons with energy in the range 16 to 30 GeV are reported. The results have uncertainties of few percent. They were compared to the predictions of the Geant… ▽ More Three spare modules of the ATLAS Tile Calorimeter were exposed to test beams from the Super Proton Synchrotron accelerator at CERN in 2017. The measurements of the energy response and resolution of the detector to positive pions and kaons and protons with energy in the range 16 to 30 GeV are reported. The results have uncertainties of few percent. They were compared to the predictions of the Geant4-based simulation program used in ATLAS to estimate the response of the detector to proton-proton events at Large Hadron Collider. The determinations obtained using experimental and simulated data agree within the uncertainties. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Rooted prism-minors and disjoint cycles containing a specified edge

Authors: João Paulo Costalonga, Talmage James Reid, Haindong Wu

Abstract: Dirac and Lovász independently characterized the $3$-connected graphs with no pair of vertex-disjoint cycles. Equivalently, they characterized all $3$-connected graphs with no prism-minors. In this paper, we completely characterize the $3$-connected graphs with an edge that is contained in the union of no pair of vertex-disjoint cycles. As applications, we answer the analogous questions for edge-d… ▽ More Dirac and Lovász independently characterized the $3$-connected graphs with no pair of vertex-disjoint cycles. Equivalently, they characterized all $3$-connected graphs with no prism-minors. In this paper, we completely characterize the $3$-connected graphs with an edge that is contained in the union of no pair of vertex-disjoint cycles. As applications, we answer the analogous questions for edge-disjoint cycles and for $4$-connected graphs and we completely characterize the $3$-connected graphs with no prism-minor using a specified edge. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Probabilistic Iterative Methods for Linear Systems

Authors: Jon Cockayne, Ilse C. F. Ipsen, Chris J. Oates, Tim W. Reid

Abstract: This paper presents a probabilistic perspective on iterative methods for approximating the solution $\mathbf{x}_* \in \mathbb{R}^d$ of a nonsingular linear system $\mathbf{A} \mathbf{x}_* = \mathbf{b}$. In the approach a standard iterative method on $\mathbb{R}^d$ is lifted to act on the space of probability distributions $\mathcal{P}(\mathbb{R}^d)$. Classically, an iterative method produces a seq… ▽ More This paper presents a probabilistic perspective on iterative methods for approximating the solution $\mathbf{x}_* \in \mathbb{R}^d$ of a nonsingular linear system $\mathbf{A} \mathbf{x}_* = \mathbf{b}$. In the approach a standard iterative method on $\mathbb{R}^d$ is lifted to act on the space of probability distributions $\mathcal{P}(\mathbb{R}^d)$. Classically, an iterative method produces a sequence $\mathbf{x}_m$ of approximations that converge to $\mathbf{x}_*$. The output of the iterative methods proposed in this paper is, instead, a sequence of probability distributions $μみゅー_m \in \mathcal{P}(\mathbb{R}^d)$. The distributional output both provides a "best guess" for $\mathbf{x}_*$, for example as the mean of $μみゅー_m$, and also probabilistic uncertainty quantification for the value of $\mathbf{x}_*$ when it has not been exactly determined. Theoretical analysis is provided in the prototypical case of a stationary linear iterative method. In this setting we characterise both the rate of contraction of $μみゅー_m$ to an atomic measure on $\mathbf{x}_*$ and the nature of the uncertainty quantification being provided. We conclude with an empirical illustration that highlights the insight into solution uncertainty that can be provided by probabilistic iterative methods. △ Less

Submitted 11 January, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Ford Highway Driving RTK Dataset: 30,000 km of North American Highways

Authors: Sarah E. Houts, Nahid Pervez, Umair Ibrahim, Gaurav Pandey, Tyler G. R. Reid

Abstract: There is a growing need for vehicle positioning information to support Advanced Driver Assistance Systems (ADAS), Connectivity (V2X), and Autonomous Driving (AD) features. These range from a need for road determination ($<$5 meters), lane determination ($<$1.5 meters), and determining where the vehicle is within the lane ($<$0.3 meters). This paper presents the Ford Highway Driving RTK (Ford-HDR)… ▽ More There is a growing need for vehicle positioning information to support Advanced Driver Assistance Systems (ADAS), Connectivity (V2X), and Autonomous Driving (AD) features. These range from a need for road determination ($<$5 meters), lane determination ($<$1.5 meters), and determining where the vehicle is within the lane ($<$0.3 meters). This paper presents the Ford Highway Driving RTK (Ford-HDR) dataset. This dataset includes nearly 30,000 km of data collected primarily on North American highways during a driving campaign designed to validate driver assistance features in 2018. This includes data from a representative automotive production GNSS used primarily for turn-by-turn navigation as well as an Inertial Navigation System (INS) which couples two survey-grade GNSS receivers with a tactical grade Inertial Measurement Unit (IMU) to act as ground truth. The latter utilized networked Real-Time Kinematic (RTK) GNSS corrections delivered over a cellular modem in real-time. This dataset is being released into the public domain to spark further research in the community. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: 8 pages, 4 figures, ION GNSS+ 2020

Journal ref: Proceedings of the 33rd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2020)

BayesCG As An Uncertainty Aware Version of CG

Authors: Tim W. Reid, Ilse C. F. Ipsen, Jon Cockayne, Chris J. Oates

Abstract: The Bayesian Conjugate Gradient method (BayesCG) is a probabilistic generalization of the Conjugate Gradient method (CG) for solving linear systems with real symmetric positive definite coefficient matrices. Our CG-based implementation of BayesCG under a structure-exploiting prior distribution represents an 'uncertainty-aware' version of CG. Its output consists of CG iterates and posterior covaria… ▽ More The Bayesian Conjugate Gradient method (BayesCG) is a probabilistic generalization of the Conjugate Gradient method (CG) for solving linear systems with real symmetric positive definite coefficient matrices. Our CG-based implementation of BayesCG under a structure-exploiting prior distribution represents an 'uncertainty-aware' version of CG. Its output consists of CG iterates and posterior covariances that can be propagated to subsequent computations. The covariances have low-rank and are maintained in factored form. This allows easy generation of accurate samples to probe uncertainty in downstream computations. Numerical experiments confirm the effectiveness of the low-rank posterior covariances. △ Less

Submitted 3 October, 2022; v1 submitted 7 August, 2020; originally announced August 2020.

Comments: 34 Pages including supplementary material (main paper is 23 pages, supplement is 11 pages). Computer codes are available at https://github.com/treid5/ProbNumCG_Supp

MSC Class: 65F10; 62F15; 65F50; 15A06; 15A10

Human Trust-based Feedback Control: Dynamically varying automation transparency to optimize human-machine interactions

Authors: Kumar Akash, Griffon McMahon, Tahira Reid, Neera Jain

Abstract: Human trust in automation plays an essential role in interactions between humans and automation. While a lack of trust can lead to a human's disuse of automation, over-trust can result in a human trusting a faulty autonomous system which could have negative consequences for the human. Therefore, human trust should be calibrated to optimize human-machine interactions with respect to context-specifi… ▽ More Human trust in automation plays an essential role in interactions between humans and automation. While a lack of trust can lead to a human's disuse of automation, over-trust can result in a human trusting a faulty autonomous system which could have negative consequences for the human. Therefore, human trust should be calibrated to optimize human-machine interactions with respect to context-specific performance objectives. In this article, we present a probabilistic framework to model and calibrate a human's trust and workload dynamics during his/her interaction with an intelligent decision-aid system. This calibration is achieved by varying the automation's transparency---the amount and utility of information provided to the human. The parameterization of the model is conducted using behavioral data collected through human-subject experiments, and three feedback control policies are experimentally validated and compared against a non-adaptive decision-aid system. The results show that human-automation team performance can be optimized when the transparency is dynamically updated based on the proposed control policy. This framework is a first step toward widespread design and implementation of real-time adaptive automation for use in human-machine interactions. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: 21 pages

ACM Class: H.1.2; H.5

Journal ref: IEEE Control Systems Magazine, vol. 40, no. 6, pp. 98-116, Dec. 2020

Satellite Navigation for the Age of Autonomy

Authors: Tyler G. R. Reid, Bryan Chan, Ashish Goel, Kazuma Gunning, Brian Manning, Jerami Martin, Andrew Neish, Adrien Perkins, Paul Tarantino

Abstract: Global Navigation Satellite Systems (GNSS) brought navigation to the masses. Coupled with smartphones, the blue dot in the palm of our hands has forever changed the way we interact with the world. Looking forward, cyber-physical systems such as self-driving cars and aerial mobility are pushing the limits of what localization technologies including GNSS can provide. This autonomous revolution requi… ▽ More Global Navigation Satellite Systems (GNSS) brought navigation to the masses. Coupled with smartphones, the blue dot in the palm of our hands has forever changed the way we interact with the world. Looking forward, cyber-physical systems such as self-driving cars and aerial mobility are pushing the limits of what localization technologies including GNSS can provide. This autonomous revolution requires a solution that supports safety-critical operation, centimeter positioning, and cyber-security for millions of users. To meet these demands, we propose a navigation service from Low Earth Orbiting (LEO) satellites which deliver precision in-part through faster motion, higher power signals for added robustness to interference, constellation autonomous integrity monitoring for integrity, and encryption / authentication for resistance to spoofing attacks. This paradigm is enabled by the 'New Space' movement, where highly capable satellites and components are now built on assembly lines and launch costs have decreased by more than tenfold. Such a ubiquitous positioning service enables a consistent and secure standard where trustworthy information can be validated and shared, extending the electronic horizon from sensor line of sight to an entire city. This enables the situational awareness needed for true safe operation to support autonomy at scale. △ Less

Submitted 18 May, 2020; originally announced May 2020.

Comments: 11 pages, 8 figures, 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS)

Journal ref: IEEE/ION PLANS, Portland, OR, pp. 342-352, 2020

Factorized Machine Learning for Performance Modeling of Massively Parallel Heterogeneous Physical Simulations

Authors: Ardavan Oskooi, Christopher Hogan, Alec M. Hammond, M. T. Homer Reid, Steven G. Johnson

Abstract: We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part… ▽ More We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part strategy. First, we perform data-driven static load balancing using regression coefficients extracted from small simulations, which both improves parallel performance and reduces the dependency of the runtime on the precise spatial layout of the heterogeneous physics. Second, we divide the execution time of these load-balanced simulations into computation and communication, factoring crude asymptotic scalings out of each term, and training neural nets for the remaining factor coefficients. This strategy is implemented for Meep, a popular and complex open-source electrodynamics simulation package, and are validated for heterogeneous simulations drawn from published engineering models. △ Less

Submitted 6 October, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 19 pages, 8 figures

Developments in Modern GNSS and Its Impact on Autonomous Vehicle Architectures

Authors: Niels Joubert, Tyler GR Reid, Fergus Noble

Abstract: This paper surveys a number of recent developments in modern Global Navigation Satellite Systems (GNSS) and investigates the possible impact on autonomous driving architectures. Modern GNSS now consist of four independent global satellite constellations delivering modernized signals at multiple civil frequencies. New ground monitoring infrastructure, mathematical models, and internet services corr… ▽ More This paper surveys a number of recent developments in modern Global Navigation Satellite Systems (GNSS) and investigates the possible impact on autonomous driving architectures. Modern GNSS now consist of four independent global satellite constellations delivering modernized signals at multiple civil frequencies. New ground monitoring infrastructure, mathematical models, and internet services correct for errors in the GNSS signals at continent scale. Mass-market automotive-grade receiver chipsets are available at low Cost, Size, Weight, and Power (CSWaP). The result is that GNSS in 2020 delivers better than lane-level accurate localization with 99.99999% integrity guarantees at over 95% availability. In autonomous driving, SAE Level 2 partially autonomous vehicles are now available to consumers, capable of autonomously following lanes and performing basic maneuvers under human supervision. Furthermore, the first pilot programs of SAE Level 4 driverless vehicles are being demonstrated on public roads. However, autonomous driving is not a solved problem. GNSS can help. Specifically, incorporating high-integrity GNSS lane determination into vision-based architectures can unlock lane-level maneuvers and provide oversight to guarantee safety. Incorporating precision GNSS into LiDAR-based systems can unlock robustness and additional fallbacks for safety and utility. Lastly, GNSS provides interoperability through consistent timing and reference frames for future V2X scenarios. △ Less

Submitted 19 April, 2020; v1 submitted 2 February, 2020; originally announced February 2020.

Comments: 2020 IEEE Intelligent Vehicles Symposium. IEEE, 2020

Report number: its.IV2020.412.8d033de9

First Principles Modeling of Topological Insulators: Structural Optimization and Exchange Correlation Functionals

Authors: Thomas K. Reid, S. Pamir Alpay, Alexander V. Balatsky, Sanjeev K. Nayak

Abstract: Topological insulators (TIs) are materials that are insulating in the bulk but have zero band gap surface states with linear dispersion and are protected by time reversal symmetry. These unique characteristics could pave the way for many promising applications that include spintronic devices and quantum computations. It is important to understand and theoretically describe TIs as accurately as pos… ▽ More Topological insulators (TIs) are materials that are insulating in the bulk but have zero band gap surface states with linear dispersion and are protected by time reversal symmetry. These unique characteristics could pave the way for many promising applications that include spintronic devices and quantum computations. It is important to understand and theoretically describe TIs as accurately as possible in order to predict properties. Quantum mechanical approaches, specifically first principles density functional theory (DFT) based methods, have been used extensively to model electronic properties of TIs. Here, we provide a comprehensive assessment of a variety of DFT formalisms and how these capture the electronic structure of TIs. We concentrate on Bi$_2$Se$_3$ and Bi$_2$Te$_3$ as examples of prototypical TI materials. We find that the generalized gradient (GGA) and kinetic density functional (metaGGA) produce displacements increasing the thickness of the TI slab, whereas we see an opposite behavior in DFT computations using LDA. Accounting for van der Waals (vdW) interactions overcomes the apparent over-relaxations and retraces the atomic positions towards the bulk. Based on an intensive computational study, we show that GGA with vdW treatment is the most appropriate method for structural optimization. Electronic structures derived from GGA or metaGGA employing experimental lattice parameters are also acceptable. In this regard, we express a slight preference for metaGGA in terms of accuracy, but an overall preference for GGA due to compensatory improvements in computability in capturing TI behavior. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Journal ref: Phys. Rev. B 101, 085140 (2020)

Application of the Biot-Savart Law to Parabolic Vortex Segments using Elliptic Integrals

Authors: Andreas Malmendier, Jackson T. Reid

Abstract: The Biot-Savart law is used in aerodynamic theory to calculate the velocity induced by curved vortex lines. Explicit formulas are developed, using multivariate Appell hypergeometric functions, for the velocity induced by a general parabolic vortex segment. The formulas are derived by constructing a particular pencil of elliptic curves whose period integrals provide the solution to the induced velo… ▽ More The Biot-Savart law is used in aerodynamic theory to calculate the velocity induced by curved vortex lines. Explicit formulas are developed, using multivariate Appell hypergeometric functions, for the velocity induced by a general parabolic vortex segment. The formulas are derived by constructing a particular pencil of elliptic curves whose period integrals provide the solution to the induced velocity. We use numerical integration and a perturbation expansion to evaluate the validity of our formulas. △ Less

Submitted 2 July, 2020; v1 submitted 23 July, 2019; originally announced July 2019.

Comments: 18 pages, 5 figures; minor typos corrected in version 2

MSC Class: 76G25; 33E05

Journal ref: Z. Angew. Math. Phys. 71 (2020), no. 4, Paper No. 104, 20 pp

Standalone and RTK GNSS on 30,000 km of North American Highways

Authors: Tyler G. R. Reid, Nahid Pervez, Umair Ibrahim, Sarah E. Houts, Gaurav Pandey, Naveen K. R. Alla, Andy Hsia

Abstract: There is a growing need for vehicle positioning information to support Advanced Driver Assistance Systems (ADAS), Connectivity (V2X), and Automated Driving (AD) features. These range from a need for road determination (<5 meters), lane determination (<1.5 meters), and determining where the vehicle is within the lane (<0.3 meters). This work examines the performance of Global Navigation Satellite S… ▽ More There is a growing need for vehicle positioning information to support Advanced Driver Assistance Systems (ADAS), Connectivity (V2X), and Automated Driving (AD) features. These range from a need for road determination (<5 meters), lane determination (<1.5 meters), and determining where the vehicle is within the lane (<0.3 meters). This work examines the performance of Global Navigation Satellite Systems (GNSS) on 30,000 km of North American highways to better understand the automotive positioning needs it meets today and what might be possible in the near future with wide area GNSS correction services and multi-frequency receivers. This includes data from a representative automotive production GNSS used primarily for turn-by-turn navigation as well as an Inertial Navigation System which couples two survey grade GNSS receivers with a tactical grade Inertial Measurement Unit (IMU) to act as ground truth. The latter utilized networked Real-Time Kinematic (RTK) GNSS corrections delivered over a cellular modem in real-time. We assess on-road GNSS accuracy, availability, and continuity. Availability and continuity are broken down in terms of satellite visibility, satellite geometry, position type (RTK fixed, RTK float, or standard positioning), and RTK correction latency over the network. Results show that current automotive solutions are best suited to meet road determination requirements at 98% availability but are less suitable for lane determination at 57%. Multi-frequency receivers with RTK corrections were found more capable with road determination at 99.5%, lane determination at 98%, and highway-level lane departure protection at 91%. △ Less

Submitted 28 August, 2019; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: Accepted for the 32nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2019), Miami, Florida, September 2019

Journal ref: Proceedings of the 32nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2019), Miami, Florida, September 2019, pp. 2135-2158

Localization Requirements for Autonomous Vehicles

Authors: Tyler G. R. Reid, Sarah E. Houts, Robert Cammarata, Graham Mills, Siddharth Agarwal, Ankit Vora, Gaurav Pandey

Abstract: Autonomous vehicles require precise knowledge of their position and orientation in all weather and traffic conditions for path planning, perception, control, and general safe operation. Here we derive these requirements for autonomous vehicles based on first principles. We begin with the safety integrity level, defining the allowable probability of failure per hour of operation based on desired im… ▽ More Autonomous vehicles require precise knowledge of their position and orientation in all weather and traffic conditions for path planning, perception, control, and general safe operation. Here we derive these requirements for autonomous vehicles based on first principles. We begin with the safety integrity level, defining the allowable probability of failure per hour of operation based on desired improvements on road safety today. This draws comparisons with the localization integrity levels required in aviation and rail where similar numbers are derived at 10^-8 probability of failure per hour of operation. We then define the geometry of the problem, where the aim is to maintain knowledge that the vehicle is within its lane and to determine what road level it is on. Longitudinal, lateral, and vertical localization error bounds (alert limits) and 95% accuracy requirements are derived based on US road geometry standards (lane width, curvature, and vertical clearance) and allowable vehicle dimensions. For passenger vehicles operating on freeway roads, the result is a required lateral error bound of 0.57 m (0.20 m, 95%), a longitudinal bound of 1.40 m (0.48 m, 95%), a vertical bound of 1.30 m (0.43 m, 95%), and an attitude bound in each direction of 1.50 deg (0.51 deg, 95%). On local streets, the road geometry makes requirements more stringent where lateral and longitudinal error bounds of 0.29 m (0.10 m, 95%) are needed with an orientation requirement of 0.50 deg (0.17 deg, 95%). △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: Under review with the SAE Journal of Connected and Automated Vehicles

Journal ref: SAE Intl. J CAV 2(3):2019

A Classification Model for Sensing Human Trust in Machines Using EEG and GSR

Authors: Kumar Akash, Wan-Lin Hu, Neera Jain, Tahira Reid

Abstract: Today, intelligent machines \emph{interact and collaborate} with humans in a way that demands a greater level of trust between human and machine. A first step towards building intelligent machines that are capable of building and maintaining trust with humans is the design of a sensor that will enable machines to estimate human trust level in real-time. In this paper, two approaches for developing… ▽ More Today, intelligent machines \emph{interact and collaborate} with humans in a way that demands a greater level of trust between human and machine. A first step towards building intelligent machines that are capable of building and maintaining trust with humans is the design of a sensor that will enable machines to estimate human trust level in real-time. In this paper, two approaches for developing classifier-based empirical trust sensor models are presented that specifically use electroencephalography (EEG) and galvanic skin response (GSR) measurements. Human subject data collected from 45 participants is used for feature extraction, feature selection, classifier training, and model validation. The first approach considers a general set of psychophysiological features across all participants as the input variables and trains a classifier-based model for each participant, resulting in a trust sensor model based on the general feature set (i.e., a "general trust sensor model"). The second approach considers a customized feature set for each individual and trains a classifier-based model using that feature set, resulting in improved mean accuracy but at the expense of an increase in training time. This work represents the first use of real-time psychophysiological measurements for the development of a human trust sensor. Implications of the work, in the context of trust management algorithm design for intelligent machines, are also discussed. △ Less

Submitted 26 March, 2018; originally announced March 2018.

Comments: 20 pages

ACM Class: H.5; H.1.2

Journal ref: ACM Trans. Interact. Intell. Syst. 8, 4, Article 27 (December 2018)

Taylor-Duffy Method for Singular Tetrahedron-Product Integrals: Efficient Evaluation of Galerkin Integrals for VIE Solvers

Authors: M. T. Homer Reid

Abstract: I present an accurate and efficient technique for numerical evaluation of singular 6-dimensional integrals over tetrahedon-product domains, with applications to calculation of Galerkin matrix elements for discretized volume-integral-equation (VIE) solvers using Schaubert-Wilton-Glisson (SWG) and other tetrahedral basis functions. My method extends the generalized Taylor-Duffy strategy---used to ha… ▽ More I present an accurate and efficient technique for numerical evaluation of singular 6-dimensional integrals over tetrahedon-product domains, with applications to calculation of Galerkin matrix elements for discretized volume-integral-equation (VIE) solvers using Schaubert-Wilton-Glisson (SWG) and other tetrahedral basis functions. My method extends the generalized Taylor-Duffy strategy---used to handle the singular \textit{triangle}-product integrals arising in discretized surface-integral-equation (SIE) formulations---to the tetrahedron-product case; it effects an exact transformation of a singular 6-dimensional integral to an nonsingular lower-dimensional integral that may be evaluated by simple numerical cubature The method is highly general and may---with the aid of automatic code generation facilitated by computer-algebra systems---be applied to a wide variety of singular integrals arising in various VIE formulations with various types of tetrahedral basis function, of which I present several examples. To demonstrate the accuracy and efficiency of my method, I apply it to the calculation of matrix elements for the volume electric-field integral equation (VEFIE) discretized with SWG basis functions, where the method yields 12-digit or higher accuracy with low computational cost---an improvement of many orders of magnitude compared to existing techniques. △ Less

Submitted 7 August, 2017; originally announced August 2017.

Comments: 9 pages, 2 figures

Photon Torpedoes and Rytov Pinwheels: Integral-Equation Modeling of Non-Equilibrium Fluctuation-Induced Forces and Torques on Nanoparticles

Authors: M. T. Homer Reid, O. D. Miller, A. G. Polimeridis, A. W. Rodriguez, E. M. Tomlinson, S. G. Johnson

Abstract: We present new theoretical tools, based on fluctuational electrodynamics and the integral-equation approach to computational electromagnetism, for numerical modeling of forces and torques on bodies of complex shapes and materials due to emission of thermal radiation out of thermal equilibrium. This extends our recently-developed fluctuating-surface-current (FSC) and fluctuating-volume-current (FVC… ▽ More We present new theoretical tools, based on fluctuational electrodynamics and the integral-equation approach to computational electromagnetism, for numerical modeling of forces and torques on bodies of complex shapes and materials due to emission of thermal radiation out of thermal equilibrium. This extends our recently-developed fluctuating-surface-current (FSC) and fluctuating-volume-current (FVC) techniques for radiative heat transfer to the computation of non-equilibrium fluctuation-induced forces and torques; as we show, the extension is non-trivial due to the greater computational cost of modeling radiative momentum transfer, including new singularities that must be carefully neutralized. We introduce a new analytical cancellation technique that addresses these challenges and allows, for the first time, accurate and efficient prediction of non-equilibrium forces and torques on bodies of essentially arbitrary shapes---including asymmetric and chiral particles---and complex material properties, including continuously-varying and anisotropic dielectrics. We validate our approach by showing that it reproduces known results, then present new numerical predictions of non-equilibrium self-propulsion, self-rotation, and momentum-transfer phenomena in complex geometries that would be difficult or impossible to study with existing methods. Our findings indicate that the fluctuation-induced dynamics of micron-size room-temperature bodies in cold environments involve microscopic length scales but macroscopic time scales, with typical linear and angular velocities on the order of microns/second and radians/second; For a micron-scale gear driven by thermal radiation from a nearby chiral emitter, we find a strong and non-monotonic dependence of the magnitude and even the \textit{sign} of the induced torque on the temperature of the emitter. △ Less

Submitted 7 August, 2017; originally announced August 2017.

Comments: 18 pages, 9 figures

Broadband Dual-Comb Spectroscopy in the Spectral Fingerprint Region

Authors: O. Kara, L. Maidment, T. Gardiner, P. G. Schunemann, D. T. Reid

Abstract: Infrared spectroscopy in the spectral fingerprint region from 6-12 um accesses the largest molecular absorption cross-sections, permitting sensitive, quantitative and species-specific measurements. Here, we show how dual-comb spectroscopy--a form of high-speed Fourier-transform spectroscopy involving no moving parts and capable of very high resolutions--can be extended to the 6-8um wavelength band… ▽ More Infrared spectroscopy in the spectral fingerprint region from 6-12 um accesses the largest molecular absorption cross-sections, permitting sensitive, quantitative and species-specific measurements. Here, we show how dual-comb spectroscopy--a form of high-speed Fourier-transform spectroscopy involving no moving parts and capable of very high resolutions--can be extended to the 6-8um wavelength band using femtosecond optical parametric oscillators (OPOs). By acquiring dual-comb interferograms faster than the mutual decoherence time of the OPO combs we implement line-shape-preserving averaging to obtain low-noise, high-fidelity spectra of H2O and CH4 at approximately 0.3 cm-1 resolutions from 1285-1585 cm-1. △ Less

Submitted 26 June, 2017; originally announced June 2017.

Journal ref: Opt. Express 25, 32713-32721 (2017)

Limits to the Optical Response of Graphene and 2D Materials

Authors: Owen D. Miller, Ognjen Ilic, Thomas Christensen, M. T. Homer Reid, Harry A. Atwater, John D. Joannopoulos, Marin Soljacic, Steven G. Johnson

Abstract: 2D materials provide a platform for strong light--matter interactions, creating wide-ranging design opportunities via new-material discoveries and new methods for geometrical structuring. We derive general upper bounds to the strength of such light--matter interactions, given only the optical conductivity of the material, including spatial nonlocality, and otherwise independent of shape and config… ▽ More 2D materials provide a platform for strong light--matter interactions, creating wide-ranging design opportunities via new-material discoveries and new methods for geometrical structuring. We derive general upper bounds to the strength of such light--matter interactions, given only the optical conductivity of the material, including spatial nonlocality, and otherwise independent of shape and configuration. Our material figure of merit shows that highly doped graphene is an optimal material at infrared frequencies, whereas single-atomic-layer silver is optimal in the visible. For quantities ranging from absorption and scattering to near-field spontaneous-emission enhancements and radiative heat transfer, we consider canonical geometrical structures and show that in certain cases the bounds can be approached, while in others there may be significant opportunity for design improvement. The bounds can encourage systematic improvements in the design of ultrathin broadband absorbers, 2D antennas, and near-field energy harvesters. △ Less

Submitted 2 August, 2017; v1 submitted 9 May, 2017; originally announced May 2017.

Comments: 7 pages, 4 figures, plus references and supplementary material

Pure down-conversion photons through sub-coherence length domain engineering

Authors: Francesco Graffitti, Dmytro Kundys, Derryck T. Reid, Agata M. Brańczyk, Alessandro Fedrizzi

Abstract: Photonic quantum technology relies on efficient sources of coherent single photons, the ideal carriers of quantum information. Heralded single photons from parametric down-conversion can approximate on-demand single photons to a desired degree, with high spectral purities achieved through group-velocity matching and tailored crystal nonlinearities. Here we propose crystal nonlinearity engineerin… ▽ More Photonic quantum technology relies on efficient sources of coherent single photons, the ideal carriers of quantum information. Heralded single photons from parametric down-conversion can approximate on-demand single photons to a desired degree, with high spectral purities achieved through group-velocity matching and tailored crystal nonlinearities. Here we propose crystal nonlinearity engineering techniques with sub-coherence-length domains. We first introduce a combination of two existing methods: a deterministic approach with coherence-length domains and probabilistic domain-width annealing. We then show how the same deterministic domain-flip approach can be implemented with sub-coherence length domains. Both of these complementary techniques create highly pure photons, outperforming previous methods, in particular for short nonlinear crystals matched to femtosecond lasers. △ Less

Submitted 20 April, 2017; v1 submitted 12 April, 2017; originally announced April 2017.

Comments: 12 pages, 4 figures. Minor update to Fig. 3

Journal ref: Quantum Sci. Technol. 2, 035001 (2017)

Non-intuitive Computational Optimization of Illumination Patterns for Maximum Optical Force and Torque

Authors: Yoonkyung E. Lee, Owen D. Miller, M. T. Homer Reid, Steven G. Johnson, Nicholas X. Fang

Abstract: This paper aims to maximize optical force and torque on arbitrary micro- and nano-scale objects using numerically optimized structured illumination. By developing a numerical framework for computer-automated design of 3d vector-field illumination, we demonstrate a 20-fold enhancement in optical torque per intensity over circularly polarized plane wave on a model plasmonic particle. The nonconvex o… ▽ More This paper aims to maximize optical force and torque on arbitrary micro- and nano-scale objects using numerically optimized structured illumination. By developing a numerical framework for computer-automated design of 3d vector-field illumination, we demonstrate a 20-fold enhancement in optical torque per intensity over circularly polarized plane wave on a model plasmonic particle. The nonconvex optimization is efficiently performed by combining a compact cylindrical Bessel basis representation with a fast boundary element method and a standard derivative-free, local optimization algorithm. We analyze the optimization results for 2000 random initial configurations, discuss the tradeoff between robustness and enhancement, and compare the different effects of multipolar plasmon resonances on enhancing force and torque. All results are obtained using open-source computational software available online. △ Less

Submitted 26 January, 2017; originally announced January 2017.

Comments: 12 pages, 6 figures

EELT-HIRES the high-resolution spectrograph for the E-ELT

Authors: A. Marconi, P. Di Marcantonio, V. D'Odorico, S. Cristiani, R. Maiolino, E. Oliva, L. Origlia, M. Riva, L. Valenziano, F. M. Zerbi, M. Abreu, V. Adibekyan, C. Allende Prieto, P. J. Amado, W. Benz, I. Boisse, X. Bonfils, F. Bouchy, L. Buchhave, D. Buscher, A. Cabral, B. L. Canto Martins, A. Chiavassa, J. Coelho, L. B. Christensen , et al. (48 additional authors not shown)

Abstract: The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and… ▽ More The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study. △ Less

Submitted 2 September, 2016; originally announced September 2016.

Comments: 12 pages, in Ground-based and Airborne Instrumentation for Astronomy VI, 2016, Proc. SPIE 9908, 23

Molecular fingerprint-region spectroscopy from 5-12 μみゅーm using an orientation-patterned gallium phosphide optical parametric oscillator

Authors: Luke Maidment, Peter G. Schunemann, Derryck T. Reid

Abstract: We report a femtosecond optical parametric oscillator (OPO) based on the new semiconductor gain material orientation patterned gallium phosphide (OP-GaP), which enables the production of high-repetition-rate femtosecond pulses spanning 5-12 μみゅーm with average powers in the few to tens of milliwatts range. This is the first example of a broadband OPO operating across the molecular fingerprint region,… ▽ More We report a femtosecond optical parametric oscillator (OPO) based on the new semiconductor gain material orientation patterned gallium phosphide (OP-GaP), which enables the production of high-repetition-rate femtosecond pulses spanning 5-12 μみゅーm with average powers in the few to tens of milliwatts range. This is the first example of a broadband OPO operating across the molecular fingerprint region, and we demonstrate its potential by conducting broadband Fourier-transform spectroscopy using water vapor and a polystyrene reference standard. △ Less

Submitted 30 June, 2016; originally announced June 2016.

Comments: 4 pages, 5 figures

Fluctuating volume-current formulation of electromagnetic fluctuations in inhomogeneous media: incandecence and luminescence in arbitrary geometries

Authors: Athanasios G. Polimeridis, M. T. Homer Reid, Weiliang Jin, Steven G. Johnson, Jacob K. White, Alejandro W. Rodriguezz

Abstract: We describe a fluctuating volume--current formulation of electromagnetic fluctuations that extends our recent work on heat exchange and Casimir interactions between arbitrarily shaped homogeneous bodies [Phys. Rev. B. 88, 054305] to situations involving incandescence and luminescence problems, including thermal radiation, heat transfer, Casimir forces, spontaneous emission, fluorescence, and Raman… ▽ More We describe a fluctuating volume--current formulation of electromagnetic fluctuations that extends our recent work on heat exchange and Casimir interactions between arbitrarily shaped homogeneous bodies [Phys. Rev. B. 88, 054305] to situations involving incandescence and luminescence problems, including thermal radiation, heat transfer, Casimir forces, spontaneous emission, fluorescence, and Raman scattering, in inhomogeneous media. Unlike previous scattering formulations based on field and/or surface unknowns, our work exploits powerful techniques from the volume--integral equation (VIE) method, in which electromagnetic scattering is described in terms of volumetric, current unknowns throughout the bodies. The resulting trace formulas (boxed equations) involve products of well-studied VIE matrices and describe power and momentum transfer between objects with spatially varying material properties and fluctuation characteristics. We demonstrate that thanks to the low-rank properties of the associatedmatrices, these formulas are susceptible to fast-trace computations based on iterative methods, making practical calculations tractable. We apply our techniques to study thermal radiation, heat transfer, and fluorescence in complicated geometries, checking our method against established techniques best suited for homogeneous bodies as well as applying it to obtain predictions of radiation from complex bodies with spatially varying permittivities and/or temperature profiles. △ Less

Submitted 19 May, 2015; originally announced May 2015.

Journal ref: Phys. Rev. B 92, 134202 (2015)

Fundamental limits to optical response in absorptive systems

Authors: Owen D. Miller, Athanasios G. Polimeridis, M. T. Homer Reid, Chia Wei Hsu, Brendan G. DeLacy, John D. Joannopoulos, Marin Soljačić, Steven G. Johnson

Abstract: At visible and infrared frequencies, metals show tantalizing promise for strong subwavelength resonances, but material loss typically dampens the response. We derive fundamental limits to the optical response of absorptive systems, bounding the largest enhancements possible given intrinsic material losses. Through basic conservation-of-energy principles, we derive geometry-independent limits to pe… ▽ More At visible and infrared frequencies, metals show tantalizing promise for strong subwavelength resonances, but material loss typically dampens the response. We derive fundamental limits to the optical response of absorptive systems, bounding the largest enhancements possible given intrinsic material losses. Through basic conservation-of-energy principles, we derive geometry-independent limits to per-volume absorption and scattering rates, and to local-density-of-states enhancements that represent the power radiated or expended by a dipole near a material body. We provide examples of structures that approach our absorption and scattering limits at any frequency, by contrast, we find that common "antenna" structures fall far short of our radiative LDOS bounds, suggesting the possibility for significant further improvement. Underlying the limits is a simple metric, $|χかい|^2 / \operatorname{Im} χかい$ for a material with susceptibility $χかい$, that enables broad technological evaluation of lossy materials across optical frequencies. △ Less

Submitted 6 January, 2016; v1 submitted 12 March, 2015; originally announced March 2015.

Comments: 21 pages and 6 figures (excluding appendices, references)

Journal ref: Opt. Express 24, 3329-3364 (2016)

On the Computation of Power in Volume Integral Equation Formulations

Authors: A. G. Polimeridis, M. T. H. Reid, S. G. Johnson, J. K. White, A. W. Rodriguez

Abstract: We present simple and stable formulas for computing power (including absorbed/radiated, scattered and extinction power) in current-based volume integral equation formulations. The proposed formulas are given in terms of vector-matrix-vector products of quantities found solely in the associated linear system. In addition to their efficiency, the derived expressions can guarantee the positivity of t… ▽ More We present simple and stable formulas for computing power (including absorbed/radiated, scattered and extinction power) in current-based volume integral equation formulations. The proposed formulas are given in terms of vector-matrix-vector products of quantities found solely in the associated linear system. In addition to their efficiency, the derived expressions can guarantee the positivity of the computed power. We also discuss the application of Poynting's theorem for the case of sources immersed in dissipative materials. The formulas are validated against results obtained both with analytical and numerical methods for scattering and radiation benchmark cases. △ Less

Submitted 12 September, 2014; v1 submitted 27 June, 2014; originally announced June 2014.

Transformation Optics scheme for two-dimensional materials

Authors: Anshuman Kumar, Kin Hung Fung, M. T. Homer Reid, Nicholas Fang

Abstract: Two dimensional optical materials, such as graphene can be characterized by a surface conductivity. So far, the transformation optics schemes have focused on three dimensional properties such as permittivity $εいぷしろん$ and permeability $μみゅー$. In this paper, we use a scheme for transforming surface currents to highlight that the surface conductivity transforms in a way different from $εいぷしろん$ and $μみゅー$. We use thi… ▽ More Two dimensional optical materials, such as graphene can be characterized by a surface conductivity. So far, the transformation optics schemes have focused on three dimensional properties such as permittivity $εいぷしろん$ and permeability $μみゅー$. In this paper, we use a scheme for transforming surface currents to highlight that the surface conductivity transforms in a way different from $εいぷしろん$ and $μみゅー$. We use this surface conductivity transformation to demonstrate an example problem of reducing scattering of plasmon mode from sharp protrusions in graphene. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Journal ref: Opt. Lett. 39, 2113-2116 (2014)

Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations

Authors: M. T. Homer Reid, Steven G. Johnson, Jacob K. White

Abstract: We present a generic technique, automated by computer-algebra systems and available as open-source software \cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solv… ▽ More We present a generic technique, automated by computer-algebra systems and available as open-source software \cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solvers has often required the aggregation of multiple disparate integral-evaluation schemes to treat all of the distinct types of integrals needed for a given BEM formulation; in contrast, our technique allows many different types of integrals to be handled by the \emph{same} algorithm and the same code implementation. Our method is a significant generalization of the Taylor--Duffy approach \cite{Taylor2003,Duffy1982}, which was originally presented for just a single type of integrand; in addition to generalizing this technique to a broad class of integrands, we also achieve a significant improvement in its efficiency by showing how the \emph{dimension} of the final numerical integral may often be reduced by one. In particular, if $n$ is the number of common vertices between the two triangles, in many cases we can reduce the dimension of the integral from $4-n$ to $3-n$, obtaining a closed-form analytical result for $n=3$ (the common-triangle case). △ Less

Submitted 19 July, 2014; v1 submitted 5 December, 2013; originally announced December 2013.

Comments: Revised version containing a rearrangement of the text and updates to figures

Journal ref: IEEE Transactions on Antennas and Propagation 63 1 (2015) 195-209