Researcher profile

Xudong Sun

Xudong Sun contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Convergence of Consensus-Based Particle Methods for Nonconvex Bi-Level Optimization

In this paper, we study a consensus-based optimization method for nonconvex bi-level optimization, where the objective is to minimize an upper-level function over the set of global minimizers of a lower-level problem. The proposed approach is derivative-free, and constructs its consensus point via smooth quantile selection combined with a Gibbs-type Laplace approximation. We establish convergence guarantees for both the associated \textit{mean-field} dynamics and its \textit{finite-particle} approximation. In particular, under suitable assumptions on smooth quantile localization, error bounds, and stability, we show that the mean-field law reaches any arbitrary prescribed Wasserstein neighborhood of the target bi-level solution with an explicit exponential rate up to the hitting time. Numerical experiments on a two-dimensional constrained problem and neural network training further support the theoretical results.

preprint2022arXiv

A New Magnetic Parameter of Active Regions Distinguishing Large Eruptive and Confined Solar Flares

With the aim of investigating how the magnetic field in solar active regions (ARs) controls flare activity, i.e., whether a confined or eruptive flare occurs, we analyze 106 flares of Geostationary Operational Environmental Satellite (GOES) class $\geq$M1.0 during 2010$-$2019. We calculate mean characteristic twist parameters $α$$_{FPIL}$ within the "flaring polarity inversion line" region and $α$$_\mathrm{HFED}$ within the area of high photospheric magnetic free energy density, which both provide measures of the nonpotentiality of AR core region. Magnetic twist is thought to be related to the driving force of electric current-driven instabilities, such as the helical kink instability. We also calculate total unsigned magnetic flux ($Φ$$_\mathrm{AR}$) of ARs producing the flare, which describes the strength of the background field confinement. By considering both the constraining effect of background magnetic fields and the magnetic non-potentiality of ARs, we propose a new parameter $α$/$Φ$$_\mathrm{AR}$ to measure the probability for a large flare to be associated with a coronal mass ejection (CME). We find that in about 90\% of eruptive flares, $α$$_\mathrm{FPIL}$/$Φ$$_\mathrm{AR}$ and $α$$_\mathrm{HFED}$/$Φ$$_\mathrm{AR}$ are beyond critical values (2.2$\times$$10^{-24}$ and 3.2$\times$$10^{-24}$ Mm$^{-1}$ Mx$^{-1}$), whereas they are less than critical values in $\sim$ 80\% of confined flares. This indicates that the new parameter $α$/$Φ$$_\mathrm{AR}$ is well able to distinguish eruptive flares from confined flares. Our investigation suggests that the relative measure of magnetic nonpotentiality within the AR core over the restriction of the background field largely controls the capability of ARs to produce eruptive flares.

preprint2022arXiv

Exploring the Solar Poles: The Last Great Frontier of the Sun

Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere.

preprint2022arXiv

On the Coordinate System of Space-Weather HMI Active Region Patches (SHARPs): A Technical Note

We describe the coordinate systems of two streams of HMI active region vector data. A distinction is made between (a) the 2D grid on which the field vector is measured (or sampled), and (b) the 3D coordinate established at each grid point, in which the field vector is presented. The HMI data reduction can involve coordinate changes on both, with those performed on the former termed "remapping", the latter "vector transformation". Relevant pipeline procedures are described. Useful examples are given for data analysis.

preprint2022arXiv

On the nature of photospheric horizontal magnetic field increase in major solar flares

Rapid increase of horizontal magnetic field ($B_h$) around the flaring polarity inversion line is the most prominent photospheric field change during flares. It is considered to be caused by the contraction of flare loops, the details behind which is still not fully understood. Here we investigate the $B_h$-increase in 35 major flares using HMI high-cadence vector magnetograms. We find that $B_h$-increase is always accompanied by the increase of field inclination. It usually initiates near the flare ribbons, showing step-like change in between the ribbons. In particular, its evolution in early flare phase shows close spatio-temporal correlation to flare ribbons. We further find that $B_h$-increase tends to have similar intensity in confined and eruptive flares, but larger spatial-extent in eruptive flares in a statistical sense. Its intensity and timescale have inverse and positive correlations to the initial ribbon separations, respectively. The results altogether are well consistent with a recent proposed scenario which suggests that the reconnection-driven contraction of flare loops enhances photospheric $B_h$ according to the ideal induction equation, providing statistical evidence to the reconnection-driven origin for $B_h$-increase for the first time.

preprint2022arXiv

The CGEM Lorentz Force Data from HMI Vector Magnetograms

We describe a new data product from the CGEM (Coronal Global Evolutionary Model) collaboration that estimates the Lorentz force in active regions (ARs) based on HMI vector magnetogram patches. Following Fisher et al. (2012), we compute three components of the integrated Lorentz force over the outer solar atmosphere every 12 minutes throughout an AR's disk passage. These estimates, differenced during solar eruptive events, can provide valuable diagnostics on dynamic processes. We describe the pipeline modules, provide data retrieval examples, and document some systematic uncertainties that users should be aware of. Finally we document the formal uncertainty propagation procedures.

preprint2021arXiv

Benchmarking time series classification -- Functional data vs machine learning approaches

Time series classification problems have drawn increasing attention in the machine learning and statistical community. Closely related is the field of functional data analysis (FDA): it refers to the range of problems that deal with the analysis of data that is continuously indexed over some domain. While often employing different methods, both fields strive to answer similar questions, a common example being classification or regression problems with functional covariates. We study methods from functional data analysis, such as functional generalized additive models, as well as functionality to concatenate (functional-) feature extraction or basis representations with traditional machine learning algorithms like support vector machines or classification trees. In order to assess the methods and implementations, we run a benchmark on a wide variety of representative (time series) data sets, with in-depth analysis of empirical results, and strive to provide a reference ranking for which method(s) to use for non-expert practitioners. Additionally, we provide a software framework in R for functional data analysis for supervised learning, including machine learning and more linear approaches from statistics. This allows convenient access, and in connection with the machine-learning toolbox mlr, those methods can now also be tuned and benchmarked.

preprint2021arXiv

Toward Improved Understanding of Magnetic Fields Participating in Solar Flares: Statistical Analysis of Magnetic Field within Flare Ribbons

Violent solar flares and coronal mass ejections (CMEs) are magnetic phenomena. However, how magnetic fields reconnecting in the flare differ from non-flaring magnetic fields remains unclear owing to the lack of studies of the flare magnetic properties. Here we present a first statistical study of flaring (highlighted by flare-ribbons) vector magnetic fields in the photosphere. Our systematic approach allows us to describe key physical properties of solar flare magnetism, including distributions of magnetic flux, magnetic shear, vertical current and net current over flaring versus non-flaring parts of the active region, and compare these with flare/CME properties. Our analysis suggests that while flares are guided by the physical properties that scale with AR size, like the total amount of magnetic flux that participates in the reconnection process and the total current (extensive properties), CMEs are guided by mean properties, like the fraction of the AR magnetic flux that participates (intensive property), with little dependence on the amount of shear at polarity inversion line (PIL) or the net current. We find that the non-neutralized current is proportional to the amount of shear at PIL, providing direct evidence that net vertical currents are formed as a result of any mechanism that could generate magnetic shear along PIL. We also find that eruptive events tend to have smaller PIL fluxes and larger magnetic shears than confined events. Our analysis provides a reference for more realistic solar and stellar flare models. The database is available online and can be used for future quantitative studies of flare magnetism.

preprint2020arXiv

Electric Current Neutralization in Solar Active Regions and Its Relation to Eruptive Activity

It is well established that magnetic free energy associated with electric currents powers solar flares and coronal mass ejections (CMEs) from solar active regions (ARs). However, the conditions that determine whether an AR will produce an eruption are not well understood. Previous work suggests that the degree to which the driving electric currents, or the sum of all currents within a single magnetic polarity, are neutralized may serve as a good proxy for assessing the ability of ARs to produce eruptions. Here, we investigate the relationship between current neutralization and flare/CME production using a sample of 15 flare-active and 15 flare-quiet ARs. All flare-quiet and 4 flare-active ARs are also CME-quiet. We additionally test the relation of current neutralization to the degree of shear along polarity inversion lines (PILs) in an AR. We find that flare-productive ARs are more likely to exhibit non-neutralized currents, specifically those that also produce a CME. We find that flare/CME-active ARs also exhibit higher degrees of PIL shear than flare/CME-quiet ARs. We additionally observe that currents become more neutralized during magnetic flux emergence in flare-quiet ARs. Our investigation suggests that current neutralization in ARs is indicative of their eruptive potential.

preprint2020arXiv

Maximum Entropy-Regularized Multi-Goal Reinforcement Learning

In Multi-Goal Reinforcement Learning, an agent learns to achieve multiple goals with a goal-conditioned policy. During learning, the agent first collects the trajectories into a replay buffer, and later these trajectories are selected randomly for replay. However, the achieved goals in the replay buffer are often biased towards the behavior policies. From a Bayesian perspective, when there is no prior knowledge about the target goal distribution, the agent should learn uniformly from diverse achieved goals. Therefore, we first propose a novel multi-goal RL objective based on weighted entropy. This objective encourages the agent to maximize the expected return, as well as to achieve more diverse goals. Secondly, we developed a maximum entropy-based prioritization framework to optimize the proposed objective. For evaluation of this framework, we combine it with Deep Deterministic Policy Gradient, both with or without Hindsight Experience Replay. On a set of multi-goal robotic tasks of OpenAI Gym, we compare our method with other baselines and show promising improvements in both performance and sample-efficiency.

preprint2020arXiv

Non-Neutralized Electric Current of Active Regions Explained as a Projection Effect

Active regions (ARs) often possess an observed net electric current in a single magnetic polarity. We show that such "non-neutralized" currents can arise from a geometric projection effect when a twisted flux tube obliquely intersects the photosphere. To this end, we emulate surface maps of an emerging AR by sampling horizontal slices of a semi-torus flux tube at various heights. Although the tube has no net toroidal current, its poloidal current, when projected along the vertical direction, amounts to a significant non-neutralized component on the surface. If the tube emerges only partially as in realistic settings, the non-neutralized current will 1) develop as double ribbons near the sheared polarity inversion line, (2) positively correlate with the twist, and 3) reach its maximum before the magnetic flux. The projection effect may be important to the photospheric current distribution, in particular during the early stages of flux emergence.

preprint2020arXiv

The PDFI_SS Electric Field Inversion Software

We describe the PDFI_SS software library, which is designed to find the electric field at the Sun's photosphere from a sequence of vector magnetogram and Doppler velocity measurements, and estimates of horizontal velocities obtained from local correlation tracking using the recently upgraded FLCT code. The library, a collection of Fortran subroutines, uses the "PDFI" technique described by Kazachenko et al. (2014), but modified for use in spherical, Plate-Carrée geometry on a staggered grid. The domain over which solutions are found is a subset of the global spherical surface, defined by user-specified limits of colatitude and longitude. Our staggered-grid approach, based on that of Yee (1966), is more conservative and self-consistent compared to the centered, Cartesian grid used by Kazachenko et al. (2014). The library can be used to compute an end-to-end solution for electric fields from data taken by the HMI instrument aboard NASA's SDO Mission. This capability has been incorporated into the HMI pipeline processing system operating at SDO's JSOC. The library is written in a general and modular way so that the calculations can be customized to modify or delete electric field contributions, or used with other data sets. Other applications include "nudging" numerical models of the solar atmosphere to facilitate assimilative simulations. The library includes an ability to compute "global" (whole-Sun) electric field solutions. The library also includes an ability to compute Potential Magnetic Field solutions in spherical coordinates. This distribution includes a number of test programs which allow the user to test the software.