Postdoctoral Researcher — University of Illinois Urbana-Champaign
Gravitational waves · Neutron star mergers · Modified gravity
Contact
I come from Hsinchu, a county in northern Taiwan, where I completed my Bachelor's degree with double majors in Physics and Mathematics. I then joined the group of Prof. Chao-Qiang Geng in Hsinchu for doctoral training. In 2019, I pursued a joint PhD between Nat'l Tsing Hua University (Taiwan) and the University of Tübingen (Germany), working in Prof. Kostas Kokkotas' group. I defended my thesis at the end of 2022 and received the degree in early 2023. I subsequently joined the Albert Einstein Institute in Golm under Prof. Masaru Shibata. I am currently a postdoctoral researcher at the University of Illinois Urbana-Champaign, working with Prof. Nicolas Yunes.
I study what neutron stars are made of, and whether gravity behaves as general relativity predicts, through the signals these stars emit before, during, and after coalescence. The work spans relativistic stellar oscillations, quasiequilibrium initial data, fully relativistic merger simulations, waveform construction, and the interpretation of high-energy transients. In each case the aim is the same: identify a physical effect, predict its observable signature, quantify the modelling uncertainty that would obscure it, and establish what a detection — or a principled nondetection — would teach us.
During inspiral the orbital frequency sweeps through a neutron star's oscillation spectrum, and resonant excitation of g-, f-, and shear modes can stress the crust to failure, redistribute angular momentum, and amplify the magnetic field. I developed the general-relativistic framework for these resonances and use precursor flare timings — including quasi-periodic features observed in short gamma-ray burst precursors — to perform asteroseismology on the flaring star, constraining its spin and the nuclear equation of state.
Recent nonlinear simulations further reveal mode-driven differential rotation, parametric instability, and spin-up, strengthening the case for premerger magnetic amplification.
Numerical accuracy is part of the physics: a waveform cannot reveal a subtle tide or a modification of gravity unless its error budget is smaller than the effect. I quantify uncertainties from grid resolution, wave extraction, initial-data solver, and evolution code in binary neutron star simulations, and use fully relativistic evolutions of tidal resonance to test where linear-response descriptions of the dynamical tide cease to hold. Residuals against analytic approximants then become diagnostics of missing physics rather than of numerical noise.
Dynamical compact stars probe gravity in regimes inaccessible to weak-field experiments. I study scalarization and descalarization in massive scalar-tensor and scalar-Gauss-Bonnet theories: their formation through core collapse and accretion, their nonlinear stability and nonuniqueness, and their imprint on quasiequilibrium binaries, merger waveforms, and postmerger remnants. Because a scalar field and an uncertain equation of state can shift the same observables, I emphasize combining inspiral, collapse threshold, and postmerger spectrum over any single anomalous frequency.
The equation of state enters as a set of measurable consequences rather than a label: composition gradients set buoyancy and the g-mode spectrum, elasticity determines whether motion penetrates or fractures the crust, tidal deformability shapes the inspiral, and thermal pressure and differential rotation govern whether a massive remnant survives. I follow these observables through a common stellar model so that degeneracies with spin, temperature, magnetic field, and the underlying theory of gravity can be cross-checked rather than assumed away.
A Kadath-based numerical-relativity code and sibling of FUKa, optimized for binary neutron star configurations. Features a 10× speedup, 70% memory reduction via a sparse matrix solver, and a new 2D axisymmetric solver.
Learn moreGeneral Relativistic Axisymmetric Spacetime Solver. Computes uniformly and differentially rotating neutron-star equilibria on a compactified meridional grid under General Relativity or scalar-tensor gravity, with tabulated EOS handling and adaptive field relaxation.
Learn more