research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

Velocity gradient dependent structures of 12–2–12 wormlike micelles: insights from small-angle neutron scattering in the 1–2 shear plane

crossmark logo

aNeutron Science and Technology Center, Comprehensive Research Organization for Science and Society (CROSS), 162-1 Shirakata, Tokai, Ibaraki 319-1106, Japan, bMaterials and Life Science Division, J-PARC Center, Japan Atomic Energy Agency, 2-4 Shirakata, Tokai, Ibaraki 319-1195, Japan, and cDepartment of Chemistry, Faculty of Science, Nara Women's University, Kitauoyanishi-machi, Nara 630-8506, Japan
*Correspondence e-mail: [email protected]

Edited by J. Houston, European Spallation Source, Sweden (Received 24 January 2025; accepted 21 July 2025; online 29 July 2025)

This article is part of a collection of articles related to the 19th International Small-Angle Scattering Conference (SAS2024) in Taipei, Taiwan.

The velocity gradient dependent shear-induced structural evolution of wormlike micelles (WLMs) formed by a gemini-type cationic surfactant (12–2–12) in an aqueous solution was investigated using small-angle neutron scattering (SANS). To achieve this, a cell for observation in the velocity–velocity gradient (1–2) shear plane was constructed. SANS measurements were conducted at multiple positions across the gap of this 1–2 shear cell under various shear rates, covering both shear-thickening and shear-thinning regimes. During shear thickening at apparent shear rates of 70–200 s−1, a clear position-dependent alignment of WLMs was observed, with higher orientation near the rotor. In contrast, the intermicellar distance remained largely uniform across the gap. Neutron transmission measurements indicated no significant macroscopic concentration fluctuations. These results suggest that shear thickening is primarily driven by shear-induced micellar elongation and orientation, rather than the formation of concentration fluctuations. This study indicates the utility of SANS measurements in the 1–2 shear plane in revealing spatially resolved structural details, providing crucial insights into the 3D behavior of complex fluids under shear and complementing previous Rheo-SANS findings.

1. Introduction

Rheological measurements are fundamental to various industrial applications, including food, cosmetics, pharmaceuticals and polymeric materials, where controlling the material deformation and flow properties is crucial for product quality and performance. Complex fluids, such as colloidal dispersions and polymer solutions, exhibit diverse rheological behaviors under shear flow. A thorough understanding of these behaviors is crucial for optimizing industrial processes and enhancing product quality (Wagner & Mewis, 2021[Wagner, N. J. & Mewis, J. (2021). Theory and applications of colloidal suspension rheology. Cambridge University Press.]).

Wormlike micelles (WLMs), formed by surfactants, are notable for their nonlinear rheological behaviors (shear thinning and shear thickening) under shear flow (Dreiss & Feng, 2017[Dreiss, C. A. & Feng, Y. (2017). Wormlike micelles: advances in systems, characterization and applications. Royal Society of Chemistry.]). This distinct rheology is closely linked to hierarchical structural changes spanning the nanoscale to the mesoscale under shear stress, which include micellar elongation, orientation, branching, the formation of higher-order structures and concentration inhomogeneities. Consequently, nanoscale structural analysis under shear flow is crucial for elucidating these phenomena. Small-angle neutron scattering (SANS) is a powerful technique for these investigations. The Rheo-SANS method, which combines rheological measurements with SANS, is particularly effective in revealing the relationship between macroscopic physical properties and nanostructures.

Numerous Rheo-SANS studies on WLMs formed by cationic surfactants [e.g. cetylpyridinium chloride (CPyCl) (Mütze et al., 2014[Mütze, A., Heunemann, P. & Fischer, P. (2014). J. Rheol. 58, 1647-1672.]; Gurnon, Lopez-Barron et al., 2014[Gurnon, A. K., Lopez-Barron, C. R., Eberle, A. P., Porcar, L. & Wagner, N. J. (2014). Soft Matter 10, 2889-2898.]; Herle et al., 2007[Herle, V., Kohlbrecher, J., Pfister, B., Fischer, P. & Windhab, E. J. (2007). Phys. Rev. Lett. 99, 158302.]), cetyltri­methyl­ammonium bromide (Lutz-Bueno et al., 2013[Lutz-Bueno, V., Kohlbrecher, J. & Fischer, P. (2013). Rheol. Acta 52, 297-312.]; Takeda et al., 2011[Takeda, M., Kusano, T., Matsunaga, T., Endo, H., Shibayama, M. & Shikata, T. (2011). Langmuir 27, 1731-1738.]; Dehmoune et al., 2009[Dehmoune, J., Decruppe, J. P., Greffier, O., Xu, H. & Lindner, P. (2009). Langmuir 25, 7271-7278.]), cetyltri­methyl­ammonium tosyl­ate (Truong & Walker, 2002[Truong, M. & Walker, L. (2002). Langmuir 18, 2024-2031.]; Berret et al., 2001[Berret, J. F., Gamez-Corrales, R., Séréro, Y., Molino, F. & Lindner, P. (2001). Europhys. Lett. 54, 605-611.]; Truong & Walker, 2000[Truong, M. T. & Walker, L. M. (2000). Langmuir 16, 7991-7998.]; Berret et al., 1998[Berret, J.-F., Gamez-Corrales, R., Oberdisse, J., Walker, L. & Lindner, P. (1998). Europhys. Lett. 41, 677-682.])] with organic salts have demonstrated shear-induced orientations and phase transitions (Helgeson et al., 2010[Helgeson, M. E., Porcar, L., Lopez-Barron, C. & Wagner, N. J. (2010). Phys. Rev. Lett. 105, 084501.]). Additionally, CPyCl/sodium salicylate systems have been shown to have shear-induced structures (Helgeson et al., 2009[Helgeson, M. E., Vasquez, P. A., Kaler, E. W. & Wagner, N. J. (2009). J. Rheol. 53, 727-756. ]). Previously, we constructed a Rheo-SANS system (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]) to study the interplay between rheological behavior and nano­structures in fluid materials, such as surfactant solutions (Nemoto et al., 2024[Nemoto, F., Takabatake, F., Yamada, N. L., Takata, S. & Seto, H. (2024). J. Chem. Phys. 161, 164902.  ]; Uyama et al., 2022[Uyama, M., Takahara, A., Higaki, Y., Yamada, N. L. & Iwase, H. (2022). J. Oleo Sci. 71, 1625-1637.]; Iwase, Kawai et al., 2021[Iwase, H., Kawai, R., Yoshimura, T., Takata, S.-i. (2021). JPS Conf. Proc. 33, 011078.]; Sugahara et al., 2020[Sugahara, T., Akamatsu, M., Iwase, H., Takamatsu, Y., Sakai, K. & Sakai, H. (2020). Langmuir 36, 4695-4701.]), polymer solutions (Gupit et al., 2020[Gupit, C. I., Li, X., Maekawa, R., Hasegawa, N., Iwase, H., Takata, S. & Shibayama, M. (2020). Macromolecules 53, 1464-1473.]), inorganic suspensions (Kajiyama et al., 2020[Kajiyama, S., Iwase, H., Nakayama, M., Ichikawa, R., Yamaguchi, D., Seto, H. & Kato, T. (2020). Nanoscale 12, 11468-11479.]) and clay suspensions (Iwase, Kubota et al., 2021[Iwase, H., Kubota, M., Itoh, T., Ogura, T., Ebina, T., Ohtani, H., Kurosaka, K. & Fukushima, Y. (2021). Langmuir 37, 6435-6441.]; Iwase et al., 2018[Iwase, H., Ogura, T., Sakuma, H., Tamura, K. & Fukushima, Y. (2018). Appl. Clay Sci. 157, 24-30.]).

To understand the rheological behavior of materials under shear flow, structural characterization in three different shear planes – velocity, velocity gradient and vorticity – is required (Fig. 1[link]). Rheo-SANS enables us to observe structures in the velocity–vorticity (1–3) and velocity gradient–vorticity (2–3) shear planes. However, structures in the velocity–velocity gradient (1–2) shear plane have remained inaccessible, limiting our ability to observe the structural changes influenced by the velocity gradients.

[Figure 1]
Figure 1
Schematic diagram of a coordinate system defined by (1) velocity, (2) velocity gradient and (3) vorticity for SANS.

It has recently been reported that particle orientation and concentration fluctuations depend on the velocity gradient and are involved in peculiar rheological behavior. In particular, macroscopic concentration fluctuations, referred to as shear bands, have been observed in various systems and are a key area of interest (Huang et al., 2025[Huang, G. R., Murphy, R. P., Porcar, L., Tung, C. H., Do, C. & Chen, W. R. (2025). J. Colloid Interface Sci. 686, 1125-1134.]; Bassu et al., 2024[Bassu, G., Houston, J. E., Lara-Peña, M. A., Kriegs, H., Lettinga, M. P., Porcar, L., Scotti, A. & Laurati, M. (2024). Phys. Fluids 36, 113116. ]; Korculanin et al., 2021[Korculanin, O., Westermeier, F., Hirsemann, H., Struth, B., Hermida-Merino, D., Wagner, U. H., Donley, G. J., Rogers, S. A. & Lettinga, M. P. (2021). Phys. Fluids 33, 123104. ]; Lang, Kohlbrecher et al., 2019[Lang, C., Kohlbrecher, J., Porcar, L., Radulescu, A., Sellinghoff, K., Dhont, J. K. G. & Lettinga, M. P. (2019). Macromolecules 52, 9604-9612.]; Calabrese et al., 2016[Calabrese, M. A., Rogers, S. A., Porcar, L. & Wagner, N. J. (2016). J. Rheol. 60, 1001-1017.]; Sambasivam et al., 2015[Sambasivam, A., Sangwai, A. V. & Sureshkumar, R. (2015). Phys. Rev. Lett. 114, 158302.]; Min Kim et al., 2014[Min Kim, J., Eberle, A. P., Kate Gurnon, A., Porcar, L. & Wagner, N. J. (2014). J. Rheol. 58, 1301-1328.]; López-Barrón et al., 2014[López-Barrón, C. R., Gurnon, A. K., Eberle, A. P., Porcar, L. & Wagner, N. J. (2014). Phys. Rev. E 89, 042301.]; Fielding, 2014[Fielding, S. M. (2014). Rep. Prog. Phys. 77, 102601.]; Kunita et al., 2012[Kunita, I., Sato, K., Tanaka, Y., Takikawa, Y., Orihara, H. & Nakagaki, T. (2012). Phys. Rev. Lett. 109, 248303.]; Helgeson et al., 2010[Helgeson, M. E., Porcar, L., Lopez-Barron, C. & Wagner, N. J. (2010). Phys. Rev. Lett. 105, 084501.]; Furukawa & Tanaka, 2009[Furukawa, A. & Tanaka, H. (2009). Nat. Mater. 8, 601-609.]; Miller & Rothstein, 2007[Miller, E. & Rothstein, J. P. (2007). J. Non-Newt. Fluid Mech. 143, 22-37.]; Hu & Lips, 2005[Hu, Y. T. & Lips, A. (2005). J. Rheol. 49, 1001-1027.]; Lettinga & Dhont, 2004[Lettinga, M. & Dhont, J. (2004). J. Phys. Condens. Matter 16, S3929-S3939.]; Salmon et al., 2003[Salmon, J.-B., Colin, A., Manneville, S. & Molino, F. (2003). Phys. Rev. Lett. 90, 228303.]; Decruppe et al., 2001[Decruppe, J. P., Lerouge, S. & Berret, J. F. (2001). Phys. Rev. E 63, 022501.]; Liu & Pine, 1996[Liu, C. & Pine, D. J. (1996). Phys. Rev. Lett. 77, 2121-2124.]). Consequently, experimental setups capable of observing structures in the 1–2 shear plane are essential for a complete understanding of material behavior under shear flow (Lang, Porcar et al., 2019[Lang, C., Porcar, L., Kriegs, H. & Lettinga, M. P. (2019). J. Phys. D Appl. Phys. 52, 074003.]; Eberle & Porcar, 2012[Eberle, A. P. R. & Porcar, L. (2012). Curr. Opin. Colloid Interface Sci. 17, 33-43.]). While Rheo-SANS systems have been developed in multiple neutron scattering facilities, cells specifically designed to observe structure in the 1–2 shear plane remain relatively rare (Velichko et al., 2019[Velichko, E., Tian, B., Nikolaeva, T., Koning, J., van Duynhoven, J. & Bouwman, W. G. (2019). Colloids Surf. A Physicochem. Eng. Asp. 566, 21-28.]; Gurnon, Godfrin et al., 2014[Gurnon, A. K., Godfrin, P. D., Wagner, N. J., Eberle, A. P., Butler, P. & Porcar, L. (2014). J. Vis. Exp. 84, e51068.]; Eberle & Porcar, 2012[Eberle, A. P. R. & Porcar, L. (2012). Curr. Opin. Colloid Interface Sci. 17, 33-43.]; Liberatore et al., 2006[Liberatore, M. W., Nettesheim, F., Wagner, N. J. & Porcar, L. (2006). Phys. Rev. E 73, 020504.]; Van der Werff et al., 1990[Van der Werff, J. C., Ackerson, B. J., May, R. P. & de Kruif, C. G. (1990). Physica A 165, 375-398.]).

In this study, a 1–2 shear cell was utilized to investigate the shear-induced structural behavior of the gemini-type cationic surfactant 12–2–12, with a particular focus on elucidating spatially resolved structural behaviors in the 1–2 shear plane. The aim of this study is to provide new insights into the 3D response to shear flow by complementing our previous Rheo-SANS investigations (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]), which characterized structures in the 1–3 and 2–3 shear planes. The 12–2–12 system, which has been observed to form WLMs and exhibit both shear-thinning and shear-thickening behaviors, serves as an ideal model to explore the structural origins of these complex rheological phenomena through spatially resolved SANS measurements in the 1–2 shear plane.

2. Materials and methods

2.1. Samples

The gemini surfactant 12–2–12 consists of two cationic amphiphilic units [C12H25N+(CH3)2Br] connected by a spacer, as shown in Fig. 2[link]. The compound was synthesized following a previously reported procedure (Esumi et al., 1996[Esumi, K., Taguma, K. & Koide, Y. (1996). Langmuir 12, 4039-4041.]; Zana & Talmon, 1993[Zana, R. & Talmon, Y. (1993). Nature 362, 228-230.]). A 12–2–12 solution was prepared in D2O with a volume fraction (φ) of 0.0135. In salt-free solutions, 12–2–12 undergoes micellar transitions, forming spherical micelles at low concentrations and transitioning to rod-like and wormlike micelles as concentration increases (In et al., 2010[In, M., Bendjeriou, B., Noirez, L. & Grillo, I. (2010). Langmuir 26, 10411-10414.]; Oda et al., 1997[Oda, R., Panizza, P., Schmutz, M. & Lequeux, F. (1997). Langmuir 13, 6407-6412.]). Further increases in concentration extend micellar lengths. Previous Rheo-SANS studies have demonstrated shear-thinning and -thickening behavior in this system at 30°C (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]).

[Figure 2]
Figure 2
Chemical structure of the gemini-type cationic surfactant (12–2–12).

2.2. Rheological measurements

Rheological measurements were conducted using an MCR302 stress-controlled rheometer (Anton Paar, Graz, Austria). A cone-plate geometry (e.g. 50 mm diameter, 1° angle, 0.100 mm truncation) was utilized. The shear rate was varied from 10 to 900 s1, while the sample temperature was maintained at 30°C using a Peltier temperature controller.

2.3. SANS

SANS measurements were conducted using the small- and wide-angle neutron scattering instrument (TAIKAN) installed at the Materials and Life Science Experimental Facility (MLF) of the Japan Proton Accelerator Research Complex (J-PARC), Tokai, Japan (Takata et al., 2015[Takata, S., Suzuki, J., Shinohara, T., Oku, T., Tominaga, T., Ohishi, K., Iwase, H., Nakatani, T., Inamura, Y., Ito, T., Suzuya, K., Aizawa, K., Arai, M., Otomo, T. & Sugiyama, M. (2015). JPS Conf. Proc. 8, 036020.]). TAIKAN has four detector banks (small-angle, middle-angle, high-angle and backward). The small-angle detector bank, located 5.65 m from the sample position, had a large effective area of ∼2.2 × 2.3 mm, excluding a 64 × 64 mm central area. A low-efficiency N2 monitor positioned behind the center of the small-angle detector bank measured transmission simultaneously with scattering. Using white neutron beams with a wavelength (λ) range of 1.0–7.5 Å based on the time-of-flight (TOF) method, a Q range from 5 × 10−3 to 17 Å−1 could be covered in a single measurement [the magnitude of the scattering vector Q is defined as Q = 4π sin(θ/2)/λ, with θ the scattering angle]. However, in this study, the observed Q range was constrained to 0.01 < Q < 1 Å−1 due to instrumental background at lower Q and low accuracy of the scattering intensity from the sample at higher Q.

The 12–2–12 solution was investigated under shear flow using a recently constructed 1–2 shear cell. As illustrated in Fig. 3[link], the 1–2 shear cell allows for structural observation in the velocity–velocity gradient (1–2) shear plane. The cell was fabricated from a Ti–6 Al–4 V (Ti-64) titanium alloy, a material selected for its exceptional corrosion resistance. The system design was informed by prior successful implementations in the field (Gurnon, Godfrin et al., 2014[Gurnon, A. K., Godfrin, P. D., Wagner, N. J., Eberle, A. P., Butler, P. & Porcar, L. (2014). J. Vis. Exp. 84, e51068.]; Liberatore et al., 2006[Liberatore, M. W., Nettesheim, F., Wagner, N. J. & Porcar, L. (2006). Phys. Rev. E 73, 020504.]), notably incorporating a liquid seal concept. The rotor, with an inner diameter of 50.0 mm, featured a flow channel width (gap) of 1.5 mm and a path length (sample thickness) of 6 mm. Shear flow was generated using a brushless motor (BLM015HK-5CS, Oriental Motor Co. Ltd) operated with a 1:5 gear ratio. Brushless motors are known for their high operational speeds and low levels of vibration. These motors are characterized by a specified rotation speed range of 100–3000 r min−1, enabling access to a shear rate (Mathematical equation) range of 40–1000 s−1 for the current experiments. The apparent shear rate was related to the motor angular velocity (Ω) via equation (1[link]), based on assumptions consistent with double-cylindrical Couette cell configurations (JIS Z 8803, 1991[JIS Z 8803 (1991). Viscosity of liquid - methods of measurement. Japanese Standards Association.]):

Mathematical equation

where Ri (mm) is the rotor radius and tf (mm) is the flow channel width. Note that the shear cell is designed with the potential to achieve higher shear rates, up to approximately 5000 s−1, through the implementation of an alternative motor or gear configuration. However, such capabilities were not explored in the present study. The 1–2 shear cell cannot directly measure shear stress because it lacks torque sensors. Therefore, Rheo-SANS measurements are necessary to investigate the relationship between rheological behavior and nanostructure. The study employed martensitic stainless-steel bearings (TBN-6H, NSK Micro Precision Co. Ltd, Tokyo, Japan). Liquid seals, fabricated from Viton O-rings, were strategically positioned between the bearings and the donut-shaped quartz glass windows (3 mm thickness). The sample was inserted into the cell by means of a tubing pump (Peristaltic Pump, Atto Corporation, Tokyo, Japan).

[Figure 3]
Figure 3
(a), (b) Schematic drawing of the 1–2 shear cell. (1) Sample inlets, (2) a brushless motor, (3) donut-shaped quartz glass windows (thickness of 3 mm), (4) a rotor, (5) a bearing, (6) O-rings. Photographs of (c) the rotor and neutron-irradiated area (dashed red window), (d) the Cd aperture, and (e) the Cd shield.

A sample aperture (w0.5 × h10 mm) was positioned on an XZ movable bench 150 mm upstream of the sample [Fig. 3[link](d)]. The XZ movable bench is usually used to adjust the standard sample aperture. Additionally, a manually operated Cd shield (w1.5 × h15 mm) was placed directly in front of the sample window [Fig. 3[link](e)]. Parasitic scattering is influenced by the sample aperture size, first aperture size and collimation length (CL), accounting for beam divergence. For measurements on the TAIKAN instrument, the first aperture size was set to 8 × 23 mm, with a CL of 7.55 m (for a 1–2 shear SANS measurement). The aperture size was similar to that previously used in Rheo-SANS measurement (tangential position). Due to the significantly reduced sample width required for these 1–2 plane measurements, a proven optical system designed for a 0.5 mm beam width was initially employed to ensure reliable performance.

The shear cell operation and data acquisition on the TAIKAN instrument were managed using the instrument software framework IROHA2 developed by J-PARC MLF (Hasemi et al., 2024[Hasemi, H., Takahashi, R., Yamauchi, Y., Ishikado, M., Ohira-Kawamura, S. & Komine, R. (2024). JPS Conf. Proc. 41, 011003.]). This framework supports automated measurements and system status monitoring, ensuring reliable and reproducible data collection. The apparent shear rates were incrementally increased from 0 to 900 s−1, with simultaneous SANS and transmission measurements at each step. The exposure time for each measurement was 15 min, and the sample temperature was maintained at 30°C using temperature-controlled circulating water.

TOF data were corrected for wavelength dependence using measurements from an empty cell. The scattering intensity I(Q) was obtained by sector averaging data along the vertical and horizontal directions with an azimuthal angle of ±15°. Parasitic scattering from the cell wall was removed analytically by pixel masking. Data reduction and analysis were performed using the UTSUSEMI software package (Inamura et al., 2013[Inamura, Y., Nakatani, T., Suzuki, J. & Otomo, T. (2013). J. Phys. Soc. Jpn 82, SA031.]). Absolute scattering intensity calibration was performed in the following manner. First, SANS measurements of a glassy carbon secondary standard (Zhang et al., 2010[Zhang, F., Ilavsky, J., Long, G. G., Quintana, J. P., Allen, A. J. & Jemian, P. R. (2010). Metall. Mater. Trans. A 41, 1151-1158.]) and the 12–2–12 solution inserted into a quartz cell were conducted using standard pinhole collimation on TAIKAN (10 mm-diameter sample aperture) to obtain their absolute scattering intensities per unit sample volume. Subsequently, SANS measurements of the 12–2–12 solution in a stationary state were performed at each of the four positions within the 1–2 shear cell using the 0.5 mm-wide neutron beam. The Q dependence of these SANS profiles was found to be identical to that obtained with the standard pinhole collimation. Therefore, by comparing these datasets, the data from the 1–2 shear cell measurements were scaled to absolute units. This procedure also allowed for an estimation of the effective irradiated area at each measurement position within the 1–2 shear cell.

3. Results and discussion

3.1. Shear rate dependence

3.1.1. Rheological measurements

Fig. 4[link] illustrates the Mathematical equation dependence of the viscosity of 12–2–12 in an aqueous solution at φ = 0.0135 at 30°C, as measured using a cone-plate geometry. At Mathematical equation < 40 s−1, the apparent viscosity was relatively low. The viscosity increased within the range of 40–150 s−1, showing shear-thickening behavior. At higher shear rates (Mathematical equation > 300 s−1), a decrease in viscosity, known as shear thinning, was observed. These rheological behaviors observed for 12–2–12 in an aqueous solution with the cone-plate geometry are qualitatively consistent with previously reported results, including our previous Rheo-SANS results (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]). While the general trend is similar, quantitative comparisons of critical shear rates or the magnitude of viscosity changes should consider the differences in measurement geometries (cone-plate here versus Couette or parallel-plate in other studies or in SANS cells). For instance, slight variations in the onset of shear thickening or the rate of viscosity increase compared with measurements in different geometries might be expected. Consequently, these rheological characteristics, determined under well defined shear conditions, provide the necessary context for interpreting the SANS results obtained under shear.

[Figure 4]
Figure 4
Shear rate dependence of the viscosity of 12–2–12 in an aqueous solution (D2O) at 30°C. The volume fraction (φ) is 0.0135.
3.1.2. SANS

Fig. 5[link] shows 2D SANS profiles for 12–2–12 in an aqueous solution at φ = 0.0135 under varying shear rates ranging from 0 to 900 s−1. The 2D SANS profiles were approximately isotropic at lower shear rates (Mathematical equation < 40 s−1), indicating a random orientation of WLMs. In contrast, at relatively high shear rates (Mathematical equation > 100 s−1), the 2D SANS profiles were anisotropic, with strong scattering observed perpendicular to the shear flow. This anisotropy is attributed to the orientation of WLMs in the flow direction, consistent with previous Rheo-SANS studies (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]).

[Figure 5]
Figure 5
2D SANS profiles I(Qx, Qy) for 12–2–12 in an aqueous solution at a volume fraction (φ) of 0.0135. (a) Mathematical equation = 0 s−1, (b) Mathematical equation = 40 s−1, (c) Mathematical equation = 60 s−1, (d) Mathematical equation = 80 s−1, (e) Mathematical equation = 100 s−1, (f) Mathematical equation = 200 s−1, (g) Mathematical equation = 600 s−1, (h) Mathematical equation = 900 s−1.

To further analyze these structural changes, 1D SANS profiles were obtained by sector averaging the 2D SANS profile. Fig. 6[link] shows the sector-averaged SANS profiles for 12–2–12 in an aqueous solution. The horizontal [Ix(Q)] and vertical [Iy(Q)] SANS profiles were perpendicular and parallel to the flow, respectively. Under stationary conditions (Mathematical equation = 0 s−1), Ix(Q) and Iy(Q) were identical within experimental uncertainty (see inset of Fig. 6[link]), indicating isotropic scattering. All SANS profiles for 12–2–12 in an aqueous solution exhibited broad peaks in the Q range of 0.015–0.04 Å−1, arising from electrostatic repulsions between the surface charges of the micelles. The peak position indicates the intermicellar distance. In the horizontal direction, the SANS profiles Ix(Q) in the Q range of 0.01 < Q < 0.2 Å−1 showed slight changes near Mathematical equation = 200 s−1, with an increase in peak intensity. In contrast, the vertical SANS intensity Iy(Q) decreased, and its peak profile broadened within the same shear rate range. The observed changes in the horizontal and vertical SANS profiles with increasing shear rates were attributed to the orientation of the WLMs along the flow direction (corresponding to the vertical direction in the profiles). A more quantitative analysis of this orientation is provided by the alignment factor discussed subsequently.

[Figure 6]
Figure 6
Shear rate dependence of the sector-averaged SANS profiles for 12–2–12 in an aqueous solution at a volume fraction (φ) of 0.0135. (a) Horizontal and (b) vertical SANS profiles which were, respectively, perpendicular and parallel to the flow. The 2D SANS profiles (in Fig. 5[link]) were sector-averaged with the azimuthal angle of ±15°. The inset of Fig. 6[link](a) presents the sector-average SANS profiles for 12–2–12 in an aqueous solution under stationary states (Mathematical equation = 0 s−1).

To quantitatively evaluate the degree of orientation of the WLMs, the alignment factor (Mathematical equation) was estimated using the following equation (Angelico et al., 2010[Angelico, R., Rossi, C. O., Ambrosone, L., Palazzo, G., Mortensen, K. & Olsson, U. (2010). Phys. Chem. Chem. Phys. 12, 8856-8862.]; Borse et al., 2004[Borse, M., Sharma, V., Aswal, V. K., Pokhriyal, N. K., Joshi, J. V., Goyal, P. S. & Devi, S. (2004). Phys. Chem. Chem. Phys. 6, 3508-3514.]; Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]):

Mathematical equation

where Ix(Q) and Iy(Q) are the SANS intensities perpendicular and parallel to the flow, respectively. The Mathematical equation values were evaluated using Qmin = 0.015 Å−1 and Qmax = 0.06 Å−1.

Fig. 7[link] shows the shear rate dependence of Mathematical equation. At Mathematical equation < 60 s−1, the Mathematical equation values remained consistently low, indicating that the WLMs were randomly distributed throughout the solution. When the shear rate was increased from 60 to 200 s−1, the Mathematical equation values gradually increased from 0 to ∼0.55. This behavior was observed near the shear thickening, indicating that the micelle orientation and shear thickening occurred in the same shear rate region. Furthermore, in the range of Mathematical equation > 200 s−1, where shear thinning was observed, the Mathematical equation values remained nearly constant. The orientational behaviors of the WLMs formed by the 12–2–12 solution in the Mathematical equation = 70–200 s−1 (shear-thickening) and >200 s−1 (shear-thinning) regions were consistent with the results obtained in previous Rheo-SANS studies on the same 12–2–12 in an aqueous solution at φ = 0.0135, suggesting that the 1–2 shear cell is capable of applying the intended shear field and detecting the expected structural response.

[Figure 7]
Figure 7
Shear rate dependence of the alignment factor Mathematical equation. The Mathematical equation value was obtained by the 1D SANS profiles.

3.2. Position (velocity gradient) dependence

3.2.1. SANS

To investigate the position-dependent structural behavior of WLMs formed by 12–2–12 under shear flow within the 1–2 shear plane, SANS measurements were performed at four different positions across the 1.5 mm gap of the flow channel. The centers of these positions were separated by approximately 0.23 mm. Fig. 8[link] shows the position dependence of the 2D SANS profiles for the 12–2–12 solution at φ = 0.0135 at Mathematical equation = (a) 0, (b) 100 and (c) 900 s−1. The 2D SANS for D2O is also shown in Fig. 8[link](d), indicating that the parasitic scattering from the 1–2 shear cell wall has been removed. Note that the profiles for position 2 are identical to those shown in Fig. 5[link].

[Figure 8]
Figure 8
Position dependence of 2D SANS profiles I(Qx, Qy) for 12–2–12 in an aqueous solution at a volume fraction (φ) of 0.0135 at (a) Mathematical equation = 0 s−1, (b) Mathematical equation = 100 s−1, (c) Mathematical equation = 900 s−1, (d) D2O (Mathematical equation = 0 s−1).

Fig. 9[link] shows the position dependence of circular-averaged SANS profiles for the 12–2–12 in an aqueous solution at Mathematical equation = 0 s−1, normalized by sample thickness. While the Q dependence of the SANS profiles demonstrated consistency across all measurement positions, significant variations were observed in SANS intensities. Specifically, the SANS intensity at position 4 was very weak compared with those at other measurement positions. When the intensity was normalized to position 2 (as indicated in the figure legend), the resulting scattering intensity values, which reflect the irradiated area, revealed that the SANS profiles at position 4 are equivalent to 41% of the intensity at position 2. This was attributed to the proximity of position 4 to the Cd shielding in front of the 1–2 shear cell. The inset of Fig. 9[link] illustrates the expected gradient of the shear flow with respect to the beam irradiation positions within the 1–2 shear cell gap, as determined from the results in Fig. 9[link]. The sector-averaged SANS profiles of the 12–2–12 solution are shown in Fig. 10[link]. At Mathematical equation = 0 s−1, all SANS profiles were isotropic, and Ix(Q) and Iy(Q) were identical, irrespective of the measurement position within the gap. This finding is consistent with a randomly oriented micellar structure under stationary conditions.

[Figure 9]
Figure 9
Position dependence of circular-averaged SANS profiles for 12–2–12 in an aqueous solution at Mathematical equation = 0 s−1. The SANS profiles were normalized by sample thickness. The inset illustrates the expected gradient of the shear flow with respect to the beam irradiation positions within the 1–2 shear cell gap. The beam irradiation areas are shifted systematically in the vertical direction to avoid overlap. The distance between the centers of each position in the horizontal direction is 0.23 mm.
[Figure 10]
Figure 10
Position dependence of the sector-averaged SANS profiles for 12–2–12 in an aqueous solution at a volume fraction (φ) of 0.013 at (a) Mathematical equation = 0 s−1, (b) Mathematical equation = 100 s−1, (c) Mathematical equation = 900 s−1. The 2D SANS profiles (in Fig. 8[link]) were sector-averaged with the azimuthal angle of ±15°.

In contrast, at Mathematical equation = 100 s−1 [Figs. 8[link](b) and 10[link](b)], the SANS profiles exhibited a clear position dependence. The SANS intensity showed strong anisotropy near the rotor and gradually became isotropic towards the periphery of the flow channel. This behavior is attributed to the velocity gradient across the gap, where the local shear rate is highest near the moving rotor (position 1) and decreases towards the stationary wall (position 4). The observed WLM alignment varies according to this local shear rate gradient. This spatially resolved observation of varying structural anisotropy within the shear-thickening regime is a key finding enabled by the 1–2 shear plane measurement geometry. At Mathematical equation = 900 s−1 [Fig. 8[link](c)], well into the shear-thinning region, the SANS profiles at all positions exhibited strong anisotropy, indicating a highly oriented state of WLMs throughout the gap. As shown in the inset of Fig. 10[link](c), the sector-averaged SANS profiles were nearly identical across all positions, suggesting no position dependence in the Mathematical equation range where shear thinning was observed.

Fig. 11[link](a) shows the shear rate dependence of the Af value, evaluated using equation (2[link]) at each of the four measurement positions. The same behavior as observed in position 2 was seen across all positions. Consistent with the 2D patterns (Fig. 8[link]), at low shear rates (Mathematical equation < 40 s−1), the Af values were low and similar for all positions, reflecting the random orientation of WLMs throughout the gap. As the shear rate increased into the shear-thickening regime (Mathematical equation ≃ 40–200 s−1), Af increased significantly, indicative of developing micellar alignment. Crucially, within this regime, Af also exhibited a strong dependence on position: higher Af values were observed closer to the moving rotor (position 1), and these values progressively decreased towards the stationary wall (position 4). This observation clearly demonstrates that the degree of shear-induced micellar orientation is non-uniform across the gap during shear thickening. In contrast, during the shear-thinning regime (Mathematical equation > 200 s−1), Af values were high and nearly constant across all positions, indicating a uniformly highly oriented state of micelles throughout the gap.

[Figure 11]
Figure 11
(a) Shear rate dependence of the alignment factor Mathematical equation and (b) peak position (Qpeak) for 12–2–12 in an aqueous solution at four positions. The Mathematical equation values were obtained by equation (2).

The position dependence of the characteristic intermicellar distance was further investigated by evaluating the peak position (Qpeak) from the Ix(Q) and Iy(Q) profiles at each of the four positions across the gap as a function of shear rate [Fig. 11[link](b)]. As can be observed, the Qpeak values remain remarkably constant across the four measurement positions. Furthermore, the Qpeak values themselves remain largely constant with increasing shear rate, with only a slight tendency to decrease (indicating a slight increase in intermicellar spacing) in the high shear rate shear-thinning regime. The lack of notable position dependence of Qpeak, even in the shear-thickening regime where Af exhibits substantial spatial variation [Fig. 11[link](a)], indicates that the average intermicellar spacing is relatively consistent across the gap. This indicates that local changes in micelle concentration are not the predominant cause of position-dependent alignment.

3.2.2. Transmission

Macroscopic concentration fluctuations, such as shear bands, are often associated with complex rheological behaviors in surfactant solutions (Liu & Pine, 1996[Liu, C. & Pine, D. J. (1996). Phys. Rev. Lett. 77, 2121-2124.]; Decruppe et al., 2001[Decruppe, J. P., Lerouge, S. & Berret, J. F. (2001). Phys. Rev. E 63, 022501.]; Salmon et al., 2003[Salmon, J.-B., Colin, A., Manneville, S. & Molino, F. (2003). Phys. Rev. Lett. 90, 228303.]; Miller & Rothstein, 2007[Miller, E. & Rothstein, J. P. (2007). J. Non-Newt. Fluid Mech. 143, 22-37.]; Sambasivam et al., 2015[Sambasivam, A., Sangwai, A. V. & Sureshkumar, R. (2015). Phys. Rev. Lett. 114, 158302.]). To investigate the possibility of macroscopic concentration fluctuations in the 12–2–12 solution, neutron transmission through the sample was measured. The transmission of the sample solution is predominantly determined by the number of hydrogen atoms, as hydrogen has the largest scattering cross section among all atoms in the 12–2–12 solution. It is known that the H2O ratio in D2O/H2O mixtures can be evaluated from transmission values (May et al., 1982[May, R. P., Ibel, K. & Haas, J. (1982). J. Appl. Cryst. 15, 15-19.]). At TAIKAN, the transmission is detected simultaneously with the scattering.

Fig. 12[link](a) shows the wavelength dependence of the transmission measured simultaneously with the SANS at a shear rate of 100 s−1. Consequently, the transmission value was directly related to the surfactant volume fraction. The wavelength dependence of the transmission of the 12–2–12 solutions at φ = 0, 0.0135 and 0.027 is shown in Fig. 12[link]. The transmission gradually decreased as the sample volume fraction increased. A clear position dependence was observed in the SANS profile for the 12–2–12 solutions, whereas the transmission remained nearly uniform. To further investigate this with improved statistical precision, the mean transmission was calculated over the wavelength range of 3 to 6 Å. Fig. 12[link](b) shows the shear rate dependence of the mean value of transmission of the 12–2–12 solutions at each of the four positions across the gap. These mean transmission values are largely uniform across all four positions for both the quiescent state and under shear at Mathematical equation = 100 s−1, within experimental uncertainty. No significant or systematic variation in transmission with position was detected that would suggest the formation of distinct shear bands with notably different surfactant concentrations.

[Figure 12]
Figure 12
Wavelength dependence of neutron transmission for 12–2–12 in an aqueous solution at four positions. The transmission and SANS were detected simultaneously during the measurement.

The transmission results were consistent with the largely uniform Qpeak values across the gap [Fig. 11[link](b)], suggesting that the shear thickening observed around 100 s−1 is primarily associated with the shear-induced changes in wormlike micellar structure (elongation and orientation) and the spatial variation of the WLMs, as revealed by SANS, rather than being dominated by macroscopic concentration fluctuations. However, it is important to recognize that, at the relatively low surfactant concentration used in this study (φ = 0.0135), the sensitivity of transmission measurements to subtle concentration changes is limited. While our data do not support strong macroscopic shear banding, the possibility of more subtle concentration variations cannot be entirely excluded on the basis of the transmission data alone. Nonetheless, the consistency of Qpeak across the gap further supports the idea that large variations in local micelle density are unlikely.

Therefore, considering the SANS results together with the transmission data, the distinct position dependence of WLM orientation [Fig. 11[link](a)], in the absence of clear macroscopic concentration fluctuations, suggests that the structural response of the micelles to the local shear rate within the velocity gradient is the primary factor driving the observed spatial inhomogeneity during shear thickening. The micellar structuring process, likely involving elongation and subsequent orientation, appears to be highly sensitive to these local shear conditions, leading to the observed non-uniform alignment across the gap.

4. Conclusions

In this study, we investigated the shear-induced structural evolution of WLMs formed by the gemini-type cationic surfactant 12–2–12 in an aqueous solution, focusing specifically on the velocity–velocity gradient (1–2) shear plane. By performing SANS measurements at four different positions across the flow channel of a 1–2 shear cell, we successfully characterized the spatial dependence of WLM structures under varying shear rates, encompassing both shear-thickening and shear-thinning regimes.

The key findings demonstrate that during the shear-thickening process (Mathematical equation ≃ 70–200 s−1) a significant position-dependent alignment of WLMs occurs across the gap, with higher orientation near the moving rotor, indicating a strong sensitivity of micellar structural response to the local shear rate. Despite this pronounced spatial variation in alignment, the intermicellar distance, inferred from the SANS peak position (Qpeak), remained largely uniform across the gap, even under conditions of strong orientational gradients (e.g. at Mathematical equation ≃ 100 s−1). Furthermore, simultaneous neutron transmission measurements, corroborated by the consistent Qpeak values, did not reveal significant macroscopic concentration fluctuations within the shear-thickening regime. These collective results strongly suggest that the observed shear thickening is primarily driven by shear-induced micellar elongation and orientation, and the spatial inhomogeneity of the WLMs, rather than by distinct shear banding involving significant concentration differences.

These observations in the 1–2 plane are broadly consistent with and complementary to our previous Rheo-SANS investigations of the same 12–2–12 system in the 1–3 and 2–3 shear planes (Iwase et al., 2019[Iwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357-366.]). However, this new study uniquely provides direct experimental evidence of the spatial inhomogeneity of micellar alignment within the velocity gradient, offering crucial insights into how local shear conditions govern the structural transitions underlying macroscopic rheological behavior. The capacity to address these position-dependent structural intricacies underscores the efficacy of SANS measurements in the 1–2 shear plane, facilitating a more comprehensive, 3D understanding of complex fluids under shear. This approach has significant potential for elucidating the intricate relationships between nanostructure, local flow kinematics and macroscopic rheology in a wide range of soft matter systems. Future investigations could extend this methodology to other shear-thickening systems or explore shear-thinning phenomena in more concentrated micellar solutions, further leveraging the capabilities of SANS measurements in the 1–2 shear plane to unravel complex flow behaviors.

Acknowledgements

The authors thank Professor Mikio Kataoka (NAIST/CROSS) for the helpful discussions on the construction of the 1–2 shear cell and Misaki Ueda (CROSS) for support in rheological measurements. The rheological measurements were carried out in the laboratory of the J-PARC research building. The SANS experiment at the Materials and Life Science Experimental Facility of J-PARC was performed under a user program (proposal No. 2023C0002). We thank Enago for the English language review.

Funding information

This work was supported by JSPS KAKENHI (grant Nos. JP19K05413, JP22K05058).

References

First citationAngelico, R., Rossi, C. O., Ambrosone, L., Palazzo, G., Mortensen, K. & Olsson, U. (2010). Phys. Chem. Chem. Phys. 12, 8856–8862.  CrossRef CAS PubMed Google Scholar
First citationBassu, G., Houston, J. E., Lara-Peña, M. A., Kriegs, H., Lettinga, M. P., Porcar, L., Scotti, A. & Laurati, M. (2024). Phys. Fluids 36, 113116.   Google Scholar
First citationBerret, J.-F., Gamez-Corrales, R., Oberdisse, J., Walker, L. & Lindner, P. (1998). Europhys. Lett. 41, 677–682.  CrossRef CAS Google Scholar
First citationBerret, J. F., Gamez-Corrales, R., Séréro, Y., Molino, F. & Lindner, P. (2001). Europhys. Lett. 54, 605–611.  CrossRef CAS Google Scholar
First citationBorse, M., Sharma, V., Aswal, V. K., Pokhriyal, N. K., Joshi, J. V., Goyal, P. S. & Devi, S. (2004). Phys. Chem. Chem. Phys. 6, 3508–3514.  CrossRef CAS Google Scholar
First citationCalabrese, M. A., Rogers, S. A., Porcar, L. & Wagner, N. J. (2016). J. Rheol. 60, 1001–1017.  CrossRef CAS Google Scholar
First citationDecruppe, J. P., Lerouge, S. & Berret, J. F. (2001). Phys. Rev. E 63, 022501.  CrossRef Google Scholar
First citationDehmoune, J., Decruppe, J. P., Greffier, O., Xu, H. & Lindner, P. (2009). Langmuir 25, 7271–7278.  CrossRef PubMed CAS Google Scholar
First citationDreiss, C. A. & Feng, Y. (2017). Wormlike micelles: advances in systems, characterization and applications. Royal Society of Chemistry.  Google Scholar
First citationEberle, A. P. R. & Porcar, L. (2012). Curr. Opin. Colloid Interface Sci. 17, 33–43.  Web of Science CrossRef CAS Google Scholar
First citationEsumi, K., Taguma, K. & Koide, Y. (1996). Langmuir 12, 4039–4041.  CrossRef CAS Google Scholar
First citationFielding, S. M. (2014). Rep. Prog. Phys. 77, 102601.  CrossRef PubMed Google Scholar
First citationFurukawa, A. & Tanaka, H. (2009). Nat. Mater. 8, 601–609.  Web of Science CrossRef PubMed CAS Google Scholar
First citationGupit, C. I., Li, X., Maekawa, R., Hasegawa, N., Iwase, H., Takata, S. & Shibayama, M. (2020). Macromolecules 53, 1464–1473.  CrossRef CAS Google Scholar
First citationGurnon, A. K., Godfrin, P. D., Wagner, N. J., Eberle, A. P., Butler, P. & Porcar, L. (2014). J. Vis. Exp. 84, e51068.  Google Scholar
First citationGurnon, A. K., Lopez-Barron, C. R., Eberle, A. P., Porcar, L. & Wagner, N. J. (2014). Soft Matter 10, 2889–2898.  CrossRef CAS PubMed Google Scholar
First citationHasemi, H., Takahashi, R., Yamauchi, Y., Ishikado, M., Ohira-Kawamura, S. & Komine, R. (2024). JPS Conf. Proc. 41, 011003.  Google Scholar
First citationHelgeson, M. E., Porcar, L., Lopez-Barron, C. & Wagner, N. J. (2010). Phys. Rev. Lett. 105, 084501.  CrossRef PubMed Google Scholar
First citationHelgeson, M. E., Vasquez, P. A., Kaler, E. W. & Wagner, N. J. (2009). J. Rheol. 53, 727–756.   CrossRef CAS Google Scholar
First citationHerle, V., Kohlbrecher, J., Pfister, B., Fischer, P. & Windhab, E. J. (2007). Phys. Rev. Lett. 99, 158302.  CrossRef PubMed Google Scholar
First citationHu, Y. T. & Lips, A. (2005). J. Rheol. 49, 1001–1027.  CrossRef CAS Google Scholar
First citationHuang, G. R., Murphy, R. P., Porcar, L., Tung, C. H., Do, C. & Chen, W. R. (2025). J. Colloid Interface Sci. 686, 1125–1134.  CrossRef CAS PubMed Google Scholar
First citationIn, M., Bendjeriou, B., Noirez, L. & Grillo, I. (2010). Langmuir 26, 10411–10414.  CrossRef CAS PubMed Google Scholar
First citationInamura, Y., Nakatani, T., Suzuki, J. & Otomo, T. (2013). J. Phys. Soc. Jpn 82, SA031.  CrossRef Google Scholar
First citationIwase, H., Kawai, R., Morishima, K., Takata, S., Yoshimura, T. & Shibayama, M. (2019). J. Colloid Interface Sci. 538, 357–366.  Web of Science CrossRef CAS PubMed Google Scholar
First citationIwase, H., Kawai, R., Yoshimura, T., Takata, S.-i. (2021). JPS Conf. Proc. 33, 011078.  Google Scholar
First citationIwase, H., Kubota, M., Itoh, T., Ogura, T., Ebina, T., Ohtani, H., Kurosaka, K. & Fukushima, Y. (2021). Langmuir 37, 6435–6441.  CrossRef CAS PubMed Google Scholar
First citationIwase, H., Ogura, T., Sakuma, H., Tamura, K. & Fukushima, Y. (2018). Appl. Clay Sci. 157, 24–30.  Web of Science CrossRef CAS Google Scholar
First citationJIS Z 8803 (1991). Viscosity of liquid – methods of measurement. Japanese Standards Association.  Google Scholar
First citationKajiyama, S., Iwase, H., Nakayama, M., Ichikawa, R., Yamaguchi, D., Seto, H. & Kato, T. (2020). Nanoscale 12, 11468–11479.  CrossRef CAS PubMed Google Scholar
First citationKorculanin, O., Westermeier, F., Hirsemann, H., Struth, B., Hermida-Merino, D., Wagner, U. H., Donley, G. J., Rogers, S. A. & Lettinga, M. P. (2021). Phys. Fluids 33, 123104.   Google Scholar
First citationKunita, I., Sato, K., Tanaka, Y., Takikawa, Y., Orihara, H. & Nakagaki, T. (2012). Phys. Rev. Lett. 109, 248303.  CrossRef PubMed Google Scholar
First citationLang, C., Kohlbrecher, J., Porcar, L., Radulescu, A., Sellinghoff, K., Dhont, J. K. G. & Lettinga, M. P. (2019). Macromolecules 52, 9604–9612.  CrossRef CAS Google Scholar
First citationLang, C., Porcar, L., Kriegs, H. & Lettinga, M. P. (2019). J. Phys. D Appl. Phys. 52, 074003.  CrossRef Google Scholar
First citationLettinga, M. & Dhont, J. (2004). J. Phys. Condens. Matter 16, S3929–S3939.  CrossRef CAS Google Scholar
First citationLiberatore, M. W., Nettesheim, F., Wagner, N. J. & Porcar, L. (2006). Phys. Rev. E 73, 020504.  CrossRef Google Scholar
First citationLiu, C. & Pine, D. J. (1996). Phys. Rev. Lett. 77, 2121–2124.  CrossRef PubMed CAS Google Scholar
First citationLópez-Barrón, C. R., Gurnon, A. K., Eberle, A. P., Porcar, L. & Wagner, N. J. (2014). Phys. Rev. E 89, 042301.  Google Scholar
First citationLutz-Bueno, V., Kohlbrecher, J. & Fischer, P. (2013). Rheol. Acta 52, 297–312.  CAS Google Scholar
First citationMay, R. P., Ibel, K. & Haas, J. (1982). J. Appl. Cryst. 15, 15–19.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationMiller, E. & Rothstein, J. P. (2007). J. Non-Newt. Fluid Mech. 143, 22–37.  CrossRef CAS Google Scholar
First citationMin Kim, J., Eberle, A. P., Kate Gurnon, A., Porcar, L. & Wagner, N. J. (2014). J. Rheol. 58, 1301–1328.  CrossRef CAS Google Scholar
First citationMütze, A., Heunemann, P. & Fischer, P. (2014). J. Rheol. 58, 1647–1672.  Google Scholar
First citationNemoto, F., Takabatake, F., Yamada, N. L., Takata, S. & Seto, H. (2024). J. Chem. Phys. 161, 164902.   Google Scholar
First citationOda, R., Panizza, P., Schmutz, M. & Lequeux, F. (1997). Langmuir 13, 6407–6412.  CrossRef CAS Google Scholar
First citationSalmon, J.-B., Colin, A., Manneville, S. & Molino, F. (2003). Phys. Rev. Lett. 90, 228303.  CrossRef PubMed Google Scholar
First citationSambasivam, A., Sangwai, A. V. & Sureshkumar, R. (2015). Phys. Rev. Lett. 114, 158302.  CrossRef PubMed Google Scholar
First citationSugahara, T., Akamatsu, M., Iwase, H., Takamatsu, Y., Sakai, K. & Sakai, H. (2020). Langmuir 36, 4695–4701.  CrossRef CAS PubMed Google Scholar
First citationTakata, S., Suzuki, J., Shinohara, T., Oku, T., Tominaga, T., Ohishi, K., Iwase, H., Nakatani, T., Inamura, Y., Ito, T., Suzuya, K., Aizawa, K., Arai, M., Otomo, T. & Sugiyama, M. (2015). JPS Conf. Proc. 8, 036020.  Google Scholar
First citationTakeda, M., Kusano, T., Matsunaga, T., Endo, H., Shibayama, M. & Shikata, T. (2011). Langmuir 27, 1731–1738.  CrossRef CAS PubMed Google Scholar
First citationTruong, M. & Walker, L. (2002). Langmuir 18, 2024–2031.  CrossRef CAS Google Scholar
First citationTruong, M. T. & Walker, L. M. (2000). Langmuir 16, 7991–7998.  CrossRef CAS Google Scholar
First citationUyama, M., Takahara, A., Higaki, Y., Yamada, N. L. & Iwase, H. (2022). J. Oleo Sci. 71, 1625–1637.  CrossRef CAS PubMed Google Scholar
First citationVan der Werff, J. C., Ackerson, B. J., May, R. P. & de Kruif, C. G. (1990). Physica A 165, 375–398.  CrossRef CAS Google Scholar
First citationVelichko, E., Tian, B., Nikolaeva, T., Koning, J., van Duynhoven, J. & Bouwman, W. G. (2019). Colloids Surf. A Physicochem. Eng. Asp. 566, 21–28.  CrossRef CAS Google Scholar
First citationWagner, N. J. & Mewis, J. (2021). Theory and applications of colloidal suspension rheology. Cambridge University Press.  Google Scholar
First citationZana, R. & Talmon, Y. (1993). Nature 362, 228–230.  CrossRef CAS Google Scholar
First citationZhang, F., Ilavsky, J., Long, G. G., Quintana, J. P., Allen, A. J. & Jemian, P. R. (2010). Metall. Mater. Trans. A 41, 1151–1158.  CrossRef Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds