Astrophysicist · AI / Agent Engineer

Jiamao Lin林佳茂

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.

  • 10MG — 240 MG fields measured
  • 4DESI/LAMOST polars modelled
  • 6first / co-author papers
  • 10⁸ ZTF light curves searched

Ph.D. candidate, Physics & Astronomy, Sun Yat-sen University (Zhuhai) · advised by Prof. Chengyuan Li · CSST Science Center for the Greater Bay Area.

What I work on

How do you weigh the magnetic field of a star you cannot resolve? You read it off the light.

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.

Cyclotron spectroscopy

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.

Automated spectral fitting

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.

Time-domain follow-up

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.

Machine learning on light curves

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.

Disk-eclipsing binaries & clusters

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.

AI agents for astronomy

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 →

Cyclotron light-curve modelling

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.

Fig. 1 A synchronous polar, solved in the co-rotating frame drag the scene or the curve to scrub φ
  • White dwarf
  • Roche-lobe donor
  • Ballistic stream
  • Dipole field
  • Cyclotron beam

Phase-folded cyclotron light curve

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.

Beaming

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°.

The donor

Φ = −µ₁/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.

The eclipse

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 white dwarf is drawn ~30× oversized so the pole geometry is legible; the field-line shells, the threading radius Rµ ≈ 0.17 a and the curtain width are schematic. Everything else — lobe, L1, stream, dipole geometry, projection, beaming and occultation — is computed. Illustrative of the AM Her mechanism, not a fit to one system.

Cyclotron spectrum lab

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.

Fig. 2 Three components, one emergent spectrum hover for λ, F and τn
Harmonic ladder

λ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.

Radiative transfer

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.

Why B is ambiguous

τ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.

Leading-order Chanmugam & Dulk (1981) harmonic scaling through a homogeneous isothermal slab, with Planck photosphere and hot-spot components; hydrogen lines and full angle-dependent transfer are omitted. Built to make the hump spacing and the optical-depth turnover legible, not to reproduce a published fit.

Square-wave eclipse: a disk-eclipsing binary

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.

Fig. 3 A knife edge crossing a star, against real ZTF photometry drag to scrub · inset is ×4
  • Occulted star
  • Companion (visible)
  • Circumbinary disk
  • Occulting edge

Folded light curve · real ZTF data

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.

The ingress

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.

Where the disk sits

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.

Why the floor is 60%

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.

Light curve = real ZTF g/r data folded at P = 36.710719 d with a χ²-fitted trapezoid (~40% depth). The scene is a knife-edge (zero-thickness) model driven by two angles; stellar radii are exaggerated ~15× so the crossing is visible, which is what the magnified inset compensates for.

Binary evolution on a cluster clock

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.

Fig. 4 A cluster clock separates three evolutionary outcomes select a stage, then play
NS Per · a surviving accreting binaryOrbital-plane view
Model chronology; stage durations and orbital motion are compressed.
    One geometry for the flow

    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.

    What is physical: Roche equipotentials and the pressure-free ballistic trajectory are calculated in the circular, co-rotating point-mass approximation. NS Per's chronology and orbital period come from the revised COSMIC production track. Its αCE = 0.20 is the representative model choice, not the posterior centre. What is schematic: common-envelope morphology, disk size, and unresolved stellar markers. Orbital speeds and stage durations are compressed. The orbital-plane view does not reproduce the measured viewing inclination. During a common envelope, equilibrium Roche surfaces are hidden.
    Model provenance and observational limits

    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.

    Interactive figure 05 · General relativity

    When a black hole gives energy back.

    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.

    Lin et al. · Physics Letters B, 2021 ↗
    Light around a black holeSchwarzschild · a = 0
    Light bending · thin disk

    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.

    Near-horizon wave · schematica/M = 1 · Q = 0
    Local scalar waveEnergy transfer from the black hole

    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.

    Covered by the stability proof

    The superradiant and bound-state conditions hold, but the paper excludes an exterior trapping well in this parameter range.

    01 / ENERGY EXTRACTION

    Can the horizon supply energy?

    0 < Mω < m/2Yes · 0.30 < 0.50

    The threshold fixes the sign of the horizon energy flux. It does not give an amplification percentage.

    02 / CONFINEMENT

    Can the field decay at infinity?

    Mω < MμYes · 0.30 < 0.80

    This is required for a massive bound state. By itself, it does not establish a trapping well or a growing mode.

    03 / THE PAPER'S RESULT

    Is a trapping well excluded?

    Mω < Mμ / √3Yes · sufficient condition met

    Within the superradiant window, this sufficient condition rules out the feedback needed for this instability.

    The stability map

    Click or drag to choose parameters
    Covered by the proofNot settled by this proofOutside joint window

    Read the boundary

    Amplification does not guarantee an instability.

    The green region meets the paper's sufficient stability condition within the superradiant and bound-state window. Moving outside green does not establish instability.

    a = M, Q = 0
    0 < Mω < m/2
    Mω < Mμ / √3

    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.

    Physics, assumptions & visual references

    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.

    1. Lin, Luo, Zheng, Yin & Huang (2021) — original paper: the analytic results implemented here.
    2. William East — superradiance simulations: a reference for showing energy transfer and separating scattering from bound-field growth. Those numerical simulations are not reproduced here.
    3. NASA / Jeremy Schnittman — a black hole's warped world: a reference for clear spatial annotations; accretion-disk lensing is a different phenomenon.
    4. ScienceClic / Alessandro Roussel — Black Hole #4: visual and interaction reference. No code or imagery is reused.
    5. Eric Bruneton — real-time black-hole rendering: a reference for Schwarzschild light bending and browser interaction. Our renderer uses direct numerical integration rather than the precomputed beam-tracing method. No shader code or imagery is reused.

    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.

    Publications

    Selected work, newest first.

    1. 2026
      First & corresponding author

      A Candidate Low-Mass Disk-Eclipsing Binary in the ~316 Myr Open Cluster UPK 13

      J. Lin et al.

      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.

      disk eclipseopen clusterscircumbinary diskZTF
      See Fig. 3 →
    2. 2026
      First & corresponding author

      Constraining Magnetic Fields of Polars through Automated Cyclotron Spectral Fitting

      J. Lin, W. Guo, Y. Sun, C. Li, W. Hou

      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.

      polarscyclotronwhite dwarfsautomated fitting
      See Fig. 2 →
    3. 2025
      First & corresponding author

      Discovery and characterization of ZTF J0112+5827: an 80.9-minute polar with strong cyclotron features

      J. Lin et al.

      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.

      polarcyclotronX-rayDoppler tomography
      See Fig. 1 →
    4. 2025
      Co-author

      A magnetic white dwarf formed through a binary merger within 35 million years

      H. Yan, J. Lin, et al.

      The Astrophysical Journal Letters (accepted)

      A young, strongly magnetic white dwarf whose field and rotation point to a recent double-degenerate merger origin.

      magnetic white dwarfbinary merger
    5. 2025
      Co-author

      Spectroscopic Identification and Population Properties of Cataclysmic Variables from DESI DR1

      W. Hou, J. Lin, et al.

      The Astrophysical Journal Supplement Series (accepted)

      A spectroscopically identified DESI DR1 sample of cataclysmic variables and their population properties.

      cataclysmic variablesDESIspectroscopy
    6. 2021
      First author

      Extremal rotating black holes, scalar perturbation and superradiant stability

      J.-M. Lin, M.-J. Luo, Z.-H. Zheng, L. Yin, J.-H. Huang

      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 →
      black holessuperradianceKerr–Newmangeneral relativity
    7. +

      Additional work

      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.

    About

    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.

    • Invited talks at Tsinghua University, the Chinese Academy of Sciences and the National Astronomical Observatory
    • Independent referee for RAA (Research in Astronomy and Astrophysics)

    Contact

    Education

    2022 — 2027 (expected)

    Ph.D., Physics — Sun Yat-sen University

    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.

    2018 — 2022

    B.Sc., Physics — South China Normal University

    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.

    Experience

    2025.10 — present

    Astro Agent — astronomy knowledge-graph & analysis agent open source

    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.

    2026.02 — 2026.04

    AI Application Engineer — knowledge graph for PKU Materials Science

    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.

    2026.04 — 2026.05

    AI Application Engineer — Blue Moon Group (蓝月亮)

    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%+.

    2023.03 — 2025.10

    Core algorithms — magnetic-field modelling & large-scale time-series mining

    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.

    Skills

    LLM · Agents · RAG

    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.

    Programming & engineering

    Python (expert) and PyTorch; MATLAB, Mathematica; parallel computing (Numba, multiprocessing); FastAPI services; Git / GitHub.

    Algorithms & data

    Global optimization (differential evolution, L-BFGS-B); CNN, SVM, DBSCAN, BP networks; preprocessing & feature engineering for massive time-series and high-dimensional spectra.