Ion aggregation in high salt solutions. VI. Spectral graph analysis of chaotropic ion aggregates

Graph Theoretical Description of Phase Transitions in Complex Multiscale Phases with Supramolecular Assemblies

Phase transitions are typically quantified using order parameters, such as crystal lattice distances and radial distribution functions, which can identify subtle changes in crystalline materials or high-contrast phases with large structural differences. However, the identification of phases with high complexity, multiscale organization and of complex patterns during the structural fluctuations preceding phase transitions, which are essential for understanding the system pathways between phases, is challenging for those traditional analyses. Here, it is shown that for two model systems- thermotropic liquid crystals and a lyotropic water/surfactant mixtures-graph theoretical (GT) descriptors can successfully identify complex phases combining molecular and nanoscale levels of organization that are hard to characterize with traditional methodologies. Furthermore, the GT descriptors also reveal the pathways between the different phases. Specifically, centrality parameters and node-based fractal dimension quantify the system behavior preceding the transitions, capturing fluctuation-induced breakup of aggregates and their long-range cooperative interactions. GT parameterization can be generalized for a wide range of chemical systems and be instrumental for the growth mechanisms of complex nanostructures.

On the Elusive Links between Solution Microstructure, Dynamics, and Solvation Thermodynamics: Demystifying the Path through a Bridge over Troubled Conjectures and Misinterpretations

We build a fundamentally based bridge between the solute-induced microstructural perturbation of the species environment and the dynamic as well as thermodynamic responses of the fluid system, regardless of the state conditions, composition, nature of the solvent, and either the magnitude or the type of solute-solvent intermolecular-interaction asymmetries. For that purpose, we advance a fluctuation-based solvation formalism of fluid mixtures to provide meaningful descriptors of solvation phenomena, the microstructural signatures of their solute-solvent intermolecular interaction asymmetry, and the thermodynamic manifestations linked to the solution nonideality. The rigorous foundations afford us to address some crucial issues frequently invoked in the literature including the microstructural perturbation domain, its proper identification and molecular-based meaning toward the interpretation of the solvation process, and the potential impact of the local differential behavior between anions and cations on the actual salt-induced perturbation of the solvent microstructure. Indeed, we link the precisely characterized species solvation behavior to fundamental thermodynamic residual-property relations, and the dynamics associated with either the viscous flow or diffusive behavior of the solvent, to finally illustrate their outcome with experimental data of aqueous electrolyte solutions from the available literature. Ultimately, this effort provides a highly desirable unambiguous identification of the cause-effect connections between the microstructurally perturbed domains and the experimentally measured macroscopic solvation properties, including their effect on the dynamics of the solvent environment. More importantly, it lends a well-established solvation framework to bridge rigorously the microstructural details of the mixture, its dynamics, and its solvation thermodynamics to enhance our understanding of well-defined ranked Hofmeister series, i.e., by avoiding ad hoc conjectures and unsupported microscopic interpretations of solvation phenomena.

Isolation, process optimisation and characterisation of the protein from the de-oiled cake flour of Madhuca latifolia

This work reports the isolation of the protein from the flour of an underutilised agro waste, a de-oiled cake of Madhuca latifolia using the bis (2-ethylehexyl) sodium sulfosuccinate salt reverse micelle and the characterisation of the protein through various techniques. The experimental conditions for the extraction were optimised using Box-Behnken design. The highest yield of the protein was achieved when the extraction parameters, i.e. KCl concentration, KCl amount, and pH of the medium, were 0.5 M, 1.25 ml, and 9.02, respectively. The experimental yield (75.56%) obtained under the optimised conditions matched extremely well with the predicted yield (75.19%). The analysis of the biochemical composition envisaged the occurrence of 2S albumin, 7S globulin, and 11S globulin as the major components in the protein. The X-ray diffraction pattern supported the β-sheets structure of the protein. The imaging of the protein through a scanning electron microscope revealed the shape and surface of the protein to be spherical and smooth, respectively. Thus, the protein isolate of the de-oiled cake flour of Madhuca latifolia could be utilised towards food product development and relevant fields.

Insights into the hydrogen bond network topology of phosphoric acid and water systems

Phosphoric acid and its mixtures with water are some of the best proton conducting materials known to science. Although the proton conductivity in pure phosphoric acid decreases upon external doping with excess H+ or OH-, the addition of water improves it substantially. A number of experimental and theoretical studies indicate that these systems form a very special case of hydrogen bond networks which not only facilitate fast proton transport but also show a number of other interesting properties such as glass forming ability. In this work, we present the molecular dynamics simulation results of the H3PO4-H2O system over the entire concentration range. The hydrogen bond networks were analyzed in terms of conventional microscopic as well as topological properties based on graph and network theory. The results show that the hydrogen bond network of H3PO4 is fundamentally different from that of H2O. On average, each phosphoric acid molecule tends to form more and stronger hydrogen bonds than water which leads to a much more connected and clustered network showing small-world properties which are absent in pure water. Moreover, these hydrogen bond network properties persist in the H3PO4-H2O mixtures as well, even at relatively high water contents. Finally, many of the physical properties such as molecular diffusion coefficients seem to be also intimately related to the network topological properties and follow similar trends with respect to system content. These results strongly indicate that many important properties such as proton transport in phosphoric acid and its aqueous systems are fundamentally related to their hydrogen bond network topology and might hold the key for their ultimate molecular understanding.

Weighted persistent homology for osmolyte molecular aggregation and hydrogen-bonding network analysis

It has long been observed that trimethylamine N-oxide (TMAO) and urea demonstrate dramatically different properties in a protein folding process. Even with the enormous theoretical and experimental research work on these two osmolytes, various aspects of their underlying mechanisms still remain largely elusive. In this paper, we propose to use the weighted persistent homology to systematically study the osmolytes molecular aggregation and their hydrogen-bonding network from a local topological perspective. We consider two weighted models, i.e., localized persistent homology (LPH) and interactive persistent homology (IPH). Boltzmann persistent entropy (BPE) is proposed to quantitatively characterize the topological features from LPH and IPH, together with persistent Betti number (PBN). More specifically, from the localized persistent homology models, we have found that TMAO and urea have very different local topology. TMAO is found to exhibit a local network structure. With the concentration increase, the circle elements in these networks show a clear increase in their total numbers and a decrease in their relative sizes. In contrast, urea shows two types of local topological patterns, i.e., local clusters around 6 Å and a few global circle elements at around 12 Å. From the interactive persistent homology models, it has been found that our persistent radial distribution function (PRDF) from the global-scale IPH has same physical properties as the traditional radial distribution function. Moreover, PRDFs from the local-scale IPH can also be generated and used to characterize the local interaction information. Other than the clear difference of the first peak value of PRDFs at filtration size 4 Å, TMAO and urea also shows very different behaviors at the second peak region from filtration size 5 Å to 10 Å. These differences are also reflected in the PBNs and BPEs of the local-scale IPH. These localized topological information has never been revealed before. Since graphs can be transferred into simplicial complexes by the clique complex, our weighted persistent homology models can be used in the analysis of various networks and graphs from any molecular structures and aggregation systems.

Persistent homology analysis of osmolyte molecular aggregation and their hydrogen-bonding networks

Dramatically different patterns can be observed in the topological fingerprints for hydrogen-bonding networks from two types of osmolyte systems.

Understanding the Structure of the Hydrogen Bond Network and Its Influence on Vibrational Spectra in a Prototypical Aprotic Ionic Liquid

Analysis of the hydrogen bond network in aprotic ionic liquid 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) has been performed based on structures obtained from ab initio or classical molecular dynamics simulations. Statistics of different donor and acceptor atoms and the amount of chelating or bifurcated bonds has been presented. Most of the hydrogen bonds in EMIM-TFSI are formed with oxygen atoms as hydrogen acceptors; and the most probable bifurcated bonds are those with a mixed pair of oxygen and nitrogen acceptors. Spectral graph analysis has shown that the cations may form hydrogen bonds with up to five different anions and the connectivity of the whole hydrogen bond network is supported mainly by H-O bonds. In the structures of the liquid simulated via force field-based dynamics, the number of hydrogen bonds is smaller and fluorine atoms are the most favored hydrogen acceptors. One-dimensional potential energy profiles for hydrogen atom displacements and corresponding vibrational frequencies have been calculated for selected C-H bonds. Individual C-H stretching frequencies vary by 200-300 cm^-1, indicating differences in local environment of hydrogen atoms forming C-H···O hydrogen bonds.

Persistent homology analysis of ion aggregations and hydrogen-bonding networks

Despite the great advancement of experimental tools and theoretical models, a quantitative characterization of the microscopic structures of ion aggregates and their associated water hydrogen-bonding networks still remains a challenging problem. In this paper, a newly-invented mathematical method called persistent homology is introduced, for the first time, to quantitatively analyze the intrinsic topological properties of ion aggregation systems and hydrogen-bonding networks. The two most distinguishable properties of persistent homology analysis of assembly systems are as follows. First, it does not require a predefined bond length to construct the ion or hydrogen-bonding network. Persistent homology results are determined by the morphological structure of the data only. Second, it can directly measure the size of circles or holes in ion aggregates and hydrogen-bonding networks. To validate our model, we consider two well-studied systems, i.e., NaCl and KSCN solutions, generated from molecular dynamics simulations. They are believed to represent two morphological types of aggregation, i.e., local clusters and extended ion networks. It has been found that the two aggregation types have distinguishable topological features and can be characterized by our topological model very well. Further, we construct two types of networks, i.e., O-networks and H2O-networks, for analyzing the topological properties of hydrogen-bonding networks. It is found that for both models, KSCN systems demonstrate much more dramatic variations in their local circle structures with a concentration increase. A consistent increase of large-sized local circle structures is observed and the sizes of these circles become more and more diverse. In contrast, NaCl systems show no obvious increase of large-sized circles. Instead a consistent decline of the average size of the circle structures is observed and the sizes of these circles become more and more uniform with a concentration increase. As far as we know, these unique intrinsic topological features in ion aggregation systems have never been pointed out before. More importantly, our models can be directly used to quantitatively analyze the intrinsic topological invariants, including circles, loops, holes, and cavities, of any network-like structures, such as nanomaterials, colloidal systems, biomolecular assemblies, among others. These topological invariants cannot be described by traditional graph and network models.

The effect of Hofmeister anions on water structure at protein surfaces

To understand the effects of specific ions on protein-water interactions and the thermodynamic stability of proteins in salt solutions, we use a molecular dynamics (MD) simulation to examine the water structure, orientational distribution, and dynamics near the surface of ubiquitin. In particular, we consider NaCl, NaBF₄, NaSCN, and NaClO₄ salt solutions containing ubiquitin, where the anions of the latter three salts are well-known chaotropic ions in the Hofmiester anion series. The number of hydrogen bonds (H-bonds) per water molecule is found to decrease significantly at the ubiquitin-water interface, indicating a significant disruption of the water H-bonding network. The distribution of the water H-bond numbers near the protein surface is modulated by dissolved ions, and the extent of the ion effect on the H-bonding network structure follows the order of the Hofmeister anion series, while there are no specific ion effects on water properties at distances larger than 5 Å from the protein surface. From detailed analyses of the surface area, volume, and root-mean-square deviation (RMSD) of ubiquitin, we show that changes in the properties of the protein could originate from the disruption of the water H-bond network induced by ions with a higher affinity for the protein surface instead of direct protein residue-ion interactions. An interesting observation made here is that the orientational distribution of water molecules at the protein-water interface is close to random, but there is a slight preference for interfacial water molecules with a straddle structure within 2.5 Å of the protein surface, where one of the two OH groups points away from the protein surface and the other points toward the surface. In addition, comparing the MD simulation results for ubiquitin solutions with dissolved NaSCN and KSCN, we show that Na⁺ affects the water H-bonding structure at the protein surface more than K⁺. It is clear that the H-bonding network structure of water more than one water layer away from the protein surface is not distinguishably different from that of neat water. We thus anticipate that the present work will provide insights into the scale of specific ion effects on the H-bonding structure and orientational distribution of water in the vicinity of protein surfaces in aqueous solutions.

Dynamic heterogeneity in aqueous ionic solutions

It is well known that supercooled liquids have heterogeneous dynamics, but it is still unclear whether dynamic heterogeneity also exists in aqueous ionic solutions at room or even higher temperatures.

Ultrafast fluxional exchange dynamics in electrolyte solvation sheath of lithium ion battery

Lithium cation is the charge carrier in lithium-ion battery. Electrolyte solution in lithium-ion battery is usually based on mixed solvents consisting of polar carbonates with different aliphatic chains. Despite various experimental evidences indicating that lithium ion forms a rigid and stable solvation sheath through electrostatic interactions with polar carbonates, both the lithium solvation structure and more importantly fluctuation dynamics and functional role of carbonate solvent molecules have not been fully elucidated yet with femtosecond vibrational spectroscopic methods. Here we investigate the ultrafast carbonate solvent exchange dynamics around lithium ions in electrolyte solutions with coherent two-dimensional infrared spectroscopy and find that the time constants of the formation and dissociation of lithium-ion···carbonate complex in solvation sheaths are on a picosecond timescale. We anticipate that such ultrafast microscopic fluxional processes in lithium-solvent complexes could provide an important clue to understanding macroscopic mobility of lithium cation in lithium-ion battery on a molecular level.