The ubiquitous presence of humic substances (HS), combined with their ability to provide multiple sites for chemical reaction, makes them relevant to numerous biogeochemical processes such as mineral weathering, nutrient bioavailability, and contaminant transport. The reactivity of HS depends on their functional group chemistry and microstructure, which are in turn influenced by the composition of the surrounding media. In order to help towards an understanding of structure conformations and aggregation process of HS in soils and waters and to get a better knowledge of these kinds of materials, a fulvic acid (FA) has been modelled as a function of its ionic state under different conditions. Our proposed theoretical model based on the Temple-Northeastern-Birmingham (TNB) monomer fits well with experimental observations on the solubility (dipolar moment) and electronic and vibrational spectra of FAs. The presence of water molecules has a great stabilization effect on the electrostatic energy; this effect is greater as ionized rate increases. In vacuum, the non-ionized aggregated species are more stable than monomers because of the increase in their interaction due to H-bonding and non-bonding forces. When the molecules are ionized, no aggregation process takes place. In solution, the FA concentration is a critical factor for the aggregation. The system containing two FA molecules probably did not form aggregates because its equivalent concentration was too low. When the concentration was increased, the system gave rise to the formation of aggregates. The ionic state is another critical factor in the aggregation process. The ionized FA has a higher electric negative charge, which increases the energetic barriers and inhibits the approximation of FA caused by the Brownian movement.