Haoran Wang
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Haoran Wang has a total of 67 co-authors. Here are the top co-authors by number of shared publications:
- Huck Beng Chew
Department of Aerospace Engineering, University of Illinois, Urbana-Champaign , Urbana, Illinois 61801, USA
5 shared publications Xueju Wang
Department of Materials Science and Engineering, University of Connecticut, Storrs, CT 06269, United States of America
4 shared publicationsShuman Xia
Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 30332, United States of America
3 shared publicationsXinnan Wang
3 shared publications- Byeongdu Lee
Argonne National Laboratory
2 shared publications Huck Beng Chew
Department of Aerospace Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
2 shared publications
Latest Publications
- 17 September 2025
- Published by The Royal Society in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 481 (2322)
Learning the governing partial differential equation of solid diffusion from atomistic simulations
Governing partial differential equations (PDEs) play a critical role in materials research and applications, as they describe essential physics underlying materials behaviour. Traditionally, these equations are developed through phenomenological modelling of experimental results or first principle analysis based on conservation laws. In addition, molecular dynamics (MD) simulations capture atomistic-scale behaviour with detailed physics. However, translating atomistic insights into continuum-scale governing equations remains a significant challenge. Empowered by recent advances in data-driven modelling, we develop a computational framework to learn governing PDEs directly from atomistic simulation data. The framework integrates numerical differentiation of MD data with the identification of constitutive relationships. It proves effective and efficient in learning governing PDEs from noisy and limited MD datasets, without requiring prior knowledge of the final PDEs. Using this framework, we identify a nonlinear PDE governing solid-state diffusion in nickel–hydrogen alloys. This PDE reveals a highly concentration-dependent diffusivity that varies over an order of magnitude. Our data-driven computational framework paves the way for cross-scale constitutive modelling.
- 17 September 2025
- Published by The Royal Society in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 481 (2322)
Learning the governing partial differential equation of solid diffusion from atomistic simulations
Governing partial differential equations (PDEs) play a critical role in materials research and applications, as they describe essential physics underlying materials behaviour. Traditionally, these equations are developed through phenomenological modelling of experimental results or first principle analysis based on conservation laws. In addition, molecular dynamics (MD) simulations capture atomistic-scale behaviour with detailed physics. However, translating atomistic insights into continuum-scale governing equations remains a significant challenge. Empowered by recent advances in data-driven modelling, we develop a computational framework to learn governing PDEs directly from atomistic simulation data. The framework integrates numerical differentiation of MD data with the identification of constitutive relationships. It proves effective and efficient in learning governing PDEs from noisy and limited MD datasets, without requiring prior knowledge of the final PDEs. Using this framework, we identify a nonlinear PDE governing solid-state diffusion in nickel–hydrogen alloys. This PDE reveals a highly concentration-dependent diffusivity that varies over an order of magnitude. Our data-driven computational framework paves the way for cross-scale constitutive modelling.
- 9 August 2024
- Published by Royal Society of Chemistry (RSC) in Journal of Materials Chemistry A
- Vol. 12 (37) , 24886-24895
Microscale mechanical property variations of Al-substituted LLZO: insights from compression testing and molecular dynamics simulations
An organophosphine-controlled diversity-oriented synthesis of chromone inden-1-one-fused cyclopentadienylides and C-acylated 2-((chromone-3-yl)methylene)-indandiones is reported. Key attributes of the methodology are the in situ generation of an allylic P-ylide and subsequent regio- and chemoselective intramolecular cyclization reactions that preferentially result in the aforementioned chromone adducts.
- 9 August 2024
- Published by Royal Society of Chemistry (RSC) in Journal of Materials Chemistry A
- Vol. 12 (37) , 24886-24895
Microscale mechanical property variations of Al-substituted LLZO: insights from compression testing and molecular dynamics simulations
An organophosphine-controlled diversity-oriented synthesis of chromone inden-1-one-fused cyclopentadienylides and C-acylated 2-((chromone-3-yl)methylene)-indandiones is reported. Key attributes of the methodology are the in situ generation of an allylic P-ylide and subsequent regio- and chemoselective intramolecular cyclization reactions that preferentially result in the aforementioned chromone adducts.
- 12 February 2024
- Published by Elsevier in International Journal of Solids and Structures
- Vol. 291
All-atom molecular dynamics simulations of nematic liquid crystal elastomers
- 12 February 2024
- Published by Elsevier in International Journal of Solids and Structures
- Vol. 291
All-atom molecular dynamics simulations of nematic liquid crystal elastomers
- Open Access
Tolerance For the Pixelation Effect in Shear Measurement
Images taken by space telescopes typically have a superb spatial resolution, but a relatively poor sampling rate due to the finite CCD pixel size. Beyond the Nyquist limit, it becomes uncertain how much the pixelation effect may affect the accuracy of galaxy shape measurement. It is timely to study this issue given that a number of space-based large-scale weak lensing surveys are planned. Using the Fourier_Quad method, we quantify the shear recovery error as a function of the sampling factor Q, i.e., the ratio between the FWHM of the point-spread-function (PSF) and the pixel size of the CCD, for different PSFs and galaxies of different sizes and noise levels. We show that sub-percent-level accuracy in shear recovery is achievable with single-exposure images for $Q\lesssim 2$. The conclusion holds for galaxies much smaller than the PSF, and those with a significant level of noise.Open AccessTolerance For the Pixelation Effect in Shear Measurement
Images taken by space telescopes typically have a superb spatial resolution, but a relatively poor sampling rate due to the finite CCD pixel size. Beyond the Nyquist limit, it becomes uncertain how much the pixelation effect may affect the accuracy of galaxy shape measurement. It is timely to study this issue given that a number of space-based large-scale weak lensing surveys are planned. Using the Fourier_Quad method, we quantify the shear recovery error as a function of the sampling factor Q, i.e., the ratio between the FWHM of the point-spread-function (PSF) and the pixel size of the CCD, for different PSFs and galaxies of different sizes and noise levels. We show that sub-percent-level accuracy in shear recovery is achievable with single-exposure images for $Q\lesssim 2$. The conclusion holds for galaxies much smaller than the PSF, and those with a significant level of noise. - 22 October 2021
- Published by IOP Publishing in Journal of Physics D: Applied Physics
- Vol. 55 (6) , 063001
A review of the multiscale mechanics of silicon electrodes in high-capacity lithium-ion batteries
Over the past decade, there has been a significant advancement in understanding the mechanics of silicon (Si) electrodes in lithium (Li)-ion batteries. Much of this interest in Si electrodes as ideal anode materials for high-capacity Li-ion batteries stems from its theoretical specific capacity of 4200 mAh g−1, which is an order-of-magnitude higher than that of conventional graphite electrodes (372 mAh g−1). However, the high capacity of Li ions is also accompanied by a ∼300% volume expansion of the Si electrode during Li intercalation, which results in massive cracking of the electrode and capacity fade. In this review article, we summarize recent progress in elucidating the underlying fracture and failure mechanics of Si electrodes using multiscale computations and experiments, spanning the quantum, atomistic, microscopic, and macroscopic length scales. We focus on four fundamental mechanics issues: (i) the mechanical properties and fracture behavior of lithiated Si electrodes; (ii) the interfacial mechanics between Si thin-film electrodes and current collectors; (iii) the deformation and failure mechanics of the solid electrolyte interphase; and (iv) the design of Si electrodes for improved mechanical performance. Current challenges and possible future directions for the field of mechanics of materials in pursuit of high-capacity rechargeable batteries are also discussed.
- 22 October 2021
- Published by IOP Publishing in Journal of Physics D: Applied Physics
- Vol. 55 (6) , 063001
A review of the multiscale mechanics of silicon electrodes in high-capacity lithium-ion batteries
Over the past decade, there has been a significant advancement in understanding the mechanics of silicon (Si) electrodes in lithium (Li)-ion batteries. Much of this interest in Si electrodes as ideal anode materials for high-capacity Li-ion batteries stems from its theoretical specific capacity of 4200 mAh g−1, which is an order-of-magnitude higher than that of conventional graphite electrodes (372 mAh g−1). However, the high capacity of Li ions is also accompanied by a ∼300% volume expansion of the Si electrode during Li intercalation, which results in massive cracking of the electrode and capacity fade. In this review article, we summarize recent progress in elucidating the underlying fracture and failure mechanics of Si electrodes using multiscale computations and experiments, spanning the quantum, atomistic, microscopic, and macroscopic length scales. We focus on four fundamental mechanics issues: (i) the mechanical properties and fracture behavior of lithiated Si electrodes; (ii) the interfacial mechanics between Si thin-film electrodes and current collectors; (iii) the deformation and failure mechanics of the solid electrolyte interphase; and (iv) the design of Si electrodes for improved mechanical performance. Current challenges and possible future directions for the field of mechanics of materials in pursuit of high-capacity rechargeable batteries are also discussed.
- Open Access
- 8 September 2021
- Published by The Royal Society in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 477 (2253) , 20210259
Nematic liquid crystalline elastomers are aeolotropic materials
Continuum models describing ideal nematic solids are widely used in theoretical studies of liquid crystal elastomers. However, experiments on nematic elastomers show a type of anisotropic response that is not predicted by the ideal models. Therefore, their description requires an additional term coupling elastic and nematic responses, to account for aeolotropic effects. In order to better understand the observed elastic response of liquid crystal elastomers, we analyse theoretically and computationally different stretch and shear deformations. We then compare the elastic moduli in the infinitesimal elastic strain limit obtained from the molecular dynamics simulations with the ones derived theoretically, and show that they are better explained by including nematic order effects within the continuum framework.
- 8 September 2021
- Published by The Royal Society in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
- Vol. 477 (2253) , 20210259
Open AccessNematic liquid crystalline elastomers are aeolotropic materials
Continuum models describing ideal nematic solids are widely used in theoretical studies of liquid crystal elastomers. However, experiments on nematic elastomers show a type of anisotropic response that is not predicted by the ideal models. Therefore, their description requires an additional term coupling elastic and nematic responses, to account for aeolotropic effects. In order to better understand the observed elastic response of liquid crystal elastomers, we analyse theoretically and computationally different stretch and shear deformations. We then compare the elastic moduli in the infinitesimal elastic strain limit obtained from the molecular dynamics simulations with the ones derived theoretically, and show that they are better explained by including nematic order effects within the continuum framework.
- 1 September 2019
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 354, 37-55
Modeling uncertainties in molecular dynamics simulations using a stochastic reduced-order basis
- 1 September 2019
- Published by Elsevier in Computer Methods in Applied Mechanics and Engineering
- Vol. 354, 37-55
Modeling uncertainties in molecular dynamics simulations using a stochastic reduced-order basis
- 12 March 2019
- Published by American Chemical Society (ACS) in ACS Applied Materials & Interfaces
- Vol. 11 (12) , 11715-11721
Nanofibrillar Si Helices for Low-Stress, High-Capacity Li+ Anodes with Large Affine Deformations
We report on the chemical lithiation of long microscale helices composed of densely packed amorphous silicon (aSi) nanofibrils, fabricated by glancing angle deposition (GLAD) through e-beam evaporation. In situ electron microscopy and companion finite element modeling demonstrate that the nanofibrillar structure of the aSi helices allows for 2 orders of magnitude faster effective rates for Li diffusion (D0 = 10–10 cm2/s) compared to solid aSi nanowires, while also averting fragmentation during lithiation. More importantly, it is shown that specific helical geometries can accommodate large, lithium-induced, volumetric expansions without shape distortion. A major advantage of the helical nanostructures is that the compressive force generated due to lithiation-induced expansion is an order of magnitude smaller than in straight nanocolumns that permanently buckle during lithiation. Thus, GLAD-fabricated films composed of dense periodic microscale helices with properly designed coil geometries are highly suitable for robust, high-capacity Li+ anodes.
- 12 March 2019
- Published by American Chemical Society (ACS) in ACS Applied Materials & Interfaces
- Vol. 11 (12) , 11715-11721
Nanofibrillar Si Helices for Low-Stress, High-Capacity Li+ Anodes with Large Affine Deformations
We report on the chemical lithiation of long microscale helices composed of densely packed amorphous silicon (aSi) nanofibrils, fabricated by glancing angle deposition (GLAD) through e-beam evaporation. In situ electron microscopy and companion finite element modeling demonstrate that the nanofibrillar structure of the aSi helices allows for 2 orders of magnitude faster effective rates for Li diffusion (D0 = 10–10 cm2/s) compared to solid aSi nanowires, while also averting fragmentation during lithiation. More importantly, it is shown that specific helical geometries can accommodate large, lithium-induced, volumetric expansions without shape distortion. A major advantage of the helical nanostructures is that the compressive force generated due to lithiation-induced expansion is an order of magnitude smaller than in straight nanocolumns that permanently buckle during lithiation. Thus, GLAD-fabricated films composed of dense periodic microscale helices with properly designed coil geometries are highly suitable for robust, high-capacity Li+ anodes.
- 1 September 2018
- Published by Institute of Electrical and Electronics Engineers (IEEE)
- Vol. 2, 403-406
Two Dimensional Distribution of Metallic (Cu) Vapor in a Forced Vacuum Arc Extinction by LIF
It is important to observe the metallic vapor dissipation after interrupting a vacuum arc, for getting a deeper understanding of dielectric recovery processes, especially for DC interruptions. The object of this paper is to obtain two dimensional density distribution of copper vapor in a forced vacuum arc extinction. The method of laser-induced fluorescence (LIF) using a pulse laser with a central wavelength at 324.8 nm was adopted here. Meanwhile, an ICCD was used to capture the 510.6 nm fluorescence signal with an exposure time of 50 ns. The fluorescence can be distinguished from 125μs (the moment of arc extinction) to 375μs with a time interval of 25μs. As a result, we could obtain the relative density distribution of ground-state copper metal vapor after the arc extinction. The results show that on the condition of CuCr25 contacts the metal vapor was evenly distributed in the gap and rapidly spread over time. The density of the copper vapor between the contacts dropped to less than 10% of the initial stage 200μs after arc extinction. On the contrary, on the condition of CuCr50, higher copper metal vapor density at the moment of arc extinction was observed. It seems that the higher Cr content in the Cu-Cr contacts will lead to producing more metal vapor around current zero.
- 1 September 2018
- Published by Institute of Electrical and Electronics Engineers (IEEE)
- Vol. 2, 403-406
Two Dimensional Distribution of Metallic (Cu) Vapor in a Forced Vacuum Arc Extinction by LIF
It is important to observe the metallic vapor dissipation after interrupting a vacuum arc, for getting a deeper understanding of dielectric recovery processes, especially for DC interruptions. The object of this paper is to obtain two dimensional density distribution of copper vapor in a forced vacuum arc extinction. The method of laser-induced fluorescence (LIF) using a pulse laser with a central wavelength at 324.8 nm was adopted here. Meanwhile, an ICCD was used to capture the 510.6 nm fluorescence signal with an exposure time of 50 ns. The fluorescence can be distinguished from 125μs (the moment of arc extinction) to 375μs with a time interval of 25μs. As a result, we could obtain the relative density distribution of ground-state copper metal vapor after the arc extinction. The results show that on the condition of CuCr25 contacts the metal vapor was evenly distributed in the gap and rapidly spread over time. The density of the copper vapor between the contacts dropped to less than 10% of the initial stage 200μs after arc extinction. On the contrary, on the condition of CuCr50, higher copper metal vapor density at the moment of arc extinction was observed. It seems that the higher Cr content in the Cu-Cr contacts will lead to producing more metal vapor around current zero.