{"id":"268601de-21bd-4704-ba03-060482384e73","arxiv_id":"2606.19854","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"M dwarf faculae appear dark because of shallower flux tubes and reduced vertical temperature gradients, as modeled with MURaM and MPS-ATLAS.","lead":"Simulations indicate that faculae on M dwarf stars appear dark rather than bright, unlike on the Sun, due to shallower magnetic flux tubes and weaker temperature gradients. This affects how stellar activity contaminates exoplanet transmission spectra observed by JWST.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Accuracy of MURaM/MPS-ATLAS in producing shallower flux tubes and reduced temperature gradients on M dwarfs without modeling bias","rationale":"The reader's weakest assumption is precisely the load-bearing step required for the qualitative explanation to be reliable. Because the manuscript presents only a qualitative interpretation derived from those specific codes, any unquantified bias in geometry or temperature structure directly undermines the claimed cause. No stronger internal inconsistency is visible from the abstract-level description.","tokens_in":1718,"tokens_out":314,"duration_ms":12313,"concrete_test":"Recompute the facular intensity contrast and flux-tube depth diagnostics from the existing MURaM snapshots using an independent radiative-transfer code (e.g., RH or MULTI3D) with the same 3D cubes but updated M-dwarf opacity tables; if the sign of the contrast or the relative tube shallowness reverses, the attribution does not survive the change in post-processing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a qualitative attribution of the dark-to-bright facular transition to shallower flux tubes and weaker vertical temperature gradients on M dwarfs. This rests entirely on the 3D MHD runs and subsequent radiative transfer correctly capturing the geometry and thermal structure; any systematic offset in how MURaM sets the Wilson depression depth, convective overshoot, or how MPS-ATLAS treats the wavelength-dependent opacity and line formation at Teff ~ 3000–3500 K would make the reported cause an artifact rather than a physical explanation.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper uses 3D MURaM radiative MHD simulations of faculae on M dwarfs (combined with MPS-ATLAS radiative transfer) to provide a qualitative explanation for the observed transition from bright faculae on G/K dwarfs to dark faculae on M dwarfs. The attribution is to shallower magnetic flux tubes and reduced vertical temperature gradients at M-dwarf surfaces relative to the Sun, with implications for stellar contamination in exoplanet transmission spectra.","tokens_in":1824,"tokens_out":492,"duration_ms":28423,"significance":"If the attribution holds, the result would be significant for modeling stellar activity signals in JWST-era exoplanet observations, where M-dwarf hosts are common. The adoption of 3D RMHD simulations is a methodological strength that allows direct inspection of flux-tube geometry and thermal structure, going beyond 1D models.","major_comments":[{"comment":"Abstract and results section: The central claim is presented purely as a qualitative interpretation of simulation outputs with no reported quantitative metrics (e.g., wavelength-dependent intensity contrasts, Wilson depression depths, or vertical temperature gradient values with uncertainties). This is load-bearing because the attribution to 'shallower flux tubes and reduced vertical temperature gradients' cannot be assessed for robustness or compared to observations without such numbers.","section":"Abstract"},{"comment":"Methods (MURaM and MPS-ATLAS description): The paper does not address or quantify potential systematic biases in the codes for Teff ~ 3000–3500 K regimes, such as treatment of wavelength-dependent opacities, convective overshoot, or line formation that could alter the reported flux-tube geometry and temperature structure. This directly affects whether the shallower tubes are a physical result or a modeling artifact.","section":"Methods"}],"minor_comments":[{"comment":"The abstract states the explanation is 'qualitative' but does not clarify what specific simulation diagnostics (e.g., maps of temperature or magnetic field at tau=1) were used to identify the shallower tubes.","section":"Abstract"},{"comment":"No direct comparison is shown between the simulated facular contrasts and existing observational constraints on M-dwarf faculae (e.g., from Kepler or TESS).","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their constructive comments and positive assessment of the work's significance. We address each major comment below, indicating where revisions will be made to strengthen the manuscript.","responses":[{"response":"We agree that quantitative metrics would strengthen the attribution and enable direct comparison with observations. Although the manuscript is framed as providing a qualitative explanation, we will revise the results section (and update the abstract accordingly) to report specific values extracted from the MURaM snapshots and MPS-ATLAS calculations. These will include Wilson depression depths, vertical temperature gradients at selected optical depths, and sample wavelength-dependent intensity contrasts, each accompanied by notes on their derivation and any associated uncertainties. This addition will support the interpretation without changing the overall qualitative focus.","revision_made":"yes","referee_comment":"[Abstract] Abstract and results section: The central claim is presented purely as a qualitative interpretation of simulation outputs with no reported quantitative metrics (e.g., wavelength-dependent intensity contrasts, Wilson depression depths, or vertical temperature gradient values with uncertainties). This is load-bearing because the attribution to 'shallower flux tubes and reduced vertical temperature gradients' cannot be assessed for robustness or compared to observations without such numbers."},{"response":"The referee correctly identifies a gap in the current Methods description. We will add a new paragraph discussing the relevant aspects of MURaM and MPS-ATLAS for the Teff ~3000–3500 K range, including opacity handling, convective overshoot, and line formation. We will reference existing validation studies for these codes in cool-star regimes and note that the shallower flux-tube geometry is a consistent outcome across our simulation set. However, a full, quantitative assessment of all possible systematic biases would require an extensive additional parameter study that lies outside the scope of this paper; the revision will therefore be limited to an explicit discussion of these issues and their potential impact.","revision_made":"partial","referee_comment":"[Methods] Methods (MURaM and MPS-ATLAS description): The paper does not address or quantify potential systematic biases in the codes for Teff ~ 3000–3500 K regimes, such as treatment of wavelength-dependent opacities, convective overshoot, or line formation that could alter the reported flux-tube geometry and temperature structure. This directly affects whether the shallower tubes are a physical result or a modeling artifact."}],"tokens_in":1364,"tokens_out":503,"duration_ms":29229,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point here is that MURaM 3D MHD runs plus MPS-ATLAS radiative transfer produce dark faculae on M dwarfs by making the flux tubes shallower and flattening the vertical temperature gradients compared to the Sun. That attribution is framed as the new result.\n\nThe work does a reasonable job framing the issue for exoplanet transmission spectra and choosing 3D simulations to capture the magnetic structures. Running the codes across spectral types to show the bright-to-dark shift is a straightforward way to test the idea.\n\nThe soft spots are the lack of any quantitative output in the abstract—no contrast values, no wavelength curves, no error bars, and no comparison to actual M-dwarf observations. The whole explanation rests on the simulations getting the Wilson depression and temperature structure right at 3000–3500 K, where opacities and convection are harder to model. If those inputs have systematic offsets, the reported cause becomes an artifact. The stress-test concern about modeling bias holds up on the given material.\n\nThis is aimed at groups correcting stellar contamination in JWST M-dwarf planet data. A reader already working on that would pick up the simulation approach, but the paper needs the missing metrics and validation to be useful. It deserves peer review so referees can check whether the full text adds those pieces and whether the literature overlap is properly handled.","headline":"Simulations suggest shallower flux tubes explain dark faculae on M dwarfs, but the support stays qualitative with no numbers or direct data checks.","tokens_in":2358,"tokens_out":345,"would_cite":false,"duration_ms":26480,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"M dwarf faculae appear dark because their magnetic flux tubes are shallower and have weaker vertical temperature gradients than on the Sun.","keywords":["M dwarfs","faculae","stellar magnetic activity","exoplanet transmission spectra","radiative MHD simulations","flux tubes","temperature gradients","stellar variability"],"falsifier":"Resolved spectroscopy or imaging of an M dwarf that measures the actual vertical extent of a facular magnetic field concentration or the temperature gradient inside it.","tokens_in":2642,"feed_emoji":"🌟","tokens_out":651,"duration_ms":12980,"temperature":0.7,"pith_summary":"The paper examines the unexpected darkening of faculae on M dwarfs, in contrast to their brightening effect on the Sun and G/K dwarfs. This matters for exoplanet studies because magnetic features on the host star alter transit depths in a wavelength-dependent manner and contaminate atmospheric transmission spectra. The authors use 3D radiative MHD simulations to model facular regions and radiative transfer calculations to produce spectra, then attribute the darkening to structural differences at M dwarf surfaces. A sympathetic reader would conclude that the net brightness contribution from small-scale magnetic fields reverses sign for cooler stars. This reversal changes how stellar activity must be subtracted when analyzing planetary signals.","feed_headline":"M dwarf faculae darken from shallower magnetic flux tubes","feed_subtitle":"Reduced temperature gradients flip the usual brightening, requiring revised corrections for exoplanet atmosphere signals from M-dwarf hosts.","key_machinery":"shallower magnetic flux tubes combined with reduced vertical temperature gradients at M dwarf surfaces","core_discovery":"The central claim is that faculae transition from bright to dark across the G to M spectral sequence because magnetic flux tubes are shallower on M dwarfs and the vertical temperature gradients near the surface are reduced relative to the Sun. The simulations demonstrate that these geometric and thermal differences cause the facular regions to appear darker than the surrounding quiet photosphere at visible wavelengths, reversing the conventional brightening seen on hotter stars.","pith_inferences":["Similar darkening could appear in even cooler objects such as L dwarfs if the same structural trends continue.","Direct measurements of facular magnetic field depth on an M dwarf would provide an independent test of the proposed mechanism.","The result suggests that assumptions about facular brightening in stellar atmosphere codes need re-examination for stars cooler than about 4000 K."],"forward_implications":["The contribution of faculae to stellar disk-integrated brightness becomes negative on M dwarfs rather than positive.","Wavelength-dependent corrections for stellar contamination in exoplanet transit observations must include darkening from faculae on M-dwarf hosts.","The contrast sign between spots and faculae reverses at a specific point in the stellar temperature sequence.","Models of stellar variability for late-type stars require separate treatment of facular contrast based on spectral type."],"fun_headline_variants":["M-dwarf faculae darken via shallower flux tubes","Reduced gradients darken faculae on M dwarfs","Shallow flux tubes flip faculae dark on M stars","Lower temperature gradients cause dark M-dwarf faculae","Magnetic shallowness darkens faculae in M dwarfs"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The 3D radiative MHD simulations accurately capture the geometry of magnetic flux tubes and the temperature structure without major unaccounted biases.","fun_headline_variants_meta":{"raw":{"variants":["M-dwarf faculae darken via shallower flux tubes","Reduced gradients darken faculae on M dwarfs","Shallow flux tubes flip faculae dark on M stars","Lower temperature gradients cause dark M-dwarf faculae","Magnetic shallowness darkens faculae in M dwarfs"]},"model":"grok-4.3","cost_usd":0.005652,"raw_usage":{"total_tokens":2620,"prompt_tokens":665,"num_sources_used":0,"completion_tokens":73,"cost_in_usd_ticks":56515500,"prompt_tokens_details":{"text_tokens":665,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":1882,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":665,"tokens_out":73,"duration_ms":19443,"temperature":1.0,"reasoning_tokens":1882,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T16:00:53.455888+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Resolved spectroscopy or imaging of an M dwarf that measures the actual vertical extent of a facular magnetic field concentration or the temperature gradient inside it.","supporting_citations":[],"review_version":1}