REVIEW 2 major objections 1 cited by
Elongated drying droplets create axial and transverse deposit inhomogeneities; longer nanowires improve connectivity and uniformity while stronger attraction boosts conductivity at the cost of clustering.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-10 21:10 UTC pith:4IKXOACW
load-bearing objection Solid CG-LBM extension to elongated nanowire lines: geometry sets axial/transverse inhomogeneity, length helps connectivity and uniformity while cohesion trades them off—useful process insight, but conductivity is geometric and self-pinning is neglected. the 2 major comments →
Evaporation-Driven Nanowire Self-Assembly in an Elongated Droplet
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The elongated droplet geometry intrinsically induces distinct axial and transverse inhomogeneities in the final nanowire deposit through anisotropic contraction and capillary flow. Increasing effective inter-nanowire attraction improves electrical connectivity via aggregation and contact formation, yet promotes clustering and local ordering that reduce structural uniformity. Increasing nanowire length yields a dual benefit: it improves long-range connectivity through more stable percolating backbones while enhancing deposit homogeneity by bridging gaps and suppressing excessive local densification.
What carries the argument
Mesoscale lattice Boltzmann color-gradient fluid dynamics two-way coupled to bead–spring filament models of nanowires, with final bead configurations mapped onto resistor networks that quantify relative conductivity and current pathways. This machinery links the two-stage (axial-then-radial) drying sequence and filament parameters to deposit morphology and transport.
Load-bearing premise
The model assumes filaments are dilute enough that deposited nanowires do not pin the contact line; if real inks self-pin, the two-stage flow sequence and the dual-benefit claim for length would change.
What would settle it
Dry elongated nanowire droplets of short versus long filaments under dilute conditions matching the simulations and measure both spatial uniformity and long-range conductivity of the deposits; if longer filaments fail to raise both connectivity and homogeneity relative to short ones, the dual-benefit claim is false. Separately, if contact-line self-pinning appears at the simulated concentrations, the dilute no-self-pinning premise fails.
If this is right
- Substrate receding contact angle can be used to tune axial deposit extent, coffee-ring strength, and whether material accumulates at edges or in the center.
- When both percolation and spatial uniformity matter for printed conductive lines, longer nanowires are preferable to stronger inter-wire attraction.
- Solvent or surface-chemistry changes that increase nanowire attraction can raise conductivity, but only by accepting more clustered, less uniform deposits.
- Geometry-imposed two-stage drying sets the global deposition pathway that microscopic filament parameters only modulate.
Where Pith is reading between the lines
- If practical inks are not dilute, self-pinning by deposited nanowires would likely lock in stronger edge-aligned structures and weaken or reverse the dual-benefit claim for length.
- The same axial-then-radial sequence should appear for other anisotropic footprints common in aerosol-jet or blade coating, making footprint geometry a general control knob beyond the rectangular patches studied.
- Resistor-network maps of dried deposits could be extended to predict anisotropic sheet resistance along versus across printed lines for circuit design.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses multi-component color-gradient lattice Boltzmann simulations coupled to bead–spring filaments to study evaporation-driven self-assembly of nanowires in elongated droplets on line-shaped hydrophilic patches. It establishes a two-stage drying pathway (axial contraction then radial recession) that produces distinct axial and transverse deposit inhomogeneities, then maps final morphologies onto a geometric resistor network to quantify relative conductivity κ*. Systematic sweeps of receding contact angle θ_r, filament length L, and normalized cohesion ε* show that stronger inter-filament attraction raises connectivity at the cost of clustering and reduced uniformity, whereas longer filaments improve long-range connectivity while also enhancing spatial homogeneity. The authors present this as design guidance for balancing transport and uniformity in printed-electronics line deposits.
Significance. If the reported trends hold under the stated idealizations, the work fills a genuine gap between coffee-ring literature on circular drops and the elongated footprints typical of inkjet/aerosol-jet printing. Strengths include a well-documented mesoscale framework (prior method papers plus open Zenodo data), multi-run averaging (20–40 seeds), explicit structural metrics (segment-wise nematic order S_x^b and density g^b), and a transparent link from drying pathway to an effective network conductivity. The dual-benefit claim for filament length and the geometry-imposed axial/transverse anisotropy are concrete, falsifiable design rules that the printed-electronics community can test. The geometric resistor model and dilute no-self-pinning assumption limit direct transfer to junction-dominated metallic nanowire inks, but the morphological results remain useful even if absolute κ* values are idealized.
major comments (2)
- Methods, Electrical network model (and Figs. 6–7, Conclusion): The dual-benefit length claim and the cohesion–uniformity trade-off rest on a purely geometric resistor network (nodes = beads, edges for r_ij < r_c = 2.1 with G_ij = 1/r_ij) that “does not explicitly distinguish between bulk filament conduction and inter-filament contact resistance.” In many metallic nanowire/CNT inks junction resistance dominates; under that regime length-driven bridging may raise effective conductivity less, and cohesion-driven clustering (which multiplies junctions) could reverse the reported trade-off. The design guidance should either (i) restate κ* as a geometric connectivity proxy rather than electrical transport, or (ii) add a sensitivity study with elevated contact resistance (or a two-parameter bulk/junction model) to show which qualitative trends survive.
- Methods, Problem Definition and Assumptions: Self-pinning of the contact line by previously deposited filaments is neglected under a dilute-filament assumption (~4%). If deposited nanowires pin the contact line in real inks, the two-stage axial-then-radial flow sequence (Fig. 2–3) and the resulting dual-benefit length claim would change. A short discussion or a limited higher-concentration/self-pinning test is needed to bound the regime of validity of the geometry-imposed pathway.
Circularity Check
No load-bearing circularity: claims arise from new forward LB simulations of elongated filament-laden droplets; self-citations supply reusable methods infrastructure and the single saturation curve is explicitly labeled a descriptive phenomenological fit.
specific steps
-
self citation load bearing
[Methods §I.A (Fluid-filament simulation method), first paragraph]
"We employ the multi-component lattice Boltzmann color-gradient (CG) method [25] to model fluid dynamics. … The same framework is adopted here with minor adjustments, and the relevant elements are summarized below. … refer to our earlier publications for a detailed derivation and validation [11, 27–29]."
The numerical engine is taken from the authors’ own prior papers. This is ordinary methods reuse and is not load-bearing for the scientific claims: the elongated-geometry, cohesion, and length results are generated by new simulation runs, not by re-invoking a uniqueness or existence theorem from those citations. Flagged only as the mildest self-citation pattern; it does not force the reported deposit morphologies or dual-benefit conclusion.
full rationale
The paper’s central claims (geometry-imposed axial/transverse inhomogeneity; dual benefit of longer filaments for connectivity plus homogeneity; cohesion improves κ* at the cost of clustering) are obtained by running new color-gradient lattice-Boltzmann simulations of evaporating elongated droplets containing bead–spring filaments, then post-processing the final bead configurations with a geometric resistor network. The simulation framework is taken from the authors’ prior methods papers, but those citations supply only the numerical machinery (CG collision, CSF surface tension, FENE/WCA filaments, solvation force, height-dependent evaporation flux); they do not encode or presuppose the elongated-deposit morphology or the dual-benefit length result. The sole fitted expression (Eq. 19) is introduced after the simulation data are shown, is called a “phenomenological exponential saturation model” and “descriptive fit rather than as a microscopic transport theory,” and is never used to generate a claimed first-principles prediction. No uniqueness theorem, self-definitional identity, or ansatz smuggled via self-citation appears in the derivation chain. The geometric idealization of the resistor network (Gij = 1/rij, no bulk/junction distinction) is a modeling assumption that may affect physical correctness, but it is not circular: the measured κ* is simply the output of that defined network applied to the simulated configurations. Hence the circularity score is at most 1.
Axiom & Free-Parameter Ledger
free parameters (6)
- Normalized cohesion ε* = ε/ε_s (and absolute LJ well depths ε, ε_s=0.006)
- Evaporation kinetic resistance K=0.01 and scale J0 in height-dependent flux J≈J0/(K+h̃)
- Filament mechanical parameters ks=0.3, Rm=2.4, kb=8.0, d0=2, WCA ε=0.03
- Resistor cutoff rc=2.1 and Gij=1/rij; electrode span fixed at 120 lu
- Saturation-fit parameters κ0*, κ∞*, A in κ*(ε*)
- Initial geometry rx0=116, ry0=42, h0≈30; filament lengths L=12,32; concentration ~4%
axioms (7)
- domain assumption Isothermal Stokes/quasi-static regime: small Ca and Bo; interface relaxes much faster than evaporation; inertial effects negligible.
- domain assumption Brownian motion omitted (high Péclet); advection dominates diffusion.
- domain assumption Marangoni stresses negligible for water droplets.
- domain assumption Self-pinning by deposited filaments neglected (dilute concentration).
- ad hoc to paper Height-dependent evaporative flux (Eq. 8) adequately represents diffusion-limited drying of the elongated drop (with elliptic first-order correction vanishing for n=2).
- ad hoc to paper Effective resistor network with geometric conductance Gij=1/rij captures relative electrical transport without separating bulk vs contact resistance.
- domain assumption Color-gradient multiphase LB with CSF tension and geometry-aware wetting is a valid continuum description of the free-surface flow.
read the original abstract
Drying of nanowire-laden elongated droplets is a ubiquitous process in printed electronics fabrication, where the resulting deposition pattern critically determines device performance by controlling nanowire alignment, connectivity, and percolating charge-transport pathways. However, the physical understanding of evaporation-driven deposition is still largely derived from studies of spherical droplets on homogeneous substrates. This gap limits the ability to predict and control deposit morphology in realistic printing scenarios. Here, we use mesoscale lattice Boltzmann simulations to investigate the drying of nanowire-laden elongated droplets on wettability-patterned substrates, focusing on the effects of droplet geometry, nanowire interactions, and nanowire length. The elongated droplet geometry is found to intrinsically induce distinct axial and transverse inhomogeneities in the final deposit. Increasing the effective attraction between nanowires, which mimics changes in surface chemistry or solvent conditions, can improve electrical connectivity but also promotes clustering and local ordering, reducing structural uniformity. In contrast, increasing nanowire length yields a dual benefit by improving long-range connectivity while simultaneously enhancing deposit homogeneity. Our findings provide design guidance for balancing electrical transport and structural uniformity in evaporation-driven printed electronics.
Figures
Forward citations
Cited by 1 Pith paper
-
Coalescence-induced alignment of anisotropic particles in drying sessile droplets
Coalescence of sessile droplets aligns suspended dumbbells along the coalescence direction, and the drying stage either preserves or weakens that alignment depending on contact angle and substrate friction.
Reference graph
Works this paper leans on
-
[1]
Nematic order and density profile We study the effect of the receding contact angle on the final deposition pattern. Fig. 4 shows top-view snap- shots of the fully deposited non-cohesive filaments for receding contact angles ofθ r = 20 ◦ andθ r = 30 ◦. The color scale represents the intramolecular nematic order Sx, where red (Sx = 1) indicates axial align...
-
[2]
Conductivity Here, we relate the deposit morphology to the result- ing relative electrical conductivity, focusing on the effects of filament cohesion and length. Fig. 6 shows the final deposition morphologies for different normalized filament cohesion strengths,ε ∗ :=ε/ε s ∈[0.3,1.3], together with the purely repulsive reference case described by the WCA ...
-
[3]
M. Steinberger, Q. Xie, O. J. J. Ronsin, P. Maisch, K. C. Tam, A. Distler, J. Harting, C. J. Brabec, and H. Egel- haaf, Challenges and opportunities in upscaling inkjet- printing of OPV, Flex. Print. Electron.9, 043001 (2024)
work page 2024
-
[4]
Q. Huang and Y. Zhu, Printing conductive nanomaterials for flexible and stretchable electronics: A review of mate- rials, processes, and applications, Adv. Mater. Technol. 4, 1800546 (2019)
work page 2019
-
[5]
B. I. Guyll, B. L. Sanford, C. L. Pint, and E. B. Secor, Controlling droplet evaporation in aerosol jet printing to understand and mitigate overspray, Small Sci.5, 2500069 (2025)
work page 2025
-
[6]
R. D. Deegan, O. Bakajin, T. F. Dupont, G. Huber, S. R. Nagel, and T. A. Witten, Capillary flow as the cause of ring stains from dried liquid drops, Nature389, 827 (1997)
work page 1997
-
[7]
R. G. Larson, Transport and deposition patterns in dry- ing sessile droplets, AIChE J.60, 1538 (2014)
work page 2014
-
[8]
M. Anyfantakis and D. Baigl, Manipulating the coffee- ring effect: Interactions at work, ChemPhysChem16, 2726 (2015)
work page 2015
- [9]
-
[10]
P. J. Yunker, T. Still, M. A. Lohr, and A. G. Yodh, Sup- pression of the coffee-ring effect by shape-dependent cap- illary interactions, Nature476, 308 (2011)
work page 2011
-
[11]
Y. Wang, M. Jian, H. Liu, and X. Zhang, Anisotropic wetting of droplets on stripe-patterned chemically het- erogeneous surfaces: Effect of length ratio and deposition position, Langmuir35, 4387 (2019)
work page 2019
-
[12]
A. Aboubakri, Y. Akkus, A. K. Sadaghiani, K. Sefiane, and A. Ko¸ sar, Computational and experimental investi- gations on the evaporation of single and multiple elon- gated droplets, Chem. Eng. J. Adv.10, 100255 (2022)
work page 2022
-
[13]
Q. Xie, T. Du, C. J. Brabec, and J. Harting, Effect of particle and substrate wettability on evaporation-driven assembly of colloidal monolayers, Langmuir41, 14995 (2025)
work page 2025
-
[14]
Y. Mino, C. Tanaka, H. Tanaka, K. Nakaso, and K. Go- toh, Numerical simulation of a drying colloidal suspen- sion on a wettable substrate using the lattice boltzmann method, Chem. Eng. Sci.263, 118050 (2022)
work page 2022
-
[15]
M. Hartmann and S. Hardt, Stability of evaporating droplets on chemically patterned surfaces, Langmuir35, 4868 (2019)
work page 2019
-
[16]
P. Kabi, R. Pal, and S. Basu, Moses effect: Splitting a sessile droplet using a vapor-mediated marangoni effect leading to designer surface patterns, Langmuir36, 1279 (2020)
work page 2020
-
[17]
N. T. Dinh, E. Sowade, T. Blaudeck, S. Hermann, R. D. Rodriguez, D. R. T. Zahn, S. E. Schulz, R. R. Baumann, and O. Kanoun, High-resolution inkjet printing of con- ductive carbon nanotube twin lines utilizing evaporation- driven self-assembly, Carbon96, 382 (2016)
work page 2016
-
[18]
P. L. Kumar, S. P. Thampi, and M. G. Basavaraj, Par- ticle size and substrate wettability dependent patterns in dried pendant drops, J. Phys.: Condens. Matter33, 024003 (2021). 12
work page 2021
-
[19]
C. P. Whitby and A. Hermant, Concentration of de- posit patterns by nanoparticles modified with short am- phiphiles, Colloids Surf. A594, 124648 (2020)
work page 2020
-
[20]
R. Bhardwaj, X. Fang, P. Somasundaran, and D. At- tinger, Self-assembly of colloidal particles from evaporat- ing droplets: Role of DLVO interactions and proposition of a phase diagram, Langmuir26, 7833 (2010)
work page 2010
-
[21]
G. L. Goh, S. Agarwala, and W. Y. Yeong, Aerosol-jet- printed preferentially aligned carbon nanotube twin-lines for printed electronics, ACS Appl. Mater. Interfaces11, 43719 (2019)
work page 2019
-
[22]
L. Tu, S. Yuan, H. Zhang, P. Wang, X. Cui, J. Wang, Y.-Q. Zhan, and L.-R. Zheng, Aerosol jet printed silver nanowire transparent electrode for flexible electronic ap- plication, J. Appl. Phys.123, 174905 (2018)
work page 2018
-
[23]
C. Belgardt, E. Sowade, T. Blaudeck, T. Baumg¨ artel, H. Graaf, C. von Borczyskowski, and R. R. Baumann, Inkjet printing as a tool for the patterned deposition of octadecylsiloxane monolayers on silicon oxide surfaces, Phys. Chem. Chem. Phys.15, 7494 (2013)
work page 2013
-
[24]
S. M. Park and D. K. Yoon, Evaporation-induced self- assembly of liquid crystal biopolymers, Mater. Horiz.11, 1843 (2024)
work page 2024
-
[25]
Y. Xia, B. P. Yalagala, A. S. Karimullah, H. Heidari, and R. Ghannam, Beyond flexibility: Transparent silver nanowire electrodes on patterned surfaces for reconfig- urable devices, Adv. Eng. Mater.26, 2301165 (2024)
work page 2024
-
[26]
X. Ye, L. Fei, L. Lu, and C. Li, Influence of anisotropic nanoparticles on the deposition pattern of an evaporating droplet, Eur. Phys. J. E42, 17 (2019)
work page 2019
-
[27]
A. K. Gunstensen, D. H. Rothman, S. Zaleski, and G. Zanetti, Lattice boltzmann model of immiscible fluids, Phys. Rev. A43, 4320 (1991)
work page 1991
-
[28]
S. Leclaire, A. Parmigiani, O. Malaspinas, B. Chopard, and J. Latt, Generalized three-dimensional lattice boltz- mann color-gradient method for immiscible two-phase pore-scale imbibition and drainage in porous media, Phys. Rev. E95, 033306 (2017)
work page 2017
-
[29]
G. Nath, O. Aouane, and J. Harting, Reaction-limited evaporation for the color-gradient lattice Boltzmann model, J. Chem. Phys.162, 114110 (2025)
work page 2025
- [30]
-
[31]
J. Sch¨ ottner, Q. Xie, G. Nath, and J. Harting, Self- assembled filament layers in drying sessile droplets: From morphology to electrical conductivity, Langmuir42, 8592 (2026)
work page 2026
-
[32]
J. U. Brackbill, D. B. Kothe, and C. Zemach, A contin- uum method for modeling surface tension, J. Comput. Phys.100, 335 (1992)
work page 1992
-
[33]
Q. Gu, J. Zhang, H. Liu, and L. Wu, Numerical study of droplet behavior passing through a constricted square channel, Phys. Fluids35, 073316 (2023)
work page 2023
- [34]
-
[35]
T. Akai, B. Bijeljic, and M. J. Blunt, Wetting bound- ary condition for the color-gradient lattice boltzmann method: Validation with analytical and experimental data, Adv. Water Resour.116, 56 (2018)
work page 2018
-
[36]
A. W. Wray and M. R. Moore, Evaporation of non- circular droplets, J. Fluid Mech.961, A11 (2023)
work page 2023
-
[37]
N. Murisic and L. Kondic, On evaporation of sessile drops with moving contact lines, J. Fluid Mech.679, 219 (2011)
work page 2011
-
[38]
C. Larsson and S. Kumar, Comparison of one-sided and diffusion-limited evaporation models for thin liquid droplets, J. Fluid Mech.976, A30 (2023)
work page 2023
-
[39]
P. Ahlrichs and B. D¨ unweg, Simulation of a single poly- mer chain in solution by combining lattice boltzmann and molecular dynamics, J. Chem. Phys.111, 8225 (1999)
work page 1999
-
[40]
M. Sega, M. Sbragaglia, S. S. Kantorovich, and A. O. Ivanov, Mesoscale structures at complex fluid–fluid in- terfaces: A novel lattice boltzmann/molecular dynamics coupling, Soft Matter9, 10092 (2013)
work page 2013
-
[41]
G. Ghigliotti, C. Zhou, and J. J. Feng, Simulations of the breakup of liquid filaments on a partially wetting solid substrate, Phys. Fluids25, 072102 (2013)
work page 2013
-
[42]
A. Dziedzic, M. Nakrani, B. Ezra, M. Syed, S. Popinet, and S. Afkhami, Breakup of finite-size liquid filaments: Transition from no-breakup to breakup including sub- strate effects, Eur. Phys. J. E42, 18 (2019)
work page 2019
- [43]
- [44]
- [45]
-
[46]
S. G. S. Ka, W. Wang, H. Giddens, Z. Chen, A. N. Khan, Y. Shui, A. S. Andy, S. Lyu, T. Hasan, Y. Hao, and Y. Y. S. Huang, Adaptive printing of conductive microfibers for seamless functional enhancement across diverse surfaces and shapes, Adv. Fiber Mater.7, 1274 (2025)
work page 2025
- [47]
-
[48]
Dai, Carbon nanotubes: Synthesis, integration, and properties, Acc
H. Dai, Carbon nanotubes: Synthesis, integration, and properties, Acc. Chem. Res.35, 1035 (2002)
work page 2002
-
[49]
B. Bessaire, M. Mathieu, V. Salles, T. Yeghoyan, C. Celle, J.-P. Simonato, and A. Brioude, Synthesis of continuous conductive PEDOT:PSS nanofibers by elec- trospinning: A conformal coating for optoelectronics, ACS Appl. Mater. Interfaces9, 950 (2017)
work page 2017
-
[50]
A. A. Hagberg, D. A. Schult, and P. J. Swart, Explor- ing network structure, dynamics, and function using net- workx, inProceedings of the 7th Python in Science Con- ference, edited by G. Varoquaux, T. Vaught, and J. Mill- man (Pasadena, CA, USA, 2008) pp. 11–15
work page 2008
-
[51]
P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Po- lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henr...
work page 2020
-
[52]
´A. G. Mar´ ın, H. Gelderblom, D. Lohse, and J. H. Snoei- jer, Order-to-disorder transition in ring-shaped colloidal stains, Phys. Rev. Lett.107, 085502 (2011)
work page 2011
-
[53]
M. Kr¨ oger, C. Luap, and P. Ilg, Ultra-slow self-similar coarsening of physical fibrillar gels formed by semiflexible polymers, Soft Matter21, 2803 (2025)
work page 2025
-
[54]
Azizian, Kinetic models of sorption: A theoretical analysis, J
S. Azizian, Kinetic models of sorption: A theoretical analysis, J. Colloid Interface Sci.276, 47 (2004)
work page 2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.