-
The Complex Interplay Between Risk Tolerance and the Spread of Infectious Diseases
Authors:
Maximilian Nguyen,
Ari Freedman,
Matthew Cheung,
Chadi Saad-Roy,
Baltazar Espinoza,
Bryan Grenfell,
Simon Levin
Abstract:
Risk-driven behavior provides a feedback mechanism through which individuals both shape and are collectively affected by an epidemic. We introduce a general and flexible compartmental model to study the effect of heterogeneity in the population with regards to risk tolerance. The interplay between behavior and epidemiology leads to a rich set of possible epidemic dynamics. Depending on the behavio…
▽ More
Risk-driven behavior provides a feedback mechanism through which individuals both shape and are collectively affected by an epidemic. We introduce a general and flexible compartmental model to study the effect of heterogeneity in the population with regards to risk tolerance. The interplay between behavior and epidemiology leads to a rich set of possible epidemic dynamics. Depending on the behavioral composition of the population, we find that increasing heterogeneity in risk tolerance can either increase or decrease the epidemic size. We find that multiple waves of infection can arise due to the interplay between transmission and behavior, even without the replenishment of susceptibles. We find that increasing protective mechanisms such as the effectiveness of interventions, the number of risk-averse people in the population, and the duration of intervention usage reduces the epidemic overshoot. When the protection is pushed past a critical threshold, the epidemic dynamics enter an underdamped regime where the epidemic size exactly equals the herd immunity threshold. Lastly, we can find regimes where epidemic size does not monotonically decrease with a population that becomes increasingly risk-averse.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Metric-guided Image Reconstruction Bounds via Conformal Prediction
Authors:
Matt Y Cheung,
Tucker J Netherton,
Laurence E Court,
Ashok Veeraraghavan,
Guha Balakrishnan
Abstract:
Recent advancements in machine learning have led to the development of novel medical imaging systems and algorithms that address ill-posed problems. Assessing their trustworthiness and understanding how to deploy them safely at test time remains an important and open problem. In this work, we propose using conformal prediction to compute valid and distribution-free bounds on downstream metrics giv…
▽ More
Recent advancements in machine learning have led to the development of novel medical imaging systems and algorithms that address ill-posed problems. Assessing their trustworthiness and understanding how to deploy them safely at test time remains an important and open problem. In this work, we propose using conformal prediction to compute valid and distribution-free bounds on downstream metrics given reconstructions generated by one algorithm, and retrieve upper/lower bounds and inlier/outlier reconstructions according to the adjusted bounds. Our work offers 1) test time image reconstruction evaluation without ground truth, 2) downstream performance guarantees, 3) meaningful upper/lower bound reconstructions, and 4) meaningful statistical inliers/outlier reconstructions. We demonstrate our method on post-mastectomy radiotherapy planning using 3D breast CT reconstructions, and show 1) that metric-guided bounds have valid coverage for downstream metrics while conventional pixel-wise bounds do not and 2) anatomical differences of upper/lower bounds between metric-guided and pixel-wise methods. Our work paves way for more meaningful and trustworthy test-time evaluation of medical image reconstructions. Code available at https://github.com/matthewyccheung/conformal-metric
△ Less
Submitted 1 July, 2024; v1 submitted 23 April, 2024;
originally announced April 2024.
-
Resolving moving heliospheric structures using interplanetary scintillation observations with the Murchison Widefield Array
Authors:
A. Waszewski,
J. S. Morgan,
R. Chhetri,
R. Ekers,
M. C. M. Cheung,
N. D. R Bhat,
M. Johnston-Hollitt
Abstract:
We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations…
▽ More
We have conducted a blind search in 49 consecutive days of interplanetary scintillation observations made by the Murchison Widefield Array from mid-2019, with overlapping daily observations approximately East and South-East of the Sun at an elongation of $\sim$30 degrees and a field of view of 30 degrees. These observations detect an unprecedented density of sources. In spite of these observations being taken at sunspot minimum, this search has revealed several interesting transitory features characterised by elevated scintillation levels. One solar wind enhancement is captured in two observations several hours apart, allowing its radial movement away from the Sun to be measured. We present here a methodology for measuring the plane-of-sky velocity for the moving heliospheric structure. The plane-of-sky velocity was inferred as $0.66\pm0.147\,^{\text{o}}\text{hr}^{-1}$, or $480\pm106\,\text{km}\,\text{s}^{-1}$ assuming a distance of 1AU. After cross-referencing our observed structure with multiple catalogues of heliospheric events, we propose that the likely source of our observed structure is a stream-interaction region originating from a low-latitude coronal hole. This work demonstrates the power of widefield interplanetary scintillation observations to capture detailed features in the heliosphere which are otherwise unresolvable and go undetected.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
Particle acceleration in solar flares with imaging-spectroscopy in soft X-rays
Authors:
Mitsuo Oka,
Amir Caspi,
Bin Chen,
Mark Cheung,
James Drake,
Dale Gary,
Lindsay Glesener,
Fan Guo,
Hantao Ji,
Xiaocan Li,
Takuma Nakamura,
Noriyuki Narukage,
Katharine Reeves,
Pascal Saint-Hilaire,
Taro Sakao,
Chengcai Shen,
Amy Winebarger,
Tom Woods
Abstract:
Particles are accelerated to very high, non-thermal energies during explosive energy-release phenomena in space, solar, and astrophysical plasma environments. In the case of solar flares, it has been established that magnetic reconnection plays an important role for releasing the magnetic energy, but it remains unclear if or how magnetic reconnection can further explain particle acceleration durin…
▽ More
Particles are accelerated to very high, non-thermal energies during explosive energy-release phenomena in space, solar, and astrophysical plasma environments. In the case of solar flares, it has been established that magnetic reconnection plays an important role for releasing the magnetic energy, but it remains unclear if or how magnetic reconnection can further explain particle acceleration during flares. Here we argue that the key issue is the lack of understanding of the precise context of particle acceleration but it can be overcome, in the near future, by performing imaging-spectroscopy in soft X-rays (SXRs). Such observations should be complemented by observations in other wavelengths such as extreme-ultraviolets (EUVs), microwaves, hard X-rays (HXRs), and gamma-rays. Also, numerical simulations will be crucial for further narrowing down the particle acceleration mechanism in the context revealed by the observations. Of all these efforts, imaging-spectroscopy in SXRs, if successfully applied to large limb flares, will be a milestone in our challenge of understanding electron acceleration in solar flares and beyond, i.e. the Plasma Universe.
△ Less
Submitted 7 June, 2023;
originally announced June 2023.
-
Potential quantum advantage for simulation of fluid dynamics
Authors:
Xiangyu Li,
Xiaolong Yin,
Nathan Wiebe,
Jaehun Chun,
Gregory K. Schenter,
Margaret S. Cheung,
Johannes Mülmenstädt
Abstract:
Numerical simulation of turbulent fluid dynamics needs to either parameterize turbulence-which introduces large uncertainties-or explicitly resolve the smallest scales-which is prohibitively expensive. Here we provide evidence through analytic bounds and numerical studies that a potential quantum exponential speedup can be achieved to simulate the Navier-Stokes equations governing turbulence using…
▽ More
Numerical simulation of turbulent fluid dynamics needs to either parameterize turbulence-which introduces large uncertainties-or explicitly resolve the smallest scales-which is prohibitively expensive. Here we provide evidence through analytic bounds and numerical studies that a potential quantum exponential speedup can be achieved to simulate the Navier-Stokes equations governing turbulence using quantum computing. Specifically, we provide a formulation of the lattice Boltzmann equation for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems and that for computationally interesting examples. This is achieved via a combination of reformulating the nonlinearity and accurately linearizing the dynamical equations, effectively trading nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size, compared to polynomial scaling in the best known classical algorithms. This work suggests that an exponential quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.
△ Less
Submitted 28 March, 2024; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling
Authors:
Graham Barnes,
Marc L. DeRosa,
Shaela I. Jones,
Charles N. Arge,
Carl J. Henney,
Mark C. M. Cheung
Abstract:
The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open…
▽ More
The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
Science Platforms for Heliophysics Data Analysis
Authors:
Monica G. Bobra,
Will T. Barnes,
Thomas Y. Chen,
Mark C. M. Cheung,
Laura A. Hayes,
Jack Ireland,
Miho Janvier,
Michael S. F. Kirk,
James P. Mason,
Stuart J. Mumford,
Paul J. Wright
Abstract:
We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments.
We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments.
△ Less
Submitted 2 January, 2023;
originally announced January 2023.
-
Towards data-driven modeling and real-time prediction of solar flares and coronal mass ejections
Authors:
M. Rempel,
Y. Fan,
M. Dikpati,
A. Malanushenko,
M. D. Kazachenko,
M. C. M. Cheung,
G. Chintzoglou,
X. Sun,
G. H. Fisher,
T. Y. Chen
Abstract:
Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientatio…
▽ More
Modeling of transient events in the solar atmosphere requires the confluence of 3 critical elements: (1) model sophistication, (2) data availability, and (3) data assimilation. This white paper describes required advances that will enable statistical flare and CME forecasting (e.g. eruption probability and timing, estimation of strength, and CME details, such as speed and magnetic field orientation) similar to weather prediction on Earth.
△ Less
Submitted 29 December, 2022;
originally announced December 2022.
-
Solar coronal heating from small-scale magnetic braids
Authors:
L. P. Chitta,
H. Peter,
S. Parenti,
D. Berghmans,
F. Auchère,
S. K. Solanki,
R. Aznar Cuadrado,
U. Schühle,
L. Teriaca,
S. Mandal,
K. Barczynski,
É. Buchlin,
L. Harra,
E. Kraaikamp,
D. M. Long,
L. Rodriguez,
C. Schwanitz,
P. J. Smith,
C. Verbeeck,
A. N. Zhukov,
W. Liu,
M. C. M. Cheung
Abstract:
Relaxation of braided coronal magnetic fields through reconnection is thought to be a source of energy to heat plasma in active region coronal loops. However, observations of active region coronal heating associated with an untangling of magnetic braids remain sparse. One reason for this paucity could be the lack of coronal observations with a sufficiently high spatial and temporal resolution to c…
▽ More
Relaxation of braided coronal magnetic fields through reconnection is thought to be a source of energy to heat plasma in active region coronal loops. However, observations of active region coronal heating associated with an untangling of magnetic braids remain sparse. One reason for this paucity could be the lack of coronal observations with a sufficiently high spatial and temporal resolution to capture this process in action. Using new observations with high spatial resolution (250-270 km on the Sun) and high cadence (3-10 s) from the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter, we observed the untangling of small-scale coronal braids in different active regions. The untangling is associated with impulsive heating of the gas in these braided loops. We assess that coronal magnetic braids overlying cooler chromospheric filamentary structures are perhaps more common. Furthermore, our observations show signatures of spatially coherent and intermittent coronal heating during the relaxation of the magnetic braids. Our study reveals the operation of gentle and impulsive modes of magnetic reconnection in the solar corona.
△ Less
Submitted 26 November, 2022; v1 submitted 25 September, 2022;
originally announced September 2022.
-
Global geomagnetic perturbation forecasting using Deep Learning
Authors:
Vishal Upendran,
Panagiotis Tigas,
Banafsheh Ferdousi,
Teo Bloch,
Mark C. M. Cheung,
Siddha Ganju,
Asti Bhatt,
Ryan M. McGranaghan,
Yarin Gal
Abstract:
Geomagnetically Induced Currents (GICs) arise from spatio-temporal changes to Earth's magnetic field which arise from the interaction of the solar wind with Earth's magnetosphere, and drive catastrophic destruction to our technologically dependent society. Hence, computational models to forecast GICs globally with large forecast horizon, high spatial resolution and temporal cadence are of increasi…
▽ More
Geomagnetically Induced Currents (GICs) arise from spatio-temporal changes to Earth's magnetic field which arise from the interaction of the solar wind with Earth's magnetosphere, and drive catastrophic destruction to our technologically dependent society. Hence, computational models to forecast GICs globally with large forecast horizon, high spatial resolution and temporal cadence are of increasing importance to perform prompt necessary mitigation. Since GIC data is proprietary, the time variability of horizontal component of the magnetic field perturbation (dB/dt) is used as a proxy for GICs. In this work, we develop a fast, global dB/dt forecasting model, which forecasts 30 minutes into the future using only solar wind measurements as input. The model summarizes 2 hours of solar wind measurement using a Gated Recurrent Unit, and generates forecasts of coefficients which are folded with a spherical harmonic basis to enable global forecasts. When deployed, our model produces results in under a second, and generates global forecasts for horizontal magnetic perturbation components at 1-minute cadence. We evaluate our model across models in literature for two specific storms of 5 August 2011 and 17 March 2015, while having a self-consistent benchmark model set. Our model outperforms, or has consistent performance with state-of-the-practice high time cadence local and low time cadence global models, while also outperforming/having comparable performance with the benchmark models. Such quick inferences at high temporal cadence and arbitrary spatial resolutions may ultimately enable accurate forewarning of dB/dt for any place on Earth, resulting in precautionary measures to be taken in an informed manner.
△ Less
Submitted 12 May, 2022;
originally announced May 2022.
-
Plasma Heating Induced by Tadpole-Like Downflows in the Flaring Solar Corona
Authors:
Tanmoy Samanta,
Hui Tian,
Bin Chen,
Katharine K. Reeves,
Mark C. M. Cheung,
Angelos Vourlidas,
Dipankar Banerjee
Abstract:
As one of the most spectacular energy release events in the solar system, solar flares are generally powered by magnetic reconnection in the solar corona. As a result of the re-arrangement of magnetic field topology after the reconnection process, a series of new loop-like magnetic structures are often formed and are known as flare loops. A hot diffuse region, consisting of around 5-10 MK plasma,…
▽ More
As one of the most spectacular energy release events in the solar system, solar flares are generally powered by magnetic reconnection in the solar corona. As a result of the re-arrangement of magnetic field topology after the reconnection process, a series of new loop-like magnetic structures are often formed and are known as flare loops. A hot diffuse region, consisting of around 5-10 MK plasma, is also observed above the loops and is called a supra-arcade fan. Often, dark, tadpole-like structures are seen to descend through the bright supra-arcade fans. It remains unclear what role these so-called supra-arcade downflows (SADs) play in heating the flaring coronal plasma. Here we show a unique flare observation, where many SADs collide with the flare loops and strongly heat the loops to a temperature of 10-20 MK. Several of these interactions generate clear signatures of quasi-periodic enhancement in the full-Sun-integrated soft X-ray emission, providing an alternative interpretation for quasi-periodic pulsations that are commonly observed during solar and stellar flares.
△ Less
Submitted 26 March, 2021;
originally announced March 2021.
-
Global Earth Magnetic Field Modeling and Forecasting with Spherical Harmonics Decomposition
Authors:
Panagiotis Tigas,
Téo Bloch,
Vishal Upendran,
Banafsheh Ferdoushi,
Mark C. M. Cheung,
Siddha Ganju,
Ryan M. McGranaghan,
Yarin Gal,
Asti Bhatt
Abstract:
Modeling and forecasting the solar wind-driven global magnetic field perturbations is an open challenge. Current approaches depend on simulations of computationally demanding models like the Magnetohydrodynamics (MHD) model or sampling spatially and temporally through sparse ground-based stations (SuperMAG). In this paper, we develop a Deep Learning model that forecasts in Spherical Harmonics spac…
▽ More
Modeling and forecasting the solar wind-driven global magnetic field perturbations is an open challenge. Current approaches depend on simulations of computationally demanding models like the Magnetohydrodynamics (MHD) model or sampling spatially and temporally through sparse ground-based stations (SuperMAG). In this paper, we develop a Deep Learning model that forecasts in Spherical Harmonics space 2, replacing reliance on MHD models and providing global coverage at one minute cadence, improving over the current state-of-the-art which relies on feature engineering. We evaluate the performance in SuperMAG dataset (improved by 14.53%) and MHD simulations (improved by 24.35%). Additionally, we evaluate the extrapolation performance of the spherical harmonics reconstruction based on sparse ground-based stations (SuperMAG), showing that spherical harmonics can reliably reconstruct the global magnetic field as evaluated on MHD simulation.
△ Less
Submitted 2 February, 2021;
originally announced February 2021.
-
Multi-Channel Auto-Calibration for the Atmospheric Imaging Assembly using Machine Learning
Authors:
Luiz F. G. dos Santos,
Souvik Bose,
Valentina Salvatelli,
Brad Neuberg,
Mark C. M. Cheung,
Miho Janvier,
Meng Jin,
Yarin Gal,
Paul Boerner,
Atılım Güneş Baydin
Abstract:
Solar activity plays a quintessential role in influencing the interplanetary medium and space-weather around the Earth. Remote sensing instruments onboard heliophysics space missions provide a pool of information about the Sun's activity via the measurement of its magnetic field and the emission of light from the multi-layered, multi-thermal, and dynamic solar atmosphere. Extreme UV (EUV) waveleng…
▽ More
Solar activity plays a quintessential role in influencing the interplanetary medium and space-weather around the Earth. Remote sensing instruments onboard heliophysics space missions provide a pool of information about the Sun's activity via the measurement of its magnetic field and the emission of light from the multi-layered, multi-thermal, and dynamic solar atmosphere. Extreme UV (EUV) wavelength observations from space help in understanding the subtleties of the outer layers of the Sun, namely the chromosphere and the corona. Unfortunately, such instruments, like the Atmospheric Imaging Assembly (AIA) onboard NASA's Solar Dynamics Observatory (SDO), suffer from time-dependent degradation, reducing their sensitivity. Current state-of-the-art calibration techniques rely on periodic sounding rockets, which can be infrequent and rather unfeasible for deep-space missions. We present an alternative calibration approach based on convolutional neural networks (CNNs). We use SDO-AIA data for our analysis. Our results show that CNN-based models could comprehensively reproduce the sounding rocket experiments' outcomes within a reasonable degree of accuracy, indicating that it performs equally well compared with the current techniques. Furthermore, a comparison with a standard "astronomer's technique" baseline model reveals that the CNN approach significantly outperforms this baseline. Our approach establishes the framework for a novel technique to calibrate EUV instruments and advance our understanding of the cross-channel relation between different EUV channels.
△ Less
Submitted 1 February, 2021; v1 submitted 27 December, 2020;
originally announced December 2020.
-
Insights from Graph Theory on the Morphologies of Actomyosin Networks with Multilinkers
Authors:
Yossi Eliaz,
Francois Nedelec,
Greg Morrison,
Herbert Levine,
Margaret S. Cheung
Abstract:
Quantifying the influence of microscopic details on the dynamics of development of the overall structure of a filamentous network is important in a number of biologically relevant contexts, but it is not obvious what order parameters can be used to adequately describe this complex process. In this paper, we investigated the role of multivalent actin-binding proteins (ABPs) in reorganizing actin fi…
▽ More
Quantifying the influence of microscopic details on the dynamics of development of the overall structure of a filamentous network is important in a number of biologically relevant contexts, but it is not obvious what order parameters can be used to adequately describe this complex process. In this paper, we investigated the role of multivalent actin-binding proteins (ABPs) in reorganizing actin filaments into higher-order complex networks via a computer model of semiflexible filaments. We characterize the importance of local connectivity among actin filaments as well as the global features of actomyosin networks. We first map the networks into local graph representations and then, using principles from network-theory order parameters, combine properties from these representations to gain insight on the heterogeneous morphologies of actomyosin networks at a global level. We find that ABPs with a valency greater than two promote filament bundles and large filament clusters to a much greater extent than bivalent multilinkers. We also show that active myosin-like motor proteins promote the formation of dendritic branches from a stalk of actin bundles. Our work motivates future studies to embrace network theory as a tool to characterize complex morphologies of actomyosin detected by experiments, leading to a quantitative understanding of the role of ABPs in manipulating the self-assembly of actin filaments into unique architectures that underlie the structural scaffold of a cell relating to its mobility and shape.
△ Less
Submitted 16 October, 2020; v1 submitted 11 June, 2020;
originally announced June 2020.
-
Coronal Dimming as a Proxy for Stellar Coronal Mass Ejections
Authors:
Meng Jin,
Mark C. M. Cheung,
Marc L. DeRosa,
Nariaki V. Nitta,
Carolus J. Schrijver,
Kevin France,
Adam Kowalski,
James P. Mason,
Rachel Osten
Abstract:
Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like…
▽ More
Solar coronal dimmings have been observed extensively in the past two decades and are believed to have close association with coronal mass ejections (CMEs). Recent study found that coronal dimming is the only signature that could differentiate powerful ares that have CMEs from those that do not. Therefore, dimming might be one of the best candidates to observe the stellar CMEs on distant Sun-like stars. In this study, we investigate the possibility of using coronal dimming as a proxy to diagnose stellar CMEs. By simulating a realistic solar CME event and corresponding coronal dimming using a global magnetohydrodynamics model (AWSoM: Alfven-wave Solar Model), we first demonstrate the capability of the model to reproduce solar observations. We then extend the model for simulating stellar CMEs by modifying the input magnetic flux density as well as the initial magnetic energy of the CME flux rope. Our result suggests that with improved instrument sensitivity, it is possible to detect the coronal dimming signals induced by the stellar CMEs.
△ Less
Submitted 14 February, 2020;
originally announced February 2020.
-
Auto-Calibration of Remote Sensing Solar Telescopes with Deep Learning
Authors:
Brad Neuberg,
Souvik Bose,
Valentina Salvatelli,
Luiz F. G. dos Santos,
Mark Cheung,
Miho Janvier,
Atilim Gunes Baydin,
Yarin Gal,
Meng Jin
Abstract:
As a part of NASA's Heliophysics System Observatory (HSO) fleet of satellites,the Solar Dynamics Observatory (SDO) has continuously monitored the Sun since2010. Ultraviolet (UV) and Extreme UV (EUV) instruments in orbit, such asSDO's Atmospheric Imaging Assembly (AIA) instrument, suffer time-dependent degradation which reduces instrument sensitivity. Accurate calibration for (E)UV instruments curr…
▽ More
As a part of NASA's Heliophysics System Observatory (HSO) fleet of satellites,the Solar Dynamics Observatory (SDO) has continuously monitored the Sun since2010. Ultraviolet (UV) and Extreme UV (EUV) instruments in orbit, such asSDO's Atmospheric Imaging Assembly (AIA) instrument, suffer time-dependent degradation which reduces instrument sensitivity. Accurate calibration for (E)UV instruments currently depends on periodic sounding rockets, which are infrequent and not practical for heliophysics missions in deep space. In the present work, we develop a Convolutional Neural Network (CNN) that auto-calibrates SDO/AIA channels and corrects sensitivity degradation by exploiting spatial patterns in multi-wavelength observations to arrive at a self-calibration of (E)UV imaging instruments. Our results remove a major impediment to developing future HSOmissions of the same scientific caliber as SDO but in deep space, able to observe the Sun from more vantage points than just SDO's current geosynchronous orbit.This approach can be adopted to perform autocalibration of other imaging systems exhibiting similar forms of degradation
△ Less
Submitted 10 November, 2019;
originally announced November 2019.
-
Using U-Nets to Create High-Fidelity Virtual Observations of the Solar Corona
Authors:
Valentina Salvatelli,
Souvik Bose,
Brad Neuberg,
Luiz F. G. dos Santos,
Mark Cheung,
Miho Janvier,
Atilim Gunes Baydin,
Yarin Gal,
Meng Jin
Abstract:
Understanding and monitoring the complex and dynamic processes of the Sun is important for a number of human activities on Earth and in space. For this reason, NASA's Solar Dynamics Observatory (SDO) has been continuously monitoring the multi-layered Sun's atmosphere in high-resolution since its launch in 2010, generating terabytes of observational data every day. The synergy between machine learn…
▽ More
Understanding and monitoring the complex and dynamic processes of the Sun is important for a number of human activities on Earth and in space. For this reason, NASA's Solar Dynamics Observatory (SDO) has been continuously monitoring the multi-layered Sun's atmosphere in high-resolution since its launch in 2010, generating terabytes of observational data every day. The synergy between machine learning and this enormous amount of data has the potential, still largely unexploited, to advance our understanding of the Sun and extend the capabilities of heliophysics missions. In the present work, we show that deep learning applied to SDO data can be successfully used to create a high-fidelity virtual telescope that generates synthetic observations of the solar corona by image translation. Towards this end we developed a deep neural network, structured as an encoder-decoder with skip connections (U-Net), that reconstructs the Sun's image of one instrument channel given temporally aligned images in three other channels. The approach we present has the potential to reduce the telemetry needs of SDO, enhance the capabilities of missions that have less observing channels, and transform the concept development of future missions.
△ Less
Submitted 10 November, 2019;
originally announced November 2019.
-
Critical phenomena in the temperature-pressure-crowding phase diagram of a protein
Authors:
Andrei G. Gasic,
Mayank M. Boob,
Maxim B. Prigozhin,
Dirar Homouz,
Caleb M. Daugherty,
Martin Gruebele,
Margaret S. Cheung
Abstract:
In the cell, proteins fold and perform complex functions through global structural rearrangements. Function requires a protein to be at the brink of stability to be susceptible to small environmental fluctuations, yet stable enough to maintain structural integrity. These apparently conflicting behaviors are exhibited by systems near a critical point, where distinct phases merge $-$ a concept beyon…
▽ More
In the cell, proteins fold and perform complex functions through global structural rearrangements. Function requires a protein to be at the brink of stability to be susceptible to small environmental fluctuations, yet stable enough to maintain structural integrity. These apparently conflicting behaviors are exhibited by systems near a critical point, where distinct phases merge $-$ a concept beyond previous studies indicating proteins have a well-defined folded/unfolded phase boundary in the pressure-temperature plane. Here, by modeling the protein phosphoglycerate kinase (PGK) on the temperature (T), pressure (P), and crowding volume-fraction ($φ$) phase diagram, we demonstrate a critical transition where phases merge, and PGK exhibits large structural fluctuations. Above the critical temperature (Tc), the difference between the intermediate and unfolded phases disappears. When $φ$ increases, the Tc moves to a lower T. We verify the calculations with experiments mapping the T-P-$φ$ space, which likewise reveal a critical point at 305 K and 170 MPa that moves to a lower T as $φ$ increases. Crowding places PGK near a critical line in its natural parameter space, where large conformational changes can occur without costly free energy barriers. Specific structures are proposed for each phase based on simulation.
△ Less
Submitted 9 June, 2019;
originally announced June 2019.
-
The impact of hydrodynamic interactions on protein folding rates depends on temperature
Authors:
Fabio C. Zegarra,
Dirar Homouz,
Yossi Eliaz,
Andrei G. Gasic,
Margaret S. Cheung
Abstract:
We investigated the impact of hydrodynamic interactions (HI) on protein folding using a coarse-grained model. The extent of the impact of hydrodynamic interactions, whether it accelerates, retards, or has no effect on protein folding, has been controversial. Together with a theoretical framework of the energy landscape theory (ELT) for protein folding that describes the dynamics of the collective…
▽ More
We investigated the impact of hydrodynamic interactions (HI) on protein folding using a coarse-grained model. The extent of the impact of hydrodynamic interactions, whether it accelerates, retards, or has no effect on protein folding, has been controversial. Together with a theoretical framework of the energy landscape theory (ELT) for protein folding that describes the dynamics of the collective motion with a single reaction coordinate across a folding barrier, we compared the kinetic effects of HI on the folding rates of two protein models that use a chain of single beads with distinctive topologies: a 64-residue alpha/beta chymotrypsin inhibitor 2 (CI2) protein, and a 57-residue beta-barrel alpha-spectrin src-Homology 3 domain (SH3) protein. When comparing the protein folding kinetics simulated with Brownian dynamics in the presence of HI to that in the absence of HI, we find that the effect of HI on protein folding appears to have a crossover behavior about the folding temperature. Meaning that at a temperature greater than the folding temperature, the enhanced friction from the hydrodynamic solvents between the beads in an unfolded configuration results in lowered folding rate; conversely, at a temperature lower than the folding temperature, HI accelerates folding by the backflow of solvent toward the native folded state. Additionally, the extent of acceleration depends on the topology of a protein: for a protein like CI2, where its folding nucleus is rather diffuse in a transition state, HI channels the formation of contacts by favoring a major folding pathway in a complex free energy landscape, thus accelerating folding. For a protein like SH3, where its folding nucleus is already specific and less diffuse, HI matters less at a temperature lower than the folding temperature. Our findings provide further theoretical insight to protein folding kinetic experiments and simulations.
△ Less
Submitted 12 December, 2017; v1 submitted 11 December, 2017;
originally announced December 2017.
-
Multi-Parametric Study of Rising 3D Buoyant Flux Tubes in an Adiabatic Stratification Using AMR
Authors:
Juan Martinez-Sykora,
Fernando Moreno-Insertis,
Mark C. M. Cheung
Abstract:
We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh…
▽ More
We study the buoyant rise of magnetic flux tubes embedded in an adiabatic stratification using two-and three-dimensional, MHD simulations. We analyze the dependence of the tube evolution on the field line twist and on the curvature of the tube axis in different diffusion regimes. To be able to achieve a comparatively high spatial resolution we use the FLASH code, which has a built-in Adaptive Mesh Refinement (AMR) capability. Our 3D experiments reach Reynolds numbers that permit a reasonable comparison of the results with those of previous 2D simulations. When the experiments are run without AMR, hence with a comparatively large diffusivity, the amount of longitudinal magnetic flux retained inside the tube increases with the curvature of the tube axis. However, when a low-diffusion regime is reached by using the AMR algorithms, the magnetic twist is able to prevent the splitting of the magnetic loop into vortex tubes and the loop curvature does not play any significant role. We detect the generation of vorticity in the main body of the tube of opposite sign on the opposite sides of the apex. This is a consequence of the inhomogeneity of the azimuthal component of the field on the flux surfaces. The lift force associated with this global vorticity makes the flanks of the tube move away from their initial vertical plane in an antisymmetric fashion. The trajectories have an oscillatory motion superimposed, due to the shedding of vortex rolls to the wake, which creates a Von Karman street.
△ Less
Submitted 6 July, 2015;
originally announced July 2015.
-
160 MeV laser-accelerated protons from CH2 nano-targets for proton cancer therapy
Authors:
B. M. Hegelich,
D. Jung,
B. J. Albright,
M. Cheung,
B. Dromey,
D. C. Gautier,
C. Hamilton,
S. Letzring,
R. Munchhausen,
S. Palaniyappan,
R. Shah,
H. -C. Wu,
L. Yin,
J. C. Fernández
Abstract:
Proton (and ion) cancer therapy has proven to be an extremely effective even supe-rior method of treatment for some tumors 1-4. A major problem, however, lies in the cost of the particle accelerator facilities; high procurement costs severely limit the availability of ion radiation therapy, with only ~26 centers worldwide. Moreover, high operating costs often prevent economic operation without sta…
▽ More
Proton (and ion) cancer therapy has proven to be an extremely effective even supe-rior method of treatment for some tumors 1-4. A major problem, however, lies in the cost of the particle accelerator facilities; high procurement costs severely limit the availability of ion radiation therapy, with only ~26 centers worldwide. Moreover, high operating costs often prevent economic operation without state subsidies and have led to a shutdown of existing facilities 5,6. Laser-accelerated proton and ion beams have long been thought of as a way out of this dilemma, with the potential to provide the required ion beams at lower cost and smaller facility footprint 7-14. The biggest challenge has been the achievement of sufficient particle energy for therapy, in the 150-250 MeV range for protons 15,16. For the last decade, the maximum exper-imentally observed energy of laser-accelerated protons has remained at ~60 MeV 17. Here we the experimental demonstration of laser-accelerated protons to energies exceeding 150 MeV, reaching the therapy window. This was achieved through a dif-ferent acceleration regime rather than a larger laser, specifically a 150 TW laser with CH2 nano-targets in the relativistically transparent regime 18,19. We also demonstrate a clear scaling law with laser intensity based on analytical theory, computer simulations and experimental validation that will enable design of a pro-totype system spanning the full range of therapeutically desirable energies.
△ Less
Submitted 31 October, 2013;
originally announced October 2013.
-
A Method for Data-Driven Simulations of Evolving Solar Active Regions
Authors:
Mark C. M. Cheung,
Marc L. DeRosa
Abstract:
We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager…
▽ More
We present a method for performing data-driven simulations of solar active region formation and evolution. The approach is based on magnetofriction, which evolves the induction equation assuming the plasma velocity is proportional to the Lorentz force. The simulations of active region coronal field are driven by temporal sequences of photospheric magnetograms from the Helioseismic Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO). Under certain conditions, the data-driven simulations produce flux ropes that are ejected from the modeled active region due to loss of equilibrium. Following the ejection of flux ropes, we find an enhancement of the photospheric horizontal field near the polarity inversion line. We also present a method for the synthesis of mock coronal images based on a proxy emissivity calculated from the current density distribution in the model. This method yields mock coronal images that are somewhat reminiscent of images of active regions taken by instruments such as SDO's Atmospheric Imaging Assembly (AIA) at extreme ultraviolet wavelengths.
△ Less
Submitted 14 August, 2012;
originally announced August 2012.
-
Simulation of the Formation of a Solar Active Region
Authors:
M. C. M. Cheung,
M. Rempel,
A. M. Title,
M. Schüssler
Abstract:
We present a radiative magnetohydrodynamics simulation of the formation of an Active Region on the solar surface. The simulation models the rise of a buoyant magnetic flux bundle from a depth of 7.5 Mm in the convection zone up into the solar photosphere. The rise of the magnetic plasma in the convection zone is accompanied by predominantly horizontal expansion. Such an expansion leads to a scalin…
▽ More
We present a radiative magnetohydrodynamics simulation of the formation of an Active Region on the solar surface. The simulation models the rise of a buoyant magnetic flux bundle from a depth of 7.5 Mm in the convection zone up into the solar photosphere. The rise of the magnetic plasma in the convection zone is accompanied by predominantly horizontal expansion. Such an expansion leads to a scaling relation between the plasma density and the magnetic field strength such that $B\propto\varrho^{1/2}$. The emergence of magnetic flux into the photosphere appears as a complex magnetic pattern, which results from the interaction of the rising magnetic field with the turbulent convective flows. Small-scale magnetic elements at the surface first appear, followed by their gradual coalescence into larger magnetic concentrations, which eventually results in the formation of a pair of opposite polarity spots. Although the mean flow pattern in the vicinity of the developing spots is directed radially outward, correlations between the magnetic field and velocity field fluctuations allow the spots to accumulate flux. Such correlations result from the Lorentz-force driven, counter-streaming motion of opposite-polarity fragments. The formation of the simulated Active Region is accompanied by transient light bridges between umbrae and umbral dots. Together with recent sunspot modeling, this work highlights the common magnetoconvective origin of umbral dots, light bridges and penumbral filaments.
△ Less
Submitted 22 June, 2010; v1 submitted 21 June, 2010;
originally announced June 2010.
-
Multiscale Investigation of Chemical Interference in Proteins
Authors:
Antonios Samiotakis,
Dirar Homouz,
Margaret S. Cheung
Abstract:
We developed a multiscale approach (MultiSCAAL) that integrates the potential of mean force (PMF) obtained from all-atomistic molecular dynamics simulations with a knowledge-based energy function for coarse-grained molecular simulations in better exploring the energy landscape of a small protein under chemical interference such as chemical denaturation. An excessive amount of water molecules in al…
▽ More
We developed a multiscale approach (MultiSCAAL) that integrates the potential of mean force (PMF) obtained from all-atomistic molecular dynamics simulations with a knowledge-based energy function for coarse-grained molecular simulations in better exploring the energy landscape of a small protein under chemical interference such as chemical denaturation. An excessive amount of water molecules in all-atomistic molecular dynamics simulations often negatively impacts the sampling efficiency of some advanced sampling techniques such as the replica exchange method and it makes the investigation of chemical interferences on protein dynamics difficult. Thus, there is a need to develop an effective strategy that focuses on sampling structural changes in protein conformations rather than solvent molecule fluctuations. In this work, we address this issue by devising a multiscale simulation scheme (MultiSCAAL) that bridges the gap between all-atomistic molecular dynamics simulation and coarse-grained molecular simulation. The two key features of this scheme are the Boltzmann inversion and a protein atomistic reconstruction method we previously developed (SCAAL). Using MultiSCAAL, we were able to enhance the sampling efficiency of proteins solvated by explicit water molecules. Our method has been tested on the folding energy landscape of a small protein Trp-cage with explicit solvent under 8M urea using both the all-atomistic replica exchange molecular dynamics (AA-REMD) and MultiSCAAL. We compared computational analyses on ensemble conformations of Trp-cage with its available experimental NOE distances. The analysis demonstrated that conformations explored by MultiSCAAL better agree with the ones probed in the experiments because it can effectively capture the changes in side chain orientations that can flip out of the hydrophobic pocket in the presence of urea and water molecules.
△ Less
Submitted 9 April, 2010;
originally announced April 2010.
-
Manipulating Biopolymer Dynamics by Anisotropic Nanoconfinement
Authors:
Shao-Qing Zhang,
Margaret S. Cheung
Abstract:
How the geometry of nano-sized confinement affects dynamics of biomaterials is interesting yet poorly understood. An elucidation of structural details upon nano-sized confinement may benefit manufacturing pharmaceuticals in biomaterial sciences and medicine. The behavior of biopolymers in nano-sized confinement is investigated using coarse-grained models and molecular simulations. Particularly,…
▽ More
How the geometry of nano-sized confinement affects dynamics of biomaterials is interesting yet poorly understood. An elucidation of structural details upon nano-sized confinement may benefit manufacturing pharmaceuticals in biomaterial sciences and medicine. The behavior of biopolymers in nano-sized confinement is investigated using coarse-grained models and molecular simulations. Particularly, we address the effects of shapes of a confinement on protein folding dynamics by measuring folding rates and dissecting structural properties of the transition states in nano-sized spheres and ellipsoids. We find that when the form of a confinement resembles the geometrical properties of the transition states, the rates of folding kinetics are most enhanced. This knowledge of shape selectivity in identifying optimal conditions for reactions will have a broad impact in nanotechnology and pharmaceutical sciences.
△ Less
Submitted 9 October, 2007; v1 submitted 8 October, 2007;
originally announced October 2007.
