Showing 1–50 of 374 results for author: Hansen, S

Statistical inference of collision frequencies from x-ray Thomson scattering spectra

Authors: Thomas W. Hentschel, Alina Kononov, Andrew D. Baczewski, Stephanie B. Hansen

Abstract: Thomson scattering spectra measure the response of plasma particles to incident radiation. In warm dense matter, which is opaque to visible light, x-ray Thomson scattering (XRTS) enables a detailed probe of the electron distribution and has been used as a diagnostic for electron temperature, density, and plasma ionization. In this work, we examine the sensitivities of inelastic XRTS signatures to… ▽ More Thomson scattering spectra measure the response of plasma particles to incident radiation. In warm dense matter, which is opaque to visible light, x-ray Thomson scattering (XRTS) enables a detailed probe of the electron distribution and has been used as a diagnostic for electron temperature, density, and plasma ionization. In this work, we examine the sensitivities of inelastic XRTS signatures to modeling details including the dynamic collision frequency and the electronic density of states. Applying verified Monte Carlo inversion methods to dynamic structure factors obtained from time-dependent density functional theory, we assess the utility of XRTS signals as a way to inform the dynamic collision frequency, especially its direct-current (DC) limit, which is directly related to the electrical conductivity. △ Less

Submitted 27 August, 2024; originally announced August 2024.

Travel time and energy dissipation minima in heterogeneous subsurface flows

Authors: Scott K. Hansen, Daniel O'Malley

Abstract: We establish a number of results concerning conditions for minimum energy dissipation and advective travel time in porous and fractured media. First, we establish a pair of converse results concerning fluid motion along a streamline between two points of fixed head: the minimal advective time is achieved under conditions of constant energy dissipation, and minimal energy dissipation is achieved un… ▽ More We establish a number of results concerning conditions for minimum energy dissipation and advective travel time in porous and fractured media. First, we establish a pair of converse results concerning fluid motion along a streamline between two points of fixed head: the minimal advective time is achieved under conditions of constant energy dissipation, and minimal energy dissipation is achieved under conditions of constant velocity along the streamline (implying homogeneous conductivity in the vicinity of the streamline). We also show directly by means of variational methods that minimum advection time along a streamline with a given average conductivity is achieved when the conductivity is constant. Finally, we turn our attention to minimum advection time and energy dissipation in parallel and sequential fracture systems governed by the cubic law: for which fracture cross-section and conductivity are intimately linked. We show that, as in porous domains, flow partitioning between different pathways always acts to minimize system energy dissipation. Finally, we consider minimum advection time as a function of aperture distribution in a sequence of fracture segments. We show that, for a fixed average aperture, a uniform-aperture system displays the shortest advection time. However, we also show that any sufficiently small small perturbations in aperture away from uniformity always act to reduce advection time. △ Less

Submitted 25 July, 2024; originally announced July 2024.

Upscaling transport in heterogeneous media featuring local-scale dispersion: flow channeling, macro-retardation and parameter prediction

Authors: Lian Zhou, Scott K. Hansen

Abstract: Many theoretical treatments of transport in heterogeneous Darcy flows consider advection only. When local-scale dispersion is neglected, flux-weighting persists over time; mean Lagrangian and Eulerian flow velocity distributions relate simply to each other and to the variance of the underlying hydraulic conductivity field. Local-scale dispersion complicates this relationship, potentially causing i… ▽ More Many theoretical treatments of transport in heterogeneous Darcy flows consider advection only. When local-scale dispersion is neglected, flux-weighting persists over time; mean Lagrangian and Eulerian flow velocity distributions relate simply to each other and to the variance of the underlying hydraulic conductivity field. Local-scale dispersion complicates this relationship, potentially causing initially flux-weighted solute to experience lower-velocity regions as well as Taylor-type macrodispersion due to transverse solute movement between adjacent streamlines. To investigate the interplay of local-scale dispersion with conductivity log-variance, correlation length, and anisotropy, we perform a Monte Carlo study of flow and advective-dispersive transport in spatially-periodic 2D Darcy flows in large-scale, high-resolution multivariate Gaussian random conductivity fields. We observe flow channeling at all heterogeneity levels and quantify its extent. We find evidence for substantial effective retardation in the upscaled system, associated with increased flow channeling, and observe limited Taylor-type macrodispersion, which we physically explain. A quasi-constant Lagrangian velocity is achieved within a short distance of release, allowing usage of a simplified continuous-time random walk (CTRW) model we previously proposed in which the transition time distribution is understood as a temporal mapping of unit time in an equivalent system with no flow heterogeneity. The numerical data set is modeled with such a CTRW; we show how dimensionless parameters defining the CTRW transition time distribution are predicted by dimensionless heterogeneity statistics and provide empirical equations for this purpose. △ Less

Submitted 24 July, 2024; originally announced July 2024.

PyDDC: An Eulerian-Lagrangian simulator for density driven convection of $\mathrm{CO_2}$--brine systems in saturated porous media

Authors: Sayan Sen, Scott K. Hansen

Abstract: PyDDC is a particle tracking reservoir simulator capable of solving non-linear density driven convection of single phase carbon-dioxide ($\mathrm{CO_2}$)--brine fluid mixture in saturated porous media at the continuum scale. In contrast to the sate-of-the-art Eulerian models, PyDDC uses a Lagrangian approach to simulate the Fickian transport of single phase solute mixtures. This introduces additio… ▽ More PyDDC is a particle tracking reservoir simulator capable of solving non-linear density driven convection of single phase carbon-dioxide ($\mathrm{CO_2}$)--brine fluid mixture in saturated porous media at the continuum scale. In contrast to the sate-of-the-art Eulerian models, PyDDC uses a Lagrangian approach to simulate the Fickian transport of single phase solute mixtures. This introduces additional flexibility of incorporating anisotropic dispersion and benefits from having no numerical artifacts in its implementation. It also includes $\mathrm{CO_2}$--brine phase equilibrium models, developed by other researchers, to study the overall dynamics in the presence of electrolyte brine at different pressure and temperatures above the critical point of $\mathrm{CO_2}$. We demonstrate the implementation procedure in depth, outlining the overall structure of the numerical solver and its different attributes that can be used for solving specific tasks. △ Less

Submitted 20 July, 2024; originally announced July 2024.

Once a Triple, Not Always a Triple: The Evolution of Hierarchical Triples that Yield Merged Inner Binaries

Authors: Cheyanne Shariat, Smadar Naoz, Kareem El-Badry, Antonio C. Rodriguez, Bradley M. S. Hansen, Isabel Angelo, Alexander P. Stephan

Abstract: More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and… ▽ More More than half of all main-sequence (MS) stars have one or more companions, and many of those with initial masses <8 M$_\odot$ are born in hierarchical triples. These systems feature two stars in a close orbit (the inner binary) while a tertiary star orbits them on a wider orbit (the outer binary). In hierarchical triples, three-body dynamics combined with stellar evolution drives interactions and, in many cases, merges the inner binary entirely to create a renovated `Post-Merger Binary' (PMB). By leveraging dynamical simulations and tracking binary interactions, we explore the outcomes of merged triples and investigate whether PMBs preserve signatures of their three-body history. Our findings indicate that in 26-54% of wide double WD binaries (s>100 au), the more massive white dwarf (WD) is a merger product, implying that these DWD binaries were previously triples. Overall, we estimate that $44\pm14\%$ of observed wide DWDs originated in triple star systems and thereby have rich dynamical histories. Additionally, our results suggest that the separations of inner and outer binaries are uncorrelated at birth, providing insights into stellar formation processes. We also examine MS+MS and MS+Red Giant mergers manifesting as Blue Straggler stars (BSSs). These PMBs have orbital configurations and ages similar to most observed BSS binaries. While the triple+merger formation channel can explain the observed chemical abundances, moderate eccentricities, and companion masses in BSS binaries, it likely only accounts for $\sim$20-25% of BSSs. Meanwhile, we predict that the majority of observed single BSSs formed as collisions in triples and harbor long-period (>10 yr) companions. Furthermore, both BSS binaries and DWDs exhibit signatures of WD birth kicks. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ. Comments are welcome

Existence of a maximum flow rate in electro-osmotic systems

Authors: Sleeba Varghese, Billy D. Todd, Jesper S. Hansen

Abstract: In this work, we investigate the effect of the hydrodynamic wall-fluid friction in electro-osmotic flows. First, we present the solution to the electro-hydrodynamic equation for the electro-osmotic velocity profile, which is derived for an ionic system composed of cations immersed in uncharged solvent particles. The system (solution and walls) is kept electrically neutral using negatively charged… ▽ More In this work, we investigate the effect of the hydrodynamic wall-fluid friction in electro-osmotic flows. First, we present the solution to the electro-hydrodynamic equation for the electro-osmotic velocity profile, which is derived for an ionic system composed of cations immersed in uncharged solvent particles. The system (solution and walls) is kept electrically neutral using negatively charged walls and will here be referred to as a "counterion-only" system. The theory predicts the existence of a counterion concentration that results in a maximum electro-osmotic flow rate, but only if the wall-fluid friction, or equivalently the slip length, is correlated with the system electrostatic screening length. Through equilibrium molecular dynamics simulations we precisely determine the hydrodynamic slip from the wall-fluid friction, and this is then used as input to the theoretical predictions. Comparison between the theory and independent non-equilibrium molecular dynamics simulation data confirms the existence of the maximum. Also, we find that standard hydrodynamic theory quantitatively agrees with simulation results for charged nanoscale systems for sufficiently small charge densities and ion charges, if the correct slip boundaries are applied. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: Main manuscript: 21 pages, 8 figures. Supporting information: 6 pages, 2 figures

The unpaved road towards efficient selective breeding in insects for food and feed

Authors: Laura Skrubbeltrang Hansen, Stine Frey Laursen, Simon Bahrndorff, Jesper Givskov Sørensen, Goutam Sahana, Torsten Nygaard Kristensen, Hanne Marie Nielsen

Abstract: Insect production for food and feed presents a promising supplement to ensure food safety and address the adverse impacts of agriculture on climate and environment in the future. However, optimisation is required for insect production to realise its full potential. This can be by targeted improvement of traits of interest through selective breeding, an approach which has so far been underexplored… ▽ More Insect production for food and feed presents a promising supplement to ensure food safety and address the adverse impacts of agriculture on climate and environment in the future. However, optimisation is required for insect production to realise its full potential. This can be by targeted improvement of traits of interest through selective breeding, an approach which has so far been underexplored and underutilised in insect farming. Here we present a comprehensive review of the selective breeding framework in the context of insect production. We systematically evaluate adjustments of selective breeding techniques to the realm of insects and highlight the essential components integral to the breeding process. The discussion covers every step of a conventional breeding scheme, such as formulation of breeding objectives, phenotyping, estimation of genetic parameters and breeding values, selection of appropriate breeding strategies, and mitigation of issues associated with genetic diversity depletion and inbreeding. This review combines knowledge from diverse disciplines, bridging the gap between animal breeding, quantitative genetics, evolutionary biology, and entomology, offering an integrated view of the insect breeding research area and uniting knowledge which has previously remained scattered across diverse fields of expertise. △ Less

Submitted 26 June, 2024; v1 submitted 24 June, 2024; originally announced June 2024.

Cosmological test of an ultraviolet origin of Dark Energy

Authors: Hans Christiansen, Bence Takács, Steen H. Hansen

Abstract: The accelerated expansion of the Universe is impressively well described by a cosmological constant. However, the observed value of the cosmological constant is much smaller than expected based on quantum field theories. Recent efforts to achieve consistency in these theories have proposed a relationship between Dark Energy and the most compact objects, such as black holes (BH). However, experimen… ▽ More The accelerated expansion of the Universe is impressively well described by a cosmological constant. However, the observed value of the cosmological constant is much smaller than expected based on quantum field theories. Recent efforts to achieve consistency in these theories have proposed a relationship between Dark Energy and the most compact objects, such as black holes (BH). However, experimental tests are very challenging to devise and perform. In this article, we present a testable model with no cosmological constant, in which the accelerated expansion can be driven by black holes. The model couples the expansion of the Universe (the Friedmann equation) with the mass-function of cosmological haloes (using the Press-Schechter formalism). Through the observed link between halo-masses and BH-masses one thus gets a coupling between the expansion rate of the Universe and the BHs. We compare the predictions of this simple BH model with SN1a data and find a poor agreement with observations. Our method is sufficiently general that it allows us to also test a fundamentally different model, also without a cosmological constant, where the accelerated expansion is driven by a new force proportional to the internal velocity dispersion of galaxies. Surprisingly enough this model cannot be excluded using the SN1a data. △ Less

Submitted 25 April, 2024; originally announced June 2024.

Comments: 11 pages, 2 figures, Accepted by Universe, special issue "The Nature of Dark Energy", https://www.mdpi.com/2218-1997/10/5/193

Analytical models for secular descents in hierarchical triple systems

Authors: Grant C. Weldon, Smadar Naoz, Bradley M. S. Hansen

Abstract: Triple body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit's eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the Eccentric Kozai-Lido… ▽ More Triple body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit's eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the Eccentric Kozai-Lidov (EKL) mechanism. The overall decrease in the inner orbit's pericenter distance has potentially dramatic effects for realistic systems, such as tidal disruption. We present an analytical approximation in the test particle limit to describe the step-wise eccentricity evolution of the inner orbit. We also integrate the equations of motion and present approximations in both the test particle limit and in the general case to describe the overall octupole-level time evolution of the eccentricity. The analytical approximations are compared to numerical simulations to show that the models accurately describe the form and timescale of the secular descent from large distances to a close-encounter distance (e.g., the Roche limit). By circumventing the need for numerical simulations to obtain the long-term behavior, these approximations can be used to readily estimate descent timescales for populations of systems. We demonstrate by calculating rates of EKL-driven migration for Hot Jupiters in stellar binaries. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures

Widespread disruption of resonant chains during protoplanetary disk dispersal

Authors: Brad M. S. Hansen, Tze-Yeung Yu, Yasuhiro Hasegawa

Abstract: We describe the evolution of low mass planets in a dispersing protoplanetary disk around a Solar mass star. The disk model is based on the results of Yu, Hansen & Hasegawa (2023), which describes a region of the inner disk where the direction of the migration torque is outwards due to the diffusion of the stellar magnetic field into the disk and the resultant gradual increase in surface density ou… ▽ More We describe the evolution of low mass planets in a dispersing protoplanetary disk around a Solar mass star. The disk model is based on the results of Yu, Hansen & Hasegawa (2023), which describes a region of the inner disk where the direction of the migration torque is outwards due to the diffusion of the stellar magnetic field into the disk and the resultant gradual increase in surface density outwards. We demonstrate that the magnetospheric rebound phase in such a disk leads to diverging orbits for double and triple planet systems, and the disruption of a high fraction of the initial resonant chains. We present simulations of three planet systems with masses based on the observed triple planet systems observed by the Kepler satellite within the context of this model. The final distribution of nearest neighbour period ratios provides an excellent fit to the observations, provided that the initial resonant configurations are compact. The occurrence rate of planets as a function of orbital period also provides a good match to the observations, for final orbital periods P<20 days. These results suggest that the period and period ratio distributions of low mass planets are primarily set in place during the disk dispersal epoch, and may not require significant dynamical evolution thereafter. △ Less

Submitted 18 July, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 17 pages, 18 figures, to appear in Open Journal of Astrophysics

First principles simulations of dense hydrogen

Authors: Michael Bonitz, Jan Vorberger, Mandy Bethkenhagen, Maximilian Böhme, David Ceperley, Alexey Filinov, Thomas Gawne, Frank Graziani, Gianluca Gregori, Paul Hamann, Stephanie Hansen, Markus Holzmann, S. X. Hu, Hanno Kählert, Valentin Karasiev, Uwe Kleinschmidt, Linda Kordts, Christopher Makait, Burkhard Militzer, Zhandos Moldabekov, Carlo Pierleoni, Martin Preising, Kushal Ramakrishna, Ronald Redmer, Sebastian Schwalbe , et al. (2 additional authors not shown)

Abstract: Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extre… ▽ More Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g. planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data. rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange-correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principle PIMC and DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Modelling the galaxy radio continuum from star formation and active galactic nuclei in the Shark semi-analytic model

Authors: Samuel P. Hansen, Claudia D. P. Lagos, Matteo Bonato, Robin H. W. Cook, Luke J. M. Davies, Ivan Delvecchio, Scott A. Tompkins

Abstract: We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGN) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission… ▽ More We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGN) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission is modelled based on the jet production rate, which depends on the black hole mass, accretion rate and spin, and includes synchrotron self-absorption. Shark reproduces radio luminosity functions (RLFs) at 1.4 GHz and 150 MHz for 0 $\leq$ z $\leq$ 4, and scaling relations between radio luminosity, star formation rate and infrared luminosity of galaxies in the local and distant universe in good agreement with observations. The model also reproduces observed number counts of radio sources from 150 MHz to 8.4 GHz to within a factor of two on average, though larger discrepancies are seen at the very bright fluxes at higher frequencies. We use this model to understand how the radio continuum emission from radio-quiet AGNs can affect the measured RLFs of galaxies. We find current methods to exclude AGNs from observational samples result in large fractions of radio-quiet AGNs contaminating the "star-forming galaxies" selection and a brighter end to the resulting RLFs. We investigate how this effects the infrared-radio correlation (IRRC) and show that AGN contamination can lead to evolution of the IRRC with redshift. Without this contamination our model predicts a redshift- and stellar mass-independent IRRC, except at the dwarf-galaxy regime. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: Accepted for publication in MNRAS. 17 pages, 7 figures

Productivity and quality-adjusted life years: QALYs, PALYs and beyond

Authors: Kristian S. Hansen, Juan D. Moreno-Ternero, Lars P. Østerdal

Abstract: We develop a unified framework for the measurement and valuation of health and productivity. Within this framework, we characterize evaluation functions allowing for compromises between the classical quality-adjusted life years (QALYs) and its polar productivity-adjusted life years (PALYs). Our framework and characterization results provide a new normative basis for the economic evaluation of heal… ▽ More We develop a unified framework for the measurement and valuation of health and productivity. Within this framework, we characterize evaluation functions allowing for compromises between the classical quality-adjusted life years (QALYs) and its polar productivity-adjusted life years (PALYs). Our framework and characterization results provide a new normative basis for the economic evaluation of health care interventions, as well as occupational health and safety policies, aimed to impact both health and productivity of individuals. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Gemini 1.5: Unlocking multimodal understanding across millions of tokens of context

Authors: Gemini Team, Petko Georgiev, Ving Ian Lei, Ryan Burnell, Libin Bai, Anmol Gulati, Garrett Tanzer, Damien Vincent, Zhufeng Pan, Shibo Wang, Soroosh Mariooryad, Yifan Ding, Xinyang Geng, Fred Alcober, Roy Frostig, Mark Omernick, Lexi Walker, Cosmin Paduraru, Christina Sorokin, Andrea Tacchetti, Colin Gaffney, Samira Daruki, Olcan Sercinoglu, Zach Gleicher, Juliette Love , et al. (1110 additional authors not shown)

Abstract: In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February… ▽ More In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February version on the great majority of capabilities and benchmarks; (2) Gemini 1.5 Flash, a more lightweight variant designed for efficiency with minimal regression in quality. Gemini 1.5 models achieve near-perfect recall on long-context retrieval tasks across modalities, improve the state-of-the-art in long-document QA, long-video QA and long-context ASR, and match or surpass Gemini 1.0 Ultra's state-of-the-art performance across a broad set of benchmarks. Studying the limits of Gemini 1.5's long-context ability, we find continued improvement in next-token prediction and near-perfect retrieval (>99%) up to at least 10M tokens, a generational leap over existing models such as Claude 3.0 (200k) and GPT-4 Turbo (128k). Finally, we highlight real-world use cases, such as Gemini 1.5 collaborating with professionals on completing their tasks achieving 26 to 75% time savings across 10 different job categories, as well as surprising new capabilities of large language models at the frontier; when given a grammar manual for Kalamang, a language with fewer than 200 speakers worldwide, the model learns to translate English to Kalamang at a similar level to a person who learned from the same content. △ Less

Submitted 8 August, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Inference for Regression with Variables Generated from Unstructured Data

Authors: Laura Battaglia, Timothy Christensen, Stephen Hansen, Szymon Sacher

Abstract: The leading strategy for analyzing unstructured data uses two steps. First, latent variables of economic interest are estimated with an upstream information retrieval model. Second, the estimates are treated as "data" in a downstream econometric model. We establish theoretical arguments for why this two-step strategy leads to biased inference in empirically plausible settings. More constructively,… ▽ More The leading strategy for analyzing unstructured data uses two steps. First, latent variables of economic interest are estimated with an upstream information retrieval model. Second, the estimates are treated as "data" in a downstream econometric model. We establish theoretical arguments for why this two-step strategy leads to biased inference in empirically plausible settings. More constructively, we propose a one-step strategy for valid inference that uses the upstream and downstream models jointly. The one-step strategy (i) substantially reduces bias in simulations; (ii) has quantitatively important effects in a leading application using CEO time-use data; and (iii) can be readily adapted by applied researchers. △ Less

Submitted 9 May, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Predicting the Dominant Formation Mechanism of Multi-Planetary Systems

Authors: Cheyanne Shariat, Yasuhiro Hasegawa, Bradley M. S. Hansen, Tze Yeung Mathew Yu, Renyu Hu

Abstract: Most, if not all, sun-like stars host one or more planets, making multi-planetary systems commonplace in our galaxy. We utilize hundreds of multi-planet simulations to explore the origin of such systems, focusing on their orbital architecture. The first set of simulations assumes in-situ assembly of planetary embryos, while the second explores planetary migration. After applying observational bias… ▽ More Most, if not all, sun-like stars host one or more planets, making multi-planetary systems commonplace in our galaxy. We utilize hundreds of multi-planet simulations to explore the origin of such systems, focusing on their orbital architecture. The first set of simulations assumes in-situ assembly of planetary embryos, while the second explores planetary migration. After applying observational biases to the simulations, we compare them to 250+ observed multi-planetary systems, including 13 systems with planets in the habitable zone. For all of the systems, we calculate two of the so-called statistical measures: the mass concentration ($S_{c}$) and orbital spacing ($S_{s}$). After analytic and empirical analyses, we find that the measures are related to first-order with a power law: $S_{c} \sim S_{s}^β$. The in-situ systems exhibit steeper power-law relations relative to the migration systems. We show that different formation scenarios cover different regions in the $S_{s} - S_{c}$ diagram with some overlap. Furthermore, we discover that observed systems with $S_{s} < 30$ are likely dominated by the migration scenario, while those with $S_{s} \geq 30$ are likely dominated by the in-situ scenario. We apply these criteria to determine that a majority (62%) of observed multi-planetary systems formed via migration, whereas most systems with currently observed habitable planets formed via in-situ assembly. This work provides methods of leveraging the statistical measures ($S_{s}$ and $S_{c}$) to disentangle the formation history of observed multi-planetary systems based on their present-day architectures. △ Less

Submitted 29 March, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJ Letters

Radiatively Cooled Magnetic Reconnection Experiments Driven by Pulsed Power

Authors: R Datta, K Chandler, C E Myers, J P Chittenden, A J Crilly, C Aragon, D J Ampleford, J T Banasek, A Edens, W R Fox, S B Hansen, E C Harding, C A Jennings, H Ji, C C Kuranz, S V Lebedev, Q Looker, S G Patel, A J Porwitzky, G A Shipley, D A Uzdensky, D A Yager-Elorriaga, J D Hare

Abstract: We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly-driven plasma flows ($M_A \approx 7$) with anti-parallel magnetic fields, which form a reconnection layer ($S_L \approx 120$) at the… ▽ More We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly-driven plasma flows ($M_A \approx 7$) with anti-parallel magnetic fields, which form a reconnection layer ($S_L \approx 120$) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate ($τ_{\text{cool}}^{-1}/τ_{\text{A}}^{-1} > 100$), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in X-ray emission from the reconnection layer, measured using filtered diodes and time-gated X-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong X-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer ($\leq$ 75 eV) and inflow temperatures (about 2 eV), and that these hotspots generate the majority of the high-energy (> 1 keV) emission. △ Less

Submitted 31 January, 2024; originally announced January 2024.

A corrected single element Maxwell visco-elastic model

Authors: Jesper Schmidt Hansen

Abstract: The original Maxwell visco-elastic constitutive model cannot predict the correct mechanical properties for most fluids. In this work, the model is generalized with respect to wave-vector and extended with a correction function. This new model has only two free parameters and avoids the attenuation-frequency locking present in the original model. Through molecular simulations it is shown that the m… ▽ More The original Maxwell visco-elastic constitutive model cannot predict the correct mechanical properties for most fluids. In this work, the model is generalized with respect to wave-vector and extended with a correction function. This new model has only two free parameters and avoids the attenuation-frequency locking present in the original model. Through molecular simulations it is shown that the model satisfactory predicts the transverse dynamics of the binary Lennard-Jones system at different temperatures, as well as water and toluene at ambient conditions. From the correction function it is shown that the viscous response is significantly reduced compared to the predictions from the original model when taking the wave-vector dependency into account. Morover, a temperature dependent characteristic length scale of maximum reduced viscous response is identified. △ Less

Submitted 24 January, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Plasmoid formation and strong radiative cooling in a driven magnetic reconnection experiment

Authors: R. Datta, K. Chandler, C. E. Myers, J. P. Chittenden, A. J. Crilly, C. Aragon, D. J. Ampleford, J. T. Banasek, A. Edens, W. R. Fox, S. B. Hansen, E. C. Harding, C. A. Jennings, H. Ji, C. C. Kuranz, S. V. Lebedev, Q. Looker, S. G. Patel, A. Porwitzky, G. A. Shipley, D. A. Uzdensky, D. A. Yager-Elorriaga, J. D. Hare

Abstract: We present results from the first experimental study of strongly radiatively-cooled magnetic reconnection. Two exploding aluminum wire arrays, driven simultaneously by the Z machine ($I_{max} = 20 \, \text{MA}$, $t_{\text{rise}} = 300 \, \text{ns}$), generate a radiatively-cooled reconnection layer ($S_L \approx 120$) in which the total cooling rate exceeds the hydrodynamic transit rate (… ▽ More We present results from the first experimental study of strongly radiatively-cooled magnetic reconnection. Two exploding aluminum wire arrays, driven simultaneously by the Z machine ($I_{max} = 20 \, \text{MA}$, $t_{\text{rise}} = 300 \, \text{ns}$), generate a radiatively-cooled reconnection layer ($S_L \approx 120$) in which the total cooling rate exceeds the hydrodynamic transit rate ($τ_{\text{hydro}}/τ_{\text{cool}} > 100$). Measurements of X-ray emission from the reconnection layer using a filtered diode ($>1$ keV) show a narrow (50 ns FWHM) burst of emission at 220 ns after current start, consistent with the formation and subsequent rapid cooling of the reconnection layer. Time-gated X-ray images of the reconnection layer show fast-moving (up to 50 km/s) hotspots inside the layer, consistent with the presence of plasmoids observed in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (at around 1.6 keV) prior to the onset of cooling, and exhibit temperatures of 170 eV, much greater than the temperature of the plasma inflows and the rest of the reconnection layer. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Simulations of Radiatively Cooled Magnetic Reconnection Driven by Pulsed Power

Authors: Rishabh Datta, Aidan J. Crilly, Jeremy P. Chittenden, Simran Chowdhry, Katherine Chandler, Nikita Chaturvedi, Clayton E. Myers, William R. Fox, Stephanie B. Hansen, Christopher A. Jennings, Hantao Ji, Carolyn C. Kuranz, Sergey V. Lebedev, Dmitri A. Uzdensky, Jack D. Hare

Abstract: Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of th… ▽ More Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the layer. In this paper, we perform 2D & 3D simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively-cooled regime. These simulations are performed in GORGON, an Eulerian resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P$_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfvénic $(M_{A} \approx 1.5)$, supersonic $(M_S \approx 4-5)$ plasma flows which generate a reconnection layer ($L/δ \approx 100, S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed Ohmic and compressional heating $τ_{cool}^{-1}/τ_A^{-1} \approx 100$; this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer. △ Less

Submitted 3 January, 2024; originally announced January 2024.

A phenomenological model for dark matter phase space distribution

Authors: Zhen Li, Steen H. Hansen

Abstract: Understanding the nature of dark matter is among the top priorities of modern physics. However, due to its inertness, detecting and studying it directly in terrestrial experiments is extremely challenging. Numerical N-body simulations currently represent the best approach for studying the particle properties and phase space distribution, assuming the collisionless nature of dark matter. These simu… ▽ More Understanding the nature of dark matter is among the top priorities of modern physics. However, due to its inertness, detecting and studying it directly in terrestrial experiments is extremely challenging. Numerical N-body simulations currently represent the best approach for studying the particle properties and phase space distribution, assuming the collisionless nature of dark matter. These simulations also address the lack of a satisfactory theory for predicting the universal properties of dark matter halos, including the density profile and velocity distribution. In this work, we propose a new phenomenological model for the dark matter phase space distribution. This model aims to provide an NFW-like density profile, velocity magnitude distribution, and velocity component distributions that align closely with simulation data. Our model is relevant both for theoretical modeling of dark matter distributions, as well as for underground detector experiments that rely on the dark matter velocity distribution for experimental analysis. △ Less

Submitted 5 March, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 9 pages, 13 figures, 2 tables

Gemini: A Family of Highly Capable Multimodal Models

Authors: Gemini Team, Rohan Anil, Sebastian Borgeaud, Jean-Baptiste Alayrac, Jiahui Yu, Radu Soricut, Johan Schalkwyk, Andrew M. Dai, Anja Hauth, Katie Millican, David Silver, Melvin Johnson, Ioannis Antonoglou, Julian Schrittwieser, Amelia Glaese, Jilin Chen, Emily Pitler, Timothy Lillicrap, Angeliki Lazaridou, Orhan Firat, James Molloy, Michael Isard, Paul R. Barham, Tom Hennigan, Benjamin Lee , et al. (1325 additional authors not shown)

Abstract: This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultr… ▽ More This report introduces a new family of multimodal models, Gemini, that exhibit remarkable capabilities across image, audio, video, and text understanding. The Gemini family consists of Ultra, Pro, and Nano sizes, suitable for applications ranging from complex reasoning tasks to on-device memory-constrained use-cases. Evaluation on a broad range of benchmarks shows that our most-capable Gemini Ultra model advances the state of the art in 30 of 32 of these benchmarks - notably being the first model to achieve human-expert performance on the well-studied exam benchmark MMLU, and improving the state of the art in every one of the 20 multimodal benchmarks we examined. We believe that the new capabilities of the Gemini family in cross-modal reasoning and language understanding will enable a wide variety of use cases. We discuss our approach toward post-training and deploying Gemini models responsibly to users through services including Gemini, Gemini Advanced, Google AI Studio, and Cloud Vertex AI. △ Less

Submitted 17 June, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Better, Not Just More: Data-Centric Machine Learning for Earth Observation

Authors: Ribana Roscher, Marc Rußwurm, Caroline Gevaert, Michael Kampffmeyer, Jefersson A. dos Santos, Maria Vakalopoulou, Ronny Hänsch, Stine Hansen, Keiller Nogueira, Jonathan Prexl, Devis Tuia

Abstract: Recent developments and research in modern machine learning have led to substantial improvements in the geospatial field. Although numerous deep learning architectures and models have been proposed, the majority of them have been solely developed on benchmark datasets that lack strong real-world relevance. Furthermore, the performance of many methods has already saturated on these datasets. We arg… ▽ More Recent developments and research in modern machine learning have led to substantial improvements in the geospatial field. Although numerous deep learning architectures and models have been proposed, the majority of them have been solely developed on benchmark datasets that lack strong real-world relevance. Furthermore, the performance of many methods has already saturated on these datasets. We argue that a shift from a model-centric view to a complementary data-centric perspective is necessary for further improvements in accuracy, generalization ability, and real impact on end-user applications. Furthermore, considering the entire machine learning cycle - from problem definition to model deployment with feedback - is crucial for enhancing machine learning models that can be reliable in unforeseen situations. This work presents a definition as well as a precise categorization and overview of automated data-centric learning approaches for geospatial data. It highlights the complementary role of data-centric learning with respect to model-centric in the larger machine learning deployment cycle. We review papers across the entire geospatial field and categorize them into different groups. A set of representative experiments shows concrete implementation examples. These examples provide concrete steps to act on geospatial data with data-centric machine learning approaches. △ Less

Submitted 22 June, 2024; v1 submitted 8 December, 2023; originally announced December 2023.

Monitoring the evolution of relative product populations at early times during a photochemical reaction

Authors: Joao Pedro Figueira Nunes, Lea Maria Ibele, Shashank Pathak, Andrew R. Attar, Surjendu Bhattacharyya, Rebecca Boll, Kurtis Borne, Martin Centurion, Benjamin Erk, Ming-Fu Lin, Ruaridh J. G. Forbes, Nate Goff, Christopher S. Hansen, Matthias Hoffmann, David M. P. Holland, Rebecca A. Ingle, Duan Luo, Sri Bhavya Muvva, Alex Reid, Arnaud Rouzée, Artem Rudenko, Sajib Kumar Saha, Xiaozhe Shen, Anbu Selvam Venkatachalam, Xijie Wang , et al. (9 additional authors not shown)

Abstract: Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging re… ▽ More Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (~50%) yield of an episulfide isomer containing a strained 3-membered ring within ~1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Quenching massive galaxies across cosmic time with the semi-analytic model SHARK v2.0

Authors: Claudia D. P. Lagos, Matias Bravo, Rodrigo Tobar, Danail Obreschkow, Chris Power, Aaron S. G. Robotham, Katy L. Proctor, Samuel Hansen, Angel Chandro-Gomez, Julian Carrivick

Abstract: We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy product… ▽ More We introduce version 2.0 of the SHARK semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being: (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy production; (iii) a model tracking the development of black hole (BH) spins; (iv) more sophisticated modelling of environmental effects on satellite galaxies; and (v) automatic parameter exploration using Particle Swarm Optimisation. We focus on two timely research topics: the structural properties of galaxies and the quenching of massive galaxies. For the former, SHARK v2.0 is capable of producing a more realistic stellar size-mass relation with a plateau marking the transition from disk- to bulge-dominated galaxies, and scaling relations between specific AM and mass that agree well with observations. For the quenching of massive galaxies, SHARK v2.0 produces massive galaxies that are more quenched than the previous version, reproducing well the observed relations between star formation rate (SFR) and stellar mass, and specific SFR and BH mass at $z=0$. SHARK v2.0 produces a number density of massive-quiescent galaxies >1dex higher than the previous version, in good agreement with JWST observations at $z\le 5$; predicts a stellar mass function of passive galaxies in reasonably good agreement with observations at $0.5<z<5$; and environmental quenching to already be effective at $z=5$. △ Less

Submitted 15 April, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS (change of terminology from radio, QSO mode to jet and wind mode, and a new section 5.3 discussing results in the context of other simulations). Supplementary material with additional comparisons with observations can be found here https://clagos.com/files/Shark_v2_SupplementaryMaterial.pdf

The destiny of Dark Matter

Authors: Fabiano Tracanna, Steen H. Hansen

Abstract: The majority of baryons, which account for $15\%$ of the matter in the Universe, will end their lives as carbon and oxygen inside cold black dwarfs. Dark matter (DM) makes up the remaining $85\%$ of the matter in the universe, however, the fate of DM is unknown. Here we show that the destiny of purely gravitationally interacting DM particles follows one of two possible routes. The first possible r… ▽ More The majority of baryons, which account for $15\%$ of the matter in the Universe, will end their lives as carbon and oxygen inside cold black dwarfs. Dark matter (DM) makes up the remaining $85\%$ of the matter in the universe, however, the fate of DM is unknown. Here we show that the destiny of purely gravitationally interacting DM particles follows one of two possible routes. The first possible route, the "radiation-destiny" scenario, is that massive DM particles lose sufficient energy through gravitational radiation causing them to spiral into a supermassive black hole that ultimately disappears through Hawking radiation. The second possible route, the "drifting-alone" destiny, applies to lighter DM particles, where only the central DM halo region spirals into the central BH which is then Hawking radiated away. The rest of the DM halo is ripped apart by the accelerated expansion of the Universe. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Comments: 15 pages, 5 figures, accepted by ApJ

State of the Art Report: Verified Computation

Authors: Jim Woodcock, Mikkel Schmidt Andersen, Diego F. Aranha, Stefan Hallerstede, Simon Thrane Hansen, Nikolaj Kuhne Jakobsen, Tomas Kulik, Peter Gorm Larsen, Hugo Daniel Macedo, Carlos Ignacio Isasa Martin, Victor Alexander Mtsimbe Norrild

Abstract: This report describes the state of the art in verifiable computation. The problem being solved is the following: The Verifiable Computation Problem (Verifiable Computing Problem) Suppose we have two computing agents. The first agent is the verifier, and the second agent is the prover. The verifier wants the prover to perform a computation. The verifier sends a description of the computation to t… ▽ More This report describes the state of the art in verifiable computation. The problem being solved is the following: The Verifiable Computation Problem (Verifiable Computing Problem) Suppose we have two computing agents. The first agent is the verifier, and the second agent is the prover. The verifier wants the prover to perform a computation. The verifier sends a description of the computation to the prover. Once the prover has completed the task, the prover returns the output to the verifier. The output will contain proof. The verifier can use this proof to check if the prover computed the output correctly. The check is not required to verify the algorithm used in the computation. Instead, it is a check that the prover computed the output using the computation specified by the verifier. The effort required for the check should be much less than that required to perform the computation. This state-of-the-art report surveys 128 papers from the literature comprising more than 4,000 pages. Other papers and books were surveyed but were omitted. The papers surveyed were overwhelmingly mathematical. We have summarised the major concepts that form the foundations for verifiable computation. The report contains two main sections. The first, larger section covers the theoretical foundations for probabilistically checkable and zero-knowledge proofs. The second section contains a description of the current practice in verifiable computation. Two further reports will cover (i) military applications of verifiable computation and (ii) a collection of technical demonstrators. The first of these is intended to be read by those who want to know what applications are enabled by the current state of the art in verifiable computation. The second is for those who want to see practical tools and conduct experiments themselves. △ Less

Submitted 16 February, 2024; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 54 pages

Bound circumplanetary orbits under the influence of radiation pressure: Application to dust in directly imaged exoplanet systems

Authors: Brad M. S. Hansen, Kevin Hayakawa

Abstract: We examine the population of simply periodic orbits in the Hill problem with radiation pressure included, in order to understand the distribution of gravitationally bound dust in orbit around a planet. We study a wide range of radiation pressure strengths, which requires the inclusion of additional terms beyond those discussed in previous analyses of this problem. In particular, our solutions reve… ▽ More We examine the population of simply periodic orbits in the Hill problem with radiation pressure included, in order to understand the distribution of gravitationally bound dust in orbit around a planet. We study a wide range of radiation pressure strengths, which requires the inclusion of additional terms beyond those discussed in previous analyses of this problem. In particular, our solutions reveal two distinct populations of stable wide, retrograde, orbits, as opposed to the single family that exists in the purely gravitational problem. We use the result of these calculations to study the observational shape of dust populations bound to extrasolar planets, that might be observable in scattered or reradiated light. In particular, we find that such dusty clouds should be elongated along the star--planet axis, and that the elongation of the bound population increases with $β$, a measure of the strength of the radiation pressure. As an application of this model, we consider the properties of the Fomalhaut system. The unusual orbital properties of the object Fomalhaut~b can be explained if the observed light was scattered by dust that was released from an object in a quasi-satellite orbit about a planet located in, or near, the observed debris ring. Within the context of the model of Hayakawa \& Hansen (2023), we find that the dust cloud around such a planet is still approximately an order of magnitude fainter than the limits set by current JWST data. △ Less

Submitted 28 February, 2024; v1 submitted 15 August, 2023; originally announced August 2023.

Comments: 21 pages, 22 figures

Journal ref: Open Journal of Astrophysics, 02/29/2024

Advanced Radiation Panel design for applications in National Security and Food Safety

Authors: A. Bross, E. C. Dukes, S. Hansen, A. Pla-Dalmau, P. Rubinov

Abstract: We describe a new concept for a basic radiation detection panel based on conventional scintillator technology and commercially available solid-state photo-detectors. The panels are simple in construction, robust, very efficient and cost-effective and are easily scalable in size, from tens of cm$^2$ to tens of m$^2$. We describe two possible applications: flagging radioactive food coontamination an… ▽ More We describe a new concept for a basic radiation detection panel based on conventional scintillator technology and commercially available solid-state photo-detectors. The panels are simple in construction, robust, very efficient and cost-effective and are easily scalable in size, from tens of cm$^2$ to tens of m$^2$. We describe two possible applications: flagging radioactive food coontamination and detection of illicit radio nucleides, such as those potentially used in a terrorist attack with a dirty bomb. △ Less

Submitted 2 April, 2024; v1 submitted 27 June, 2023; originally announced July 2023.

Comments: 10 pages, 11 figures

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Authors: Alina Kononov, Thomas Hentschel, Stephanie B. Hansen, Andrew D. Baczewski

Abstract: Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for ev… ▽ More Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminium crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via "ouroboros" effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1-2 order of magnitude cost reductions for obtaining representatively averaged and converged results. △ Less

Submitted 6 July, 2023; originally announced July 2023.

Procedural content generation of puzzle games using conditional generative adversarial networks

Authors: Andreas Hald, Jens Struckmann Hansen, Jeppe Kristensen, Paolo Burelli

Abstract: In this article, we present an experimental approach to using parameterized Generative Adversarial Networks (GANs) to produce levels for the puzzle game Lily's Garden. We extract two condition vectors from the real levels in an effort to control the details of the GAN's outputs. While the GANs perform well in approximating the first condition (map shape), they struggle to approximate the second co… ▽ More In this article, we present an experimental approach to using parameterized Generative Adversarial Networks (GANs) to produce levels for the puzzle game Lily's Garden. We extract two condition vectors from the real levels in an effort to control the details of the GAN's outputs. While the GANs perform well in approximating the first condition (map shape), they struggle to approximate the second condition (piece distribution). We hypothesize that this might be improved by trying out alternative architectures for both the Generator and Discriminator of the GANs. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: Proceedings of the 15th International Conference on the Foundations of Digital Games 2020

Dynamical Evolution of White Dwarfs in Triples in the Era of Gaia

Authors: Cheyanne Shariat, Smadar Naoz, Bradley M. S. Hansen, Isabel Angelo, Erez Michaely, Alexander P. Stephan

Abstract: The $Gaia$ mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orb… ▽ More The $Gaia$ mission has detected many white dwarfs (WDs) in binary and triple configurations, and while observations suggest that triple stellar systems are common in our Galaxy, not much attention was devoted to WDs in triples. For stability reasons, these triples must have hierarchical configurations, i.e., two stars are on a tight orbit (the inner binary), with the third companion on a wider orbit about the inner binary. In such a system, the two orbits torque each other via the eccentric Kozai-Lidov mechanism (EKL), which can alter the orbital configuration of the inner binary. We simulate thousands of triple stellar systems for over $10$ Gyr, tracking gravitational interactions, tides, general relativity, and stellar evolution up to their WD fate. As demonstrated here, three-body dynamics coupled with stellar evolution is a critical channel to form tight WD binaries or merge a WD binary. Amongst these triples, we explore their manifestations as cataclysmic variables, Type Ia supernovae, and gravitational-wave events. The simulated systems are then compared to a sample of WD triples selected from the $Gaia$ catalog. We find that including the effect of mass loss-induced kicks is crucial for producing a distribution of the inner binary-tertiary separations that is consistent with $Gaia$ observations. Lastly, we leverage this consistency to estimate that, at minimum, $30\%$ of solar-type stars in the local $200$ parsecs were born in triples. △ Less

Submitted 7 September, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 18 pages, 7 figures, revised version accepted for publication in ApJ

Impacts of permeability heterogeneity and background flow on supercritical CO2 dissolution in the deep subsurface

Authors: Scott K. Hansen, Yichen Tao, Satish Karra

Abstract: Motivated by CO2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi… ▽ More Motivated by CO2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO2 dissolution in a kilometer-scale aquifer over 1000 y. Predictive dimensionless scaling relationships relating CO2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Peclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, $\rm Pe\ = Ra^{3/4}$, was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow. △ Less

Submitted 21 May, 2023; originally announced May 2023.

MSC Class: 76R99

Debris Rings from Extrasolar Irregular Satellites

Authors: Kevin T. Hayakawa, Bradley M. S. Hansen

Abstract: Irregular satellites are the minor bodies found orbiting all four Solar System giant planets, with large semi-major axes, eccentricities, and inclinations. Previous studies have determined that the Solar System's irregular satellites are extremely collisionally evolved populations today, having lost $\sim$99 per cent of their initial mass over the course of hundreds of Myr. Such an evolution impli… ▽ More Irregular satellites are the minor bodies found orbiting all four Solar System giant planets, with large semi-major axes, eccentricities, and inclinations. Previous studies have determined that the Solar System's irregular satellites are extremely collisionally evolved populations today, having lost $\sim$99 per cent of their initial mass over the course of hundreds of Myr. Such an evolution implies that the irregular satellites must have produced a population of dusty collisional debris in the past, which is potentially observable due to the resulting reprocessing of stellar light. In this paper we examine the signatures of the debris discs produced by extrasolar analogues of this process. Radiation pressure, quantified by the parameter $β$, is the driving force behind the liberation of dust grains from the planetary Hill sphere, and results in the formation of circumstellar dust rings, even in the absence of an underlying belt of asteroids in the system. Our simulated discs reproduce many of the same features seen in some classes of observed debris discs, such as thin ring morphology, a large blowout size, and azimuthal symmetry. We compare our simulated discs' radial profiles to those of the narrow dust rings observed around Fomalhaut and HR 4796A, and show that they can broadly reproduce the observed radial distribution of dust. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 19 pages, 17 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 522, 2 (2023), 2115-2132

On inequalities involving counts of the prime factors of an odd perfect number

Authors: Graeme Clayton, Cody S. Hansen

Abstract: Let $N$ be an odd perfect number. Let $ω(N)$ be the number of distinct prime factors of $N$ and let $Ω(N)$ be the total number (counting multiplicity) of prime factors of $N$. We prove that $\frac{99}{37}ω(N) - \frac{187}{37} \leq Ω(N)$ and that if $3\nmid N$, then $\frac{51}{19}ω(N)-\frac{46}{19} \leq Ω(N)$. Let $N$ be an odd perfect number. Let $ω(N)$ be the number of distinct prime factors of $N$ and let $Ω(N)$ be the total number (counting multiplicity) of prime factors of $N$. We prove that $\frac{99}{37}ω(N) - \frac{187}{37} \leq Ω(N)$ and that if $3\nmid N$, then $\frac{51}{19}ω(N)-\frac{46}{19} \leq Ω(N)$. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 17 pages

MSC Class: 11A25 (Primary) 11A51; 11N32 (Secondary)

Causal Dependence Plots

Authors: Joshua R. Loftus, Lucius E. J. Bynum, Sakina Hansen

Abstract: Explaining artificial intelligence or machine learning models is increasingly important. To use such data-driven systems wisely we must understand how they interact with the world, including how they depend causally on data inputs. In this work we develop Causal Dependence Plots (CDPs) to visualize how one variable--an outcome--depends on changes in another variable--a predictor--… ▽ More Explaining artificial intelligence or machine learning models is increasingly important. To use such data-driven systems wisely we must understand how they interact with the world, including how they depend causally on data inputs. In this work we develop Causal Dependence Plots (CDPs) to visualize how one variable--an outcome--depends on changes in another variable--a predictor--$\textit{along with any consequent causal changes in other predictor variables}$. Crucially, CDPs differ from standard methods based on holding other predictors constant or assuming they are independent. CDPs make use of an auxiliary causal model because causal conclusions require causal assumptions. With simulations and real data experiments, we show CDPs can be combined in a modular way with methods for causal learning or sensitivity analysis. Since people often think causally about input-output dependence, CDPs can be powerful tools in the xAI or interpretable machine learning toolkit and contribute to applications like scientific machine learning and algorithmic fairness. △ Less

Submitted 5 July, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Self-Supervised Few-Shot Learning for Ischemic Stroke Lesion Segmentation

Authors: Luca Tomasetti, Stine Hansen, Mahdieh Khanmohammadi, Kjersti Engan, Liv Jorunn Høllesli, Kathinka Dæhli Kurz, Michael Kampffmeyer

Abstract: Precise ischemic lesion segmentation plays an essential role in improving diagnosis and treatment planning for ischemic stroke, one of the prevalent diseases with the highest mortality rate. While numerous deep neural network approaches have recently been proposed to tackle this problem, these methods require large amounts of annotated regions during training, which can be impractical in the medic… ▽ More Precise ischemic lesion segmentation plays an essential role in improving diagnosis and treatment planning for ischemic stroke, one of the prevalent diseases with the highest mortality rate. While numerous deep neural network approaches have recently been proposed to tackle this problem, these methods require large amounts of annotated regions during training, which can be impractical in the medical domain where annotated data is scarce. As a remedy, we present a prototypical few-shot segmentation approach for ischemic lesion segmentation using only one annotated sample during training. The proposed approach leverages a novel self-supervised training mechanism that is tailored to the task of ischemic stroke lesion segmentation by exploiting color-coded parametric maps generated from Computed Tomography Perfusion scans. We illustrate the benefits of our proposed training mechanism, leading to considerable improvements in performance in the few-shot setting. Given a single annotated patient, an average Dice score of 0.58 is achieved for the segmentation of ischemic lesions. △ Less

Submitted 16 March, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Improving dynamic collision frequencies: impacts on dynamic structure factors and stopping powers in warm dense matter

Authors: Thomas W. Hentschel, Alina Kononov, Alexandra Olmstead, Attila Cangi, Andrew D. Baczewski, Stephanie B. Hansen

Abstract: Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, we systematically investigate variations in dynamic electron-ion collision frequencies $ν(ω)$ in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum den… ▽ More Simulations and diagnostics of high-energy-density plasmas and warm dense matter rely on models of material response properties, both static and dynamic (frequency-dependent). Here, we systematically investigate variations in dynamic electron-ion collision frequencies $ν(ω)$ in warm dense matter using data from a self-consistent-field average-atom model. We show that including the full quantum density of states, strong collisions, and inelastic collisions lead to significant changes in $ν(ω)$. These changes result in red shifts and broadening of the plasmon peak in the dynamic structure factor, an effect observable in x-ray Thomson scattering spectra, and modify stopping powers around the Bragg peak. These changes improve the agreement of computationally efficient average-atom models with first-principles time-dependent density functional theory in warm dense aluminum, carbon, and deuterium. △ Less

Submitted 23 January, 2023; originally announced January 2023.

In-context Reinforcement Learning with Algorithm Distillation

Authors: Michael Laskin, Luyu Wang, Junhyuk Oh, Emilio Parisotto, Stephen Spencer, Richie Steigerwald, DJ Strouse, Steven Hansen, Angelos Filos, Ethan Brooks, Maxime Gazeau, Himanshu Sahni, Satinder Singh, Volodymyr Mnih

Abstract: We propose Algorithm Distillation (AD), a method for distilling reinforcement learning (RL) algorithms into neural networks by modeling their training histories with a causal sequence model. Algorithm Distillation treats learning to reinforcement learn as an across-episode sequential prediction problem. A dataset of learning histories is generated by a source RL algorithm, and then a causal transf… ▽ More We propose Algorithm Distillation (AD), a method for distilling reinforcement learning (RL) algorithms into neural networks by modeling their training histories with a causal sequence model. Algorithm Distillation treats learning to reinforcement learn as an across-episode sequential prediction problem. A dataset of learning histories is generated by a source RL algorithm, and then a causal transformer is trained by autoregressively predicting actions given their preceding learning histories as context. Unlike sequential policy prediction architectures that distill post-learning or expert sequences, AD is able to improve its policy entirely in-context without updating its network parameters. We demonstrate that AD can reinforcement learn in-context in a variety of environments with sparse rewards, combinatorial task structure, and pixel-based observations, and find that AD learns a more data-efficient RL algorithm than the one that generated the source data. △ Less

Submitted 25 October, 2022; originally announced October 2022.

Graphical model inference with external network data

Authors: Jack Jewson, Li Li, Laura Battaglia, Stephen Hansen, David Rossell, Piotr Zwiernik

Abstract: We consider two applications where we study how dependence structure between many variables is linked to external network data. We first study the interplay between social media connectedness and the co-evolution of the COVID-19 pandemic across USA counties. We next study study how the dependence between stock market returns across firms relates to similarities in economic and policy indicators fr… ▽ More We consider two applications where we study how dependence structure between many variables is linked to external network data. We first study the interplay between social media connectedness and the co-evolution of the COVID-19 pandemic across USA counties. We next study study how the dependence between stock market returns across firms relates to similarities in economic and policy indicators from text regulatory filings. Both applications are modelled via Gaussian graphical models where one has external network data. We develop spike-and-slab and graphical LASSO frameworks to integrate the network data, both facilitating the interpretation of the graphical model and improving inference. The goal is to detect when the network data relates to the graphical model and, if so, explain how. We found that counties strongly connected on Facebook are more likely to have similar COVID-19 evolution (positive partial correlations), accounting for various factors driving the mean. We also found that the association in stock market returns depends in a stronger fashion on economic than on policy indicators. The examples show that data integration can improve interpretation, statistical accuracy, and out-of-sample prediction, in some instances using significantly sparser graphical models. △ Less

Submitted 13 November, 2023; v1 submitted 20 October, 2022; originally announced October 2022.

ProtoVAE: A Trustworthy Self-Explainable Prototypical Variational Model

Authors: Srishti Gautam, Ahcene Boubekki, Stine Hansen, Suaiba Amina Salahuddin, Robert Jenssen, Marina MC Höhne, Michael Kampffmeyer

Abstract: The need for interpretable models has fostered the development of self-explainable classifiers. Prior approaches are either based on multi-stage optimization schemes, impacting the predictive performance of the model, or produce explanations that are not transparent, trustworthy or do not capture the diversity of the data. To address these shortcomings, we propose ProtoVAE, a variational autoencod… ▽ More The need for interpretable models has fostered the development of self-explainable classifiers. Prior approaches are either based on multi-stage optimization schemes, impacting the predictive performance of the model, or produce explanations that are not transparent, trustworthy or do not capture the diversity of the data. To address these shortcomings, we propose ProtoVAE, a variational autoencoder-based framework that learns class-specific prototypes in an end-to-end manner and enforces trustworthiness and diversity by regularizing the representation space and introducing an orthonormality constraint. Finally, the model is designed to be transparent by directly incorporating the prototypes into the decision process. Extensive comparisons with previous self-explainable approaches demonstrate the superiority of ProtoVAE, highlighting its ability to generate trustworthy and diverse explanations, while not degrading predictive performance. △ Less

Submitted 14 October, 2022; originally announced October 2022.

Consequences of dynamically unstable moons in extrasolar systems

Authors: Bradley M. S. Hansen

Abstract: Moons orbiting rocky exoplanets in compact orbits about other stars experience an accelerated tidal evolution, and can either merge with their parent planet or reach the limit of dynamical instability within a Hubble time. We review the parameter space over which moons become unbound, including the effects of atmospheric tides on the planetary spin. We find that such tides can change the final o… ▽ More Moons orbiting rocky exoplanets in compact orbits about other stars experience an accelerated tidal evolution, and can either merge with their parent planet or reach the limit of dynamical instability within a Hubble time. We review the parameter space over which moons become unbound, including the effects of atmospheric tides on the planetary spin. We find that such tides can change the final outcome from merger to escape, albeit over a limited parameter space. We also follow the further evolution of unbound moons, and demonstrate that the overwhelmingly most likely long-term outcome is that the unbound moon returns to collide with its original parent planet. The dust released by such a collision is estimated to reach optical depths approximately 0.001, exhibit characteristic temperatures of a few hundred degrees Kelvin, and last for a few thousand years. These properties make such events an attractive model for the emerging class of middle-aged main sequence stars that are observed to show transient clouds of warm dust. Furthermore, a late collision between a planet and a returning moon on a hyperbolic orbit may sterilise an otherwise habitable planet. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: to appear in Monthly Notices of the Royal Astronomical Society

A Path Towards Clinical Adaptation of Accelerated MRI

Authors: Michael S. Yao, Michael S. Hansen

Abstract: Accelerated MRI reconstructs images of clinical anatomies from sparsely sampled signal data to reduce patient scan times. While recent works have leveraged deep learning to accomplish this task, such approaches have often only been explored in simulated environments where there is no signal corruption or resource limitations. In this work, we explore augmentations to neural network MRI image recon… ▽ More Accelerated MRI reconstructs images of clinical anatomies from sparsely sampled signal data to reduce patient scan times. While recent works have leveraged deep learning to accomplish this task, such approaches have often only been explored in simulated environments where there is no signal corruption or resource limitations. In this work, we explore augmentations to neural network MRI image reconstructors to enhance their clinical relevancy. Namely, we propose a ConvNet model for detecting sources of image artifacts that achieves a classifier $F_2$ score of 79.1%. We also demonstrate that training reconstructors on MR signal data with variable acceleration factors can improve their average performance during a clinical patient scan by up to 2%. We offer a loss function to overcome catastrophic forgetting when models learn to reconstruct MR images of multiple anatomies and orientations. Finally, we propose a method for using simulated phantom data to pre-train reconstructors in situations with limited clinically acquired datasets and compute capabilities. Our results provide a potential path forward for clinical adaptation of accelerated MRI. △ Less

Submitted 28 November, 2022; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: Accepted to ML4H 2022

Journal ref: In Proceedings of the 2nd Machine Learning for Health Symposium 193:489-511, 2022

Uniqueness and Complexity of Inverse MDP Models

Authors: Marcus Hutter, Steven Hansen

Abstract: What is the action sequence aa'a" that was likely responsible for reaching state s"' (from state s) in 3 steps? Addressing such questions is important in causal reasoning and in reinforcement learning. Inverse "MDP" models p(aa'a"|ss"') can be used to answer them. In the traditional "forward" view, transition "matrix" p(s'|sa) and policy π(a|s) uniquely determine "everything": the whole dynamics p… ▽ More What is the action sequence aa'a" that was likely responsible for reaching state s"' (from state s) in 3 steps? Addressing such questions is important in causal reasoning and in reinforcement learning. Inverse "MDP" models p(aa'a"|ss"') can be used to answer them. In the traditional "forward" view, transition "matrix" p(s'|sa) and policy π(a|s) uniquely determine "everything": the whole dynamics p(as'a's"a"...|s), and with it, the action-conditional state process p(s's"...|saa'a"), the multi-step inverse models p(aa'a"...|ss^i), etc. If the latter is our primary concern, a natural question, analogous to the forward case is to which extent 1-step inverse model p(a|ss') plus policy π(a|s) determine the multi-step inverse models or even the whole dynamics. In other words, can forward models be inferred from inverse models or even be side-stepped. This work addresses this question and variations thereof, and also whether there are efficient decision/inference algorithms for this. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: 34 pages, 4 fgures

The level-1 trigger for the SuperCDMS experiment at SNOLAB

Authors: Jonathan S. Wilson, Hanno Meyer zu Theenhausen, Belina von Krosigk, Elham Azadbakht, Ray Bunker, Jeter Hall, Sten Hansen, Bruce Hines, Ben Loer, Jamieson T. Olsen, Scott M. Oser, Richard Partridge, Matthew Pyle, Joel Sander, Bruno Serfass, David Toback, Samuel L. Watkins, Xuji Zhao

Abstract: The SuperCDMS SNOLAB dark matter search experiment aims to be sensitive to energy depositions down to O(1 eV). This imposes requirements on the resolution, signal efficiency, and noise rejection of the trigger system. To accomplish this, the SuperCDMS level-1 trigger system is implemented in an FPGA on a custom PCB. A time-domain optimal filter algorithm realized as a finite impulse response filte… ▽ More The SuperCDMS SNOLAB dark matter search experiment aims to be sensitive to energy depositions down to O(1 eV). This imposes requirements on the resolution, signal efficiency, and noise rejection of the trigger system. To accomplish this, the SuperCDMS level-1 trigger system is implemented in an FPGA on a custom PCB. A time-domain optimal filter algorithm realized as a finite impulse response filter provides a baseline resolution of 0.38 times the standard deviation of the noise, $σ_n$, and a 99.9% trigger efficiency for signal amplitudes of 1.1 $σ_n$ in typical noise conditions. Embedded in a modular architecture, flexible trigger logic enables reliable triggering and vetoing in a dead-time-free manner for a variety of purposes and run conditions. The trigger architecture and performance are detailed in this article. △ Less

Submitted 23 May, 2022; v1 submitted 27 April, 2022; originally announced April 2022.

Comments: Accepted for publication in JINST

Journal ref: 2022 JINST 17 P07010

Enhancement or damping of fast neutrino flavor conversions due to collisions

Authors: Rasmus S. L. Hansen, Shashank Shalgar, Irene Tamborra

Abstract: Fast neutrino flavor conversion can occur in core-collapse supernovae or compact binary merger remnants when non-forward collisions are also at play, and neutrinos are not fully decoupled from matter. This work aims to shed light on the conditions under which fast flavor conversion is enhanced or suppressed by collisions. By relying on a neutrino toy model with three angular bins in the absence of… ▽ More Fast neutrino flavor conversion can occur in core-collapse supernovae or compact binary merger remnants when non-forward collisions are also at play, and neutrinos are not fully decoupled from matter. This work aims to shed light on the conditions under which fast flavor conversion is enhanced or suppressed by collisions. By relying on a neutrino toy model with three angular bins in the absence of spatial inhomogeneities, we consider two angular configurations: The first one with angular distributions of $ν_e$ and $\barν_e$ that are almost isotropic as expected before complete neutrino decoupling and showing little flavor conversion when collisions are absent. The second one with angular distributions of $ν_e$ and $\barν_e$ that are forward peaked as expected in the free-streaming regime and showing significant flavor conversion in the absence of collisions. By including angle-independent, direction-changing collisions, we find that collisions are responsible for an overall enhancement (damping) of flavor conversion in the former (latter) angular configuration. These opposite outcomes are due to the non-trivial interplay between collisions, flavor conversion, and the initial angular distributions of the electron type neutrinos. The enhancement in neutrino flavor conversion is found to be anticorrelated with the magnitude of flavor conversions in the absence of collisions. △ Less

Submitted 9 June, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: 12 pages, 7 figures, matches published version

Journal ref: Phys. Rev. D 105, 123003, 2022

Prime factors of $Φ_3(x)$ of the same form

Authors: Cody S. Hansen, Pace P. Nielsen

Abstract: We parameterize solutions to the equality $Φ_3(x)=Φ_3(a_1)Φ_3(a_2)\cdotsΦ_3(a_n)$ when each $Φ_3(a_i)$ is prime. Our focus is on the special cases when $n=2,3,4$, as this analysis simplifies and extends bounds on the total number of prime factors of an odd perfect number. We parameterize solutions to the equality $Φ_3(x)=Φ_3(a_1)Φ_3(a_2)\cdotsΦ_3(a_n)$ when each $Φ_3(a_i)$ is prime. Our focus is on the special cases when $n=2,3,4$, as this analysis simplifies and extends bounds on the total number of prime factors of an odd perfect number. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Journal ref: Integers 22:A71, 9 pages, 2022, available at: http://math.colgate.edu/~integers/w71/w71.pdf

When Do Stars Go BOOM?

Authors: Harvey B. Richer, Roger E. Cohen, Matteo Correnti, Ilaria Caiazzo, Jeffrey Cummings, Paul Goudfrooij, Bradley M. S. Hansen, Jeremy Heyl, Molly Peeples, Jason Kalirai, Elena Sabbi, Pier-Emmanuel Tremblay, Benjamin Williams

Abstract: The maximum mass of a star that can produce a white dwarf (WD) is an important astrophysical quantity. One of the best approaches to establishing this limit is to search for WDs in young star clusters in which only massive stars have had time to evolve and where the mass of the progenitor can be established from the cooling time of the WD together with the age of the cluster. Searches in young Mil… ▽ More The maximum mass of a star that can produce a white dwarf (WD) is an important astrophysical quantity. One of the best approaches to establishing this limit is to search for WDs in young star clusters in which only massive stars have had time to evolve and where the mass of the progenitor can be established from the cooling time of the WD together with the age of the cluster. Searches in young Milky Way clusters have not thus far yielded WD members more massive than about 1.1$~M_{\odot}$, well below the Chandrasekhar mass of $1.38~M_{\odot}$, nor progenitors with masses in excess of about $6~M_{\odot}$. However, the hunt for potentially massive WDs that escaped their cluster environs is yielding interesting candidates. To expand the cluster sample further, we used HST to survey four young and massive star clusters in the Magellanic Clouds for bright WDs that could have evolved from stars as massive as 10$~M_{\odot}$. We located five potential WD candidates in the oldest of the four clusters examined, the first extragalactic single WDs thus far discovered. As these hot WDs are very faint at optical wavelengths, final confirmation will likely have to await spectroscopy with 30-metre class telescopes. △ Less

Submitted 10 April, 2022; v1 submitted 21 March, 2022; originally announced March 2022.

Comments: 10 pages, 5 figures, accepted to the Astrophysical Journal Letters

Anomaly Detection-Inspired Few-Shot Medical Image Segmentation Through Self-Supervision With Supervoxels

Authors: Stine Hansen, Srishti Gautam, Robert Jenssen, Michael Kampffmeyer

Abstract: Recent work has shown that label-efficient few-shot learning through self-supervision can achieve promising medical image segmentation results. However, few-shot segmentation models typically rely on prototype representations of the semantic classes, resulting in a loss of local information that can degrade performance. This is particularly problematic for the typically large and highly heterogene… ▽ More Recent work has shown that label-efficient few-shot learning through self-supervision can achieve promising medical image segmentation results. However, few-shot segmentation models typically rely on prototype representations of the semantic classes, resulting in a loss of local information that can degrade performance. This is particularly problematic for the typically large and highly heterogeneous background class in medical image segmentation problems. Previous works have attempted to address this issue by learning additional prototypes for each class, but since the prototypes are based on a limited number of slices, we argue that this ad-hoc solution is insufficient to capture the background properties. Motivated by this, and the observation that the foreground class (e.g., one organ) is relatively homogeneous, we propose a novel anomaly detection-inspired approach to few-shot medical image segmentation in which we refrain from modeling the background explicitly. Instead, we rely solely on a single foreground prototype to compute anomaly scores for all query pixels. The segmentation is then performed by thresholding these anomaly scores using a learned threshold. Assisted by a novel self-supervision task that exploits the 3D structure of medical images through supervoxels, our proposed anomaly detection-inspired few-shot medical image segmentation model outperforms previous state-of-the-art approaches on two representative MRI datasets for the tasks of abdominal organ segmentation and cardiac segmentation. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: Accepted in Medical Image Analysis

Demonstrating The Risk of Imbalanced Datasets in Chest X-ray Image-based Diagnostics by Prototypical Relevance Propagation

Authors: Srishti Gautam, Marina M. -C. Höhne, Stine Hansen, Robert Jenssen, Michael Kampffmeyer

Abstract: The recent trend of integrating multi-source Chest X-Ray datasets to improve automated diagnostics raises concerns that models learn to exploit source-specific correlations to improve performance by recognizing the source domain of an image rather than the medical pathology. We hypothesize that this effect is enforced by and leverages label-imbalance across the source domains, i.e, prevalence of a… ▽ More The recent trend of integrating multi-source Chest X-Ray datasets to improve automated diagnostics raises concerns that models learn to exploit source-specific correlations to improve performance by recognizing the source domain of an image rather than the medical pathology. We hypothesize that this effect is enforced by and leverages label-imbalance across the source domains, i.e, prevalence of a disease corresponding to a source. Therefore, in this work, we perform a thorough study of the effect of label-imbalance in multi-source training for the task of pneumonia detection on the widely used ChestX-ray14 and CheXpert datasets. The results highlight and stress the importance of using more faithful and transparent self-explaining models for automated diagnosis, thus enabling the inherent detection of spurious learning. They further illustrate that this undesirable effect of learning spurious correlations can be reduced considerably when ensuring label-balanced source domain datasets. △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: To appear in ISBI 2022