
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38594372
58854
10.1038/s41598-024-58854-6
Article
Designing the structure of cationic star-shaped trimeric surfactants most active in micelle formation using molecular connectivity indices
Mozrzymas Anna anna.mozrzymas@upwr.edu.pl

https://ror.org/05cs8k179 grid.411200.6 0000 0001 0694 6014 Department of Physics and Biophysics, Wrocław University of Environmental and Life Sciences, ul. Norwida 25, 50-375 Wrocław, Poland
9 4 2024
9 4 2024
2024
14 832316 1 2024
3 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
A model expressing the relationship between the logarithm of critical micelle concentration of cationic star-shaped trimeric surfactants and topological indices was obtained using only molecular connectivity indices. Based on the obtained model, the exemplary compound most active in micelle formation was designed. The analysis of the influence of various structural factors on the value of the critical micelle concentration was supported by atomic charge studies. The obtained model will be used to design new star-shaped trimeric surfactants that are more active in formation of micelle as well as to predict their critical micelle concentration.

Subject terms

Polymer chemistry
Chemistry
Computational chemistry
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Oligomeric surfactants with multiple hydrophobic and hydrophilic groups are a new class of surfactants which become attractive in recent years. Dimeric surfactant, also called gemini surfactants, are simplest oligomeric surfactants which contain two hydrophobic and two hydrophilic groups connected by the spacer group. Trimeric surfactants are natural extension of the gemini surfactants. Their structure represent intermediate structure between dimeric surfactants and higher oligomeric ones. The first report on trimeric surfactants was by Raoul Zana et al.1–3. Generally, the oligomeric surfactants are constructed by two or more hydrophobic tails and polar head groups linked by the spacer groups. The spacer groups of trimeric and higher oligomerization degree surfactants can be linear, ring-like or star-shaped therefore the oligomeric surfactants are categorized into linear, ring-like and star-shaped, and their aggregation behavior strongly depends on their topological structures4. The star-shaped trimeric surfactants exhibit more unique self-aggregation behavior in aqueous solution compared to their linear analogues4.

Generally, the trimeric surfactants show an excellent surface-active and self-aggregation properties in aqueous solutions. The experimental studies5–9 and also the theoretical considerations10 show that the critical micelle concentration decreases with the increasing the degree of oligomerization. Thus, the critical micelle concentration of trimeric surfactants is much lower than those of the corresponding monomeric and dimeric analogues. Also, the cationic trimeric compounds exhibit strong antimicrobial activity. They are active against board range of microorganism such as bacteria and fungi11,12 and probably can be used against some viruses.

The main aim of this work is designing the structure of cationic star-shaped trimeric surfactants most active in micelle formation otherwise to study the effect of the structure modifications on the critical micelle concentration value. This work is a continuation of studies on structure–property relationship (QSPR) of cationic surfactants13–15. The subject of the previous studies was the critical micelle concentration (cmc) of cationic monomeric13 and cationic gemini14,15 surfactants. In the papers13–15, the relationship of the critical micelle concentration and the structure of cationic surfactants was investigated using the molecular connectivity indices only. In the present work these topological indices were also used to study the influence of the chains structure on the value of critical micelle concentration of star-shaped cationic trimeric surfactants. As was suggested in paper4 the topological structure of the oligomeric surfactants strongly affect the aggregation behavior of these compounds. Therefore, the topological descriptors like the molecular connectivity indices16 can be very good representation of the oligomeric surfactant structure in studies of self-aggregation properties, and in particular of the critical micelle concentration, as was shown in the previous papers concerning the gemini surfactants14,15.

The critical micelle concentration depends not only on geometrical structure but also on a number of other parameters, among them a kind of counterion. Therefore, in order to minimize the influence of factors other than geometrical on cmc value, only the star-shaped cationic surfactants with bromides as counterions were taken into account. Although, the obtained QSPR model has been derived for compounds with fixed counterions, it can be used to design the structure of cationic trimeric star-shaped surfactants with any kind of counterion because the impact of modifying the geometrical structure on changes in the cmc value, its decrease or increase, should be the same.

The semi-empirical calculations of the atomic charges were also performed to investigate the effect of branches and heteroatoms contained in the spacer group on the critical micelle concentration of studied trimeric surfactants.

Results and discussion

All investigated cationic trimeric surfactants are star-shaped type surfactants. The geometrical structures of investigated compounds are significantly different. The entire data set includes thirteen training set compounds (compounds 1–13) and five test compounds (compounds T1–T5). The structures of the molecules taken from literature along with the logarithms of the literature cmc values are shown in Methods section.

Molecular connectivity indices were calculated basing on the graphic structural formula of the molecules 1–13 using the expressions contained in Methods section (Eqs. 3–5). The values of five molecular connectivity indices from zero to fourth order and five valence molecular connectivity indices from zero to fourth order of training compounds are given in Table 1.Table 1 Values of molecular connectivity indices.

No	0χ	1χ	2χ	3χc	4χpc	0χν	1χν	2χν	3χcν	4χpcν	
1	36.53319	23.46739	19.70997	3.82544	4.05433	36.2447	22.65065	18.47576	3.39712	3.57442	
2	40.77583	26.46739	21.83129	3.82544	4.05433	40.48734	25.65065	20.59708	3.39712	3.57442	
3	45.01847	29.46739	23.95261	3.82544	4.05433	44.72998	28.65065	22.71840	3.39712	3.57442	
4	43.38657	27.64893	23.75389	4.69147	4.76746	40.19419	24.30208	19.38003	3.57390	3.30413	
5	51.87185	33.64893	27.99653	4.69147	4.76746	48.67947	30.30208	23.62267	3.57390	3.30413	
6	56.11449	36.64893	30.11785	4.69147	4.76746	52.92211	33.30208	25.74399	3.57390	3.30413	
7	44.54127	28.59842	25.07538	5.09972	5.98618	40.23638	23.89611	19.00814	3.65964	3.82895	
8	53.02655	34.59842	29.31802	5.09972	5.98618	48.72167	29.89611	23.25078	3.65964	3.82895	
9	44.88657	28.89157	27.21624	6.81279	7.88938	41.82432	26.00647	22.85695	5.06880	5.50002	
10	49.12921	31.89157	29.33757	6.81279	7.88938	46.06697	29.00647	24.97828	5.06880	5.50002	
11	32.49730	22.69453	17.90271	1.42887	2.87217	30.41853	19.56576	13.81274	0.79057	1.43061	
12	33.65200	23.64402	19.31693	1.83712	3.42036	31.08205	19.90978	14.43237	0.98601	1.54138	
13	50.17281	34.23886	27.51366	2.80888	5.40758	44.53661	28.39050	20.17867	1.22422	1.86717	

Based on the values of connectivity indices (Table 1) and the logarithms of literature cmc values of training compounds 1–13, using the polynomial regression analysis and stepwise method, the two-variable equation (Eq. 1), expressing the relation between logarithm of cmc and the molecular connectivity indices, have been obtained:1 Log10cmc=-(1.10878∓0.18486)+(0.00428±0.00058)·(2χ)2-(0.00709∓0.00055)·(1χν)2

The statistical characteristics of the equation variables are given in Table 2.Table 2 Statistical characteristics of variables included in Eq. (1).

Variable	Coefficient	Standard error	t-value	p-value	
Constant	− 1.10878	0.18486	− 5.9979	0.00013	
(2χ)2	0.00428	0.00058	7.3276	0.00003	
(1χν)2	− 0.00709	0.00055	− 12.9841	0.00000	

The obtained model (Eq. 1) contains first-order valence molecular connectivity index 1χν and second-order molecular connectivity index 2χ. The first-order valence molecular connectivity index 1χν is path-type index and it represents one-bonds fragments in a molecule. The values of this index depend on the isomers of the compound16 and decrease with increasing in branching. The molecular valence connectivity index 1χν is a valence connectivity index thus differentiates heteroatoms and multiple bonds. The second-order 2χ molecular connectivity index appearing in the model (Eq. 1) represents two-bonds fragments in the molecule. The values of this index also depend on the isomers of the compound16, and in this case the values increase with increasing in branching. Also, the values of 2χ and 1χν indices increase with the increase in the number of atoms in the molecule by extending the hydrocarbon chains or adding atoms to the chains through branching.

From Table 1 it can be concluded that the values of 2χ index are smaller compared to the values of 1χν index. The obtained model (Eq. 1) contains the (1χν)2 variable with a negative coefficient and the (2χ)2 variable with a positive one, and the absolute value of these coefficients are 0.00709 and 0.00428, respectively. Therefore, the analysis of equation (Eq. 1) allows us to conclude that as the 1χν index increases, the cmc decreases.

The graphical comparison of the calculated Log10cmc values using Eq. (1) and the experimental Log10cmc values is shown in Fig. 1.Figure 1 Scatter plot of the calculated Log10cmc versus the experimental Log10cmc for compounds of training set (r = 0.981, F = 280.480, s = 0.173, r2adj = 0.959, QLOO2 = 0.925) (rhombus) and test set (Rpred2 = 0.531) (triangle).

The Log10cmc values calculated using Eq. (1) and the experimental Log10cmc values for studied trimeric surfactants, along with the values of residuals, are contained in Table 3.Table 3 Calculated and literature values Log10cmc of compounds from training and test set.

Compound	Experimental Log10cmc	Calculated Log10cmc	Residual	
1	− 2.9956817	− 3.08207	0.08639	
2	− 3.6840317	− 3.73189	0.04786	
3	− 4.5114517	− 4.47077	− 0.04068	
4	− 2.7619518	− 2.87906	0.11711	
5	− 4.2941418	− 4.26129	− 0.03285	
6	− 5.2076118	− 5.08599	− 0.12162	
7	− 2.3555619	− 2.46403	0.10847	
8	− 3.4894519	− 3.76371	0.27426	
9	− 3.0000020	− 2.73117	− 0.26883	
10	− 3.3010320	− 3.38734	0.08631	
11	− 2.2898821	− 2.44998	0.16010	
12	− 2.6382721	− 2.32084	− 0.31743	
13	− 3.6798522	− 3.58075	− 0.09910	
T1	− 2.7958823	− 2.79056	− 0.00532	
T2	− 3.7520323	− 3.39585	− 0.35618	
T3	− 3.2924324	− 3.45880	0.16637	
T4	− 3.9065818	− 3.52565	0.38093	
T5	− 2.6989720	− 2.16405	0.53492	

The plot of residuals versus the experimental values of Log10cmc is shown in Fig. 2.Figure 2 Plot of residuals versus the experimental Log10cmc values for training set (rhombus) and test compounds (triangle).

As can be seen from Table 3 and Figs. 1–2, the calculated values of Log10cmc using Eq. (1) are very close to the from literature ones.

As shown in Table 6 given in Methods section, the Eq. (1) has been obtained based on molecules with a different topological structure. These compounds differ in the length of the hydrophobic chains as well as the structure of spacer groups. Therefore, it is worth analyzing, based on the obtained equation, the influence of hydrophobic tails and spacer group on the values of logarithm of cmc.

Hydrophobic tails effect

The experimental cmc values and thus Log10cmc of star-shaped cationic trimeric surfactants decrease when the tails are lengthened between ten and fourteen (compounds 1–3) or ten and sixteen (compounds 4–6 and T4) carbon atoms. The corresponding plots of the calculated, using Eq. (1), and the experimental Log10cmc versus the alkyl chain carbon number of compounds 1–3 and compounds 4–6 and T4 are shown in Figs. 3–4.Figure 3 Calculated (square) and experimental17 (triangle) Log10cmc values of star-shaped trimeric compounds 1–3 versus the number of alkyl chain carbon atoms (n).

Figure 4 Calculated (square) and experimental18 (triangle) Log10cmc values of star-shaped trimeric compounds 5–6 and T4 versus the number of alkyl chain carbon atoms (n).

As shown in Figs. 3–4, the Log10cmc values calculated using Eq. (1) also decrease with tails lengthening from ten to eighteen (Figs. 3–4) carbon atoms, and as shown in Figs. 3–4, with this range the dependence of both the calculated and also experimental logarithm values of cmc on the number of carbon atoms in the alkyl chains is linear. As can be seen in Fig. 3–4, the calculated Log10cmc values are very close to experimental values. For experimental Log10cmc values of compounds in Fig. 4 there is a little deviation from linearity for the twelfth carbon atoms chain, but this deviation is probably within the margin of experimental error. Thus, it can be concluded that in the range between the ten and the eighteen carbon atoms the dependence of Log10cmc on the number of alkyl chain carbon atoms (n) of star-shaped cationic trimeric surfactants can be described by linear function.

However, in the case of longer chains, above eighteen carbon atoms, the dependence of the Log10cmc values calculated using Eq. (1) on the number of carbon atoms will be non-linear due to the non-linear nature of this equation. The non-linear dependence of logarithm of cmc on the chains length has been observed for various cationic gemini surfactants of different nature and flexibility spacer groups or tails25–27. Also, the non-linear relationship between Log10cmc and the number of tail carbon atoms for cationic gemini surfactants was obtained basing on theoretical model28.

Spacer group effect

The structure and nature of the spacer group play an important role in micelle formation of oligomeric surfactants. The studies on gemini surfactants29–33 show that such features of the spacer group as flexibility and hydrophobicity have significant impact on the aggregation behavior in aqueous solution of these compounds. Therefore, the effect of such factors which influence the flexibility and hydrophobicity as spacer’s length, branching and heteroatoms on the cmc of star-shaped trimeric surfactants have been investigated using Eq. (1) (Table 4). In Table 4 the star-shaped type trimeric surfactants with different spacers groups structures and corresponding Log10cmc values calculated using Eq. (1) are shown.Table 4 Calculated Log10cmc values of star-shaped trimeric surfactants with different spacer group (S1-S10).

	
Compound	X	Y	Calculated Log10cmc	
S1		–CH2–	− 4.80429	
S2		–CH2–	− 4.43451	
S3	CH	–CH2–	− 4.09246	
S4	N	–CH2–	− 4.01400	
S5	N	–NH–	− 3.75740	
S6	N	–O–	− 3.61353	
S7	N	–CH2–CH2–	− 4.47315	
	
S8	N	–CH3	− 3.78889	
S9	N	–NH2	− 3.51421	
S10	N	–OH	− 3.43116	

The inspection of data contained in Table 4 shows that the greatest value of critical micelle concentration is for compound S6 and the lowest for compound S1.

In the case of compounds with different central groups in the spacer group (compounds S1–S4), the lowest cmc value is for S1 compound with the cyclohexane-based spacer group, and the highest for compound S4 with the central nitrogen atom in the spacer group. This lowest cmc value of compound S1 may be due to the fact that the spacer is more hydrophobic which may cause its incorporation into the interior of the formed micelle. Also, as shown in Table 4, the cmc value of compound S2 is lower than cmc value of compound S4. Similar results was obtained by Nacham et al.21. The authors of paper21 have studied among others the ionic liquid-based star-shaped trimeric surfactants containing triethylamine and triethylbenzene spacer groups. As was shown in21 the cmc value of IL-based trimeric surfactant containing triethylbenzene spacer groups is lower than the cmc of that containing triethylamine spacer group with identical alkyl tails and head groups. The authors of this paper also suggest that this difference of cmc values may be due to the higher hydrophobicity imparted by the benzyl core in investigated IL-trimeric surfactant. Table 4 also shows that cmc value of compound S1 is lower compared to the cmc value of compound S2. This difference may result from the fact that substituted cyclohexane ring is more flexible compared to the benzene one, and thus the cyclohexane-based spacer group easily incorporates into the interior of the formed micelle.

Among the compounds having nitrogen as the central group and the same number of non-hydrogen atoms in the spacer group (compounds S4–S6 and S8–S10), the lowest critical micelle concentration is for compound S4 and the greatest for compound S10. In the case of the compounds with straight spacer chains (compounds S4–S6), the cmc value increases in the order: –CH2– < –NH– < –O–. The similar order is for compounds with branched spacer group (compounds S8–S10), the order of increasing cmc value is: –CH3 < –NH2 < –OH. In addition, the comparison of the critical micelle concentration values of compounds with straight chains of star-shaped spacers (compounds S4–S6) with the corresponding compounds having branched spacer group (compounds S8–S10) shows that the critical micelle concentration values of compounds having the branched spacer groups are greater compared to the cmc values of the corresponding compounds with straight spacer’s chains.

In summary, Table 4 shows that in the case of star-shaped spacers having the same number of non-hydrogen atoms, the branches and the heteroatoms cause the increase the critical micelle concentration value.

The effect of spacer’s branches and heteroatoms on cmc value of compounds presented in Table 4 (compounds S4–S6 and S8–S10) was also analyzed using the atomic charges. The total charges of different spacer’s functional groups found in the Y position of compounds S4–S6 and S8–S10 are presented in Table 5.Table 5 Total charges of the different functional groups found in the Y position (see Table 4).

Straight spacer (S4–S6)	Branched spacer (S8–S10)	
Group	Charge	Group	Charge	
Y1	Y2	Y3	Y1	Y2	Y3	
–CH2–	0.0057	0.0179	0.0435	–CH3	0.0599	0.0514	0.0550	
–NH–	− 0.211	− 0.1377	− 0.1263	–NH2	− 0.0065	− 0.0293	− 0.030	
–O–	− 0.3271	− 0.2867	− 0.2987	–OH	− 0.099	− 0.0961	− 0.1037	

The AM1 semi-empirical calculations (Table 5) show that the atoms of high electronegativity such oxygen and nitrogen introduce a great negative charge to the spacer group, greater than the carbon atom thus the total charges of functional groups containing heteroatoms are negative. Moreover, the total charges of the various functional groups shown in Table 5 show that for compounds both straight and branched spacer group the total charge changes from positive for functional group with carbon atom to negative for functional groups with nitrogen or oxygen. Additionally, the negative charge is greater in the case of a group with an oxygen atom than with a nitrogen atom. Comparing these results with the data from Table 4, it can be concluded that the cmc values increase with the increase in the negative charge of the functional group.

The data in Table 5 also show that the total (negative) charge of the functional group, in the branched spacer chain compared to that of the corresponding functional group in the straight chain, is lower for branched spacer when the non-hydrogen atom in the functional group in Y position is oxygen or nitrogen, and greater (positive) charge also for branched spacer group when the non-hydrogen atom in the functional group in Y position is carbon atom. Comparing these results with the data from Table 4, it can be seen that when the non-hydrogen atom in the functional group is carbon atom, the results from Table 5 agree very well with those in Table 4 discussed above, i.e. the greater the positive charge, the greater the critical micelle concentration. In the case of functional groups containing heteroatoms (Table 4), it can be seen that the cmc value is greater for compounds with branched spacer group for which the negative charge of the functional group is lower than the negative charge of analogous functional group in straight spacer chain (Table 5).

The obtained model also shows that the cmc value of star-shaped trimeric surfactants decreases with the increase in the chains length of the spacer group (compounds S4 and S7). This is due to the fact that changes in the length of hydrocarbon tails have a significant impact on cmc values, and thus on obtained model. However, the experimental data34 shows the opposite conclusion. As shown in34, the cmc value of 3C12tris-s-Q cationic trimeric surfactants is greater for s = 6 than for s = 3. This suggests that the cmc increases with increasing spacer chains length (s) from three to about six methylene groups. A similar effect of spacer chains length is observed for dimeric surfactants29,30. The cmc values of 12–s–12 gemini surfactants increase with increasing number of spacer carbon atoms, up to four or five methylene groups, and then decrease with further elongation of the spacer chain29,30. Therefore, probably the increase in cmc value of trimeric surfactants with spacer chains length may have maximum at about six carbon atoms and then, as for gemini surfactants, the cmc may decrease with further chains length elongation. To the best of the author's knowledge, spacer chains lengthening above six methylene groups for star-shaped trimeric surfactants has not been studied. Therefore, the obtained model should be used probably for compounds with more than six carbon atoms in each chain of the trimeric surfactant spacer group when examining the effect of the spacer chains elongation on the cmc value.

Finally, if we examine only the changes in cmc value due to changes in the geometrical structure, the obtained equation (Eq. 1), although derived for bromide compounds, can also be used for compounds with another counterion, because the changes in cmc, its increase or decrease, should be the same. In other words, when examining changes in the cmc of bromide star-shaped trimeric surfactants, we can expect that in the case of chloride analogous compounds will be similar. For example, for chloride compounds 1,1,1-Tris[2-hydroxy-3-(dodecyldimethylammonio)-propoxymethyl]ethane Trichloride35 and the so-called III-12-4 trimeric surfactant36, for which the experimental cmc values are 0.223 (mM) and 0.15 (mM), respectively, the calculated values using Eq. (1) are 0.198 (mM) and 0.187 (mM), respectively. So, the calculated cmc value decreases as the experiment shows.

Methods

Experimental data

All investigated molecules are trimeric quaternary ammonium bromide surfactants. The structures of all considered compounds along with the logarithms of the literature cmc values17–24 are shown in Table 6.Table 6 Molecular structures and logarithms of experimental cmc values of training (1–13) and test (T1–T5) compounds.

No	Molecular structure	n	Log10cmc	Ref	
1

2

3

		10

12

14

	− 2.99568

− 3.68403

− 4.51145

	17	
4

5

6

		10

14

16

	− 2.76195

− 4.29414

− 5.20761

	18	
7

8

		8

12

	− 2.35556

− 3.48945

	19	
9

10

		10

12

	− 3.00000

− 3.30103

	20	
11		8	− 2.28988	21	
12		8	− 2.63827	21	
13		12	− 3.67985	22	
T1

T2

		10

12

	− 2.79588

− 3.75203

	23	
T3		12	− 3.29243	24	
T4		12	− 3.90658	18	
T5		8	− 2.69897	20	

The literature values of critical micelle concentration were given in (mM)17–24. The cmc values expressed in molar units (M) have been converted to logarithms of cmc (Table 6).

Most values of cmc were measured at 25 °C in aqueous solution. The cmc of compound T3 was measured at 22 oC.

Molecular connectivity indices

Molecular connectivity indices, some of the topological descriptors to characterize molecules in structure–property and structure–activity studies, are calculated from the molecular graph. The molecular graph is a graphical representation of the structural formula of the chemical compound, in which vertices represents atoms and edges symbolize covalence bonds.

The first connectivity index was proposed by Randić37 and was defined as:2 χ=∑kδiδjk-0.5

where δi is a connectivity degree i.e. the number of non-hydrogen atoms to which the i-th non-hydrogen atom is bonded. The Kier and Hall molecular connectivity indices16 are generalizations of Randić’s connectivity index and the m-th order Kier-Hall molecular connectivity index is defined as16:3 mχk=∑j=1nm∏i=1m+1δij-0.5

where δi is a connectivity degree, m is the order of the connectivity index, k denotes type of the fragment of the molecule for example: path (p), cluster (c) and path-cluster (pc), nm is the number of fragments of type k and order m.

The m-th order valence molecular connectivity index is defined16:4 mχkν=∑j=1nm∏i=1m+1δiνj-0.5

where the valence connectivity degree δiν is defined as:5 δiν=Ziν-hiZi-Ziν-1

where Ziν is the number of valence electrons in the i-th atom, hi is the number of hydrogen atoms connected to the i-th atom and Zi is the number of all electrons in the i-th atom.

The values of the molecular and valence molecular connectivity indices of training compounds are contained in Table 1.

Atomic charges

Atomic charges were calculated using the semi-empirical molecular orbital package MOPAC 7 included in the VEGA program38,39, employing the semi-empirical AM1 method. All calculated atomic charges are expressed in atomic units (a.u.).

Statistics

Each formula expressing the relationship between the Log10cmc and molecular connectivity indices was generated using the least-squares method and the final equation was obtained using the stepwise method. Pearson correlation coefficient (r), the adjusted coefficient of determination (r2adj), the standard deviation of the fit (s) and the Fisher ratio (F) were used to select the best model. The model obtained was selected according to following principles: highest correlation coefficient, adjusted coefficient of determination and Fisher value, the lowest standard deviation of the fit and also the smallest possible number of significant descriptors in the model. For good quality QSPR model the positive value of correlation coefficient (r) should be closer to 1. High values of F-test indicate that the model is statistically significant. The number of variables in the model should not exceed the number of compounds divided by five40. The statistical characteristics of the equation variables includes standard error, t-value and p-value. High absolute Student t value of the variable expresses that the coefficient of the variable is significantly larger than the standard error. Variable with p value below 0.05 is considered statistically significant41.

The leave-one-out cross-validated correlation coefficient (QLOO2) and predictive correlation coefficient (Rpred2) were also used to indicate the internal and external validation of derived model. The model is considered to be excellent, if QLOO2 is equal or more than 0.942. For acceptable QSPR model, the values of QLOO2 and Rpred2 should be more than 0.540,43.

The statistical calculations were made using the program STATISTICA 1244.

Conclusions

The aim of the presented work was to find a simple equation expressing the dependence of the critical micelle concentration of cationic star-shaped trimeric surfactants on their geometrical structures represented by topological indices, that will allow for examining the impact of structure modifications on the cmc values and thus it can be helpful in the design of novel star-shaped cationic trimeric surfactants. Employing the polynomial regression analysis, the equation with two molecular connectivity indices was obtained. Using the obtained model, the influence of hydrophobic tails length and structure of star-shaped spacer group on the critical micelle concentration was examined.

The analysis of the influence of tails length in the range between ten and eighteen carbon atoms confirmed the experimental results and the Klevens equation45 which expresses a linear dependence of logarithm of cmc on the number of carbon atoms in the alkyl chains.

The analysis of star-shaped spacers with the same number of non-hydrogen atoms allows us to conclude that heteroatoms cause the increase the cmc. These results were confirmed by atomic charge analysis. Functional groups containing highly electronegativity atoms, such as oxygen or nitrogen, introduce into the spacer group the negative charges, which cause an increase in the critical micelle concentration. The spacer group analysis also reveals that the branches cause the increase the critical micelle concentration too.

Acknowledgements

This paper was financially supported by the Department of Physics and Biophysics at Wrocław University of Environmental and Life Sciences. The article processing charge (APC) is financed/co-financed by Wrocław University of Environmental and Life Sciences. The statistical calculations were performed using the program Statistica 12 provided by Wrocław University of Environmental and Life Sciences.

Author contributions

A.M. conceived the topic of the study, made the calculations, performed the analysis and wrote the paper.

Funding

This paper is financially supported by the Department of Physics and Biophysics at Wrocław University of Environmental and Life Sciences. The article processing charge is financed/co-financed by Wrocław University of Environmental and Life Sciences.

Data availability

All data generated or analysed during study are included in this published article. The experimental data were taken from literature.

Competing interests

The author declares no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Zana R Levy H Papoutsi D Beinert G Micellization of two triquaternary ammonium surfactants in aqueous solution Langmuir 1995 11 3694 3698 10.1021/la00010a018
Zana, R., Levy, H., Papoutsi, D. & Beinert, G. Micellization of two triquaternary ammonium surfactants in aqueous solution. Langmuir 11, 3694–3698. 10.1021/la00010a018 (1995).10.1021/la00010a018
2. Danino D Talmon Y Levy H Beinert G Zana R Branched threadlike micelles in an aqueous solution of a trimeric surfactants Science 1995 269 1420 1421 10.1126/science.269.5229.1420 17731153
Danino, D., Talmon, Y., Levy, H., Beinert, G. & Zana, R. Branched threadlike micelles in an aqueous solution of a trimeric surfactants. Science 269, 1420–1421. 10.1126/science.269.5229.1420 (1995).17731153 10.1126/science.269.5229.1420
3. Zana R Dimeric and oligomeric surfactants. Behavior at interfaces and in aqueous solution: A review Adv. Colloid. Interf. Sci. 2002 97 205 253 10.1016/S0001-8686(01)00069-0
Zana, R. Dimeric and oligomeric surfactants. Behavior at interfaces and in aqueous solution: A review. Adv. Colloid. Interf. Sci. 97, 205–253. 10.1016/S0001-8686(01)00069-0 (2002).10.1016/S0001-8686(01)00069-0
4. Fan Y Wang Y Self-assembly and functions of star-shaped oligomeric surfactants Langmuir 2018 34 11220 11241 10.1021/acs.langmuir.8b00290 29616549
Fan, Y. & Wang, Y. Self-assembly and functions of star-shaped oligomeric surfactants. Langmuir 34, 11220–11241. 10.1021/acs.langmuir.8b00290 (2018).29616549 10.1021/acs.langmuir.8b00290
5. Laschewsky A Wattebled L Arotcaréna M Habib-Jiwan J-L Rakotoaly RH Synthesis and properties of cationic oligomeric surfactants Langmuir 2005 21 7170 7179 10.1021/la050952o 16042438
Laschewsky, A., Wattebled, L., Arotcaréna, M., Habib-Jiwan, J.-L. & Rakotoaly, R. H. Synthesis and properties of cationic oligomeric surfactants. Langmuir 21, 7170–7179. 10.1021/la050952o (2005).16042438 10.1021/la050952o
6. Yoshimura T Yoshida H Ohno A Esumi K Physicochemical properties of quaternary ammonium bromide-type trimeric surfactants J. Colloid. Interf. Sci. 2003 267 167 172 10.1016/S0021-9797(03)00694-5
Yoshimura, T., Yoshida, H., Ohno, A. & Esumi, K. Physicochemical properties of quaternary ammonium bromide-type trimeric surfactants. J. Colloid. Interf. Sci. 267, 167–172. 10.1016/S0021-9797(03)00694-5 (2003).10.1016/S0021-9797(03)00694-5
7. In M Bec V Aguerre-Chariol O Zana R Quaternary ammonium bromide surfactant oligomers in aqueous solution: Self-association and microstructure Langmuir 2000 16 141 148 10.1021/la990645g
In, M., Bec, V., Aguerre-Chariol, O. & Zana, R. Quaternary ammonium bromide surfactant oligomers in aqueous solution: Self-association and microstructure. Langmuir 16, 141–148. 10.1021/la990645g (2000).10.1021/la990645g
8. Chen J Qiao M Gao N Ran Q Wu J Shan G Qi S Wu S Cationic oligomeric surfactants as novel air entraining agents for concrete Colloid. Surf. A Physicochem. Eng. Asp. 2018 538 686 693 10.1016/j.colsurfa.2017.11.065
Chen, J. et al. Cationic oligomeric surfactants as novel air entraining agents for concrete. Colloid. Surf. A Physicochem. Eng. Asp. 538, 686–693. 10.1016/j.colsurfa.2017.11.065 (2018).10.1016/j.colsurfa.2017.11.065
9. Jurasin D Habus I Filipowić-Vinceković N Role of the alkyl chain number and head groups location on surfactants self-assembly in aqueous solution Colloid. Surf. A Physicochem. Eng. Asp. 2010 368 119 128 10.1016/j.colsurfa.2010.07.02
Jurasin, D., Habus, I. & Filipowić-Vinceković, N. Role of the alkyl chain number and head groups location on surfactants self-assembly in aqueous solution. Colloid. Surf. A Physicochem. Eng. Asp. 368, 119–128. 10.1016/j.colsurfa.2010.07.02 (2010).10.1016/j.colsurfa.2010.07.02
10. Maiti PK Lansac Y Glaser MA Clark NA Rouault Y Self-assembly in surfactant oligomers: A coarse-grained description through molecular dynamic simulation Langmuir 2002 18 1908 1918 10.1021/la0111203
Maiti, P. K., Lansac, Y., Glaser, M. A., Clark, N. A. & Rouault, Y. Self-assembly in surfactant oligomers: A coarse-grained description through molecular dynamic simulation. Langmuir 18, 1908–1918. 10.1021/la0111203 (2002).10.1021/la0111203
11. Murguia MC Cristaldi MD Porto A Conza JD Grau RJ Synthesis, surface-active properties, and antimicrobial activities of new neutral and cationic trimeric surfactants J. Surfact. Deterg. 2008 11 41 48 10.1007/s11743-007-1052-4
Murguia, M. C., Cristaldi, M. D., Porto, A., Conza, J. D. & Grau, R. J. Synthesis, surface-active properties, and antimicrobial activities of new neutral and cationic trimeric surfactants. J. Surfact. Deterg. 11, 41–48. 10.1007/s11743-007-1052-4 (2008).10.1007/s11743-007-1052-4
12. Zhou C Wang F Chen H Li M Qiao F Liu Z Hou Y Wu C Fan Y Liu L Wang S Wang Y Selective antimicrobial activities and action mechanism of micelles self-assembled by cationic oligomeric surfactants ACS Appl. Mater. Interfaces 2016 8 4242 4249 10.1021/acsami.5b12688 26820390
Zhou, C. et al. Selective antimicrobial activities and action mechanism of micelles self-assembled by cationic oligomeric surfactants. ACS Appl. Mater. Interfaces 8, 4242–4249. 10.1021/acsami.5b12688 (2016).26820390 10.1021/acsami.5b12688
13. Mozrzymas A Różycka-Roszak B Prediction of critical micelle concentration of cationic surfactants using connectivity indices J. Math. Chem. 2011 49 276 289 10.1007/s10910-010-9738-7
Mozrzymas, A. & Różycka-Roszak, B. Prediction of critical micelle concentration of cationic surfactants using connectivity indices. J. Math. Chem. 49, 276–289. 10.1007/s10910-010-9738-7 (2011).10.1007/s10910-010-9738-7
14. Mozrzymas A Molecular connectivity indices for modeling of the critical micelle concentration of cationic (chloride) gemini surfactants Colloid Polym. Sci. 2017 295 75 87 10.1007/s00396-016-3979-3 28111493
Mozrzymas, A. Molecular connectivity indices for modeling of the critical micelle concentration of cationic (chloride) gemini surfactants. Colloid Polym. Sci. 295, 75–87. 10.1007/s00396-016-3979-3 (2017).28111493 10.1007/s00396-016-3979-3
15. Mozrzymas A Modelling of the critical micelle concentration of cationic gemini surfactants using molecular connectivity indices J. Solution Chem. 2013 42 2187 2199 10.1007/s10953-013-0095-6 24319301
Mozrzymas, A. Modelling of the critical micelle concentration of cationic gemini surfactants using molecular connectivity indices. J. Solution Chem. 42, 2187–2199. 10.1007/s10953-013-0095-6 (2013).24319301 10.1007/s10953-013-0095-6
16. Kier LB Hall LH Molecular connectivity in structure- activity analysis 1986 Research Studies Press Ltd
Kier, L. B. & Hall, L. H. Molecular connectivity in structure- activity analysis (Research Studies Press Ltd, 1986).
17. Morita T Yada S Yoshimura T Linear- and star-type quaternary ammonium salt-based trimeric surfactants: Effect of structure on adsorption and aggregation properties Colloid. Surf. A Physicochem. Eng. Asp. 2023 656 130364 1 9 10.1016/j.colsurfa.2022.130364
Morita, T., Yada, S. & Yoshimura, T. Linear- and star-type quaternary ammonium salt-based trimeric surfactants: Effect of structure on adsorption and aggregation properties. Colloid. Surf. A Physicochem. Eng. Asp. 656(130364), 1–9. 10.1016/j.colsurfa.2022.130364 (2023).10.1016/j.colsurfa.2022.130364
18. Liu XG Xing XJ Gao ZN Synthesis and physicochemical properties of star-like cationic trimeric surfactants Colloid. Surf. A: Physicochem. Eng. Asp. 2014 457 374 381 10.1016/j.colsurfa.2014.06.019
Liu, X. G., Xing, X. J. & Gao, Z. N. Synthesis and physicochemical properties of star-like cationic trimeric surfactants. Colloid. Surf. A: Physicochem. Eng. Asp. 457, 374–381. 10.1016/j.colsurfa.2014.06.019 (2014).10.1016/j.colsurfa.2014.06.019
19. El-Said WA Moharram AS Hussein EM El-Khawaga AM Synthesis, characterization, and applications of some new trimeric-type cationic surfactants J. Surfact. Deterg. 2018 21 343 353 10.1002/jsde.12041
El-Said, W. A., Moharram, A. S., Hussein, E. M. & El-Khawaga, A. M. Synthesis, characterization, and applications of some new trimeric-type cationic surfactants. J. Surfact. Deterg. 21, 343–353. 10.1002/jsde.12041 (2018).10.1002/jsde.12041
20. Zhu H Guo J Yang C Liu S Cui Y Zhong X Synthesis of adamantine-based trimeric cationic surfactants Synthetic Commun. 2013 43 1161 1167 10.1080/00397911.2011.625134
Zhu, H. et al. Synthesis of adamantine-based trimeric cationic surfactants. Synthetic Commun. 43, 1161–1167. 10.1080/00397911.2011.625134 (2013).10.1080/00397911.2011.625134
21. Nacham O Martin-Perez A Steyer DJ Trujillo-Rodriguez MJ Anderson JL Pino V Afonso AM Interfacial and aggregation behaviour of dicationic and tricationic ionic liquid-based surfactants in aqueous solution Colloid. Surf. A: Physicochem. Eng. Asp. 2015 469 224 234 10.1016/j.colsurfa.2015.01.026
Nacham, O. et al. Interfacial and aggregation behaviour of dicationic and tricationic ionic liquid-based surfactants in aqueous solution. Colloid. Surf. A: Physicochem. Eng. Asp. 469, 224–234. 10.1016/j.colsurfa.2015.01.026 (2015).10.1016/j.colsurfa.2015.01.026
22. Shaban MM Eid AM Farag RK Negm NA Fadda AA Migahed MA Novel trimeric cationic pyridinium surfactants as bi-functional corrosion inhibitors and antiscalants for API 5L × 70 carbon steel against oilfield formation water J. Mol. Liquids 2020 305 112817 112839 10.1016/j.molliq.2020.112817
Shaban, M. M. et al. Novel trimeric cationic pyridinium surfactants as bi-functional corrosion inhibitors and antiscalants for API 5L × 70 carbon steel against oilfield formation water. J. Mol. Liquids 305, 112817–112839. 10.1016/j.molliq.2020.112817 (2020).10.1016/j.molliq.2020.112817
23. Yoshimura T Kusano T Iwase H Shibayama M Ogawa T Kurata H Star-shaped trimeric quaternary ammonium bromide surfactants: Adsorption and aggregation properties Langmuir 2012 28 9322 9331 10.1021/la301220y 22642550
Yoshimura, T. et al. Star-shaped trimeric quaternary ammonium bromide surfactants: Adsorption and aggregation properties. Langmuir 28, 9322–9331. 10.1021/la301220y (2012).22642550 10.1021/la301220y
24. El-Tabei AS Hegazy MA Bedair AH Sadeq MA Synthesis and inhibition effect of a novel tri-cationic surfactant on carbon steel corrosion in 0.5M H2SO4 solution J. Surfact. Deterg. 2014 17 341 352 10.1007/s11743-013-1524-7
El-Tabei, A. S., Hegazy, M. A., Bedair, A. H. & Sadeq, M. A. Synthesis and inhibition effect of a novel tri-cationic surfactant on carbon steel corrosion in 0.5M H2SO4 solution. J. Surfact. Deterg. 17, 341–352. 10.1007/s11743-013-1524-7 (2014).10.1007/s11743-013-1524-7
25. Rosen MJ Surfactants and Interfacial Phenomena 2004 3 Wiley
Rosen, M. J. Surfactants and Interfacial Phenomena 3rd edn. (Wiley, 2004).
26. Taleb K Mohamed-Benkada M Benhamed N Saidi-Besbes S Grohens Y Derdour A Benzene ring containing cationic gemini surfactants: Synthesis, surface properties and antimicrobial activity J. Mol. Liquids 2017 241 81 90 10.1016/j.molliq.2017.06.008
Taleb, K. et al. Benzene ring containing cationic gemini surfactants: Synthesis, surface properties and antimicrobial activity. J. Mol. Liquids 241, 81–90 (2017).10.1016/j.molliq.2017.06.008
27. Ghumare AK Pawar BV Bhagwat SS Synthesis and antibacterial activity of novel amido-amine-cationic gemini surfactants J. Surfact. Deterg. 2013 16 85 93 10.1007/s11743-012-1406-4
Ghumare, A. K., Pawar, B. V. & Bhagwat, S. S. Synthesis and antibacterial activity of novel amido-amine-cationic gemini surfactants. J. Surfact. Deterg. 16, 85–93 (2013).10.1007/s11743-012-1406-4
28. Mozrzymas A On the hydrophobic chains effect on critical micelle concentration of cationic gemini surfactants using molecular connectivity indices Monatshefte Chem. Chem. Mon. 2020 151 525 531 10.1007/s00706-020-02581-x
Mozrzymas, A. On the hydrophobic chains effect on critical micelle concentration of cationic gemini surfactants using molecular connectivity indices. Monatshefte Chem. Chem. Mon. 151, 525–531. 10.1007/s00706-020-02581-x (2020).10.1007/s00706-020-02581-x
29. Zana R Dimeric (gemini) surfactants: Effect of the spacer group on the association behavior in aqueous solution J. Colloid Interface Sci. 2002 248 203 220 10.1006/jcis.2001.8104 16290524
Zana, R. Dimeric (gemini) surfactants: Effect of the spacer group on the association behavior in aqueous solution. J. Colloid Interface Sci. 248, 203–220. 10.1006/jcis.2001.8104 (2002).16290524 10.1006/jcis.2001.8104
30. Wettig SD Verrall RE Thermodynamic studies of aqueous m-s-m gemini surfactant system J. Colloid. Interf. Sci. 2001 235 310 316 10.1006/jcis.2000.7348
Wettig, S. D. & Verrall, R. E. Thermodynamic studies of aqueous m-s-m gemini surfactant system. J. Colloid. Interf. Sci. 235, 310–316. 10.1006/jcis.2000.7348 (2001).10.1006/jcis.2000.7348
31. Zhang Z Zheng P Guo Y Yang Y Chen Z Wang X An X Shen W The effect of the spacer rigidity on the aggregation behaviour of two ester-containing Gemini surfactants J. Colloid Interface Sci. 2012 379 64 71 10.1016/j.jcis.2012.04.052 22608148
Zhang, Z. et al. The effect of the spacer rigidity on the aggregation behaviour of two ester-containing Gemini surfactants. J. Colloid Interface Sci. 379, 64–71. 10.1016/j.jcis.2012.04.052 (2012).22608148 10.1016/j.jcis.2012.04.052
32. Laschewsky A Lunkenheimer K Rakotoaly RH Wattebled L Spacer effects in dimeric cationic surfactants Colloid Polym. Sci. 2005 283 469 467 10.1007/s00396-004-1219-8
Laschewsky, A., Lunkenheimer, K., Rakotoaly, R. H. & Wattebled, L. Spacer effects in dimeric cationic surfactants. Colloid Polym. Sci. 283, 469–467. 10.1007/s00396-004-1219-8 (2005).10.1007/s00396-004-1219-8
33. Wettig SD Nowak P Verrall RE Thermodynamic and aggregation properties of gemini surfactants with hydroxyl substituted spacers in aqueous solution Langmuir 2002 18 5354 5359 10.1021/la011782s
Wettig, S. D., Nowak, P. & Verrall, R. E. Thermodynamic and aggregation properties of gemini surfactants with hydroxyl substituted spacers in aqueous solution. Langmuir 18, 5354–5359. 10.1021/la011782s (2002).10.1021/la011782s
34. Morita T Yada S Yoshimura T Adsorption dynamics of quaternary-ammonium-salt-based gemini and trimeric surfactants with different spacer structures at air/water interface Colloid. Surf. A: Physicochem. Eng. Asp. 2023 671 131589 1 11 10.1016/j.colsurfa.2023.131589
Morita, T., Yada, S. & Yoshimura, T. Adsorption dynamics of quaternary-ammonium-salt-based gemini and trimeric surfactants with different spacer structures at air/water interface. Colloid. Surf. A: Physicochem. Eng. Asp. 671(131589), 1–11. 10.1016/j.colsurfa.2023.131589 (2023).10.1016/j.colsurfa.2023.131589
35. Murguia MC Cristaldi MD Porto A Conza J Grau RJ Synthesis, surface-active properties and antimicrobial activities of new neutral and cationic trimeric surfactants J. Surfact. Deterg. 2008 11 41 48 10.1007/s11743-007-1052-4
Murguia, M. C., Cristaldi, M. D., Porto, A., Conza, J. & Grau, R. J. Synthesis, surface-active properties and antimicrobial activities of new neutral and cationic trimeric surfactants. J. Surfact. Deterg. 11, 41–48. 10.1007/s11743-007-1052-4 (2008).10.1007/s11743-007-1052-4
36. Xu H Chen D Cui Z Study on the synthesis and surface-active properties of a novel surfactant with triple quaternary ammonium groups and triple dodecyl chains derived from glycerin J. Surfact. Deterg. 2011 14 167 172 10.1007/s11743-010-1213-8
Xu, H., Chen, D. & Cui, Z. Study on the synthesis and surface-active properties of a novel surfactant with triple quaternary ammonium groups and triple dodecyl chains derived from glycerin. J. Surfact. Deterg. 14, 167–172. 10.1007/s11743-010-1213-8 (2011).10.1007/s11743-010-1213-8
37. Randić M On characterization of molecular branching J. Am. Chem. Soc. 1975 97 6609 6615 10.1021/ja00856a001
Randić, M. On characterization of molecular branching. J. Am. Chem. Soc. 97, 6609–6615. 10.1021/ja00856a001 (1975).10.1021/ja00856a001
38. Pedretti A Villa L Vistoli G VEGA—An open platform to develop chemo-bio-informatics applications using plug-in architecture and script programming J. Comput. Aided Mol. Des. 2004 18 167 173 10.1023/B:JCAM.0000035186.90683.f2 15368917
Pedretti, A., Villa, L. & Vistoli, G. VEGA—An open platform to develop chemo-bio-informatics applications using plug-in architecture and script programming. J. Comput. Aided Mol. Des. 18, 167–173. 10.1023/B:JCAM.0000035186.90683.f2 (2004).15368917 10.1023/B:JCAM.0000035186.90683.f2
39. Pedretti A Mazzolari A Gervasoni S Fumagalli L Vistoli G The VEGA suite of programs: An versatile platform for cheminformatics and drug design project Bioinformatics 2021 37 8 1174 1175 10.1093/bioinformatics/btaa774 33289523
Pedretti, A., Mazzolari, A., Gervasoni, S., Fumagalli, L. & Vistoli, G. The VEGA suite of programs: An versatile platform for cheminformatics and drug design project. Bioinformatics 37(8), 1174–1175. 10.1093/bioinformatics/btaa774 (2021).33289523 10.1093/bioinformatics/btaa774
40. Yousefinejad S Hemmateenejad B Chemometrics tools in QSAR/QSPR studies: A historical perspective Chemom. Intell. Lab. Syst. 2015 149 177 204 10.1016/j.chemolab.2015.06.016
Yousefinejad, S. & Hemmateenejad, B. Chemometrics tools in QSAR/QSPR studies: A historical perspective. Chemom. Intell. Lab. Syst. 149, 177–204. 10.1016/j.chemolab.2015.06.016 (2015).10.1016/j.chemolab.2015.06.016
41. Xu J Zhu L Fang D Liu L Wang L Xu W Prediction of dielectric dissipation factors of polymers from cyclic dimmer structure using multiple linear regression and support vector machine Colloid Polym. Sci. 2013 291 551 561 10.1007/s00396-012-2743-6
Xu, J. et al. Prediction of dielectric dissipation factors of polymers from cyclic dimmer structure using multiple linear regression and support vector machine. Colloid Polym. Sci. 291, 551–561. 10.1007/s00396-012-2743-6 (2013).10.1007/s00396-012-2743-6
42. Singh J Shaik B Singh S Agraval VK Khadikar PV Deeb O Supuran CT Comparative QSAR study on para-substituted aromatic sulphonamides as CAII inhibitors: Information versus topological (distance-based and connectivity) indices Chem. Biol. Drug. Des. 2008 71 244 259 10.1111/j.1747-0285.2007.00625.x 18221308
Singh, J. et al. Comparative QSAR study on para-substituted aromatic sulphonamides as CAII inhibitors: Information versus topological (distance-based and connectivity) indices. Chem. Biol. Drug. Des. 71, 244–259. 10.1111/j.1747-0285.2007.00625.x (2008).18221308 10.1111/j.1747-0285.2007.00625.x
43. Roy PP Roy K On some aspects of variable selection for partial least squares regression models QSAR Comb. Sci. 2008 27 302 313 10.1002/qsar.200710043
Roy, P. P. & Roy, K. On some aspects of variable selection for partial least squares regression models. QSAR Comb. Sci. 27, 302–313. 10.1002/qsar.200710043 (2008).10.1002/qsar.200710043
44. STATISTICA (data analysis software system). version 12 StatSoft. Inc.; software. http://www.statsoft.com (2014).
45. Klevens HB Structure and aggregation in dilute solutions of surface active agents J. Am. Oil Chem. Soc. 1953 30 74 80 10.1007/BF02635002
Klevens, H. B. Structure and aggregation in dilute solutions of surface active agents. J. Am. Oil Chem. Soc. 30, 74–80. 10.1007/BF02635002 (1953).10.1007/BF02635002
