Grand Canonical Ensemble Research Articles

Riemann surfaces and winding numbers of Rényi phase structure of charged-flat black holes

It’s widely recognized that the free energy landscape captures the essentials of thermodynamic phase transitions. In this work, we extend the findings of [1] by incorporating the nonextensive nature of black hole entropy. Specifically, the connection between black hole phase transitions and the winding number of Riemann surfaces derived through complex analysis is extended to the Rényi entropy framework. This new geometrical and nonextensive formalism is employed to predict the phase portraits of charged-flat black holes within both the canonical and grand canonical ensembles. Furthermore, we elucidate novel relations between the number of sheets comprising the Riemann surface of the Hawking–Page and Van der Waals transitions and the dimensionality of black hole spacetimes. Notably, these new numbers are consistent with those found for charged-AdS black holes in Gibbs–Boltzmann statistics, providing another significant example of the potential connection between the cosmological constant and the nonextensive Rényi parameter.

Investigation of hybrid modeling and its transferability in building load prediction used for district heating systems

In the district heating systems, the historical operation data of the buildings in those areas would be partially or entirely missing. The traditional data-driven model is hard to predict the ground truth results because the historical data is not available for model training. However, utilizing the physics-based methods for load calculation takes a long time to process and encounters low accuracy issues. This paper investigates several hybrid models that integrate the data-driven model and the physics-based models with different fusion methods. The physics-based models calculate envelope load and infiltration load, based on Fourier's law and the grand canonical ensemble theory, respectively. After undergoing load processing, features fusion, and residual connection, the best advanced hybrid models generate 21.35%, 16.35%, and 12.73% better prediction results compared with the data-driven model. Moreover, the advanced hybride models also perform strong transferability across all the data quantity groups. In terms of practical application, the advanced hybrid models could be deployed with effective generalization in limited data scenarios and robust transfer capabilities. The selected best model constructed by hybrid modeling displays the highest performance and saves the total training costs with strong transferability.

Topological properties of black holes in five-dimensional gauged supergravity

In this paper, we study the topological properties of five-dimensional rotating charged black holes in different ensembles.The topological numbers for the black holes are gotten in the grand canonical and canonical ensembles, which are 1. When the charge and cosmological constant disappear, their topological numbers are 0. When the pressure is lower than the critical pressure, the phase transition exists in the canonical ensemble. However, the phase transition also exists in the grand canonical ensemble when the pressure is higher than the critical pressure and two independent rotational parameters are a=1 and b=−1. This may be due to the fact that the values of the rotational parameters change the supersymmetry of the black hole and lead to the phase transition.

Phantom BTZ black holes

Motivated by the impact of the phantom field (or anti-Maxwell field) on the structure of three-dimensional black holes in the presence of the cosmological constant, we present the first extraction of solutions for the phantom BTZ (A)dS black hole. In this study, we analyze the effect of the phantom field on the horizon structure. Furthermore, we compare the BTZ black holes in the presence of both the phantom and Maxwell fields. Additionally, we calculate the conserved and thermodynamic quantities of the phantom BTZ black holes, demonstrating their compliance with the first law of thermodynamics. Subsequently, we assess the effects of the electrical charge and the cosmological constant on the local stability in the canonical ensemble by considering these fields with respect to the heat capacity. We then investigate the global stability area of the BTZ black holes with phantom and Maxwell fields within the grand canonical ensemble using Gibbs free energy. In this analysis, we evaluate the influence of the electrical charge and the cosmological constant on this area.

Fleeting-Active-Site-Thrust Oxygen Evolution Reaction by Iron Cations from the Electrolyte.

Oxygen evolution reaction (OER) is key to sustainable energy and environmental engineering, thus necessitating rational design of high-performing electrocatalysts that requires understanding the structure-performance relationship with a possible dynamic nature under working conditions. Herein, we uncover a novel type of OER mechanisms thrust by the fleeting active sites (FASs) dynamically formed on Ni-based layered double hydroxides (Ni-LDHs) by Fe cations from the electrolyte under OER potentials. We employ grand-canonical ensemble methods and microkinetic modeling to elucidate the potential-dependent structures of FASs on Ni-LDHs and demonstrate that the fleeting-active-site-thrust (FAST) mechanism delivers superior OER activity via the FAST intramolecular oxygen coupling pathway, which also suppresses the lattice oxygen mechanism, leading to improved operando stability of Ni-LDHs. We further reveal that introducing only trace-level loadings (10-100 ppm) of FASs on Ni-LDHs can significantly boost and govern the catalytic performance for OER. This underscores the crucial importance of considering the novel FAST mechanism in OER and also suggests the electrolyte as a key part of the structure-performance relationship as well as an effective design strategy via engineering the electrolyte.

Selective methanol production via CO2 reduction on Cu2O revealed by micro-kinetic study combined with constant potential model

High efficiency electrocatalysis of greenhouse gas CO2 to liquid fuel methanol attracts great attentions that allows energy transformation from renewable electric energy to storable chemical energy. Catalytically converting CO2 to methanol using green hydrogen is a promising pathway to achieve carbon neutrality. In this work, the catalytic activity and product selectivity of the electrochemical CO2 reduction reaction on Cu2O were studied using a micro-kinetic model based on Marcus charge transfer theory. The constant potential model under the grand-canonical ensemble enables the accurate description of binding energies and ground-state charge transfer upon surface adsorption as a function of electrode potential. Firstly, the potential-dependent activation energies and rate constants of all possible elementary reactions involved in CO2 reduction are calculated. Then, the surface coverage and its derivative can be obtained in a steady-state approximation. Finally, the total reaction rate and partial reaction rate can be deduced, allowing for direct comparison with catalytic activity (current density) and product selectivity (Faradaic efficiency). It is found that CH3OH is the dominant product on the Cu2O(111) surface, while CO is the main product on the Cu2O(110) surface. Based on the derived total and partial reaction rates under the steady-state approximation, it is demonstrated that Cu2O(111) exhibits the highest electrocatalytic activity and methanol selectivity at U = −0.6 V vs. SHE. This research presents a methodology for predicting apparent electrocatalytic performance from microscopic reaction energy calculations, which can be applied to other important electrocatalytic reaction mechanisms and performance predictions.

Semi-empirical supported, Ab initio derived thermodynamic properties for ClO2 and its sub and extended species, applied in water treatment cycles

This paper describes a group of sixty (60) sub and extended chlorine oxide species with the general formulae of ClxOy (with x ≤ 2, y ≤ 8). Their role in water treatment cycles, behaving as key reactive species, is represented by a complex sequence of chemical inter-dependencies, exposed as a cohesive set of chemical reactions to demonstrate their cyclic role in aqueous media. An empirical/semi-empirical computational approach, supported by Ab Initio simulations, in accordance with open-shell character, has been followed to determine their optimum molecular geometries, to obtain their thermochemical properties. Besides a single molecular analysis, Grand Canonical Ensemble simulations, supported by a revised library of force field parameters, constituted a core component of the computational approach and proved to be invaluable in confirming thermochemical properties. This approach also offered finite estimates of optimum model sizes, a benefit with wider modelling application. Extended molecular species of ClO2 display a complex sequence of bonding character, with a variable charge dissipation (reported as partial charges), which complicates selection of basis sets in optimizing molecular geometries, during Ab Initio analyses.Optimum molecular geometries were obtained using Gaussian and MOPAC, which in turn resulted in reliable Heats of Formation. These correlated well with energies extracted from the open literature. Thermodynamic Analysis of the reaction of selected chlorine oxides with water using FactSage, predicted the production of known and two previously undetected chlorine species, [ClO4]- and [ClOH2]+.

Statistical Thermodynamic Description of Self-Assembly of Large Inclusions in Biological Membranes.

Recent studies revealed anomalous underscreening in concentrated electrolytes, and we suggest that the underscreened electrostatic forces between membrane proteins play a significant role in the process of self-assembly. In this work, we assumed that the underscreened electrostatic forces compete with the thermodynamic Casimir forces induced by concentration fluctuations in the lipid bilayer, and developed a simplified model for a binary mixture of oppositely charged membrane proteins with different preference to liquid-ordered and liquid-disordered domains in the membrane. In the model, like macromolecules interact with short-range Casimir attraction and long-range electrostatic repulsion, and the cross-interaction is of the opposite sign. We determine energetically favored patterns in a system in equilibrium with a bulk reservoir of the macromolecules. Different patterns consisting of clusters and stripes of the two components and of vacancies are energetically favorable for different values of the chemical potentials. Effects of thermal flutuations at low temperature are studied using Monte Carlo simulations in grand canonical and canonical ensembles. For fixed numbers of the macromolecules, a single two-component cluster with a regular pattern coexists with dispersed small one-component clusters, and the number of small clusters depends on the ratio of the numbers of the molecules of the two components. Our results show that the pattern formation is controlled by the shape of the interactions, the density of the proteins, and the proportion of the components.

Discharging of Ramsdellite MnO2 Cathode in a Lithium-Ion Battery.

We report an application of our unbiased Monte Carlo approach to investigate thermodynamic and electrochemical properties of lithiated manganese oxide in the ramsdellite phase (R-MnO2) to uncover the mechanism of lithium intercalation and understand charging/discharging of R-MnO2 as a cathode material in lithium-ion batteries. The lithium intercalation reaction was computationally explored by modeling thermodynamically significant distributions of lithium and reduced manganese in the R-MnO2 framework for a realistic range of lithium molar fractions 0 < x < 1 in Li x MnO2. We employed interatomic potentials and analyzed the thermodynamics of the resultant grand canonical ensemble. We found ordered or semiordered phases at x = 0.5 and 1.0 in Li x MnO2, verified by configurational entropy changes and simulated X-ray diffraction patterns of partially lithiated R-MnO2. The radial distribution functions show the preference of lithium for homogeneous distributions across the one-dimensional 2 × 1 ramsdellite channels accompanied by alternating reduced/oxidized manganese ions. The occupation of the interstitial sites in the channels is correlated with the calculated voltage profile, showing a sharp voltage drop at x = 0.5, which is explained by the energy penalty of shifting lithium ions from stable tetrahedral to unstable octahedral sites. To facilitate this work, our in-house software, Knowledge Led Master Code (KLMC) was extended to support massive parallelism on high-performance computers.

Partition function zeros of zeta-urns

We discuss the distribution of partition function zeros for the grand-canonical ensemble of the zeta-urn model, where tuning a single parameter can give a first or any higher order condensation transition. We compute the locus of zeros for finite-size systems and test scaling relations describing the accumulation of zeros near the critical point against theoretical predictions for both the first and higher order transition regimes.

Plasma-Plasma Third Order Phase Transition from Type IIB Supergravity.

We show that the planar, charged black hole in asymptotically anti-de Sitter spacetime, dual to the strongly coupled quark-gluon plasma thermal state of large N, SU(N), N=4 super Yang-Mills at finite chemical potential undergoes a third-order phase transition in the grand canonical ensemble to a hairy black hole of type IIB supergravity. The hairy phase is another strongly coupled fluid with a conformal equation of state and can be interpreted as another kind of quark-gluon plasma. This new quark-gluon plasma has less entropy and, therefore, seems to characterize some form of smooth hadronization. The locus of the transition in terms of the critical values of the "baryon" chemical potential μ_{c} and the temperature T_{c} is μ_{c}=2πT_{c}.

RASPA3: A Monte Carlo code for computing adsorption and diffusion in nanoporous materials and thermodynamics properties of fluids.

We present RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. In this article, we introduce the new additions and changes compared to RASPA2. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub. Using this code, we compare four Monte Carlo (MC) insertion/deletion techniques: unbiased Metropolis MC, Configurational-Bias Monte Carlo (CBMC), Continuous Fractional Component MC (CFCMC), and CB/CFCMC. We compare particle distribution shapes, acceptance ratios, accuracy and speed of isotherm computation, enthalpies of adsorption, and chemical potentials, over a wide range of loadings and systems, for the grand canonical ensemble and for the Gibbs ensemble.

Accelerated convergence via adiabatic sampling for adsorption and desorption processes.

Under isothermal conditions, phase transitions occur through a nucleation event when conditions are sufficiently close to coexistence. The formation of a nucleus of the new phase requires the system to overcome a free energy barrier of formation, whose height rapidly rises as supersaturation decreases. This phenomenon occurs both in the bulk and under confinement and leads to a very slow kinetics for the transition, ultimately resulting in hysteresis, where the system can remain in a metastable state for a long time. This has broad implications, for instance, when using simulations to predict phase diagrams or screen porous materials for gas storage applications. Here, we leverage simulations in an adiabatic statistical ensemble, known as adiabatic grand-isochoric ensemble (μ, V, L) ensemble, to reach equilibrium states with a greater efficiency than its isothermal counterpart, i.e., simulations in the grand-canonical ensemble. For the bulk, we show that at low supersaturation, isothermal simulations converge slowly, while adiabatic simulations exhibit a fast convergence over a wide range of supersaturation. We then focus on adsorption and desorption processes in nanoporous materials, assess the reliability of (μ, V, L) simulations on the adsorption of argon in IRMOF-1, and demonstrate the efficiency of adiabatic simulations to predict efficiently the equilibrium loading during the adsorption and desorption of argon in MCM-41, a system that exhibits significant hysteresis. We provide quantitative measures of the increased rate of convergence when using adiabatic simulations. Adiabatic simulations explore a wide temperature range, leading to a more efficient exploration of the configuration space.

Exploring the UV and IR of a type-II holographic superconductor using a dyonic black hole

In this study, we investigate a type-II holographic superconductor with a perturbative scalar field over a (3 + 1)-dimensional electric and magnetically charged planar AdS black hole. We review the thermodynamical properties of the background relevant for the dual description of the Ginzburg–Landau density of superconducting states, showing that the adoption of a London gauge allows a consistent description of the Abrikosov vortex lattices, typical of type-II superconductors. We further derive a new expression for the upper critical magnetic field as a function of temperature in the grand canonical ensemble, supplementary to the one previously obtained in the canonical ensemble setup. These results confirm that our perturbative scalar field model consistently reproduces the well-known temperature behavior of the upper critical magnetic field according to the Ginzburg–Landau theory and other Abelian–Higgs holographic developments for type-II superconductors. We then continue with an original analysis of the scalar field equation in terms of a Schrödinger potential that led to the observation of bound states, a signal for the condensation of Cooper pairs. These results provide robust evidence for the existence of an IR order parameter near extremality. In view of this, we performed a closer inspection of the IR effective scalar equation in which the geometry adopts a Schwarzschild AdS2×R2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {AdS}_{2}\ imes {\\mathbb {R}}^{2}$$\\end{document} structure, finding novel condensed bound states in the IR regime.

An Inverse Cluster Expansion for the Chemical Potential

Interacting particle systems in a finite-volume in equilibrium are often described by a grand-canonical ensemble induced by the corresponding Hamiltonian, i.e. a finite-volume Gibbs measure. However, in practice, directly measuring this Hamiltonian is not possible, as such, methods need to be developed to calculate the Hamiltonian potentials from measurable data. In this work, we give an expansion of the chemical potential in terms of the correlation functions of such a system in the thermodynamic limit. This is a justification of a formal approach of Nettleton and Green from the 50’s, that can be seen as an inverse cluster expansion.

Topology of charged AdS black hole in restricted phase space* *Supported by the National Natural Science Foundation of China (12075143, 12375050), and the Natural Science Foundation of Shanxi Province, China (202203021221209, 202303021211180)

The local topological properties of black hole systems can be expressed through winding numbers as defects. To date, AdS black hole thermodynamics are often depicted by the dual parameters of in the extended phase space, while there have been several studies on the black hole thermodynamics in the restricted phase space. In this paper, we analyze the topological properties of charged AdS black holes in the restricted phase space under the higher-dimension and higher-order curvature gravity frame. The results show that the topological number of the charged black hole in the same canonical ensembles is a constant and is independent of the concrete dual thermodynamical parameters. However, the topological number in the grand canonical ensemble is different from that in the canonical ensemble for the same black hole system. Furthermore, these results are independent of dimension d and highest order k of the Lanczos-Lovelock densities.

Hyperdensity Functional Theory of Soft Matter.

We present a scheme for investigating arbitrary thermal observables in spatially inhomogeneous equilibrium many-body systems. Extending the grand canonical ensemble yields any given observable as an explicit hyperdensity functional. Associated local fluctuation profiles follow from an exact hyper-Ornstein-Zernike equation. While the local compressibility and simple observables permit analytic treatment, complex order parameters are accessible via simulation-based supervised machine learning of neural hyperdirect correlation functionals. We exemplify efficient and accurate neural predictions for the cluster statistics of hard rods, square-well rods, and hard spheres. The theory allows one to treat complex observables, as is impossible in standard density functional theory.

Type II string theory on AdS3×S3×T4 and symmetric orbifolds

We discuss in detail the 1+1-dimensional superconformal field theory dual to type II string theory on AdS3×S3×T4, emphasizing the string theoretic aspects of this duality. For one unit of Neveu-Schwarz (NS-NS) 5-brane flux (Q5=1), this string theory has been suggested to be dual to a grand-canonical ensemble of T4N/SN free symmetric orbifold conformal field theories (CFTs). We show how the string genus expansion emerges to all orders for the free orbifold grand-canonical correlation functions. We also discuss how the strong coupling limit of the NS-NS string theory arises (even at large N) in the free orbifold description, and argue why this limit does not have a weakly coupled RR description. The dual conformal field theory (CFT) includes (for all values of Q5) an extra T4 factor that is decoupled from perturbative string theory. We discuss the exactly marginal deformations that relate the different values of Q5, including the precise JJ¯ deformations mixing this extra T4 with the symmetric orbifold. Published by the American Physical Society 2024

Theory and Monte Carlo simulation of the ideal gas with shell particles in the canonical, isothermal-isobaric, grand canonical, and Gibbs ensembles.

Theories of small systems play an important role in the fundamental understanding of finite size effects in statistical mechanics, as well as the validation of molecular simulation results as no computer can simulate fluids in the thermodynamic limit. Previously, a shell particle was included in the isothermal-isobaric ensemble in order to resolve an ambiguity in the resulting partition function. The shell particle removed either redundant volume states or redundant translational degrees of freedom of the system and yielded quantitative differences from traditional simulations in this ensemble. In this work, we investigate the effect of including a shell particle in the canonical, grand canonical, and Gibbs ensembles. For systems comprised of a pure component ideal gas, analytical expressions for various thermodynamic properties are obtained. We also derive the Metropolis Monte Carlo simulation acceptance criteria for these ensembles with shell particles, and the results of the simulations of an ideal gas are in excellent agreement with the theoretical predictions. The system size dependence of various important ensemble averages is also analyzed.

Corrections of Electron–Phonon Coupling for Second‐Order Structural Phase Transitions

Structural phase transitions are accompanied by a movement of one nucleus (or a few) in the crystallographic unit cell. If the nucleus movement is continuous, a second‐order phase transition without latent heat results, whereas an abrupt nucleus displacement indicates a first‐order phase transition with accompanying latent heat. Herein, a Hamiltonian including electron–phonon coupling (EPC) as proposed by Kristoffel and Konsel is taken. Contrary to their treatment, both the kinetic energy of the nucleus and its position are treated. The interaction of the many‐electron system with the single nucleus is taken into account by the Born–Oppenheimer approximation and perturbative expressions for the free energies are derived. The nuclei corrections due to the entangled electrons are found to be minor, but highlight the importance of the symmetry breaking at low temperature. Furthermore the free energy for a canonical ensemble is computed, whereas Kristoffel and Konsel use a grand canonical ensemble, which allows to derive more stringent bounds on the free energy. For the zero‐order nucleus correction, the shift of the phase transition temperature by evaluating the free energy is deduced.