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

Wenrui Hao

Wenrui Hao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical AnalysisMachine Learningmath.APInformation Theorymath.ITmath.OCQuantitative Methods

Trust snapshot

Quick read

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

preprint2026arXiv

A Structure-Preserving Scheme for the Time-Dependent Ginzburg-Landau Model with BCS Gap Coupling

We propose a structure-preserving scheme for a hybrid model that couples the time-dependent Ginzburg-Landau (TDGL) equation of superconducting vortex dynamics and the nonlinear Bardeen-Cooper-Schrieffer (BCS) gap equation. This formulation is consistent with the classical TDGL equation in the near-critical temperature, while extending the applicability of the existing TDGL model to regimes beyond the critical temperature. The resulting system poses significant computational challenges due to its nonlinear and coupled structure. To achieve stable and reliable simulations of the vortex dynamics and accompanying morphological transitions, we develop a maximum bound preserving, energy-stable implicit-explicit (IMEX) scheme. The structure-preserving properties of the scheme are rigorously established, ensuring long-time stability and physical consistency. Through two- and three-dimensional simulations, the hybrid model successfully captures the temporal and spatial formation and alignment of vortices and the suppression of superconductivity under increasing magnetic fields, demonstrating both the accuracy and robustness of the proposed computational approach.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

A Systematic Computational Framework for Practical Identifiability Analysis in Mathematical Models Arising from Biology

Practical identifiability is a critical concern in data-driven modeling of mathematical systems. In this paper, we propose a novel framework for practical identifiability analysis to evaluate parameter identifiability in mathematical models of biological systems. Starting with a rigorous mathematical definition of practical identifiability, we demonstrate its equivalence to the invertibility of the Fisher Information Matrix. Our framework establishes the relationship between practical identifiability and coordinate identifiability, introducing a novel metric that simplifies and accelerates the evaluation of parameter identifiability compared to the profile likelihood method. Additionally, we introduce new regularization terms to address non-identifiable parameters, enabling uncertainty quantification and improving model reliability. To guide experimental design, we present an optimal data collection algorithm that ensures all model parameters are practically identifiable. Applications to Hill functions, neural networks, and dynamic biological models demonstrate the feasibility and efficiency of the proposed computational framework in uncovering critical biological processes and identifying key observable variables.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

ZENN: A Thermodynamics-Inspired Computational Framework for Heterogeneous Data-Driven Modeling

Traditional entropy-based methods - such as cross-entropy loss in classification problems - have long been essential tools for representing the information uncertainty and physical disorder in data and for developing artificial intelligence algorithms. However, the rapid growth of data across various domains has introduced new challenges, particularly the integration of heterogeneous datasets with intrinsic disparities. To address this, we introduce a zentropy-enhanced neural network (ZENN), extending zentropy theory into the data science domain via intrinsic entropy, enabling more effective learning from heterogeneous data sources. ZENN simultaneously learns both energy and intrinsic entropy components, capturing the underlying structure of multi-source data. To support this, we redesign the neural network architecture to better reflect the intrinsic properties and variability inherent in diverse datasets. We demonstrate the effectiveness of ZENN on classification tasks and energy landscape reconstructions, showing its superior generalization capabilities and robustness-particularly in predicting high-order derivatives. ZENN demonstrates superior generalization by introducing a learnable temperature variable that models latent multi-source heterogeneity, allowing it to surpass state-of-the-art models on CIFAR-10/100, BBCNews, and AGNews. As a practical application in materials science, we employ ZENN to reconstruct the Helmholtz energy landscape of Fe$_3$Pt using data generated from density functional theory (DFT) and capture key material behaviors, including negative thermal expansion and the critical point in the temperature-pressure space. Overall, this work presents a zentropy-grounded framework for data-driven machine learning, positioning ZENN as a versatile and robust approach for scientific problems involving complex, heterogeneous datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nonlinear simulation of vascular tumor growth with chemotaxis and the control of necrosis

In this paper, we develop a sharp interface tumor growth model in two dimensions to study the effect of both the intratumoral structure using a controlled necrotic core and the extratumoral nutrient supply from vasculature on tumor morphology. We first show that our model extends the benchmark results in the literature using linear stability analysis. Then we solve this generalized model numerically using a spectrally accurate boundary integral method in an evolving annular domain, not only with a Robin boundary condition on the outer boundary for the nutrient field which models tumor vasculature, but also with a static boundary condition on the inner boundary for pressure field which models the control of tumor necrosis. The discretized linear systems for both pressure and nutrient fields are shown to be well-conditioned through tracing GMRES iteration numbers. Our nonlinear simulations reveal the stabilizing effects of angiogenesis and the destabilizing ones of chemotaxis and necrosis in the development of tumor morphological instabilities if the necrotic core is fixed in a circular shape. When the necrotic core is controlled in a non-circular shape, the stabilizing effects of proliferation and the destabilizing ones of apoptosis are observed. Finally, the values of the nutrient concentration with its fluxes and the pressure level with its normal derivatives, which are solved accurately at the boundaries, help us to characterize the corresponding tumor morphology and the level of the biophysical quantities on interfaces required in keeping various shapes of the necrotic region of the tumor. Interestingly, we notice that when the necrotic region is fixed in a 3-fold non-circular shape, even if the initial shape of the tumor is circular, the tumor will evolve into a shape corresponding to the 3-fold symmetry of the shape of the fixed necrotic region.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Weight Initialization Based on the Linear Product Structure for Neural Networks

Weight initialization plays an important role in training neural networks and also affects tremendous deep learning applications. Various weight initialization strategies have already been developed for different activation functions with different neural networks. These initialization algorithms are based on minimizing the variance of the parameters between layers and might still fail when neural networks are deep, e.g., dying ReLU. To address this challenge, we study neural networks from a nonlinear computation point of view and propose a novel weight initialization strategy that is based on the linear product structure (LPS) of neural networks. The proposed strategy is derived from the polynomial approximation of activation functions by using theories of numerical algebraic geometry to guarantee to find all the local minima. We also provide a theoretical analysis that the LPS initialization has a lower probability of dying ReLU comparing to other existing initialization strategies. Finally, we test the LPS initialization algorithm on both fully connected neural networks and convolutional neural networks to show its feasibility, efficiency, and robustness on public datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Learn bifurcations of nonlinear parametric systems via equation-driven neural networks

Nonlinear parametric systems have been widely used in modeling nonlinear dynamics in science and engineering. Bifurcation analysis of these nonlinear systems on the parameter space are usually used to study the solution structure such as the number of solutions and the stability. In this paper, we develop a new machine learning approach to compute the bifurcations via so-called equation-driven neural networks (EDNNs). The EDNNs consist of a two-step optimization: the first step is to approximate the solution function of the parameter by training empirical solution data; the second step is to compute bifurcations by using the approximated neural network obtained in the first step. Both theoretical convergence analysis and numerical implementation on several examples have been performed to demonstrate the feasibility of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Power Series Expansion Neural Network

In this paper, we develop a new neural network family based on power series expansion, which is proved to achieve a better approximation accuracy in comparison with existing neural networks. This new set of neural networks embeds the power series expansion (PSE) into the neural network structure. Then it can improve the representation ability while preserving comparable computational cost by increasing the degree of PSE instead of increasing the depth or width. Both theoretical approximation and numerical results show the advantages of this new neural network.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Two neural-network-based methods for solving obstacle problems

Two neural-network-based numerical schemes are proposed to solve the classical obstacle problems. The schemes are based on the universal approximation property of neural networks, and the cost functions are taken as the energy minimization of the obstacle problems. We rigorously prove the convergence of the two schemes and derive the convergence rates with the number of neurons $N$. In the simulations, we use two example problems (1-D & 2-D) to verify the convergence rate of the methods and the quality of the results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An adaptive homotopy method for computing bifurcations of nonlinear parametric systems

In this paper, we present an adaptive step-size homotopy tracking method for computing bifurcation points of nonlinear systems. There are four components in this new method: 1) an adaptive tracking technique is developed near bifurcation points; 2) an inflation technique is backed up when the adaptive tracking fails; 3) Puiseux series interpolation is used to compute bifurcation points; and 4)the tangent cone structure of the bifurcation point is approximated numerically to compute solutions on different branches. Various numerical examples of nonlinear systems are given to illustrate the efficiency of this new approach. This new adaptive homotopy tracking method is also applied to a system of nonlinear PDEs and shows robustness and efficiency for large-scale nonlinear discretized systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bifurcation analysis of a free boundary model of plaque formation associated with the cholesterol ratio

The low-density lipoprotein (LDL)/high-density lipoprotein (HDL)-cholesterol ratio has been shown a high correlation with the cardiovascular risk assessment. Is it possible to quantify the correlation mathematically? In this paper, we develop a bifurcation analysis for a mathematical model of the plaque formation with a free boundary in the early stage of atherosclerosis. This bifurcation analysis, to the ratio of LDL/HDL, is based on the explicit formulation of radially symmetric steady-state solutions. By performing the perturbation analysis to these solutions, we establish the existence of bifurcation branches and derive a theoretical condition that a bifurcation occurs for different modes. Then we also analyze the stability of radially symmetric steady-state solutions and conduct numerical simulations to verify all the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A homotopy training algorithm for fully connected neural networks

In this paper, we present a Homotopy Training Algorithm (HTA) to solve optimization problems arising from fully connected neural networks with complicated structures. The HTA dynamically builds the neural network starting from a simplified version to the fully connected network via adding layers and nodes adaptively. Therefore, the corresponding optimization problem is easy to solve at the beginning and connects to the original model via a continuous path guided by the HTA, which provides a high probability to get a global minimum. By gradually increasing the complexity of model along the continuous path, the HTA gets a rather good solution to the original loss function. This is confirmed by various numerical results including VGG models on CIFAR-10. For example, on the VGG13 model with batch normalization, HTA reduces the error rate by 11.86% on test dataset comparing with the traditional method. Moreover, the HTA also allows us to find the optimal structure for a fully connected neural network by building the neutral network adaptively.

0 citations0 reviewsNo code linkedOpen access