Search NASA⌕ Search

SEARCH · Search NASA

Results for “Debye Model”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Synthesis of cerium oxide (CeO 2 ) by co-precipitation for application as a reference material for X-ray powder diffraction peak widths

Several methods can be used to obtain, from powder diffraction patterns, crystallite size and lattice strain of polycrystalline samples. Some examples are the Scherrer equation, Williamson–Hall plots, Warren/Averbach Fourier decomposition, Whole Powder Pattern Modeling, and Debye function analysis. To apply some of these methods, it is necessary to remove the contribution of the instrument to the widths of the diffraction peaks. Nowadays, one of the main samples used for this purpose is the LaB6 SRM660b commercialized by the National Institute of Standard Technology; the width of the diffraction peak of this sample is caused only by the instrumental apparatus. However, this sample can be expensive for researchers in developing countries. In this work, the authors present a simple route to obtain micron-sized polycrystalline CeO 2 that have a full width at half maximum comparable with the SRM660b and therefore it can be used to remove instrumental broadening.

de Lima Batista, Anderson Márcio↗
DOE OSTI2021

Thermal Conductivity of β-Phase Ga 2 O 3 and (Al x Ga 1– x ) 2 O 3 Heteroepitaxial Thin Films

Heteroepitaxy of β-phase gallium oxide (β-Ga 2 O 3 ) thin films on foreign substrates shows promise for the development of next-generation deep ultraviolet solar blind photodetectors and power electronic devices. In this work, the influences of the film thickness and crystallinity on the thermal conductivity of ($\bar{2}01$)-oriented β-Ga 2 O 3 heteroepitaxial thin films were investigated. Unintentionally doped β-Ga 2 O 3 thin films were grown on c-plane sapphire substrates with off-axis angles of 0° and 6° toward $\langle$$11\bar{2}0$$\rangle$ via metal–organic vapor phase epitaxy (MOVPE) and low-pressure chemical vapor deposition. The surface morphology and crystal quality of the β-Ga 2 O 3 thin films were characterized using scanning electron microscopy, X-ray diffraction, and Raman spectroscopy. The thermal conductivities of the β-Ga 2 O 3 films were measured via time-domain thermoreflectance. The interface quality was studied using scanning transmission electron microscopy. The measured thermal conductivities of the submicron-thick β-Ga 2 O 3 thin films were relatively low as compared to the intrinsic bulk value. The measured thin film thermal conductivities were compared with the Debye–Callaway model incorporating phononic parameters derived from first-principles calculations. The comparison suggests that the reduction in the thin film thermal conductivity can be partially attributed to the enhanced phonon-boundary scattering when the film thickness decreases. They were found to be a strong function of not only the layer thickness but also the film quality, resulting from growth on substrates with different offcut angles. Growth of β-Ga 2 O 3 films on 6° offcut sapphire substrates was found to result in higher crystallinity and thermal conductivity than films grown on on-axis c-plane sapphire. However, the β-Ga 2 O 3 films grown on 6° offcut sapphire exhibit a lower thermal boundary conductance at the β-Ga 2 O 3 /sapphire heterointerface. In addition, the thermal conductivity of MOVPE-grown ($\bar{2}01$)-oriented β-(Al x Ga 1–x ) 2 O 3 thin films with Al compositions ranging from 2% to 43% was characterized. Because of phonon-alloy disorder scattering, the β-(Al x Ga 1–x ) 2 O 3 films exhibit lower thermal conductivities (2.8–4.7 W/m∙K) than the β-Ga 2 O 3 thin films. The dominance of the alloy disorder scattering in β-(Al x Ga 1–x ) 2 O 3 is further evidenced by the weak temperature dependence of the thermal conductivity. This work provides fundamental insight into the physical interactions that govern phonon transport within heteroepitaxially grown β-phase Ga 2 O 3 and (Al x Ga 1–x ) 2 O 3 thin films and lays the groundwork for the thermal modeling and design of β-Ga 2 O 3 electronic and optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A study of electroactive polyvinyl chloride (PVC) gel actuators through the use of the electric modulus formalism and cyclic linear voltage sweeps

The conductivity and dielectric properties are integral to the function of polyvinyl chloride (PVC) gel actuators. The frequency-dependent properties of PVC gel actuators are investigated here in terms of their impedance, permittivity, and for the first time, electric modulus. The data shows that PVC gels’ conductive properties are just as, if not more, important as their dielectric properties in electromechanical transduction applications. The electrode polarization (EP) and its impact on the impedance and dielectric spectra of PVC gels, as well as the developed asymmetric space charge at the anode, are discussed. The electric modulus and tan δ spectra are used for the fitting of Cole-Cole (CC) and Debye relaxation models for gels of varying plasticizer content. The electrostatic adhesive force for PVC gels of varying plasticizer is also measured, indicating large electrostatic adhesion (>2 N cm -2 ). A cyclic linear voltage sweep is used to clarify the dynamics of space charge within the gels. The peak current (associated with space charge development) is seen to be concurrent with the onset of mechanical deformation, showing the asymmetric charge as the origin of electromechanical transduction. Additionally, the maximum charge transferred (as measured by the integration of current over time) before space charge development is found to correlate with the electrostatic adhesive force measured for the gels, pointing to a new method of characterizing PVC gels for actuation.

Instruments & Instrumentation↗
DOE OSTI2024

Computational design of high-entropy rare earth aluminum garnets for advanced thermal and environmental barrier coatings

To enhance the protection of Ni-based superalloys in gas turbine engines’ high-temperature environments, it’s crucial to develop advanced thermal/environmental barrier coating (T/EBC) materials with a balanced combination of thermal and mechanical properties. This optimization is essential to safeguard against chemical and thermal challenges. Here, in this study, we harness the power of density functional theory (DFT) in conjunction with combinatorial chemistry methodologies to engineer high-performance high-entropy rare earth disilicates of the RE 3 Al 5 O 12 family (where RE denotes Y, Gd, Er, and Yb). These materials are meticulously designed to exhibit superior phase stability, a targeted coefficient of thermal expansion (CTE), low lattice thermal conductivity, and robust mechanical properties. The determination of CTE values is accomplished through phonon calculations at various volume settings within the quasi-harmonic approximation, while lattice thermal conductivities are rigorously assessed employing the Debye-Callaway model, accounting for three distinct phonon processes. Our findings highlight the remarkable attributes of the solid solution (Y 1/4 Gd 1/4 Er 1/4 Yb 1/4 ) 3 Al 5 O 12 , which displays a reduction in lattice thermal conductivity compared to its individual constituents while maintaining a favorable range of CTE values. The novel T/EBC material, distinguished by their multifaceted functionalities, are poised to use in substantial enhancements in the performance of engines.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles

It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Ground-State and Thermodynamical Properties of Uranium Mononitride from Anharmonic First-Principles Theory

We report on an advanced density-functional theory (DFT) approach for investigating the ground-state and thermodynamical properties of uranium mononitride (UN). The electronic structure for UN at zero temperature is obtained from DFT that utilizes the generalized gradient approximation (GGA) for the electron exchange and correlation functional and includes spin-orbit interaction and an extension with orbital polarization. Thermodynamical properties are computed within the quasi-harmonic approximation in the Debye–Grüneisen model while anharmonicity is captured in the self-consistent ab initio lattice dynamics (SCAILD) scheme. Anharmonic phonons have heretofore never been modeled from first-principles for UN but they turn out to be important. The computed free energy compares well with that of a CALPHAD (CALculation of PHAse Diagrams) assessment of available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Pressure effects on the electronic structure and superconductivity of ( TaNb ) 0.67 ( HfZrTi ) 0.33 high entropy alloy

Effects of pressure on the electronic structure, electron-phonon interaction and superconductivity of the high entropy alloy (TaNb) 0.67 (HfZrTi) 0.33 , is studied in the range from 0 to 100 GPa. Electronic structure is calculated using the Korringa-Kohn-Rostoker method and the coherent potential approximation. Effects of pressure on the lattice dynamics were simulated using the Debye- Gru ¨neisen model, and the Gruneisen parameter at ambient conditions, Debye temperature and Sommerfeld electronic heat capacity coefficient were additionally experimentally determined. The electron-phonon coupling parameter λ is calculated using the McMillan-Hopfield parameters, computed in the rigid muffin tin approximation. We find, that the system undergoes the Lifshitz transition, as one of the bands cross the Fermi level at elevated pressures. The electron-phonon coupling parameter ? decreases above 10 GPa, however, calculated superconducting Tc increases up to 40 - 50 GPa and later is stabilized at the larger value than for the ambient conditions, in agreement with the experimental findings. Our results show, that the evolution of Tc with pressure in (TaNb) 0.67 (HfZrTi) 0.33 can be very well explained by the classical electron-phonon mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

The structure of boron in boron fibres

The structure of noncrystalline, chemically vapour-deposited boron fibres was investigated by computer modelling the experimentally obtained X-ray diffraction patterns. The diffraction patterns from the models were computed using the Debye scattering equation. The modelling was done utilizing the minimum nearest-neighbour distance, the density of the model, and the broadening and relative intensity of the various peaks as boundary conditions. The results suggest that the fibres consist of a continuous network of randomly oriented regions of local atomic order, about 2 nm in diameter, containing boron atoms arranged in icosahedra. Approximately half of these regions have a tetragonal structure and the remaining half a distorted rhombohedral structure. The model also indicates the presence of many partial icosahedra and loose atoms not associated with any icosahedra. The partial icosahedra and loose atoms indicated in the present model are in agreement with the relaxing sub-units which have been suggested to explain the anelastic behavior of fibre boron and the loosely bound boron atoms which have been postulated to explain the strengthening mechanism in boron fibres during thermal treatment.

Bhardwaj, J.↗
NASA NTRS1989

LEO high voltage solar array arcing response model, continuation 5

The modeling of the Debye Approximation electron sheaths in the edge and strip geometries was completed. Electrostatic potentials in these sheaths were compared to NASCAP/LEO solutions for similar geometries. Velocity fields, charge densities and particle fluxes to the biased surfaces were calculated for all cases. The major conclusion to be drawn from the comparisons of our Debye Approximation calculations with NASCAP-LEO output is that, where comparable biased structures can be defined and sufficient resolution obtained, these results are in general agreement. Numerical models for the Child-Langmuir, high-voltage electron sheaths in the edge and strip geometries were constructed. Electrostatic potentials were calculated for several cases in each of both geometries. Velocity fields and particle fluxes were calculated. The self-consistent solution process was carried through one cycle and output electrostatic potentials compared to NASCAP-type input potentials.

Metz, Roger N.↗
DOE OSTI2019

Thermal and Elastic Properties of Ce{sub x}Th{sub 1 –x}O{sub 2} Mixed Oxides: A Self-Consistent Thermodynamic Approach

A self-consistent thermodynamic model of thorium dioxide and its dilute alloy Ce{sub 0.05}Th{sub 0.95}O{sub 2}, which are of great importance for atomic power engineering, has been constructed based on generalized Debye and Einstein models. In terms of this model, the temperature dependences of the bulk modulus, density, and heat capacity of the systems have been obtained. The anharmonicity of the acoustical and optical modes of the crystal lattice and its change in passing from ThO{sub 2} to Ce{sub 0.05}Th{sub 0.95}O{sub 2} has been discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Temperature-dependent Structural and Spectroscopic Studies of (Bi 1– x Fe x )FeO 3

We report on temperature-dependent structural and spectroscopic features of (Bi 1–x Fe x )FeO 3 perovskite for x = 0.15 and 0.25. Samples were synthesized by heating quantum crystalline precursors obtained by the polyol method. Crystal structures of each composition were obtained from in-house X-ray, synchrotron X-ray, and time-of-flight neutron powder diffraction data Rietveld refinements. Partial replacement of the Bi site by the Fe 3+ cation significantly changes the crystal physicochemical properties, such as thermal expansion, polyhedral distortion, Debye temperature, and vibrational and magnetic properties. Whereas BiFeO 3 is multiferroic, both Bi 0.85 Fe 0.15 FeO 3 and Bi 0.75 Fe 0.25 FeO 3 are found to be superparamagnetic, as observed by temperature-dependent Mössbauer and SQUID measurements. Lattice thermal expansion was modeled using the Debye–Einstein-anharmonicity approach. Debye temperatures obtained from the mean-squared atomic displacement parameter and lattice thermal expansion are compared. Furthermore, temperature dependence of selective Raman modes is also analyzed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Black carbon aerosol number and mass concentration measurements by picosecond short-range elastic backscatter lidar

Abstract Black carbon aerosol emissions are recognized as contributors to global warming and air pollution. There remains, however, a lack of techniques to remotely measure black carbon aerosol particles with high range and time resolution. This article presents a direct and contact-free remote technique to estimate the black carbon aerosol number and mass concentration at a few meters from the emission source. This is done using the Colibri instrument based on a novel technique, referred to here as Picosecond Short-Range Elastic Backscatter Lidar (PSR-EBL). To address the complexity of retrieving lidar products at short measurement ranges, we apply a forward inversion method featuring radiometric lidar calibration. Our method is based on an extension of a well-established light-scattering model, the Rayleigh–Debye–Gans for Fractal-Aggregates (RDG-FA) theory, which computes an analytical expression of lidar parameters. These parameters are the backscattering cross-sections and the lidar ratio for black carbon fractal aggregates. Using a small-scale Jet A-1 kerosene pool fire, we demonstrate the ability of the technique to quantify the aerosol number and mass concentration with centimetre range-resolution and millisecond time-resolution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Local structure, bonding, and asymmetry of $\mathrm{((NH_2)_2CH)PbBr_3}$, $\mathrm{CsPbBr_3}$, and $\mathrm{(CH_3MH_3)PbBr_3}$

Here we report local structure measurements for CsPbBr 3 and ((NH 2 ) 2 CH)PbBr 3 (FAPbBr 3 ) and compare them with recent results for (CH 3 NH 3 )PbBr 3 (MAPbBr 3 ). The Pb-Br bonding is similar for all three systems; the effective spring constants, κ, are comparable (ranging from 1.20 to 1.95 eV/Å 2 ), but small in magnitude indicating very soft materials. However, there are also important differences between the three systems. Static disorder is very small for CsPbBr3 but increases somewhat with the size of the organic molecular ions MA + and FA + . At room temperature, dynamic disorder dominates in all compounds. The thermal disorder of the Pb-Br pair distribution function (PDF), i.e., the Debye-Waller factor σ 2 follows a correlated Debye or Einstein model up to 300 K in CsPbBr 3 (orthorhombic phase), but for FAPbBr 3 and MAPbBr 3 , there is a break in the σ 2 (T) curve at the orthorhombic-tetragonal transition (o-t) near 150 K, indicating a small change in the spring constant κ. κ increases for MAPbBr 3 but decreases in FAPbBr 3 at this transition. These changes are attributed to changes in the H-bonding between Br - and MA + or FA + at this transition, as a result of librations or rotations of these molecular cations. In addition, the Pb-Br PDF becomes asymmetric at a relatively low temperature for FAPbBr 3 and MAPbBr 3 , while this effect is significantly smaller for CsPbBr 3 . Finally, we address the question of a model to explain the asymmetric PDF. Two main models are discussed in the literature, an anharmonic pair potential and a split-pair distribution, possibly driven by the presence of a lone pair on the Pb ion. We show that the fourth cumulant C 4 can differentiate between these two models and other possible models. Experimentally C 4 is positive at 250 K and above, for all three systems and that is inconsistent with a split-peak model, for which C 4 is negative for splittings larger than 0.12 Å.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2012

Ocean Surface Emissivity at L-band (1.4 GHz): The Dependence on Salinity and Roughness

A characterization of the emissivity of sea water at L-band is important for the remote sensing of sea surface salinity. Measurements of salinity are currently being made in the radio astronomy band at 1.413 GHz by ESA's Soil Moisture and Ocean Salinity (SMOS) mission and NASA's Aquarius instrument aboard the Aquarius/SAC-D observatory. The goal of both missions is accuracy on the order of 0.1 psu. This requires accurate knowledge of the dielectric constant of sea water as a function of salinity and temperature and also the effect of waves (roughness). The former determines the emissivity of an ideal (i.e. flat) surface and the later is the major source of error from predictions based on a flat surface. These two aspects of the problem of characterizing the emissivity are being addressed in the context of the Aquarius mission. First, laboratory measurements are being made of the dielectric constant of sea water. This is being done at the George Washington University using a resonant cavity. In this technique, sea water of known salinity and temperature is fed into the cavity along its axis through a narrow tube. The sea water changes the resonant frequency and Q of the cavity which, if the sample is small enough, can be related to the dielectric constant of the sample. An extensive set of measurements have been conducted at 1.413 GHz to develop a model for the real and imaginary part of the dielectric constant as a function of salinity and temperature. The results are compared to the predictions of models based on parameterization of the Debye resonance of the water molecule. The models and measurements are close; however, the differences are significant for remote sensing of salinity. This is especially true at low temperatures where the sensitivity to salinity is lowest.

LeVine, D. M.↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Electrostatic shape control of a charged molecular membrane from ribbon to scroll

Bilayers of amphiphiles can organize into spherical vesicles, nanotubes, planar, undulating, and helical nanoribbons, and scroll-like cochleates. These bilayer-related architectures interconvert under suitable conditions. Here, a charged, chiral amphiphile (palmitoyl-lysine, C 16 -K 1 ) is used to elucidate the pathway for planar nanoribbon to cochleate transition induced by salt (NaCl) concentration. In situ small- and wide-angle X-ray scattering (SAXS/WAXS), atomic force and cryogenic transmission electron microscopies (AFM and cryo-TEM) tracked these transformations over angstrom to micrometer length scales. AFM reveals that the large length (L) to width (W) ratio nanoribbons (L/W > 10) convert to sheets (L/W → 1) before rolling into cochleates. Furthermore, a theoretical model based on electrostatic and surface energies shows that the nanoribbons convert to sheets via a first-order transition, at a critical Debye length, with 2 shallow minima of the order of thermal energy at L/W >> 1 and at L/W = 1. SAXS shows that interbilayer spacing ( D ) in the cochleates scales linearly with the Debye length, and ranges from 13 to 35 nm for NaCl concentrations from 100 to 5 mM. Theoretical arguments that include electrostatic and elastic energies explain the membrane rolling and the bilayer separation–Debye length relationship. These models suggest that the salt-induced ribbon to cochleate transition should be common to all charged bilayers possessing an intrinsic curvature, which in the present case originates from molecular chirality. Our studies show how electrostatic interactions can be tuned to attain and control cochleate structures, which have potential for encapsulating, and releasing macromolecules in a size-selective manner.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Characteristic temperatures of the Mossbauer fraction and thermal-shift measurements in iron and iron salts.

Measurement of Mossbauer spectra in metallic iron, sodium nitroprusside, sodium ferrocyanide, and potassium ferrocyanide absorbers between 78 and 293 K. The temperature dependences of the Mossbauer fraction and the resonant velocity were fitted to Einstein and Debye lattice-vibration models. The characteristic temperatures of the models fitted to the Mossbauer fraction are consistently lower than those fitted to the resonant velocity showing the sensitivity of the Mossbauer fraction to low-frequency modes of vibration. The characteristic temperatures obtained from the resonant velocity are higher for the salts than for the metal, indicating the presence of higher-frequency modes in the salts. This interpretation is verified semiquantitatively by comparing the thermal-shift Debye temperatures of the salts to their infrared absorption frequencies. The Mossbauer fraction of potassium ferrocyanide shows a weaker temperature dependence than that expected for a harmonic solid, which suggests that potassium ferrocyanide is anharmonic in the temperature range studied.

Lafleur, L. D.↗
DOE OSTI2017

A molecular Debye-Huckel theory of solvation in polar fluids: An extension of the Born model

A dielectric response theory of solvation beyond the conventional Born model for polar fluids is presented. The dielectric response of a polar fluid is described by a Born response mode and a linear combination of Debye-Hückel-like response modes that capture the nonlocal response of polar fluids. The Born mode is characterized by a bulk dielectric constant, while a Debye-Hückel mode is characterized by its corresponding Debye screening length. Both the bulk dielectric constant and the Debye screening lengths are determined from the bulk dielectric function of the polar fluid. The linear combination coefficients of the response modes are evaluated in a self-consistent way and can be used to evaluate the electrostatic contribution to the thermodynamic properties of a polar fluid. In conclusion, our theory is applied to a dipolar hard sphere fluid as well as interaction site models of polar fluids such as water, where the electrostatic contribution to their thermodynamic properties can be obtained accurately.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗