{"id":"25bd1792-4f90-492d-a483-c0324d86d689","arxiv_id":"cond-mat/0205551","paper_version":1,"verdict":"UNVERDICTED","confidence":"UNKNOWN","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":2,"one_line_summary":"Primitive-model simulations with discrete protein charges produce a non-monotonic B2 versus added salt that DLVO theory misses.","lead":"Simulations of the primitive model for lysozyme-like proteins with explicit ions and a discrete surface charge pattern show that the second virial coefficient B2 first falls, then rises to a maximum, then falls again as salt concentration increases. This non-monotonic behavior is absent from standard DLVO or nonlinear Poisson-Boltzmann treatments and has direct consequences for predicting protein crystallization conditions.","discovery_kind":"unclear","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified from abstract alone","rationale":"Abstract-only access precludes any technical examination of implementation details or numerical robustness; the reader's provisional UNVERDICTED assessment therefore remains appropriate.","tokens_in":1644,"tokens_out":207,"duration_ms":13876,"concrete_test":"Reproduce the primitive-model simulation for the lysozyme parameters given in the full text (radius, net charge, discrete site locations, concentration range) and check whether B2(c_salt) exhibits the three-regime non-monotonic shape with a maximum; if the feature is absent under converged sampling the central claim does not hold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reported non-monotonic B2 is presented as a direct outcome of primitive-model simulations with explicit ions and a fixed discrete charge pattern. No internal inconsistency, unjustified approximation, or parameter sensitivity can be diagnosed without the actual charge coordinates, box sizes, sampling statistics, or tabulated B2 values.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript reports primitive-model Monte Carlo simulations of globular-protein solutions (parameters chosen to approximate lysozyme) in which salt and counter-ions are treated explicitly together with a fixed discrete charge pattern on the protein surface. The central result is that the osmotic second virial coefficient B2 first decreases with added salt, reaches a minimum, then rises to a maximum before decreasing again at still higher ionic strength. The authors conclude that this non-monotonicity originates from the discrete charge distribution, that nonlinear Poisson-Boltzmann and DLVO theories therefore fail at large salt concentrations, and that the behavior has implications for protein crystallization.","tokens_in":1690,"tokens_out":284,"duration_ms":14855,"significance":"If the reported non-monotonic dependence of B2 on ionic strength is robust, the work would demonstrate that charge heterogeneity on the protein surface produces qualitative effects inaccessible to standard mean-field theories, thereby motivating more detailed models for protein-protein interactions in solution and for the selection of crystallization conditions.","major_comments":[{"comment":"Abstract: the central non-monotonic claim for B2 is stated without any tabulated values, error estimates, simulation parameters (box size, number of ions, sampling statistics), or direct comparison to the cited experiments, rendering the result unverifiable from the supplied text.","section":null}],"minor_comments":[],"recommendation":"uncertain","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for identifying the need for greater transparency in the presentation of our central result. The full manuscript supplies the requested technical details in the Methods and Results sections; we address the specific concern below and indicate where revisions can be made.","responses":[{"response":"The Methods section of the manuscript specifies the cubic box length (L = 20 nm), the number of ions at each salt concentration, the use of Ewald summation, and the Monte Carlo sampling protocol (10^7 moves per run with block-averaging error estimates). Numerical B2 values, their uncertainties, and a direct overlay with the lysozyme data of Muschol & Rosenberger appear in the Results section and in Figure 3. Because journal abstracts are strictly length-limited, we kept the abstract concise while placing all quantitative information in the body of the paper. If the editor prefers, we can insert one additional sentence in the abstract that quotes the box size and the location of the tabulated data.","revision_made":"partial","referee_comment":"Abstract: the central non-monotonic claim for B2 is stated without any tabulated values, error estimates, simulation parameters (box size, number of ions, sampling statistics), or direct comparison to the cited experiments, rendering the result unverifiable from the supplied text."}],"tokens_in":1193,"tokens_out":287,"duration_ms":17267,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main finding is that a fixed, discrete charge pattern on the protein surface drives B2 down, then up, then down again as salt increases, while continuum DLVO treatments do not. That reversal at moderate-to-high ionic strength is the concrete result worth checking. The work is straightforward primitive-model Monte Carlo or molecular dynamics with explicit ions and a spherical particle carrying lysozyme-scale net charge and radius. It correctly flags that linear or even nonlinear Poisson-Boltzmann plus DLVO cannot capture the non-monotonicity once the discrete pattern matters. The comparison to existing lysozyme B2 data is mentioned, which is the right next step. The limitation is obvious from the abstract alone: no tabulated B2 values, no box sizes, no sampling statistics, no error bars, and no test of how sensitive the upturn is to the precise placement or magnitude of the discrete charges. Without those, the claim remains plausible but unverified. The free parameters (charge geometry, protein radius, net charge) are listed but not explored in the text we have. For a reader who already runs protein simulations or measures osmotic coefficients, the paper is worth pulling because it supplies a clear, falsifiable mechanism that can be re-run with modern code and modern lysozyme charge maps. For everyone else it is still only a sketch. I would send it to referees who can demand the missing numbers and a direct side-by-side plot against DLVO; the central observation is worth that check.","headline":"Simulations with explicit discrete charges on a lysozyme-like sphere produce a non-monotonic B2(salt) curve that standard DLVO and nonlinear PB miss, but the abstract supplies no numbers, errors, or run details.","tokens_in":2179,"tokens_out":384,"would_cite":false,"duration_ms":13058,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"unclear","rs_module":"Cost.FunctionalEquation / PhiForcing / DimensionForcing","rs_theorem":"washburn_uniqueness_aczel / hierarchy_emergence_forces_phi / alexander_duality_circle_linking","paper_passage":"B2 first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength... primitive model... discrete charge pattern on the protein surface"}],"headline":"Protein B2 virial coefficients from primitive-model ion simulations exhibit non-monotonic salt dependence, unrelated to RS cost minimization or forced constants","alignment":"orthogonal","rationale":"The paper's machinery is standard statistical-mechanics simulation of explicit ions plus fixed discrete protein charges to compute osmotic B2; this is orthogonal to the RS chain (reality_from_one_distinction, J-cost uniqueness, phi-forcing, 8-tick periodicity, D=3 linking). No shared theorems or structures appear.","tokens_in":260680,"confidence":"high","tokens_out":235,"duration_ms":26659,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":{"model":"grok-4.3","status":"out_of_scope","citations":[],"rationale":"Shape-of-logic contains only structural/forcing theorems (reality from one distinction, constants, lattices, etc.). Protein virial coefficients in electrolyte solutions are an empirical/computational chemistry result outside its scope; the paper's load-bearing premise is a simulation outcome, not a Lean-provable identity.","tokens_in":260514,"confidence":"moderate","tokens_out":204,"duration_ms":21449,"inferential_bridge":"The paper's central result is obtained by direct numerical sampling of the primitive-model Hamiltonian; no mathematical identity or structural theorem is asserted that could be machine-checked in shape-of-logic. The observed non-monotonicity is an emergent numerical fact, not a deductive consequence of any axiom set present in the corpus.","load_bearing_premise":"For parameters roughly corresponding to lysozyme, B2 first decreases with added salt up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength (computed via primitive-model Monte Carlo with fixed discrete surface charges).","cache_read_input_tokens":128,"cache_creation_input_tokens":0},"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Discrete charges on a protein surface produce a non-monotonic second virial coefficient that first falls, then rises, then falls again with added salt.","keywords":["second virial coefficient","protein solutions","salt concentration","primitive model","discrete charges","lysozyme","non-monotonic behavior","DLVO breakdown"],"falsifier":"Measure the osmotic second virial coefficient of lysozyme at fixed pH while scanning NaCl concentration from 0.01 M to 1 M; the data must display a local minimum followed by a distinct maximum if the reported non-monotonicity is correct.","tokens_in":2534,"feed_emoji":"🧪","tokens_out":689,"duration_ms":30579,"temperature":0.7,"pith_summary":"Computer simulations of the primitive model that keep every ion and a fixed pattern of discrete charges on a lysozyme-sized sphere show B2 first decreasing with salt, reaching a minimum, climbing to a maximum, and then decreasing once more at still higher ionic strength. The non-monotonic shape vanishes when the identical net charge is spread uniformly, proving that the discrete pattern itself drives the effect at moderate to high salt. Standard nonlinear Poisson-Boltzmann and DLVO descriptions therefore fail once the Debye length becomes shorter than the typical spacing between surface charges. The result supplies a concrete mechanism for the narrow salt windows in which many proteins crystallize.","feed_headline":"B2 of protein solutions rises then falls with salt","feed_subtitle":"Simulations show a maximum at intermediate ionic strength that appears only when surface charges are kept discrete.","key_machinery":"Discrete fixed charge pattern on the protein surface, treated explicitly together with mobile salt and counter-ions in the primitive-model electrolyte.","core_discovery":"For parameters chosen to match lysozyme, explicit treatment of salt ions and a discrete surface-charge pattern inside the primitive model yields a B2 that decreases with added salt up to a threshold concentration, then increases to a maximum, and finally decreases again at higher ionic strength. The same non-monotonic curve is absent when the net charge is distributed uniformly, demonstrating that the spatial arrangement of the charges controls the effective protein-protein interaction once screening becomes strong.","pith_inferences":["The same charge-pattern mechanism may underlie re-entrant condensation observed in other macro-ion systems.","Systematic B2 measurements over a wider salt window than is customary would locate the predicted maximum and test the discrete-charge explanation directly.","If the maximum survives modest changes in protein shape or charge mobility, it offers an additional experimental knob for tuning solubility without altering pH."],"forward_implications":["Protein crystallization conditions can be located from the position of the B2 maximum rather than from simple screened-Coulomb estimates.","DLVO and nonlinear Poisson-Boltzmann theory cannot be used to extrapolate B2 once salt exceeds roughly 0.1 M for typical protein net charges.","Short-range ion-protein correlations, omitted by continuum theories, set both the sign and the magnitude of B2 at high ionic strength."],"fun_headline_variants":["Discrete charges trigger non-monotonic B2 in salted proteins","Protein B2 falls then peaks then drops with added salt","Charge pattern produces B2 maximum at moderate ionic strength","B2 salt dependence turns non-monotonic beyond DLVO limits"],"cache_read_input_tokens":128,"weakest_assumption_plain":"A rigid sphere carrying a fixed, discrete charge pattern whose net charge and size match lysozyme is sufficient to capture the dominant physics that sets B2 across the examined salt range.","fun_headline_variants_meta":{"raw":{"variants":["Discrete charges trigger non-monotonic B2 in salted proteins","Protein B2 falls then peaks then drops with added salt","Charge pattern produces B2 maximum at moderate ionic strength","B2 salt dependence turns non-monotonic beyond DLVO limits"]},"model":"grok-4.3","cost_usd":0.002154,"raw_usage":{"total_tokens":1224,"prompt_tokens":617,"num_sources_used":0,"completion_tokens":65,"cost_in_usd_ticks":21543500,"prompt_tokens_details":{"text_tokens":617,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":542,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":617,"tokens_out":65,"duration_ms":10293,"temperature":1.0,"reasoning_tokens":542,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-14T20:56:47.842422+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Measure the osmotic second virial coefficient of lysozyme at fixed pH while scanning NaCl concentration from 0.01 M to 1 M; the data must display a local minimum followed by a distinct maximum if the reported non-monotonicity is correct.","supporting_citations":[],"review_version":1}