SEARCH · Search NASA

Results for “COHERENCE COEFFICIENT”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Linear polarization–direction correlations in $\gamma $-ray scattering experiments

Scattering measurements with incident linearly polarized γ rays provide information on spins, parities, and γ-ray multipolarity mixing coefficients, and, therefore, on the nuclear matrix elements involved in the transitions. Here, we present the general formalism for analyzing the observed angular correlations. The expressions are used to compute three-dimensional radiation patterns, which are important tools for optimizing experimental setups. Frequently, γ-ray transitions can proceed via two multipolarities that mix coherently. In such cases, the relative phases of the nuclear matrix elements are important when comparing results from different measurements. We discuss different phase conventions that have been used in the literature and present their relationships. Finally, we propose a basic experimental geometry consisting of detectors located at four different spatial locations. For this geometry, we present the measured anisotropies of the emitted γ rays in graphical format as an aid in the data analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Extremely Slow Diffusion of Gold Nanoparticles under Confinement in Mesoporous Silica

Here, we report an X-ray photon correlation spectroscopy (XPCS) study of the mobility of colloidal gold nanoparticles with a diameter of approximately 4 nm in glycerol under confinement within the channels of SBA-15 mesoporous silica with pore diameters of approximately 6 and 12 nm. The XPCS correlation functions result from an effective heterodyne signal due to the mixing of the coherent scattering from the nanoparticles and the mesoporous silica. Over the range of wave vectors and hence length scales probed, the nanoparticle dynamics are well described by one-dimensional diffusion. The nanoparticle diffusivity varies with temperature in a manner expected based on the temperature-dependent viscosity of glycerol; however, the magnitudes of the diffusion coefficients are several orders of magnitude smaller than those of the nanoparticles in bulk glycerol. We consider mechanisms that might contribute to this reduction in diffusivity including enhancement in hydrodynamic drag under confinement, effects of nanoparticle adsorption to the pore walls, and slowing structural dynamics in glycerol due to the combination of the nanometer-scale confinement and the presence of the nanoparticles.

36 MATERIALS SCIENCE↗

NASA NTRS1996

Structure Identification Within a Transitioning Swept-Wing Boundary Layer

Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.

Chapman, Keith↗

DOE OSTI2023

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1985

Advances in sputtered and ion plated solid film lubrication

The glow discharge or ion assisted vacuum deposition techniques, primarily sputtering and ion plating, have rapidly emerged and offer great potential to deposit solid lubricants. The increased energizing of these deposition processes lead to improved adherence and coherence, favorable morphological growth, higher density, and reduced residual stresses in the film. These techniques are of invaluable importance where high precision machines tribo-components require very thin, uniform lubricating films (0.2 m), which do not interface with component tolerances. The performance of sputtered MoS2 films and ion plated Au and Pb films are described in terms of film thickness, coefficient of friction, and wear lives.

Spalvins, T.↗

DOE OSTI2024

Dielectric loss due to charged-defect acoustic phonon emission

The coherence times of state-of-the-art superconducting qubits are limited by bulk dielectric loss, yet the microscopic mechanism leading to this loss is unclear. Here, we propose that the experimentally observed loss can be attributed to the presence of charged defects that enable the absorption of electromagnetic radiation by the emission of acoustic phonons. Our explicit derivation of the absorption coefficient for this mechanism allows us to derive a loss tangent of 7.2 × 10−9 for Al2O3, in good agreement with recent high-precision measurements [Read et al., Phys. Rev. Appl. 19, 034064 (2023)]. We also find that for temperatures well below ∼0.2 K, the loss should be independent of temperature, which is also in agreement with observations. Our investigations show that the loss per defect depends mainly on properties of the host material, and a high-throughput search suggests that diamond, cubic BN, AlN, and SiC are optimal in this respect.

Turiansky, Mark E. (ORCID:0000000291543582)↗

NASA NTRSDate not supplied

Development of the Ames Global Hyperspectral Synthetic Dataset

This study develops the surface BRDF (bidirectional reflectance distribution function) product of the Ames Global Hyperspectral Synthetic Dataset (AGHSD), based on the corresponding MODIS products, to support the NASA Surface Biology and Geology mission development. A main challenge in deriving a hyperspectral dataset from the multi-band satellite products is how to identify a succinct yet robust algorithm that allow us to infer BRDF at unobserved wavelengths based on the few observed bands. Using the theories of radiative transfer in vegetation canopies, we arrive at a simple equation that accurately approximates hyperspectral surface BRDF as the weighted sum of components from the soil and the vegetation. Each of the components is modeled by the product of the spectrally-dependent optical properties of a surface element (the spectra of the soil surface reflectance, the leaf single albedo, or the canopy scattering coefficient) and a spectrally-independent bidirectional scattering function. The optical properties of the soil and the vegetation can be obtained from existing spectral libraries or model simulations. The bidirectional scattering functions are represented by the Ross-Thick-Li-Sparse BRDF model, where the linear coefficients are estimated with regression analysis from the multi-band MODIS data. We validate the algorithm with simulations by Monte Carlo Ray Tracing model experiments, and the results are highly consistent with the theoretic derivation. We apply the algorithm to generate the AGHSD BRDF product at 1km and 8-day resolutions for the year of 2019. The results are biogeochemically and physically coherent and consistent, and thus serve the goal to support the science and application development of the SBG community.

Hyperspectral↗

DOE OSTI2024

Cyclic superconducting refrigerators using guided fluxon propagation

We propose cyclic refrigeration in solid state, employing a gas of magnetic field vortices in a type-II superconductor—also known as fluxons—as the cooling agent. Refrigeration cycles are realized by envisioning a racetrack geometry consisting of both adiabatic and isothermal arms, etched into a type-II superconductor. The guided propagation of fluxons in the racetrack is achieved by applying an external electrical current, in a Corbino geometry, through the sample. A gradient of magnetic field is set across the racetrack allowing one to adiabatically cool down and heat up the fluxons, which subsequently exchange heat with the cold and hot reservoirs, respectively. We characterize the steady state of refrigeration cycles thermodynamically for both s -wave and d -wave pairing symmetries, and present their figures of merit such as the cooling power delivered, and the coefficient of performance. Our cooling principle can offer significant cooling for on-chip microrefrigeration purposes, by locally cooling below the base temperatures achievable in a conventional dilution refrigerator. We estimate n W / m m 2 of cooling power per unit area assuming a tunnel coupling with ∼ M Ω µ m 2 specific resistance. Integrating the fluxon fridge to quantum circuits can enhance their coherence time by locally suppressing thermal fluctuations, and improve the efficiency of single photon detectors and charge sensors. Published by the American Physical Society 2024

Karmakar, Tathagata (ORCID:0000000187812017)↗

NASA NTRS1992

The NDSC Primary Site at the Observatoire De Haute-provence: Ozone and Aerosols Observations (1986-1992)

The primary site of the newly implemented Network for Detection of Stratospheric Changes (NDSC) includes three main sites in the alpine regions of southern Europe. At one of these observatories, a coherent set of experimental systems has been established to monitor on a long-term basis ozone, aerosols, minor constituents, and temperature vertical distributions in the troposphere and stratosphere. Such techniques include several lidar systems, balloon-borne sondes and ultraviolet visible spectrometry. Lidar measurements of stratospheric aerosols are performed on a routine basis, since 1980. These include the vertical distribution of the extinction and backscattering coefficients from the tropopause. Recently, measurements performed following the eruption of Mount Pinatubo, have allowed researchers to follow the time evolution of the aerosol cloud in middle latitude regions.

Megie, G.↗

NASA NTRS1996

Macroscopic Violation of Three Cauchy-Schwarz Inequalities Using Correlated Light Beams From an Infra-Red Emitting Semiconductor Diode Array

We briefly review quantum mechanical and semi-classical descriptions of experiments which demonstrate the macroscopic violation of the three Cauchy-Schwarz inequalities: g(sup 2)(sub 11)(0) greater than or equal to 1; g(sup 2)(sub 11)(0) greater than or equal to g(sup 2)(sub 11)(t), (t approaches infinity); (the absolute value of g(sup 2)(sub 11)(0))(exp 2) less than or equal to g(sup 2)(sub 11)(0) g(sup 2)(sub 11)(0). Our measurements demonstrate the violation, at macroscopic intensities, of each of these inequalities. We show that their violation, although weak, can be demonstrated through photodetector current covariance measurements on correlated sub-Poissonian Poissonian, and super Poissonian light beams. Such beams are readily generated by a tandem array of infrared-emitting semiconductor junction diodes. Our measurements utilize an electrically coupled array of one or more infrared-emitting diodes, optically coupled to a detector array. The emitting array is operated in such a way as to generate highly correlated beams of variable photon Fano Factor. Because the measurements are made on time scales long compared with the first order coherence time and with detector areas large compared with the corresponding coherence areas, first order interference effects are negligible. The first and second inequalities are violated, as expected, when a sub-Poissonian light beam is split and the intensity fluctuations of the two split beams are measured by two photodetectors and subsequently cross-correlated. The third inequality is violated by bunched (as well as anti-bunched) beams of equal intensity provided the measured cross correlation coefficient exceeds (F - 1)/F, where F is the measured Fano Factor of each beam. We also investigate the violation for the case of unequal beams.

Edwards, P. J.↗

DOE OSTI2022

The calculation of light element impurity (α, n) yield curves in a PuO 2 matrix and associated specific yield coefficients: Influence of the reaction cross sections

Most of the Pu separated from irradiated commercial nuclear fuel is stored as PuO 2 . The primary quantitative nondestructive measurement technique used to verify the amount of Pu in storage containers is passive neutron correlation counting. An important physical property of the oxide material is the ratio, α, of the rate of (α, n) neutrons produced inside the item to the rate of neutrons produced by spontaneous fission. This ratio influences the precision of the correlated counting method and affects the interpretation of the data because of how it changes both the primary total neutron production rate and the rate of induced fission events taking place inside the item. In addition to the main O(α, n) contribution, additional contributions come from α-particle interactions with light element impurities that are inevitably present. In this work, we calculate specific (α, n) yield coefficients, expressed in units of neutrons per second per gram of α-emitting nuclide per part per million by mass of the specified impurity element distributed in a pure PuO 2 matrix, for some key α-emitting actinides commonly present in reprocessed Pu ( 238–242 Pu+ 241 Am). These coefficients are directly applicable to nuclear safeguards verification work in which the α ratio is often calculated from the Pu-isotopic composition and chemical information obtained by other means. Further, they also provide a convenient up-to-date reference set against which values generated by other methods can be compared. Results are presented for impurities with atomic number from 3 to 17 inclusive, plus K and Fe. In most cases, these coefficients are not expected to change by more than 5%–10% at any time in the future. However, as new data become available, changes as large as 20% may be needed for some targets (e.g., F). The present yield calculations are limited by the general shortage of quality experimental total (α, n) reaction cross section data, which, together with unexplained variation between determinations, means that an objective and coherent evaluation is not possible. The situation is even less satisfactory for the partial differential cross section needed to calculate neutron spectra.

(𝛼, n) reactions↗

DOE OSTI2024

Underlying mechanism of structural transformation between GaSb and GaAs response to intense electronic excitation

Ion irradiation of semiconductors has emerged as a promising approach for fabricating self-organized nanosystems with high atomic precision, despite often being accompanied by undesirable phenomena. Exploring the mechanisms underlying structural transformations is crucial for assessing nanostructure array types under complex irradiation environments. By quantitatively calculating the thermodynamically driven processes and analyzing the impact of intrinsic structural parameters, distinct structural transformations in response to intense electronic excitation are systematically investigated in gallium antimonide (GaSb) and gallium arsenide (GaAs) systems. In high-energy regimes, the nanofibers layer of GaSb exhibits intriguing structural discrepancy, characterized by partial nanofibers with coherent boundaries, interspersed nanopores accompanied by antisite defects and Ga precipitates, distinguishing to a series of discontinuous latent tracks that emerged within cylindrical trajectories in GaAs. Furthermore, significant diffusion behaviors of the nanohillocks are discovered in GaAs, with higher average roughness than GaSb, driven by the gradient stress distribution influenced by the free-surface effects. The deposition energy for melting phase formation, Gibbs free energy, and Ga diffusion coefficients contribute to the distinctive structural features, evidencing relatively stable morphological configurations and higher irradiation resistance in GaAs. Consequently, special optoelectronic properties associated with structural discrepancies facilitate the design and optimization of material functionalities by irradiation technologies.

36 MATERIALS SCIENCE↗

NASA NTRS1985

Ladder and cross terms in second-order distorted Born approximation

In the strong fluctuation theory for a bounded layer of random discrete scatterers, the second moments of the fields in the second-order distorted Born approximation are obtained for copolarized and cross-polarized fields. The backscattering cross sections per unit area are calculated by including the mutual coherence of the fields due to the coincidental ray paths, and that due to the opposite ray paths, corresponding to the ladder and cross terms in the Feynman diagramatic representation. It is proved that the contributions from ladder and cross terms for the copolarized backscattering cross sections are the same, while the contributions for the cross-polarized backscattering cross sections are of the same order. The bistatic scattering coefficients in the second-order approximation for both the ladder and cross terms are also obtained. The contributions from the cross terms explain the enhancement in the backscattering direction.

Jin, Y. Q.↗

NASA NTRS1991

The dynamics of dense and dilute clusters of drops evaporating in large, coherent vortices

The behavior of evaporating clusters of drops embedded into large, coherent vortices is described using a formulation which is valid for both dense and dilute clusters. Drops and gas interact both dynamically and thermodynamically. Dynamic coupling occurs through a force on the drops due to drag resulting from a slip velocity between the two phases. The net interaction force on the gas with drops is due to a source thrust from evaporation plus drag on each drop. The drag coefficient accounts for blowing from the drop surface. Thermodynamic coupling is a result of drop heating and evaporation. Limitations due to drop proximity on heating and evaporation are taken into account. The vortical motion of the drops in the cluster results in the formation of a core region devoid of drops at the center of the vortex, and a shell region containing the drops and surrounding the inner core. Results are presented showing the dependence of the evaporation time, the final to initial volume ratio and the final to initial shell thickness ratio upon the initial air/fuel mass ratio and as a function of the initial tangential velocities, upon the initial Stokes number, initial drop radius and initial outer cluster radius. Differences in behavior between and control parameters of dense and dilute clusters are pointed out by these new results. It is found that for dense clusters the final to initial volume ratio and final to initial shell thickness scale with the initial Stokes number, a new result which must be validated experimentally.

Bellan, J.↗

NASA NTRS1989

Redistribution in astrophysically important hydrogen lines

Under typical solar chromospheric conditions for hydrogen radiators, strong collisions due to both electrons and ions are well separated in time, so that a binary collision theory for collisional redistribution is applicable. However, a simple impact approximation may not be used, but rather a unified type theory is required in which frequency dependent line shape parameters are used to describe both impact and quasi-static regions of the spectrum. In addition, correlated terms which describe absorption and emission during a collision are important, and, in fact, without correlated terms describing both transfer of excitation and emission during the same collision unphysical predictions (such as negative intensities) would be obtained. In this paper theory is specifically developed for the coupled Lyman-alpha, Lyman-beta, Hydrogen-alpha system, and equations of statistical equilibrium and absorption and emission coefficients are given. All correlated events are examined and emission during a collision is found to be important in the line wings. Stimulated emission and absorption is also included within a broadband approximation. The major approximation is to ignore lower state interaction. It is found that for Lyman-beta Raman-coupling with Hydrogen-alpha occurs and the overall scattering of radiation in the line wings is mostly coherent. In contrast, for Hydrogen-alpha, incoherent redistribution due to lower state radiative decay (which occurs even in the absence of collisions) is found to dominate the coherent scattering. Finally, in the Lyman series the dominant incoherent contribution is associated with cascade transitions and inelastic collisions between different principal quantum states.

Cooper, J.↗

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Shebalin, John V↗

DOE OSTI2020

Computational thermodynamics and its applications

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable or un- stable, when interacting with the surroundings. In this overview, fundamentals of thermodynamics are briefly reviewed through the combination of first and second laws of thermodynamics for open and nonequilibrium systems to demonstrate that the reversible equilibrium and irreversible nonequilibrium thermodynamics can be integrated to enhance the power and utilities of thermodynamics. The recent progresses in computational thermodynamics, the remaining challenges, and potential impacts in broad scientific fields are discussed here. It is shown that computational thermodynamics enables the modeling of thermodynamics of a state as a function of both external and internal variables and quantitative calculations of a broad range of properties of a multicomponent system in terms of first and second derivatives of energy, including not only equilibrium states when there are no driving forces for any internal processes and but also non-equilibrium states with driving forces for internal processes. Consequently, external constraints such as fixed strain and internal degree of freedoms such as ordering and defects can be described in a coherent framework and applied to materials design. Furthermore, two important but largely overlooked aspects in thermodynamics will be discussed, i.e. the rigorous application of statistical thermodynamics with the probability of configurations and their contributions to system properties, and the applications of second derivatives of energy with respect to either two extensive variables or two potentials or a mixture of them in terms of understanding and predicting emergent behaviors, critical phenomena, kinetic coefficients, and mechanical properties.

36 MATERIALS SCIENCE↗