{"id":"784a9d98-f1c0-4d60-b6af-d42558dff2d4","arxiv_id":"2607.17144","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"First numerical ray-shooting simulation of compound (two-plane) quasar microlensing in the double-source-plane lens J1721+8842 shows new caustic morphologies — convexity violations and lip caustics — that do not occur in single-plane microlensing.","lead":"This paper simulates, for the first time, what happens when quasar light is microlensed by stars in two separate foreground galaxies instead of one — 'microlensing of microlensing.' It shows the predicted extra layer of caustic features and brightness peaks appears in realistic maps of the recently discovered Einstein zig-zag lens J1721+8842, which matters for precision Hubble-constant measurements with such systems.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ray-shooting applies two single-plane lens equations sequentially, omitting the η factor in the multi-plane lens equation; smooth-limit maps would give μ_lens×μ_s1 instead of μ_eff, so the compound maps are not valid two-plane predictions.","rationale":"The Reader's weakest_assumption focuses on the κ_s = 0 matter model and the simplistic macromodel. Those are important physical limitations, but they are explicitly acknowledged by the authors and do not by themselves invalidate the demonstration under the assumed (maximal-stellar-density) model. The concern I identify is more fundamental: the ray-shooting algorithm, as described, appears to implement two independent single-plane lens equations in sequence, which is not equivalent to the multi-plane lens equation the paper itself derives in Equation (5). The missing η factor means the maps are not correct two-plane lensing predictions for J1721. This directly undermines the central claim of 'first recursive two-plane ray-shooting' and the practical demonstration of compound microlensing in this system. If the code is indeed sequential, the paper requires a major correction—either re-implementing the ray-shooting with the correct multi-plane equation or repositioning the work as an approximate study. The smooth-limit test is decisive and easy to run; the product-vs-effective magnification discrepancy provides a crisp numerical check. I therefore recommend REJECT for the current version, though the scientific idea may be salvageable with a corrected implementation.","tokens_in":14579,"tokens_out":49606,"duration_ms":416222,"concrete_test":"Re-run the compound ray-shooting for image A with the point-mass microlenses switched off (set all m_i = 0) while retaining the smooth M matrices from Tables 2 and 3. Measure the resulting large-scale magnification of the map. If it is approximately μ_lens × μ_s1 = −16.6 rather than the macro-model value μ_eff = −14.5 (Table 1), the code is implementing two sequential single-plane lens equations instead of the correct multi-plane equation (5). A complementary test is to compare the maps directly against a reference implementation that evaluates β = θ − η α_lens(θ) − α_s1(θ − α_lens(θ)) for the same microlens fields; any difference in the caustic networks would confirm the omission.","verdict_should_be":"REJECT","load_bearing_attack":"The central claim—that these are the first recursive two-plane ray-shooting simulations of compound microlensing—requires the code to solve the correct multi-plane lens equation (5): β = θ_lens − η α_lens(θ_lens) − α_s1(θ_s1), with θ_s1 = θ_lens − α_lens(θ_lens). However, the code description (§3.2) states that rays are 'sequentially lensed at multiple planes' and that the compound map 'takes rays generated on the lens plane, maps them to the s1 plane, and then further maps these onto the s2 plane.' This is implemented by applying Equation (8) twice: first with the foreground-lens parameters from Table 2 (which are defined for source s1, not s2), then with the s1-plane parameters from Table 3. This yields β = M_s1(M_lens θ − Σ_lens) − Σ_s1, which is missing the η factor that multiplies α_lens in the exact multi-plane equation. The omission is not a minor detail: in the smooth-only limit, the sequential composition gives magnification det(M_s1 M_lens)^{-1} = μ_lens × μ_s1, whereas the correct effective magnification from Table 1 is μ_eff. For image A, μ_lens × μ_s1 = (−6.178)(2.687) = −16.6, while μ_eff = −14.5, a ~14% discrepancy. Thus the compound maps and lightcurves do not correspond to the stated DSPL geometry of J1721+8842. The fact that η ≈ 1.02 for this system makes the error numerically small, but the method is not a valid multi-plane ray-shooting scheme and can fail more severely for other DSPL configurations where η differs from unity.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the first claimed numerical simulations of compound microlensing in the double-source-plane lens J1721+8842. The authors fit a macromodel of the system, extract convergence and shear at the six foreground-lens image positions and at two intermediate s1-plane image clusters, and then ray-shoot through two successive planes populated by 1 M⊙ point masses using a GPU-accelerated JAX code. They produce six compound magnification maps, lightcurves for point-like and extended quasar sources, and identify caustic morphologies—convexity-violating 'second-kind' cusps and a lip caustic—that are forbidden in single-plane lensing. The central claim is that these features are not theoretical curiosities but persist in realistic high-density microlensing environments, with observable signatures in the lightcurves.","tokens_in":14996,"tokens_out":7874,"duration_ms":87493,"significance":"If the simulations are faithful to the multi-plane lens equation, this is a genuinely new demonstration of a predicted effect. The paper is timely given the recent discovery of the Einstein zig-zag lens and the growing catalogue of DSPLs, and it opens a numerical pathway for studying compound microlensing. Strengths include the first bespoke recursive two-plane ray-shooting implementation, the use of a concrete astrophysical system with a fitted macromodel, and the explicit comparison of single-plane and compound maps and lightcurves. The authors are also candid about several limitations. However, the central result is currently compromised by an apparent mismatch between the stated DSPL lens equation and the implemented sequential two-plane ray-shooting scheme, and by the assumption that all convergence is in compact point masses. The circularity concern raised in the stress-test note does not, in my reading, land: Petters & Wicklin's classification is used to interpret the maps, not to generate them.","major_comments":[{"comment":"The code description states that rays are sequentially lensed: first by the foreground lens using the Table 2 parameters, then by the s1 plane using the Table 3 parameters. This implements a composition of two single-plane lens equations, β = M_s1(M_lens θ − Σ_lens) − Σ_s1. The exact DSPL equation (5), however, contains η α_lens(θ_lens) as the foreground contribution to the s2 source plane. No η factor appears anywhere in §3.2. In the smooth-only limit the sequential composition predicts μ = μ_lens μ_s1; for image A this gives (−6.178)(2.687) ≈ −16.6, whereas Table 1 lists μ_eff = −14.5, a ~14% discrepancy. Thus the maps and lightcurves as described do not solve the stated multi-plane lens equation. If η is in fact incorporated inside the code through a rescaling of Table 2 parameters, that needs to be stated explicitly; otherwise the 'first recursive ray-shooting' claim is not supportab","section":"§3.2, Eq. (5)"},{"comment":"All convergence on both planes is assigned to static 1 M⊙ point masses (κ* = κ, κs = 0). The conclusion that the new caustic morphologies 'persist in realistic, high-density microlensing environments' is stronger than the simulations support: a high-density stellar field with zero smooth matter is an extreme end of the parameter space, not a demonstration of persistence under realistic smooth+compact decompositions. The authors themselves note in §5 that 'It is currently unclear how much convexity violation will be suppressed under a more realistic assumption that the convergence is formed of a smooth distribution of matter as well as compact masses.' This is a load-bearing limitation for the headline claim. I ask the authors to run at least a few maps with a nonzero smooth component (e.g., κs/κ = 0.5 and 0.8) and quantify how the frequency of convexity-violating cusps and lip caustics c","section":"§3.2, §5"},{"comment":"The s1-plane macro-parameters, which drive half of the compound effect, come from a deliberately simplified model in which the s1 quasar images are described as 'not as well focused' and are evaluated only at two mean positions. Tables 1–3 report single median values without credible intervals from the SVI posterior. Since the compound maps depend directly on these values, it is difficult to judge whether the reported morphologies and lightcurve differences are robust against plausible variations in the macromodel. At minimum, the authors should report some measure of uncertainty in κ, γ, and μ for the s1-plane parameters and test whether the second-kind cusps survive within that range.","section":"Appendix A, Tables 1–3"}],"minor_comments":[{"comment":"The identification of convexity-violating cusps and the lip caustic is currently visual. A small automated criterion (e.g., computing the curvature of the two fold segments on either side of each cusp and checking the sign relative to the caustic interior) would make the central morphological claim more quantitative and less vulnerable to selection effects.","section":"§4, Figure 7"},{"comment":"The y-axis labels on the lightcurve panels appear to be magnification values but are unlabelled. Please add axis labels or clarify the units in the captions.","section":"Figures 4–5"},{"comment":"The code and data are described as available 'on reasonable request'. Given that the paper's novelty is a new numerical pipeline, a public release of the code would materially help reproducibility and future use.","section":"Data Availability"},{"comment":"The relationship between η and the alternative β parameter is explained in the footnote, but it may be helpful to state explicitly that Tables 2 and 3 are defined with s1 as the source and therefore that the η factor in Eq. (5) is required when these tables are used for s2 maps.","section":"§2.2.1"}],"recommendation":"major_revision","confidential_remarks":"The missing η factor is the key issue. If the authors can show that the code already rescales the foreground-lens parameters by η (e.g., inside the construction of the matrix M), then the paper may be close to acceptance after a clarification. If not, the simulation pipeline and all downstream results need to be redone. The κs=0 assumption is also a substantive overreach relative to the 'realistic environment' language, though it is fixable with additional runs. I would not reject outright because the underlying compound-microlensing idea is sound and the qualitative effect is likely robust, but the manuscript in its current form does not yet demonstrate what it claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what's new: the paper builds the first recursive two-plane microlensing ray-shooting code and produces the first compound microlensing magnification maps and lightcurves for a real double-source-plane lens (J1721+8842). The formalism in §2 is standard and clearly laid out, and the authors are upfront about their simplifications—κ_s=0, a simplistic macromodel with poorly focused s1 images, and no propagated uncertainties. That honesty is real, and the qualitative observation that a second plane of compact masses adds a new caustic network and extra peaks in the lightcurves is credible.\n\nThe problem is in the ray-shooting implementation. The code applies Equation (8) twice: first from the foreground lens to the s1 plane, then from s1 to s2. That sequential composition gives β = θ − α_lens(θ) − α_s1(θ − α_lens(θ)), whereas the correct multi-plane lens equation is θ_s2 = θ − η α_lens(θ) − α_s1(θ − α_lens(θ)). The η factor is missing. The stress-test note is right: in the smooth-limit the sequential maps give μ_lens × μ_s1, not μ_eff. For image A that's a 14% discrepancy in magnification. For J1721 η ≈ 1.02, so the error is numerically small, but the method is not a valid multi-plane ray-shooting scheme, and the maps do not describe the stated geometry. This is a load-bearing flaw for the quantitative claims, though the qualitative phenomenon—two planes of deflectors produce richer caustic structure—almost certainly survives, since the sequential composition is itself a two-lens system with η=1.\n\nThe other soft spots are secondary. The κ_s=0 assumption is acknowledged and could suppress the convexity violations; the macromodel uncertainties are not propagated into the maps. But the η issue is the one a referee should demand be fixed before any quantitative use.\n\nOverall: worth a serious referee, but not as is. The concept is timely and the authors are thinking clearly. I'd send it to review with a request to correct the multi-plane equation, rerun the maps, and test a smooth dark-matter component.","headline":"First compound-microlensing simulation for a real DSPL, but the ray shooting appears to drop the η factor in the multi-plane lens equation, so the maps are for η=1 rather than the claimed J1721 geometry.","tokens_in":15566,"tokens_out":11266,"would_cite":false,"duration_ms":107743,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Sequential microlensing by two lens planes produces caustic shapes that a single plane forbids, as shown by first recursive ray-shooting simulations of the Einstein zig-zag lens J1721+8842.","keywords":["compound microlensing","double-source-plane lens","Einstein zig-zag","J1721+8842","caustics","ray shooting","quasar microlensing","gravitational lensing"],"falsifier":"Compute the same six J1721 magnification maps with, say, half of the convergence placed in a smooth component (κ_s/κ = 0.5) while keeping the macro parameters fixed; if convexity-violating cusps and lip caustics become too rare to distinguish from single-plane maps, the claim that they persist in realistic high-density microlensing environments would be refuted.","tokens_in":14385,"feed_emoji":"🔭","tokens_out":8259,"duration_ms":71134,"temperature":0.7,"pith_summary":"The paper argues that when quasar light is deflected by compact stars on two successive lens planes, the combined microlensing produces magnification maps with caustic morphologies that are strictly forbidden in single-plane lensing: convexity-violating 'second-kind' cusps and lip caustics. Using the recently discovered double-source-plane system J1721+8842 as a fiducial case, the authors build the first recursive ray-shooting code that traces rays through two planes of point masses, and show these structures persist in realistic high-density star fields. The simulated lightcurves correspondingly show both an overall amplification of peaks from the first lens plane and an added population of small-scale peaks from new caustic crossings. This matters because compound lenses are promising cosmological probes, and microlensing noise must be understood if they are to deliver percent-level Hubble constant measurements.","feed_headline":"Two lens planes create caustics a single lens cannot","feed_subtitle":"First recursive ray-shooting maps of the Einstein zig-zag lens J1721+8842 show extra caustics and new peaks.","key_machinery":"The load-bearing device is a recursive, GPU-accelerated ray-shooting pipeline that traces rays through two successive planes of point-mass microlenses. The lens equation (a one-to-many mapping) is applied twice: first from the foreground lens to the intermediate source plane, then from that plane to the background quasar source plane, with each ray's final position determined by summing deflections from all microlenses on each plane. The code tiles and vectorises the computation to manage the large number of lenses and rays, and it produces 10^4 × 10^4 pixel magnification maps for each of the six images. The singularity-theory classification of multi-plane lensing is what identifies the extr","core_discovery":"The central claim is that compound microlensing—sequential deflection of quasar light by compact masses on two lens planes—produces microcaustic networks that are qualitatively richer than single-plane microlensing. Building on the singularity-theory classification of multi-plane lensing, the authors demonstrate by simulation that the double-plane maps contain caustics of the 'second kind', where one fold at a cusp is concave relative to the inside of the caustic, as well as lip caustics, both of which are absent from single-plane maps because a single lens mapping cannot produce such convexity violations. These features appear in all six simulated images of the Einstein zig-zag lens J1721+8","pith_inferences":["The paper sets every bit of surface density in both planes to point masses, so the signature caustics may be weaker in real galaxies; if a smooth dark-matter component is added, the convexity violations could disappear or become too rare to detect.","The two 'zig-zag' images D and F showed more convexity-violating cusps in the simulations; if this holds over a wider parameter range, zig-zag images could serve as a targeted diagnostic for compound microlensing.","Typical double-source-plane systems have late-type star-forming galaxies on the intermediate plane; their stellar mass functions differ from J1721's early-type deflector, so the visibility of compound signatures may vary across the population.","The recursive ray-shooting approach could be extended to more than two lens planes; if so, the most complex caustic networks may appear in the rare systems with multiple background sources."],"forward_implications":["If compound microlensing produces these additional caustic structures, then flux ratios and lightcurves of double-source-plane quasars encode information about the mass distribution on both lens planes, not just one.","Lightcurves of the six J1721 images should show both enhanced peak magnifications and an extra population of short-duration peaks relative to a single-plane model; this is directly testable with monitoring data.","Compound microlensing adds a new source of brightness noise for double-source-plane time-delay cosmography, so statistical characterization of this effect is needed before percent-level Hubble constant measurements are attempted.","The persistence of the distinctive caustics in dense star fields validates the use of static point-mass fields for such simulations and motivates extending them to moving microlenses and realistic mass functions."],"fun_headline_variants":["Double-plane microlensing adds caustics single-plane can't make","Caustics of the second kind emerge in double-plane microlensing","Two-plane lensing spawns caustic networks absent in single-plane","Compound microlensing shows new caustics in Einstein zig-zag"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper assumes that all of the convergence on both lens planes is made of static solar-mass stars (no smooth dark matter), so the distinctive caustic shapes could be suppressed in a realistic mix of stars and dark matter—a possibility the authors explicitly flag.","fun_headline_variants_meta":{"raw":{"variants":["Double-plane microlensing adds caustics single-plane can't make","Caustics of the second kind emerge in double-plane microlensing","Two-plane lensing spawns caustic networks absent in single-plane","Compound microlensing shows new caustics in Einstein zig-zag"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001286,"raw_usage":{"total_tokens":5079,"prompt_tokens":722,"completion_tokens":4357,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":4290}},"tokens_in":466,"tokens_out":4357,"duration_ms":28814,"temperature":1.0,"reasoning_tokens":4290,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T18:53:41.447384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same six J1721 magnification maps with, say, half of the convergence placed in a smooth component (κ_s/κ = 0.5) while keeping the macro parameters fixed; if convexity-violating cusps and lip caustics become too rare to distinguish from single-plane maps, the claim that they persist in realistic high-density microlensing environments would be refuted.","supporting_citations":[],"review_version":1}