The evolution of solution aggregates of the anionic form of the native monorhamnolipid (mRL) mixture produced by Pseudomonas aeruginosa ATCC 9027 is explored at pH 8.0 using both experimental and computational approaches. Experiments utilizing surface tension measurements, dynamic light scattering, and both steady-state and time-resolved fluorescence spectroscopy reveal solution aggregation properties. All-atom molecular dynamics simulations on self-assemblies of the most abundant monorhamnolipid molecule, L-rhamnosyl-beta-hydroxydecanoyl-beta-hydroxydecanoate (Rha-C10-C10), in its anionic state explore the formation of aggregates and the role of hydrogen bonding, substantiating the experimental results. At pH 8.0, at concentrations above the critical aggregation concentration of 201 mu M but below similar to 7.5 mM, small premicelles exist in solution; above similar to 75 mM, micelles with hydrodynamic radii of similar to 2.5 nm dominate, although two discrete populations of larger lamellar aggregates (hydrodynamic radii of similar to 10 and 90 nm) are also present in solution in much smaller number densities. The critical aggregation number for the micelles is determined to be similar to 26 monomers/micelle using fluorescence quenching measurements, with micelles gradually increasing in size with monorhamnolipid concentration. Molecular dynamics simulations on systems with between 10 and 100 molecules of Rha-C10-C10 indicate the presence of stable premicelles of seven monomers with the most prevalent micelle being similar to 25 monomers and relatively spherical. A range of slightly larger micelles of comparable stability can also exist that become increasing elliptical with increasing monomer number. Intermolecular hydrogen bonding is shown to play a significant role in stabilization of these aggregates. In total, the computational results are in excellent agreement with the experimental results.