Background {#Sec1}
==========

Diesel fuel has a vital function in the transportation sector, yet its combustion emits greenhouse gases and it is a finite resource; a cost-effective renewable substitute should have equivalent fuel efficiency, small net carbon emission, and be readily available worldwide. A potentially attractive alternative to fossil fuel is the use of plant oils---biodiesel \[[@CR1]\]. Biodiesel is defined as mono-alkyl esters of long-chain fatty acids (FAs), offering a viable alternative to petroleum-based diesel fuel. Biodiesel is non-toxic, degrades four times faster than diesel, and its blending with diesel fuel increases engine efficiency. It also does not produce net greenhouse effects and is safer in storage due to its high flash point \[[@CR2]\]. For these and other reasons, biodiesel production has gradually grown in recent years, raising the need for new rapid and cost-effective analytical tools and technologies for developmental characterization and quality control. ^1^H low-field nuclear magnetic resonance (LF-NMR) holds good potential in the fuel industry with many applications including determination of physical, chemical, structural, and dynamic properties.

In our previous work \[[@CR3]\], the molecular packing of methyl oleate (18:1) in its liquid phase was studied. The results from X-ray, LF-NMR diffusiometry, and high-field (HF)-NMR measurements were rationalized for the first time by a given liquid-phase packing model of fatty acid methyl esters (FAMEs; biodiesel). The developed model for 18:1 proposed that liquid FAMEs have short range order, where molecules arrange in a head-to-head conformation due to polar interactions, and their aggregate morphology retains a quasi-smectic liquid order mainly due to weak intermolecular interactions between two adjacent backbone molecules. The two molecule chains arrange longitudinally, and alternately to make an interdigitated structure, where in the same lateral plane, the ester groups of one molecule and the terminal methyl groups of the neighboring molecule are aligned side by side, similar to FAs \[[@CR4], [@CR5]\]. The translational movement of FAMEs, on the other hand, was found to differ from that of FAs, which are mostly dimerized due to head group hydrogen bonding. While the FA dimer is the basic unit of inter- and intra-molecular movements \[[@CR4], [@CR5]\], FAME molecules diffuse as monomers.

Biodiesel is a mixture of different FAMEs, with 18:1 as one of the main components, and each component affects the biodiesel properties as a function of its relative concentration. Therefore, to fully characterize biodiesel properties, it is necessary to further investigate different FAMEs that vary in chain lengths and number of double bonds.

Very little research has been performed on the liquid-phase molecular organization of FAMEs. In general, the physical properties of FAs and their biodiesel derivatives are largely determined by the length of the hydrocarbon chain, the degree of unsaturation, and the effect of molecular packing. In the fully saturated compounds, the hydrocarbon chain is highly flexible with free rotation around each carbon--carbon bond. The most stable conformation is the completely stretched arrangement, wherein the steric hindrance of neighboring atoms is minimized. This conformation allows tight packing in nearly crystalline arrays, stabilized by van der Waals force interactions between the atoms of neighboring molecules. In unsaturated hydrocarbons, chain bending occurs due to *cis* double bonds. The intermolecular interactions of molecules with alkyl chains having several *cis* bond bends are weaker than molecules with only one *cis*, because of less compact packing compared to one *cis* bond bend or a fully saturated hydrocarbon. These loosely ordered arrays of unsaturated molecules have lower melting points than monounsaturated molecules of the same chain length, because less thermal energy is needed for overcoming molecular interactions \[[@CR6]\].

Matsuzawa et al. \[[@CR7]\] studied different molecular packing densities and suggested that the existence of aggregate clusters most likely determines the liquid properties of FAs such as density and fluidity. It is feasible to assume that this is similar for FAMEs; hence, the liquid morphological structure of these materials will affect the physicochemical properties of the biodiesel including viscosity, density, fluid dynamics, and low temperature operability. These properties are of great importance to the field of biodiesel and each is based on different liquid characteristics. For example, viscosity, defined as a liquid's resistance to flow, is a function of the intermolecular forces of attraction within a liquid. Density, on the other hand, defined as the mass per unit volume of a material, depends on how tightly the molecules are packed together. The former rely on interactions between one molecule to its neighbor, while the latter depends on the conformation of a molecule and its packing density properties.

It has been shown that ^1^H LF-NMR spin-spin (*T*~2~) relaxometry can be applied to differentiate between morphological populations in complex systems \[[@CR8]--[@CR12]\]. Still, there isn't a certainty about the origin of triacylglycerol peaks in ^1^H LF-NMR relaxation time distributions. It has been suggested that the bimodal *T*~2~ distribution of liquid tricaprin is due to inhomogeneous relaxation rates for the protons along the side chains or inhomogeneous organization of the triacylglycerols in the liquid with intermolecular interactions \[[@CR13]\]. In our previous work \[[@CR3]\], these hypotheses were considered regarding oleic acid (18:1 acid) and methyl oleate (18:1), two similar alkyl chain molecules with different head groups, and it was suggested that the two peaks are the result of two distinct mobility populations of the protons on the chain.

In the present study, in order to characterize FAME aggregate structures and how they influence viscosity, liquid density, and temperature effects, we focused on further exploring the assignment of the bimodal peaks in ^1^H LF-NMR *T*~2~ distributions for different FAME molecules, with similar head groups but different alkyl chain lengths (10 to 20 carbons) and degrees of unsaturation (0, 1, 2, and 3 double bonds). Our objective was to study the relationship between molecular organization versus segmental and translational movements of different FAMEs in their pure liquid phase using ^1^H LF-NMR relaxometry and supporting advanced technologies, including X-ray diffraction, ^1^H LF-NMR diffusiometry, and ^13^C HF-NMR.

Results and discussion {#Sec2}
======================

^1^H LF-NMR *T*~2~ distributions of FAMEs at 313 K {#Sec3}
--------------------------------------------------

Biodiesel is a complex mixture of FAMEs with different lengths, degrees of unsaturation, and composition. It is therefore important to understand the molecular organization versus segmental and translational movements of the separate FAMEs, in order to explain their behavior in the biodiesel blend. The FAMEs evaluated in this study and their literature data of melting temperatures (denoted as *T*~m~) are summarized in Table [1](#Tab1){ref-type="table"} \[[@CR14]--[@CR17]\]. Throughout this manuscript, FAMEs are identified by their structures.Table 1Summary of the FAMEs evaluated in this study and literature data of melting temperatures (*T* ~m~s)IUPAC nameCommon nameStructure*T* ~m~ \[K\]Methyl decanoateMethyl caprate10:0259.5 \[[@CR14]\]Methyl dodecanoateMethyl laurate12:0278.0 \[[@CR15]\]Methyl tetradecanoateMethyl myristate14:0292.0 \[[@CR16]\]Methyl hexadecanoateMethyl palmitate16:0303.0 \[[@CR16]\]Methyl octadecanoateMethyl stearate18:0312.1 \[[@CR16]\]Methyl eicosanoateMethyl arachidate20:0318.8 \[[@CR17]\]Methyl Z-9-hexadecenoateMethyl palmitoleate16:1238.9 \[[@CR14]\]Methyl Z-9-octadecenoateMethyl oleate18:1253.0 \[[@CR15]\]Methyl Z,Z-9,12-octadecadienoateMethyl linoleate18:2238.0 \[[@CR15]\]Methyl Z,Z,Z-9,12,15-octadecatrienoateMethyl linolenate18:3227.5 \[[@CR21]\]*T* ~m~ melting temperature

The combined ^1^H LF-NMR *T*~2~ distributions of some saturated and unsaturated FAMEs at 313 K are presented in Fig. [1](#Fig1){ref-type="fig"}. *T*~2~ distributions are arranged by increasing chain length and degree of unsaturation. Intrinsic *T*~2~ values and percentage relative contribution (RC) of the peaks are marked on each plot. All *T*~2~ distributions exhibit two distinct peaks at different *T*~2~ values and RC. These will be denoted as *P*~1~ and *P*~2~ with increasing *T*~2~ values (*T*~2,1~ and RC~1~ will therefore stand for intrinsic *T*~2~ value and RC of *P*~1~). FAME 18:3 was found to be very prone to oxidation. Almost immediately following the first measurements an additional peak at low *T*~2~ values appeared that increased over time (Additional file [1](#MOESM1){ref-type="media"}). This same event occurred with different fresh samples and at different temperatures. In this study, therefore, fresh samples were used whenever the oxidation peak exceeded an RC of 5 %, and only the two main peaks, not related to oxidation, are discussed.Fig. 1Combined ^1^H LF-NMR *T* ~2~ distributions of FAMEs at 313 K. Plots are arranged by increasing chain length (*bottom* to *top*) and increasing degrees of unsaturation. The relative contributions of each peak in relation to the other peak and intrinsic *T* ~2~ values are shown on each plot. FAMEs are referred by their structures

As shown, each FAME exhibits a slightly different *T*~2~ distribution, which can be explained by the differences in chemical structures. For the saturated esters, with increasing chain length RC~1~ increases in relation to RC~2~ (0, 4, 15, 33, and 48 % RC~1~ for 10 to 18 carbons) and both *T*~2,1~ and *T*~2,2~ values decrease. For the unsaturated esters, with increasing number of double bonds, again RC~1~ increases in relation to RC~2~ (18, 46, and 65 % RC~1~ for 1 to 3 double bonds, respectively) whereas *T*~2~ values of both peaks increase.

Following the peak assignments we previously suggested \[[@CR3]\]; it is proposed that FAMEs have reduced mobility with increasing chain length and saturation, and intra-molecularly more rigid parts with increasing chain length and degree of unsaturation. Interestingly, for the 18 carbon esters, when comparing the fully saturated to monounsaturated, RC~1~ decreases to a greater extent (48 versus 18 % RC~1~ for 18:0 and 18:1, respectively). This may suggest different molecular organizations and/or type of interactions for the saturated versus unsaturated FAMEs, since the morphology of the alkyl chains, intermolecular interactions, and absolute temperature of measurement contribute to the ^1^H LF-NMR *T*~2~ distributions of FAMEs.

In our previous work, we compared *T*~2~ distributions of 18:1 and 18:1 acid relative to the temperature of melting of each compound, and several similarities were suggested for the peaks of each standard, which strengthened their assignment to two distinct mobility populations of the protons on the chain. In the case of 18:1 at 288 K, the short *T*~2~ peak was assigned as the protons on the chain backbone and the larger *T*~2~ peak as the protons on both ends of the chain. Adam-Berret et al. \[[@CR13]\] suggested that liquid triglycerides with three similar fully saturated alkyl chains have comparable behavior for the same difference from their melting temperature and that the effect of chain length was attenuated with this representation. A similar observation was suggested for different chain length liquid alkanes \[[@CR18]\]. This similar behavior is the result of comparable structures and type of interactions between alkyl chains.

Prior to melting, the atoms in a solid have restricted molecular motion and are confined to vibrating about their mean positions within the lattice/morphology structure of the solid. By increasing the temperature of the solid matter, the amplitude of the molecular vibrations increases, until at a certain temperature, intermolecular bonds within the solid break, allowing for bond rotation, and the molecules become free to rotate and translate within the liquid volume. The temperature of transition for a specific substance is the melting point. This is determined by the strength of a crystal lattice, which in turn is controlled primarily by four factors: the nature and number of intermolecular forces, molecular symmetry and packing efficiency, and the conformational degrees of freedom of a molecule \[[@CR19]\].

Following this approach, in order to gain a deeper understanding of the molecular organization and rotational and translational motions of saturated FAMEs in the liquid phase, additional measurements in this work include comparisons according to melting temperatures of each FAME material. Unsaturated FAMEs, on the other hand, have different intermolecular interactions, molecular freedom, and packing, leading to more complicated molecular structures, and rotational and translational motions, which cannot be standardized comprehensively by melting point. These will therefore be compared at absolute temperatures of measurement.

X-ray measurements {#Sec4}
------------------

In order to explore the molecular arrangement of the saturated FAMEs in this present study, the long- and short-range spacings were determined using small angle X-ray scattering (SAXS) and X-ray diffraction (XRD), respectively, close to their melting points at *T*~m~ + 15 K (Table [2](#Tab2){ref-type="table"}). The X-ray diffraction spectra acquired resemble those of 18:1 acid and 18:1 methyl ester, as previously reported \[[@CR3]\]. Two bands at around 0.14 and 0.03 nm^−1^ were observed, which give a measure of the spacing between adjacent molecules (short-range spacing) and long-range spacing of the plane as a result of head groups packing of the aligned molecules, respectively. Due to technical specifications of the XRD instrument, short-range spacing measurements were only available above ambient temperature; therefore measurements of short (10 and 12 carbons) saturated FAMEs, whose *T*~m~ + 15 K are below 298 K were performed at 298 K. SAXS showed no peaks for 10:0 at *T*~m~ + 15 K. It was thus measured at *T*~m~ + 5 K.Table 2Short- and long-range spacing of saturated FAMEs measured at *T* ~m~ + 15 KFAME structureShort spacing \[nm\]Long spacing \[nm\]10:00.459^a^1.70^b^12:00.461^a^2.0814:00.4592.2716:00.4622.5318:00.4552.6420:00.4572.83*T* ~m~ melting temperature^a^Measurements performed at ambient temperature^b^Measurement performed at *T* ~m~ + 5 K

All the FAMEs measured exhibited similar short-range spacing. As anticipated, for the saturated FAMEs, longer long-range spacings were found with increasing chain length. The long-range spacing for 10:0, however, was expected to be approximately half that of 20:0. According to Table [2](#Tab2){ref-type="table"}, this was not the case. This same occurrence was previously presented for the long-range spacing of 8:0 acid \[[@CR7]\] and was explained as 8:0 acid molecules being only partly interdigitated compared to 18:0 acid molecules, which are completely interdigitated. We believe that the same arrangement applies also for the FAMEs in hand, since even though the head group is different, the structure of the tails is similar.

Above melting temperature, and in accordance with our previous work on 18:1 \[[@CR3]\], the results are explained in that the studied FAMEs possess a liquid crystal-like structure, consisting of randomly aggregated liquid clusters with void spaces between them. This kind of microstructural organization can be the result of structural memory coming from the solid structure of FAMEs molecules, where polar interactions exist between heads and van der Waals forces are between two adjacent backbones forming an interdigitated structure. Evidence for this microstructural arrangement can be viewed by differences in density (Table [3](#Tab3){ref-type="table"}). The reported data was compiled from the studies by Pratas et al. \[[@CR20], [@CR21]\] for an absolute temperature of 313 K. With increasing temperature, the density of saturated and monounsaturated FAMEs decreases linearly, as suggested by Knothe and Steidley \[[@CR22]\].Table 3Densities, *ρ*, of some FAMEs measured at 313 K (compiled from the literature)FAME structure*ρ* \[kg/m^3^\]10:0856.0 \[[@CR20]\]12:0853.9 \[[@CR20]\]14:0852.2 \[[@CR20]\]16:0850.8 \[[@CR20]\]18:0849.8 \[[@CR20]\]16:1853.8 \[[@CR21]\]18:1859.5 \[[@CR20]\]18:2871.5 \[[@CR20]\]18:3887.0 \[[@CR20]\]*ρ* density

According to this data, for saturated FAMEs, density slightly decreases with increasing chain length. For the same number of carbons, density increases with increasing number of double bonds. These density measurements suggest, contrary to our expectation from molecular translational studies, that 18:3 molecules are the most closely packed, while saturated FAMEs are the most loosely packed. Iwahashi and Kasahara \[[@CR23]\] found similar behavior for saturated and monounsaturated FAs and have explained this discrepancy by the form of cluster aggregation. They suggested that saturated molecules aggregate tightly to make rigid clusters. These clusters form a liquid morphology with many void spaces between them leading to an overall lower liquid density and, consequently, a large apparent molar volume. Unsaturated molecules, on the other hand, form similar but softer clusters and can aggregate closer together to form smaller and fewer void spaces, leading to increased density and molar volume.

The same holds for the FAMEs in this study. Apparently, chain length can have the same effect on density, where the shorter the chain, the less number of interactions and the softer the cluster leading to increased packing density of the cluster and liquid density. Interestingly, the decrease in density with increasing chain length is not uniform. For the saturated compounds at 313 K, density decreases by 2.1, 1.7, 1.4, and 1.0 kg/m^3^ for 10 to 18 carbons, respectively. This is true also for other temperatures according to Knothe and Steidley \[[@CR22]\] and may be explained by the degree of interdigitation of rods inside the clusters as previously discussed.

Self-diffusion coefficients {#Sec5}
---------------------------

The effect of the chemical structure of FAMEs on their translational movement is shown in Fig. [2a, b](#Fig2){ref-type="fig"} for measurement at a single temperature that was above the melt points for all the FAMEs and for normalized temperature differences according to the *T*~m~ of each FAME, respectively. Interestingly, these results show conflicting trends that can be traced back to the change in morphology and translational motion above the melting point. The self-diffusion coefficient, *D*, is related to a fluid's viscosity from the Stokes-Einstein equation. From the physicochemical point of view, viscosity is related to the resistance of a molecule to move/slide relative to another molecule. Therefore, viscosity must be closely correlated with the structural parameters of the fluid particles \[[@CR24]\].Fig. 2Self-diffusion coefficients, *D*, of FAMEs in response to temperature. *D* measurements are compared at **a** absolute temperature and **b** temperature distances from each sample's melting point (*T* ~m~) for all FAMEs. FAMEs are referred to by their structures. For the absolute temperatures (**a**), diffusion increases with decreasing chain length and for the same chain length increasing number of double bonds. For the normalized temperatures (**b**), saturated FAMEs exhibit the highest and almost similar translational movement, followed by the monounsaturated FAMEs (16:1, 18:1), the di- (18:2), and tri-unsaturated (18:3) ones. The conflicting trends can be traced back to the change in morphology and translational motion above the melting point

In Fig. [2a](#Fig2){ref-type="fig"} (the same temperature for all FAMEs), the shorter the chain of the fully saturated FAME, the faster it diffuses. This can be explained by fewer intermolecular bonds and by the partly interdigitated structure of the shorter FAMEs. For the same chain length, the more double bonds, the larger the *D*. These results correlate very well with the *T*~2~ relaxation distributions presented in Fig. [1](#Fig1){ref-type="fig"}. It is well known that kinematic viscosity increases with chain length and with increasing degree of saturation \[[@CR15]\]. The effect of unsaturation has been attributed to interference of the double bonds with the molecules' ordered structure by adding kinks to the chain. Ramirez-Verduzco et al. \[[@CR25]\] suggested that coil-like *cis* configuration hinders the interactive approach of the double bond carbon atoms with the double bond carbons of neighboring molecules. This means that the translational movement of the unsaturated FAMEs increases with the number of double bonds.

When comparing the translational motion of the FAMEs at similar distances from their melting point (Fig. [2b](#Fig2){ref-type="fig"}) different trends are observed. Apparently all saturated FAMEs, regardless of chain length, exhibit almost similar translational movement. This can be rationalized by their similar structures and interactions, resulting in similar aggregate morphologies. The increase in melting points with chain length (Table [1](#Tab1){ref-type="table"}) is accounted for by the increase in the number of van der Waals interactions. Therefore, at similar distances from the melting point, the differences in translational motion are attenuated by the relative quantities of thermal energy. In addition, saturated FAMEs exhibit the highest translational movement, followed by the monounsaturated FAMEs (16:1, 18:1) and then by the di- (18:2) and tri-unsaturated (18:3) ones. As suggested before, when comparing FAMsEs with different degrees of unsaturation, several additional variables need to be considered that can explain the opposite trend in *D*, including temperature, type, and number of interactions, and aggregate morphologies for different alkyl chain configurations.

Segmental motion {#Sec6}
----------------

Segmental motion through the reciprocal of the effective correlation time, 1/*τ*~c~, of each carbon can be calculated from the longitudinal relaxation time, *T*~1~, measured by ^13^C HF-NMR. The 1/*τ*~c~ values of some saturated and unsaturated FAMEs were measured at *T*~m~ + 15 K and 298 K, respectively (Fig. [3a, b](#Fig3){ref-type="fig"}). *T*~1~, is likely to be correlated to the movement of the carbon atoms, i.e., segmental motion (specifically, rotational tumbling and to a lesser extent translational and internal motion) in the molecule. Assignment of ^13^C chemical shifts to the appropriate peaks was performed according to the literature \[[@CR26]\]. The 1/*τ*~c~ value for C1 is not shown since the carbonyl carbon has no covalent bonded protons. According to the decrease in 1/*τ*~c~ values from both ends towards the center of the chain, it is suggested that both the tail and head in all FAMEs have a higher freedom of movement compared to the interior chain. Although some degree of polar interactions exists between two heads, 1/*τ*~c~ values suggest that this interaction is very weak.Fig. 3Segmental motion of some FAMEs. Segmental motion through the reciprocal of the effective correlation time, 1/*τ* ~c~, of the carbon atoms at different positions of **a** saturated FAMEs (10:0, 16:0, and 18:0) at *T* ~m~ (melting point) +15 K and **b** unsaturated FAMEs (18:1, 18:2, and 18:3) at 298 K. FAMEs are referred to by their structures. The structures of 18:0 and 18:1 in **a** and **b**, respectively, along with designation of carbon numbers are shown for reference

According to Fig. [3a](#Fig3){ref-type="fig"}, the longer the chain length, the more rigid carbons there are in the interior part of the FAME molecule chain. These results fit very well with the data presented in the literature \[[@CR27]\] for the segmental motion of 9:0 acid and 18:0 acid. This can be rationalized by an increase in van der Waals interactions with increasing chain length, as suggested also by the molecular organization as viewed from the X-ray results (partial versus full interdigitated structure for short- and long-chain FAMEs, respectively). The lower 1/*τ*~c~ values for all carbons in 10:0 compared to 16:0 and 18:0 can be attributed to the large difference in the absolute temperature of measurement of the three (273 K for 10:0 versus 318 and 328 K for 16:0 and 18:0, respectively).

For the unsaturated FAMEs at an absolute temperature of 298 K, the rotational movement of both ends of the chain increases with number of double bonds. This correlates with the self-diffusion coefficient of unsaturated FAMEs at absolute temperatures (Fig. [2a](#Fig2){ref-type="fig"}), since as suggested in the literature, the segmental movements at the end and near the end of the molecule are probably most important for the FAME molecules to find the spaces for their translational diffusion \[[@CR3], [@CR28]\]. Segmental motion of the head and tail is facilitated due to the bending of the molecule, which leads to reduced van der Waals interactions and higher degree of freedom. The motion of the double bond carbons, on the other hand, is considerably restricted (Fig. [3b](#Fig3){ref-type="fig"}). Pi et al. \[[@CR29]\] stated that the movement and bending of alkyl chain of 18:2 acid from C9 to C13 is more regulated compared with that in the alkyl chain of 18:1 acid due to the presence of an additional C = C bond. The decrease in segmental motion at the double bond position can therefore be attributed to the larger energy barrier to bond rotation as a result of the double bonds and stronger intermolecular interactions between the *π* electrons. This implies increased rigidity with an increased number of double bonds as previously suggested according to the relative contribution of the peaks in *T*~2~ distributions (Fig. [1](#Fig1){ref-type="fig"}).

^1^H LF-NMR *T*~2~ distributions at different temperatures {#Sec7}
----------------------------------------------------------

The combined ^1^H LF-NMR *T*~2~ distributions of saturated and unsaturated FAMEs are presented in Figs. [4a--e](#Fig4){ref-type="fig"} and [5a--c](#Fig5){ref-type="fig"} with increasing chain length and degree of unsaturation, respectively. For the saturated FAMEs, measurements are compared at similar temperature distances from their *T*~m~, whereas the results for the unsaturated FAMEs are shown at absolute temperatures of measurement. Intrinsic *T*~2~ values and percentage RC of each peak are marked on each plot.Fig. 4Combined ^1^H LF-NMR *T* ~2~ distributions of saturated FAMEs at different temperatures. Comparison between the *T* ~2~ distributions of the saturated **a** 10:0, **b** 12:0, **c** 14:0, **d** 16:0, and **e** 18:0 FAMEs at different temperature distances from their melting points (*T* ~m~). FAMEs are referred to by their structures. The relative contributions of each peak in relation to the other peak and intrinsic *T* ~2~ values are shown on each plotFig. 5Combined ^1^H LF-NMR *T* ~2~ distributions of unsaturated FAMEs at different temperatures. Comparison between the *T* ~2~ distributions of the unsaturated **a** 18:1, **b** 18:2, and **c** 18:3 FAMEs at different temperatures. FAMEs are referred to by their structures. The relative contributions of each peak in relation to the other peak and intrinsic *T* ~2~ values are shown on each plot

In accordance with our previous work \[[@CR3]\], for the saturated FAMEs, as temperature increases *T*~2,1~ and *T*~2,2~ shift to higher values and RC~1~ decreases in favor of RC~2~. For the saturated FAMEs, the effect of presenting results at specific distances from melting point versus one absolute temperature can be seen by comparing *T*~2~ distributions in Fig. [4a--e](#Fig4){ref-type="fig"}, at any one of the temperature distances, to Fig. [1](#Fig1){ref-type="fig"}. According to this normalization, certain constancy can be observed in *T*~2,1~ and *T*~2,2~, especially for 10:0 and 12:0, and for 14:0--18:0. RC~1~, on the other hand, increases with chain length in relation to RC~2~, as stated for the absolute single temperature comparison (Fig. [1](#Fig1){ref-type="fig"}). This increase can be explained by the addition of rigid parts of the molecule and the number of van der Waals interactions, as manifested from the segmental motion (Fig. [3a](#Fig3){ref-type="fig"}).

It was established in the past that mono-exponential *T*~2~ values of FAMEs correlate with their viscosity for the same temperature \[[@CR30], [@CR31]\]. Table [4](#Tab4){ref-type="table"} presents the mono-exponential *T*~2~s for the saturated FAMEs and dynamic viscosities, *η*, according to \[[@CR20]\], at *T*~m~ + 5 K and at 313 K for reference. Mono-exponential *T*~2~s at 313 K were calculated from the relaxation data used for the analyses presented in Fig. [1](#Fig1){ref-type="fig"}. As shown, saturated FAMEs exhibit very similar translational motion at similar distances from melting point compared to absolute temperatures, according to the dynamic viscosity values and mono-exponential *T*~2~s, even though their melting points differ to a great extent (Table [1](#Tab1){ref-type="table"}). These results are comparable with the self-diffusion coefficients presented in Fig. [2a, b](#Fig2){ref-type="fig"} for absolute and normalized temperatures, respectively. In addition, the mono-exponential *T*~2~s at *T*~m~ + 5 K in Table [4](#Tab4){ref-type="table"} resemble each other to a greater degree than *T*~2,1~ and *T*~2,2~ in Fig. [4](#Fig4){ref-type="fig"} for the same temperatures (for 10:0 to 18:0 a *T*~2~ range of 589--617 versus 644--761 ms for mono-exponential *T*~2~ and *T*~2,2~, respectively). This suggests that the overall translational motion of saturated FAMEs is not affected solely by the mobility of the different parts of the molecules, and that the RCs of the peaks indicate other molecular parameters (e.g., molecular interactions) that also play an important part.Table 4Mono-exponential, *T* ~2~, and dynamic viscosity, *η*, of saturated FAMEs at *T* ~m~ + 5 K and 313 K*T* ~m~ + 5 K313 KFAME structure*T* ~2~ \[ms\]*η* ^a^ \[mPa · s\]*T* ~2~ \[ms\]*η* ^a^ \[mPa · s\]10:0589na12061.4812:06014.079742.0814:06063.987742.8416:06124.216623.7518:06174.435604.99*na* data not available, *T* ~m~ melting temperature, *T* ~2~ spin-spin relaxation time, *η* dynamic viscosity^a^Dynamic viscosity values were taken from \[[@CR20]\]

The unsaturated FAMEs (Fig. [5a--c](#Fig5){ref-type="fig"}) also display an increase in *T*~2,1~ and *T*~2,2~ with temperature. The RCs of the peaks, on the other hand, exhibit different trends with temperature and number of double bonds. In general for the unsaturated FAMEs, RC~2~ increases with temperature in relation to RC~1~; however, this response is attenuated with increasing number of double bonds. For the temperature range 258--308 K, RC~2~ increases in the ranges 24--77, 42--53, and 33--35 % for 1--3 double bonds, respectively. The change in RC of the peaks with increasing temperature and unsaturation may reflect the number and type of weak intermolecular interactions.

In this study, our characterization of different FAME molecules by ^1^H LF-NMR relaxometry, diffusion, ^13^C HF-NMR, and X-ray methods gives a clear indication of their molecular morphology and intermolecular interactions. For the saturated molecules, the longer the chain length, the higher the melting point. This is due to their linear chain configuration, which allows for molecules to pack closely together, maximizing the number of van der Waals contacts. This is manifested by increasing RC~1~, reduced diffusion coefficient, and a decrease in the segmental motion of the backbone carbons. For the same chain length and increasing number of double bonds, melting temperatures decrease. Unsaturated molecules of *cis* configuration cannot pack as close due to bending of the chain at the double bond position, minimizing secondary interactions. As seen from their melting temperatures, they require much less energy for disordering the crystal structure and breaking intermolecular interactions to achieve molecular motion above the melt point. This leads to an increase in diffusion coefficients and segmental motion of the tail.

When going from the fully saturated to the monounsaturated 18-carbon molecule, a reduction in RC~1~ occurs in the ^1^H LF-NMR *T*~2~ distributions (Fig. [1](#Fig1){ref-type="fig"}). This can be explained as a sharp decrease in the number of van der Waals interactions due to bending of the 18:1 chain, as previously discussed. When increasing the number of double bonds, increase in π-π interactions takes place as revealed by an increase in RC~1~ as shown when going from 18:1 → 18:2 → 18:3. This can also explain the decrease with unsaturation of the RC~1~ response to increasing temperature, since it is well known that less energy (temperature) is required to break van der Waals interactions compared to π-π interactions.

These results strengthen the peak assignment suggested for the ^1^H LF-NMR *T*~2~ distributions of FAMEs, where the two peaks are the result of two distinct mobility populations of the protons on the chain affected by the molecular structure and weak intermolecular interactions, especially those between the backbone of two adjacent chains. In a previous paper \[[@CR32]\], we presented the ^1^H LF-NMR *T*~2~ distribution of a rapeseed biodiesel sample measured at 313 K. In that study, the biodiesel sample exhibited three peaks with intrinsic *T*~2~ values of 338, 671, and 1141 ms. Given that the main constituents in this biodiesel at decreasing ratios are 18:1, 18:2, 18:3, 16:0, and 18:0, the peaks can be assigned as the average contribution according to the three regions designated a, b, and c in Additional file [2](#MOESM2){ref-type="media"}. An interesting study would be to explore the change in ^1^H LF-NMR RC of the peaks in biodiesels from different sources, storage conditions and shelf-life, and different temperatures, to study dynamic processes, melting mechanisms, and structural organizations of alkyl chains, with important applications in the development of biodiesel fuels.

Conclusions {#Sec8}
===========

Both the peaks' assignments for ^1^H LF-NMR *T*~2~ distributions of FAMEs and the model for their liquid crystal-like structure/morphology in the liquid phase, used to rationalize the assignment, were confirmed in the present work. This morphology along with the number and type of interactions and temperature effects generated differences in translational and rotational movements of the molecules, which were monitored using the presented ^1^H LF-NMR methodology. The study of morphological structures within liquids and their response to temperature changes by ^1^H LF-NMR is a powerful tool and is supported by traditional methods of characterization, such as X-ray and ^13^C HF-NMR. This new application of ^1^H LF-NMR is of potentially great interest to the field of biodiesel and to other research and applied disciplines with the potential of studying numerous physicochemical- and molecular organizational-based properties, processes, and mechanisms of alkyl chains.

Materials and methods {#Sec9}
=====================

Materials {#Sec10}
---------

Pure samples (≥99 %) of methyl ester standards (Table [1](#Tab1){ref-type="table"}) were purchased from Sigma Aldrich and used without further purification. These FAMEs exhibit a wide range of melting points and consequently are in different states for the same absolute temperature \[[@CR33]\]. Some of the measurements were therefore compared according to a given temperature distance from melting point (*T* = *T*~m~ + *d*, where *T* is the actual temperature of measurement, *T*~m~ is the melting point, and *d* the temperature distance from *T*~m~). The melting temperatures used for the materials in this study are within ±2 K from the melting temperatures reported in the literature (Table [1](#Tab1){ref-type="table"}).

^1^H LF-NMR relaxometry and diffusiometry {#Sec11}
-----------------------------------------

Measurements were carried out on a 20-MHz minispec bench-top pulsed NMR analyzer (Bruker Analytic GmbH, Germany), equipped with a permanent magnet, and a 10-mm temperature controlled probe head. Prior to measurement, samples were heated from 193 K for a minimum 3 h and then allowed to equilibrate inside the instrument for 5 min. All measurements were performed on liquid standards (above melting point). Receiver gain was optimized for each temperature and sample.

Determination of spin-spin relaxation time constant (*T*~2~) was performed using a Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence \[[@CR34], [@CR35]\]. *τ* and recycle delay were between 0.4 and 1.5 s and 4 and 8 s, respectively. For all the analyses, 32 scans were accumulated and 8000 echoes were acquired. Data was acquired in magnitude mode due to better repeatability and stability of results and further analyzed using the PDCO inverse Laplace transform optimization algorithm with *a*~2~ = 0.5 as described in the literature \[[@CR36]\].

Mono-exponential fitting of the acquired CPMG raw data was performed by SPSS software (version 15.0, SPSS Inc.) using Eq. [1](#Equ1){ref-type=""}:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ s(t)=w{\mathrm{e}}^{-t/{T}_2} $$\end{document}$$where *s*(*t*) is the acquired signal at *t* time, *w* is the pre-exponential weighting factor, and *T*~2~ is the mono-exponential relaxation time constant for transverse relaxation.

The self-diffusion coefficient, *D*, was determined by the pulsed-field gradient spin echo (PFGSE) method \[[@CR37]\]. The pulse sequence was used with 16 scans, *τ* of 7.5 ms, and a recycle delay of 6 s. Typical gradient parameters were *Δ* of 7.5 ms, *δ* of 0.5 ms, time between the 90° pulse to the first gradient pulse of 1 ms, and *G* of 1.6 T/m. A double distilled water sample was used for calibration. *D* values of water were taken from \[[@CR38]\]. Each reported value is the average of a minimum of ten measurements.

High-field (HF) ^13^C-NMR relaxometry {#Sec12}
-------------------------------------

Measurements were performed on a BRUKER AVANCE III operating at 150 MHz for ^13^C. Prior to measurement, samples were heated for minimum 10 min and added to a 5-mm NMR tube. For lock signal, a closed 1-mm capillary tube, filled with D~2~O, was added to the sample. The non-spinning samples were allowed to equilibrate inside the instrument for 15 min after reaching the set temperature. Before each measurement, shimming was optimized using automated and manual procedures. The chemical shifts, in parts per millions (ppm), were obtained without reference signal. The spectra were obtained using eight scans and recycle delay of 120 s.

The longitudinal relaxation times, *T*~1~, were measured using the INVREC method \[[@CR39]\]. The calculations of *T*~1~ were carried out with the subroutine included in the TOPSPIN 3.2 software package.

^13^C HF-NMR spin--lattice relaxation of a protonated carbon is overwhelmingly dominated by dipole-dipole interactions with the attached protons \[[@CR28]\]. *T*~1~ is therefore related to the number of directly bonded hydrogen, N, and the effective correlation time, *τ*~c~, for the rotational movement of the carbon atoms in the object molecule. Thus, *T*~1~ is approximately given in terms of N and 1/*τ*~c~:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {\mathrm{T}}_1=\frac{{\mathsf{r}}_{\mathsf{CH}}^{\mathsf{6}}}{N{\hslash}^2{\gamma}_{\mathrm{C}}^2{\gamma}_{\mathrm{H}}^2}\cdot \frac{1}{\tau_{\mathrm{C}}} $$\end{document}$$where *ħ* is Planck's constant and *γ*~C~ and *γ*~H~ are the gyromagnetic ratios of ^13^C and ^1^H, respectively. Here, *r*~CH~ is the C--H distance, usually about 0.109 nm, and the reciprocal of the effective correlation time, 1/*τ*~c~, represents the magnitude of the segmental rotation for the carbon atom at a different position.

X-ray methods {#Sec13}
-------------

XRD and SAXS techniques were used for measuring the short- and long-range spacing between adjacent molecules, respectively.

XRD data was collected on Panalytical Empyrean Powder Diffractometer equipped with position sensitive (PSD) X'Celerator detector using Cu K~α~ radiation (*λ* = 0.154 nm) and operated at 40 kV and 30 mA. The usual Bragg-Brentano *θ*/2θ geometry was employed. *θ*/2θ scans were run during 15 min in a 2θ range of 2--35° with step equal to \~0.0167°.

SAXS measurements were performed on a SAXSLAB GANESHA 300-XL (Skovlunde, Denmark) instrument. Cu K~α~ (*λ* = 0.154 nm) radiation was generated by Genix 3D Cu-source (operated at 47 mV and 0.55 mA) with integrated Monochromator, three pinholes collimation, and two-dimensional Pilatus 300 K detector. The distance between the sample and detector was 350 mm. *q* range was between 0.0012 and 0.067 nm^−1^.

Additional files {#Sec14}
================

Additional file 1:**Combined** ^**1**^ **H LF-NMR** ***T*** ~**2**~ **distributions of 18:3 at increasing degrees of oxidation.** Measurements were performed at 293 K. Plots are denoted A, B, and C by increasing time of natural oxidation. The relative contribution (RC) of the additional peak was found to increase over time, as marked on each plot, while the ratio between RC~1~ and RC~2~ (RC of *P* ~1~ and *P* ~2~, respectively) was kept unchanged.Additional file 2:**Combined** ^**1**^ **H LF-NMR** ***T*** ~**2**~ **distributions of a rapeseed biodiesel sample and its main FAMEs at 313 K.** Plots are arranged by increasing chain length (bottom to top) and increasing degrees of unsaturation. Peaks are assigned to three regions (a, b, c) according to intrinsic *T* ~2~s. FAMEs are referred to by their structures.

FA

:   Fatty acid

FAME

:   Fatty acid methyl ester

HF-NMR

:   High-field nuclear magnetic resonance

LF-NMR

:   Low-field nuclear magnetic resonance

PFGSE

:   Pulsed-field gradient spin echo

RC

:   Relative contribution

SAXS

:   Small angle X-ray scattering

XRD

:   X-ray diffraction

**Competing interests**

The authors declare that they have no competing interests.

**Authors' contributions**

NM performed the ^1^H LF-NMR measurements and analyzed all the results. PB drafted the manuscript and contributed to the analysis of results. LAC and TBM performed the HF-NMR measurements. NM, PB, LAC, CL, and ZW contributed to the development of the experimental design and proofread the manuscript. ZW led and coordinated the overall project. All authors read and approved the final manuscript.

NM acknowledges support from the Substitutes for Oil Transportation fellowship of the Ministry of Science, Technology and Space, Israel. PB acknowledges support from the Women in Science scholarship of the Israel Ministry of Science and Technology. The authors would like to thank Dr. Dimitri Mogiliansk and Dr. Sharon Hazan from the Ilse Katz Institute for Nanoscale Science and Technology at Ben-Gurion University, for performing the X-ray measurements; and the Phyto-Lipid Biotechnology Lab (PLBL) members at Ben-Gurion University of the Negev for their contribution to this work.
