Showing 1–50 of 287 results for author: Robinson, D

Microphysical Prescriptions for Parameterized Water Cloud Formation on Ultra-cool Substellar Objects

Authors: James Mang, Caroline V. Morley, Tyler D. Robinson, Peter Gao

Abstract: Water must condense into ice clouds in the coldest brown dwarfs and exoplanets. When they form, these icy clouds change the emergent spectra, temperature structure, and albedo of the substellar atmosphere. The properties of clouds are governed by complex microphysics but these complexities are often not captured by the simpler parameterized cloud models used in climate models or retrieval models.… ▽ More Water must condense into ice clouds in the coldest brown dwarfs and exoplanets. When they form, these icy clouds change the emergent spectra, temperature structure, and albedo of the substellar atmosphere. The properties of clouds are governed by complex microphysics but these complexities are often not captured by the simpler parameterized cloud models used in climate models or retrieval models. Here, we combine microphysical cloud modeling and 1D climate modeling to incorporate insights from microphysical models into a self-consistent, parameterized cloud model. Using the 1D Community Aerosol and Radiation Model for Atmospheres (CARMA), we generate microphysical water clouds and compare their properties with those from the widely-used EddySed cloud model (Ackerman & Marley 2001) for a grid of Y dwarfs. We find that the mass of water condensate in our CARMA water clouds is significantly limited by available condensation nuclei; in models without additional seed particles for clouds added, the atmosphere becomes supersaturated. We incorporate water latent heat release in the convective and radiative parts of the atmosphere and find no significant impact on water-ice cloud formation for typical gas giant compositions. Our analysis reveals the CARMA cloud profiles have a gradual decrease in opacity of approximately 4% per bar below the cloud base. Incorporating this gradual cloud base falloff and a variable $f_{sed}$ parameter allows spectra generated from the parameterized Eddysed model to better match those of the microphysical CARMA model. This work provides recommendations for efficiently generating microphysically-informed water clouds for future models of cold substellar objects with H/He atmospheres. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: 20 pages, 19 figures. Accepted for publication in ApJ

Imagen 3

Authors: Imagen-Team-Google, :, Jason Baldridge, Jakob Bauer, Mukul Bhutani, Nicole Brichtova, Andrew Bunner, Kelvin Chan, Yichang Chen, Sander Dieleman, Yuqing Du, Zach Eaton-Rosen, Hongliang Fei, Nando de Freitas, Yilin Gao, Evgeny Gladchenko, Sergio Gómez Colmenarejo, Mandy Guo, Alex Haig, Will Hawkins, Hexiang Hu, Huilian Huang, Tobenna Peter Igwe, Christos Kaplanis, Siavash Khodadadeh , et al. (227 additional authors not shown)

Abstract: We introduce Imagen 3, a latent diffusion model that generates high quality images from text prompts. We describe our quality and responsibility evaluations. Imagen 3 is preferred over other state-of-the-art (SOTA) models at the time of evaluation. In addition, we discuss issues around safety and representation, as well as methods we used to minimize the potential harm of our models. We introduce Imagen 3, a latent diffusion model that generates high quality images from text prompts. We describe our quality and responsibility evaluations. Imagen 3 is preferred over other state-of-the-art (SOTA) models at the time of evaluation. In addition, we discuss issues around safety and representation, as well as methods we used to minimize the potential harm of our models. △ Less

Submitted 13 August, 2024; originally announced August 2024.

MapTune: Advancing ASIC Technology Mapping via Reinforcement Learning Guided Library Tuning

Authors: Mingju Liu, Daniel Robinson, Yingjie Li, Cunxi Yu

Abstract: Technology mapping involves mapping logical circuits to a library of cells. Traditionally, the full technology library is used, leading to a large search space and potential overhead. Motivated by randomly sampled technology mapping case studies, we propose MapTune framework that addresses this challenge by utilizing reinforcement learning to make design-specific choices during cell selection. By… ▽ More Technology mapping involves mapping logical circuits to a library of cells. Traditionally, the full technology library is used, leading to a large search space and potential overhead. Motivated by randomly sampled technology mapping case studies, we propose MapTune framework that addresses this challenge by utilizing reinforcement learning to make design-specific choices during cell selection. By learning from the environment, MapTune refines the cell selection process, resulting in a reduced search space and potentially improved mapping quality. The effectiveness of MapTune is evaluated on a wide range of benchmarks, different technology libraries and technology mappers. The experimental results demonstrate that MapTune achieves higher mapping accuracy and reducing delay/area across diverse circuit designs, technology libraries and mappers. The paper also discusses the Pareto-Optimal exploration and confirms the perpetual delay-area trade-off. Conducted on benchmark suites ISCAS 85/89, ITC/ISCAS 99, VTR8.0 and EPFL benchmarks, the post-technology mapping and post-sizing quality-of-results (QoR) have been significantly improved, with average Area-Delay Product (ADP) improvement of 22.54\% among all different exploration settings in MapTune. The improvements are consistently remained for four different technologies (7nm, 45nm, 130nm, and 180 nm) and two different mappers. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Comments: IEEE/ACM International Conference on Computer-Aided Design (ICCAD '24), October 27--31, 2024

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Authors: David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used i… ▽ More Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHAP importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere. We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 12 + 5 pages, 7 + 4 figures. To be submitted to OJA. Comments welcome!

Technical design report for the CODEX-$β$ demonstrator

Authors: CODEX-b collaboration, :, Giulio Aielli, Juliette Alimena, James Beacham, Eli Ben Haim, Andras Burucs, Roberto Cardarelli, Matthew Charles, Xabier Cid Vidal, Albert De Roeck, Biplab Dey, Silviu Dobrescu, Ozgur Durmus, Mohamed Elashri, Vladimir Gligorov, Rebeca Gonzalez Suarez, Thomas Gorordo, Zarria Gray, Conor Henderson, Louis Henry, Philip Ilten, Daniel Johnson, Jacob Kautz, Simon Knapen , et al. (28 additional authors not shown)

Abstract: The CODEX-$β$ apparatus is a demonstrator for the proposed future CODEX-b experiment, a long-lived-particle detector foreseen for operation at IP8 during HL-LHC data-taking. The demonstrator project, intended to collect data in 2025, is described, with a particular focus on the design, construction, and installation of the new apparatus. The CODEX-$β$ apparatus is a demonstrator for the proposed future CODEX-b experiment, a long-lived-particle detector foreseen for operation at IP8 during HL-LHC data-taking. The demonstrator project, intended to collect data in 2025, is described, with a particular focus on the design, construction, and installation of the new apparatus. △ Less

Submitted 22 May, 2024; originally announced June 2024.

Influence of Orbit and Mass Constraints on Reflected Light Characterization of Directly Imaged Rocky Exoplanets

Authors: Arnaud Salvador, Tyler D. Robinson, Jonathan J. Fortney, Mark S. Marley

Abstract: Survey strategies for upcoming exoplanet direct imaging missions have considered varying assumptions of prior knowledge. Precursor radial velocity surveys could have detected nearby exo-Earths and provided prior orbit and mass constraints. Alternatively, a direct imaging mission performing astrometry could yield constraints on orbit and phase angle of target planets. Understanding the impact of pr… ▽ More Survey strategies for upcoming exoplanet direct imaging missions have considered varying assumptions of prior knowledge. Precursor radial velocity surveys could have detected nearby exo-Earths and provided prior orbit and mass constraints. Alternatively, a direct imaging mission performing astrometry could yield constraints on orbit and phase angle of target planets. Understanding the impact of prior mass and orbit information on planetary characterization is crucial for efficiently recognizing habitable exoplanets. To address this question, we use a reflected-light retrieval tool to infer the atmospheric and bulk properties of directly imaged Earth-analogs while considering varying levels of prior information and signal-to-noise ratio (SNR). Because of the strong correlation between the orbit-related parameters and the planetary radius, prior information on the orbital distance and planetary phase yield tight constraints on the planetary radius: from $R_{\rm{p}}=2.95^{+2.69}_{-1.95}~R_{\oplus}$ without prior knowledge, to $R_{\rm{p}}=1.01^{+0.33}_{-0.19}~R_{\oplus}$ with prior determination of the orbit for $\rm{SNR}=20$ in the visible/near-infrared spectral range, thus allowing size determination from reflected light observations. However, additional knowledge of planet mass does not notably enhance radius ($R_{\rm{p}}=0.98^{+0.17}_{-0.14}~R_{\oplus}$) or atmospheric characterization. Also, prior knowledge of the mass alone does not yield a tight radius constraint ($R_{\rm{p}}=1.64^{+1.29}_{-0.80}~R_{\oplus}$) nor improves atmospheric composition inference. By contrast, because of its sensitivity to gas column abundance, detecting a Rayleigh scattering slope or bounding Rayleigh opacity helps to refine gas mixing ratio inferences without requiring prior mass knowledge. Overall, apart from radius determination, increasing the SNR is more beneficial than additional prior observations. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJL

Loss-Versus-Fair: Efficiency of Dutch Auctions on Blockchains

Authors: Ciamac C. Moallemi, Dan Robinson

Abstract: Milionis et al.(2023) studied the rate at which automated market makers leak value to arbitrageurs when block times are discrete and follow a Poisson process, and where the risky asset price follows a geometric Brownian motion. We extend their model to analyze another popular mechanism in decentralized finance for onchain trading: Dutch auctions. We compute the expected losses that a seller incurs… ▽ More Milionis et al.(2023) studied the rate at which automated market makers leak value to arbitrageurs when block times are discrete and follow a Poisson process, and where the risky asset price follows a geometric Brownian motion. We extend their model to analyze another popular mechanism in decentralized finance for onchain trading: Dutch auctions. We compute the expected losses that a seller incurs to arbitrageurs and expected time-to-fill for Dutch auctions as a function of starting price, volatility, decay rate, and average interblock time. We also extend the analysis to gradual Dutch auctions, a variation on Dutch auctions for selling tokens over time at a continuous rate. We use these models to explore the tradeoff between speed of execution and quality of execution, which could help inform practitioners in setting parameters for starting price and decay rate on Dutch auctions, or help platform designers determine performance parameters like block times. △ Less

Submitted 26 July, 2024; v1 submitted 31 May, 2024; originally announced June 2024.

Paths to Robust Exoplanet Science Yield Margin for the Habitable Worlds Observatory

Authors: Christopher C. Stark, Bertrand Mennesson, Steve Bryson, Eric B. Ford, Tyler D. Robinson, Ruslan Belikov, Matthew R. Bolcar, Lee D. Feinberg, Olivier Guyon, Natasha Latouf, Avi M. Mandell, Bernard J. Rauscher, Dan Sirbu, Noah W. Tuchow

Abstract: The Habitable Worlds Observatory (HWO) will seek to detect and characterize potentially Earth-like planets around other stars. To ensure that the mission achieves the Astro2020 Decadal's recommended goal of 25 exoEarth candidates (EECs), we must take into account the probabilistic nature of exoplanet detections and provide "science margin" to budget for astrophysical uncertainties with a reasonabl… ▽ More The Habitable Worlds Observatory (HWO) will seek to detect and characterize potentially Earth-like planets around other stars. To ensure that the mission achieves the Astro2020 Decadal's recommended goal of 25 exoEarth candidates (EECs), we must take into account the probabilistic nature of exoplanet detections and provide "science margin" to budget for astrophysical uncertainties with a reasonable level of confidence. In this study, we explore the probabilistic distributions of yields to be expected from a blind exoEarth survey conducted by such a mission. We identify and estimate the impact of all major known sources of astrophysical uncertainty on the exoEarth candidate yield. As expected, eta_Earth uncertainties dominate the uncertainty in EEC yield, but we show that sampling uncertainties inherent to a blind survey are another important source of uncertainty that should be budgeted for during mission design. We adopt the Large UV/Optical/IR Surveyor Design B (LUVOIR-B) as a baseline and modify the telescope diameter to estimate the science margin provided by a larger telescope. We then depart from the LUVOIR-B baseline design and identify six possible design changes that, when compiled, provide large gains in exoEarth candidate yield and more than an order of magnitude reduction in exposure times for the highest priority targets. We conclude that a combination of telescope diameter increase and design improvements could provide robust exoplanet science margins for HWO. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 90 pages, 26 figures, 4 tables, Submitted to JATIS

Requirements are All You Need: The Final Frontier for End-User Software Engineering

Authors: Diana Robinson, Christian Cabrera, Andrew D. Gordon, Neil D. Lawrence, Lars Mennen

Abstract: What if end users could own the software development lifecycle from conception to deployment using only requirements expressed in language, images, video or audio? We explore this idea, building on the capabilities that generative Artificial Intelligence brings to software generation and maintenance techniques. How could designing software in this way better serve end users? What are the implicati… ▽ More What if end users could own the software development lifecycle from conception to deployment using only requirements expressed in language, images, video or audio? We explore this idea, building on the capabilities that generative Artificial Intelligence brings to software generation and maintenance techniques. How could designing software in this way better serve end users? What are the implications of this process for the future of end-user software engineering and the software development lifecycle? We discuss the research needed to bridge the gap between where we are today and these imagined systems of the future. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Accepted at International Workshop on Software Engineering 2030 in Porto de Galinhas, Brazil (July 2024)

Impact of Planetary Parameters on Water Clouds Microphysics

Authors: Huanzhou Yang, Thaddeus D. Komacek, Owen B. Toon, Eric T. Wolf, Tyler D. Robinson, Caroline Chael, Dorian S. Abbot

Abstract: Potentially habitable exoplanets are targets of great interest for the James Webb Space Telescope and upcoming mission concepts such as the Habitable Worlds Observatory. Clouds strongly affect climate and habitability, but predicting their properties is difficult. In Global Climate Models (GCMs), especially those aiming at simulating Earth, cloud microphysics is often crudely approximated by assum… ▽ More Potentially habitable exoplanets are targets of great interest for the James Webb Space Telescope and upcoming mission concepts such as the Habitable Worlds Observatory. Clouds strongly affect climate and habitability, but predicting their properties is difficult. In Global Climate Models (GCMs), especially those aiming at simulating Earth, cloud microphysics is often crudely approximated by assuming that all cloud particles have a single, constant size or a prescribed size distribution and that all clouds in a grid cell are identical. For exoplanets that range over a large phase space of planetary properties, this method could result in large errors. In this work, our goal is to determine how cloud microphysics on terrestrial exoplanets, whose condensable is mainly water vapor, depend on aerosol properties and planetary parameters such as surface pressure, surface gravity, and incident stellar radiation. We use the Community Aerosol and Radiation Model for Atmospheres as a 1D microphysical model to simulate the formation and evolution of clouds including the processes of nucleation, condensation, evaporation, coagulation, and vertical transfer. In these 1D idealized experiments, we find that the parameters that determine the macrophysical thermal structure of the atmospheres, including surface pressure and stellar flux, impact cloud radiative effect (CRE) most significantly. Parameters such as gravity and number density of aerosols working as cloud condensation nuclei affect the microphysical processes of cloud formation, including activation and vertical transfer. They also have a significant, though weaker effect on CRE. This work motivates the development of more accurate GCM cloud schemes and should aid in the interpretation of future observations. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 17 pages, 10 figures

Journal ref: Huanzhou Yang et al 2024 ApJ 966 152

Exomoons & Exorings with the Habitable Worlds Observatory I: On the Detection of Earth-Moon Analog Shadows & Eclipses

Authors: Mary Anne Limbach, Jacob Lustig-Yaeger, Andrew Vanderburg, Johanna M. Vos, Rene Heller, Tyler D. Robinson

Abstract: The highest priority recommendation of the Astro2020 Decadal Survey for space-based astronomy was the construction of an observatory capable of characterizing habitable worlds. In this paper series we explore the detectability of and interference from exomoons and exorings serendipitously observed with the proposed Habitable Worlds Observatory (HWO) as it seeks to characterize exoplanets, starting… ▽ More The highest priority recommendation of the Astro2020 Decadal Survey for space-based astronomy was the construction of an observatory capable of characterizing habitable worlds. In this paper series we explore the detectability of and interference from exomoons and exorings serendipitously observed with the proposed Habitable Worlds Observatory (HWO) as it seeks to characterize exoplanets, starting in this manuscript with Earth-Moon analog mutual events. Unlike transits, which only occur in systems viewed near edge-on, shadow (i.e., solar eclipse) and lunar eclipse mutual events occur in almost every star-planet-moon system. The cadence of these events can vary widely from ~yearly to multiple events per day, as was the case in our younger Earth-Moon system. Leveraging previous space-based (EPOXI) lightcurves of a Moon transit and performance predictions from the LUVOIR-B concept, we derive the detectability of Moon analogs with HWO. We determine that Earth-Moon analogs are detectable with observation of ~2-20 mutual events for systems within 10pc, and larger moons should remain detectable out to 20pc. We explore the extent to which exomoon mutual events can mimic planet features and weather. We find that HWO wavelength coverage in the near-IR, specifically in the 1.4 micron water band where large moons can outshine their host planet, will aid in differentiating exomoon signals from exoplanet variability. Finally, we predict that exomoons formed through collision processes akin to our Moon are more likely to be detected in younger systems, where shorter orbital periods and favorable geometry enhance the probability and frequency of mutual events. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 16 pages, 11 figures, in review at AJ, includes revisions from first round of review, comments welcome

Detectability of Surface Biosignatures for Directly-Imaged Rocky Exoplanets

Authors: Schuyler R. Borges, Gabrielle G. Jones, Tyler D. Robinson

Abstract: Modeling the detection of life has never been more opportune. With next generation space telescopes, like the currently developing Habitable Worlds Observatory (HWO) concept, we will begin to characterize rocky exoplanets potentially similar to Earth. However, currently, few realistic planetary spectra containing surface biosignatures have been paired with direct imaging telescope instrument model… ▽ More Modeling the detection of life has never been more opportune. With next generation space telescopes, like the currently developing Habitable Worlds Observatory (HWO) concept, we will begin to characterize rocky exoplanets potentially similar to Earth. However, currently, few realistic planetary spectra containing surface biosignatures have been paired with direct imaging telescope instrument models. Therefore, we use a HWO instrument noise model to assess the detection of surface biosignatures affiliated with oxygenic, anoxygenic, and nonphotosynthetic extremophiles. We pair the HWO telescope model to a 1-D radiative transfer model to estimate the required exposure times necessary for detecting each biosignature on planets with global microbial coverage and varying atmospheric water vapor concentrations. For modeled planets with 0% - 50% cloud coverage, we determine pigments and the red edge could be detected within 1,000 hours (100 hours) at distances within 15 pc (11 pc). However, tighter telescope inner working angles (2.5 lambda/D) would allow surface biosignature detection at further distances. Anoxygenic photosynthetic biosignatures could also be more easily detectable than nonphotosynthetic pigments and the photosynthetic red edge when compared against a false positive iron oxide slope. Future life detection missions should evaluate the influence of false positives on the detection of multiple surface biosignatures. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: in press; as-accepted version, prior to journal editing

A Proximal-Gradient Method for Constrained Optimization

Authors: Yutong Dai, Xiaoyi Qu, Daniel P. Robinson

Abstract: We present a new algorithm for solving optimization problems with objective functions that are the sum of a smooth function and a (potentially) nonsmooth regularization function, and nonlinear equality constraints. The algorithm may be viewed as an extension of the well-known proximal-gradient method that is applicable when constraints are not present. To account for nonlinear equality constraints… ▽ More We present a new algorithm for solving optimization problems with objective functions that are the sum of a smooth function and a (potentially) nonsmooth regularization function, and nonlinear equality constraints. The algorithm may be viewed as an extension of the well-known proximal-gradient method that is applicable when constraints are not present. To account for nonlinear equality constraints, we combine a decomposition procedure for computing trial steps with an exact merit function for determining trial step acceptance. Under common assumptions, we show that both the proximal parameter and merit function parameter eventually remain fixed, and then prove a worst-case complexity result for the maximum number of iterations before an iterate satisfying approximate first-order optimality conditions for a given tolerance is computed. Our preliminary numerical results indicate that our approach has great promise, especially in terms of returning approximate solutions that are structured (e.g., sparse solutions when a one-norm regularizer is used). △ Less

Submitted 10 April, 2024; originally announced April 2024.

MSC Class: 49M37; 65K05; 65K10; 65Y20; 68Q25; 90C30; 90C60

am-AMM: An Auction-Managed Automated Market Maker

Authors: Austin Adams, Ciamac C. Moallemi, Sara Reynolds, Dan Robinson

Abstract: Automated market makers (AMMs) have emerged as the dominant market mechanism for trading on decentralized exchanges implemented on blockchains. This paper presents a single mechanism that targets two important unsolved problems for AMMs: reducing losses to informed orderflow, and maximizing revenue from uninformed orderflow. The ``auction-managed AMM'' works by running a censorship-resistant oncha… ▽ More Automated market makers (AMMs) have emerged as the dominant market mechanism for trading on decentralized exchanges implemented on blockchains. This paper presents a single mechanism that targets two important unsolved problems for AMMs: reducing losses to informed orderflow, and maximizing revenue from uninformed orderflow. The ``auction-managed AMM'' works by running a censorship-resistant onchain auction for the right to temporarily act as ``pool manager'' for a constant-product AMM. The pool manager sets the swap fee rate on the pool, and also receives the accrued fees from swaps. The pool manager can exclusively capture some arbitrage by trading against the pool in response to small price movements, and also can set swap fees incorporating price sensitivity of retail orderflow and adapting to changing market conditions, with the benefits from both ultimately accruing to liquidity providers. Liquidity providers can enter and exit the pool freely in response to changing rent, though they must pay a small fee on withdrawal. We prove that under certain assumptions, this AMM should have higher liquidity in equilibrium than any standard, fixed-fee AMM. △ Less

Submitted 23 May, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Ten simple rules for collaborating with wet lab researchers for computational researchers

Authors: Mark D. Robinson, Peiying Cai, Martin Emons, Reto Gerber, Pierre-Luc Germain, Samuel Gunz, Siyuan Luo, Giulia Moro, Emanuel Sonder, Anthony Sonrel, Jiayi Wang, David Wissel, Izaskun Mallona

Abstract: Computational biologists are frequently engaged in collaborative data analysis with wet lab researchers. These interdisciplinary projects, as necessary as they are to the scientific endeavour, can be surprisingly challenging due to cultural differences in operations and values. In these Ten Simple Rules guide we aim to help dry lab researchers identify sources of friction; and provide actionable t… ▽ More Computational biologists are frequently engaged in collaborative data analysis with wet lab researchers. These interdisciplinary projects, as necessary as they are to the scientific endeavour, can be surprisingly challenging due to cultural differences in operations and values. In these Ten Simple Rules guide we aim to help dry lab researchers identify sources of friction; and provide actionable tools to facilitate respectful, open, transparent and rewarding collaborations. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 8 pages, 1 figure

On-shell recursion and holomorphic HQET for heavy quark hadronic resonances

Authors: Claudio Andrea Manzari, Dean J. Robinson

Abstract: We develop a new theoretical framework for the treatment of heavy quark (HQ) resonances within heavy quark effective theory (HQET). This framework uses on-shell recursion techniques to express the resonant amplitude as a product of on-shell subamplitudes, which allows one to employ a form-factor representation of the hadronic matrix elements and to obtain an HQ expansion, but at the price of intro… ▽ More We develop a new theoretical framework for the treatment of heavy quark (HQ) resonances within heavy quark effective theory (HQET). This framework uses on-shell recursion techniques to express the resonant amplitude as a product of on-shell subamplitudes, which allows one to employ a form-factor representation of the hadronic matrix elements and to obtain an HQ expansion, but at the price of introducing complex momenta. We construct a generalized "holomorphic HQET" onto which such complex-momentum matrix elements can be matched, and we show that $PT$ symmetry ensures the Isgur-Wise functions (and the perturbative corrections) become holomorphic functions of the complex recoil parameter with real coefficients. They are thus an analytic continuation of the standard HQET description. This framework admits a HQ hadron (strong decay) width expansion. At second order, we show it is in good agreement with data for the $B_{1(2)}^{(*)}$ and $D_{1(2)}^{(*)}$ HQ doublets. Taking the $\bar{B} \to (D_1^*(1^-) \to Dπ)lν$ system as an example, we compute the holomorphic HQET expansion to first order, as well as the complex-momentum on-shell subamplitudes. A toy numerical study of the resulting differential rates demonstrates that this framework generates HQ resonance lineshapes with large tails, resembling those seen in data. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 65 pages, 6 figures, 1 table

Exoplanet Analog Observations of Earth from Galileo Disk-integrated Photometry

Authors: Ryder H. Strauss, Tyler D. Robinson, David E. Trilling, Ryan Cummings, Christopher J. Smith

Abstract: The Galileo spacecraft had distant encounters with Earth in 1990 and 1992. Limited Solid State Imager (SSI) data acquired during these encounters has been previously presented, but the majority of the data from these Earth flybys have not been presented in the literature. Observations of Earth taken from afar are both rare and directly relevant to the development of any future exo-Earth direct ima… ▽ More The Galileo spacecraft had distant encounters with Earth in 1990 and 1992. Limited Solid State Imager (SSI) data acquired during these encounters has been previously presented, but the majority of the data from these Earth flybys have not been presented in the literature. Observations of Earth taken from afar are both rare and directly relevant to the development of any future exo-Earth direct imaging mission. Here we present a pipeline that vets, calibrates, and measures the disk-integrated brightness of the Earth, in multiple filters, from the complete SSI data sets from both the 1990 and 1992 Galileo flybys. The result is over 1500 usable photometric measurements for Earth as an analog for an exoplanet. The 1990 data set includes full rotational lightcurves in six bandpasses spanning the optical range. The 1992 data set is more limited, with lightcurves only spanning 14 hr. Time-averaged photometry for both encounters is presented while variability and color are discussed relative to findings from NASA's EPOXI mission (which also provided photometric lighturves for Earth). The new Galileo/SSI data are used to further validate the Virtual Planetary Laboratory 3D spectral Earth model, which often serves as a stand-in for true disk-integrated observations of our planet. The revived Galileo/SSI data for Earth is a testament to the ability of NASA's Planetary Data System to maintain data over decades-long timescales. The disk-integrated products derived from these data add to a very short list of calibrated and published whole-disk observations of the Pale Blue Dot. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Journal ref: AJ 167 87 (2024)

Inner Edge Habitable Zone Limits Around Main Sequence Stars: Cloudy Estimates

Authors: James D. Windsor, Tyler D. Robinson, Ravi kumar Kopparapu, Arnaud Salvador, Amber V. Young, Victoria S. Meadows

Abstract: Understanding the limits of rocky planet habitability is one of the key goals of current and future exoplanet characterization efforts. An intrinsic concept of rocky planet habitability is the Habitable Zone. To date, the most widely used estimates of the Habitable Zone are based on cloud-free, one-dimensional (vertical) radiative-convective climate model calculations. However, recent three-dimens… ▽ More Understanding the limits of rocky planet habitability is one of the key goals of current and future exoplanet characterization efforts. An intrinsic concept of rocky planet habitability is the Habitable Zone. To date, the most widely used estimates of the Habitable Zone are based on cloud-free, one-dimensional (vertical) radiative-convective climate model calculations. However, recent three-dimensional global climate modeling efforts have revealed that rocky planet habitability is strongly impacted by radiative cloud feedbacks, where computational expense and model limitations can prevent these tools from exploring the limits of habitability across the full range of parameter space. We leverage a patchy cloud one-dimensional radiative-convective climate model with parameterized cloud microphysics to investigate Inner Edge limits to the Habitable Zone for main sequence stars ($T_{\rm eff}$ = 2600 -7200K). We find that Inner Edge limits to the Habitable Zone can be 3.3 and 4.7 times closer than previous cloud-free estimates for Earth- and super-Earth-sized worlds, respectively, depending on bulk cloud parameters (e.g., fractional cloudiness and sedimentation efficiency). These warm, moist Inner Edge climates are expected to have extensive cloud decks that could mute deep atmosphere spectral features. To aid in rocky planet characterization studies, we identify the potential of using $\rm{CO_{\rm 2}}$ absorption features in transmission spectroscopy as a means of quantifying cloud deck height and cloud sedimentation efficiency. Moist greenhouse climates may represent key yet poorly understood states of habitable planets for which continued study will uncover new insights into the search and characterization of habitable worlds. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 33 pages, 12 figures, submitted, community comments welcome

Perceptual and technical barriers in sharing and formatting metadata accompanying omics studies

Authors: Yu-Ning Huang, Michael I. Love, Cynthia Flaire Ronkowski, Dhrithi Deshpande, Lynn M. Schriml, Annie Wong-Beringer, Barend Mons, Russell Corbett-Detig, Christopher I Hunter, Jason H. Moore, Lana X. Garmire, T. B. K. Reddy, Winston A. Hide, Atul J. Butte, Mark D. Robinson, Serghei Mangul

Abstract: Metadata, often termed "data about data," is crucial for organizing, understanding, and managing vast omics datasets. It aids in efficient data discovery, integration, and interpretation, enabling users to access, comprehend, and utilize data effectively. Its significance spans the domains of scientific research, facilitating data reproducibility, reusability, and secondary analysis. However, nume… ▽ More Metadata, often termed "data about data," is crucial for organizing, understanding, and managing vast omics datasets. It aids in efficient data discovery, integration, and interpretation, enabling users to access, comprehend, and utilize data effectively. Its significance spans the domains of scientific research, facilitating data reproducibility, reusability, and secondary analysis. However, numerous perceptual and technical barriers hinder the sharing of metadata among researchers. These barriers compromise the reliability of research results and hinder integrative meta-analyses of omics studies . This study highlights the key barriers to metadata sharing, including the lack of uniform standards, privacy and legal concerns, limitations in study design, limited incentives, inadequate infrastructure, and the dearth of well-trained personnel for metadata management and reuse. Proposed solutions include emphasizing the promotion of standardization, educational efforts, the role of journals and funding agencies, incentives and rewards, and the improvement of infrastructure. More accurate, reliable, and impactful research outcomes are achievable if the scientific community addresses these barriers, facilitating more accurate, reliable, and impactful research outcomes. △ Less

Submitted 22 November, 2023; originally announced January 2024.

Study Duration Prediction for Clinical Trials with Time-to-Event Endpoints Using Mixture Distributions Accounting for Heterogeneous Population

Authors: Hong Zhang, Jie Pu, Shibing Deng, Satrajit Roychoudhury, Haitao Chu, Douglas Robinson

Abstract: In the era of precision medicine, more and more clinical trials are now driven or guided by biomarkers, which are patient characteristics objectively measured and evaluated as indicators of normal biological processes, pathogenic processes, or pharmacologic responses to therapeutic interventions. With the overarching objective to optimize and personalize disease management, biomarker-guided clinic… ▽ More In the era of precision medicine, more and more clinical trials are now driven or guided by biomarkers, which are patient characteristics objectively measured and evaluated as indicators of normal biological processes, pathogenic processes, or pharmacologic responses to therapeutic interventions. With the overarching objective to optimize and personalize disease management, biomarker-guided clinical trials increase the efficiency by appropriately utilizing prognostic or predictive biomarkers in the design. However, the efficiency gain is often not quantitatively compared to the traditional all-comers design, in which a faster enrollment rate is expected (e.g. due to no restriction to biomarker positive patients) potentially leading to a shorter duration. To accurately predict biomarker-guided trial duration, we propose a general framework using mixture distributions accounting for heterogeneous population. Extensive simulations are performed to evaluate the impact of heterogeneous population and the dynamics of biomarker characteristics and disease on the study duration. Several influential parameters including median survival time, enrollment rate, biomarker prevalence and effect size are identitied. Re-assessments of two publicly available trials are conducted to empirically validate the prediction accuracy and to demonstrate the practical utility. The R package \emph{detest} is developed to implement the proposed method and is publicly available on CRAN. △ Less

Submitted 31 December, 2023; originally announced January 2024.

Appraising constrained second-order power corrections in HQET with $Λ_b \to Λ_c l ν$

Authors: Florian U. Bernlochner, Michele Papucci, Dean J. Robinson

Abstract: We derive the $Λ_b \to Λ_c$ form factors for the Standard Model and beyond at second order in Heavy Quark Effective Theory (HQET), applying the recently-proposed Residual Chiral Expansion (RCE) to reduce the set of unknown subsubleading hadronic functions to a single, highly-constrained function, that is fully determined by hadron mass parameters at zero recoil. We fit a form factor parametrizatio… ▽ More We derive the $Λ_b \to Λ_c$ form factors for the Standard Model and beyond at second order in Heavy Quark Effective Theory (HQET), applying the recently-proposed Residual Chiral Expansion (RCE) to reduce the set of unknown subsubleading hadronic functions to a single, highly-constrained function, that is fully determined by hadron mass parameters at zero recoil. We fit a form factor parametrization based on these results to all available Lattice QCD (LQCD) predictions and experimental data. We find that the constrained and predictive structure of the form factors under the RCE is in excellent agreement with LQCD predictions and experimental data, as well as prior HQET-based fits. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 37 pages, 5 figures, 10 tables

Report number: CALT-TH-2023-049

Retrievals Applied To A Decision Tree Framework Can Characterize Earth-like Exoplanet Analogs

Authors: Amber V. Young, Jaime Crouse, Giada Arney, Shawn Domagal-Goldman, Tyler D. Robinson, Sandra T. Bastelberger

Abstract: Exoplanet characterization missions planned for the future will soon enable searches for life beyond our solar system. Critical to the search will be the development of life detection strategies that can search for biosignatures while maintaining observational efficiency. In this work, we adopted a newly developed biosignature decision tree strategy for remote characterization of Earth-like exopla… ▽ More Exoplanet characterization missions planned for the future will soon enable searches for life beyond our solar system. Critical to the search will be the development of life detection strategies that can search for biosignatures while maintaining observational efficiency. In this work, we adopted a newly developed biosignature decision tree strategy for remote characterization of Earth-like exoplanets. The decision tree offers a step-by-step roadmap for detecting exoplanet biosignatures and excluding false positives based on Earth's biosphere and its evolution over time. We followed the pathways for characterizing a modern Earth-like planet and an Archean Earth-like planet and evaluated the observational trades associated with coronagraph bandpass combinations of designs consistent with The Habitable Worlds Observatory (HWO) precursor studies. With retrieval analyses of each bandpass (or combination), we demonstrate the utility of the decision tree and evaluated the uncertainty on a suite of biosignature chemical species and habitability indicators (i.e., the gas abundances of H$_2$O, O$_2$, O$_3$, CH$_4$, and CO$_2$). Notably for modern Earth, less than an order of magnitude spread in the 1-$σ$ uncertainties were achieved for the abundances of H$_2$O and O$_2$, planetary surface pressure, and atmospheric temperature with three strategically placed bandpasses (two in the visible and one in the near-infrared). For the Archean, CH$_4$ and H$_2$O were detectable in the visible with a single bandpass. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: Planetary Science Journal. For Source Data see https://zenodo.org/doi/10.5281/zenodo.10105573

Inferring Chemical Disequilibrium Biosignatures for Proterozoic Earth-Like Exoplanets

Authors: Amber V. Young, Tyler D. Robinson, Joshua Krissansen-Totton, Edward W. Schwieterman, Nicholas F. Wogan, Michael J. Way, Linda E. Sohl, Giada N. Arney, Christopher T. Reinhard, Michael R. Line, David C. Catling, James D. Windsor

Abstract: Chemical disequilibrium quantified via available free energy has previously been proposed as a potential biosignature. However, exoplanet biosignature remote sensing work has not yet investigated how observational uncertainties impact the ability to infer a life-generated available free energy. We pair an atmospheric retrieval tool to a thermodynamics model to assess the detectability of chemical… ▽ More Chemical disequilibrium quantified via available free energy has previously been proposed as a potential biosignature. However, exoplanet biosignature remote sensing work has not yet investigated how observational uncertainties impact the ability to infer a life-generated available free energy. We pair an atmospheric retrieval tool to a thermodynamics model to assess the detectability of chemical disequilibrium signatures of Earth-like exoplanets, emphasizing the Proterozoic Eon where atmospheric abundances of oxygen-methane disequilibrium pairs may have been relatively high. Retrieval model studies applied across a range of gas abundances revealed that order-of-magnitude constraints on disequilibrium energy are achieved with simulated reflected-light observations at the high abundance scenario and signal-to-noise ratios (50) while weak constraints are found at moderate SNRs (20\,--\,30) for med\,--\,low abundance cases. Furthermore, the disequilibrium energy constraints are improved by modest thermal information encoded in water vapor opacities at optical and near-infrared wavelengths. These results highlight how remotely detecting chemical disequilibrium biosignatures can be a useful and metabolism-agnostic approach to biosignature detection. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: Nature Astronomy. Supplementary information see https://zenodo.org/records/10093798 For Source Data see https://zenodo.org/records/8335447

Emergence of the Temperature-Density Relation in the Low-Density Intergalactic Medium

Authors: Alexandra Wells, David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The PDF of gas temperature at fixed density exhibits two clear modes: a warm and cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We fin… ▽ More We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The PDF of gas temperature at fixed density exhibits two clear modes: a warm and cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is "universal" in the sense that its timing is accurately parameterized by the value of the volume-weighted neutral fraction for any reionization history. This "universality" is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history and the evolving gas density PDF. Furthermore, the "emergence" of the tight temperature-density relation in the warm mode is also approximately "universally" controlled by the volume-weighted neutral fraction for any reionization history. In particular, the "emergence" of the temperature-density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is $10^{-4}\lesssim X_\mathrm{HI} \lesssim10^{-3}$ for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature-density relation, regardless of reionization history. △ Less

Submitted 22 August, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 7 pages, 10 figures. Updated to match accepted version in MNRAS

Exploring the Dependence of Gas Cooling and Heating Functions on the Incident Radiation Field with Machine Learning

Authors: David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capa… ▽ More Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code Cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (SHAP) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behavior of cooling and heating functions. △ Less

Submitted 18 January, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 12+4 pages, 6+2 figures. Published in MNRAS. Updated to match MNRAS accepted version

Journal ref: MNRAS, Volume 528, Issue 1, February 2024

The Iberian Exception: An overview of its effects over its first 100 days

Authors: David Robinson, Angel Arcos-Vargas, Micheael Tennican, Fernando Núñez

Abstract: This paper offers an independent assessment of certain key economic effects of the Iberian Exception (IE). Their stated aim was to reduce the major component of electricity prices for most Iberian consumers, a component which was indexed to Iberian wholesale power market spot prices power market prices that were rising alarmingly due to extremely tight international markets for natural gas. The Sp… ▽ More This paper offers an independent assessment of certain key economic effects of the Iberian Exception (IE). Their stated aim was to reduce the major component of electricity prices for most Iberian consumers, a component which was indexed to Iberian wholesale power market spot prices power market prices that were rising alarmingly due to extremely tight international markets for natural gas. The Spanish Government estimates that, during its first 100 days, the IE provided substantial benefits for consumers affected by the IE, which included over 10 million small consumers as well as many large ones, but the authors of this study question that estimate. The authors of this paper argue that the estimated effect of the IE on retail prices depends critically on the assumptions about what would have occurred in the absence of the IE, i.e., in a counterfactual scenario. Although counterfactuals are always difficult to construct, the government s counterfactual ignores demand elasticity, and this inflates their estimate of immediate consumer benefits. Using hourly data on the wholesale electricity market for the first 100 days of the IE, this paper s analysis of alternative counterfactuals that reflect the effects of demand elasticity shows substantially lower benefits of the IE for consumers than the Spanish government estimates. Indeed, this paper s analysis suggests that affected consumers would have paid somewhat less for electricity in the first 100 days of the IE had it not been introduced. The authors identify several other potential short and longterm effects of the IE that deserve further study. These include increased margins for fossil fired generators, reduced margins for some decarbonized inframarginal plant, heightened investor perceptions of regulatory risk, weakened incentives for efficient consumption, and higher carbon emissions and gas prices. △ Less

Submitted 25 September, 2023; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: 34 pages, 9 figures and 4 tables

MSC Class: Economics

Transitioning ECP Software Technology into a Foundation for Sustainable Research Software

Authors: Gregory R. Watson, Addi Malviya-Thakur, Daniel S. Katz, Elaine M. Raybourn, Bill Hoffman, Dana Robinson, John Kellerman, Clark Roundy

Abstract: Research software plays a crucial role in advancing scientific knowledge, but ensuring its sustainability, maintainability, and long-term viability is an ongoing challenge. The Sustainable Research Software Institute (SRSI) Model has been designed to address the concerns, and presents a comprehensive framework designed to promote sustainable practices in the research software community. However th… ▽ More Research software plays a crucial role in advancing scientific knowledge, but ensuring its sustainability, maintainability, and long-term viability is an ongoing challenge. The Sustainable Research Software Institute (SRSI) Model has been designed to address the concerns, and presents a comprehensive framework designed to promote sustainable practices in the research software community. However the SRSI Model does not address the transitional requirements for the Exascale Computing Project (ECP) Software Technology (ECP-ST) focus area specifically. This white paper provides an overview and detailed description of how ECP-ST will transition into the SRSI in a compressed time frame that a) meets the needs of the ECP end-of-technical-activities deadline; and b) ensures the continuity of the sustainability efforts that are already underway. △ Less

Submitted 30 August, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 7 pages, 1 figure

Report number: 200366

An Open Community-Driven Model For Sustainable Research Software: Sustainable Research Software Institute

Authors: Gregory R. Watson, Addi Malviya-Thakur, Daniel S. Katz, Elaine M. Raybourn, Bill Hoffman, Dana Robinson, John Kellerman, Clark Roundy

Abstract: Research software plays a crucial role in advancing scientific knowledge, but ensuring its sustainability, maintainability, and long-term viability is an ongoing challenge. To address these concerns, the Sustainable Research Software Institute (SRSI) Model presents a comprehensive framework designed to promote sustainable practices in the research software community. This white paper provides an i… ▽ More Research software plays a crucial role in advancing scientific knowledge, but ensuring its sustainability, maintainability, and long-term viability is an ongoing challenge. To address these concerns, the Sustainable Research Software Institute (SRSI) Model presents a comprehensive framework designed to promote sustainable practices in the research software community. This white paper provides an in-depth overview of the SRSI Model, outlining its objectives, services, funding mechanisms, collaborations, and the significant potential impact it could have on the research software community. It explores the wide range of services offered, diverse funding sources, extensive collaboration opportunities, and the transformative influence of the SRSI Model on the research software landscape △ Less

Submitted 30 August, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 13 pages, 1 figure

Report number: 200363

Earth as a Transiting Exoplanet: A Validation of Transmission Spectroscopy and Atmospheric Retrieval Methodologies for Terrestrial Exoplanets

Authors: Jacob Lustig-Yaeger, Victoria S. Meadows, David Crisp, Michael R. Line, Tyler D. Robinson

Abstract: The James Webb Space Telescope (JWST) will enable the search for and characterization of terrestrial exoplanet atmospheres in the habitable zone via transmission spectroscopy. However, relatively little work has been done to use solar system data, where ground truth is known, to validate spectroscopic retrieval codes intended for exoplanet studies, particularly in the limit of high resolution and… ▽ More The James Webb Space Telescope (JWST) will enable the search for and characterization of terrestrial exoplanet atmospheres in the habitable zone via transmission spectroscopy. However, relatively little work has been done to use solar system data, where ground truth is known, to validate spectroscopic retrieval codes intended for exoplanet studies, particularly in the limit of high resolution and high signal-to-noise (S/N). In this work, we perform such a validation by analyzing a high S/N empirical transmission spectrum of Earth using a new terrestrial exoplanet atmospheric retrieval model with heritage in Solar System remote sensing and gaseous exoplanet retrievals. We fit the Earth's 2-14 um transmission spectrum in low resolution (R=250 at 5 um) and high resolution (R=100,000 at 5 um) under a variety of assumptions about the 1D vertical atmospheric structure. In the limit of noiseless transmission spectra, we find excellent agreement between model and data (deviations < 10%) that enable the robust detection of H2O, CO2, O3, CH4, N2, N2O, NO2, HNO3, CFC-11, and CFC-12 thereby providing compelling support for the detection of habitability, biosignature, and technosignature gases in the atmosphere of the planet using an exoplanet-analog transmission spectrum. Our retrievals at high spectral resolution show a marked sensitivity to the thermal structure of the atmosphere, trace gas abundances, density-dependent effects, such as collision-induced absorption and refraction, and even hint at 3D spatial effects. However, we used synthetic observations of TRAPPIST-1e to verify that the use of simple 1D vertically homogeneous atmospheric models will likely suffice for JWST observations of terrestrial exoplanets transiting M dwarfs. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 37 pages, 14 figures, 9 tables. Accepted for publication in PSJ

Transferable Models for Bioacoustics with Human Language Supervision

Authors: David Robinson, Adelaide Robinson, Lily Akrapongpisak

Abstract: Passive acoustic monitoring offers a scalable, non-invasive method for tracking global biodiversity and anthropogenic impacts on species. Although deep learning has become a vital tool for processing this data, current models are inflexible, typically cover only a handful of species, and are limited by data scarcity. In this work, we propose BioLingual, a new model for bioacoustics based on contra… ▽ More Passive acoustic monitoring offers a scalable, non-invasive method for tracking global biodiversity and anthropogenic impacts on species. Although deep learning has become a vital tool for processing this data, current models are inflexible, typically cover only a handful of species, and are limited by data scarcity. In this work, we propose BioLingual, a new model for bioacoustics based on contrastive language-audio pretraining. We first aggregate bioacoustic archives into a language-audio dataset, called AnimalSpeak, with over a million audio-caption pairs holding information on species, vocalization context, and animal behavior. After training on this dataset to connect language and audio representations, our model can identify over a thousand species' calls across taxa, complete bioacoustic tasks zero-shot, and retrieve animal vocalization recordings from natural text queries. When fine-tuned, BioLingual sets a new state-of-the-art on nine tasks in the Benchmark of Animal Sounds. Given its broad taxa coverage and ability to be flexibly queried in human language, we believe this model opens new paradigms in ecological monitoring and research, including free-text search on the world's acoustic monitoring archives. We open-source our models, dataset, and code. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Recommended Implementation of Quantitative Susceptibility Mapping for Clinical Research in The Brain: A Consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group

Authors: Berkin Bilgic, Mauro Costagli, Kwok-Shing Chan, Jeff Duyn, Christian Langkammer, Jongho Lee, Xu Li, Chunlei Liu, José P. Marques, Carlos Milovic, Simon Daniel Robinson, Ferdinand Schweser, Karin Shmueli, Pascal Spincemaille, Sina Straub, Peter van Zijl, Yi Wang, ISMRM Electro-Magnetic Tissue Properties Study Group

Abstract: This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic suscepti… ▽ More This article provides recommendations for implementing quantitative susceptibility mapping (QSM) for clinical brain research. It is a consensus of the ISMRM Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available give rise to the need in the neuroimaging community for guidelines on implementation. This article describes relevant considerations and provides specific implementation recommendations for all steps in QSM data acquisition, processing, analysis, and presentation in scientific publications. We recommend that data be acquired using a monopolar 3D multi-echo GRE sequence, that phase images be saved and exported in DICOM format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields should be removed within the brain mask using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of whole brain as a region of interest in the analysis, and QSM results should be reported with - as a minimum - the acquisition and processing specifications listed in the last section of the article. These recommendations should facilitate clinical QSM research and lead to increased harmonization in data acquisition, analysis, and reporting. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: Authors are listed in alphabetical order. * indicates corresponding authors

Scattering Transparency of Clouds in Exoplanet Transit Spectra

Authors: Bhavesh Jaiswal, Tyler D. Robinson

Abstract: The presence of aerosols in an exoplanet atmosphere can veil the underlying material and can lead to a flat transmission spectrum during primary transit observations. In this work, we explore forward scattering effects from super-micron sized aerosol particles present in the atmosphere of a transiting exoplanet. We find that the impacts of forward scattering from larger aerosols can significantly… ▽ More The presence of aerosols in an exoplanet atmosphere can veil the underlying material and can lead to a flat transmission spectrum during primary transit observations. In this work, we explore forward scattering effects from super-micron sized aerosol particles present in the atmosphere of a transiting exoplanet. We find that the impacts of forward scattering from larger aerosols can significantly impact exoplanet transits and the strength of these effects can be dependent on wavelength. In certain cloud configurations, the forward-scattered light can effectively pass through the clouds unhindered, thus rendering the clouds transparent. The dependence of the aerosol scattering properties on wavelength can then lead to a positive slope in the transit spectrum. These slopes are characteristically different from both Rayleigh and aerosol absorption slopes. As examples, we demonstrate scattering effects for both a rocky world and a hot Jupiter. In these models, the predicted spectral slopes due to forward scattering effects can manifest in the transit spectrum at the level of $\sim$10s to $\sim$100s of parts per million and, hence, could be observable with NASA's James Webb Space Telescope. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 9 pages, 7 figures, published in MNRAS

Plasma flows during the ablation stage of an over-massed pulsed-power-driven exploding planar wire array

Authors: R. Datta, J. Angel, J. B. Greenly, S. N. Bland, J. P. Chittenden, E. S. Lavine, W. M. Potter, D. Robinson, T. W. O. Varnish, E. Wong, D. A. Hammer, B. R. Kusse, J. D. Hare

Abstract: We characterize the plasma flows generated during the ablation stage of an over-massed exploding planar wire array, fielded on the COBRA pulsed-power facility (1 MA peak current, 250 ns rise time). The planar wire array is designed to provide a driving magnetic field (80-100 T) and current per wire distribution (about 60 kA), similar to that in a 10 MA cylindrical exploding wire array fielded on t… ▽ More We characterize the plasma flows generated during the ablation stage of an over-massed exploding planar wire array, fielded on the COBRA pulsed-power facility (1 MA peak current, 250 ns rise time). The planar wire array is designed to provide a driving magnetic field (80-100 T) and current per wire distribution (about 60 kA), similar to that in a 10 MA cylindrical exploding wire array fielded on the Z machine. Over-massing the arrays enables continuous plasma ablation over the duration of the experiment. The requirement to over-mass on the Z machine necessitates wires with diameters of 75-100 $μ$m, which are thicker than wires usually fielded on wire array experiments. To test ablation with thicker wires, we perform a parametric study by varying the initial wire diameter between 33-100 $μ$m. The largest wire diameter (100 $μ$m) array exhibits early closure of the AK gap, while the gap remains open during the duration of the experiment for wire diameters between 33-75 $μ$m. Laser plasma interferometry and time-gated XUV imaging are used to probe the plasma flows ablating from the wires. The plasma flows from the wires converge to generate a pinch, which appears as a fast-moving ($V \approx {100}$ kms$^{-1}$) column of increased plasma density ($\bar{n}_e \approx 2 \times 10^{18}$ cm$^{-3}$) and strong XUV emission. Finally, we compare the results with three-dimensional resistive-magnetohydrodynamic (MHD) simulations performed using the code GORGON, the results of which reproduce the dynamics of the experiment reasonably well. △ Less

Submitted 12 July, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: 14 pages; 14 figures

Journal ref: Phys. Plasmas 30, 092104 (2023)

A Stochastic-Gradient-based Interior-Point Algorithm for Solving Smooth Bound-Constrained Optimization Problems

Authors: Frank E. Curtis, Vyacheslav Kungurtsev, Daniel P. Robinson, Qi Wang

Abstract: A stochastic-gradient-based interior-point algorithm for minimizing a continuously differentiable objective function (that may be nonconvex) subject to bound constraints is presented, analyzed, and demonstrated through experimental results. The algorithm is unique from other interior-point methods for solving smooth nonconvex optimization problems since the search directions are computed using sto… ▽ More A stochastic-gradient-based interior-point algorithm for minimizing a continuously differentiable objective function (that may be nonconvex) subject to bound constraints is presented, analyzed, and demonstrated through experimental results. The algorithm is unique from other interior-point methods for solving smooth nonconvex optimization problems since the search directions are computed using stochastic gradient estimates. It is also unique in its use of inner neighborhoods of the feasible region -- defined by a positive and vanishing neighborhood-parameter sequence -- in which the iterates are forced to remain. It is shown that with a careful balance between the barrier, step-size, and neighborhood sequences, the proposed algorithm satisfies convergence guarantees in both deterministic and stochastic settings. The results of numerical experiments show that in both settings the algorithm can outperform projection-based methods. △ Less

Submitted 13 March, 2024; v1 submitted 28 April, 2023; originally announced April 2023.

Imaging 3D Chemistry at 1 nm Resolution with Fused Multi-Modal Electron Tomography

Authors: Jonathan Schwartz, Zichao Wendy Di, Yi Jiang, Jason Manassa, Jacob Pietryga, Yiwen Qian, Min Gee Cho, Jonathan L. Rowell, Huihuo Zheng, Richard D. Robinson, Junsi Gu, Alexey Kirilin, Steve Rozeveld, Peter Ercius, Jeffrey A. Fessler, Ting Xu, Mary Scott, Robert Hovden

Abstract: Measuring the three-dimensional (3D) distribution of chemistry in nanoscale matter is a longstanding challenge for metrological science. The inelastic scattering events required for 3D chemical imaging are too rare, requiring high beam exposure that destroys the specimen before an experiment completes. Even larger doses are required to achieve high resolution. Thus, chemical mapping in 3D has been… ▽ More Measuring the three-dimensional (3D) distribution of chemistry in nanoscale matter is a longstanding challenge for metrological science. The inelastic scattering events required for 3D chemical imaging are too rare, requiring high beam exposure that destroys the specimen before an experiment completes. Even larger doses are required to achieve high resolution. Thus, chemical mapping in 3D has been unachievable except at lower resolution with the most radiation-hard materials. Here, high-resolution 3D chemical imaging is achieved near or below one nanometer resolution in a Au-Fe$_3$O$_4$ metamaterial, Co$_3$O$_4$ - Mn$_3$O$_4$ core-shell nanocrystals, and ZnS-Cu$_{0.64}$S$_{0.36}$ nanomaterial using fused multi-modal electron tomography. Multi-modal data fusion enables high-resolution chemical tomography often with 99\% less dose by linking information encoded within both elastic (HAADF) and inelastic (EDX / EELS) signals. Now sub-nanometer 3D resolution of chemistry is measurable for a broad class of geometrically and compositionally complex materials. △ Less

Submitted 18 June, 2024; v1 submitted 24 April, 2023; originally announced April 2023.

Journal ref: Nat Commun 15, 3555 (2024)

Predicting dynamic, motion-related changes in B0 field in the brain at a 7 T MRI using a subject-specific fine-tuned U-net

Authors: Stanislav Motyka, Paul Weiser, Beata Bachrata, Lukas Hingerl, Bernhard Strasser, Gilbert Hangel, Eva Niess, Dario Goranovic, Fabian Niess, Maxim Zaitsev, Simon Daniel Robinson, Georg Langs, Siegfried Trattnig, Wolfgang Bogner

Abstract: Subject movement during the magnetic resonance examination is inevitable and causes not only image artefacts but also deteriorates the homogeneity of the main magnetic field (B0), which is a prerequisite for high quality data. Thus, characterization of changes to B0, e.g. induced by patient movement, is important for MR applications that are prone to B0 inhomogeneities. We propose a deep learning… ▽ More Subject movement during the magnetic resonance examination is inevitable and causes not only image artefacts but also deteriorates the homogeneity of the main magnetic field (B0), which is a prerequisite for high quality data. Thus, characterization of changes to B0, e.g. induced by patient movement, is important for MR applications that are prone to B0 inhomogeneities. We propose a deep learning based method to predict such changes within the brain from the change of the head position to facilitate retrospective or even real-time correction. A 3D U-net was trained on in vivo brain 7T MRI data. The input consisted of B0 maps and anatomical images at an initial position, and anatomical images at a different head position (obtained by applying a rigid-body transformation on the initial anatomical image). The output consisted of B0 maps at the new head positions. We further fine-tuned the network weights to each subject by measuring a limited number of head positions of the given subject, and trained the U-net with these data. Our approach was compared to established dynamic B0 field mapping via interleaved navigators, which suffer from limited spatial resolution and the need for undesirable sequence modifications. Qualitative and quantitative comparison showed similar performance between an interleaved navigator-equivalent method and proposed method. We therefore conclude that it is feasible to predict B0 maps from rigid subject movement and, when combined with external tracking hardware, this information could be used to improve the quality of magnetic resonance acquisitions without the use of navigators. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 9 pages, 6 figures

RESPECT: Reinforcement Learning based Edge Scheduling on Pipelined Coral Edge TPUs

Authors: Jiaqi Yin, Yingjie Li, Daniel Robinson, Cunxi Yu

Abstract: Deep neural networks (DNNs) have substantial computational and memory requirements, and the compilation of its computational graphs has a great impact on the performance of resource-constrained (e.g., computation, I/O, and memory-bound) edge computing systems. While efficient execution of their computational graph requires an effective scheduling algorithm, generating the optimal scheduling soluti… ▽ More Deep neural networks (DNNs) have substantial computational and memory requirements, and the compilation of its computational graphs has a great impact on the performance of resource-constrained (e.g., computation, I/O, and memory-bound) edge computing systems. While efficient execution of their computational graph requires an effective scheduling algorithm, generating the optimal scheduling solution is a challenging NP-hard problem. Furthermore, the complexity of scheduling DNN computational graphs will further increase on pipelined multi-core systems considering memory communication cost, as well as the increasing size of DNNs. Using the synthetic graph for the training dataset, this work presents a reinforcement learning (RL) based scheduling framework RESPECT, which learns the behaviors of optimal optimization algorithms and generates near-optimal scheduling results with short solving runtime overhead. Our framework has demonstrated up to $\sim2.5\times$ real-world on-chip inference runtime speedups over the commercial compiler with ten popular ImageNet models deployed on the physical Coral Edge TPUs system. Moreover, compared to the exact optimization methods, the proposed RL scheduling improves the scheduling optimization runtime by up to 683$\times$ speedups compared to the commercial compiler and matches the exact optimal solutions with up to 930$\times$ speedups. Finally, we perform a comprehensive generalizability test, which demonstrates RESPECT successfully imitates optimal solving behaviors from small synthetic graphs to large real-world DNNs computational graphs. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: 6 pages, ACM/IEEE Design Automation Conference (DAC'23)

Deep Learning Pipeline for Preprocessing and Segmenting Cardiac Magnetic Resonance of Single Ventricle Patients from an Image Registry

Authors: Tina Yao, Nicole St. Clair, Gabriel F. Miller, Adam L. Dorfman, Mark A. Fogel, Sunil Ghelani, Rajesh Krishnamurthy, Christopher Z. Lam, Joshua D. Robinson, David Schidlow, Timothy C. Slesnick, Justin Weigand, Michael Quail, Rahul Rathod, Jennifer A. Steeden, Vivek Muthurangu

Abstract: Purpose: To develop and evaluate an end-to-end deep learning pipeline for segmentation and analysis of cardiac magnetic resonance images to provide core-lab processing for a multi-centre registry of Fontan patients. Materials and Methods: This retrospective study used training (n = 175), validation (n = 25) and testing (n = 50) cardiac magnetic resonance image exams collected from 13 institution… ▽ More Purpose: To develop and evaluate an end-to-end deep learning pipeline for segmentation and analysis of cardiac magnetic resonance images to provide core-lab processing for a multi-centre registry of Fontan patients. Materials and Methods: This retrospective study used training (n = 175), validation (n = 25) and testing (n = 50) cardiac magnetic resonance image exams collected from 13 institutions in the UK, US and Canada. The data was used to train and evaluate a pipeline containing three deep-learning models. The pipeline's performance was assessed on the Dice and IoU score between the automated and reference standard manual segmentation. Cardiac function values were calculated from both the automated and manual segmentation and evaluated using Bland-Altman analysis and paired t-tests. The overall pipeline was further evaluated qualitatively on 475 unseen patient exams. Results: For the 50 testing dataset, the pipeline achieved a median Dice score of 0.91 (0.89-0.94) for end-diastolic volume, 0.86 (0.82-0.89) for end-systolic volume, and 0.74 (0.70-0.77) for myocardial mass. The deep learning-derived end-diastolic volume, end-systolic volume, myocardial mass, stroke volume and ejection fraction had no statistical difference compared to the same values derived from manual segmentation with p values all greater than 0.05. For the 475 unseen patient exams, the pipeline achieved 68% adequate segmentation in both systole and diastole, 26% needed minor adjustments in either systole or diastole, 5% needed major adjustments, and the cropping model only failed in 0.4%. Conclusion: Deep learning pipeline can provide standardised 'core-lab' segmentation for Fontan patients. This pipeline can now be applied to the >4500 cardiac magnetic resonance exams currently in the FORCE registry as well as any new patients that are recruited. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 17 pages, 6 figures

Anomaly Search Over Many Sequences With Switching Costs

Authors: Matthew Ubl, Benjamin D. Robinson, Matthew T. Hale

Abstract: This paper considers the quickest search problem to identify anomalies among large numbers of data streams. These streams can model, for example, disjoint regions monitored by a mobile robot. A particular challenge is a version of the problem in which the experimenter must suffer a cost each time the data stream being sampled changes, such as the time the robot must spend moving between regions. I… ▽ More This paper considers the quickest search problem to identify anomalies among large numbers of data streams. These streams can model, for example, disjoint regions monitored by a mobile robot. A particular challenge is a version of the problem in which the experimenter must suffer a cost each time the data stream being sampled changes, such as the time the robot must spend moving between regions. In this paper, we propose an algorithm which accounts for switching costs by varying a confidence threshold that governs when the algorithm switches to a new data stream. Our main contributions are easily computable approximations for both the optimal value of this threshold and the optimal value of the parameter that determines when a stream must be re-sampled. Further, we empirically show (i) a uniform improvement for switching costs of interest and (ii) roughly equivalent performance for small switching costs when comparing to the closest available algorithm. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 6 pages, 4 figures

Sequential Quadratic Optimization for Stochastic Optimization with Deterministic Nonlinear Inequality and Equality Constraints

Authors: Frank E. Curtis, Daniel P. Robinson, Baoyu Zhou

Abstract: A sequential quadratic optimization algorithm for minimizing an objective function defined by an expectation subject to nonlinear inequality and equality constraints is proposed, analyzed, and tested. The context of interest is when it is tractable to evaluate constraint function and derivative values in each iteration, but it is intractable to evaluate the objective function or its derivatives in… ▽ More A sequential quadratic optimization algorithm for minimizing an objective function defined by an expectation subject to nonlinear inequality and equality constraints is proposed, analyzed, and tested. The context of interest is when it is tractable to evaluate constraint function and derivative values in each iteration, but it is intractable to evaluate the objective function or its derivatives in any iteration, and instead an algorithm can only make use of stochastic objective gradient estimates. Under loose assumptions, including that the gradient estimates are unbiased, the algorithm is proved to possess convergence guarantees in expectation. The results of numerical experiments are presented to demonstrate that the proposed algorithm can outperform an alternative approach that relies on the ability to compute more accurate gradient estimates. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 27 pages, 1 figure

Report number: Report-no: 23T-017 MSC Class: 49M05; 49M37; 65K05; 65K10; 90C15; 90C30; 90C55

The Local Ledoit-Peche Law

Authors: Van Latimer, Benjamin D. Robinson

Abstract: Ledoit and Peche proved convergence of certain functions of a random covariance matrix's resolvent; we refer to this as the Ledoit-Peche law. One important application of their result is shrinkage covariance estimation with respect to so-called Minimum Variance (MV) loss, discussed in the work of Ledoit and Wolf. We provide an essentially optimal rate of convergence and hypothesize it to be the sm… ▽ More Ledoit and Peche proved convergence of certain functions of a random covariance matrix's resolvent; we refer to this as the Ledoit-Peche law. One important application of their result is shrinkage covariance estimation with respect to so-called Minimum Variance (MV) loss, discussed in the work of Ledoit and Wolf. We provide an essentially optimal rate of convergence and hypothesize it to be the smallest possible rate of excess MV loss within the shrinkage class. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 12 pages

MSC Class: 60B20 ACM Class: G.3

A Variance-Reduced and Stabilized Proximal Stochastic Gradient Method with Support Identification Guarantees for Structured Optimization

Authors: Yutong Dai, Guanyi Wang, Frank E. Curtis, Daniel P. Robinson

Abstract: This paper introduces a new proximal stochastic gradient method with variance reduction and stabilization for minimizing the sum of a convex stochastic function and a group sparsity-inducing regularization function. Since the method may be viewed as a stabilized version of the recently proposed algorithm PStorm, we call our algorithm S-PStorm. Our analysis shows that S-PStorm has strong convergenc… ▽ More This paper introduces a new proximal stochastic gradient method with variance reduction and stabilization for minimizing the sum of a convex stochastic function and a group sparsity-inducing regularization function. Since the method may be viewed as a stabilized version of the recently proposed algorithm PStorm, we call our algorithm S-PStorm. Our analysis shows that S-PStorm has strong convergence results. In particular, we prove an upper bound on the number of iterations required by S-PStorm before its iterates correctly identify (with high probability) an optimal support (i.e., the zero and nonzero structure of an optimal solution). Most algorithms in the literature with such a support identification property use variance reduction techniques that require either periodically evaluating an exact gradient or storing a history of stochastic gradients. Unlike these methods, S-PStorm achieves variance reduction without requiring either of these, which is advantageous. Moreover, our support-identification result for S-PStorm shows that, with high probability, an optimal support will be identified correctly in all iterations with the index above a threshold. We believe that this type of result is new to the literature since the few existing other results prove that the optimal support is identified with high probability at each iteration with a sufficiently large index (meaning that the optimal support might be identified in some iterations, but not in others). Numerical experiments on regularized logistic loss problems show that S-PStorm outperforms existing methods in various metrics that measure how efficiently and robustly iterates of an algorithm identify an optimal support. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: The work is accepted for presentation at AISTATS 2023. This is a technical report version

MSC Class: 49M29; 49M37; 65K05; 65K10; 65Y20; 68Q25; 90C30; 90C60

Autonomy Is An Acquired Taste: Exploring Developer Preferences for GitHub Bots

Authors: Amir Ghorbani, Nathan Cassee, Derek Robinson, Adam Alami, Neil A. Ernst, Alexander Serebrenik, Andrzej Wasowski

Abstract: Software bots fulfill an important role in collective software development, and their adoption by developers promises increased productivity. Past research has identified that bots that communicate too often can irritate developers, which affects the utility of the bot. However, it is not clear what other properties of human-bot collaboration affect developers' preferences, or what impact these pr… ▽ More Software bots fulfill an important role in collective software development, and their adoption by developers promises increased productivity. Past research has identified that bots that communicate too often can irritate developers, which affects the utility of the bot. However, it is not clear what other properties of human-bot collaboration affect developers' preferences, or what impact these properties might have. The main idea of this paper is to explore characteristics affecting developer preferences for interactions between humans and bots, in the context of GitHub pull requests. We carried out an exploratory sequential study with interviews and a subsequent vignette-based survey. We find developers generally prefer bots that are personable but show little autonomy, however, more experienced developers tend to prefer more autonomous bots. Based on this empirical evidence, we recommend bot developers increase configuration options for bots so that individual developers and projects can configure bots to best align with their own preferences and project cultures. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: appears in proceedings of ICSE, 2023

Exploring the $τ$ polarization in $B\to Xτ\barν$ along different axes

Authors: Florian U. Bernlochner, Zoltan Ligeti, Michele Papucci, Dean J. Robinson

Abstract: The $τ$ polarization in semileptonic $B$ decays provides probes of new physics complementary to decay rate distributions of the three-body final state. Prior calculations for inclusive decays used a definition for the polarization axis that is different from the choice used in calculations (and the only measurement) for exclusive channels. To compare inclusive and exclusive predictions, we calcula… ▽ More The $τ$ polarization in semileptonic $B$ decays provides probes of new physics complementary to decay rate distributions of the three-body final state. Prior calculations for inclusive decays used a definition for the polarization axis that is different from the choice used in calculations (and the only measurement) for exclusive channels. To compare inclusive and exclusive predictions, we calculate the $τ$ polarization in inclusive $B\to Xτ\barν$ using the same choice as in the exclusive decays, and construct a sum rule relating the inclusive $τ$ polarization to a weighted sum of exclusive decay polarizations. We use this relation, experimental data, and theoretical predictions for the decays to the lightest charm or up-type hadrons to make predictions for excited channels. △ Less

Submitted 30 March, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Comments: 6 pages, 2 figures

Integrated distance sampling models for simple point counts

Authors: Marc Kery, J. Andrew Royle, Tyler Hallman, W. Douglas Robinson, Nicolas Strebel, Kenneth F. Kellner

Abstract: Point counts (PCs) are widely used in biodiversity surveys, but despite numerous advantages, simple PCs suffer from several problems: detectability, and therefore abundance, is unknown; systematic spatiotemporal variation in detectability produces biased inferences, and unknown survey area prevents formal density estimation and scaling-up to the landscape level. We introduce integrated distance sa… ▽ More Point counts (PCs) are widely used in biodiversity surveys, but despite numerous advantages, simple PCs suffer from several problems: detectability, and therefore abundance, is unknown; systematic spatiotemporal variation in detectability produces biased inferences, and unknown survey area prevents formal density estimation and scaling-up to the landscape level. We introduce integrated distance sampling (IDS) models that combine distance sampling (DS) with simple PC or detection/nondetection (DND) data and capitalize on the strengths and mitigate the weaknesses of each data type. Key to IDS models is the view of simple PC and DND data as aggregations of latent DS surveys that observe the same underlying density process. This enables estimation of separate detection functions, along with distinct covariate effects, for all data types. Additional information from repeat or time-removal surveys, or variable survey duration, enables separate estimation of the availability and perceptibility components of detectability. IDS models reconcile spatial and temporal mismatches among data sets and solve the above-mentioned problems of simple PC and DND data. To fit IDS models, we provide JAGS code and the new IDS() function in the R package unmarked. Extant citizen-science data generally lack adjustments for detection biases, but IDS models address this shortcoming, thus greatly extending the utility and reach of these data. In addition, they enable formal density estimation in hybrid designs, which efficiently combine distance sampling with distance-free, point-based PC or DND surveys. We believe that IDS models have considerable scope in ecology, management, and monitoring. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 24 pages, plus 5 figures, plus 2 Appendices

Geometry Optimization for Long-lived Particle Detectors

Authors: Thomas Gorordo, Simon Knapen, Benjamin Nachman, Dean J. Robinson, Adi Suresh

Abstract: The proposed designs of many auxiliary long-lived particle (LLP) detectors at the LHC call for the instrumentation of a large surface area inside the detector volume, in order to reliably reconstruct tracks and LLP decay vertices. Taking the CODEX-b detector as an example, we provide a proof-of-concept optimization analysis that demonstrates the required instrumented surface area can be substantia… ▽ More The proposed designs of many auxiliary long-lived particle (LLP) detectors at the LHC call for the instrumentation of a large surface area inside the detector volume, in order to reliably reconstruct tracks and LLP decay vertices. Taking the CODEX-b detector as an example, we provide a proof-of-concept optimization analysis that demonstrates the required instrumented surface area can be substantially reduced for many LLP models, while only marginally affecting the LLP signal efficiency. This optimization permits a significant reduction in cost and installation time, and may also inform the installation order for modular detector elements. We derive a branch-and-bound based optimization algorithm that permits highly computationally efficient determination of optimal detector configurations, subject to any specified LLP vertex and track reconstruction requirements. We outline the features of a newly-developed generalized simulation framework, for the computation of LLP signal efficiencies across a range of LLP models and detector geometries. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 46 pages, 11 figures, 3 tables

Inexact Proximal-Gradient Methods with Support Identification

Authors: Yutong Dai, Daniel P. Robinson

Abstract: We consider the proximal-gradient method for minimizing an objective function that is the sum of a smooth function and a non-smooth convex function. A feature that distinguishes our work from most in the literature is that we assume that the associated proximal operator does not admit a closed-form solution. To address this challenge, we study two adaptive and implementable termination conditions… ▽ More We consider the proximal-gradient method for minimizing an objective function that is the sum of a smooth function and a non-smooth convex function. A feature that distinguishes our work from most in the literature is that we assume that the associated proximal operator does not admit a closed-form solution. To address this challenge, we study two adaptive and implementable termination conditions that dictate how accurately the proximal-gradient subproblem is solved. We prove that the number of iterations required for the inexact proximal-gradient method to reach a $τ> 0$ approximate first-order stationary point is $\mathcal{O}(τ^{-2})$, which matches the similar result that holds when exact subproblem solutions are computed. Also, by focusing on the overlapping group $\ell_1$ regularizer, we propose an algorithm for approximately solving the proximal-gradient subproblem, and then prove that its iterates identify (asymptotically) the support of an optimal solution. If one imposes additional control over the accuracy to which each subproblem is solved, we give an upper bound on the maximum number of iterations before the support of an optimal solution is obtained. △ Less

Submitted 3 November, 2022; originally announced November 2022.

MSC Class: 49M29; 49M37; 65K05; 65K10; 65Y20; 68Q25; 90C30; 90C60

A Radiative-Convective Model for Terrestrial Planets with Self-Consistent Patchy Clouds

Authors: James D. Windsor, Tyler D. Robinson, Ravi kumar Kopparapu, David E. Trilling, Joe LLama, Amber Young

Abstract: Clouds are ubiquitous\, -- \,they arise for every solar system planet that possesses an atmosphere and have also been suggested as a leading mechanism for obscuring spectral features in exoplanet observations. As exoplanet observations continue to improve, there is a need for efficient and general planetary climate models that appropriately handle the possible cloudy atmospheric environments that… ▽ More Clouds are ubiquitous\, -- \,they arise for every solar system planet that possesses an atmosphere and have also been suggested as a leading mechanism for obscuring spectral features in exoplanet observations. As exoplanet observations continue to improve, there is a need for efficient and general planetary climate models that appropriately handle the possible cloudy atmospheric environments that arise on these worlds. We generate a new 1D radiative-convective terrestrial planet climate model that self-consistently handles patchy clouds through a parameterized microphysical treatment of condensation and sedimentation processes. Our model is general enough to recreate Earth's atmospheric radiative environment without over-parameterization, while also maintaining a simple implementation that is applicable to a wide range of atmospheric compositions and physical planetary properties. We first validate this new 1D patchy cloud radiative-convective climate model by comparing it to Earth thermal structure data and to existing climate and radiative transfer tools. We produce partially-clouded Earth-like climates with cloud structures that are representative of deep tropospheric convection and are adequate 1D representations of clouds within rocky planet atmospheres. After validation against Earth, we then use our partially clouded climate model and explore the potential climates of super-Earth exoplanets with secondary nitrogen-dominated atmospheres which we assume are abiotic. We also couple the partially clouded climate model to a full-physics, line-by-line radiative transfer model and generate high-resolution spectra of simulated climates. These self-consistent climate-to-spectral models bridge the gap between climate modeling efforts and observational studies of rocky worlds. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 28 pages, 15 figures, submitted; community comments welcome

Tunable intervalence charge transfer in ruthenium Prussian blue analogue enables stable and efficient biocompatible artificial synapses

Authors: Donald A. Robinson, Michael E. Foster, Christopher H. Bennett, Austin Bhandarkar, Elizabeth R. Webster, Aleyna Celebi, Nisa Celebi, Elliot J. Fuller, Vitalie Stavila, Catalin D. Spataru, David S. Ashby, Matthew J. Marinella, Raga Krishnakumar, Mark D. Allendorf, A. Alec Talin

Abstract: Emerging concepts for neuromorphic computing, bioelectronics, and brain-computer interfacing inspire new research avenues aimed at understanding the relationship between oxidation state and conductivity in unexplored materials. Here, we present ruthenium Prussian blue analogue (RuPBA), a mixed valence coordination compound with an open framework structure and ability to conduct both ionic and elec… ▽ More Emerging concepts for neuromorphic computing, bioelectronics, and brain-computer interfacing inspire new research avenues aimed at understanding the relationship between oxidation state and conductivity in unexplored materials. Here, we present ruthenium Prussian blue analogue (RuPBA), a mixed valence coordination compound with an open framework structure and ability to conduct both ionic and electronic charge, for flexible artificial synapses that reversibly switch conductance by more than four orders of magnitude based on electrochemically tunable oxidation state. Retention of programmed states is improved by nearly two orders of magnitude compared to the extensively studied organic polymers, thus reducing the frequency, complexity and energy costs associated with error correction schemes. We demonstrate dopamine detection using RuPBA synapses and biocompatibility with neuronal cells, evoking prospective application for brain-computer interfacing. By application of electron transfer theory to in-situ spectroscopic probing of intervalence charge transfer, we elucidate a switching mechanism whereby the degree of mixed valency between N-coordinated Ru sites controls the carrier concentration and mobility, as supported by DFT. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Interpreting LHCb's $Λ_b\to Λ_cτ\barν$ measurement and puzzles in semileptonic $Λ_b$ decays

Authors: Florian U. Bernlochner, Zoltan Ligeti, Michele Papucci, Dean J. Robinson

Abstract: Normalizing the recent LHCb measurement of $Λ_b \to Λ_c τ\barν$ to the standard model (SM) prediction for the $Λ_b \to Λ_c μ\barν$ rate, instead of a LEP measurement, provides a more consistent comparison with the SM prediction for the lepton flavor universality ratio $R(Λ_c)$. This modestly increases $R(Λ_c)$ compared to the quoted LHCb result, such that it no longer hints at a suppression compar… ▽ More Normalizing the recent LHCb measurement of $Λ_b \to Λ_c τ\barν$ to the standard model (SM) prediction for the $Λ_b \to Λ_c μ\barν$ rate, instead of a LEP measurement, provides a more consistent comparison with the SM prediction for the lepton flavor universality ratio $R(Λ_c)$. This modestly increases $R(Λ_c)$ compared to the quoted LHCb result, such that it no longer hints at a suppression compared to the SM, which would be hard to accommodate in new physics scenarios that enhance $R(D^{(*)})$. We point out that the fraction of excited states in inclusive semileptonic $Λ_b$ decay may be significantly greater than the corresponding fraction in $B$ decays. Possible implications are speculated upon. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: 3 pages

Report number: CALT-TH-2022-023