Solvation Energy Research Articles

Bound Charge Layering Near a Spherical Ion: Moving Beyond Born Theory.

The response of aqueous solvent to a dissolved ion is analyzed in terms of the bound charge, the net charge of the solvent in the vicinity of the solute. The total amount of bound charge is , where Qion is the charge of the ion and ϵ is the solvent dielectric constant, in both continuum and molecular theory. Aqueous solvation involves an inner layer of bound charge way over this value, followed by another layer that almost or over compensates the first layer, and so on. We demonstrate how layering of charge explains the strong solvation response of aqueous solvent. Born theory, in which the ion resides in a cavity within a dielectric continuum, places all the bound charge on the cavity surface. Without accounting for bound charge layering, it cannot describe the strong aqueous solvation response. The only adjustable parameter in Born theory is the cavity radius, and unphysically small values of the radius are required to match the Born prediction to the actual solvation free energy. We propose a simple analytical model for aqueous solvation of a spherical ion that incorporates bound charge layering. We point out which parameters are expected to be solvent-specific and transferable between different solutes, while other parameters should depend upon ion and solvent size. The solvation energy from finite, periodically replicated simulations must be corrected to describe an ion at infinitely dilution. We present a very simple correction, and demonstrate its accuracy.

Use of green RPLC methods as an alternative to conventional RPLC methods in the determination of the pKa values of some H2 receptor antagonists

In this study, the chromatographic retention behavior and dissociation constants (pKa ) of H2 receptor antagonists ranitidine and famotidine were determined by the green reverse phase liquid chromatography (RPLC) and hybrid micellar chromatography methods, which are widely used today. Replicate analyses were performed on a Gemini NX C18 (250 × 3.0 mm ID, 5 μm) column at constant temperature (37 °C) and flow rate (1 mL/min). The obtained data analyses were evaluated with the linear solvation energy relationship (LSER) approach and the consistency of the data was assessed by calculating the theoretical retention time values. The pK a 1 s s and pK a 2 s s values obtained in ethanol–water binary mixtures containing 11%, 13%, 15%, and 20% (v/v) ethanol were also calculated with this approach. The pKa values of these compounds in water and pure micellar environment were calculated with pK a s s values and solvent. It was observed that the p K a w w data obtained by three different methods were compatible. Accordingly, with the hybrid micelles method, p K a 2 w w and p K a 3 w w values for famotidine were calculated as 7.305 and 9.869, respectively, and p K a 1 w w and p K a 2 w w values for ranitidine were calculated as 2.952 and 8.589. The environmental impacts of the methods applied in the study were evaluated using green metric tools.

Synergistic Modeling of Liquid Properties: Integrating Neural Network-Derived Molecular Features with Modified Kernel Models.

A significant challenge in applying machine learning to computational chemistry, particularly considering the growing complexity of contemporary machine learning models, is the scarcity of available experimental data. To address this issue, we introduce an approach that derives molecular features from an intricate neural network-based model and applies them to a simpler conventional machine learning model that is robust to overfitting. This method can be applied to predict various properties of a liquid system, including viscosity or surface tension, based on molecular features drawn from the ab initio calculated free energy of solvation. Furthermore, we propose a modified kernel model that includes Arrhenius temperature dependence to incorporate theoretical principles and diminish extreme nonlinearity in the model. The modified kernel model demonstrated significant improvements in certain scenarios and possible extensions to various theoretical concepts of molecular systems.

Computational evaluation of phosphonium ILs as CO2 absorbents for electrochemistry

Since phosphonium-based ionic liquids (ILs) are suitable electrolytes for electrochemical reduction, a computationally informed analysis of these ILs was conducted to assess their ability to absorb CO2. In this work, theoretical and experimental approaches were mixed with a focus on gas solubility, transport properties, and structural characterization. The calculations included free energy of solvation (ΔG), the Henry constant (Kh), and short-range energies (Coulomb plus Lennard–Jones). Two ILs, [P1444][TFSI] and [P1444][FSI], were selected for further evaluation because of their lower Kh values. The enthalpy and entropy of the solvation process for these selected ILs were computationally determined; both ILs were found to have a favourable enthalpy of solvation for gas dissolution, and [P1444][TFSI] had a lower penalty for solvation (less negative entropy). Moreover, the analysis revealed competition among the ions interacting with the gas. The structural characterization of the [P1444][TFSI] and [P1444][FSI] with the CO2 systems involved radial and spatial distribution functions. The gas was likely enveloped by an anion cage, and oxygen and fluorine exhibited greater interactions with the carbon atom of CO2 than nitrogen. The viscosity and density were also calculated for the two systems, and the addition of CO2 only caused a slight change on these values. Experiments on viscosity and density confirmed the computational results; here, compared with those of neat ionic liquids, an approximate 4 % decrease in the viscosity of CO2-saturated systems was observed. Finally, the actual CO2 solubility in [P1444][TFSI] was experimentally determined to be 120 mM; this value was nearly four times greater than those of typical aqueous bicarbonate solutions.

5-methyl-N'-(thiophen-2-ylmethylene)-1H-pyrazole-3-carbohydrazide as potent anticancer agent: Synthesis, spectroscopic characterization, anticancer activity and DFT studies

A new pyrazole-hydrazide hydrazone derivative, named 5-methyl-N'-(thiophen-2-ylmethylene)-1H-pyrazole-3-carbohydrazide (TMPC) has been synthesized and characterized by FT-IR, UV–Vis, 1H NMR, 13C NMR, and HRMS-ESI. Its cytotoxic activity was tested on various cancer and healthy cell lines using the XTT method. In addition, the advanced anticancer mechanism was investigated by flow cytometry analysis. Hybrid B3LYP/6–311++G** calculations reveal two stable C1 and C2 conformers of TMPC with energy differences of 28.96 kJ/mol in the gas phase and 8.50 kJ/mol in ethanol solution. Comparisons among the experimental NMR, UV and infrared spectra with the predicted ones support the presence of both forms in solution. The presence of dimer C2 justify the intensities and numerous bands observed in the IR spectrum in the 2000–10 cm−1 region. Probably, the repulsion between the N7 and O1 atoms in C1 explain its lower stability and the higher solvation energy (-108.20 kJ/mol) observed for this form in ethanol than C2 (-88.01 kJ/mol). NBO calculations suggest that both forms are stable in the two media while the AIM analyses evidence the higher stability of C2 than C1 due to the N9-H22···N7 interaction observed only for this conformer in the gas phase. Gap values suggest high reactivities of both forms, as compared with similar pyrazole-carbohydrazide species. Full vibrational assignments are reported for both forms of TMPC together with the main scaled force constants. The XTT results showed that TMPC exhibited a selective cytotoxic effect on the tested cancer cell lines in a dose- and time-dependent manner, and also displayed the most potent cytotoxic activity against breast cancer MCF-7 cells with an IC50 value of 22.79 µM after 48 h Flow cytometric analyses showed that TMPC induced apoptosis and suppressed PI3K/Akt signaling in MCF-7 cells.

On delivering polar solvation free energy of proteins from energy minimized structures using a regularized super-Gaussian Poisson-Boltzmann model.

The biomolecules interact with their partners in an aqueous media; thus, their solvation energy is an important thermodynamics quantity. In previous works (J. Chem. Theory Comput. 14(2): 1020-1032), we demonstrated that the Poisson-Boltzmann (PB) approach reproduces solvation energy calculated via thermodynamic integration (TI) protocol if the structures of proteins are kept rigid. However, proteins are not rigid bodies and computing their solvation energy must account for their flexibility. Typically, in the framework of PB calculations, this is done by collecting snapshots from molecular dynamics (MD) simulations, computing their solvation energies, and averaging to obtain the ensemble-averaged solvation energy, which is computationally demanding. To reduce the computational cost, we have proposed Gaussian/super-Gaussian-based methods for the dielectric function that use the atomic packing to deliver smooth dielectric function for the entire computational space, the protein and water phase, which allows the ensemble-averaged solvation energy to be computed from a single structure. One of the technical difficulties associated with the smooth dielectric function presentation with respect to polar solvation energy is the absence of a dielectric border between the protein and water where induced charges should be positioned. This motivated the present work, where we report a super-Gaussian regularized Poisson-Boltzmann method and use it for computing the polar solvation energy from single energy minimized structures and assess its ability to reproduce the ensemble-averaged polar solvation on a dataset of 74 high-resolution monomeric proteins.

Solvatochromism, electric dipole moments, optical properties and electronic structure of biphenylcarbonitrile liquid crystals

The UV–Vis. absorbance and fluorescence spectra of 4HC (4′-heptyl-4-biphenyl carbonitrile) and 4OB (4′-Octyl-4-biphenylcarbonitrile) liquid crystals were measured in 23 solvents with different polarities. As dependent on changing of solvents, while the change in absorbance spectrum of liquid crystals have less, change in fluorescence spectrum have more. Molecular interactions were quantitavely investigated by using Linear solvation energy relationships. Optical forbidden energy gap has been determined with using Tauc plot. The refractive index has been calculated with semi-empirical methods. The relationships between refractive index and Eg has been commented in n-hexane, acetic acid, water and 1,4-dioxane. The electric dipole moment has been calculated by using the Bilot-Kawski, Lippert-Mataga, Bakshiev, Kawski-Chamma-Viallet and Reichardt methods. Solvatochromic shifts showed to having large dipole moments of 4OB and 4HC LCs. Their molecular orbital energies, electric dipole moments, MEP, SAS, and reactivity properties of 4HC and 4OB LCs has been found by quantum chemical calculations. It was observed that the dipole moment of 4HC was higher than 4OB.

Recognition of halide anions by N,N’-bis(4-nitrophenyl)urea in the gas-phase and solution: a computational study

ABSTRACT This study investigates the recognition of halide anions by N,N’-bis(4-nitrophenyl)urea using computational methods in both gas-phase and solution. It examines interaction, stabilization, and deformation energies, revealing a gas-phase selectivity trend of [L⋯F]−>[L⋯Cl]−>[L⋯Br]−>[L⋯I]−. NBO, QTAIM, and EDA-NOCV analyses highlight N-H∙∙∙X hydrogen bonds, with electrostatic forces contributing around 60% to the total interaction. Despite unfavorable solvation energy changes, the intrinsic affinity remains the main factor for selectivity, which aligns with experimental data. The strong correlation between theoretical and experimental formation constants supports the reliability of the computational findings.

Bandgap tuning for transition metal oxides via PEGylation

Abstract Bandgap engineering is controlled manipulation of the bandgap of materials/meta-materials to achieve desired properties. The electrical and optical properties of materials are significantly affected by bandgap tuning; therefore, bandgap engineering is a powerful technique for designing electronic and optoelectronic devices. Compositional engineering, strain engineering, and nanoscience and technology are the three major fields associated with bandgap engineering. Any unique combination of this engineering can provide novel strategies to produce novel band-structured devices. In this method article, we have demonstrated how solvation energy can alter the bandgap energy, a fact that is generally ignored due to misconceptions about quantum/size confinement. Here, we prepare nanostructured transition metal oxides (Co3O4, CuO, and ZnO) with polyethylene glycol (PEG), and the method is termed PEGylation. We investigate the influence of PEGylation on the structural, electrochemical, and electronic nature of these oxides. It is observed that the bandgap tunability (7.33%) is maximum for ZnO. Our study suggests that band alteration is significantly correlated with the change in lattice parameters; however, it is orientation dependent as the correlation coefficient reduces to 0.85 from 1 for the change in lattice parameter b along the y-axis compared to the other two lattice parameters. Similarly, band alteration is also known to have some correlation with the electrochemical potential, but is surprisingly almost independent of size confinement.

Designing Alkylammonium Cations for Enhanced Solubility of Anionic Active Materials in Redox Flow Batteries: The Role of Bulk and Chain Length.

Advancing grid-scale energy storage technologies is crucial for realizing a fully renewable energy landscape, with non-aqueous redox flow batteries (NRFBs) presenting a promising solution. One of the current challenges in NRFBs stems from the low energy density of redox active materials, primarily due to their limited solubility in non-aqueous solvents. Herein, this study explores the solubility of vanadium(IV/V) bis-hydroxyiminodiacetate (VBH) crystals in acetonitrile, aiming to use them as anionic catholytes in NRFBs. We focused on enhancing VBH solubility by modifying the structure of the alkylammonium cation. Employing periodic density functional theory and a solvation model, we calculated the dissolution free energy ([[EQUATION]]), which includes sublimation ([[EQUATION]]) and solvation ([[EQUATION]]) energies. Our results indicate that neither elongating straight-chain alkyl groups beyond a tetrabutylammonium baseline nor introducing bulky substituents at the nitrogen center significantly enhances solubility. However, the introduction of carbon spacers combined with terminal bulky substituents markedly improves solubility by favorably altering both [[EQUATION]] and [[EQUATION]]. These findings underline the nuanced impact of cation structure on solubility and suggest a viable approach to optimize VBH-based anionic catholytes. This advancement promises to enhance NRFB efficiency and sustainability, marking a significant step forward in energy storage technology.

Functionalization of single-walled aluminum nitride nanotube with amino acids using the first principle's study

In this study, we investigate the functionalization of aluminum nitride nanotube (AlNNT) with amino acids (Phenylalanine, Arginine, and Glutamine) using first-principle density functional theory (DFT) calculations. DFT methods were applied to model the interactions between amino acids and AlNNT, exploring the structural and electronic features of the functionalized systems, and examining various aspects of the structures, including their geometry and electronic properties, such as bond lengths, bond angles, total energy, formation energy, band structure, chemical potential, density of states (DOS), projected density of states (PDOS), charge transfer, and electron difference density. The (5, 5) AlNNT exhibits an energy bandgap value of 3.31 eV, implying its semiconducting nature. Furthermore, the pure system's chemical potential was -4.12 eV. The DOS examination shows that the pristine system's band structure profile and energy states match well. Notably, the valence band displays a more pronounced level of energy hybridization compared to the conduction band, which is reflected in the band structure. Functionalization of AlNNT with amino acids enhances the p character and polar covalency, thereby increasing the percentage ionic character and Gibb's free energy of solvation which are essential parameters to determine the aqueous solubility.

Quantum Chemical (QC) Calculations and Linear Solvation Energy Relationships (LSER): Hydrogen-Bonding Calculations with New QC-LSER Molecular Descriptors

Quantum Chemical (QC) Calculations and Linear Solvation Energy Relationships (LSER): Hydrogen-Bonding Calculations with New QC-LSER Molecular Descriptors

Transfer learning for molecular property predictions from small datasets

Machine learning has emerged as a new tool in chemistry to bypass expensive experiments or quantum-chemical calculations, for example, in high-throughput screening applications. However, many machine learning studies rely on small datasets, making it difficult to efficiently implement powerful deep learning architectures such as message passing neural networks. In this study, we benchmark common machine learning models for the prediction of molecular properties on two small datasets, for which the best results are obtained with the message passing neural network PaiNN as well as SOAP molecular descriptors concatenated to a set of simple molecular descriptors tailored to gradient boosting with regression trees. To further improve the predictive capabilities of PaiNN, we present a transfer learning strategy that uses large datasets to pre-train the respective models and allows us to obtain more accurate models after fine-tuning on the original datasets. The pre-training labels are obtained from computationally cheap ab initio or semi-empirical models, and both datasets are normalized to mean zero and standard deviation one to align the labels’ distributions. This study covers two small chemistry datasets, the Harvard Organic Photovoltaics dataset (HOPV, HOMO–LUMO-gaps), for which excellent results are obtained, and the FreeSolv dataset (solvation energies), where this method is less successful, probably due to a complex underlying learning task and the dissimilar methods used to obtain pre-training and fine-tuning labels. Finally, we find that for the HOPV dataset, the final training results do not improve monotonically with the size of the pre-training dataset, but pre-training with fewer data points can lead to more biased pre-trained models and higher accuracy after fine-tuning.

Predicting pKa Values of Para-Substituted Aniline Radical Cations vs. Stable Anilinium Ions in Aqueous Media

The focus of pKa calculations has primarily been on stable molecules, with limited studies comparing radical cations and stable cations. In this study, we comprehensively investigate models with implicit solvent and explicit water molecules, direct and indirect calculation approaches, as well as methods for calculating free energy, solvation energy, and quasi-harmonic oscillator approximation for para-substituted aniline radical cations (R-PhNH2•+) and anilinium cations (R-PhNH3+) in the aqueous phase. Properly including and positioning explicit H2O molecules in the models is important for reliable pKa predictions. For R-PhNH2•+, precise pKa values were obtained using models with one or two explicit H2O molecules, resulting in a root mean square error (RMSE) of 0.563 and 0.384, respectively, for both the CBS-QB3 and M062X(D3)/ma-def2QZVP methods. Further improvement was achieved by adding H2O near oxygen-containing substituents, leading to the lowest RMSE of 0.310. Predicting pKa values for R-PhNH3+ was more challenging. CBS-QB3 provided an RMSE of 0.349 and the M062X(D3)/ma-def2QZVP method failed to calculate pKa accurately (RMSE > 1). However, by adopting the double-hybrid functional method and adding H2O near the R substituent group, the calculations were significantly improved with an average absolute difference (ΔpKa) of 0.357 between the calculated and experimental pKa values. Our study offers efficient and reliable methods for pKa calculations of R-PhNH2•+ (especially) and R-PhNH3+ based on currently mature quantum chemistry software.

Quantitative Analysis of the Li-Ion Solvation Structure in Li-Ion Battery Electrolytes Using ATR-FTIR Spectroscopy.

Attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy is widely used to study condensed materials due to its convenient sample preparation and ability to avoid absorption saturation. Recently, it has been applied to in situ and in operando observations of chemical reactions within electrochemical devices, such as lithium-ion batteries. However, because ATR-FTIR spectroscopy relies on frequency-dependent attenuated reflectance, quantitative concentration measurements of chemical species using the Beer-Lambert law are challenging. Despite the availability of several correction methods, discrepancies remain in the solvation structures around Li+ ions when comparing transmission-type FTIR and ATR-FTIR spectroscopy results, which complicate the determination of solvation and desolvation energies. In this study, we investigate ATR-FTIR correction algorithms, elucidate the reasons for the discrepancies between ATR-FTIR and transmission FTIR spectroscopy results, and develop a method to correct the ATR-FTIR spectrum to accurately determine the solvation structures around Li+ ions.

Solvation free energy in governing equations for DNA hybridization, protein-ligand binding, and protein folding.

This work examines the thermodynamics of model biomolecular interactions using a governing equation that accounts for the participation of bulk water in the equilibria. In the first example, the binding affinities of two DNA duplexes, one of nine and one of 10 base pairs in length, are measured and characterized by isothermal titration calorimetry (ITC) as a function of concentration. The results indicate that the change in solvation free energy that accompanies duplex formation (ΔGS) is large and unfavorable. The duplex with the larger number of G:C pairings yields the largest change in solvation free energy, ΔGS = +460 kcal·mol-1per base pair at 25 °C. A van't Hoff analysis of the data is complicated by the varying degree of intramolecular base stacking within each DNA strand as a function of temperature. A modeling study reveals how the solvation free energy alters the output of a typical ITC experiment and leads to a good, though misleading, fit to the classical equilibrium equation. The same thermodynamic framework is applied to a model protein-ligand interaction, the binding of ribonuclease A with the nucleotide inhibitor 3'-UMP, and to a conformational equilibrium, the change in tertiary structure of α-lactalbumin in molar guanidinium chloride solutions. The ribonuclease study yields a value of ΔGS = +160 kcal·mol-1, whereas the folding equilibrium yields ΔGS ≈ 0, an apparent characteristic of hydrophobic interactions. These examples provide insight on the role of solvation energy in binding equilibria and suggest a pivot in the fundamental application of thermodynamics to solution chemistry.

Preparation and evaluation of a pyridine sulfonate betaine-based zwitterionic stationary phase for hydrophilic interaction chromatography

A zwitterionic stationary phase comprising pyridinium cations and sulfonate anions was successfully developed through thiol-ene click chemistry. Using seven polar small molecules as probes, the zwitterionic stationary phase showed high separation selectivity and excellent column efficiency (35,200–54,800 plates/m) compared with two commercial columns. The influence of water proportion, salt concentration, and pH in the mobile phase, and column temperature, on the retention of six polar compounds was examined. The retention mechanism was explored by three hydrophilic retention models, Tanaka test and linear solvation energy relationship analysis. For the analysis of sample dairy products (milk powder, milk, and yogurt), the stationary phase was operated in hydrophilic interaction chromatography mode without the addition of buffer salts, facilitating rapid and efficient detection and quantification of melamine. The LOD and LOQ are 0.04 mg⋅g−1 and 0.13 mg⋅g−1, respectively, and the recovery rate is 90.3 − 102.8 %. The zwitterionic stationary phase has the advantages of simple preparation, good method reproducibility, good selectivity and high precision.

Effect of the alkyl chain length of α,ω-dichloroalkane on the Gibbs energy of transfer for functional groups.

Ion-transfer reactions of alkyl and perfluoroalkyl carboxylate ions (CH3(CH2)n-2COO- with n = 8-12 and CF3(CF2)n-2COO- with n = 3-9) were investigated at the polarized Cl-(CH2)m-Cl with m = 2, 4, 6, and 8 (O) | water (W) interface to evaluate the effects of n and m on the solvation energy of the ions, as well as on their methylene and terminal groups. These ions exhibited reversible or quasi-reversible voltammetric waves due to their transfer across the O | W interfaces, enabling the determination of formal potentials and formal Gibbs transfer energies from O to W, . The values for CH3(CH2)n-2COO- and CF3(CF2)n-2COO- increased linearly with n, allowing the estimation of for methylene and difluoromethylene groups, (CH2) and (CF2), and for their terminal groups, (COO- + CH3) and (COO- + CF3). Whereas the (CH2) and (CF2) hardly changed with the variation in m, the (COO- + CH3) and (COO- + CF3) decreased noticeably. These results suggest that the solvation energy for ions in Cl-(CH2)m-Cl increases with m, regardless of hydrophilic or lipophilic nature of the ions. Based on these findings, the advantage of using Cl-(CH2)m-Cl with a large m as a non-aqueous solvent for ion-transfer voltammetry was discussed.

Exploring the molecular interactions of dorzolamide hydrochloride in aqueous solutions: a comprehensive study using thermophysical techniques and molecular simulations across various temperatures

This research begins with the initial assessment of key physical properties—density ( ρ), speed of sound ( u), and refractive index ( nD ) in binary aqueous solutions of dorzolamide hydrochloride. These properties were measured in solutions with varying molalities, ranging from 0.0 to 0.0276, under standard atmospheric conditions and at temperatures between 300.15 and 320.15 K. From these measurements, we were able to calculate a range of important parameters, including the apparent molar volume ( Vϕ), its limiting values ([Formula: see text]), limiting apparent molar expansivity ([Formula: see text]), Helper’s constant ([Formula: see text]), apparent molar isentropic compression ( κϕ), and its limiting values ([Formula: see text]), specific refraction ( RD), molar refraction ( RM), and hydration number ( nh). These calculated values provided insights into the strong interactions between the solute and solvent in the liquid system. They also shed light on the presence of various molecular interactions, such as inter- and intramolecular hydrogen bonds, dipole–dipole, and dipole-induced-dipole interactions, as well as the ability of the solute to influence the structure of the surrounding water molecules. Additionally, we derived the densities of dorzolamide hydrochloride solutions from molecular dynamics (MD) simulations, which showed a strong correlation with our experimental findings, underscoring the reliability of MD simulations in such studies. In the end, our combined study on dorzolamide hydrochloride in water demonstrates that it shapes structures because it has lower volumetric properties and stronger hydrogen bonds in MD simulations. More alchemical free energy calculations show that there are big changes in solvation energies. This is important for understanding how dorzolamide hydrochloride works in water.

New approach: Solvent effects of benzaldehyde in aliphatic alcohol solvents with FT-IR spectroscopy and augmented vertex-adjacency matrices

In this study of benzaldehyde (BNZ) in 7 aliphatic alcohol solvents, were undertaken to investigate the solvent/solute interaction with Fouirer Transform Infrared Spectroscopy (FT-IR), BNZ produced carbonyl group stretching vibration frequencies (ν(C=O)1 and ν(C=O)2). In the present paper, ν(C=O)1 and ν(C=O)2 of BNZ are observed and correlated with the solvent acceptor number (AN), dipolarity/polarizability (π*), hydrogen-bond donor acidity (α), hydrogen-bond acceptor basicity (β), solvent polarity parameter [ET(30)], KBM parameter ƒ(ε), boiling point (bp) and the linear solvation energy relationship (LSER) using single and multiple regression analysis. All calculations were performed for ν(C=O)1 and ν(C=O)2 and compared. ν(C=O)1 was believed to correspond to a single hydrogen bonded BNZ-alcohol complex, while ν(C=O)2 was believed to correspond to a multiple hydrogen bonded BNZ-alcohol complex and multiple hydrogen-bonded alcohol molecules. It was seen that LSER equation as multiple regression analysis was better than single regression analysis in explaining solute-solvent interactions due to the higher correlation coefficients obtained. The correlations of a set of the mathematical characteristics of the augmented vertex-adjacency matrices (AV-AM) obtained by the chemical structure of BNZ with solvent parameters and ν(C=O) were investigated. It was determined that the correlation factors (R2 > 0.999) of the multiparameter regression models obtained using the matrix results were higher than the correlation factors (R2 > 0.99) of the LSER equations including solvent parameters.