Search

O and Nb are two representative alloying elements of Ti to form high-temperature and corrosion resistance α Ti alloys. The investigation on the thermodynamic characteristics of α Ti–O and Ti–Nb has attracted much attention in recent years. However, in this regard, a satisfied experimental technique or modeling scheme is still yet to be developed due to the appearance of a variety of oxides in Ti–O and the mechanical instability present in Ti–Nb. Herein, we combined first-principles calculations with the cluster expansion method to investigate the ground-state characteristics for α Ti–O and α Ti–Nb systems. The atomic bonding interactions in these two systems were first revealed based on the calculated electronic structures. Afterward, the Debye–Grüneisen model and Monte Carlo simulations were employed together to investigate the thermodynamic properties of α phases in these two systems, and the effect of vibrational entropy on the order–disorder transition temperatures of the phases in α Ti–O system was first examined. A good agreement with experimentally reported phase boundaries is obtained in the Ti–Nb system by handling the mechanical instabilities introduced by the highly distorted structures. In addition, the cluster expansion coefficients for the Ti–O and Ti–Nb system offer a good starting point to investigate the phase equilibrium in Ti–Nb–O ternary alloy. We also believe the insights provided here would be helpful for those who would like to seek an efficient scheme they are confident with to investigate the phase thermodynamic properties of other hcp Ti-based alloys. [ABSTRACT FROM AUTHOR]

Water shows anomalous properties that are enhanced upon supercooling. The unusual behavior is observed in both H2O and D2O, however, with different temperature dependences for the two isotopes. It is often noted that comparing the properties of the isotopes at two different temperatures (i.e., a temperature shift) approximately accounts for many of the observations—with a temperature shift of 7.2 K in the temperature of maximum density being the most well-known example. However, the physical justification for such a shift is unclear. Motivated by recent work demonstrating a "corresponding-states-like" rescaling for water properties in three classical water models that all exhibit a liquid–liquid transition and critical point [Uralcan et al., J. Chem. Phys. 150, 064503 (2019)], the applicability of this approach for reconciling the differences in the temperature- and pressure-dependent thermodynamic properties of H2O and D2O is investigated here. Utilizing previously published data and equations-of-state for H2O and D2O, we show that the available data and models for these isotopes are consistent with such a low temperature correspondence. These observations provide support for the hypothesis that a liquid–liquid critical point, which is predicted to occur at low temperatures and high pressures, is the origin of many of water's anomalies. [ABSTRACT FROM AUTHOR]

Light and heavy water show similar anomalies in thermodynamic and dynamic properties, with a consistent trend of anomalies occurring at higher temperatures in heavy water. Viscosity also increases faster upon cooling in heavy water, causing a giant isotope effect, with a viscosity ratio near 2.4 at 244 K. While a simple temperature shift apparently helps in collapsing experimental data for both isotopes, it lacks a clear justification, changes value with the property considered, and requires additional ad hoc scaling factors. Here, we use a corresponding states analysis based on the possible existence of a liquid–liquid critical point in supercooled water. This provides a coherent framework that leads to the collapse of thermodynamic data. The ratio between the dynamic properties of the isotopes is strongly reduced. In particular, the decoupling between viscosity η and self-diffusion D, measured as a function of temperature T by the Stokes–Einstein ratio Dη/T, is found to collapse after applying the corresponding states analysis. Our results are consistent with simulations and suggest that the various isotope effects mirror the one on the liquid–liquid transition. [ABSTRACT FROM AUTHOR]

The long-term storage performance of energetic multilayer nanofilms is of great significance for their applications. In this paper, it is proposed to add a 10 nm Cu barrier layer between Al/CuO composite films to increase their storage stability. The Al/CuO composite film and Al/Cu/CuO composite film were aged for 14 days in an environment with a relative humidity of 40% and a temperature of 71 °C. Scanning electron microscopy and differential scanning calorimetry were used to analyze the microstructure and thermodynamic properties of the energetic films before and after aging, and the electrical detonation performance and ignition ability of energy-containing semiconductor bridges were studied. The results indicate that after aging for 14 days in an environment with a relative humidity of 40% and a temperature of 71 °C, the Al layer of the Al/CuO composite film becomes thinner, the Al2O3 interface layer increases, and the heat release decreases. The interlayer microstructure of the Al/Cu/CuO energetic multilayer nanofilms did not change significantly, and the addition of a 10 nm Cu layer formed a low-temperature Al–Cu alloy, reducing the reaction initiation temperature from 626 to 570 °C. The critical ignition time and critical ignition energy of the Al/CuO-energetic semiconductor bridge (ESCB) increased, the flame duration shortened from 440 to 300 μs, the flame size decreased by 50%, the plasma temperature decreased, and aging had no significant effect on the electrical explosion performance of Al/Cu/CuO-ESCB. After aging for 14 days in an environment with a relative humidity of 40% and a temperature of 71 °C, the maximum ignition gap of B/KNO3 for Al/CuO-ESCB decreased from 1.4 to 1.2 mm, while the maximum ignition gap for Al/Cu/CuO-ESCB remained at 1.6 mm, which significantly improved the ignition performance and long storage performance of the energetic semiconductor bridge. [ABSTRACT FROM AUTHOR]

We develop a theory to calculate structural correlations and thermodynamic properties of a fluid confined in a random porous solid medium (matrix). We used density functional formalism to derive an annealed averaged expression for the density profile and excess free energy of fluid arising due to random fields of a particular realization of the matrix. After performing the second average over the quenched-disordered variables, the excess free energy is organized to give one- and two-body potentials for fluid particles. The average over disorder reduces the system to an effective one-component system of fluid in which particles feel one-body (external) potential and interact via effective pair potential. The effective pair potential is a sum of the bare (the one in the pure fluid) and the matrix-induced potential. The resulting partition function involves only fluid variables. Equations are derived for fluid–fluid and fluid–matrix correlation functions and for free energy, pressure, and chemical potential of the fluid. The theory is applied to a model system of hard spheres and results for the effective pair potential, correlation functions, and thermodynamic properties are reported. The effective pair potential is found to be attractive at the contact and develops a repulsive peak before decaying to zero. Results for pair correlation function and structure factor are compared with simulation results for several fluid densities at two matrix densities. In all the cases, a very good agreement has been found. [ABSTRACT FROM AUTHOR]

A proof-of-concept study is undertaken to demonstrate the utility of the machine learning combined with the thermodynamic perturbation theory (MLPT) to test the accuracy of electronic structure methods in finite-temperature thermodynamic calculations. As a test example, formic acid dimer is chosen, which is one of the systems included in the popular benchmark set S22 [Jurečka et al., Phys. Chem. Chem. Phys. 8, 1985–1993 (2006)]. Starting from the explicit molecular dynamics and thermodynamic integration performed at the PBE + D2 level, the MLPT is used to obtain fully anharmonic dimerization free and internal energies at the reference quality CCSD(T) level and 19 different density functional approximations, including GGA, meta-GGA, non-local, and hybrid functionals with and without dispersion corrections. Our finite-temperature results are shown to be both qualitatively and quantitatively different from those obtained using the conventional benchmarking strategy based on fixed structures. The hybrid functional HSE06 is identified as the best performing approximate method tested, with the errors in free and internal energies of dimerization being 36 and 41 meV, respectively. [ABSTRACT FROM AUTHOR]

A fundamental equation of state in terms of the Helmholtz energy is presented for mixtures of nitrogen, oxygen, and argon at any composition. It is expressed in terms of the residual Helmholtz energy and can be used to calculate all thermodynamic equilibrium properties including vapor–liquid equilibria. The parameters of the equations for the pure-fluid and mixture contributions are fitted exclusively to molecular simulation data so that the model has a predictive character. The description of the mixture-specific reducing parameters is realized via generalized correlations of the critical parameters of the pure fluids so that an extension of the model to additional components can be implemented straightforwardly. Extensive comparisons to experimental data and the GERG-2008 reference equation of state show that the prediction of thermodynamic properties is satisfactory. [ABSTRACT FROM AUTHOR]

The ability to distinguish between stochastic systems based on their trajectories is crucial in thermodynamics, chemistry, and biophysics. The Kullback–Leibler (KL) divergence, D KL A B (0 , τ) , quantifies the distinguishability between the two ensembles of length-τ trajectories from Markov processes A and B. However, evaluating D KL A B (0 , τ) from histograms of trajectories faces sufficient sampling difficulties, and no theory explicitly reveals what dynamical features contribute to the distinguishability. This work provides a general formula that decomposes D KL A B (0 , τ) in space and time for any Markov processes, arbitrarily far from equilibrium or steady state. It circumvents the sampling difficulty of evaluating D KL A B (0 , τ). Furthermore, it explicitly connects trajectory KL divergence with individual transition events and their waiting time statistics. The results provide insights into understanding distinguishability between Markov processes, leading to new theoretical frameworks for designing biological sensors and optimizing signal transduction. [ABSTRACT FROM AUTHOR]

The focus of our study is an in-depth investigation of the quantum effects associated with the surface tension and other thermodynamic properties of nanoscopic liquid drops. The behavior of drops of quantum Lennard-Jones fluids is investigated with path-integral Monte Carlo simulations, and the test-area method is used to determine the surface tension of the spherical vapor–liquid interface. As the thermal de Broglie wavelength, λB, becomes more significant, the average density of the liquid drop decreases, with the drop becoming mechanically unstable at large wavelengths. As a consequence, the surface tension is found to decrease monotonically with λB, vanishing altogether for dominant quantum interactions. Quantum effects can be significant, leading to values that are notably lower than the classical thermodynamic limit, particularly for smaller drops. For planar interfaces (with infinite periodicity in the direction parallel to the interface), quantum effects are much less significant with the same values of λB but are, nevertheless, consequential for values representative of hydrogen or helium-4 at low temperatures corresponding to vapor–liquid coexistence. Large quantum effects are found for small drops of molecules with quantum interactions corresponding to water, ethane, methanol, and carbon dioxide, even at ambient conditions. The notable decrease in the density and tension has important consequences in reducing the Gibbs free-energy barrier of a nucleating cluster, enhancing the nucleation kinetics of liquid drops and of bubble formation. This implies that drops would form at a much greater rate than is predicted by classical nucleation theory. [ABSTRACT FROM AUTHOR]

The reversible hydrogenation/dehydrogenation of aromatic molecules, known as liquid organic hydrogen carriers (LOHCs), is considered an attractive option for the safe storage and release of elemental hydrogen. The LOHC systems based on the alkoxy-naphthalene/alkoxy-decalin studied in this work can become potentially attractive from the point of view of the thermodynamic conditions of the reversible hydrogenation/dehydrogenation processes. This work reports the results of a complex experimental investigation of the thermochemical properties of the reactants of the LOHC systems. The enthalpies of formation were measured using high-precision combustion calorimetry, the enthalpies of vaporization and sublimation were derived from the vapor pressure–temperature dependencies measured using the transpiration method, and the melting temperatures and enthalpies of fusion were measured using the differential scanning calorimetry method. The liquid-phase enthalpies of formation of methoxy- and ethoxy-substituted naphthalenes and methoxy- and ethoxy-substituted decalins were derived and used for the thermodynamic analysis of hydrogenation/dehydrogenation reactions and transferhydrogenation reactions. [ABSTRACT FROM AUTHOR]

In this study, we present a new equation of state for electrolyte solutions, integrating the statistical associating fluid theory for variable range interactions utilizing the generic Mie form and binding Debye–Hückel theories. This equation of state underscores the pivotal role of ion–ion association in determining the properties of electrolyte solutions. We propose a unified framework that simultaneously examines the thermodynamic properties of electrolyte solutions and their electrical conductivity, given the profound impact of ion pairing on this transport property. Using this equation of state, we predict the liquid density, mean ionic activity coefficient, and osmotic coefficient for binary NaCl, Na2SO4, and MgSO4 aqueous solutions at 298.15 K. Additionally, we evaluate the molar conductivity of these systems by considering the fraction of free ions derived from our equation of state in conjunction with two advanced electrical conductivity models. Our results reveal that, while ion–ion association has a minimal influence on the modification of the predicted properties of sodium chloride solutions, their impact on sodium and magnesium sulfate solutions is considerably more noticeable. [ABSTRACT FROM AUTHOR]

The influence of microstructure on the response of energetic structural materials (ESMs) under shock conditions remains inadequately quantified, and the energy release process is not thoroughly understood. In this work, taking the classical Ni/Al ESM as an example, the shock response was investigated by the shock compression theory with the microstructure-based chemical reaction kinetics model. This theory mainly refers to the equation of the state of multi-component materials with mixture rule, and the reaction at the particle contact interface is built to form the multiscale thermodynamics model. The physical states of material after shock, including relative volume, temperature, and extent of reaction, were analyzed. The results revealed the effect of the burn velocity, particle size and molar ratio on the shock response. Furthermore, the model facilitates a comprehensive understanding of energy release, the extent of the intermetallic reaction, and the oxidation reaction. Despite the involvement of only a small portion of materials in the oxidation reaction, the energy release proportion was comparable to that of the intermetallic reaction. Additionally, insights into the effect of the microstructure on the energy release revealed by the model matched the tests well. [ABSTRACT FROM AUTHOR]

Anharmonic effects play a crucial role in determining thermochemical properties of liquids and gases. For such extended phases, the inclusion of anharmonicity in reliable electronic structure methods is computationally extremely demanding, and hence, anharmonic effects are often lacking in thermochemical calculations. In this study, we apply the quantum cluster equilibrium method to transfer density functional theory calculations at the cluster level to the macroscopic, liquid, and gaseous phase of hydrogen fluoride. This allows us to include anharmonicity, either via vibrational self-consistent field calculations for smaller clusters or using a regression model for larger clusters. We obtain the structural composition of the fluid phases in terms of the population of different clusters as well as isobaric heat capacities as an example for thermodynamic properties. We study the role of anharmonicities for these analyses and observe that, in particular, the dominating structural motifs are rather sensitive to the anharmonicity in vibrational frequencies. The regression model proves to be a promising way to get access to anharmonic features, and the extension to more sophisticated machine-learning models is promising. [ABSTRACT FROM AUTHOR]

Experimental and theoretical advances for understanding the temperature of maximum density (TMD) of aqueous solutions are outlined. The main equations that relate the TMD behavior to key thermodynamic properties are stated. The experimental TMD data are classified as a function of the nature of the solute (inorganic electrolytes, non-electrolytes, organic salts and ionic liquids, and amino acids and proteins). In addition, the experimental results that explore the effect of pressure are detailed. These experimental data are rationalized by making use of qualitative and semi-quantitative arguments based on the thermodynamics of aqueous systems. The main theoretical and simulation advances in TMD for aqueous solutions are also shown—including new calculations in the context of the scaled particle theory—and their ability to reproduce the experimental data is evaluated. Finally, new experiments and theoretical and simulation developments, which could give important insights into the problem of TMD for aqueous solutions, are proposed. [ABSTRACT FROM AUTHOR]

Metal–organic frameworks (MOFs), with their unique porous structures and versatile functionality, have emerged as promising materials for the adsorption, separation, and storage of diverse molecular species. In this study, we investigate water adsorption in MOF-808, a prototypical MOF that shares the same secondary building unit (SBU) as UiO-66, and elucidate how differences in topology and connectivity between the two MOFs influence the adsorption mechanism. To this end, molecular dynamics simulations were performed to calculate several thermodynamic and dynamical properties of water in MOF-808 as a function of relative humidity (RH), from the initial adsorption step to full pore filling. At low RH, the μ3-OH groups of the SBUs form hydrogen bonds with the initial water molecules entering the pores, which triggers the filling of these pores before the μ3-OH groups in other pores become engaged in hydrogen bonding with water molecules. Our analyses indicate that the pores of MOF-808 become filled by water sequentially as the RH increases. A similar mechanism has been reported for water adsorption in UiO-66. Despite this similarity, our study highlights distinct thermodynamic properties and framework characteristics that influence the adsorption process differently in MOF-808 and UiO-66. [ABSTRACT FROM AUTHOR]

Studies of water thermodynamics have long been tied to the identification of two distinct families of local structures, whose competition could explain the origin of the many thermodynamic anomalies and the hypothesized liquid–liquid critical point in water. Despite the many successes and insights gained, the structural indicators proposed throughout the years were not able to unequivocally identify these two families over a wide range of conditions. We show that a recently introduced indicator, Ψ, which exploits information on the hydrogen bond network connectivity, can reliably identify these two distinct local environments over a wide range of thermodynamic conditions (188–300 K and 0–13 kbar) and that close to the liquid–liquid critical point, the spatial correlations of density fluctuations are identical to those of the Ψ indicator. Our results strongly support the idea that water thermodynamic properties arise from the competition between two distinct and identifiable local environments. [ABSTRACT FROM AUTHOR]

Vapor–liquid equilibria and thermodynamic properties of saturated argon and krypton were calculated by semi-classical Monte Carlo simulations with the NpT + test particle method using ab initio potentials for the two-body and nonadditive three-body interactions. The NpT + test particle method was extended to the calculation of second-order thermodynamic properties, such as the isochoric and isobaric heat capacities or the speed of sound, of the saturated liquid and vapor by using our recently developed approach for the systematic calculation of arbitrary thermodynamic properties in the isothermal–isobaric ensemble. Generally, the results for all simulated properties agree well with experimental data and the current reference equations of state for argon and krypton. In particular, the results for the vapor pressure and for the density and speed of sound of the saturated liquid and vapor agree with the most accurate experimental data for both noble gases almost within the uncertainty of these data, a level of agreement unprecedented for many-particle simulations. This study demonstrates that the vapor–liquid equilibrium and thermodynamic properties at saturation of a pure fluid can be predicted by Monte Carlo simulations with high accuracy when the intermolecular interactions are described by state-of-the-art ab initio pair and nonadditive three-body potentials and quantum effects are accounted for. [ABSTRACT FROM AUTHOR]

This study reviews available van der Waals- and Platteeuw-based hydrate models considering multiple occupancy of cavities. Small guest molecules, such as hydrogen and nitrogen, are known to occupy lattice cavities multiple times. This phenomenon has a significant impact on hydrate stability and thermodynamic properties of the hydrate phase. The objective of this work is to provide a comprehensive overview and required correlations for the implementation of a computationally sufficient cluster model that considers up to five guest molecules per cavity. Two methodologies for cluster size estimation are evaluated by existing nitrogen hydrate models showing accurate results for phase equilibria calculations. Furthermore, a preliminary hydrogen hydrate model is introduced and compared with the results of other theoretical studies, indicating that double occupancy of small sII cavities is improbable and four-molecule clusters are predominant in large sII cavities for pressures above 300 MPa. This work lays the foundation for further exploration and optimization of hydrate-based technologies for small guest molecules, e.g., storage and transportation, emphasizing their role in the future landscape of sustainable energy solutions. [ABSTRACT FROM AUTHOR]

The formation of biomolecular condensates in vivo is increasingly recognized to underlie a multitude of crucial cellular functions. Furthermore, the evolution of highly dynamic protein condensates into progressively less reversible assemblies is thought to be involved in a variety of disorders, from cancer over neurodegeneration to rare genetic disorders. There is an increasing need for efficient experimental methods to characterize the thermodynamics of condensate formation and that can be used in screening campaigns to identify and rationally design condensate modifying compounds. Theoretical advances in the field are also identifying the key parameters that need to be measured in order to obtain a comprehensive understanding of the underlying interactions and driving forces. Here, we review recent progress in the development of efficient and quantitative experimental methods to study the driving forces behind and the temporal evolution of biomolecular condensates. [ABSTRACT FROM AUTHOR]

We present a theoretical model to study the effect of counterion size on the effective charge, size, and thermodynamic behavior of a single, isolated, and flexible polyelectrolyte (PE) chain. We analyze how altering counterion size modifies the energy and entropy contributions to the system, including the ion-pair free energy, excluded volume interactions, entropy of free and condensed ions, and dipolar attraction among monomer–counterion pairs, which result in competing effects challenging intuitive predictions. The PE self-energy is calculated using the Edwards–Muthukumar Hamiltonian, considering a Gaussian monomer distribution for the PE. The condensed ions are assumed to be confined within a cylindrical volume around the PE backbone. The dipolar and excluded volume interactions are described by the second and third virial coefficients. The assumption of freely rotating dipoles results in a first-order coil–globule transition of the PE chain. A more realistic, weaker dipolar attraction, parameterized in our theory, shifts it to a second-order continuous transition. We calculate the size scaling-exponent of the PE and find exponents according to the relative dominance of the electrostatic, excluded volume, or dipolar effects. We further identify the entropy- and energy-driven regimes of the effective charge and conformation of the PE, highlighting the interplay of free ion entropy and ion-pair energy with varying electrostatic strengths. The crossover strength, dependent on the counterion size, indicates that diminishing sizes favor counterion condensation at the expense of free ion entropy. The predictions of the model are consistent with trends in simulations and generalize the findings of the point-like counterion theories. [ABSTRACT FROM AUTHOR]

Polymeric surfactants are amphiphilic molecules with two or more different types of monomers. If one type of monomer interacts favorably with a liquid, and another type of monomer interacts favorably with another, immiscible liquid, then polymeric surfactants adsorb at the interface between the two liquids and reduce the interfacial tension. The effects of polymer architecture on the structural and thermodynamic properties of the liquid–liquid interface are studied using molecular simulations. The interface is modeled with a non-additive binary Lennard-Jones fluid in the two-phase region of the phase diagram. Block and gradient copolymer surfactants are represented with coarse-grained, bead-spring models, where each component of the polymer favors one or the other liquid. Gradient copolymers have a greater concentration at the interface than do block copolymers because the gradient copolymers adopt conformations partially aligned with the interface. The interfacial tension is determined as a function of the surface excess of polymeric surfactant. Gradient copolymers are more potent surfactants than block copolymers because the gradient copolymers cross the dividing surface multiple times, effectively acting as multiple individual surfactants. For a given surface excess, the interfacial tension decreases monotonically when changing from a block to a gradient architecture. The coarse-grained simulations are complemented by all-atom simulations of acrylic-acid/styrene copolymers at the chloroform-water interface, which have been studied in experiments. The agreement between the simulations (both coarse-grained and atomistic) and experiments is shown to be excellent, and the molecular-scale structures identified in the simulations help explain the variation of surfactancy with copolymer architecture. [ABSTRACT FROM AUTHOR]

Utilization of metal tungsten (W) as the structural material or pressure scale requires accurate knowledge of the equation of state (EOS), which is far beyond the available experimental measurements. In the present work, a direct integral approach (DIA) with ultrahigh efficiency was applied to calculate the EOS of W up to 500 GPa and 3000 K with ab initio calculations. Compared with previous static compression experiments up to 150 GPa under room temperature and 35 GPa at high temperatures up to 1673 K, all the deviations of the calculated pressure are within or comparable to the uncertainty of experiments. Predictions for higher-temperature and simultaneously higher-pressure EOS up to 300 GPa and 3000 K differ slightly from the comprehensive analysis by Litasov et al. [J. Appl. Phys. 113, 133505 (2013)] via fitting available experimental data with the empirical equation. These results indicate that the EOS of crystal W obtained from DIA should be convincible, and DIA without any empirical or artificial parameters may find its wide applications for predicting thermodynamic properties of condensed matter in the future. [ABSTRACT FROM AUTHOR]

A major limitation of amine-based post-combustion carbon capture technology is the necessity to regenerate amines at high temperatures, which dramatically increases operating costs. This paper concludes the effect of solvent choice as a possible route to modify the thermodynamics and kinetics characterizing the involved amine regeneration reactions and discusses whether these modifications can be economically beneficial. We report experimentally benchmarked computational chemistry calculations of monoethanolamine regeneration reactions employing aqueous and non-aqueous solvents with a wide range of dielectric constants. Unlike previous studies, our improved computational chemistry framework could accurately reproduce the right experimental activation energy of zwitterion formation. From the thermodynamics and kinetics of the predicted reactions, the use of non-aqueous solvents with small dielectric constants led to reductions in regeneration Gibbs free energies, activation barriers, and enthalpy changes. This can reduce energy consumption and give an opportunity to run desorption columns at relatively lower temperatures, thus offering the possibility of relying on low-grade waste heat as an energy input. [ABSTRACT FROM AUTHOR]

A guided ion beam tandem mass spectrometer was employed to study the reactions of U+ + CO2, UO+ + O2, and the reverse of the former, UO+ + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U+ + CO2 proceeds to form UO+ + CO with an efficiency of 118% ± 24% as well as generating UO2+ + C and UCO+ + O. The reaction of UO+ + O2 forms UO2+ in an exothermic, barrierless process and also results in the collision-induced dissociation of UO+ to yield U+. In the UO+ + CO reaction, the formation of UO2+ in an endothermic process is the dominant reaction, but minor products of UCO+ + O and U+ + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D0(U+–O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D0(U+–CO) = 0.73 ± 0.13 eV, and D0(OU+–O) = 7.56 ± 0.12 eV. The values obtained for D0(U+–O) and D0(OU+–O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D0(U+–CO). An analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO2+ relative to ThO2+ and especially transition metal dioxide cations. [ABSTRACT FROM AUTHOR]

Chemical-ordering in metallic liquid alloys affects important structural, thermodynamic, and dynamical factors governing the kinetics of the glass formation process and the glass-forming ability. The present study on Zr50Cu50, Zr50Cu45Al5, and Zr50Cu45Ag5 metallic glass-forming liquids reveal that minor addition of Al/Ag in Zr50Cu50 leads to different chemical short-range orders due to the hetero-coordination tendency of Al in Zr50Cu45Al5 and homo-coordination tendency of Ag in Zr50Cu45Ag5. Different chemical short-range ordering causes qualitatively different topological short-range orders in the two ternary alloys. Results of inherent structure energy and excess entropy indicate modification of the potential energy landscape such that the local minima (metabasins) on the landscape of Zr50Cu45Al5 become deeper and less rugged, whereas the metabasins become shallower and more rugged in Zr50Cu45Ag5. Single-particle dynamics investigations clearly demonstrate the effect of difference in the chemical-ordering, topological short-range order, and the potential energy landscape on the atomic diffusion, structural relaxation, and dynamic heterogeneities in the ternary alloys. It signifies that the dynamics of the studied glass-forming alloys is closely linked with the structure and thermodynamics. The study also provides a very useful insight of the correlation between the chemical-ordering and the short-time dynamical features in the studied metallic glass-forming liquids. [ABSTRACT FROM AUTHOR]

Nanomaterials have excellent adsorption performance due to nano-effect and have been widely used in many fields. The adsorption performance of nanomaterials depends on their adsorption thermodynamics and kinetics. Here, the thermodynamic and kinetic equations of adsorption for spherical nanoparticles are first derived in theory. The mechanisms and regularities of influences of nano-effect on thermodynamics and kinetics, using nano-CdS adsorption of methylene blue as a probe, were investigated. The results show that there are significant influences of nano-effect (i.e., interface area effect and interface tension effect) on the thermodynamics and kinetics of nanoparticle adsorption. When the particle radius is larger and beyond the nanometer scale (∼100 nm), the nano-effect can be neglected. When the radius of the nanoparticle is within 10–50 nm, the nano-effect of adsorption is mainly the interface area effect. When the radius is less than 10 nm, the influences of both the interface area effect and the interface tension effect on adsorption thermodynamics and kinetics become more significant. The theoretic equations and experimental insights in the present work could provide a significant basis and reference for improving and perfecting the adsorption theory involving nanomaterials and provide guidelines for the explanations of adsorption mechanisms as well as the selection of adsorbents. [ABSTRACT FROM AUTHOR]

Conventional theories of weak polyelectrolytes are either computationally prohibitive to account for the multidimensional inhomogeneity of polymer ionization in a liquid environment or oversimplistic in describing the coupling effects of ion-explicit electrostatic interactions and long-range intrachain correlations. To bridge this gap, we implement the Ising density functional theory (iDFT) for ionizable polymer systems using the single-chain-in-mean-field algorithm. The single-chain-in-iDFT (sc-iDFT) shows significant improvements over conventional mean-field methods in describing segment-level dissociation equilibrium, specific ion effects, and long-range intrachain correlations. With an explicit consideration of the fluctuations of polymer configurations and the position-dependent ionization of individual polymer segments, sc-iDFT provides a faithful description of the structure and thermodynamic properties of inhomogeneous weak polyelectrolyte systems across multiple length scales. [ABSTRACT FROM AUTHOR]

Charge-stabilized colloidal suspensions display a rich variety of microstructural and thermodynamic properties, which are determined by electro-steric interactions between all ionic species. The large size asymmetry between molecular-scale microions and colloidal macroions allows the microion degrees of freedom to be integrated out, leading to an effective one-component model of microion-dressed colloidal quasi-particles. For highly charged colloids with strong macroion–microion correlations, nonlinear effects can be incorporated into effective interactions by means of charge renormalization methods. Here, we compare and partially extend several practical mean-field methods of calculating renormalized colloidal interaction parameters, including effective charges and screening constants, as functions of concentration and ionic strength. Within the one-component description, we compute structural and thermodynamic properties from the effective interactions and assess the accuracy of the different methods by comparing predictions with elaborate primitive-model simulations [P. Linse, J. Chem. Phys. 113, 4359 (2000)]. We also compare various prescriptions for the osmotic pressure of suspensions in Donnan equilibrium with a salt ion reservoir and analyze instances where the macroion effective charge becomes larger than the bare one. The methods assessed include single-center cell, jellium, and multi-center mean-field theories. The strengths and weaknesses of the various methods are critically assessed, with the aim of guiding optimal and accurate implementations. [ABSTRACT FROM AUTHOR]

Block copolymers at homopolymer interfaces are poised to play a critical role in the compatibilization of mixed plastic waste, an area of growing importance as the rate of plastic accumulation rapidly increases. Using molecular dynamics simulations of Kremer–Grest polymer chains, we have investigated how the number of blocks and block degree of polymerization in a linear multiblock copolymer impacts the interface thermodynamics of strongly segregated homopolymer blends, which is key to effective compatibilization. The second virial coefficient reveals that interface thermodynamics are more sensitive to block degree of polymerization than to the number of blocks. Moreover, we identify a strong correlation between surface pressure (reduction of interfacial tension) and the spatial uniformity of block junctions on the interface, yielding a morphological framework for interpreting the role of compatibilizer architecture (number of blocks) and block degree of polymerization. These results imply that, especially at high interfacial loading, the choice of architecture of a linear multiblock copolymer compatibilizing surfactant does not greatly affect the modification of interfacial tension. [ABSTRACT FROM AUTHOR]

Equilibrium thermodynamics describes the energy exchange of a body with its environment. Here, we describe the global energy exchange of an ideal gas in the Coutte flow in a thermodynamic-like manner. We derive a fundamental relation between internal energy as a function of parameters of state. We analyze a non-equilibrium transition in the system and postulate the extremum principle, which determines stable steady states in the system. The steady-state thermodynamic framework resembles equilibrium thermodynamics. [ABSTRACT FROM AUTHOR]

A model system of identical particles interacting via a hard-sphere potential is essential in condensed matter physics; it helps to understand in and out of equilibrium phenomena in complex fluids, such as colloidal dispersions. Yet, most of the fixed time-step algorithms to study the transport properties of those systems have drawbacks due to the mathematical nature of the interparticle potential. Because of this, mapping a hard-sphere potential onto a soft potential has been recently proposed [Báez et al., J. Chem. Phys. 149, 164907 (2018)]. More specifically, using the second virial coefficient criterion, one can set a route to estimate the parameters of the soft potential that accurately reproduces the thermodynamic properties of a monocomponent hard-sphere system. However, real colloidal dispersions are multicomponent or polydisperse, making it important to find an efficient way to extend the potential model for dealing with such kind of many-body systems. In this paper, we report on the extension and applicability of the second virial coefficient criterion to build a description that correctly captures the phenomenology of both multicomponent and polydisperse hard-sphere dispersions. To assess the accuracy of the continuous potentials, we compare the structure of soft polydisperse systems with their hard-core counterpart. We also contrast the structural and thermodynamic properties of soft binary mixtures with those obtained through mean-field approximations and the Ornstein–Zernike equation for the two-component hard-sphere dispersion. [ABSTRACT FROM AUTHOR]

There is an ever increasing use of local density dependent potentials in the mesoscale modeling of complex fluids. Questions remain, though, about the dependence of the thermodynamic and structural properties of such systems on the cutoff distance used to calculate these local densities. These questions are particularly acute when it comes to the stability and structure of the vapor/liquid interface. In this article, we consider local density dependent potentials derived from an underlying van der Waals equation of state. We use simulation and density functional theory to examine how the bulk thermodynamic and interfacial properties vary with the cutoff distance, rc, used to calculate the local densities. We show quantitatively how the simulation results for bulk thermodynamic properties and vapor–liquid equilibrium approach the van der Waals limit as rc increases and demonstrate a scaling law for the radial distribution function in the large rc limit. We show that the vapor–liquid interface is stable with a well-defined surface tension and that the interfacial density profile is oscillatory, except for temperatures close to critical. Finally, we show that in the large rc limit, the interfacial tension is proportional to rc and, therefore, unlike the bulk thermodynamic properties, does not approach a constant value as rc increases. We believe that these results give new insights into the properties of local density dependent potentials, in particular their unusual interfacial behavior, which is relevant for modeling complex fluids in soft matter. [ABSTRACT FROM AUTHOR]

A review of the structural, elastic, and thermodynamic properties of cubic and hexagonal ScxAl1−xN crystals over the range of possible random alloys is provided. Based on measured and simulated lattice and internal cell parameters of NaCl (B1), CsCl (B2), and α-ZnS (B3) type cubic ScxAl1−xN lattices as well as of β-ZnS (B4), lh-MgO (Bk), and NiAs (B81) type hexagonal ScxAl1−xN crystals, their atomic positions, distances to nearest neighbor atoms, geometric dimensions of crystal cells, mass density, as well as their average bond length and bond angles are presented in dependence on the alloy composition. The understanding gained about the crystal lattices is used to provide a model for the transitions from the β-ZnS to the lh-MgO or NaCl lattice induced by the alloying of AlN with ScN. Based on published data sets of stiffness coefficients, the compliance coefficients, Young's modulus, shear modulus, Poisson's ratio, compressibility, and the sound velocities are presented in relation to the orientation of representative crystal planes and axes for rock salt, layered hexagonal, and wurtzite ScxAl1−xN crystals. Particular attention is paid to the directional anisotropies of elastic properties of the different crystal lattices if Sc atoms substitute an increasing number of Al atoms. Based on sound velocities determined, an overview of the fundamental thermodynamic properties of cubic and hexagonal ScxAl1−xN alloys is provided, such as the Debye temperature, heat capacity, minimum heat conduction, and melting temperature. [ABSTRACT FROM AUTHOR]

This paper presents an exploration of a specific type of a generalized multicomponent solution model, which appears to be first given by Saulov in the current explicit form. The assumptions of the underlying theory and a brief derivation of the main equation have been provided preliminarily for completeness and notational consistency. The resulting formulae for the Gibbs free energy of mixing and the chemical potentials are multivariate polynomials with physically meaningful coefficients and the mole fractions of the components as variables. With one additional assumption about the relative magnitudes of the solvent–solute and solute–solute interaction exchange energies, combining rules were obtained that express the mixed coefficients of the polynomial in terms of its pure coefficients. This was done by exploiting the mathematical structure of the asymmetric form of the solvent chemical potential equation. The combining rules allow one to calculate the thermodynamic properties of the solvent with multiple solutes from binary mixture data only (i.e., each solute with the solvent), and hence, are of practical importance. Furthermore, a connection was established between the osmotic virial coefficients derived in this work and the original osmotic virial coefficients of Hill found by employing a different procedure, illustrating the equivalency of what appears to be two different theories. A validation of the combining rules derived here has been provided in a separate paper where they were successfully used to predict the freezing points of ternary salt solutions of water. [ABSTRACT FROM AUTHOR]

Thermodynamics of far-from-equilibrium systems often require measurement of effective parameters such as temperature. Whether such approach is valid for the general case of resetting protocols, active systems, or of confined systems under time-varying fields is still under investigation. We report on the effect of switching ON-OFF of an asymmetric bistable potential to the mean first passage time (MFPT) of a probed particle to go from one potential minima to the other. Experimental results coupled with numerical simulations shows the potential becoming more symmetric at slow switching. Moreover, the MFPT deviates from equilibrium condition with an effective temperature, Teff < T, at slow switching but approaches room temperature, T, at fast switching. For each switching rate, we quantify how far the system is from equilibrium by measuring deviation from a detailed balance like relation and the net circulation of flux present in phase-space. Both analysis suggest equilibrium condition are met at high switching. [ABSTRACT FROM AUTHOR]

Contrary to the case of solids and gases, where Debye theory and kinetic theory offer a good description for most of the physical properties, a complete theoretical understanding of the vibrational and thermodynamic properties of liquids is still missing. Liquids exhibit a vibrational density of states (VDOS) which does not obey Debye law, and a heat capacity which decreases monotonically with temperature, rather than growing as in solids. Despite many attempts, a simple, complete and widely accepted theoretical framework able to formally derive the aforementioned properties has not been found yet. Here, we revisit one of the theoretical proposals, and in particular we re-analyze the properties of liquids within the soft-potential model, originally formulated for glasses. We confirm that, at least at a qualitative level, many characteristic properties of liquids can be rationalized within this model. We discuss the validity of several phenomenological expressions proposed in the literature for the density of unstable modes, and in particular for its temperature and frequency dependence. We discuss the role of negative curvature regions and unstable modes as fundamental ingredients to have a linear in frequency VDOS. Finally, we compute the heat capacity within the soft potential model for liquids and we show that it decreases with temperature, in agreement with experimental and simulation data. [ABSTRACT FROM AUTHOR]

The surface charges of catalysts have intricate influences on the thermodynamics and kinetics of electrochemical reactions. Herein, we develop a grand-canonical iteration method based on density functional theory calculations to explore the effect of surface charges on reaction kinetics beyond the traditional Butler–Volmer picture. Using the hydrogen evolution reaction on S vacancies of MoS2 as an example, we show how to track the change of surface charge in a reaction and to analyze its influence on the kinetics. Protons adsorb on S vacancies in a tough and charge-insensitive water splitting manner, which explains the observed large Tafel slope. Grand-canonical calculations report an unanticipated surface charge-induced change of the desorption pathway from the Heyrovsky route to a Volmer–Tafel route. During an electrochemical reaction, a net electron inflow into the catalyst may bring two effects, i.e., stabilization of the canonical energy and destabilization of the charge-dependent grand-canonical part. On the contrary, a net outflow of electrons from the catalyst can reverse the two effects. This surface charge effect has substantial impacts on the overpotential and the Tafel slope. We suggest that the surface charge effect is universal for all electrochemical reactions and significant for those involving interfacial proton transfers. [ABSTRACT FROM AUTHOR]

The returning probability (RP) theory, a rigorous diffusion-influenced reaction theory, enables us to analyze the binding process systematically in terms of thermodynamics and kinetics using molecular dynamics (MD) simulations. Recently, the theory was extended to atomistically describe binding processes by adopting the host–guest interaction energy as the reaction coordinate. The binding rate constants can be estimated by computing the thermodynamic and kinetic properties of the reactive state existing in the binding processes. Here, we propose a methodology based on the RP theory in conjunction with the energy representation theory of solution, applicable to complex binding phenomena, such as protein–ligand binding. The derived scheme of calculating the equilibrium constant between the reactive and dissociate states, required in the RP theory, can be used for arbitrary types of reactive states. We apply the present method to the bindings of small fragment molecules [4-hydroxy-2-butanone (BUT) and methyl methylthiomethyl sulphoxide (DSS)] to FK506 binding protein (FKBP) in an aqueous solution. Estimated binding rate constants are consistent with those obtained from long-timescale MD simulations. Furthermore, by decomposing the rate constants to the thermodynamic and kinetic contributions, we clarify that the higher thermodynamic stability of the reactive state for DSS causes the faster binding kinetics compared with BUT. [ABSTRACT FROM AUTHOR]

We study the thermodynamic behavior of sodium perchlorate solutions in supercooled water through molecular dynamics numerical simulations. These solutions are of special interest because of the recent experimental results that led to hypothesize the presence of liquid water in perchlorate solutions beneath the Martian soil. We model water using the TIP4P/2005 potential. The results we obtain for solutions with concentrations 1.63 and 15.4 wt% are in agreement with those of a system undergoing a liquid–liquid phase transition where the liquid–liquid critical point shifts to slightly higher temperatures and lower pressures. The structure of the system is also analyzed, and we come to the conclusion that, even at the highest concentration considered, water retains its anomalous behavior. [ABSTRACT FROM AUTHOR]

Orientation-dependent integral equation theory (ODIET) was applied to the rose water model. Structural and thermodynamic properties of water modeled with the rose model were calculated using ODIET and compared to results from orientation-averaged integral equation theory (IET) and Monte Carlo simulations. Rose water model is a simple two-dimensional water model where molecules of water are represented as Lennard–Jones disks with explicit hydrogen bonding potential in form of rose functions. Orientational dependency significantly improves IET, as the thermodynamic results obtained using ODIET are significantly more in agreement with results calculated using MC than in the case of the orientationally averaged version. At high temperatures, the agreement between the simulation and theory is quantitative; however, when temperatures lower, a slight deviation between results obtained with different methods appear. ODIET correctly predicts the radial distribution function; moreover, ODIet also enables the calculation of angular distributions. While the angular distributions obtained with ODIET are in qualitative agreement with distributions from MC simulations, the height of the peaks in angular distributions differs between methods. Using results from ODIET, the spatial distribution of water molecules was constructed, which aids in the interpretation of other structural properties. ODIET was also used to calculate fractions of molecules with different number of hydrogen bonds, which is in the agreement with the simulations. Overall, use of ODIET significantly improves the obtained results in comparison to standard IET. [ABSTRACT FROM AUTHOR]

The computational modeling of electrochemical interfaces and their applications in electrocatalysis has attracted great attention in recent years. While tremendous progress has been made in this area, however, the accurate atomistic descriptions at the electrode/electrolyte interfaces remain a great challenge. The Computational Hydrogen Electrode (CHE) method and continuum modeling of the solvent and electrolyte interactions form the basis for most of these methodological developments. Several posterior corrections have been added to the CHE method to improve its accuracy and widen its applications. The most recently developed grand canonical potential approaches with the embedded diffuse layer models have shown considerable improvement in defining interfacial interactions at electrode/electrolyte interfaces over the state-of-the-art computational models for electrocatalysis. In this Review, we present an overview of these different computational models developed over the years to quantitatively probe the thermodynamics and kinetics of electrochemical reactions in the presence of an electrified catalyst surface under various electrochemical environments. We begin our discussion by giving a brief picture of the different continuum solvation approaches, implemented within the ab initio method to effectively model the solvent and electrolyte interactions. Next, we present the thermodynamic and kinetic modeling approaches to determine the activity and stability of the electrocatalysts. A few applications to these approaches are also discussed. We conclude by giving an outlook on the different machine learning models that have been integrated with the thermodynamic approaches to improve their efficiency and widen their applicability. [ABSTRACT FROM AUTHOR]

Synthetic multiblock copolymers are an interesting class of polymeric chains and have emerged as promising materials to mimic the function of complex biomolecules. In this work, we use Wang–Landau sampling to study sequences of multiblock ( A n B n ) m copolymers on the simple cubic lattice, where n represents the block length and m represents the number of blocks. We first compare to the thermodynamic and structural properties of four sequences previously studied in the continuum [W. Wang et al., J. Chem. Phys. 141, 244907 (2014)] to observe the differences that arise during the collapse process. We then focus on the structural transitions that occur at temperatures below the coil-to-globule transition in the lattice. Moreover, by studying additional sequences, we detail the relationship between the block length, number of blocks, and, thus, overall polymer length with respect to said structural transitions. Finally, we observe how the formation and shape of a ground state core of the more strongly interacting monomer type affect the procession of structural changes that occurs as temperature increases. [ABSTRACT FROM AUTHOR]

Recent vibrational sum frequency generation spectroscopic experiments [Sengupta et al., J. Phys. Chem. Lett. 13, 11391–11397 (2022)] demonstrated synergistic interfacial adsorption effects between the anionic dodecyl sulfate (DS−) and the polar, but charge-neutral hexaethylene glycol monododecyl ether (C12E6), surfactants. In this study, the interfacial adsorption thermodynamics underlying these synergistic effects are analyzed through free energy decompositions. A general decomposition method utilizing alchemical intermediate states is outlined. Combining free energy decompositions with the potential distribution theorem illuminates the statistical interpretations of correlated effects between different system components. This approach allows for the identification of the physical effects leading to synergistic adsorption thermodynamics of DS− binding to the air-C12E6-water interface. The binding properties are found to result from a combination of effects predominantly including energetic van der Waals stabilization between DS− and C12E6, as well as competing energetic and entropic effects due to changes in the interfacial water structure as a result of introducing a C12E6 monolayer into the bare air–water interface. [ABSTRACT FROM AUTHOR]

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. A detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology. [ABSTRACT FROM AUTHOR]

From the perspective of a bond-order-length-strength correlation, we put forward an analytical solution to describe the size, shape, and wall thickness dependency of melting temperature, entropy, and enthalpy for metallic nanostructures. Theoretical reproduction of measurements clarified that (i) when the crystal size reduces, the atomic coordination number lowers, the atomic cohesive energy decreases, and the surface-to-volume ratio increases; (ii) at the same equivalent radius, with the decrease in the number of sides for polyhedral nanoparticles and polygonal nanowires or nanotubes, the melting temperature, entropy, and enthalpy depress; and (iii) the melting temperature, entropy, and enthalpy of nanotubes are always lower than those of nanowires with the same cross-sectional radius. The present formulation is accurate and convenient, which not only shows deeper insight into the physical origins of a melting thermodynamic property response to perturbations but also provides guidance for the design and optimization of electronic nanodevices. [ABSTRACT FROM AUTHOR]

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes. [ABSTRACT FROM AUTHOR]

Coarse-grained (CG) simulation models of condensed-phase systems can be derived with well-established methods that perform coarse-graining in space and provide an effective Hamiltonian with which some of the structural and thermodynamic properties of the underlying fine-grained (FG) reference system can be represented. Coarse-graining in time potentially provides CG models that furthermore represent dynamic properties. However, systematic efforts in this direction have so far been limited, especially for moderately coarse-grained, chemistry-specific systems with complicated conservative interactions. With the aim of representing structural, thermodynamic, and dynamic properties in CG simulations of multi-component molecular systems, we investigated a recently introduced method in which the force on a CG particle originates from conservative interactions with surrounding particles and non-Markovian dissipative interactions, the latter introduced by means of a colored-noise thermostat. We examined two different methods to derive isotropic memory kernels required for integrating the corresponding generalized Langevin equation (GLE) of motion, based on the orthogonal dynamics of the FG forces and on an iterative optimization scheme. As a proof of concept, we coarse-grain single-component molecular liquids (cyclohexane, tetrachloromethane) and ideal and non-ideal binary mixtures of cyclohexane/tetrachloromethane and ethanol/tetrachloromethane, respectively. We find that for all systems, the FG single particle velocity auto-correlation functions and, consequently, both the short time and long time diffusion coefficients can be quantitatively reproduced with the CG-GLE models. We furthermore demonstrate that the present GLE-approach leads to an improved description of the rate with which the spatial correlations decay, which is artificially accelerated in the absence of dissipation. [ABSTRACT FROM AUTHOR]

In this work, we have proposed a Continuous Constant Potential Model (CCPM) based on grand canonical density functional theory for describing the electrocatalytic thermodynamics on single atom electrocatalysts dispersed on graphene support. The linearly potential-dependent capacitance is introduced to account for the net charge variation of the electrode surface and to evaluate the free energetics. We have chosen the CO2 electro-reduction reaction on single-copper atom catalysts, dispersed by nitrogen-doped graphene [CuNX@Gra (X = 2, 4)], as an example to show how our model can predict the potential-dependent free energetics. We have demonstrated that the net charges of both catalyst models are quadratically correlated with the applied potentials and, thus, the quantum capacitance is linearly dependent on the applied potentials, which allows us to continuously quantify the potential effect on the free energetics during the carbon dioxide reduction reaction instead of confining it to a specific potential. On the CuN4@Gra model, it is suggested that CO2 adsorption, coupled with an electron transfer, is a potential determining step that is energetically unfavorable even under high overpotentials. Interestingly, the hydrogen adsorption on CuN4@Gra is extremely easy to occur at both the Cu and N sites, which probably results in the reconstruction of the CuN4@Gra catalyst, as reported by many experimental observations. On CuN2@Gra, the CO2RR is found to exhibit a higher activity at the adjacent C site, and the potential determining step is shifted to the *CO formation step at a wide potential range. In general, CCPM provides a simple method for studying the free energetics for the electrocatalytic reactions under constant potential. [ABSTRACT FROM AUTHOR]

Intermetallic compounds were applied widely in the fields of automotive and aerospace because of its excellent shape memory effect. In this work, the structural, mechanical, and thermodynamic properties of Ni–Ti intermetallic compounds are studied by first-principle calculation. The results show that Ti, NiTi3, NiTi2, NiTi, Ni4Ti3, Ni3Ti, Ni4Ti, and Ni are stable based on Born stable criterion. On the mechanical aspect, the bulk modulus of Ni–Ti intermetallic compounds increases with the rising electron density. Their ductility is ranked as follows: NiTi3 > Ni4Ti3 > NiTi2 > Ni4Ti > NiTi > Ni3Ti. And when the ratio of Ni and Ti is 3:1, it is the hardest. In addition, the Ni–Ti intermetallic compounds in this work are anisotropy, except NiTi2. On the thermodynamic aspect, the Debye temperature θD, melting point Tm, and minimum thermal conductivity kmin of Ni–Ti intermetallic compounds increase first and then decrease with the rising content of Ni. The θD and kmin of NiTi3 are the minimum among the Ni–Ti ICs mentioned above, and θD of Ni3Ti and kmin of NiTi are both the maximum. The melting point of Ni3Ti is 1940 K, it is the highest temperature. [ABSTRACT FROM AUTHOR]