Showing 1–50 of 112 results for author: Stevens, L

Binding Affinity Prediction: From Conventional to Machine Learning-Based Approaches

Authors: Xuefeng Liu, Songhao Jiang, Xiaotian Duan, Archit Vasan, Chong Liu, Chih-chan Tien, Heng Ma, Thomas Brettin, Fangfang Xia, Ian T. Foster, Rick L. Stevens

Abstract: Protein-ligand binding is the process by which a small molecule (drug or inhibitor) attaches to a target protein. The binding affinity, which refers to the strength of this interaction, is central to many important problems in bioinformatics such as drug design. An extensive amount of work has been devoted to predicting binding affinity over the past decades due to its significance. In this paper,… ▽ More Protein-ligand binding is the process by which a small molecule (drug or inhibitor) attaches to a target protein. The binding affinity, which refers to the strength of this interaction, is central to many important problems in bioinformatics such as drug design. An extensive amount of work has been devoted to predicting binding affinity over the past decades due to its significance. In this paper, we review all significant recent works, focusing on the methods, features, and benchmark datasets. We have observed a rising trend in the use of traditional machine learning and deep learning models for predicting binding affinity, accompanied by an increasing amount of data on proteins and small drug-like molecules. While prediction results are constantly improving, we also identify several open questions and potential directions that remain unexplored in the field. This paper could serve as an excellent starting point for machine learning researchers who wish to engage in the study of binding affinity, or for anyone with general interests in machine learning, drug discovery, and bioinformatics. △ Less

Submitted 29 September, 2024; originally announced October 2024.

Assessing Reusability of Deep Learning-Based Monotherapy Drug Response Prediction Models Trained with Omics Data

Authors: Jamie C. Overbeek, Alexander Partin, Thomas S. Brettin, Nicholas Chia, Oleksandr Narykov, Priyanka Vasanthakumari, Andreas Wilke, Yitan Zhu, Austin Clyde, Sara Jones, Rohan Gnanaolivu, Yuanhang Liu, Jun Jiang, Chen Wang, Carter Knutson, Andrew McNaughton, Neeraj Kumar, Gayara Demini Fernando, Souparno Ghosh, Cesar Sanchez-Villalobos, Ruibo Zhang, Ranadip Pal, M. Ryan Weil, Rick L. Stevens

Abstract: Cancer drug response prediction (DRP) models present a promising approach towards precision oncology, tailoring treatments to individual patient profiles. While deep learning (DL) methods have shown great potential in this area, models that can be successfully translated into clinical practice and shed light on the molecular mechanisms underlying treatment response will likely emerge from collabor… ▽ More Cancer drug response prediction (DRP) models present a promising approach towards precision oncology, tailoring treatments to individual patient profiles. While deep learning (DL) methods have shown great potential in this area, models that can be successfully translated into clinical practice and shed light on the molecular mechanisms underlying treatment response will likely emerge from collaborative research efforts. This highlights the need for reusable and adaptable models that can be improved and tested by the wider scientific community. In this study, we present a scoring system for assessing the reusability of prediction DRP models, and apply it to 17 peer-reviewed DL-based DRP models. As part of the IMPROVE (Innovative Methodologies and New Data for Predictive Oncology Model Evaluation) project, which aims to develop methods for systematic evaluation and comparison DL models across scientific domains, we analyzed these 17 DRP models focusing on three key categories: software environment, code modularity, and data availability and preprocessing. While not the primary focus, we also attempted to reproduce key performance metrics to verify model behavior and adaptability. Our assessment of 17 DRP models reveals both strengths and shortcomings in model reusability. To promote rigorous practices and open-source sharing, we offer recommendations for developing and sharing prediction models. Following these recommendations can address many of the issues identified in this study, improving model reusability without adding significant burdens on researchers. This work offers the first comprehensive assessment of reusability and reproducibility across diverse DRP models, providing insights into current model sharing practices and promoting standards within the DRP and broader AI-enabled scientific research community. △ Less

Submitted 18 September, 2024; originally announced September 2024.

Comments: 12 pages, 2 figures

Mixed Source Region Signatures Inside Magnetic Switchback Patches Inferred by Heavy Ion Diagnostics

Authors: Yeimy J. Rivera, Samuel T. Badman, Michael L. Stevens, Jim M. Raines, Christopher J. Owen, Kristoff Paulson, Tatiana Niembro, Stefano A. Livi, Susan T. Lepri, Enrico Landi, Jasper S. Halekas, Tamar Ervin, Ryan M. Dewey, Jesse T. Coburn, Stuart D. Bale, B. L. Alterman

Abstract: Since Parker Solar Probe's (Parker's) first perihelion pass at the Sun, large amplitude Alfvén waves grouped in patches have been observed near the Sun throughout the mission. Several formation processes for these magnetic switchback patches have been suggested with no definitive consensus. To provide insight to their formation, we examine the heavy ion properties of several adjacent magnetic swit… ▽ More Since Parker Solar Probe's (Parker's) first perihelion pass at the Sun, large amplitude Alfvén waves grouped in patches have been observed near the Sun throughout the mission. Several formation processes for these magnetic switchback patches have been suggested with no definitive consensus. To provide insight to their formation, we examine the heavy ion properties of several adjacent magnetic switchback patches around Parker's 11th perihelion pass capitalizing on a spacecraft lineup with Solar Orbiter where each samples the same solar wind streams over a large range of longitudes. Heavy ion properties (Fe/O, C$^{6+}$/C$^{5+}$, O$^{7+}$/O$^{6+}$) related to the wind's coronal origin, measured with Solar Orbiter can be linked to switchback patch structures identified near the Sun with Parker. We find that switchback patches do not contain distinctive ion and elemental compositional signatures different than the surrounding non-switchback solar wind. Both the patches and ambient wind exhibit a range of fast and slow wind qualities, indicating coronal sources with open and closed field lines in close proximity. These observations and modeling indicate switchback patches form in coronal hole boundary wind and with a range of source region magnetic and thermal properties. Furthermore, the heavy ion signatures suggest interchange reconnection and/or shear driven processes may play a role in their creation. △ Less

Submitted 5 September, 2024; originally announced September 2024.

Comments: Accepted for publication in ApJ on September 4th, 2024

In situ observations of large amplitude Alfvén waves heating and accelerating the solar wind

Authors: Yeimy J. Rivera, Samuel T. Badman, Michael L. Stevens, Jaye L. Verniero, Julia E. Stawarz, Chen Shi, Jim M. Raines, Kristoff W. Paulson, Christopher J. Owen, Tatiana Niembro, Philippe Louarn, Stefano A. Livi, Susan T. Lepri, Justin C. Kasper, Timothy S. Horbury, Jasper S. Halekas, Ryan M. Dewey, Rossana De Marco, Stuart D. Bale

Abstract: After leaving the Sun's corona, the solar wind continues to accelerate and cools, but more slowly than expected for a freely expanding adiabatic gas. We use in situ measurements from the Parker Solar Probe and Solar Orbiter spacecrafts to investigate a stream of solar wind as it traverses the inner heliosphere. The observations show heating and acceleration of the the plasma between the outer edge… ▽ More After leaving the Sun's corona, the solar wind continues to accelerate and cools, but more slowly than expected for a freely expanding adiabatic gas. We use in situ measurements from the Parker Solar Probe and Solar Orbiter spacecrafts to investigate a stream of solar wind as it traverses the inner heliosphere. The observations show heating and acceleration of the the plasma between the outer edge of the corona and near the orbit of Venus, in connection to the presence of large amplitude Alfvén waves. Alfvén waves are perturbations in the interplanetary magnetic field that transport energy. Our calculations show the damping and mechanical work performed by the Alfvén waves is sufficient to power the heating and acceleration of the fast solar wind in the inner heliosphere. △ Less

Submitted 5 September, 2024; v1 submitted 30 August, 2024; originally announced September 2024.

Comments: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science on August 30 2024, DOI: 10.1126/science.adk6953

Journal ref: Science, 385, 962-966 (2024)

Observations of Kappa Distributions in Solar Energetic Protons and Derived Thermodynamic Properties

Authors: M. E. Cuesta, A. T. Cummings, G. Livadiotis, D. J. McComas, C. M. S. Cohen, L. Y. Khoo, T. Sharma, M. M. Shen, R. Bandyopadhyay, J. S. Rankin, J. R. Szalay, H. A. Farooki, Z. Xu, G. D. Muro, M. L. Stevens, S. D. Bale

Abstract: In this paper we model the high-energy tail of observed solar energetic proton energy distributions with a kappa distribution function. We employ a technique for deriving the thermodynamic parameters of solar energetic proton populations measured by the Parker Solar Probe (PSP) Integrated Science Investigation of the Sun (IS$\odot$IS) EPI-Hi high energy telescope (HET), over energies from 10 - 60… ▽ More In this paper we model the high-energy tail of observed solar energetic proton energy distributions with a kappa distribution function. We employ a technique for deriving the thermodynamic parameters of solar energetic proton populations measured by the Parker Solar Probe (PSP) Integrated Science Investigation of the Sun (IS$\odot$IS) EPI-Hi high energy telescope (HET), over energies from 10 - 60 MeV. With this technique we explore, for the first time, the characteristic thermodynamic properties of the solar energetic protons associated with an interplanetary coronal mass ejection (ICME) and its driven shock. We find that (1) the spectral index, or equivalently, the thermodynamic parameter kappa of solar energetic protons ($κ_{\rm EP}$) gradually increases starting from the pre-ICME region (upstream of the CME-driven shock), reaching a maximum in the CME ejecta ($κ_{\rm EP} \approx 3.5$), followed by a gradual decrease throughout the trailing portion of the CME; (2) solar energetic proton temperature and density ($T_{\rm EP}$ and $n_{\rm EP}$) appear anti-correlated, a behavior consistent to sub-isothermal polytropic processes; and (3) values of $T_{\rm EP}$ and $κ_{\rm EP}$ appear are positively correlated, indicating an increasing entropy with time. Therefore, these proton populations are characterized by a complex and evolving thermodynamic behavior, consisting of multiple sub-isothermal polytropic processes, and a large-scale trend of increasing temperature, kappa, and entropy. This study and its companion study by Livadiotis et al. (2024) open a new set of procedures for investigating the thermodynamic behavior of energetic particles and their shared thermal properties. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 16 pages, 6 figures, 1 table

Entropy-Reinforced Planning with Large Language Models for Drug Discovery

Authors: Xuefeng Liu, Chih-chan Tien, Peng Ding, Songhao Jiang, Rick L. Stevens

Abstract: The objective of drug discovery is to identify chemical compounds that possess specific pharmaceutical properties toward a binding target. Existing large language models (LLMS) can achieve high token matching scores in terms of likelihood for molecule generation. However, relying solely on LLM decoding often results in the generation of molecules that are either invalid due to a single misused tok… ▽ More The objective of drug discovery is to identify chemical compounds that possess specific pharmaceutical properties toward a binding target. Existing large language models (LLMS) can achieve high token matching scores in terms of likelihood for molecule generation. However, relying solely on LLM decoding often results in the generation of molecules that are either invalid due to a single misused token, or suboptimal due to unbalanced exploration and exploitation as a consequence of the LLMs prior experience. Here we propose ERP, Entropy-Reinforced Planning for Transformer Decoding, which employs an entropy-reinforced planning algorithm to enhance the Transformer decoding process and strike a balance between exploitation and exploration. ERP aims to achieve improvements in multiple properties compared to direct sampling from the Transformer. We evaluated ERP on the SARS-CoV-2 virus (3CLPro) and human cancer cell target protein (RTCB) benchmarks and demonstrated that, in both benchmarks, ERP consistently outperforms the current state-of-the-art algorithm by 1-5 percent, and baselines by 5-10 percent, respectively. Moreover, such improvement is robust across Transformer models trained with different objectives. Finally, to further illustrate the capabilities of ERP, we tested our algorithm on three code generation benchmarks and outperformed the current state-of-the-art approach as well. Our code is publicly available at: https://github.com/xuefeng-cs/ERP. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Published in ICML2024

Calibrating the WSA model in EUHFORIA based on PSP observations

Authors: Evangelia Samara, Charles N. Arge, Rui F. Pinto, Jasmina Magdalenic, Nicolas Wijsen, Michael L. Stevens, Luciano Rodriguez, Stefaan Poedts

Abstract: We employ Parker Solar Probe (PSP) observations during the latest solar minimum period (years 2018 -2021) to calibrate the version of the Wang-Sheeley-Arge (WSA) coronal model used in the European space weather forecasting tool EUHFORIA. WSA provides a set of boundary conditions at 0.1 au necessary to initiate the heliospheric part of EUHFORIA, namely, the domain extending beyond the solar Alfveni… ▽ More We employ Parker Solar Probe (PSP) observations during the latest solar minimum period (years 2018 -2021) to calibrate the version of the Wang-Sheeley-Arge (WSA) coronal model used in the European space weather forecasting tool EUHFORIA. WSA provides a set of boundary conditions at 0.1 au necessary to initiate the heliospheric part of EUHFORIA, namely, the domain extending beyond the solar Alfvenic point. To calibrate WSA, we observationally constrain four constants in the WSA semi-empirical formula based on PSP observations. We show how the updated (after the calibration) WSA boundary conditions at 0.1 au are compared to PSP observations at similar distances, and we further propagate these conditions in the heliosphere according to EUHFORIAs magnetohydrodynamic (MHD) approach. We assess the predictions at Earth based on the Dynamic Time Warping technique. Our findings suggest that, for the period of interest, the WSA configurations which resembled optimally the PSP observations close to the Sun, were different than the ones needed to provide better predictions at Earth. One reason for this discrepancy can be attributed to the scarcity of fast solar wind velocities recorded by PSP. The calibration of the model was performed based on unexpectedly slow velocities that did not allow us to achieve generally and globally improved solar wind predictions, compared to older studies. Other reasons can be attributed to missing physical processes from the heliospheric part of EUHFORIA but also the fact that the currently employed WSA relationship, as coupled to the heliospheric MHD domain, may need a global reformulation beyond that of just updating the four constant factors that were taken into account in this study. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Origin and Properties of the Near Subsonic Solar Wind Observed by Parker Solar Probe

Authors: Wenshuai Cheng, Ying D. Liu, Hao Ran, Yiming Jiao, Michael L. Stevens, Justin C. Kasper

Abstract: We identify and examine the solar wind intervals near the sonic critical point (i.e., $M_S \sim 1$) observed by the Parker Solar Probe (PSP). The near subsonic wind intervals show similar properties: a low density, an extremely low velocity, a low proton temperature, and essentially no magnetic field deflections compared with the surrounding solar wind. The extremely low velocity is the primary co… ▽ More We identify and examine the solar wind intervals near the sonic critical point (i.e., $M_S \sim 1$) observed by the Parker Solar Probe (PSP). The near subsonic wind intervals show similar properties: a low density, an extremely low velocity, a low proton temperature, and essentially no magnetic field deflections compared with the surrounding solar wind. The extremely low velocity is the primary contributor to the near crossing of the sonic critical point rather than the sound speed, which is roughly constant in these intervals. Source tracing with a potential field source surface (PFSS) model suggests that the near subsonic intervals all connect to the boundaries inside coronal holes. Heliospheric current sheet (HCS) and partial HCS crossings around the near subsonic intervals indicate that the near subsonic wind is a transition layer between the slow and fast wind. The above scenario is consistent with the nature of the near subsonic wind as a low Mach-number boundary layer (LMBL), which facilitates the crossing of the sonic critical point at 15-20 $R_S$. Moreover, we find a dependence of the amplitude of switchbacks on the radial sonic Mach number. Magnetic field deflections essentially disappear near the sonic critical point, which suggests that switchbacks originate from above the sonic critical point. △ Less

Submitted 15 April, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ

Correlation of Coronal Mass Ejection Shock Temperature with Solar Energetic Particle Intensity

Authors: Manuel Enrique Cuesta, D. J. McComas, L. Y. Khoo, R. Bandyopadhyay, T. Sharma, M. M. Shen, J. S. Rankin, A. T. Cummings, J. R. Szalay, C. M. S. Cohen, N. A. Schwadron, R. Chhiber, F. Pecora, W. H. Matthaeus, R. A. Leske, M. L. Stevens

Abstract: Solar energetic particle (SEP) events have been observed by the Parker Solar Probe (PSP) spacecraft since its launch in 2018. These events include sources from solar flares and coronal mass ejections (CMEs). Onboard PSP is the IS\(\odot\)IS instrument suite measuring ions over energies from ~ 20 keV/nucleon to 200 MeV/nucleon and electrons from ~ 20 keV to 6 MeV. Previous studies sought to group C… ▽ More Solar energetic particle (SEP) events have been observed by the Parker Solar Probe (PSP) spacecraft since its launch in 2018. These events include sources from solar flares and coronal mass ejections (CMEs). Onboard PSP is the IS\(\odot\)IS instrument suite measuring ions over energies from ~ 20 keV/nucleon to 200 MeV/nucleon and electrons from ~ 20 keV to 6 MeV. Previous studies sought to group CME characteristics based on their plasma conditions and arrived at general descriptions with large statistical errors, leaving open questions on how to properly group CMEs based solely on their plasma conditions. To help resolve these open questions, plasma properties of CMEs have been examined in relation to SEPs. Here we reexamine one plasma property, the solar wind proton temperature, and compare it to the proton SEP intensity in a region immediately downstream of a CME-driven shock for seven CMEs observed at radial distances within 1 au. We find a statistically strong correlation between proton SEP intensity and bulk proton temperature, indicating a clear relationship between SEPs and the conditions in the solar wind. Furthermore, we propose that an indirect coupling of SEP intensity to the level of turbulence and the amount of energy dissipation that results is mainly responsible for the observed correlation between SEP intensity and proton temperature. These results are key to understanding the interaction of SEPs with the bulk solar wind in CME-driven shocks and will improve our ability to model the interplay of shock evolution and particle acceleration. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 12 pages, 4 figures, and 2 tables

Direct In Situ Measurements of a Fast Coronal Mass Ejection and Associated Structures in the Corona

Authors: Ying D. Liu, Bei Zhu, Hao Ran, Huidong Hu, Mingzhe Liu, Xiaowei Zhao, Rui Wang, Michael L. Stevens, Stuart D. Bale

Abstract: We report on the first direct in situ measurements of a fast coronal mass ejection (CME) and shock in the corona, which occurred on 2022 September 5. In situ measurements from the Parker Solar Probe (PSP) spacecraft near perihelion suggest two shocks with the second one decayed, which is consistent with more than one eruptions in coronagraph images. Despite a flank crossing, the measurements indic… ▽ More We report on the first direct in situ measurements of a fast coronal mass ejection (CME) and shock in the corona, which occurred on 2022 September 5. In situ measurements from the Parker Solar Probe (PSP) spacecraft near perihelion suggest two shocks with the second one decayed, which is consistent with more than one eruptions in coronagraph images. Despite a flank crossing, the measurements indicate unique features of the young ejecta: a plasma much hotter than the ambient medium suggestive of a hot solar source, and a large plasma $β$ implying a highly non-force-free state and the importance of thermal pressure gradient for CME acceleration and expansion. Reconstruction of the global coronal magnetic fields shows a long-duration change in the heliospheric current sheet (HCS), and the observed field polarity reversals agree with a more warped HCS configuration. Reconnection signatures are observed inside an HCS crossing as deep as the sonic critical point. As the reconnection occurs in the sub-Alfvénic wind, the reconnected flux sunward of the reconnection site can close back to the Sun, which helps balance magnetic flux in the heliosphere. The nature of the sub-Alfvénic wind after the HCS crossing as a low Mach-number boundary layer (LMBL) leads to in situ measurements of the near subsonic plasma at a surprisingly large distance. Specifically, an LMBL may provide favorable conditions for the crossings of the sonic critical point in addition to the Alfvén surface. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: Accepted for publication in the The Astrophysical Journal

On the Mesoscale Structure of CMEs at Mercury's Orbit: BepiColombo and Parker Solar Probe Observations

Authors: Erika Palmerio, Fernando Carcaboso, Leng Ying Khoo, Tarik M. Salman, Beatriz Sánchez-Cano, Benjamin J. Lynch, Yeimy J. Rivera, Sanchita Pal, Teresa Nieves-Chinchilla, Andreas J. Weiss, David Lario, Johannes Z. D. Mieth, Daniel Heyner, Michael L. Stevens, Orlando M. Romeo, Andrei N. Zhukov, Luciano Rodriguez, Christina O. Lee, Christina M. S. Cohen, Laura Rodríguez-García, Phyllis L. Whittlesey, Nina Dresing, Philipp Oleynik, Immanuel C. Jebaraj, David Fischer , et al. (5 additional authors not shown)

Abstract: On 2022 February 15, an impressive filament eruption was observed off the solar eastern limb from three remote-sensing viewpoints, namely Earth, STEREO-A, and Solar Orbiter. In addition to representing the most-distant observed filament at extreme ultraviolet wavelengths -- captured by Solar Orbiter's field of view extending to above 6 $R_{\odot}$ -- this event was also associated with the release… ▽ More On 2022 February 15, an impressive filament eruption was observed off the solar eastern limb from three remote-sensing viewpoints, namely Earth, STEREO-A, and Solar Orbiter. In addition to representing the most-distant observed filament at extreme ultraviolet wavelengths -- captured by Solar Orbiter's field of view extending to above 6 $R_{\odot}$ -- this event was also associated with the release of a fast ($\sim$2200 km$\cdot$s$^{-1}$) coronal mass ejection (CME) that was directed towards BepiColombo and Parker Solar Probe. These two probes were separated by 2$^{\circ}$ in latitude, 4$^{\circ}$ in longitude, and 0.03 au in radial distance around the time of the CME-driven shock arrival in situ. The relative proximity of the two probes to each other and to the Sun ($\sim$0.35 au) allows us to study the mesoscale structure of CMEs at Mercury's orbit for the first time. We analyse similarities and differences in the main CME-related structures measured at the two locations, namely the interplanetary shock, the sheath region, and the magnetic ejecta. We find that, despite the separation between the two spacecraft being well within the typical uncertainties associated with determination of CME geometric parameters from remote-sensing observations, the two sets of in-situ measurements display some profound differences that make understanding of the overall 3D CME structure particularly challenging. Finally, we discuss our findings within the context of space weather at Mercury's distances and in terms of the need to investigate solar transients via spacecraft constellations with small separations, which has been gaining significant attention during recent years. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: 31 pages, 13 figures, 5 tables, accepted for publication in The Astrophysical Journal

Viscoelastic mechanics of tidally induced lake drainage in the Amery grounding zone

Authors: Hanwen Zhang, Richard F. Katz, Laura A. Stevens

Abstract: Drainage of supraglacial lakes to the ice-sheet bed can occur when a hydrofracture propagates downward, driven by the weight of the water in the lake. For supraglacial lakes in the grounding zones of Antarctic glaciers, the mechanics of drainage is complicated by their proximity to the grounding line. Recently, a series of supraglacial lake-drainage events through hydrofractures was observed in th… ▽ More Drainage of supraglacial lakes to the ice-sheet bed can occur when a hydrofracture propagates downward, driven by the weight of the water in the lake. For supraglacial lakes in the grounding zones of Antarctic glaciers, the mechanics of drainage is complicated by their proximity to the grounding line. Recently, a series of supraglacial lake-drainage events through hydrofractures was observed in the Amery Ice Shelf grounding zone, East Antarctica. The lake depth at drainage varied considerably between events, raising questions about the mechanisms that induce hydrofracture, even at low lake depths. Here we use a modelling approach to investigate the contribution of tidally driven flexure to hydrofracture propagation. We model the viscoelastic response of a marine ice sheet to tides, the stresses that are induced, and the contribution of tidal stresses to hydrofracture propagation. Our results show that ocean tides and lake-water pressure together control supraglacial lake drainage through hydrofractures in the grounding zone. We give a model-based criterion that predicts supraglacial lake drainage as a function of daily maximum tidal amplitude and lake depth. Our model-based criterion agrees with remotely sensed data, indicating the importance of tidal flexure to processes associated with hydrofracturing such as supraglacial lake drainage, rifting and calving. △ Less

Submitted 6 March, 2024; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 23 pages, 10 figures

Influencing factors on false positive rates when classifying tumor cell line response to drug treatment

Authors: Priyanka Vasanthakumari, Thomas Brettin, Yitan Zhu, Hyunseung Yoo, Maulik Shukla, Alexander Partin, Fangfang Xia, Oleksandr Narykov, Rick L. Stevens

Abstract: Informed selection of drug candidates for laboratory experimentation provides an efficient means of identifying suitable anti-cancer treatments. The advancement of artificial intelligence has led to the development of computational models to predict cancer cell line response to drug treatment. It is important to analyze the false positive rate (FPR) of the models, to increase the number of effecti… ▽ More Informed selection of drug candidates for laboratory experimentation provides an efficient means of identifying suitable anti-cancer treatments. The advancement of artificial intelligence has led to the development of computational models to predict cancer cell line response to drug treatment. It is important to analyze the false positive rate (FPR) of the models, to increase the number of effective treatments identified and to minimize unnecessary laboratory experimentation. Such analysis will also aid in identifying drugs or cancer types that require more data collection to improve model predictions. This work uses an attention based neural network classification model to identify responsive/non-responsive drug treatments across multiple types of cancer cell lines. Two data filtering techniques have been applied to generate 10 data subsets, including removing samples for which dose response curves are poorly fitted and removing samples whose area under the dose response curve (AUC) values are marginal around 0.5 from the training set. One hundred trials of 10-fold cross-validation analysis is performed to test the model prediction performance on all the data subsets and the subset with the best model prediction performance is selected for further analysis. Several error analysis metrics such as the false positive rate (FPR), and the prediction uncertainty are evaluated, and the results are summarized by cancer type and drug mechanism of action (MoA) category. The FPR of cancer type spans between 0.262 and 0.5189, while that of drug MoA category spans almost the full range of [0, 1]. This study identifies cancer types and drug MoAs with high FPRs. Additional drug screening data of these cancer and drug categories may improve response modeling. Our results also demonstrate that the two data filtering approaches help improve the drug response prediction performance. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Blending Imitation and Reinforcement Learning for Robust Policy Improvement

Authors: Xuefeng Liu, Takuma Yoneda, Rick L. Stevens, Matthew R. Walter, Yuxin Chen

Abstract: While reinforcement learning (RL) has shown promising performance, its sample complexity continues to be a substantial hurdle, restricting its broader application across a variety of domains. Imitation learning (IL) utilizes oracles to improve sample efficiency, yet it is often constrained by the quality of the oracles deployed. which actively interleaves between IL and RL based on an online estim… ▽ More While reinforcement learning (RL) has shown promising performance, its sample complexity continues to be a substantial hurdle, restricting its broader application across a variety of domains. Imitation learning (IL) utilizes oracles to improve sample efficiency, yet it is often constrained by the quality of the oracles deployed. which actively interleaves between IL and RL based on an online estimate of their performance. RPI draws on the strengths of IL, using oracle queries to facilitate exploration, an aspect that is notably challenging in sparse-reward RL, particularly during the early stages of learning. As learning unfolds, RPI gradually transitions to RL, effectively treating the learned policy as an improved oracle. This algorithm is capable of learning from and improving upon a diverse set of black-box oracles. Integral to RPI are Robust Active Policy Selection (RAPS) and Robust Policy Gradient (RPG), both of which reason over whether to perform state-wise imitation from the oracles or learn from its own value function when the learner's performance surpasses that of the oracles in a specific state. Empirical evaluations and theoretical analysis validate that RPI excels in comparison to existing state-of-the-art methodologies, demonstrating superior performance across various benchmark domains. △ Less

Submitted 4 October, 2023; v1 submitted 2 October, 2023; originally announced October 2023.

Estimates of Proton and Electron Heating Rates Extended to the Near-Sun Environment

Authors: R. Bandyopadhyay, C. M. Meyer, W. H. Matthaeus, D. J. McComas, S. R. Cranmer, J. S. Halekas, J. Huang, D. E. Larson, R. Livi, A. Rahmati, P. L. Whittlesey, M. L. Stevens, J. C. Kasper, S. D. Bale

Abstract: A central problem of space plasma physics is how protons and electrons are heated in a turbulent, magnetized plasma. The differential heating of charged species due to dissipation of turbulent fluctuations plays a key role in solar wind evolution. Measurements from previous heliophysics missions have provided estimates of proton and electron heating rates beyond 0.27 au. Using Parker Solar Probe (… ▽ More A central problem of space plasma physics is how protons and electrons are heated in a turbulent, magnetized plasma. The differential heating of charged species due to dissipation of turbulent fluctuations plays a key role in solar wind evolution. Measurements from previous heliophysics missions have provided estimates of proton and electron heating rates beyond 0.27 au. Using Parker Solar Probe (PSP) data accumulated during the first ten encounters, we extend the evaluation of the individual rates of heat deposition for protons and electrons in to a distance of 0.063 au (13.5 solar radii), in the newly formed solar wind. The PSP data in the near-Sun environment show different behavior of the electron heat conduction flux from what was predicted from previous fits to Helios and Ulysses data. Consequently, the empirically derived proton and electron heating rates exhibit significantly different behavior than previous reports, with the proton heating becoming increasingly dominant over electron heating at decreasing heliocentric distances. We find that the protons receive about 80% of the total plasma heating at ~ 13 solar radii, slightly higher than the near-Earth values. This empirically derived heating partition between protons and electrons will help to constrain theoretical models of solar wind heating. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: Accepted for publication in Astrophysical Journal Letters

The Effects of Non-Equilibrium Velocity Distributions on Alfvén Ion-Cyclotron Waves in the Solar Wind

Authors: Jada Walters, Kristopher G. Klein, Emily Lichko, Michael L. Stevens, Daniel Verscharen, Benjamin D. G. Chandran

Abstract: In this work, we investigate how the complex structure found in solar wind proton velocity distribution functions (VDFs), rather than the commonly assumed two-component bi-Maxwellian structure, affects the onset and evolution of parallel-propagating microinstabilities. We use the Arbitrary Linear Plasma Solver (ALPS), a numerical dispersion solver, to find the real frequencies and growth/damping r… ▽ More In this work, we investigate how the complex structure found in solar wind proton velocity distribution functions (VDFs), rather than the commonly assumed two-component bi-Maxwellian structure, affects the onset and evolution of parallel-propagating microinstabilities. We use the Arbitrary Linear Plasma Solver (ALPS), a numerical dispersion solver, to find the real frequencies and growth/damping rates of the Alfvén modes calculated for proton VDFs extracted from Wind spacecraft observations of the solar wind. We compare this wave behavior to that obtained by applying the same procedure to core-and-beam bi-Maxwellian fits of the Wind proton VDFs. We find several significant differences in the plasma waves obtained for the extracted data and bi-Maxwellian fits, including a strong dependence of the growth/damping rate on the shape of the VDF. By application of the quasilinear diffusion operator to these VDFs, we pinpoint resonantly interacting regions in velocity space where differences in VDF structure significantly affect the wave growth and damping rates. This demonstration of the sensitive dependence of Alfvén mode behavior on VDF structure may explain why the Alfvén ion-cyclotron instability thresholds predicted by linear theory for bi-Maxwellian models of solar wind proton background VDFs do not entirely constrain spacecraft observations of solar wind proton VDFs, such as those made by the Wind spacecraft. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 11 pages, 4 figures, 1 table. To be published in The Astrophysical Journal

Transferable Graph Neural Fingerprint Models for Quick Response to Future Bio-Threats

Authors: Wei Chen, Yihui Ren, Ai Kagawa, Matthew R. Carbone, Samuel Yen-Chi Chen, Xiaohui Qu, Shinjae Yoo, Austin Clyde, Arvind Ramanathan, Rick L. Stevens, Hubertus J. J. van Dam, Deyu Lu

Abstract: Fast screening of drug molecules based on the ligand binding affinity is an important step in the drug discovery pipeline. Graph neural fingerprint is a promising method for developing molecular docking surrogates with high throughput and great fidelity. In this study, we built a COVID-19 drug docking dataset of about 300,000 drug candidates on 23 coronavirus protein targets. With this dataset, we… ▽ More Fast screening of drug molecules based on the ligand binding affinity is an important step in the drug discovery pipeline. Graph neural fingerprint is a promising method for developing molecular docking surrogates with high throughput and great fidelity. In this study, we built a COVID-19 drug docking dataset of about 300,000 drug candidates on 23 coronavirus protein targets. With this dataset, we trained graph neural fingerprint docking models for high-throughput virtual COVID-19 drug screening. The graph neural fingerprint models yield high prediction accuracy on docking scores with the mean squared error lower than $0.21$ kcal/mol for most of the docking targets, showing significant improvement over conventional circular fingerprint methods. To make the neural fingerprints transferable for unknown targets, we also propose a transferable graph neural fingerprint method trained on multiple targets. With comparable accuracy to target-specific graph neural fingerprint models, the transferable model exhibits superb training and data efficiency. We highlight that the impact of this study extends beyond COVID-19 dataset, as our approach for fast virtual ligand screening can be easily adapted and integrated into a general machine learning-accelerated pipeline to battle future bio-threats. △ Less

Submitted 14 September, 2023; v1 submitted 17 July, 2023; originally announced August 2023.

Comments: 8 pages, 5 figures, 2 tables, accepted by ICLMA2023

ACM Class: I.2.1

Observations and Modeling of Unstable Proton and Alpha Particle Velocity Distributions in Sub-Alfvenic Solar Wind at PSP Perihelia

Authors: Leon Ofman, Scott A Boardsen, Lan K Jian, Parisa Mostafavi, Jaye L Verniero, Roberto Livi, Michael McManus, Ali Rahmati, Davin Larson, Michael L Stevens

Abstract: Past observations show that solar wind (SW) acceleration occurs inside the sub-Alfvenic region, reaching the local Alfven speed at typical distances ~ 10 - 20 Rs (solar radii). Recently, Parker Solar Probe (PSP) traversed regions of sub-Alfvenic SW near perihelia in encounters E8-E12 for the first time providing data in these regions. It became evident that properties of the magnetically dominated… ▽ More Past observations show that solar wind (SW) acceleration occurs inside the sub-Alfvenic region, reaching the local Alfven speed at typical distances ~ 10 - 20 Rs (solar radii). Recently, Parker Solar Probe (PSP) traversed regions of sub-Alfvenic SW near perihelia in encounters E8-E12 for the first time providing data in these regions. It became evident that properties of the magnetically dominated SW are considerably different from the super-Alfvenic wind. For example, there are changes in relative abundances and drift of alpha particles with respect to protons, as well as in the magnitude of magnetic fluctuations. We use data of the magnetic field from the FIELDS instrument, and construct ion velocity distribution functions (VDFs) from the sub-Alfvenic regions using Solar Probe Analyzer Ions (SPAN-I) data, and run 2.5D and 3D hybrid models of proton-alpha sub-Alfvenic SW plasma. We investigate the nonlinear evolution of the ion kinetic instabilities in several case studies, and quantify the transfer of energy between the protons, alpha particles, and the kinetic waves. The models provide the 3D ion VDFs at the various stages of the instability evolution in the SW frame. By combining observational analysis with the modeling results, we gain insights on the evolution of the ion instabilities, the heating and the acceleration processes of the sub-Alfvenic SW plasma and quantify the exchange of energy between the magnetic and kinetic components. The modeling results suggest that the ion kinetic instabilities are produced locally in the SW, resulting in anisotropic heating of the ions, as observed by PSP. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: Accepted for publication in The Astrophysical Journal

The Temperature, Electron, and Pressure Characteristics of Switchbacks: Parker Solar Probe Observations

Authors: Jia Huang, Justin C. Kasper, Davin E. Larson, Michael D. McManus, Phyllis Whittlesey, Roberto Livi, Ali Rahmati, Orlando M. Romeo, Mingzhe Liu, Lan K. Jian, J. L. Verniero, Marco Velli, Samuel T. Badman, Yeimy J. Rivera, Tatiana Niembro, Kristoff Paulson, Michael L. Stevens, Anthony W. Case, Trevor A. Bowen, Marc Pulupa, Stuart D. Bale, Jasper S. Halekas

Abstract: Parker Solar Probe (PSP) observes unexpectedly prevalent switchbacks, which are rapid magnetic field reversals that last from seconds to hours, in the inner heliosphere, posing new challenges to understanding their nature, origin, and evolution. In this work, we investigate the thermal states, electron pitch angle distributions, and pressure signatures of both inside and outside switchbacks, separ… ▽ More Parker Solar Probe (PSP) observes unexpectedly prevalent switchbacks, which are rapid magnetic field reversals that last from seconds to hours, in the inner heliosphere, posing new challenges to understanding their nature, origin, and evolution. In this work, we investigate the thermal states, electron pitch angle distributions, and pressure signatures of both inside and outside switchbacks, separating a switchback into spike, transition region (TR), and quiet period (QP). Based on our analysis, we find that the proton temperature anisotropies in TRs seem to show an intermediate state between spike and QP plasmas. The proton temperatures are more enhanced in spike than in TR and QP, but the alpha temperatures and alpha-to-proton temperature ratios show the opposite trends, implying that the preferential heating mechanisms of protons and alphas are competing in different regions of switchbacks. Moreover, our results suggest that the electron integrated intensities are almost the same across the switchbacks but the electron pitch angle distributions are more isotropic inside than outside switchbacks, implying switchbacks are intact structures but strong scattering of electrons happens inside switchbacks. In addition, the examination of pressures reveals that the total pressures are comparable through an individual switchback, confirming switchbacks are pressure-balanced structures. These characteristics could further our understanding of ion heating, electron scattering, and the structure of switchbacks. △ Less

Submitted 29 August, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: published in ApJ

Journal ref: ApJ 954 133 (2023)

Quantifying the Energy Budget in the Solar Wind from 13.3-100 Solar Radii

Authors: J. S. Halekas, S. D. Bale, M. Berthomier, B. D. G. Chandran, J. F. Drake, J. C. Kasper, K. G. Klein, D. E. Larson, R. Livi, M. P. Pulupa, M. L. Stevens, J. L. Verniero, P. Whittlesey

Abstract: A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric… ▽ More A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ~100 solar radii (RS). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e. "fast radial scans"). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 RS. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Total Electron Temperature Derived from Quasi-Thermal Noise Spectroscopy In the Pristine Solar Wind: Parker Solar Probe Observations

Authors: M. Liu, K. Issautier, M. Moncuquet, N. Meyer-Vernet, M. Maksimovic, J. Huang, M. Martinovic, L. Griton, N. Chrysaphi, V. K. Jagarlamudi, S. Bale, M. Pulupa, J. C. Kasper, M. L. Stevens

Abstract: The Quasi-thermal noise (QTN) technique is a reliable tool to yield accurate measurements of the electron parameters in the solar wind. We apply this method on Parker Solar Probe (PSP) observations to derive the total electron temperature ($T_e$) from the linear fit of the high-frequency part of the QTN spectra acquired by the RFS/FIELDS instrument, and present a combination of 12-day period of ob… ▽ More The Quasi-thermal noise (QTN) technique is a reliable tool to yield accurate measurements of the electron parameters in the solar wind. We apply this method on Parker Solar Probe (PSP) observations to derive the total electron temperature ($T_e$) from the linear fit of the high-frequency part of the QTN spectra acquired by the RFS/FIELDS instrument, and present a combination of 12-day period of observations around each perihelion from Encounter One (E01) to Ten (E10) (with E08 not included) with the heliocentric distance varying from about 13 to 60 solar radii ($R_\odot{}$). We find that the total electron temperature decreases with the distance as $\sim$$R^{-0.66}$, which is much slower than adiabatic. The extrapolated $T_e$ based on PSP observations is consistent with the exospheric solar wind model prediction at $\sim$10 $R_\odot{}$, Helios observations at $\sim$0.3 AU and Wind observations at 1 AU. Also, $T_e$, extrapolated back to 10 $R_\odot{}$, is almost the same as the strahl electron temperature $T_s$ (measured by SPAN-E) which is considered to be closely related to or even almost equal to the coronal electron temperature. Furthermore, the radial $T_e$ profiles in the slower solar wind (or flux tube with larger mass flux) are steeper than those in the faster solar wind (or flux tube with smaller mass flux). More pronounced anticorrelated $V_p$-$T_e$ is observed when the solar wind is slower and closer to the Sun. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 10 pages, 7 figures, and Astronomy & Astrophysics Accepted

Journal ref: A&A 674, A49 (2023)

Prediction and Verification of Parker Solar Probe Solar Wind Sources at 13.3 R$_\odot$

Authors: Samuel T. Badman, Pete Riley, Shaela I. Jones, Tae K. Kim, Robert C. Allen, C. Nick Arge, Stuart D. Bale, Carl J. Henney, Justin C. Kasper, Parisa Mostafavi, Nikolai V. Pogorelov, Nour E. Raouafi, Michael L. Stevens, J. L. Verniero

Abstract: Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing \textit{in situ} timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (PSP), is necessary for achieving useful coordinated remote observatio… ▽ More Drawing connections between heliospheric spacecraft and solar wind sources is a vital step in understanding the evolution of the solar corona into the solar wind and contextualizing \textit{in situ} timeseries. Furthermore, making advanced predictions of this linkage for ongoing heliospheric missions, such as Parker Solar Probe (PSP), is necessary for achieving useful coordinated remote observations and maximizing scientific return. The general procedure for estimating such connectivity is straightforward (i.e. magnetic field line tracing in a coronal model) but validating the resulting estimates difficult due to the lack of an independent ground truth and limited model constraints. In its most recent orbits, PSP has reached perihelia of 13.3$R_\odot$ and moreover travels extremely fast prograde relative to the solar surface, covering over 120 degrees longitude in three days. Here we present footpoint predictions and subsequent validation efforts for PSP Encounter 10, the first of the 13.3$R_\odot$ orbits, which occurred in November 2021. We show that the longitudinal dependence of \textit{in situ} plasma data from these novel orbits provides a powerful method of footpoint validation. With reference to other encounters, we also illustrate that the conditions under which source mapping is most accurate for near-ecliptic spacecraft (such as PSP) occur when solar activity is low, but also requires that the heliospheric current sheet is strongly warped by mid-latitude or equatorial coronal holes. Lastly, we comment on the large-scale coronal structure implied by the Encounter 10 mapping, highlighting an empirical equatorial cut of the Alfvèn surface consisting of localized protrusions above unipolar magnetic separatrices. △ Less

Submitted 29 March, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: 33 Pages, 7 Figures, JGR Space Physics, Published Online 2023/3/28, In Final Production

On the evolution of the Anisotropic Scaling of Magnetohydrodynamic Turbulence in the Inner Heliosphere

Authors: Nikos Sioulas, Marco Velli, Zesen Huang, Chen Shi, Trevor A. Bowen, B. D. G. Chandran, Ioannis Liodis, Nooshin Davis, Stuart D. Bale, T. S. Horbury, Thierry Dudok de Wit, Davin Larson, Justin Kasper, Christopher J. Owen, Michael L. Stevens, Anthony Case, Marc Pulupa, David M. Malaspina, J. W. Bonnell, Keith Goetz, Peter R. Harvey, Robert J. MacDowall

Abstract: We analyze a merged Parker Solar Probe ($PSP$) and Solar Orbiter ($SO$) dataset covering heliocentric distances $13 \ R_{\odot} \lesssim R \lesssim 220$ $R_{\odot}$ to investigate the radial evolution of power and spectral-index anisotropy in the wavevector space of solar wind turbulence. Our results show that anisotropic signatures of turbulence display a distinct radial evolution when fast,… ▽ More We analyze a merged Parker Solar Probe ($PSP$) and Solar Orbiter ($SO$) dataset covering heliocentric distances $13 \ R_{\odot} \lesssim R \lesssim 220$ $R_{\odot}$ to investigate the radial evolution of power and spectral-index anisotropy in the wavevector space of solar wind turbulence. Our results show that anisotropic signatures of turbulence display a distinct radial evolution when fast, $V_{sw} \geq ~ 400 ~km ~s^{-1}$, and slow, $V_{sw} \leq ~ 400 ~km ~s^{-1}$, wind streams are considered. The anisotropic properties of slow wind in Earth orbit are consistent with a ``critically balanced'' cascade, but both spectral-index anisotropy and power anisotropy diminish with decreasing heliographic distance. Fast streams are observed to roughly retain their near-Sun anisotropic properties, with the observed spectral index and power anisotropies being more consistent with a ``dynamically aligned'' type of cascade, though the lack of extended fast-wind intervals makes it difficult to accurately measure the anisotropic scaling. A high-resolution analysis during the first perihelion of PSP confirms the presence of two sub-ranges within the inertial range, which may be associated with the transition from weak to strong turbulence. The transition occurs at $κd_{i} \approx 6 \times 10^{-2}$, and signifies a shift from -5/3 to -2 and -3/2 to -1.57 scaling in parallel and perpendicular spectra, respectively. Our results provide strong observational constraints for anisotropic theories of MHD turbulence in the solar wind. △ Less

Submitted 20 March, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: Accepted to APJ

Deep learning methods for drug response prediction in cancer: predominant and emerging trends

Authors: Alexander Partin, Thomas S. Brettin, Yitan Zhu, Oleksandr Narykov, Austin Clyde, Jamie Overbeek, Rick L. Stevens

Abstract: Cancer claims millions of lives yearly worldwide. While many therapies have been made available in recent years, by in large cancer remains unsolved. Exploiting computational predictive models to study and treat cancer holds great promise in improving drug development and personalized design of treatment plans, ultimately suppressing tumors, alleviating suffering, and prolonging lives of patients.… ▽ More Cancer claims millions of lives yearly worldwide. While many therapies have been made available in recent years, by in large cancer remains unsolved. Exploiting computational predictive models to study and treat cancer holds great promise in improving drug development and personalized design of treatment plans, ultimately suppressing tumors, alleviating suffering, and prolonging lives of patients. A wave of recent papers demonstrates promising results in predicting cancer response to drug treatments while utilizing deep learning methods. These papers investigate diverse data representations, neural network architectures, learning methodologies, and evaluations schemes. However, deciphering promising predominant and emerging trends is difficult due to the variety of explored methods and lack of standardized framework for comparing drug response prediction models. To obtain a comprehensive landscape of deep learning methods, we conducted an extensive search and analysis of deep learning models that predict the response to single drug treatments. A total of 60 deep learning-based models have been curated and summary plots were generated. Based on the analysis, observable patterns and prevalence of methods have been revealed. This review allows to better understand the current state of the field and identify major challenges and promising solution paths. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Towards a Verified Tableau Prover for a Quantifier-Free Fragment of Set Theory

Authors: Lukas Stevens

Abstract: Using Isabelle/HOL, we verify the state-of-the-art decision procedure for multi-level syllogistic with singleton (MLSS for short), which is a quantifier-free fragment of set theory. We formalise its syntax and semantics as well as a sound and complete tableau calculus for it. We also provide an executable specification of a decision procedure that exhaustively applies the rules of the calculus and… ▽ More Using Isabelle/HOL, we verify the state-of-the-art decision procedure for multi-level syllogistic with singleton (MLSS for short), which is a quantifier-free fragment of set theory. We formalise its syntax and semantics as well as a sound and complete tableau calculus for it. We also provide an executable specification of a decision procedure that exhaustively applies the rules of the calculus and prove its termination. Furthermore, we extend the calculus with a lightweight type system that paves the way for an integration of the procedure into Isabelle/HOL. △ Less

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

Magnetic field spectral evolution in the inner heliosphere

Authors: Nikos Sioulas, Zesen Huang, Chen Shi, Marco Velli, Anna Tenerani, Loukas Vlahos, Trevor A. Bowen, Stuart D. Bale, J. W. Bonnell, P. R. Harvey, Davin Larson, arc Pulupa, Roberto Livi, L. D. Woodham, T. S. Horbury, Michael L. Stevens, T. Dudok de Wit, R. J. MacDowall, David M. Malaspina, K. Goetz, Jia Huang, Justin Kasper, Christopher J. Owen, Milan Maksimović, P. Louarn , et al. (1 additional authors not shown)

Abstract: Parker Solar Probe and Solar Orbiter data are used to investigate the radial evolution of magnetic turbulence between $0.06 ~ \lesssim R ~\lesssim 1$ au. The spectrum is studied as a function of scale, normalized to the ion inertial scale $d_{i}$. In the vicinity of the Sun, the inertial range is limited to a narrow range of scales and exhibits a power-law exponent of, $α_{B} = -3/2$, independent… ▽ More Parker Solar Probe and Solar Orbiter data are used to investigate the radial evolution of magnetic turbulence between $0.06 ~ \lesssim R ~\lesssim 1$ au. The spectrum is studied as a function of scale, normalized to the ion inertial scale $d_{i}$. In the vicinity of the Sun, the inertial range is limited to a narrow range of scales and exhibits a power-law exponent of, $α_{B} = -3/2$, independent of plasma parameters. The inertial range grows with distance, progressively extending to larger spatial scales, while steepening towards a $α_{B} =-5/3$ scaling. It is observed that spectra for intervals with large magnetic energy excesses and low Alfvénic content steepen significantly with distance, in contrast to highly Alfvénic intervals that retain their near-Sun scaling. The occurrence of steeper spectra in slower wind streams may be attributed to the observed positive correlation between solar wind speed and Alfvénicity. △ Less

Submitted 28 December, 2022; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: Accepted to APJ letters with minor revisions

Engaging, Large-Scale Functional Programming Education in Physical and Virtual Space

Authors: Kevin Kappelmann, Jonas Rädle, Lukas Stevens

Abstract: Worldwide, computer science departments have experienced a dramatic increase in the number of student enrolments. Moreover, the ongoing COVID-19 pandemic requires institutions to radically replace the traditional way of on-site teaching, moving interaction from physical to virtual space. We report on our strategies and experience tackling these issues as part of a Haskell-based functional programm… ▽ More Worldwide, computer science departments have experienced a dramatic increase in the number of student enrolments. Moreover, the ongoing COVID-19 pandemic requires institutions to radically replace the traditional way of on-site teaching, moving interaction from physical to virtual space. We report on our strategies and experience tackling these issues as part of a Haskell-based functional programming and verification course, accommodating over 2000 students in the course of two semesters. Among other things, we fostered engagement with weekly programming competitions and creative homework projects, workshops with industry partners, and collaborative pair-programming tutorials. To offer such an extensive programme to hundreds of students, we automated feedback for programming as well as inductive proof exercises. We explain and share our tools and exercises so that they can be reused by other educators. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: In Proceedings TFPIE 2021/22, arXiv:2207.11600

Journal ref: EPTCS 363, 2022, pp. 93-113

A Linter for Isabelle: Implementation and Evaluation

Authors: Yecine Megdiche, Fabian Huch, Lukas Stevens

Abstract: In interactive theorem proving, formalization quality is a key factor for maintainability and re-usability of developments and can also impact proof-checking performance. Commonly, anti-patterns that cause quality issues are known to experienced users. However, in many theorem prover systems, there are no automatic tools to check for their presence and make less experienced users aware of them. We… ▽ More In interactive theorem proving, formalization quality is a key factor for maintainability and re-usability of developments and can also impact proof-checking performance. Commonly, anti-patterns that cause quality issues are known to experienced users. However, in many theorem prover systems, there are no automatic tools to check for their presence and make less experienced users aware of them. We attempt to fill this gap in the Isabelle environment by developing a linter as a publicly available add-on component. The linter offers basic configurability, extensibility, Isabelle/jEdit integration, and a standalone command-line tool. We uncovered 480 potential problems in Isabelle/HOL, 14016 in other formalizations of the Isabelle distribution, and an astonishing 59573 in the AFP. With a specific lint bundle for AFP submissions, we found that submission guidelines were violated in 1595 cases. We set out to alleviate problems in Isabelle/HOL and solved 168 of them so far; we found that high-severity lints corresponded to actual problems most of the time, individual users often made the same mistakes in many places, and that solving those problems retrospectively amounts to a substantial amount of work. In contrast, solving these problems interactively for new developments usually incurs only little overhead, as we found in a quantitative user survey with 22 participants (less than a minute for more than 60% of participants). We also found that a good explanation of problems is key to the users' ease of solving these problems (correlation coefficient 0.48), and their satisfaction with the end result (correlation coefficient 0.62). △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: Isabelle Workshop 2022

The Radial Evolution of the Solar Wind as Organized by Electron Distribution Parameters

Authors: J. S. Halekas, P. Whittlesey, D. E. Larson, M. Maksimovic, R. Livi, M. Berthomier, J. C. Kasper, A. W. Case, M. L. Stevens, S. D. Bale, R. J. MacDowall, M. P. Pulupa

Abstract: We utilize observations from the Parker Solar Probe (PSP) to study the radial evolution of the solar wind in the inner heliosphere. We analyze electron velocity distribution functions observed by the Solar Wind Electrons, Alphas, and Protons suite to estimate the coronal electron temperature and the local electric potential in the solar wind. From the latter value and the local flow speed, we comp… ▽ More We utilize observations from the Parker Solar Probe (PSP) to study the radial evolution of the solar wind in the inner heliosphere. We analyze electron velocity distribution functions observed by the Solar Wind Electrons, Alphas, and Protons suite to estimate the coronal electron temperature and the local electric potential in the solar wind. From the latter value and the local flow speed, we compute the asymptotic solar wind speed. We group the PSP observations by asymptotic speed, and characterize the radial evolution of the wind speed, electron temperature, and electric potential within each group. In agreement with previous work, we find that the electron temperature (both local and coronal) and the electric potential are anti-correlated with wind speed. This implies that the electron thermal pressure and the associated electric field can provide more net acceleration in the slow wind than in the fast wind. We then utilize the inferred coronal temperature and the extrapolated electric + gravitational potential to show that both electric field driven exospheric models and the equivalent thermally driven hydrodynamic models can explain the entire observed speed of the slowest solar wind streams. On the other hand, neither class of model can explain the observed speed of the faster solar wind streams, which thus require additional acceleration mechanisms. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: Submitted to the Astrophysical Journal

Cost-Effective Online Contextual Model Selection

Authors: Xuefeng Liu, Fangfang Xia, Rick L. Stevens, Yuxin Chen

Abstract: How can we collect the most useful labels to learn a model selection policy, when presented with arbitrary heterogeneous data streams? In this paper, we formulate this task as an online contextual active model selection problem, where at each round the learner receives an unlabeled data point along with a context. The goal is to output the best model for any given context without obtaining an exce… ▽ More How can we collect the most useful labels to learn a model selection policy, when presented with arbitrary heterogeneous data streams? In this paper, we formulate this task as an online contextual active model selection problem, where at each round the learner receives an unlabeled data point along with a context. The goal is to output the best model for any given context without obtaining an excessive amount of labels. In particular, we focus on the task of selecting pre-trained classifiers, and propose a contextual active model selection algorithm (CAMS), which relies on a novel uncertainty sampling query criterion defined on a given policy class for adaptive model selection. In comparison to prior art, our algorithm does not assume a globally optimal model. We provide rigorous theoretical analysis for the regret and query complexity under both adversarial and stochastic settings. Our experiments on several benchmark classification datasets demonstrate the algorithm's effectiveness in terms of both regret and query complexity. Notably, to achieve the same accuracy, CAMS incurs less than 10% of the label cost when compared to the best online model selection baselines on CIFAR10. △ Less

Submitted 17 February, 2023; v1 submitted 13 July, 2022; originally announced July 2022.

Magnetic field intermittency in the solar wind: PSP and SolO observations ranging from the Alfven region out to 1 AU

Authors: Nikos Sioulas, Zesen Huang, Marco Velli, Rohit Chhiber, Manuel E. Cuesta, Chen Shi, William H. Matthaeus, Riddhi Bandyopadhyay, Loukas Vlahos, Trevor A. Bowen, Ramiz A. Qudsi, Stuart D. Bale, Christopher J. Owen, P. Louarn, A. Fedorov, Milan Maksimovic, Michael L. Stevens, Justin Kasper, Davin Larson, Roberto Livi

Abstract: $PSP$ and $SolO$ data are utilized to investigate magnetic field intermittency in the solar wind (SW). Small-scale intermittency $(20-100d_{i})$ is observed to radially strengthen when methods relying on higher-order moments are considered ($SF_q$, $SDK… ▽ More $PSP$ and $SolO$ data are utilized to investigate magnetic field intermittency in the solar wind (SW). Small-scale intermittency $(20-100d_{i})$ is observed to radially strengthen when methods relying on higher-order moments are considered ($SF_q$, $SDK$), but no clear trend is observed at larger scales. However, lower-order moment-based methods (e.g., PVI) are deemed more appropriate for examining the evolution of the bulk of Coherent Structures (CSs), $PVI \ge 3$. Using PVI, we observe a scale-dependent evolution in the fraction of the dataset occupied by CSs, $f_{PVI \ge 3}$. Specifically, regardless of the SW speed, a subtle increase is found in $f_{PVI\ge3}$ for $\ell =20 d_i$, in contrast to a more pronounced radial increase in CSs observed at larger scales. Intermittency is investigated in relation to plasma parameters. Though, slower SW speed intervals exhibit higher $f_{PVI \geq 6}$ and higher kurtosis maxima, no statistical differences are observed for $f_{PVI \geq 3}$. Highly Alfvénic intervals, display lower levels of intermittency. The anisotropy with respect to the angle between the magnetic field and SW flow, $Θ_{VB}$ is investigated. Intermittency is weaker at $Θ_{VB} \approx 0^{\circ}$ and is strengthened at larger angles. Considering the evolution at a constant alignment angle, a weakening of intermittency is observed with increasing advection time of the SW. Our results indicate that the strengthening of intermittency in the inner heliosphere is driven by the increase in comparatively highly intermittent perpendicular intervals sampled by the probes with increasing distance, an effect related directly to the evolution of the Parker spiral. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Eruption and Interplanetary Evolution of a Stealthy Streamer-Blowout CME Observed by PSP at ${\sim}$0.5~AU

Authors: Sanchita Pal, Benjamin J. Lynch, Simon W. Good, Erika Palmerio, Eleanna Asvestari, Jens Pomoell, Michael L. Stevens, Emilia K. J. Kilpua

Abstract: Streamer-blowout coronal mass ejections (SBO-CMEs) are the dominant CME population during solar minimum. Although they are typically slow and lack clear low-coronal signatures, they can cause geomagnetic storms. With the aid of extrapolated coronal fields and remote observations of the off-limb low corona, we study the initiation of an SBO-CME preceded by consecutive CME eruptions consistent with… ▽ More Streamer-blowout coronal mass ejections (SBO-CMEs) are the dominant CME population during solar minimum. Although they are typically slow and lack clear low-coronal signatures, they can cause geomagnetic storms. With the aid of extrapolated coronal fields and remote observations of the off-limb low corona, we study the initiation of an SBO-CME preceded by consecutive CME eruptions consistent with a multi-stage sympathetic breakout scenario. From inner-heliospheric Parker Solar Probe (PSP) observations, it is evident that the SBO-CME is interacting with the heliospheric magnetic field and plasma sheet structures draped about the CME flux rope. We estimate that $18 \, \pm \, 11\%$ of the CME's azimuthal magnetic flux has been eroded through magnetic reconnection and that this erosion began after a heliospheric distance of ${\sim}0.35$ AU from the Sun was reached. This observational study has important implications for understanding the initiation of SBO-CMEs and their interaction with the heliospheric surroundings. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 21 pages, 6 figures, 3 videos

Basal hydrofractures near sticky patches

Authors: Hanwen Zhang, Timothy Davis, Richard F. Katz, Laura A. Stevens, Dave A. May

Abstract: Basal crevasses are macroscopic structural discontinuities at the base of ice sheets and glaciers. Motivated by observations and the mechanics of elastic fracture, we hypothesise that in the presence of basal water pressure, spatial variations in basal stress can promote and localise basal crevassing. We quantify this process in the theoretical context of linear elastic fracture mechanics. We deve… ▽ More Basal crevasses are macroscopic structural discontinuities at the base of ice sheets and glaciers. Motivated by observations and the mechanics of elastic fracture, we hypothesise that in the presence of basal water pressure, spatial variations in basal stress can promote and localise basal crevassing. We quantify this process in the theoretical context of linear elastic fracture mechanics. We develop a model evaluating the effect of shear stress variation on the growth of basal crevasses. Our results indicate that sticky patches promote the initiation of basal crevasses, increase their length of propagation into the ice and, under some conditions, give them curved trajectories that incline upstream. A detailed exploration of the parameter space is conducted to gain a better understanding of the conditions under which sticky-patch-induced basal crevassing likely occurs beneath ice sheets and glaciers. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 29 pages, 9 figures

Data augmentation and multimodal learning for predicting drug response in patient-derived xenografts from gene expressions and histology images

Authors: Alexander Partin, Thomas Brettin, Yitan Zhu, James M. Dolezal, Sara Kochanny, Alexander T. Pearson, Maulik Shukla, Yvonne A. Evrard, James H. Doroshow, Rick L. Stevens

Abstract: Patient-derived xenografts (PDXs) are an appealing platform for preclinical drug studies because the in vivo environment of PDXs helps preserve tumor heterogeneity and usually better mimics drug response of patients with cancer compared to CCLs. We investigate multimodal neural network (MM-Net) and data augmentation for drug response prediction in PDXs. The MM-Net learns to predict response using… ▽ More Patient-derived xenografts (PDXs) are an appealing platform for preclinical drug studies because the in vivo environment of PDXs helps preserve tumor heterogeneity and usually better mimics drug response of patients with cancer compared to CCLs. We investigate multimodal neural network (MM-Net) and data augmentation for drug response prediction in PDXs. The MM-Net learns to predict response using drug descriptors, gene expressions (GE), and histology whole-slide images (WSIs) where the multi-modality refers to the tumor features. We explore whether the integration of WSIs with GE improves predictions as compared with models that use GE alone. We use two methods to address the limited number of response values: 1) homogenize drug representations which allows to combine single-drug and drug-pairs treatments into a single dataset, 2) augment drug-pair samples by switching the order of drug features which doubles the sample size of all drug-pair samples. These methods enable us to combine single-drug and drug-pair treatments, allowing us to train multimodal and unimodal neural networks (NNs) without changing architectures or the dataset. Prediction performance of three unimodal NNs which use GE are compared to assess the contribution of data augmentation methods. NN that uses the full dataset which includes the original and the augmented drug-pair treatments as well as single-drug treatments significantly outperforms NNs that ignore either the augmented drug-pairs or the single-drug treatments. In assessing the contribution of multimodal learning based on the MCC metric, MM-Net statistically significantly outperforms all the baselines. Our results show that data augmentation and integration of histology images with GE can improve prediction performance of drug response in PDXs. △ Less

Submitted 25 April, 2022; originally announced April 2022.

CMEs and SEPs During November-December 2020: A Challenge for Real-Time Space Weather Forecasting

Authors: Erika Palmerio, Christina O. Lee, M. Leila Mays, Janet G. Luhmann, David Lario, Beatriz Sánchez-Cano, Ian G. Richardson, Rami Vainio, Michael L. Stevens, Christina M. S. Cohen, Konrad Steinvall, Christian Möstl, Andreas J. Weiss, Teresa Nieves-Chinchilla, Yan Li, Davin E. Larson, Daniel Heyner, Stuart D. Bale, Antoinette B. Galvin, Mats Holmström, Yuri V. Khotyaintsev, Milan Maksimovic, Igor G. Mitrofanov

Abstract: Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric sy… ▽ More Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric system, and multi-point measurements of solar transients can provide useful insights and validations for prediction models. In this work, we aim to analyse the whole inner heliospheric context between two eruptive flares that took place in late 2020, i.e. the M4.4 flare of November 29 and the C7.4 flare of December 7. This period is especially interesting because the STEREO-A spacecraft was located ~60° east of the Sun-Earth line, giving us the opportunity to test the capabilities of "predictions at 360°" using remote-sensing observations from the Lagrange L1 and L5 points as input. We simulate the CMEs that were ejected during our period of interest and the SEPs accelerated by their shocks using the WSA-Enlil-SEPMOD modelling chain and four sets of input parameters, forming a "mini-ensemble". We validate our results using in-situ observations at six locations, including Earth and Mars. We find that, despite some limitations arising from the models' architecture and assumptions, CMEs and shock-accelerated SEPs can be reasonably studied and forecast in real time at least out to several tens of degrees away from the eruption site using the prediction tools employed here. △ Less

Submitted 5 May, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 50 pages, 16 figures, 3 tables, accepted for publication in Space Weather (v2 contains updated CCMC links after their website update)

Radius Constraints from Reflection Modeling of Cygnus X-2 with NuSTAR and NICER

Authors: R. M. Ludlam, E. M. Cackett, J. A. García, J. M. Miller, A. L. Stevens, A. C. Fabian, J. Homan, M. NG, S. Guillot, D. J. K. Buisson, D. Chakrabarty

Abstract: We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in t… ▽ More We present a spectral analysis of NuSTAR and NICER observations of the luminous, persistently accreting neutron star (NS) low-mass X-ray binary Cygnus X-2. The data were divided into different branches that the source traces out on the Z-track of the X-ray color-color diagram; namely the horizontal branch, normal branch, and the vertex between the two. The X-ray continuum spectrum was modeled in two different ways that produced a comparable quality fit. The spectra showed clear evidence of a reflection component in the form of a broadened Fe K line, as well as a lower energy emission feature near 1 keV likely due to an ionized plasma located far from the innermost accretion disk. We account for the reflection spectrum with two independent models (relxillns and rdblur*rfxconv). The inferred inclination is in agreement with earlier estimates from optical observations of ellipsoidal light curve modeling (relxillns: $i=67^{\circ}\pm4^{\circ}$, rdblur*rfxconv: $i=60^{\circ}\pm10^{\circ}$). The inner disk radius remains close to the NS ($R_{\rm in}\leq1.15\ R_{\mathrm{ISCO}}$) regardless of the source position along the Z-track or how the 1 keV feature is modeled. Given the optically determined NS mass of $1.71\pm0.21\ M_{\odot}$, this corresponds to a conservative upper limit of $R_{\rm in}\leq19.5$ km for $M=1.92\ M_{\odot}$ or $R_{\rm in}\leq15.3$ km for $M=1.5\ M_{\odot}$. We compare these radius constraints to those obtained from NS gravitational wave merger events and recent NICER pulsar light curve modeling measurements. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ, 13 pages, 5 tables, 5 figures

Direct First PSP Observation of the Interaction of Two Successive Interplanetary Coronal Mass Ejections in November 2020

Authors: Teresa Nieves-Chinchilla, Nathalia Alzate, Hebe Cremades, Laura Rodriguez-Garcia, Luiz F. G. Dos Santos, Ayris Narock, Hong Xie, Adam Szabo Vratislav Krupar, Marc Pulupa, David Lario, Michael L. Stevens, Erika Palmerio, Lynn B. Wilson III, Katharine K. Reeves Ryun-Young Kwon, M. Leila Mays, O. Chris St. Cyr, Phillip Hess, Daniel B. Seaton, Tatiana Niembro, Stuart D. Bale, Justin C. Kasper

Abstract: We investigate the effects of the evolutionary processes in the internal magnetic structure of two interplanetary coronal mass ejections (ICMEs) detected in situ between 2020 November 29 and December 1 by Parker Solar Probe (PSP). The sources of the ICMEs were observed remotely at the Sun in EUV and subsequently tracked to their coronal counterparts in white light. This period is of particular int… ▽ More We investigate the effects of the evolutionary processes in the internal magnetic structure of two interplanetary coronal mass ejections (ICMEs) detected in situ between 2020 November 29 and December 1 by Parker Solar Probe (PSP). The sources of the ICMEs were observed remotely at the Sun in EUV and subsequently tracked to their coronal counterparts in white light. This period is of particular interest to the community since it has been identified as the first widespread solar energetic particle event of Solar Cycle 25. The distribution of various solar and heliospheric-dedicated spacecraft throughout the inner heliosphere during PSP observations of these large-scale magnetic structures enables a comprehensive analysis of the internal evolution and topology of such structures. By assembling different models and techniques, we identify the signatures of interaction between the two consecutive ICMEs and the implications for their internal structure. We use multispacecraft observations in combination with a remote-sensing forward modeling technique, numerical propagation models, and in-situ reconstruction techniques. The outcome, from the full reconciliations, demonstrates that the two CMEs are interacting in the vicinity of PSP. Thus, we identify the in-situ observations based on the physical processes that are associated with the interaction and collision of both CMEs. We also expand the flux rope modeling and in-situ reconstruction technique to incorporate the aging and expansion effects in a distorted internal magnetic structure and explore the implications of both effects in the magnetic configuration of the ICMEs. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Journal ref: The Astrophysical Journal, 930:88, 2022

Statistical analysis of intermittency and its association with proton heating in the near Sun environment

Authors: Nikos Sioulas, Marco Velli, Rohit Chhiber, Loukas Vlahos, William H. Matthaeus, Riddhi Bandyopadhyay, Manuel E. Cuesta, Chen Shi, Trevor A. Bowen, Ramiz A. Qudsi, Michael L. Stevens, Stuart D. Bale

Abstract: We use data from the first six encounters of Parker Solar Probe and employ the Partial Variance of Increments ($PVI$) method to study the statistical properties of coherent structures in the inner heliosphere with the aim of exploring physical connections between magnetic field intermittency and observable consequences such as plasma heating and turbulence dissipation. Our results support proton h… ▽ More We use data from the first six encounters of Parker Solar Probe and employ the Partial Variance of Increments ($PVI$) method to study the statistical properties of coherent structures in the inner heliosphere with the aim of exploring physical connections between magnetic field intermittency and observable consequences such as plasma heating and turbulence dissipation. Our results support proton heating localized in the vicinity of, and strongly correlated with, magnetic structures characterized by $PVI \geq 1$. We show that on average, such events constitute $\approx 19\%$ of the dataset, though variations may occur depending on the plasma parameters. We show that the waiting time distribution ($WT$) of identified events is consistent across all six encounters following a power-law scaling at lower $WTs$. This result indicates that coherent structures are not evenly distributed in the solar wind but rather tend to be tightly correlated and form clusters. We observe that the strongest magnetic discontinuities, $PVI \geq 6$, usually associated with reconnection exhausts, are sites where magnetic energy is locally dissipated in proton heating and are associated with the most abrupt changes in proton temperature. However, due to the scarcity of such events, their relative contribution to energy dissipation is minor. Taking clustering effects into consideration, we show that smaller scale, more frequent structures with PVI between, $1\lesssim PVI \lesssim 6$, play the major role in magnetic energy dissipation. The number density of such events is strongly associated with the global solar wind temperature, with denser intervals being associated with higher $T_{p}$. △ Less

Submitted 25 January, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted for publication in APJ

Sliding across a surface: particles with fixed and mobile ligands

Authors: Janna Lowensohn, Laurie Stevens, Daniel Goldstein, Bortolo Matteo Mognetti

Abstract: A quantitative model of the mobility of functionalized particles at the interface is pivotal to understanding important systems in biology and nanotechnology. In this work, we investigate the emerging dynamics of particles anchored through ligand-receptor bridges to functionalized surfaces. We consider systems with reversible bridges in which ligand-receptor pairs bind/unbind with finite reaction… ▽ More A quantitative model of the mobility of functionalized particles at the interface is pivotal to understanding important systems in biology and nanotechnology. In this work, we investigate the emerging dynamics of particles anchored through ligand-receptor bridges to functionalized surfaces. We consider systems with reversible bridges in which ligand-receptor pairs bind/unbind with finite reaction rates. For a given set of bridges, the particle can explore a tiny fraction of the surface as the extensivity of the bridges is finite. We show how at time scales longer than the bridges' lifetime, the averaged position of the particle diffuses away from its initial value. We distill our findings into two analytic equations for the sliding diffusion constant of particles carrying mobile and fixed ligands. We quantitatively validate our theoretical predictions using reaction-diffusion simulations. Our results, along with recent literature, will allow inferring the microscopic parameters at play in complex biological systems from experimental trajectories. △ Less

Submitted 5 April, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Version accepted for publication in The Journal of Chemical Physics. The first two authors contributed equally

Spatial epidemiology and adaptive targeted sampling to manage the Chagas disease vector Triatoma dimidiata

Authors: B. K. M. Case, Jean-Gabriel Young, Daniel Penados, Carlota Monroy, Laurent Hébert-Dufresne, Lori Stevens

Abstract: Widespread application of insecticide remains the primary form of control for Chagas disease in Central America, despite only temporarily reducing domestic levels of the endemic vector Triatoma dimidiata and having little long-term impact. Recently, an approach emphasizing community feedback and housing improvements has been shown to yield lasting results. However, the additional resources and per… ▽ More Widespread application of insecticide remains the primary form of control for Chagas disease in Central America, despite only temporarily reducing domestic levels of the endemic vector Triatoma dimidiata and having little long-term impact. Recently, an approach emphasizing community feedback and housing improvements has been shown to yield lasting results. However, the additional resources and personnel required by such an intervention likely hinders its widespread adoption. One solution to this problem would be to target only a subset of houses in a community while still eliminating enough infestations to interrupt disease transfer. Here we develop a sequential sampling framework that adapts to information specific to a community as more houses are visited, thereby allowing us to efficiently find homes with domiciliary vectors while minimizing sampling bias. The method fits Bayesian geostatistical models to make spatially informed predictions, while gradually transitioning from prioritizing houses based on prediction uncertainty to targeting houses with a high risk of infestation. A key feature of the method is the use of a single exploration parameter, $α$, to control the rate of transition between these two design targets. In a simulation study using empirical data from five villages in southeastern Guatemala, we test our method using a range of values for $α$, and find it can consistently select fewer homes than random sampling, while still bringing the village infestation rate below a given threshold. We further find that when additional socioeconomic information is available, much larger savings are possible, but that meeting the target infestation rate is less consistent, particularly among the less exploratory strategies. Our results suggest new options for implementing long-term T. dimidiata control. △ Less

Submitted 10 February, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Exploring the Solar Wind from its Source on the Corona into the Inner Heliosphere during the First Solar Orbiter - Parker Solar Probe Quadrature

Authors: Daniele Telloni, Vincenzo Andretta, Ester Antonucci, Alessandro Bemporad, Giuseppe E. Capuano, Silvano Fineschi, Silvio Giordano, Shadia Habbal, Denise Perrone, Rui F. Pinto, Luca Sorriso-Valvo, Daniele Spadaro, Roberto Susino, Lloyd D. Woodham, Gary P. Zank, Marco Romoli, Stuart D. Bale, Justin C. Kasper, Frédéric Auchère, Roberto Bruno, Gerardo Capobianco, Anthony W. Case, Chiara Casini, Marta Casti, Paolo Chioetto , et al. (46 additional authors not shown)

Abstract: This Letter addresses the first Solar Orbiter (SO) -- Parker Solar Probe (PSP) quadrature, occurring on January 18, 2021, to investigate the evolution of solar wind from the extended corona to the inner heliosphere. Assuming ballistic propagation, the same plasma volume observed remotely in corona at altitudes between 3.5 and 6.3 solar radii above the solar limb with the Metis coronagraph on SO ca… ▽ More This Letter addresses the first Solar Orbiter (SO) -- Parker Solar Probe (PSP) quadrature, occurring on January 18, 2021, to investigate the evolution of solar wind from the extended corona to the inner heliosphere. Assuming ballistic propagation, the same plasma volume observed remotely in corona at altitudes between 3.5 and 6.3 solar radii above the solar limb with the Metis coronagraph on SO can be tracked to PSP, orbiting at 0.1 au, thus allowing the local properties of the solar wind to be linked to the coronal source region from where it originated. Thanks to the close approach of PSP to the Sun and the simultaneous Metis observation of the solar corona, the flow-aligned magnetic field and the bulk kinetic energy flux density can be empirically inferred along the coronal current sheet with an unprecedented accuracy, allowing in particular estimation of the Alfvén radius at 8.7 solar radii during the time of this event. This is thus the very first study of the same solar wind plasma as it expands from the sub-Alfvénic solar corona to just above the Alfvén surface. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 10 pages, 4 figures

Journal ref: Telloni, D., Andretta, V., Antonucci, E., et al. 2021, ApJL, 920, L14

Predicting the Magnetic Fields of a Stealth CME Detected by Parker Solar Probe at 0.5 AU

Authors: Erika Palmerio, Christina Kay, Nada Al-Haddad, Benjamin J. Lynch, Wenyuan Yu, Michael L. Stevens, Sanchita Pal, Christina O. Lee

Abstract: Stealth coronal mass ejection (CMEs) are eruptions from the Sun that are not associated with appreciable low-coronal signatures. Because they often cannot be linked to a well-defined source region on the Sun, analysis of their initial magnetic configuration and eruption dynamics is particularly problematic. In this manuscript, we address this issue by undertaking the first attempt at predicting th… ▽ More Stealth coronal mass ejection (CMEs) are eruptions from the Sun that are not associated with appreciable low-coronal signatures. Because they often cannot be linked to a well-defined source region on the Sun, analysis of their initial magnetic configuration and eruption dynamics is particularly problematic. In this manuscript, we address this issue by undertaking the first attempt at predicting the magnetic fields of a stealth CME that erupted in 2020 June from the Earth-facing Sun. We estimate its source region with the aid of off-limb observations from a secondary viewpoint and photospheric magnetic field extrapolations. We then employ the Open Solar Physics Rapid Ensemble Information (OSPREI) modelling suite to evaluate its early evolution and forward-model its magnetic fields up to Parker Solar Probe, which detected the CME in situ at a heliocentric distance of 0.5 AU. We compare our hindcast prediction with in-situ measurements and a set of flux rope reconstructions, obtaining encouraging agreement on arrival time, spacecraft crossing location, and magnetic field profiles. This work represents a first step towards reliable understanding and forecasting of the magnetic configuration of stealth CMEs and slow, streamer-blowout events. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 9 pages, 4 figures, 1 table, accepted for publication in The Astrophysical Journal

Characteristic scales of magnetic switchback patches near the Sun and their possible association with solar supergranulation and granulation

Authors: Naïs Fargette, Benoit Lavraud, Alexis Rouillard, Victor Réville, Thierry Dudok De Wit, Clara Froment, Jasper S. Halekas, Tai Phan, David Malaspina, Stuart D. Bale, Justin Kasper, Philippe Louarn, Anthony W. Case, Kelly E. Korreck, Davin E. Larson, Marc Pulupa, Michael L. Stevens, Phyllis L. Whittlesey, Matthieu Berthomier

Abstract: Parker Solar Probe (PSP) data recorded within a heliocentric radial distance of 0.3 AU have revealed a magnetic field dominated by Alfvénic structures that undergo large local variations or even reversals of the radial magnetic field. They are called magnetic switchbacks, they are consistent with folds in magnetic field lines within a same magnetic sector, and are associated with velocity spikes d… ▽ More Parker Solar Probe (PSP) data recorded within a heliocentric radial distance of 0.3 AU have revealed a magnetic field dominated by Alfvénic structures that undergo large local variations or even reversals of the radial magnetic field. They are called magnetic switchbacks, they are consistent with folds in magnetic field lines within a same magnetic sector, and are associated with velocity spikes during an otherwise calmer background. They are thought to originate either in the low solar atmosphere through magnetic reconnection processes, or result from the evolution of turbulence or velocity shears in the expanding solar wind. In this work, we investigate the temporal and spatial characteristic scales of magnetic switchback patches. We define switchbacks as a deviation from the nominal Parker spiral direction and detect them automatically for PSP encounters 1, 2, 4 and 5. We focus in particular on a 5.1-day interval dominated by switchbacks during E5. We perform a wavelet transform of the solid angle between the magnetic field and the Parker spiral and find periodic spatial modulations with two distinct wavelengths, respectively consistent with solar granulation and supergranulation scales. In addition we find that switchback occurrence and spectral properties seem to depend on the source region of the solar wind rather than on the radial distance of PSP. These results suggest that switchbacks are formed in the low corona and modulated by the solar surface convection pattern. △ Less

Submitted 13 May, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures

Ambipolar electric field and potential in the solar wind estimated from electron velocity distribution functions

Authors: Laura Bercic, Milan Maksimovic, Jasper S. Halekas, Smone Landi, Christopher J. Owen, Daniel Verscharen, Davin Larson, Phyllis Whittlesey, Samuel T. Badman, Stuart. D. Bale, Anthony W. Case, Keith Goetz, Peter R. Harvey, Justin C. Kasper, Kelly E. Korreck, Roberto Livi, Robert J. MacDowall, David M. Malaspina, Marc Pulupa, Michael L. Stevens

Abstract: The solar wind escapes from the solar corona and is accelerated, over a short distance, to its terminal velocity. The energy balance associated with this acceleration remains poorly understood. To quantify the global electrostatic contribution to the solar wind dynamics, we empirically estimate the ambipolar electric field ($\mathrm{E}_\parallel$) and potential ($Φ_\mathrm{r,\infty}$). We analyse… ▽ More The solar wind escapes from the solar corona and is accelerated, over a short distance, to its terminal velocity. The energy balance associated with this acceleration remains poorly understood. To quantify the global electrostatic contribution to the solar wind dynamics, we empirically estimate the ambipolar electric field ($\mathrm{E}_\parallel$) and potential ($Φ_\mathrm{r,\infty}$). We analyse electron velocity distribution functions (VDFs) measured in the near-Sun solar wind, between 20.3\,$R_S$ and 85.3\,$R_S$, by the Parker Solar Probe. We test the predictions of two different solar wind models. Close to the Sun, the VDFs exhibit a suprathermal electron deficit in the sunward, magnetic field aligned part of phase space. We argue that the sunward deficit is a remnant of the electron cutoff predicted by collisionless exospheric models (Lemaire & Sherer 1970, 1971, Jockers 1970). This cutoff energy is directly linked to $Φ_\mathrm{r,\infty}$. Competing effects of $\mathrm{E}_\parallel$ and Coulomb collisions in the solar wind are addressed by the Steady Electron Runaway Model (SERM) (Scudder 2019). In this model, electron phase space is separated into collisionally overdamped and underdamped regions. We assume that this boundary velocity at small pitch angles coincides with the strahl break-point energy, which allows us to calculate $\mathrm{E}_\parallel$. The obtained $Φ_\mathrm{r,\infty}$ and $\mathrm{E}_\parallel$ agree well with theoretical expectations. They decrease with radial distance as power law functions with indices $α_Φ= -0.66$ and $α_\mathrm{E} = -1.69$. We finally estimate the velocity gained by protons from electrostatic acceleration, which equals to 77\% calculated from the exospheric models, and to 44\% from the SERM model. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Spatial Graph Attention and Curiosity-driven Policy for Antiviral Drug Discovery

Authors: Yulun Wu, Mikaela Cashman, Nicholas Choma, Érica T. Prates, Verónica G. Melesse Vergara, Manesh Shah, Andrew Chen, Austin Clyde, Thomas S. Brettin, Wibe A. de Jong, Neeraj Kumar, Martha S. Head, Rick L. Stevens, Peter Nugent, Daniel A. Jacobson, James B. Brown

Abstract: We developed Distilled Graph Attention Policy Network (DGAPN), a reinforcement learning model to generate novel graph-structured chemical representations that optimize user-defined objectives by efficiently navigating a physically constrained domain. The framework is examined on the task of generating molecules that are designed to bind, noncovalently, to functional sites of SARS-CoV-2 proteins. W… ▽ More We developed Distilled Graph Attention Policy Network (DGAPN), a reinforcement learning model to generate novel graph-structured chemical representations that optimize user-defined objectives by efficiently navigating a physically constrained domain. The framework is examined on the task of generating molecules that are designed to bind, noncovalently, to functional sites of SARS-CoV-2 proteins. We present a spatial Graph Attention (sGAT) mechanism that leverages self-attention over both node and edge attributes as well as encoding the spatial structure -- this capability is of considerable interest in synthetic biology and drug discovery. An attentional policy network is introduced to learn the decision rules for a dynamic, fragment-based chemical environment, and state-of-the-art policy gradient techniques are employed to train the network with stability. Exploration is driven by the stochasticity of the action space design and the innovation reward bonuses learned and proposed by random network distillation. In experiments, our framework achieved outstanding results compared to state-of-the-art algorithms, while reducing the complexity of paths to chemical synthesis. △ Less

Submitted 11 May, 2022; v1 submitted 3 June, 2021; originally announced June 2021.

An Empirical Background Model for the NICER X-ray Timing Instrument

Authors: Ronald A. Remillard, Michael Loewenstein, James F. Steiner, Gregory Y. Prigozhin, Beverly LaMarr, Teruaki Enoto, Keith C. Gendreau, Zaven Arzoumanian, Craig Markwardt, Arkadip Basak, Abigail L. Stevens, Paul S. Ray, Diego Altamirano, Douglas J. K. Buisson

Abstract: NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV;… ▽ More NICER has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4-12 keV; 50 detectors) of 0.87 c/s, but in 5 percent (1 percent) of cases, the rate exceeds 10 (300) c/s. Model residuals persist at 20-30 percent of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c/s (3 sigma, single GTI) at 0.4-12 keV, equivalent to 3.6e-12 erg/cm^2/s for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c/s at 0.3-2.0 keV, or 4.3e-13 erg/cm^2/s for a 100 eV blackbody. Faint-source filtering selects 85 percent of the background GTIs, and higher rates are expected for targets scheduled more favorably. An application of the model to 1 s timescale makes it possible to distinguish source flares from possible surges in the background. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 27 pages, 14 figures; submitted to the Astronomical Journal

A powerful machine learning technique to extract proton core, beam and alpha-particle parameters from velocity distribution functions in space plasmas

Authors: Daniel Vech, Michael L. Stevens, Kristoff W. Paulson, David M. Malaspina, Anthony W. Case, Kristopher G. Klein, Justin C. Kasper

Abstract: Context: The analysis of the thermal part of velocity distribution functions (VDF) is fundamentally important for understanding the kinetic physics that governs the evolution and dynamics of space plasmas. However, calculating the proton core, beam and alpha-particle parameters for large data sets of VDFs is a time consuming and computationally demanding process that always requires supervision by… ▽ More Context: The analysis of the thermal part of velocity distribution functions (VDF) is fundamentally important for understanding the kinetic physics that governs the evolution and dynamics of space plasmas. However, calculating the proton core, beam and alpha-particle parameters for large data sets of VDFs is a time consuming and computationally demanding process that always requires supervision by a human expert. Aims: We developed a machine learning tool that can extract proton core, beam and alpha-particle parameters using images (2-D grid consisting pixel values) of VDFs. Methods: A database of synthetic VDFs is generated, which is used to train a convolutional neural network that infers bulk speed, thermal speed and density for all three particle populations. We generate a separate test data set of synthetic VDFs that we use to compare and quantify the predictive power of the neural network and a fitting algorithm. Results: The neural network achieves significantly smaller root-mean-square errors to infer proton core, beam and alpha-particle parameters than a traditional fitting algorithm. Conclusion: The developed machine learning tool has the potential to revolutionize the processing of particle measurements since it allows the computation of more accurate particle parameters than previously used fitting procedures. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: Accepted in Astronomy and Astrophysics

A Verified Decision Procedure for Orders in Isabelle/HOL

Authors: Lukas Stevens, Tobias Nipkow

Abstract: We present the first verified implementation of a decision procedure for the quantifier-free theory of partial and linear orders. We formalise the procedure in Isabelle/HOL and provide a specification that is made executable using Isabelle's code generator. The procedure is already part of the development version of Isabelle as a sub-procedure of the simplifier. We present the first verified implementation of a decision procedure for the quantifier-free theory of partial and linear orders. We formalise the procedure in Isabelle/HOL and provide a specification that is made executable using Isabelle's code generator. The procedure is already part of the development version of Isabelle as a sub-procedure of the simplifier. △ Less

Submitted 2 July, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Using Parker Solar Probe observations during the first four perihelia to constrain global magnetohydrodynamic models

Authors: Pete Riley, Roberto Lionello, Ronald M. Caplan, Cooper Downs, Jon A. Linker, Samuel T. Badman, Michael L. Stevens

Abstract: Parker Solar Probe (PSP) is providing an unprecedented view of the Sun's corona as it progressively dips closer into the solar atmosphere with each solar encounter. Each set of observations provides a unique opportunity to test and constrain global models of the solar corona and inner heliosphere and, in turn, use the model results to provide a global context for interpreting such observations. In… ▽ More Parker Solar Probe (PSP) is providing an unprecedented view of the Sun's corona as it progressively dips closer into the solar atmosphere with each solar encounter. Each set of observations provides a unique opportunity to test and constrain global models of the solar corona and inner heliosphere and, in turn, use the model results to provide a global context for interpreting such observations. In this study, we develop a set of global magnetohydrodynamic (MHD) model solutions of varying degrees of sophistication for PSP's first four encounters and compare the results with in situ measurements from PSP, Stereo-A, and Earth-based spacecraft, with the objective of assessing which models perform better or worse. All models were primarily driven by the observed photospheric magnetic field using data from Solar Dynamics Observatory's Helioseismic and Magnetic Imager (HMI) instrument. Overall, we find that there are substantial differences between the model results, both in terms of the large-scale structure of the inner heliosphere during these time periods, as well as in the inferred time-series at various spacecraft. The "thermodynamic" model, which represents the "middle ground", in terms of model complexity, appears to reproduce the observations most closely for all four encounters. Our results also contradict an earlier study that had hinted that the open flux problem may disappear nearer the Sun. Instead, our results suggest that this "missing" solar flux is still missing even at 26.9 Rs, and thus it cannot be explained by interplanetary processes. Finally, the model results were also used to provide a global context for interpreting the localized in situ measurements. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Journal ref: A&A 650, A19 (2021)

Determination of Solar Wind Angular Momentum and Alfvén Radius from Parker Solar Probe Observations

Authors: Ying D. Liu, Chong Chen, Michael L. Stevens, Mingzhe Liu

Abstract: As fundamental parameters of the Sun, the Alfvén radius and angular momentum loss determine how the solar wind changes from sub-Alfvénic to super-Alfvénic and how the Sun spins down. We present an approach to determining the solar wind angular momentum flux based on observations from Parker Solar Probe (PSP). A flux of about $0.15\times10^{30}$ dyn cm sr$^{-1}$ near the ecliptic plane and 0.7:1 pa… ▽ More As fundamental parameters of the Sun, the Alfvén radius and angular momentum loss determine how the solar wind changes from sub-Alfvénic to super-Alfvénic and how the Sun spins down. We present an approach to determining the solar wind angular momentum flux based on observations from Parker Solar Probe (PSP). A flux of about $0.15\times10^{30}$ dyn cm sr$^{-1}$ near the ecliptic plane and 0.7:1 partition of that flux between the particles and magnetic field are obtained by averaging data from the first four encounters within 0.3 au from the Sun. The angular momentum flux and its particle component decrease with the solar wind speed, while the flux in the field is remarkably constant. A speed dependence in the Alfvén radius is also observed, which suggests a "rugged" Alfvén surface around the Sun. Substantial diving below the Alfvén surface seems plausible only for relatively slow solar wind given the orbital design of PSP. Uncertainties are evaluated based on the acceleration profiles of the same solar wind streams observed at PSP and a radially aligned spacecraft near 1 au. We illustrate that the "angular momentum paradox" raised by Réville et al. can be removed by taking into account the contribution of the alpha particles. The large proton transverse velocity observed by PSP is perhaps inherent in the solar wind acceleration process, where an opposite transverse velocity is produced for the alphas with the angular momentum conserved. Preliminary analysis of some recovered alpha parameters tends to agree with the results. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: Accepted for publication in The Astrophysical Journal Letters