Source author record

Shravan Hanasoge

Shravan Hanasoge appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 published item(s)

preprint2026arXiv

Inferring Asteroseismic Parameters from Short Observations Using Deep Learning: Application to TESS and K2 Red Giants

Asteroseismology is the study of resonant oscillations of stars to infer their internal structure and dynamics. It is also a powerful tool for precisely determining stellar parameters such as mass, radius, surface gravity, and age. The ongoing TESS mission, with its nearly complete sky coverage, presents a unique opportunity to uniformly probe stellar populations across the Milky Way. TESS is estimated to have observed more than 300,000 oscillating red giants, most of which have one to two months of observations. Given the scale of this dataset, we need a fast, efficient, and robust way to analyse the data. In this work, our objective is to develop a machine learning (ML) based method to infer asteroseismic parameters from short-duration observations. Specifically, we focus on two global seismic parameters, the large frequency separation ($Δν$) and the frequency at maximum power ($ν_{\mathrm{max}}$), from one-month-long TESS observations of red giants. Meanwhile, for K2 data, our focus extends to inferring the period spacings of dipolar gravity modes ($ΔΠ_{1}$), in addition to $Δν$ and $ν_{\mathrm{max}}$. Our findings demonstrate that our machine learning algorithm can accurately infer $Δν$ and $ν_{\mathrm{max}}$ for approximately 50% of samples created by taking one-month Kepler and K2 observations. For TESS one sector data however, we recover reliable $Δν$ for only about 23% of the stars. Additionally, we get reliable $ΔΠ_{1}$ inferences for about 200 young red-giants from K2. For these $ΔΠ_{1}$ inferences, we see a good match with the well known $Δν-ΔΠ_{1}$ degenerate sequence observed in Kepler red-giants.

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

Imaging the Sun's near-surface flows using mode-coupling analysis

The technique of normal-mode coupling is a powerful tool with which to seismically image non-axisymmetric phenomena in the Sun. Here we apply mode coupling in the Cartesian approximation to probe steady, near-surface flows in the Sun. Using Doppler cubes obtained from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory, we perform inversions on mode-coupling measurements to show that the resulting divergence and radial vorticity maps at supergranular length scales ($\sim$30 Mm) near the surface compare extremely well with those obtained using the Local Correlation Tracking method. We find that the Pearson correlation coefficient is $\geq$ 0.9 for divergence flows, while $\geq$ 0.8 is obtained for the radial vorticity.

preprint2022arXiv

Measuring frequency and period separations in red-giant stars using machine learning

Asteroseismology is used to infer the interior physics of stars. The \textit{Kepler} and TESS space missions have provided a vast data set of red-giant light curves, which may be used for asteroseismic analysis. These data sets are expected to significantly grow with future missions such as \textit{PLATO}, and efficient methods are therefore required to analyze these data rapidly. Here, we describe a machine learning algorithm that identifies red giants from the raw oscillation spectra and captures \textit{p} and \textit{mixed} mode parameters from the red-giant power spectra. We report algorithmic inferences for large frequency separation ($Δν$), frequency at maximum amplitude ($ν_{max}$), and period separation ($ΔΠ$) for an ensemble of stars. In addition, we have discovered $\sim$25 new probable red giants among 151,000 \textit{Kepler} long-cadence stellar-oscillation spectra analyzed by the method, among which four are binary candidates which appear to possess red-giant counterparts. To validate the results of this method, we selected $\sim$ 3,000 \textit{Kepler} stars, at various evolutionary stages ranging from subgiants to red clumps, and compare inferences of $Δν$, $ΔΠ$, and $ν_{max}$ with estimates obtained using other techniques. The power of the machine-learning algorithm lies in its speed: it is able to accurately extract seismic parameters from 1,000 spectra in $\sim$5 seconds on a modern computer (single core of the Intel Xeon Platinum 8280 CPU).

preprint2020arXiv

Properties of solar Rossby waves from normal mode coupling and characterizing its systematics

Rossby waves play an important role in mediating the angular momentum of rotating spherical fluids, creating weather on Earth and tuning exoplanet orbits in distant stellar systems (Ogilvie 2014). Their recent discovery in the solar convection zone provides an exciting opportunity to appreciate the detailed astrophysics of Rossby waves. Large-scale Rossby waves create subtle drifts in acoustic oscillations in the convection zone, which we measure using helioseismology to image properties of Rossby waves in the interior. We analyze 20 years of space-based observations, from 1999 to 2018, to measure Rossby-mode frequencies, line-widths, and amplitudes. Spatial leakage affects the measurements of normal model coupling and complicates the analysis of separating out specific harmonic degree and the azimuthal number of features on the Sun. Here we demonstrate a novel approach to overcome this difficulty and test it by performing synthetic tests. We find that the root-mean-square velocity of the modes is of the order of 0.5 m/s at the surface.

preprint2020arXiv

Solar wind prediction using deep learning

Emanating from the base of the Sun's corona, the solar wind fills the interplanetary medium with a magnetized stream of charged particles whose interaction with the Earth's magnetosphere has space-weather consequences such as geomagnetic storms. Accurately predicting the solar wind through measurements of the spatio-temporally evolving conditions in the solar atmosphere is important but remains an unsolved problem in heliophysics and space-weather research. In this work, we use deep learning for prediction of solar wind (SW) properties. We use Extreme Ultraviolet images of the solar corona from space based observations to predict the SW speed from the NASA OMNIWEB dataset, measured at Lagragian point 1. We evaluate our model against autoregressive and naive models, and find that our model outperforms the benchmark models, obtaining a best-fit correlation of 0.55 $\pm$ 0.03 with the observed data. Upon visualization and investigation of how the model uses data to make predictions, we find higher activation at the coronal holes for fast wind prediction ($\approx$ 3 to 4 days prior to prediction), and at the active regions for slow wind prediction. These trends bear an uncanny similarity to the influence of regions potentially being the sources of fast and slow wind, as reported in literature. This suggests that our model was able to learn some of the salient associations between coronal and solar wind structure without built-in physics knowledge. Such an approach may help us discover hitherto unknown relationships in heliophysics datasets.

preprint2020arXiv

Validating inversions for toroidal flows using normal-mode coupling

Normal-mode coupling is a helioseismic technique that uses measurements of mode eigenfunctions to infer the interior structure of the Sun. This technique has led to insights into the evolution and structure of toroidal flows in the solar interior. Here, we validate an inversion algorithm for normal-mode coupling by generating synthetic seismic measurements associated with input flows and comparing the input and inverted velocities. We study four different cases of input toroidal flows and compute synthetics that take into account the partial visibility of the Sun. We invert the synthetics using Subtractive Optimally Localized Averages (SOLA) and also try to mitigate the systematics of mode leakage. We demonstrate that, ultimately, inversions are only as good as the model we assume for the correlation between flow velocities.

preprint2016arXiv

Asteroseismic determination of fundamental parameters of sun-like stars using multi-layered neural networks

The advent of space-based observatories such as CoRoT and Kepler has enabled the testing of our understanding of stellar evolution on thousands of stars. Evolutionary models typically require five input parameters, the mass, initial Helium abundance, initial metallicity, mixing- length (assumed to be constant over time), and the age to which the star must be evolved. Some of these parameters are also very useful in characterizing the associated planets and in studying galactic archaeology. How to obtain these parameters from observations rapidly and accurately, specifically in the context of surveys of thousands of stars, is an outstanding ques- tion, one that has eluded straightforward resolution. For a given star, we typically measure the effective temperature and surface metallicity spectroscopically and low-degree oscillation frequencies through space observatories. Here we demonstrate that statistical learning, using artificial neural networks, is successful in determining the evolutionary parameters based on spectroscopic and seismic measurements. Our trained networks show robustness over a broad range of parameter space, and critically, are entirely computationally inexpensive and fully automated. We analyze the observations of a few stars using this method and the results com- pare well to inferences obtained using other techniques. This method is both computationally cheap and inferentially accurate, paving the way for analyzing the vast quantities of stellar observations from past, current, and future missions.

preprint2016arXiv

Spatio-spectral concentration of convolutions

Differential equations may possess coefficients that vary on a spectrum of scales. Because coefficients are typically multiplicative in real space, they turn into convolution operators in spectral space, mixing all wavenumbers. However, in many applications, only the largest scales of the solution are of interest and so the question turns to whether it is possible to build effective coarse-scale models of the coefficients in such a manner that the large scales of the solution are left intact. Here we apply the method of numerical homogenization to deterministic linear equations to generate sub-grid-scale models of coefficients at desired frequency cutoffs. We use the Fourier basis to project, filter and compute correctors for the coefficients. The method is tested in 1D and 2D scenarios and found to reproduce the coarse scales of the solution to varying degrees of accuracy depending on the cutoff. We relate this method to mode-elimination Renormalization Group (RG) and discuss the connection between accuracy and the cutoff wavenumber. The tradeoff is governed by a form of the uncertainty principle for convolutions, which states that as the convolution operator is squeezed in the spectral domain, it broadens in real space. As a consequence, basis sparsity is a high virtue and the choice of the basis can be critical.

preprint2015arXiv

Seismic Sounding of Convection in the Sun

Our Sun, primarily composed of ionized hydrogen and helium, has a surface temperature of 5777~K and a radius $R_\odot \approx 696,000$ km. In the outer $R_\odot/3$, energy transport is accomplished primarily by convection. Using typical convective velocities $u\sim100\,\rm{m\,s^{-1}}$ and kinematic viscosities of order $10^{-4}$ m$^{2}$s$^{-1}$, we obtain a Reynolds number $Re \sim 10^{14}$. Convection is thus turbulent, causing a vast range of scales to be excited. The Prandtl number, $Pr$, of the convecting fluid is very low, of order $10^{-7}$\,--\,$10^{-4}$, so that the Rayleigh number ($\sim Re^2 Pr$) is on the order of $10^{21}\,-\,10^{24}$. Solar convection thus lies in extraordinary regime of dynamical parameters, highly untypical of fluid flows on Earth. Convective processes in the Sun drive global fluid circulations and magnetic fields, which in turn affect its visible outer layers ("solar activity") and, more broadly, the heliosphere ("space weather"). The precise determination of the depth of solar convection zone, departures from adiabaticity of the temperature gradient, and the internal rotation rate as a function of latitude and depth are among the seminal contributions of helioseismology towards understanding convection in the Sun. Contemporary helioseismology, which is focused on inferring the properties of three-dimensional convective features, suggests that transport velocities are substantially smaller than theoretical predictions. Furthermore, helioseismology provides important constraints on the anisotropic Reynolds stresses that control the global dynamics of the solar convection zone. This review discusses the state of our understanding of convection in the Sun, with a focus on helioseismic diagnostics. We present our considerations with the interests of fluid dynamicists in mind.

preprint2014arXiv

Full Waveform Inversion for Time-Distance Helioseismology

Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions, and likely for supergranular flows as well. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the sub-surface structure of large-amplitude perturbations. Defining the penalty functional as the $L_2$ norm of wave travel-time deviations, we compute the the total misfit gradient of this functional with respect to the relevant model parameters %(only sound speed in this case) at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or quasi-Newton limited-memory BFGS. Performing nonlinear iterations requires privileging pixels (such as those in the near-field of the scatterer), a practice not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a 2-D inverse problem with only a sound-speed perturbation.

preprint2013arXiv

Interpreting Cross-correlations of One-bit Filtered Seismic Noise

Seismic noise, generated by oceanic microseisms and other sources, illuminates the crust in a manner different from tectonic sources, and therefore provides independent information. The primary measurable is the two-point cross-correlation, evaluated using traces recorded at a pair of seismometers over a finite-time interval. However, raw seismic traces contain intermittent large-amplitude perturbations arising from tectonic activity and instrumental errors, which may corrupt the estimated cross-correlations of microseismic fluctuations. In order to diminish the impact of these perturbations, the recorded traces are filtered using the nonlinear one-bit digitizer, which replaces the measurement by its sign. Previous theory shows that for stationary Gaussian-distributed seismic noise fluctuations one-bit and raw correlation functions are related by a simple invertible transformation. Here we extend this to show that the simple correspondence between these two correlation techniques remains valid for {\it non-stationary} Gaussian and a very broad range of {\it non-Gaussian} processes as well. For a limited range of stationary and non-stationary Gaussian fluctuations, we find that one-bit filtering performs at least as well as spectral whitening. We therefore recommend using one-bit filtering when processing terrestrial seismic noise, with the substantial benefit that the measurements are fully compatible with current theoretical interpretation (e.g., adjoint theory). Given that seismic records are non-stationary and comprise small-amplitude fluctuations and intermittent, large-amplitude tectonic/other perturbations, we outline an algorithm to accurately retrieve the correlation function of the small-amplitude signals.

preprint2013arXiv

Propagation of seismic waves through a spatio-temporally fluctuating medium: Homogenization

Measurements of seismic wave travel times at the photosphere of the Sun have enabled inferences of its interior structure and dynamics. In interpreting these measurements, the simplifying assumption that waves propagate through a temporally stationary medium is almost universally invoked. However, the Sun is in a constant state of evolution, on a broad range of spatio-temporal scales. At the zero wavelength limit, i.e., when the wavelength is much shorter than the scale over which the medium varies, the WKBJ (ray) approximation may be applied. Here, we address the other asymptotic end of the spectrum, the infinite wavelength limit, using the technique of homogenization. We apply homogenization to scenarios where waves are propagating through rapidly varying media (spatially and temporally), and derive effective models for the media. One consequence is that a scalar sound speed becomes a tensorial wavespeed in the effective model and anisotropies can be induced depending on the nature of the perturbation. The second term in this asymptotic two-scale expansion, the so-called corrector, contains contributions due to higher-order scattering, leading to the decoherence of the wavefield. This decoherence may be causally linked to the observed wave attenuation in the Sun. Although the examples we consider here consist of periodic arrays of perturbations to the background, homogenization may be extended to ergodic and stationary random media. This method may have broad implications for the manner in which we interpret seismic measurements in the Sun and for modeling the effects of granulation on the scattering of waves and distortion of normal-mode eigenfunctions.

preprint2013arXiv

Seismic constraints on rotation of Sun-like star and mass of exoplanet

Rotation is thought to drive cyclic magnetic activity in the Sun and Sun-like stars. Stellar dynamos, however, are poorly understood owing to the scarcity of observations of rotation and magnetic fields in stars. Here, inferences are drawn on the internal rotation of a distant Sun-like star by studying its global modes of oscillation. We report asteroseismic constraints imposed on the rotation rate and the inclination of the spin axis of the Sun-like star HD 52265, a principal target observed by the CoRoT satellite that is known to host a planetary companion. These seismic inferences are remarkably consistent with an independent spectroscopic observation (rotational line broadening) and with the observed rotation period of star spots. Furthermore, asteroseismology constrains the mass of exoplanet HD 52265b. Under the standard assumption that the stellar spin axis and the axis of the planetary orbit coincide, the minimum spectroscopic mass of the planet can be converted into a true mass of 1.85 (+0.52,-0.42) M_Jupiter, which implies that it is a planet, not a brown dwarf.

preprint2012arXiv

Anomalously Weak Solar Convection

Convection in the solar interior is thought to comprise structures on a spectrum of scales. This conclusion emerges from phenomenological studies and numerical simulations, though neither covers the proper range of dynamical parameters of solar convection. Here, we analyze observations of the wavefield in the solar photosphere using techniques of time-distance helioseismology to image flows in the solar interior. We downsample and synthesize 900 billion wavefield observations to produce 3 billion cross-correlations, which we average and fit, measuring 5 million wave travel times. Using these travel times, we deduce the underlying flow systems and study their statistics to bound convective velocity magnitudes in the solar interior, as a function of depth and spherical-harmonic degree $\ell$. Within the wavenumber band $\ell<60$, Convective velocities are 20-100 times weaker than current theoretical estimates. This suggests the prevalence of a different paradigm of turbulence from that predicted by existing models, prompting the question: what mechanism transports the heat flux of a solar luminosity outwards? Advection is dominated by Coriolis forces for wavenumbers $\ell<60$, with Rossby numbers smaller than $\sim10^{-2}$ at $r/R_\odot=0.96$, suggesting that the Sun may be a much faster rotator than previously thought, and that large-scale convection may be quasi-geostrophic. The fact that iso-rotation contours in the Sun are not co-aligned with the axis of rotation suggests the presence of a latitudinal entropy gradient.

preprint2012arXiv

Detectable seismic consequences of the interaction of a primordial black hole with Earth

Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate.

preprint2012arXiv

Seismic probes of solar interior magnetic structure

Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfvén waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun.

preprint2011arXiv

The Adjoint Method Applied to Time-Distance Helioseismology

For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fréchet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed.

preprint2011arXiv

Transient solar oscillations driven by primordial black holes

Stars are transparent to the passage of primordial black holes (PBHs) and serve as seismic detectors for such objects. The gravitational field of a PBH squeezes a star and causes it to ring acoustically. We calculate the seismic signature of a PBH passing through the Sun. The background for this signal is the observed spectrum of solar oscillations excited by supersonic turbulence. We predict that PBHs more massive than 10^21 g (comparable in mass to an asteroid) are detectable by existing solar observatories. The oscillations excited by PBHs peak at large scales and high frequencies, making them potentially detectable in other stars. The discovery of PBHs would have profound implications for cosmology and high-energy physics.