This content will become publicly available on January 1, 2025
 Award ID(s):
 2305155
 NSFPAR ID:
 10503541
 Publisher / Repository:
 arXiv
 Date Published:
 Journal Name:
 Physical Chemistry Chemical Physics
 ISSN:
 14639076
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

We present Thermobar, a new opensource Python3 package for calculating pressures, temperatures, and melt compositions from mineral and mineralmelt equilibrium. Thermobar allows users to perform calculations with >100 popular parametrizations involving liquid, olivineliquid, olivinespinel, pyroxene only, pyroxeneliquid, two pyroxene, feldsparliquid, two feldspar, amphibole only, amphiboleliquid, and garnet equilibria. Thermobar is the first opensource tool which can match up all possible pairs of phases from a given region, and apply various equilibrium tests to identify pairs from which to calculate pressures and temperatures (e.g. pyroxeneliquid, two pyroxene, feldsparliquid, two feldspar, amphiboleliquid). Thermobar also contains functions allowing users to propagate analytical errors using MonteCarlo methods, convert pressures to depths using different crustal density profiles, plot mineral classification and mineralmelt equilibrium diagrams, calculate liquid viscosities, and convert between oxygen fugacity values, buffer positions and Fe speciation in a silicate melt. Thermobar can be downloaded using pip and extensive documentation is available at https://thermobar.readthedocs.io/.more » « less

null (Ed.)We present a quasipolynomial time classical algorithm that estimates the partition function of quantum manybody systems at temperatures above the thermal phase transition point. It is known that in the worst case, the same problem is NPhard below this point. Together with our work, this shows that the transition in the phase of a quantum system is also accompanied by a transition in the hardness of approximation. We also show that in a system of n particles above the phase transition point, the correlation between two observables whose distance is at least Ω(logn) decays exponentially. We can improve the factor of logn to a constant when the Hamiltonian has commuting terms or is on a 1D chain. The key to our results is a characterization of the phase transition and the critical behavior of the system in terms of the complex zeros of the partition function. Our work extends a seminal work of Dobrushin and Shlosman on the equivalence between the decay of correlations and the analyticity of the free energy in classical spin models. On the algorithmic side, our result extends the scope of a recent approach due to Barvinok for solving classical counting problems to quantum manybody systems.more » « less

Density functional theory (DFT) prescribes the existence of a groundstate configuration at 0 K for a given system. However, the groundstate configuration alone is insufficient to describe a phase at finite temperatures as symmetrybreaking nongroundstate configurations are excited statistically at temperatures above 0 K. Our multiscale entropy approach, Zentropy, postulates that the entropy of a phase is composed of the sum of the entropy of each configuration weighted by its probability plus the configurational entropy among all configurations. Consequently, the partition function of each configuration in statistical mechanics needs to be evaluated by its free energy rather than total energy. The combination of the groundstate and symmetrybreaking nongroundstate configurations represents the building blocks of materials and can be used to quantitatively predict free energy of individual phases with the free energy of each configuration predicted from DFT as well as all properties derived from free energy of individual phases.more » « less

de Jong, Bert ; Nieminen, Risto (Ed.)A kernel polynomial method is developed to calculate the random phase approximation (RPA) correlation energy. In the method, the RPA correlation energy is formulated in terms of the matrix that is the product of the Coulomb potential and the density linear response functions. The integration over the matrix's eigenvalues is calculated by expanding the density of states of the matrix in terms of the Chebyshev polynomials. The coefficients in the expansion are obtained through stochastic sampling. Since it is often the energy difference between two systems that is of much interest in practice, another focus of this work is to develop a correlated sampling scheme to accelerate the convergence of the stochastic calculations of the RPA correlation energy difference between two similar systems. The scheme is termed the atombased correlated sampling (ACS). The performance of ACS is examined by calculating the isomerization energy of acetone to 2propenol and the energy of the water–gas shift reaction. Using ACS, the convergences of these two examples are accelerated by 3.6 and 4.5 times, respectively. The methods developed in this work are expected to be useful for calculating RPAlevel reaction energies for the reactions that take place in local regions, such as calculating the adsorption energies of molecules on transition metal surfaces for modeling surface catalysis.more » « less

A theoretical treatment based on the equations of motion of an electronic reduced density matrix, and related computational modeling, is used to describe and calculate relaxation times for nanostructured TiO_{2}(110) surfaces, here for Ag and Ag_{2}adsorbates. The theoretical treatment deals with the preparation of a photoexcited system under two different conditions, by steady light absorption with a cutoff and by a light pulse, and describes the following relaxation of electronic densities. On the computational modeling, results are presented for electronic density of states, light absorbance, and relaxation dynamics, comparing results for Ag and Ag_{2}adsorbates. The aim of this work is to provide insight on the dynamics and magnitude of relaxation rates for a surface with adsorbed open and closedshell Ag species to determine whether the advantages in using them to enhance light absorbance remain valid in the presence of charge density relaxation. Different behaviors can be expected depending on whether the adsorbate particles (Ag metal clusters in our present choice) have electronic openshell or closedshell structures. Calculated electron and hole lifetimes are given for pure TiO_{2}(110), Ag/TiO_{2}(110), and Ag_{2}/TiO_{2}(110). The present results, while limited to chosen structures and photon wavelengths, show that relaxation rates are noticeably different for electrons and holes, but comparable in magnitude for pure and adsorbate surfaces. Overall, the introduction of the adsorbates does not lead to rapid loss of charge carriers, while they give large increases in light absorption. This appears to be advantageous for applications to photocatalysis.