SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The [PW 12 O 40 ] 3− [M + (H 2 O) 16 ] 3 is a “hydrogen bonded” ion-pair complex for M = Li, Na and K, but is a “contact” ion-pair complex for M = Rb and Cs, intermolecular charge transfer from the solvated counter cations M + (H 2 O) 16 to the anion [PW 12 O 40 ] 3− .
Development of calcium metal batteries has been historically frustrated by a lack of electrolytes capable of supporting reversible calcium electrodeposition. In this paper, we report the study of an electrolyte consisting of Ca(BH 4 ) 2 in tetrahydrofuran (THF) to gain important insight into the role of the liquid solvation environment in facilitating the reversible electrodeposition of this highly reactive, divalent metal. Through interrogation of the Ca 2+ solvation environment and comparison with Mg 2+ analogs, we show that an ability to reversibly electrodeposit metal at reasonable rates is strongly regulated by dication charge density and polarizability. Our results indicate that the greater polarizability of Ca 2+ over Mg 2+ confers greater configurational flexibility, enabling ionic cluster formation via neutral multimer intermediates. Increased concentration of the proposed electroactive species, CaBH 4 + , enables rapid and stable delivery of Ca 2+ to the electrode interface. Finally, this work helps set the stage for future progress in the development of electrolytes for calcium and other divalent metal batteries.
We report a Many Body Energy (MBE) analysis of aqueous ionic clusters containing anions and cations at the two opposite ends of the Hofmeister series, viz. the kosmotropes Ca2+, SO42- and chaotropes NH4+ and ClO4- with 9 water molecules to quantify the how these ions in altering the interaction between the water molecules in their immediate surrounding. The current results are contrasted to the ones reported earlier for water clusters as well as for alkali metal and halide ion aqueous clusters of the same size, which lie in the middle of the Hofmeister series. Through this analysis, noteworthy differences between the MBE of kosmotropes and chaotropes were identified. The MBE of kosmotropes is dominated by ion-water interactions that extends beyond the 4-body term, the point at which the MBE of pure water converges. The percentage contribution of the 2- B to the total cluster binding energy is noticeably larger. The disruption due to the dominant ion results in weak, unfavorable water-water interactions. The MBE for chaotropes, on the other hand, was found to converge more quickly as it more closely resembles that of pure water clusters. Chaotropes exhibit weaker overall binding energies and ion-water interactions with more favorable water-water interactions, somewhat recovering the pattern of the 2-4 body terms exemplified by pure water clusters. More importantly, both kosmotropic and chaotropic ions exhibit an anticorrelation between the 2-B ion-water (I-W) and water-water (W-W) interactions as well as between the 3-B (I-W-W) and (I-W) interactions. The consideration of two different structural arrangements (ion inside and outside of a water cluster) suggests that fully solvated (ion inside) chaotropes disrupt the hydrogen bonding network in a similar manner as partially solvated (ion outside) kosmotropes and offer useful insights into the modeling requirements of bulk vs. an interface. Finally, the 2-B contribution to the total Basis Set Superposition Error (BSSE) correction for the kosmotropic and chaotropic ions follows the previously reported erf profile vs. intermolecular distance. When scaled for the corresponding dimer energies and distances, a single profile fits the current results together with all previously reported ones for the pure water and halide water clusters.
In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.
The effects of an aqueous phase on ethanol decomposition for hydrogen production over a Ru(0001) catalyst surface model have been investigated from first principles. Solvent effects on the reaction mechanism and kinetic parameters have been quantified with the help of a microkinetic reactor model, density functional theory, and an implicit solvation scheme (iSMS). Our calculations indicate that in both vapor- and aqueous-phase reaction environments, the ethanol decomposition starts with acetaldehyde formation on the surface, some of which further dehydrogenates to ketenyl species (CHCO), where the C–C bond cleaves to form methylidyne (CH) and CO. In the vapor phase, adsorbed CH gets hydrogenated to methane, and CO desorbs or undergoes methanation reducing the amount of hydrogen produced. In contrast, under aqueous phase reaction conditions, the methanation is inhibited, and the water–gas shift (WGS) reaction is accelerated, leading to complete conversion of CO to CO 2 and H 2 . Furthermore, calculations indicate that the observed reaction behavior under aqueous phase reforming conditions originates primarily from the higher water chemical potential, and implicit solvent models predict only a small solvation effect.
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Here, we carry out Diffusion Monte Carlo simulations of up to five hydrogen molecules aggregated with two Stockmayer clusters that solvate a single lithium ion. The first one contains six point dipole solvent particles with parameters tuned to emulate nitromethane. The second cluster is a relative large system investigated recently [G. DiEmma, S. Kalette, and E. Curotto, Chem. Phys. Lett. 2019, 725, 80–86]. In both cases we find that the aggregated hydrogen molecules perturb significantly the ground state of the host cluster and form a distorted tetrahedral cage around the Li + ion. The fifth hydrogen molecule is absorbed by the larger Stockmayer cluster while remaining in the proximity of the solvated charge.
The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. Here, in order to gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li + , Na + and K + . The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.
Electrolytes in lithium-ion batteries comprise solvent mixtures, but analysis of ion transport is always based on treating the solvents as a single-entity. We combine electrophoretic NMR (eNMR) measurements and molecular dynamics (MD) simulations to quantify electric-field-induced transport in a concentrated solution containing LiPF 6 salt dissolved in an ethylene carbonate/ethyl methyl carbonate (EC/EMC) mixture. The selective transport of EC relative to EMC is reflected in the difference between two transference numbers, defined as the fraction of current carried by cations relative to the velocity of each solvent species. This difference arises from the preferential solvation of cations by EC and its dynamic consequences. The simulations reveal the presence of a large variety of transient solvent-containing clusters which migrate at different velocities. Rigorous averaging over different solvation environments is essential for comparing simulated and measured transference numbers. Our study emphasizes the necessity of acknowledging the presence of four species in mixed-solvent electrolytes.
Radiation induced damage of extractant molecules is a well-known phenomenon responsible for reducing efficiency and increasing the waste and cost of reprocessing used nuclear fuel (UNF). As such, understanding early-stage (pico- to nanoseconds) radiation-induced reaction mechanisms is essential for informing the design of next generation extractants with enhanced radiation robustness. Here, in this work, we utilized picosecond and nanosecond electron pulse radiolysis experiments to probe the early-stage radioactive environment experienced by the organic phase extractant N,N,N',N'-tetraoctyldiglycolamide (TODGA), proposed for separating highly radioactive trivalent minor actinides (specifically americium and curium) from the trivalent lanthanides. Using comparisons to the similar ionization potential (IP) solute p-xylene, this work determined the mechanism of reaction with the ionized diluent (i.e., n-dodecane radical cation, DD˙ + ) is hole transfer to produce TODGA˙ + . At high TODGA concentrations (>100 mM), the majority of this transfer occurs faster than 10 ps via the capture of DD˙ + holes prior to their solvation with a C 37 = 300 mM. The surviving solvated holes were captured with k = (2.38 ± 0.15) × 10 10 M -1 s -1 . Attempts at subsequent hole transfer to lower IP solutes found that only 10% of holes were transferred, indicating bond rupture of TODGA˙ + occurs within 2.6 ns at 200 mM TODGA. Possible reaction pathways for the rapid decomposition of TODGA˙ + were explored using a combination of experiments and density functional theory (DFT) calculations.
Lithium/sodium metal batteries (LMBs/SMBs) possess immense potential for various applications due to their high energy density. Nevertheless, LMBs/SMBs are highly susceptible to the detrimental effects of an unstable solid electrolyte interphase (SEI) and dendrites during practical applications, particularly pronounced in low-temperature environments. Furthermore, sluggish ion transportation further compromises the cycling stability of LMBs/SMBs at low temperatures. To achieve stable operation of LMBs/SMBs at low temperatures, researchers have made numerous efforts including electrolyte optimization aimed at creating stable SEIs and suppressing the metal dendrites under low temperature conditions. Despite the significant advancements made recently in the liquid electrolyte design, there remains considerable hurdle in electrolyte engineering for practical low-temperature, high energy density LMBs/SMBs, calling for a profound comprehension of the intricate interplay between the electrochemical reaction kinetics and electrolyte compositions. Here, this review provides a thorough overview of various strategies in optimizing liquid electrolytes covering weakly solvating electrolytes, concentration-designed electrolytes, and solvation structure-designed electrolytes, to address the challenges faced by LMBs/SMBs at low temperatures, including slow reaction kinetics and the difficulties in Li + /Na + solvation/desolvation. Furthermore, this review discusses future prospects for the advancement of this field, intending to provide valuable insights and support for subsequent research undertakings.
Rare earth tris(alkyl) complexes such as M(CH 2 SiMe 3 ) 3 (sol) n are widely used as precursors for many compounds and as homogeneous catalysts for alkene polymerization and alkane functionalization. However, the thermal instability of those most conveniently made from the commercially available lithium salt of the neosilyl anion, LiCH 2 SiMe 3 , Li(r), restricts their utility. We present a new range of synthetically useful, more kinetically stable rare earth neosilyl solvates, derived from a full kinetic study of the various possible decomposition mechanisms of 7 known and 12 new solvated rare earth neosilyl complexes M(CH 2 SiMe 3 ) 3 (sol) n M = Sc(III), Y(III), Lu(III), Sm(III), and sol = THF; TMEDA; DMPE; diglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O, G 2 ), triglyme ((CH 3 ) 2 (OCH 2 CH 2 ) 3 O, G 3 ). Surprisingly, simply using higher-denticity donors to sterically disfavor neosilyl γ-H elimination is not effective. While Sc(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) has a half-life, t 1/2 , of 258.1 h, six times longer than for Sc(r) 3 (C 4 H 8 O) 2 (t 1/2 = 43 h), Lu(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) and Y(r) 3 ((CH 3 ) 2 (OCH 2 CH 2 ) 2 O) do not show the expected, analogous increased t 1/2 . This is because new decomposition pathways appear for poorly fitting donors. Finally, kinetic studies demonstrate the impact of small, and increasing amounts of LiCl on the kinetics of the reactivity of the smaller alkyls Y(r) 3 (THF) 2 and Lu(r) 3 (THF) 2 ; molecules used in hydrocarbon chemistry and catalysis for fifty years. A new route to pure Y(r) 3 (THF) 2 , which avoids the traditional use of Li(r), is presented.
Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.
Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.
Mixing anions is emerging as a promising strategy for multivalent electrolyte design, allowing for adjustment of the solvation structure of bulk cations and enhancing the efficiency of electrochemical processes (e.g. metal deposition for batteries and catalysis). Further progress in electrolyte development requires a fundamental understanding of how tailored electrolyte speciation in mixed anion systems can modify the dynamic electrochemical interface during metal cycling. In this study, we present an anode-focused mechanistic study of exemplar Mg electrolytes containing three different secondary anions, correlating electrochemical behavior with bulk speciation and operando interfacial dynamics. Electrospray Ionization-Mass Spectrometry (ESI-MS) results reveal a general trend of forming mixed anion contact ion pairs (CIPs) across various anions, with the extent of ion pairing influenced by the association strength of the secondary anion. Operando multiharmonic electrochemical quartz crystal microbalance with dissipation (EQCM-D) reveals how these bulk species influence interfacial mass uptake, viscoelasticity, and solvent-coupled hydrodynamic behavior during deposition and stripping. The results indicate that Mg-containing ion pairs and solvated complexes shape adsorption, nucleation, and deposit growth, leading to distinct anion-dependent interphases ranging from more permeable, solvent-coupled layers to relatively compact and rigid deposits. This work establishes a quantitative link between bulk speciation and interfacial dynamics in divalent metal electrodeposition and provides mechanistic guidance for electrolyte design.
The calculation of a protein-ligand binding free energy based on molecular dynamics (MD) simulations generally relies on a thermodynamic cycle in which the ligand is alchemically inserted into the system, both in the solvated protein and free in solution. The corresponding ligand-insertion free energies are typically calculated in nanoscale computational boxes simulated under periodic boundary conditions and considering electrostatic interactions defined by a periodic lattice-sum. This is distinct from the ideal bulk situation of a system of macroscopic size simulated under non-periodic boundary conditions with Coulombic electrostatic interactions. This discrepancy results in finite-size effects, which affect primarily the charging component of the insertion free energy, are dependent on the box size, and can be large when the ligand bears a net charge, especially if the protein is charged as well. This article investigates finite-size effects on calculated charging free energies using as a test case the binding of the ligand 2-amino-5-methylthiazole (net charge +1 e) to a mutant form of yeast cytochrome c peroxidase in water. Considering different charge isoforms of the protein (net charges −5, 0, +3, or +9 e), either in the absence or the presence of neutralizing counter-ions, and sizes of the cubic computational box (edges ranging from 7.42 to 11.02 nm), the potentially large magnitude of finite-size effects on the raw charging free energies (up to 17.1 kJ mol{sup −1}) is demonstrated. Two correction schemes are then proposed to eliminate these effects, a numerical and an analytical one. Both schemes are based on a continuum-electrostatics analysis and require performing Poisson-Boltzmann (PB) calculations on the protein-ligand system. While the numerical scheme requires PB calculations under both non-periodic and periodic boundary conditions, the latter at the box size considered in the MD simulations, the analytical scheme only requires three non-periodic PB calculations for a given system, its dependence on the box size being analytical. The latter scheme also provides insight into the physical origin of the finite-size effects. These two schemes also encompass a correction for discrete solvent effects that persists even in the limit of infinite box sizes. Application of either scheme essentially eliminates the size dependence of the corrected charging free energies (maximal deviation of 1.5 kJ mol{sup −1}). Because it is simple to apply, the analytical correction scheme offers a general solution to the problem of finite-size effects in free-energy calculations involving charged solutes, as encountered in calculations concerning, e.g., protein-ligand binding, biomolecular association, residue mutation, pK{sub a} and redox potential estimation, substrate transformation, solvation, and solvent-solvent partitioning.
Hydrated Fe{sup +} ions are produced in a laser-vaporization cluster source of a triple quadrupole mass spectrometer. The Fe{sup +}(H{sub 2}O){sub n} (n = 3–8) complexes are mass-selected and probed with infrared (IR) photodissociation spectroscopy in the OH-stretch region. Density functional theory (DFT) calculations are also carried out for analyzing the experimental IR spectra and for evaluating thermodynamic quantities of low-lying isomers. Solvation through H-bonding instead of direct coordination to Fe{sup +} is observed already at n = 3, indicating the completion of the first hydration shell with two H{sub 2}O molecules. Size dependent variations in the spectra for n = 5–7 provide evidence for the second-shell completion at n = 6, where a linearly coordinated Fe{sup +}(H{sub 2}O){sub 2} subunit is solvated with four H{sub 2}O molecules. Overall spectral features for n = 3–8 agree well with those predicted for 2-coordinated structures. DFT calculations predict that such 2-coordinated structures are lowest in energy for smaller n. However, 4-coordinated isomers are predicted to be more stable for n = 7 and 8; the energy ordering is in conflict with the IR spectroscopic observation. Examination of free energy as a function of temperature suggests that the ordering of the isomers at warmer temperatures can be different from the ordering near 0 K. For n = 7 and 8, the 4-coordinated isomers should be observed at low temperatures because they are lowest in enthalpy. Meanwhile, outer-shell waters in the 2-coordinated structures are bound less rigidly; their contribution to entropy is rather large. The 2-coordinated structures become abundant at warmer temperatures, owing to the entropy effect.
The initial purpose of the study is to systematically investigate the solvation properties of different proteins in water solution by terahertz (THz) radiation absorption. Transmission measurements of protein water solutions have been performed using a vector network analyser-driven quasi-optical bench covering the WR-3 waveguide band (0.220–0.325 THz). The following proteins, ranging from low to high molecular weight, were chosen for this study: lysozyme, myoglobin, and bovine serum albumin (BSA). Absorption properties of solutions were studied at different concentrations of proteins ranging from 2 to 100 mg/ml. The concentration-dependent absorption of protein molecules was determined by treating the solution as a two-component model first; then, based on protein absorptivity, the extent of the hydration shell is estimated. Protein molecules are shown to possess a concentration-dependent absorptivity in water solutions. Absorption curves of all three proteins sharply peak towards a dilution-limit that is attributed to the enhanced flexibility of protein and amino acid side chains. An alternative approach to the determination of hydration shell thickness is thereby suggested, based on protein absorptivity. The proposed approach is independent of the absorption of the hydration shell. The derived estimate of hydration shell thickness for each protein supports previous findings that protein-water interaction dynamics extends beyond 2-3 water solvation-layers as predicted by molecular dynamics simulations and other techniques such as NMR, X-ray scattering, and neutron scattering. According to our estimations, the radius of the dynamic hydration shell is 16, 19, and 25 Å, respectively, for lysozyme, myoglobin, and BSA proteins and correlates with the dipole moment of the protein. It is also seen that THz radiation can serve as an initial estimate of the protein hydrophobicity.