William R. W. Welch

Effect of Asphaltene Structure on Association and Aggregation Using Molecular Dynamics

The aggregation of asphaltenes has been established for decades by numerous experimental techniques; however, very few studies have been performed on the association free energy and asphaltene aggregation in solvents. The lack of reliable and coherent data on the free energy of association and aggregation size of asphaltene has imposed severe limitations on the thermodynamic modeling of asphaltene phase behavior. Current thermodynamic models either consider asphaltenes as non-associating components or use fitting parameters to characterize the association. In this work, the relations between Gibbs free energy of asphaltene association and asphaltene molecular structure are studied using molecular dynamics (MD). The free energy of association is computed from the potential of mean force profile along the separation distance between the centers of mass of two asphaltene molecules using the umbrella sampling technique in the GROMACS simulation package. The average aggregation number for asphaltene nanoaggregates and clusters is also calculated through MD simulations of 36 asphaltene molecules in toluene and heptane in order to estimate the effects of association free energy and steric repulsion on the aggregation behavior of asphaltenes. Our simulation results confirm that the interactions between aromatic cores of asphaltene molecules are the major driving force for association as the energy of association increases substantially with the number of aromatic rings. Moreover, heteroatoms attached to the aromatic cores have much more influence on the association free energy than to ones attached to the aliphatic chains. The length and number of aliphatic chains do not seem to have a noticeable effect on asphaltene dimerization; however, they have a profound effect on asphaltene aggregation size since steric repulsion can prevent asphaltenes from forming T-shape configurations and therefore decrease the aggregation size of asphaltenes significantly. Our MD simulation results show for the first time asphaltene precipitation in heptane as an explicit solvent, and predict three distinct stages of aggregation (nanoaggregation, clustering, and flocculation) as proposed by the modified Yen model. Finally, the association free energy for asphaltenes in heptane is higher than that in toluene, which is consistent with asphaltene aggregate sizes obtained from MD simulations.

Infrared, vibrational circular dichroism, and Raman spectral simulations for β-sheet structures with various isotopic labels, interstrand, and stacking arrangements using density functional theory.

Infrared (IR), Raman, and vibrational circular dichroism (VCD) spectral variations for different β-sheet structures were studied using simulations based on density functional theory (DFT) force field and intensity computations. The DFT vibrational parameters were obtained for β-sheet fragments containing nine-amides and constrained to a variety of conformations and strand arrangements. These were subsequently transferred onto corresponding larger β-sheet models, normally consisting of five strands with ten amides each, for spectral simulations. Further extension to fibril models composed of multiple stacked β-sheets was achieved by combining the transfer of DFT parameters for each sheet with dipole coupling methods for interactions between sheets. IR spectra of the amide I show different splitting patterns for parallel and antiparallel β-sheets, and their VCD, in the absence of intersheet stacking, have distinct sign variations. Isotopic labeling by (13)C of selected residues yields spectral shifts and intensity changes uniquely sensitive to relative alignment of strands (registry) for antiparallel sheets. Stacking of multiple planar sheets maintains the qualitative spectral character of the single sheet but evidences some reduction in the exciton splitting of the amide I mode. Rotating sheets with respect to each other leads to a significant VCD enhancement, whose sign pattern and intensity is dependent on the handedness and degree of rotation. For twisted β-sheets, a significant VCD enhancement is computed even for sheets stacked with either the same or opposite alignments and the inter-sheet rotation, depending on the sense, can either further increase or weaken the enhanced VCD intensity. In twisted, stacked structures (without rotation), similar VCD amide I patterns (positive couplets) are predicted for both parallel and antiparallel sheets, but different IR intensity distributions still enable their differentiation. Our simulation results prove useful for interpreting experimental vibrational spectra in terms of β-sheet and fibril structure, as illustrated in the accompanying paper.

Insight into the packing pattern of β2 fibrils: A model study of glutamic acid rich oligomers with 13C isotopic edited vibrational spectroscopy

Polyglutamic acid at low pH forms aggregates and self-assembles into a spiral, fibril-like superstructure formed as a β2-type sheet conformation that has a more compact intersheet packing than commonly found. This is stabilized by three-centered bifurcated hydrogen bonding of the amide carbonyl involving the protonated glutamic acid side chain. We report vibrational spectroscopic results and analyses for oligopeptides rich in glutamic acid enhanced with (13)C isotope labeling in a study modeling low pH poly-Glu self-assembly. Our results indicate bifurcated H-bonding and β2 aggregation can be attained in these model decamers, confirming they have the same conformations as poly-Glu. We also prepared conventional β1-sheet aggregates by rapid precipitation from the residual peptides in the higher pH supernatant. By comparing the isotope-enhanced IR and VCD spectra with theoretical predictions, we deduced that the oligo-Glu β2 structure is based on stacked, twisted, antiparallel β-sheets. The best fit to theoretical predictions was obtained for the strands being out of register, sequentially stepped by one residue, in a ladder-like fashion. The alternate β1 conformer for this oligopeptide was similarly shown to be antiparallel but was less ordered and apparently had a different registry in its aggregate structure.

Molecular Dynamics Simulations of CO2/Water/Quartz Interfacial Properties: Impact of CO2 Dissolution in Water

The safe trapping of carbon dioxide (CO2) in deep saline aquifers is one of the major concerns of CO2 sequestration. The amount of capillary trapping is dominated by the capillary pressure of water and CO2 inside the reservoir, which in turn is controlled by the interfacial tension (IFT) and the contact angle (CA) of CO2/water/rock systems. The measurement of IFT and CA could be very challenging at reservoir conditions, especially in the presence of toxic cocontaminants. Thus, the ability to accurately predict these interfacial properties at reservoir conditions is very advantageous. Although the majority of existing molecular dynamics (MD) studies of CO2/water/mineral systems were able to capture the trends in IFT and CA variations with pressure and temperature, their predictions often deviated from experimental data, possibly due to erroneous models and/or overlooked chemical reactions. The objective of this study was to improve the MD predictions of IFT and CA of CO2/water/quartz systems at various pressure and temperature conditions by (i) considering the chemical reactions between CO2 and water and (ii) using a new molecular model for α-quartz surface. The results showed that the presence of carbonic acid at the CO2/water interface improved the predictions of IFT, especially at low temperature and high pressure where more CO2 dissolution occurs. On the other hand, the effect on CA was minor. The slight decrease in CA observed across the pressure range investigated could be attributed to an increase in the total number of H-bonds between fluid molecules and quartz surface.