Polars (AM Her stars)
Cataclysmic variables whose white-dwarf field — 10 to 240 MG — is strong enough to lock the binary into synchronous rotation and suppress the accretion disk entirely. The gas free-falls along field lines onto one magnetic pole.
Astrophysicist · AI / Agent Engineer
I study magnetic cataclysmic variables — binaries in which a strongly magnetised white dwarf strips gas from its companion and funnels it down field lines onto a single pole. I weigh those fields by modelling the cyclotron humps that the shocked, gyrating electrons stamp onto the spectrum.
Explore five interactive figures: Roche geometry and accretion streams, cyclotron spectra, disk eclipses, binary evolution on a cluster clock, and black-hole superradiance. Each figure connects the visual model to its physical assumptions, with observed constraints and schematic stages identified alongside it.
Ph.D. candidate, Physics & Astronomy, Sun Yat-sen University (Zhuhai) · advised by Prof. Chengyuan Li · CSST Science Center for the Greater Bay Area.
How do you weigh the magnetic field of a star you cannot resolve? You read it off the light.
Cataclysmic variables whose white-dwarf field — 10 to 240 MG — is strong enough to lock the binary into synchronous rotation and suppress the accretion disk entirely. The gas free-falls along field lines onto one magnetic pole.
At the accretion shock, hot electrons gyrate and radiate at harmonics of νc = eB/2πmec. The spacing of the resulting humps is a direct, nearly model-free field measurement.
A Koester photosphere, a hot-spot blackbody and an isothermal constant-Λ cyclotron slab, fitted simultaneously. Profiled χ²(B) curves expose the competing harmonic assignments instead of hiding them behind one number.
Six to seven years of ZTF photometry, mined for orbital periods and their aliases — which is how a rare short-period eclipsing-polar candidate below the period gap turned up.
CNN and SVM classifiers, with engineered period and shape features plus DBSCAN clustering, sift ~10⁸ ZTF light curves down to a short, vettable candidate list: magnetic CVs, eclipsing systems and other oddities.
Flat-bottomed “square-wave” eclipses in which a misaligned circumbinary disk occults a low-mass binary, à la KH 15D. Open-cluster membership turns these into age-dated probes of long-lived circumstellar dust.
LLM-agent pipelines that wrap retrieval, modelling and evidence-auditing into a reviewable graph — RAG plus knowledge graphs, tool-calling, and domain-tuned models that refuse to publish a claim the evidence doesn't support.
See the AI / Agents track →A polar is synchronous — the white dwarf's spin is locked to the orbit — so the whole magnetic geometry is rigid in the co-rotating frame and only the observer moves. That is exactly how the figure is built: the donor's Roche lobe, the ballistic stream out of L1 and the dipole accretion curtain are solved once, then viewed from a direction that rotates with phase. The light curve is not drawn separately; it is the same geometry, integrated.
Hover to read a phase; drag either panel to move the scene with it. The shaded band is the eclipse window implied by the current inclination — it appears on its own once i passes ic.
F ∝ sin²µ · vis, cos µ = B̂spot · n̂(φ)
Cyclotron radiation leaves perpendicular to B, so the flux peaks when the field at the footpoint lies across the line of sight. B̂ is the true dipole direction there, not the pole axis — at the footpoint the field is tilted out of the radial by roughly 14°.
Φ = −µ₁/r₁ − µ₂/r₂ − ½[(x−xcm)² + y²]
The teardrop is the equipotential through L1, found by bisection along 420 directions. Change q and L1 slides, the lobe changes shape, and the stream — integrated with RK4 in the rotating frame — takes a different Coriolis curve. Nothing here is a drawn ellipse.
ic = arccos [(RL,proj + Rwd) / a]
There is no eclipse toggle. Raise i and the projected lobe eventually covers the white dwarf on the sky; the shaded window in the light curve appears at the same moment, because both read the same projected silhouette.
The spectrum of a polar is a sum, not a curve: a white-dwarf photosphere, a heated pole cap, and an isothermal cyclotron slab. Only the third one knows about B, and it only shows resolved humps where it has gone optically thin. Drag Λ — the slab's column parameter — and watch the turnover harmonic march up the ladder. That drift is the whole reason harmonic assignment is ambiguous, and why the paper reports competing branches rather than a single field strength.
λn ≈ 1.071×10⁶ / (n · B[MG]) Å
Consecutive humps are spaced by Δλ/λ ≈ 1/n, so the pattern compresses toward the blue. Measuring the spacing, not the absolute position, is what makes B nearly model-independent — provided you have the right n.
I = Bλ(Te) [1 − e−Στn]
Not a sum of Gaussians. Low harmonics saturate at the source function and merge into a smooth blue continuum; only where τn < 1 do separate humps survive. The inset bars show log τn and where it crosses 1.
τn ∝ Λ n2n/(2ⁿn!) · βn−1 sin2n−2θ
Raise Λ or kTe and the turnover moves up, so the same observed humps can be labelled n or n+1 — and B shifts by tens of MG. My fitting code profiles χ²(B) across those branches and reports them, instead of quietly picking one.
UPK 13-c2 dims into a flat-bottomed square wave every 36.71 days. A misaligned circumbinary disk presents a sharp inner edge, and only the wider-orbit component swings behind it — that star carries ~40% of the light, so the floor lands near 60% while its companion keeps shining. The clincher is the slow, multi-day ingress: switch the occulted source to a white dwarf and the walls snap vertical, which is exactly why a WD is ruled out.
Dots = real ZTF g/r photometry, phase-folded on P = 36.710719 d. Dashed = the χ²-fitted trapezoid. Solid = the live model — it only leaves the fit when you change the occulted source or the geometry.
v⊥ = vorb sin α, ting ≈ 2R★/v⊥
A grazing geometry stretches a stellar-diameter crossing into days. Shrink α and the walls flatten further; a point-like occulted source removes R★ from the numerator entirely and the eclipse becomes a step function.
Rin ≈ 2–3 a ≈ 130 R⊙
For a ≈ 0.24 AU (Kepler, Mtot ≈ 1.4 M⊙), tidal truncation (Artymowicz & Lubow 1994) clears the disk inside ~2.5 a. The ring in the figure is drawn at that radius, not at a convenient one.
Fmin = 1 − Locc/Ltot ≈ 0.60
A single star behind an opaque edge cannot produce a 40% flat bottom unless something else is still shining. The near-achromatic depth from the optical through W1, plus a mid-IR excess, is what points to dust rather than a stellar companion.
One initial binary, three possible outcomes. Too little envelope ejection ends in a merger; a wider surviving orbit can run out of time before accretion begins. Follow the red NS Per track to a dwarf nova within the 209 Myr cluster age, then explore the complementary merger and survival channels in HSC 1224.
Reduced motion is enabled. Select individual stages or use the timeline to explore.
The representative track reaches contact before the cluster age. At 209 Myr its period is about 6.208 h; the observed period is 6.29500 ± 0.00050 h.
The critical Roche surfaces meet at L₁. A stream appears only after the donor fills this surface. Its ballistic path is traced to impact or first closest approach, with all structures rotating in the same orbital plane.
HSC 1224: object identities follow the 12 September 2026 thesis revision. BSS-1 is near Roche contact: transfer may have just ended or continue slowly. WD4 and WD6 are independent merger-remnant candidates; WD2 is the double-white-dwarf candidate in the wider branch. These are different cluster members, not the future identities of BSS-1 or BSS-2. The BSS-2 orbit and the WD2 orbit remain unconfirmed.
NS Per: N_representative_track, 11 September 2026; production grid results_big, binary 143961, Z = 0.0132. BPP event nodes preserve the instantaneous CE jump; BCM samples are interpolated consistently for the scene, readouts and cursor. The 209 Myr endpoint is interpolated between the 200 and 210 Myr samples.
Geometry: the synchronous Roche potential; ballistic motion follows Newton's equations in the rotating frame. Gas pressure, shocks, magnetism and envelope hydrodynamics are not simulated.
The machine-learning and LLM side of this work — a locally fine-tuned 8B astronomy model (QLoRA + AttnRes), a knowledge-graph & experiment-design agent, the Astro Agent research pipeline, and the ZTF CNN classifier — lives on its own page, each project fully visualised.
Interactive figure 05 · General relativity
A rotating black hole can transfer energy to a wave. Can that wave remain trapped and keep growing? Explore the distinction at the heart of my work on superradiant stability.
Preparing the black-hole view…
Camera and zoom controls remain available.Follow the far side of the disk as gravity bends its light around the shadow.
A nonrotating optical view with an illustrative disk. Explore the rotating black holes and superradiance from my paper below.
The colors trace a field and energy transfer, not visible light or an accretion disk. Wave amplitudes and animation time are illustrative.
G = c = ℏ = 1 · m = 1 · a = M. The black hole is exactly extremal; only the scalar-field parameters vary.
The superradiant and bound-state conditions hold, but the paper excludes an exterior trapping well in this parameter range.
The threshold fixes the sign of the horizon energy flux. It does not give an amplification percentage.
This is required for a massive bound state. By itself, it does not establish a trapping well or a growing mode.
Within the superradiant window, this sufficient condition rules out the feedback needed for this instability.
Read the boundary
The green region meets the paper's sufficient stability condition within the superradiant and bound-state window. Moving outside green does not establish instability.
The hatched region needs further mode analysis. This figure evaluates the analytic inequalities; it does not compute a mode spectrum, a growth rate, or the evolution of a boson cloud.
What is calculated. The parameter map uses the strict inequalities in Lin et al. (2021), Eqs. (1), (16), (35), (48), (49), (57), and (68). Frequencies and scalar masses are expressed as Mω and Mμ, with G = c = ℏ = 1. The azimuthal number is fixed at m = 1; this does not select a radial or angular eigenmode. Equality boundaries are excluded from the stated sufficient proof.
Charged case. Extremality enforces Q/M = √(1 − k²), with positive Q and 0 ≤ k < 1. The qQ slider changes the scalar charge coupling. The sufficient condition requires qQ > 0; the neutral or oppositely charged cases are not certified by that condition. The superradiant threshold is Mωc = (mk + qQ)/(1 + k²). Charge can supply energy as well as rotation.
Optical view. The large view integrates Schwarzschild null geodesics with a fourth-order Runge–Kutta scheme, in units G = c = M = 1. The horizon is at r = 2M; the thin disk extends from 6M to 22M. A static observer at 36M uses a local orthonormal camera frame. The browser uses WebGL when available and otherwise a Canvas 2D renderer with cached ray geometry. Both use the same optical model. The software view first shows a fast preview, then refines it at the display resolution after the camera stops moving. The disk texture, color mapping, exposure and animation timing are illustrative; finite integration steps limit the higher-order images. It is not a Kerr image, an accretion-flow simulation or an image predicted by the superradiance paper. The color-shift toggle includes the frequency shift for circular disk motion and gravitational redshift, with artistic tone mapping and a small photographic glow.
Wave illustration. The smaller panel is a local field schematic in a compressed radial layout, not ray tracing or a solution of the Klein–Gordon equation. No wave is drawn escaping from inside the horizon. For a bound-state choice, the blue field is local rather than an incoming propagating wave from infinity. No trapping well is invented for the unclassified region.
Figure 5 · Black-hole superradiance and stability. A Schwarzschild optical view introduces light bending. A separate schematic explains energy exchange in the rotating models; the linked map evaluates the sufficient stability bounds from my 2021 paper. Green is a proved region under the stated assumptions. Hatching means the present criterion does not decide.
Selected work, newest first.
The Astrophysical Journal (accepted) · partly assisted by Astro Agent
UPK 13-c2 shows a flat-bottomed 36.71-day eclipse with a near-achromatic ~40% decrement from the optical through W1 and a mid-IR excess. A slow 2.5-day ingress rules out a white-dwarf occulter, favouring a late-K/early-M binary behind a misaligned circumbinary disk — potentially the oldest known main-sequence disk-eclipsing binary.
See Fig. 3 →The Astrophysical Journal (under review)
Simultaneous continuum + cyclotron fitting of DESI/LAMOST polar candidates yields a robust 56.8 MG field for J0005+2941 — a rare eclipsing-polar candidate below the period gap — and reports explicit harmonic-branch constraints rather than hiding the ambiguous cases.
See Fig. 2 →Astronomy & Astrophysics (accepted)
ROSAT X-rays, ZTF gri photometry, and Palomar 200-inch time-resolved spectroscopy with Doppler tomography confirm ZTF J0112+5827 as a polar: an 80.9-min magnetic CV showing accretion streams (no disk) and a white-dwarf field of 38.7 MG measured from its cyclotron harmonics.
See Fig. 1 →The Astrophysical Journal Letters (accepted)
A young, strongly magnetic white dwarf whose field and rotation point to a recent double-degenerate merger origin.
The Astrophysical Journal Supplement Series (accepted)
A spectroscopically identified DESI DR1 sample of cataclysmic variables and their population properties.
Physics Letters B 819, 136392 (2021)
Analytic superradiant-stability conditions for four-dimensional extremal Kerr and Kerr–Newman black holes under massive (and charged-massive) scalar perturbations — e.g. an extremal Kerr black hole is superradiantly stable when ω < μ/√3.
Explore the interactive stability map →A first-author methods paper in College Physics (digital holographic interferometry for liquid-phase diffusion, accepted), plus further third/fourth-author papers in ApJ and related journals. Full list and DOIs on ORCID.
I'm a Ph.D. candidate in physics & astronomy at Sun Yat-sen University (Zhuhai), advised by Prof. Chengyuan Li and affiliated with the CSST Science Center for the Greater Bay Area — and an AI application engineer at heart.
On the science side I model accretion onto magnetised white dwarfs from wide-field surveys (DESI, LAMOST, ZTF), and I care a great deal about methods that report their own uncertainties honestly — including the ambiguous cases. On the engineering side I build retrieval-augmented, tool-calling LLM agents that carry the same evidence-first discipline into software.
Direct master–Ph.D. track, School of Physics & Astronomy; advisor Prof. Chengyuan Li. Top of cohort. Honours: Caixia Zhanchi fellowship, university first-class scholarship, outstanding student leader, first prize in the physics-experiment design competition.
GPA 3.95/5.0 (recommended for postgraduate admission). Calculus 100/100, Linear Algebra 99/100. First-class scholarship; first prize in the physics-experiment design competition.
Solo full-stack build. Domain knowledge base joining vector search, keyword search and knowledge-graph queries; a retrieval → reasoning → tool-calling → result-verification agent loop; domain adaptation of Qwen3 / Llama-7B with LoRA + Attn-Res. Helped produce a published paper.
Built a knowledge-graph pipeline (literature parsing → chunk cleaning → LLM entity/relation extraction → entity disambiguation → Neo4j triples), a dual-track RAG store (vector + multi-hop graph), and a LangGraph ReAct Q&A agent that verifies intermediate results before generating and checking materials-experiment plans.
Designed a “deterministic engine + LLM agent” hybrid for large-scale advertising-creative analysis — pushing the analytics logic into a trusted compute backend to eliminate quantitative hallucination, with self-contained interactive ECharts reports. Cut monthly review effort by 95%+.
A “differential evolution + L-BFGS-B + DBSCAN” framework for white-dwarf magnetic-field inversion; Numba- and multicore-accelerated spectral fitting for large samples; and SVM/CNN classification over ~10⁸ ZTF light curves to surface rare-object candidates. Served a national survey programme.
Fine-tuning & local deployment (Ollama, vLLM); LoRA + Attn-Res domain adaptation and alignment; dual-track RAG (vector + keyword); LangGraph agents, ReAct loops and knowledge graphs (Neo4j); prompt engineering.
Python (expert) and PyTorch; MATLAB, Mathematica; parallel computing (Numba, multiprocessing); FastAPI services; Git / GitHub.
Global optimization (differential evolution, L-BFGS-B); CNN, SVM, DBSCAN, BP networks; preprocessing & feature engineering for massive time-series and high-dimensional spectra.