BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Braxton Osting

Braxton Osting contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.OCMachine Learningmath.SPmath.APmath.DSmath.NAmath.PRNumerical AnalysisComputationComputational Geometrycond-mat.mtrl-sciDiscrete Mathematicsmath.DGmath.STphysics.flu-dynphysics.med-phStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
17topics
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

11 published item(s)

preprint2025arXiv

Cuffless, calibration-free hemodynamic monitoring with physics-informed machine learning models

Wearable technologies have the potential to transform ambulatory and at-home hemodynamic monitoring by providing continuous assessments of cardiovascular health metrics and guiding clinical management. However, existing cuffless wearable devices for blood pressure (BP) monitoring often rely on methods lacking theoretical foundations, such as pulse wave analysis or pulse arrival time, making them vulnerable to physiological and experimental confounders that undermine their accuracy and clinical utility. Here, we developed a smartwatch device with real-time electrical bioimpedance (BioZ) sensing for cuffless hemodynamic monitoring. We elucidate the biophysical relationship between BioZ and BP via a multiscale analytical and computational modeling framework, and identify physiological, anatomical, and experimental parameters that influence the pulsatile BioZ signal at the wrist. A signal-tagged physics-informed neural network incorporating fluid dynamics principles enables calibration-free estimation of BP and radial and axial blood velocity. We successfully tested our approach with healthy individuals at rest and after physical activity including physical and autonomic challenges, and with patients with hypertension and cardiovascular disease in outpatient and intensive care settings. Our findings demonstrate the feasibility of BioZ technology for cuffless BP and blood velocity monitoring, addressing critical limitations of existing cuffless technologies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A dynamical systems based framework for dimension reduction

We propose a novel framework for learning a low-dimensional representation of data based on nonlinear dynamical systems, which we call dynamical dimension reduction (DDR). In the DDR model, each point is evolved via a nonlinear flow towards a lower-dimensional subspace; the projection onto the subspace gives the low-dimensional embedding. Training the model involves identifying the nonlinear flow and the subspace. Following the equation discovery method, we represent the vector field that defines the flow using a linear combination of dictionary elements, where each element is a pre-specified linear/nonlinear candidate function. A regularization term for the average total kinetic energy is also introduced and motivated by optimal transport theory. We prove that the resulting optimization problem is well-posed and establish several properties of the DDR method. We also show how the DDR method can be trained using a gradient-based optimization method, where the gradients are computed using the adjoint method from optimal control theory. The DDR method is implemented and compared on synthetic and example datasets to other dimension reductions methods, including PCA, t-SNE, and Umap.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Isoperimetric Sloshing Problem in a Shallow Container with Surface Tension

In 1965, B. A. Troesch solved the isoperimetric sloshing problem of determining the container shape that maximizes the fundamental sloshing frequency among two classes of shallow containers: symmetric canals with a given free surface width and cross-sectional area, and radially symmetric containers with a given rim radius and volume [doi:10.1002/cpa.3160180124]. Here, we extend these results in two ways: (i) we consider surface tension effects on the fluid free surface, assuming a flat equilibrium free surface together with a pinned contact line, and (ii) we consider sinusoidal waves traveling along the canal with wavenumber $α\ge 0$ and spatial period $2π/α$; two-dimensional sloshing corresponds to the case $α= 0$. Generalizing our recent variational characterization of fluid sloshing with surface tension to the case of a pinned contact line, we derive the pinned-edge linear shallow sloshing problem, which is an eigenvalue problem for a generalized Sturm-Liouville system. In the case without surface tension, we show that the optimal shallow canal is a rectangular canal for any $α> 0$. In the presence of surface tension, we solve for the maximizing cross-section explicitly for shallow canals with any given $α\ge 0$ and shallow radially symmetric containers with $m$ azimuthal nodal lines, $m = 0, 1$. Our results reveal that the squared maximal sloshing frequency increases considerably as surface tension increases. Interestingly, both the optimal shallow canal for $α= 0$ and the optimal shallow radially symmetric container are not convex. As a consequence of our explicit solutions, we establish convergence of the maximizing cross-sections, as surface tension vanishes, to the maximizing cross-sections without surface tension.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extremal graph realizations and graph Laplacian eigenvalues

For a regular polyhedron (or polygon) centered at the origin, the coordinates of the vertices are eigenvectors of the graph Laplacian for the skeleton of that polyhedron (or polygon) associated with the first (non-trivial) eigenvalue. In this paper, we generalize this relationship. For a given graph, we study the eigenvalue optimization problem of maximizing the first (non-trivial) eigenvalue of the graph Laplacian over non-negative edge weights. We show that the spectral realization of the graph using the eigenvectors corresponding to the solution of this problem, under certain assumptions, is a centered, unit-distance graph realization that has maximal total variance. This result gives a new method for generating unit-distance graph realizations and is based on convex duality. A drawback of this method is that the dimension of the realization is given by the multiplicity of the extremal eigenvalue, which is typically unknown prior to solving the eigenvalue optimization problem. Our results are illustrated with a number of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flat tori with large Laplacian eigenvalues in dimensions up to eight

We consider the optimization problem of maximizing the $k$-th Laplacian eigenvalue, $λ_{k}$, over flat $d$-dimensional tori of fixed volume. For $k=1$, this problem is equivalent to the densest lattice sphere packing problem. For larger $k$, this is equivalent to the NP-hard problem of finding the $d$-dimensional (dual) lattice with longest $k$-th shortest lattice vector. As a result of extensive computations, for $d \leq 8$, we obtain a sequence of flat tori, $T_{k,d}$, each of volume one, such that the $k$-th Laplacian eigenvalue of $T_{k,d}$ is very large; for each (finite) $k$ the $k$-th eigenvalue exceeds the value in (the $k\to \infty$ asymptotic) Weyl's law by a factor between 1.54 and 2.01, depending on the dimension. Stationarity conditions are derived and numerically verified for $T_{k,d}$ and we describe the degeneration of the tori as $k \to \infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Probabilistic methods for approximate archetypal analysis

Archetypal analysis is an unsupervised learning method for exploratory data analysis. One major challenge that limits the applicability of archetypal analysis in practice is the inherent computational complexity of the existing algorithms. In this paper, we provide a novel approximation approach to partially address this issue. Utilizing probabilistic ideas from high-dimensional geometry, we introduce two preprocessing techniques to reduce the dimension and representation cardinality of the data, respectively. We prove that provided the data is approximately embedded in a low-dimensional linear subspace and the convex hull of the corresponding representations is well approximated by a polytope with a few vertices, our method can effectively reduce the scaling of archetypal analysis. Moreover, the solution of the reduced problem is near-optimal in terms of prediction errors. Our approach can be combined with other acceleration techniques to further mitigate the intrinsic complexity of archetypal analysis. We demonstrate the usefulness of our results by applying our method to summarize several moderately large-scale datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A continuum limit for the PageRank algorithm

Semi-supervised and unsupervised machine learning methods often rely on graphs to model data, prompting research on how theoretical properties of operators on graphs are leveraged in learning problems. While most of the existing literature focuses on undirected graphs, directed graphs are very important in practice, giving models for physical, biological, or transportation networks, among many other applications. In this paper, we propose a new framework for rigorously studying continuum limits of learning algorithms on directed graphs. We use the new framework to study the PageRank algorithm, and show how it can be interpreted as a numerical scheme on a directed graph involving a type of normalized graph Laplacian. We show that the corresponding continuum limit problem, which is taken as the number of webpages grows to infinity, is a second-order, possibly degenerate, elliptic equation that contains reaction, diffusion, and advection terms. We prove that the numerical scheme is consistent and stable and compute explicit rates of convergence of the discrete solution to the solution of the continuum limit PDE. We give applications to proving stability and asymptotic regularity of the PageRank vector. Finally, we illustrate our results with numerical experiments and explore an application to data depth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A non-autonomous equation discovery method for time signal classification

Certain neural network architectures, in the infinite-layer limit, lead to systems of nonlinear differential equations. Motivated by this idea, we develop a framework for analyzing time signals based on non-autonomous dynamical equations. We view the time signal as a forcing function for a dynamical system that governs a time-evolving hidden variable. As in equation discovery, the dynamical system is represented using a dictionary of functions and the coefficients are learned from data. This framework is applied to the time signal classification problem. We show how gradients can be efficiently computed using the adjoint method, and we apply methods from dynamical systems to establish stability of the classifier. Through a variety of experiments, on both synthetic and real datasets, we show that the proposed method uses orders of magnitude fewer parameters than competing methods, while achieving comparable accuracy. We created the synthetic datasets using dynamical systems of increasing complexity; though the ground truth vector fields are often polynomials, we find consistently that a Fourier dictionary yields the best results. We also demonstrate how the proposed method yields graphical interpretability in the form of phase portraits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computation of free boundary minimal surfaces via extremal Steklov eigenvalue problems

Recently Fraser and Schoen showed that the solution of a certain extremal Steklov eigenvalue problem on a compact surface with boundary can be used to generate a free boundary minimal surface, i.e., a surface contained in the ball that has (i) zero mean curvature and (ii) meets the boundary of the ball orthogonally (doi:10.1007/s00222-015-0604-x). In this paper, we develop numerical methods that use this connection to realize free boundary minimal surfaces. Namely, on a compact surface, $Σ$, with genus $γ$ and $b$ boundary components, we maximize $σ_j(Σ,g) \ L(\partial Σ, g)$ over a class of smooth metrics, $g$, where $σ_j(Σ,g)$ is the $j$-th nonzero Steklov eigenvalue and $L(\partial Σ, g)$ is the length of $\partial Σ$. Our numerical method involves (i) using conformal uniformization of multiply connected domains to avoid explicit parameterization for the class of metrics, (ii) accurately solving a boundary-weighted Steklov eigenvalue problem in multi-connected domains, and (iii) developing gradient-based optimization methods for this non-smooth eigenvalue optimization problem. For genus $γ=0$ and $b=2,\dots, 9, 12, 15, 20$ boundary components, we numerically solve the extremal Steklov problem for the first eigenvalue. The corresponding eigenfunctions generate a free boundary minimal surface, which we display in striking images. For higher eigenvalues, numerical evidence suggests that the maximizers are degenerate, but we compute local maximizers for the second and third eigenvalues with $b=2$ boundary components and for the third and fifth eigenvalues with $b=3$ boundary components.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Consistency of archetypal analysis

Archetypal analysis is an unsupervised learning method that uses a convex polytope to summarize multivariate data. For fixed $k$, the method finds a convex polytope with $k$ vertices, called archetype points, such that the polytope is contained in the convex hull of the data and the mean squared distance between the data and the polytope is minimal. In this paper, we prove a consistency result that shows if the data is independently sampled from a probability measure with bounded support, then the archetype points converge to a solution of the continuum version of the problem, of which we identify and establish several properties. We also obtain the convergence rate of the optimal objective values under appropriate assumptions on the distribution. If the data is independently sampled from a distribution with unbounded support, we also prove a consistency result for a modified method that penalizes the dispersion of the archetype points. Our analysis is supported by detailed computational experiments of the archetype points for data sampled from the uniform distribution in a disk, the normal distribution, an annular distribution, and a Gaussian mixture model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spectral band degeneracies of $\fracπ{2}-$rotationally invariant periodic Schrödinger operators

The dynamics of waves in periodic media is determined by the band structure of the underlying periodic Hamiltonian. Symmetries of the Hamiltonian can give rise to novel properties of the band structure. Here we consider a class of periodic Schrödinger operators, $H_V=-Δ+V$, where $V$ is periodic with respect to the lattice of translates $Λ=\mathbb{Z}^2$. The potential is also assumed to be real-valued, sufficiently regular and such that, with respect to some origin of coordinates, inversion symmetric (even) and invariant under $π/2$ rotation. The present results are the $\mathbb{Z}^2-$ analogue of results obtained for conical degenerate points (Dirac points) in honeycomb structures. Our proofs make use of the framework developed by Fefferman-Weinstein and Fefferman-Lee-Thorp-Weinstein.

0 citations0 reviewsNo code linkedOpen access