1. How do massive stars assemble against their own feedback?

    Massive stars shape galaxies even as they form — launching radiation, winds, and collimated outflows that push back on the gas feeding them. I build radiation-magnetohydrodynamic models to ask how a star keeps gaining mass despite forces strong enough to reverse its own accretion.

    • Radiation pressure
    • Protostellar outflows
    • Magnetic fields
    • Accretion
  2. How does stellar feedback sculpt clusters and the interstellar medium?

    The winds, radiation, and supernovae of young massive clusters carve bubbles, stir turbulence, and regulate further star formation. I connect simulations to multi-wavelength observations — from X-rays to the far-infrared — to trace where that feedback energy actually goes.

    • Stellar winds
    • H II regions & superbubbles
    • Turbulence
    • 30 Doradus & the LMC
  3. What can observations genuinely tell us about stellar populations?

    The initial mass function, binarity, and cluster structure are all inferred from incomplete data. I've shown that ignoring unresolved binaries can make high-mass IMF conclusions confidently wrong as samples grow. I want to know what surveys like Rubin and Gaia can — and cannot — actually recover.

    • The initial mass function
    • Unresolved binaries
    • Identifiability
    • Rubin/LSST & Gaia
    Three panels showing binaries raising the line-of-sight velocity dispersion, binary-blind mass estimates near three times the truth against a calibrated control and binary-aware recovery, and radial target allocations shifting between mass-first and degeneracy-first observing designs.Preliminary · in prep
    Unresolved binaries inflate the observed velocity dispersion, and a binary-blind fit turns that orbital motion into 2.9 times the true cluster mass — while forecasting 4.5% precision. That is a bias 42 times its own error bar. Modelling the binaries recovers both the mass and the binary fraction, and the same machinery then says where to point the telescope next.A. L. Rosen · Jaxstro (informax). Preliminary.
  4. How can differentiable models transform astrophysical inference?

    If every step from a system's birth to its observed light is a differentiable function, we can compute exactly how each observable responds to the underlying physics — and design the observations most likely to break a degeneracy. This is the idea behind the Jaxstro ecosystem, and it reaches well beyond any one problem.

    • Automatic differentiation
    • Forward modeling
    • Optimal experimental design
    • Scientific software
    Three panels. A 512-star cluster in the x-y plane with a central massive binary; cluster and binary scales versus time showing the cluster expanding while the binary semimajor axis contracts in discrete steps; and a matrix of logarithmic sensitivities of endpoint observables to birth parameters, with cluster radii responding to the initial half-mass radius and the binary responding to its own initial separation and eccentricity.Preliminary · in prep
    A differentiable model has to earn the word. Here a 512-star cluster with a hard massive binary is evolved for three crossing times — the binary hardens, the cluster expands — and the endpoint's sensitivity to every birth parameter is obtained by differentiating straight through the adaptive timestep, agreeing with complete finite-difference reruns to three parts in a hundred thousand. The result is legible: cluster size remembers where the cluster started, the binary remembers its own orbit. That is what makes the model invertible against data rather than merely runnable.A. L. Rosen · Jaxstro (gravax). Preliminary.