{"id":"48694dcb-7d71-412c-8ca1-b84ccfbd95e5","arxiv_id":"1908.03246","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Confined liquid crystal skyrmions in open traps numerically form binary pseudo-spins that relax into ice-rule states on square and hexagonal lattices.","lead":"A numerical study shows that confining liquid crystal skyrmions in open-ended traps creates binary variables that settle into ice-rule states on square and hexagonal lattices. The result proposes liquid crystals as a new, reconfigurable platform for artificial spin ice.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ice-rule result is computed in the 2D K=0 limit, while the claimed platform is a 3D confined LC cell; the main text does not demonstrate that this mapping preserves the spin-ice observables.","rationale":"The reader and I identify the same weakest assumption: the 2D K=0 model is the bridge between the simulated observable and the claimed 3D LC platform. I considered other candidates, including the overbroad 'first Skyrmion Spin Ice' claim relative to Ref. [18], the lack of repeated-run statistics, and the absence of code and data. Those are real but secondary: the novelty wording does not affect whether the ice-rule physics is correctly demonstrated, and the numerical protocol could be repeated without changing the physics. The 2D-to-3D mapping is different because it directly determines whether the proposed experimental system would behave as a spin ice. The paper's own statement that the main text is 2D, with 3D in the SM, makes this an acknowledged limitation rather than a hidden one. If the SM already contains a quantitative 3D comparison of vertex fractions for the same geometries, my proposed test would confirm it; if not, the test is the minimal check needed before the 'realized' claim can be accepted. The parameter scans and consistency with known spin-ice phenomenology (ordered square ground state, disordered hexagonal manifold) support the 2D result, so this is not a challenge to internal consistency. Nor do I claim the 2D model is wrong; I claim the main text does not yet demonstrate that it is quantitatively faithful in 3D. Hence the verdict stays CONDITIONAL.","tokens_in":8974,"tokens_out":6338,"duration_ms":66352,"concrete_test":"Run the full 3D Q-tensor relaxation of Eq. (1) on a square-lattice ice geometry matching the trap aspect ratio of Fig. 3 (e.g., a 5x5 vertex array, roughly 50 skyrmions) at Nz=0.36p with K=K0 and alpha in the skyrmion-stability range, using cylindrical obstacles with homeotropic anchoring, and apply the same swell/deswell protocol. Measure the defect ratio and compare it with the 2D K=0 result at the same omega/D. If the 3D defect ratio is more than 10 percentage points above the 2D value, or lies outside the green ice-rule region of Fig. 4(a), the 2D K=0 model is not a faithful proxy and the experimental realization claim is not established. Repeat with two independent random initial skyrmion placements to also address the missing repetition statistics.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central ice-rule result is obtained from Eq. (1) in the 'special case of K = 0' that 'yields a z-invariant structure' and is modeled in 2D, while the paper states 'Main text includes only 2D simulations and full 3D simulations are presented in SM.' The experimental cell, however, has finite thickness (Nz/p ≈ 0.36), homeotropic anchoring, and, for K ≠ 0, barrel-like 3D skyrmions. The 2D K=0 skyrmions are not distorted along z and their mutual repulsion is computed from a z-invariant texture; in a real cell the cholesteric twist, boundary layers, and obstacle anchoring can alter both the binary trap potential and the range or anisotropy of skyrmion-skyrmion forces. Because the ice-rule states in Fig. 3 emerge specifically from frustrated mutual repulsion and binary occupancy, a material change in those interactions in 3D could suppress the type-IV and type-II vertex fractions. The manuscript asserts the K=0 mapping as a computational convenience rather than demonstrating that it is quantitatively faithful for the spin-ice observables; the 3D evidence is cited but not shown in the main text. This is the load-bearing gap: if 3D simulations do not reproduce the 2D defect fractions, the proposal is not validated as a real LC platform.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a liquid-crystal realization of artificial spin ice in which chiral nematic skyrmions are confined in binary, open-ended traps and interact through mutual repulsion. The authors simulate ensembles of skyrmions on square and hexagonal lattices using a Landau–de Gennes free energy, Eq. (1), starting from random skyrmion placements and employing a swell-and-deswell relaxation protocol. They report that the relaxed states obey the ice rule: square ice organizes into an ordered array of Type IV vertices with domain walls, and hexagonal ice relaxes into a disordered manifold of Type II and Type III vertices. The paper also presents parameter studies of trap aspect ratio, obstacle strength, skyrmion size, and quenched disorder. The central claim is that this is the first skyrmion spin ice and a viable new platform for frustrated soft-matter systems.","tokens_in":9244,"tokens_out":4180,"duration_ms":46081,"significance":"If the central claim holds, this is a substantive contribution to artificial spin ice: liquid-crystal skyrmions offer reconfigurable traps, optical creation and annihilation, size tunability, and controllable interactions, which are not readily available in magnetic or colloidal platforms. The modeling has real strengths. The free energy in Eq. (1) is a standard cholesteric Landau–de Gennes energy with no parameters fitted to the ice-rule states; the simulated systems are large (288 and 192 skyrmions); the initial configurations are random; and the paper explicitly examines parameter dependence and disorder effects. The main limitation is that the central ice-rule results are computed in the 2D, K = 0 limit, while the proposed experimental platform is a finite-thickness 3D cell with homeotropic anchoring and, for K ≠ 0, barrel-like skyrmions. Establishing quantitative fidelity of the 2D mapping for the spin-ice observables is therefore load-bearing for the platform claim.","major_comments":[{"comment":"The central demonstration of ice-rule states is performed for the K = 0, z-invariant limit of Eq. (1), modeled in 2D, whereas the proposed experimental cell has thickness Nz/p ≈ 0.36, homeotropic anchoring at both boundaries, and barrel-like three-dimensional skyrmions when K ≠ 0. The effective binary trap potential and the range or anisotropy of skyrmion-skyrmion repulsion could change in a 3D texture because of cholesteric twist, boundary layers, and obstacle anchoring; the ice-rule statistics in Fig. 3 emerge specifically from frustrated mutual repulsion and binary occupancy. The manuscript therefore needs to show in the main text, or explicitly cite with quantitative values from the SM, that 3D simulations reproduce the same vertex fractions and defect ratios as the 2D K = 0 model. As written, the platform claim rests on an unverified dimensional-reduction assumption.","section":"Fig. 3 and the statement 'Main text includes only 2D simulations and full 3D simulations are presented in SM'"},{"comment":"The reported defect ratios appear to come from single realizations or a very small number of runs, with no error bars or ensemble averages. The text itself acknowledges that the system 'preserves memory of its preparation' and that cyclic swelling and deswelling changes the final state, which implies that relaxation is history-dependent. Without repeated independent initializations and a statement of run-to-run variability, the reader cannot assess whether the ice-rule behavior is robust or an artifact of a particular random starting configuration and swell-deswell schedule. Please provide at least several independent runs for the key parameter values and report mean and standard deviation of the defect fraction.","section":"Fig. 3 and Fig. 4(a,b): defect statistics"},{"comment":"The quantitative comparability of the parameter scans is compromised by the admitted differences in system size and simulation time across panels: Fig. 4(c,d) use systems four times smaller than Fig. 4(a,b), and Fig. 4(e) uses a simulation time four times longer. This makes the apparent zero-defect regions and the disorder-dependence comparison difficult to interpret as a function of the control parameters alone. The authors should either repeat the relevant scans with consistent system size, simulation time, and relaxation protocol, or provide a scaling analysis showing that these differences do not affect the conclusions.","section":"Fig. 4 and the paragraph beginning 'One last note'"}],"minor_comments":[{"comment":"The word 'gemoetries' should be 'geometries'.","section":"Conclusion"},{"comment":"The caption and text describe the light-exposure curve as modeled by q0 → q0/1.2, which 'produces an effect similar to changing K → 1.2K'; please clarify in the caption whether the plotted curve uses the q0 reduction or the equivalent K change, and define the dimensionless units used on the horizontal axis.","section":"Fig. 4(c)"},{"comment":"The phrase 'previously studies in detail' should read 'previously studied in detail'.","section":"Page 4"},{"comment":"The main text repeatedly refers to the SM for simulation parameters but does not list the numerical values of L, a, b, c, Γ, lattice size, and grid spacing needed to reproduce the figures. Including a parameter table in the main text would improve reproducibility.","section":"Eq. (1) and SM"},{"comment":"The description of Type IV vertices as '2 skyrmions in the vertex, and two out of the vertex' would be clearer if the spin convention for each lattice edge were stated explicitly, since the square and hexagonal vertex types are defined relative to the four or three incident traps respectively.","section":"Fig. 2(b,c)"}],"recommendation":"major_revision","confidential_remarks":"The paper is appropriate for cond-mat.soft and will interest the artificial-spin-ice community. The main concern is the gap between the 2D K = 0 simulations shown in the main text and the 3D experimental platform being proposed. If the SM already contains a full 3D vertex analysis that reproduces the 2D results, the revision is straightforward: bring that comparison into the main text. If it does not, the authors need to perform the 3D simulations or substantially soften the platform claim. I have no concerns about novelty disclosure; the authors cite prior work on skyrmion confinement and on particle-based ice appropriately."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does what it says: it shows numerically that confined LC skyrmions in open-ended traps relax to ice-rule states on square and hexagonal lattices. The open-end trap design is a smart fix to the closed-dumbbell geometry used in colloidal ice, and the parameter scans for aspect ratio, field strength, and disorder are genuinely informative. Square ice orders antiferromagnetically, hexagonal ice stays disordered—exactly what the ice rule predicts. That is a solid numerical demonstration of a known physics in a new platform, not a new theoretical framework, and I think the authors are appropriately measured in most of their claims.\n\nThe soft spots are real but not fatal. The main one is exactly what the stress test flags: the headline results come from 2D K=0 simulations, while the experimental cell is 3D with finite thickness and anchoring. The authors are transparent about this and point to 3D results in the SM, but the main text does not show that the K=0 limit preserves the binary trap potential and repulsion that drive the ice-rule behavior. A referee should push for that evidence in the main text or a more explicit argument for why the 3D corrections do not change the vertex energetics. This is a load-bearing gap but not a collapse: the K=0 structure is a known limit of the same Landau–de Gennes model, and the physics likely survives, but the claim 'realized via confined skyrmions' needs the 3D check to be presented rather than deferred.\n\nThe other weaknesses are minor in comparison. There are no repeated-run statistics; each parameter point appears to be a single relaxation, so defect ratios in Fig. 4 could be noisy. The 'first Skyrmion Spin Ice' claim is slightly overbroad given the cited magnetic skyrmion work [18], though the LC-specific realization is new. No code or data is provided, which limits reproducibility. The citation practice is fine: the self-citations are to the authors' own earlier work on LC skyrmion stability and confinement, and they are used appropriately.\n\nWho is this for? Researchers in artificial spin ice and soft matter topology. It is a good example of exporting spin-ice physics to a reconfigurable platform, and it suggests concrete experimental directions. I would bring it to a reading group.\n\nMy recommendation: send it to peer review. It is a serious, well-posed proposal with honest limitations. A referee should request 3D simulation results in the main text or a rigorous justification of the K=0 mapping, and repeated runs for the defect statistics. With those, it would be a solid contribution.","headline":"A credible numerical proposal for liquid-crystal skyrmion spin ice with a real 2D-to-3D gap, worth refereeing despite missing statistics.","tokens_in":9810,"tokens_out":1916,"would_cite":true,"duration_ms":21654,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that liquid-crystal skyrmions confined in open-ended traps act as binary pseudo-spins whose mutual repulsion drives the system into ice-rule states on square and hexagonal lattices.","keywords":["artificial spin ice","liquid crystal skyrmions","ice rule","geometric frustration","topological solitons","chiral nematic","particle ice","Q-tensor model"],"falsifier":"Run a full three-dimensional simulation, or a chiral-nematic cell experiment, of the same square trap array with intermediate aspect ratio and measure the vertex statistics after relaxation: if the fraction of ice-rule (type-IV) vertices is not substantially higher than random, with defect ratios near 1 rather than the low values reported in the paper, the central claim fails.","tokens_in":8737,"feed_emoji":"🌀","tokens_out":10815,"duration_ms":106191,"temperature":0.7,"pith_summary":"Liquid-crystal skyrmions are particle-like topological solitons in the director field of a chiral nematic. The paper proposes that, when confined in open-ended traps, they behave as binary pseudo-spins whose mutual elastic repulsion is frustrated at lattice vertices, producing artificial spin ice. Two large numerical relaxations--288 skyrmions on a square lattice and 192 on a hexagonal lattice--converge to ice-rule-obeying states: an ordered antiferromagnetic tessellation of type-IV vertices on the square lattice and a disordered manifold of 2-in/1-out and 1-in/2-out vertices on the hexagonal lattice. If the two-dimensional model captures real three-dimensional confinement, this would make liquid crystals a reconfigurable platform for studying frustration, with trap geometry and skyrmion size tunable in place.","feed_headline":"Skyrmion spin ice: liquid-crystal solitons obey ice rules","feed_subtitle":"Simulations show square ice ordering and hexagonal ice staying disordered, a reconfigurable frustrated-matter platform.","key_machinery":"The central object is the liquid-crystal skyrmion, a particle-like soliton in which the director rotates by 180 degrees from core to periphery, carrying unit topological charge. The key construction is the open-ended binary trap: a channel whose ends confine a skyrmion to either of two positions, with the middle left open so that skyrmions on neighboring traps can interact through their elastic fields; the paper explains that closed dumbbell traps would suppress this interaction and defeat the spin ice. The simulations evolve a Landau-de Gennes $Q$-tensor free energy with cholesteric twist, surface anchoring, and an electric field, using the $z$-invariant $K=0$ limit to make large two-dimensional systems tractable, and implement the relaxation on GPUs.","core_discovery":"The central claim is that a chiral nematic liquid crystal can realize a skyrmion spin ice: arrays of open-ended traps, each holding at most one skyrmion, define binary variables, and the elastic repulsion between skyrmions imposes the ice rule at every vertex. In overdamped relaxations, square-lattice arrays settle into an ordered 'antiferromagnetic' state of type-IV vertices (two skyrmions in, two out), with defects concentrated in domain walls, while hexagonal arrays settle into a disordered mixture of type-II and type-III vertices (2-in/1-out and 1-in/2-out) with sparse monopole defects. The paper identifies a window of trap aspect ratios where ice behavior is best, shows that cyclically swelling and deswelling skyrmions lowers defect counts, and finds that quenched disorder hurts square ice more than hexagonal ice. The main-text evidence is two-dimensional, using the $K=0$ $z$-invariant limit; full three-dimensional simulations are presented in the supplementary material. On this basis the authors state that they have demonstrated numerically that liquid crystals are a new platform for spin-ice physics on the two most common geometries.","pith_inferences":["If the 2D $z$-invariant model carries over to finite-thickness cells, the defect-ratio curves in the paper give a quantitative target for an experiment: measuring vertex statistics as a function of trap aspect ratio and field strength would test the platform directly.","The open-ended trap design suggests a general principle for particle-based ices: confinement must not screen the inter-particle interaction, so traps should localize a particle without isolating it from its neighbors.","The wide-trap regime that spoils the binary character could be a practical route to classical spin-1 ice (vertex occupations 0, 1, or 2), a model that is hard to reach with magnetic islands.","The stronger sensitivity of square ice to quenched disorder points to a design heuristic: choose hexagonal lattices when robustness to fabrication disorder matters, and square lattices when an ordered ground state is desired."],"forward_implications":["Square-lattice skyrmion ice should relax to an ordered antiferromagnetic state of type-IV vertices, with ice-rule violations appearing mainly inside domain walls.","Hexagonal-lattice skyrmion ice should remain a disordered manifold of 2-in/1-out and 1-in/2-out vertices, with sparse monopole defects.","Ice-rule behavior is conditional on trap shape and skyrmion size: traps that are too narrow freeze the particles, traps that are too wide allow centered positions, and intermediate aspect ratios are required.","Cyclically swelling and deswelling the skyrmions anneals the system and extends the range of trap parameters over which the ice rule is reached.","Because skyrmions can be created, annihilated, resized, and steered optically or electrically, the same sample can be reconfigured between different ice geometries, enabling memory, doping, and decimation experiments."],"supporting_citations":[{"why":"It supplies the original colloidal artificial ice construction: binary traps as pseudo-spins and repulsion enforcing ice rules.","marker":"[3]"},{"why":"It provides the theory of non-local frustration in particle-based ice that motivates why repelling particles should obey the ice rule.","marker":"[21]"},{"why":"It gives the pseudo-spin language and vertex-type classification for particle ice used throughout the paper.","marker":"[17]"},{"why":"It demonstrates experimentally that two-dimensional skyrmions exist in confined chiral nematics, grounding the quasi-particle picture.","marker":"[24]"},{"why":"It establishes the 2D $z$-invariant skyrmion model and the alignment-induced attraction and repulsion used to build traps.","marker":"[26]"},{"why":"It determines the electric-field stability range and size control of skyrmions used in the swelling and deswelling protocol.","marker":"[28]"},{"why":"It shows commensurate trapping of liquid-crystal skyrmions in a square array, the direct basis for the binary traps.","marker":"[41]"}],"fun_headline_variants":["Liquid-crystal skyrmions make a reconfigurable spin ice","Skyrmions in liquid crystals obey ice rules","Reconfigurable spin ice from LC skyrmions","LC skyrmion traps give square and hexagonal ice","Skyrmion spin ice: frustrated solitons in liquid crystals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the two-dimensional, $z$-invariant model (the special $K=0$ case) captures the same binary-trap confinement and mutual repulsion that real three-dimensional skyrmions would experience in a cell of finite thickness, so that ice-rule behavior predicted in 2D survives in experiments.","fun_headline_variants_meta":{"raw":{"variants":["Liquid-crystal skyrmions make a reconfigurable spin ice","Skyrmions in liquid crystals obey ice rules","Reconfigurable spin ice from LC skyrmions","LC skyrmion traps give square and hexagonal ice","Skyrmion spin ice: frustrated solitons in liquid crystals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000262,"raw_usage":{"total_tokens":1563,"prompt_tokens":877,"completion_tokens":686,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":604}},"tokens_in":493,"tokens_out":686,"duration_ms":6344,"temperature":1.0,"reasoning_tokens":604,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:19:54.488736+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a full three-dimensional simulation, or a chiral-nematic cell experiment, of the same square trap array with intermediate aspect ratio and measure the vertex statistics after relaxation: if the fraction of ice-rule (type-IV) vertices is not substantially higher than random, with defect ratios near 1 rather than the low values reported in the paper, the central claim fails.","supporting_citations":[{"cited_title":"Realizing Colloidal Artificial Ice on Arrays of Optical Traps,","cited_arxiv_id":null,"evidence_quote":"It supplies the original colloidal artificial ice construction: binary traps as pseudo-spins and repulsion enforcing ice rules."},{"cited_title":"Unexpected phenomenology in particle-based ice absent in magnetic spin ice,","cited_arxiv_id":null,"evidence_quote":"It provides the theory of non-local frustration in particle-based ice that motivates why repelling particles should obey the ice rule."},{"cited_title":"Col- loquium: Ice rule and emergent frustration in particle ice and beyond,","cited_arxiv_id":null,"evidence_quote":"It gives the pseudo-spin language and vertex-type classification for particle ice used throughout the paper."},{"cited_title":"Two- dimensional skyrmions and other solitonic structures in confinement-frustrated chiral nematics,","cited_arxiv_id":null,"evidence_quote":"It demonstrates experimentally that two-dimensional skyrmions exist in confined chiral nematics, grounding the quasi-particle picture."},{"cited_title":"Alignment induced re-configurable walls for patterning and assembly of liquid crystal skyrmions","cited_arxiv_id":"2001.09615","evidence_quote":"It establishes the 2D $z$-invariant skyrmion model and the alignment-induced attraction and repulsion used to build traps."},{"cited_title":"Comparing skyrmions and merons in chiral liquid crystals and magnets,","cited_arxiv_id":null,"evidence_quote":"It determines the electric-field stability range and size control of skyrmions used in the swelling and deswelling protocol."},{"cited_title":"Commensurate states and pat- tern switching via liquid crystal skyrmions trapped in a square lattice,","cited_arxiv_id":null,"evidence_quote":"It shows commensurate trapping of liquid-crystal skyrmions in a square array, the direct basis for the binary traps."}],"review_version":1}