We investigate the morphology and energetics of a self-associating model cationic surfactant in water using coarse-grained molecular dynamics simulations. We develop an algorithm to track micelle contours and quantify various microstructural features, such as contour length, persistence length, and mesh size. We demonstrate that branched and multiconnected structures govern the anomalous dependence of zero-shear viscosity on salt concentration. We predict reliably the end-cap energy of micelles, for the first time, directly from the simulations.In aqueous solutions, surfactant molecules spontaneously self-assemble into diverse geometrically complex and dynamically fluctuating morphologies. It has long been known that the emerging superstructures can range from spherical and elongated cylindrical to very long, flexible wormlike micelles with or without branches [1][2][3][4][5][6][7][8][9] and topologically rich knotted structures [7,8]. The diversity in microstructure and rheological properties make micellar solutions beneficial to numerous applications [10] as hydrofracking fluids in oil industry, turbulent friction drag reducing agents [8], thickening agents in consumer products, drug carriers in targeted delivery [11], and templates to create functional nanofluids with tunable mechanical or optical properties [12][13][14].Since the early work of Debye [1], the microstructural transitions in micellar solutions have been investigated both theoretically [4,[15][16][17][18] and experimentally [1,2,19,20]. It is now well-recognized that the molecular structure of co-surfactant or salt has a spectacular effect on the morphology. In particular, aromatic organic salts have stronger binding affinity to the micelles and induce enormous growth. Consequently, micelles become very long, flexible and entangle even at relatively low surfactant concentration c D . Dilute solutions with spherical or short cylindrical structures exhibit Newtonian fluid rheology [19][20][21]. In contrast, solutions above the overlap concentration φ * consist of very long-thread like structures with contour lengths that span from a few 10s of nm to several µm, and show viscoelastic behavior reminiscent of flexible polymer solutions [4,19,20]. However, unlike polymers, WLMs can merge and undergo reversible breaking at time scales that are detectable both in scattering experiments and simulations. Under non-equilibrium conditions, such as, under shear flow, the structure, dynamics, and the resultant rheological properties could change dramatically: two notable examples are the shear induced structure (SIS) formation, and shear banding [19][20][21][22][23][24][25][26][27]. Due to such dynamical complexities, a quantitative description of the microstructure of micellar fluids is incomplete. In this paper, we present a comprehensive simulation study of self-assembly, emerging structures, length scales, and the energetics of a model surfactant solution with explicit solvent, electrostatic and hydrodynamic interactions.Despite decades of research aimed ...