Date: Fri, 9 Aug 2024 | Total: 73

Calcium carbonate plays a crucial role in the global carbon cycle, and its phase diagram has always been of significant scientific interest. In this study, we used molecular dynamics (MD) to investigate several structural phase transitions of calcium carbonate. Using the Raiteri potential model, we explored the structural transitions occurring at a constant pressure of 1 bar with temperatures ranging from 300 K to 2500 K, and at a constant temperature of 1600 K with pressures ranging from 0 to 13 GPa. At increasing temperatures, the transitions calcite, CaCO$_3$-IV, and CaCO$_3$-V are observed and characterized. Within the calcite structure, CO$_3^{2-}$ ions are ordered between layers. As temperature increases, the calcite to CaCO$_3$-IV transition occurs, determining the partial disordering of CO$_3^{2-}$ ions. At a higher temperature, CaCO$_3$-IV transforms into CaCO$_3$-V. By applying free energy analysis, we have classified the last transition as a continuous order-disorder phase transition. At a temperature of 2000 K, it appears a `disordered CaCO$_3$' structure, characterized by low order within the calcium and carbonate sublattices and the free rotation of CO$_3^{2-}$ ions. At increasing pressures, two calcium carbonate transformations were observed. At $P=$ 2 GPa, the CaCO$_3$-IV phase undergoes a phase transition into CaCO$_3$-V, demonstrating that the model can describe the transition between these two phases as pressure and temperature-driven. Another phase transition was detected at $P=$ 4.25 GPa -- CaCO$_3$-V transits into the recently discovered CaCO$_3$-Vb phase. This transition is classified as a first-order phase transition by structural analysis and free energy-based arguments.

Elucidating the fundamental limitations inherent in physical systems is a central subject in physics. For important thermodynamic operations such as information erasure, cooling, and copying, resources like time and energetic cost must be expended to achieve the desired outcome within a predetermined error margin. In this study, we introduce the concept of {\it separated states}, which consist of fully unoccupied and occupied states. This concept generalizes many critical states involved in relevant thermodynamic operations. We uncover a three-way trade-off relation between {\it time}, {\it cost}, and {\it error} for a general class of thermodynamic operations aimed at creating separated states, simply expressed as $\tau{\cal C}\varepsilon_{\tau}\ge 1-\eta$. This fundamental relation is applicable to diverse thermodynamic operations, including information erasure, cooling, and copying. It provides a profound quantification of the unattainability principle in the third law of thermodynamics in a general form. Building upon this relation, we explore the quantitative limitations governing cooling operations, the preparation of separated states, and a no-go theorem for exact classical copying. Furthermore, we extend these findings to the quantum regime, encompassing both Markovian and non-Markovian dynamics. Specifically, within Lindblad dynamics, we derive a similar three-way trade-off relation that quantifies the cost of achieving a pure state with a given error. The generalization to general quantum dynamics involving a system coupled to a finite bath implies that heat dissipation becomes infinite as the quantum system is exactly cooled down to the ground state or perfectly reset to a pure state, thereby resolving an open question regarding the thermodynamic cost of information erasure.

From a data perspective, the field of materials mechanics is characterized by a sparsity of available data, mainly due to the strong microstructure-sensitivity of properties such as strength, fracture toughness, and fatigue limit. Consequently, individual tests are needed for specimens with various thermo-mechanical process histories, even if their chemical composition remains the same. Experimental data on the mechanical behavior of materials is usually rare, as mechanical testing is typically a destructive method requiring large amounts of material and effort for specimen preparation and testing. Furthermore, mechanical behavior is typically characterized in simplified tests under uniaxial loading conditions, whereas a complete characterization of mechanical material behavior requires multiaxial testing conditions. To address this data sparsity, different simulation methods, such as micromechanical modeling or even atomistic simulations, can contribute to microstructure-sensitive data collections. These methods cover a wide range of materials with different microstructures characterized under multiaxial loading conditions. In the present work, we describe a novel data schema that integrates metadata and mechanical data itself, following the workflows of the material modeling processes by which the data has been generated. Each run of this workflow results in unique data objects due to the incorporation of various elements such as user, system, and job-specific information in correlation with the resulting mechanical properties. Hence, this integrated data format provides a sustainable way of generating data objects that are Findable, Accessible, Interoperable, and Reusable (FAIR). The choice of metadata elements has centered on necessary features required to characterize microstructure-specific data objects, simplifying how purpose-specific datasets are collected by search algorithms.

In this study, we uncover the intrinsic information processes in non-Hermitian quantum systems and their thermodynamic effects. We demonstrate that these systems can exhibit negative entropy production, making them potential candidates for information engines. We also identify a key informational quantity that can characterize phase transitions beyond the reach of traditional partition functions. This work enhances our understanding of the interplay between information and thermodynamics, providing a new perspective on non-Hermitian quantum systems.

Discovering mixed state quantum orders is an on-going issue. Recently, it has been recognized that there are (at least) two kinds of symmetries in the mixed state; strong and weak symmetries. Under symmetry-respective decoherence, spontaneous strong-to-weak symmetry breaking (SSSB) can occur. This work provides a scheme to describe SSSB and other decoherence phenomena in the mixed state by employing the stabilizer formalism and the efficient numerical algorithm of Clifford circuits. We present two systematic numerical studies. In a two-dimensional (2D) circuit with a stochastic Ising type decoherence, an SSSB phase transition is clearly observed and its criticality is elucidated by the numerical methods. In particular, we calculate R\'{e}nyi-2 correlations and estimate critical exponents of the SSSB transition. For the second system, we introduce an idea of subsystem SSSB. As an example, we study a system with symmetry-protected-topological (SPT) order provided by both one-form and zero-form symmetries, and observe how the system evolves under decoherence. After displaying numerical results, we show that viewpoint of percolation is quite useful to understand the SSSB transition, which is applicable for a wide range of decohered states. Finally, we comment on SSSB of one-form-symmetry exemplifying toric code.

Modelling electrolytes accurately on both a nanoscale and cell level can contribute to improving battery chemistries.[Armand and Tarascon, Nature, 2008, 451, 652-657] We previously presented a thermodynamic continuum model for electrolytes.[arXiv:2010.14915] In this paper we include solvation interactions between the ions and solvent, which alter the structure of the electochemical double layer (EDL). We are able to combine a local solvation model -- permitting examination of the interplay between electric forces and the ion-solvent binding -- with a full electrolyte model. Using this, we can investigate double layer structures for a wide range of electrolytes, especially including highly concentrated solutions. We find that some of the parameters of our model significantly affect the solvent concentration at the electrode surface, and thereby the rate of solvent decomposition. Firstly, an increased salt concentration weakens the solvation shells, making it possible to strip the solvent in the EDL before the ions reach the surface. The strength of the ion-solvent interaction also affects at which potential difference the solvation shells removed. We are therefore able to qualitatively predict EDL structures for different electrolytes based on parameters like molecule size, solvent binding energy and salt concentration.

We prove that in strongly disordered, interacting, quantum chains, the conductance of a chain of length $L$ vanishes faster than $1/L$. This means that transport is anomalous in such chains. This phenomenon was first claimed in \cite{Basko2006,gornyi2005interacting} and a pioneering treatment appeared in \cite{Imbrie2016a}

We consider a composite defect system where a lower-dimensional defect (sub-defect) is embedded to a higher-dimensional one, and examine renormalization group (RG) flows localized on the defect. A composite defect is constructed in the $(3-\epsilon)$-dimensional free $\text{O}(N)$ vector model with line and surface interactions by triggering localized RG flows to non-trivial IR fixed points. Focusing on the case where the symmetry group $\text{O}(N)$ is broken to a subgroup $\text{O}(m)\times\text{O}(N-m)$ on the line defect, there is an $\text{O}(N)$ symmetric fixed point for all $N$, while two additional $\text{O}(N)$ symmetry breaking ones appear for $N\ge 23$. We also examine a $\mathcal{C}$-theorem for localized RG flows along the sub-defect and show that the $\mathcal{C}$-theorem holds in our model by perturbative calculations.

Water's ability to autoionize into hydroxide and hydronium ions profoundly influences surface properties, rendering interfaces either basic or acidic. While it is well-established that the water-air interface is acidic, a critical knowledge gap exists in technologically relevant surfaces like the graphene-water interface. Here we use machine learning-based simulations with first-principles accuracy to unravel the behavior of the hydroxide and hydronium ions at the graphene-water interface. Our findings reveal that the graphene-water interface is acidic, with the hydronium ion predominantly residing in the first contact layer of water. In contrast, the hydroxide ion exhibits a bimodal distribution, found both near the surface and towards the interior layers. Analysis of the underlying electronic structure reveals strong polarization effects, resulting in counterintuitive charge rearrangement. Proton propensity to the graphene-water interface challenges the interpretation of surface experiments and is expected to have far-reaching consequences for ion conductivity, interfacial reactivity, and proton-mediated processes.

Recently it was shown that the topological entanglement entropy (TEE) of a two-dimensional gapped ground state obeys the universal inequality $\gamma \geq \log \mathcal{D}$, where $\gamma$ is the TEE and $\mathcal{D}$ is the total quantum dimension of all anyon excitations, $\mathcal{D} = \sqrt{\sum_a d_a^2}$. Here we present an alternative, more direct proof of this inequality. Our proof uses only the strong subadditivity property of the von Neumann entropy together with a few physical assumptions about the ground state density operator. Our derivation naturally generalizes to a variety of systems, including spatially inhomogeneous systems with defects and boundaries, higher dimensional systems, and mixed states.

Recent years have seen substantial advances in our understanding of high-dimensional ridge regression, but existing theories assume that training examples are independent. By leveraging recent techniques from random matrix theory and free probability, we provide sharp asymptotics for the in- and out-of-sample risks of ridge regression when the data points have arbitrary correlations. We demonstrate that in this setting, the generalized cross validation estimator (GCV) fails to correctly predict the out-of-sample risk. However, in the case where the noise residuals have the same correlations as the data points, one can modify the GCV to yield an efficiently-computable unbiased estimator that concentrates in the high-dimensional limit, which we dub CorrGCV. We further extend our asymptotic analysis to the case where the test point has nontrivial correlations with the training set, a setting often encountered in time series forecasting. Assuming knowledge of the correlation structure of the time series, this again yields an extension of the GCV estimator, and sharply characterizes the degree to which such test points yield an overly optimistic prediction of long-time risk. We validate the predictions of our theory across a variety of high dimensional data.

We establish a framework for a mechanism steering state vector collapse through time evolution as described by the Schr\"odinger's equation. Under the premise of a non-Hermitian Hamiltonian dynamics, a main outcome is that collapse has to occur when the Hamiltonian completes a closed loop in the parameter space encoding the interaction with the meter. For two-level systems, we put forward the phenomenon of chiral state conversion as a mechanism effectively eliminating superpositions. Our results open a way to simulate quantum measurements in classical systems that up to now were restricted to the Schr\"odinger part of the quantum dynamics.

Orthorhombic pseudobrookites CoTi$_2$O$_5$ and FeTi$_2$O$_5$ have a low-symmetry crystal structure comprising magnetic Co$^{2+}$/Fe$^{2+}$ ions and nonmagnetic Ti$^{4+}$ ions, where the orbital-nondegenerate Co$^{2+}$/Fe$^{2+}$ ions form one-dimensional chains running along the orthorhombic $a$ axis. These compounds undergo an antiferromagnetic phase transition at $T_N \sim$ 26 K for CoTi$_2$O$_5$ and $T_N \sim$ 40 K for FeTi$_2$O$_5$. We perform ultrasound velocity measurements on single crystals of CoTi$_2$O$_5$ and FeTi$_2$O$_5$. The measurements of these compounds reveal that the symmetry-lowering elastic modes of shear elastic moduli exhibit unusual elastic softness in the paramagnetic phase above $T_N$. This elastic softness indicates the presence of spin-lattice-coupled fluctuations above $T_N$ that should be a precursor to the symmetry-lowering lattice distortion at $T_N$. Furthermore, it is revealed that the magnitude of the unusual elastic softness is larger in CoTi$_2$O$_5$ than in FeTi$_2$O$_5$, which indicates that the spin-lattice coupling is stronger in CoTi$_2$O$_5$ than in FeTi$_2$O$_5$. The present study suggests that CoTi$_2$O$_5$ and FeTi$_2$O$_5$ are unique spin Jahn--Teller systems with low crystal symmetry, where, although the nature of exchange interactions is quasi-one-dimensional, the three-dimensional spin-lattice coupling releases the frustration by further lowering the crystal symmetry.

The Hulth\'{e}n potential is a short-range potential that has been widely used in various fields of physics. In this paper, we investigate the distribution functions for the Hulth\'{e}n potential by using statistical and superstatistical methods. We first review the ordinary statistics and superstatistics methods. We then consider some distribution functions, such as uniform, 2-level, gamma, and log-normal and F distributions. Finally, we investigate the behavior of the Hulth\'{e}n potential for statistical and superstatistical methods and compare the results with each other. We use the Tsallis statistics of the superstatistical system. We conclude that the Tsallis behavior of different distribution functions for the Hulth\'{e}n potential exhibits better results than the statistical method. We examined the thermal properties of the Hulth\'{e}n potential for five different distributions: Uniform, 2-level, Gamma, Log-normal, and F. We plotted the Helmholtz free energy and the entropy as functions of temperature for various values of q. It shows that the two uniform and 2-level distributions have the same results due to the universal relationship and that the F distribution does not become ordinary statistics at q=1. It also reveals that the curves of the Helmholtz free energy and the entropy change their order and behavior as q increases and that some distributions disappear or coincide at certain values of q. One can discuss the physical implications of our results and their applications in nuclear and atomic physics in the future.

The widely used thermal Hartree-Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is defined by the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasi-particle) energies in such a way that all thermodynamic relations are obeyed. The associated density matrix is of the one-electron type, whose diagonal elements take the form of the Fermi-Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi-Dirac or thermal HF theory. The theory thus postulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green's function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperature, the theory approaches finite-temperature MBPT of the same order, but it outperforms the latter at intermediate temperature by including additional electron-correlation effects through orbital energies. A physical meaning of these thermal orbital energies (including that of thermal HF orbital energies, which has been elusive) is proposed.

In the past decade, it has been demonstrated that monolayers of metal dichalcogenides are well-suited for thermoelectric applications. ZrX2N4 (X = Si, Ge) is a reasonable choice for thermoelectric applications when considering a favorable value of the figure of merit in two-dimensional (2D) layered materials. In this study, we examined the thermoelectric characteristics of the two-dimensional monolayer of ZrX2N4 (where X can be either Si or Ge) using a combination of Density Functional Theory (DFT) and the Boltzmann Transport Equation (BTE). A thermoelectric figure of merit (ZT) of 0.90 was achieved at a temperature of 900 K for p-type ZrGe2N4, while a ZT of 0.83 was reported for n-type ZrGe2N4 at the same temperature. In addition, the ZrGe2N4 material exhibited a thermoelectric figure of merit (ZT) of around 0.7 at room temperature for the p-type. Conversely, the ZrSi2N4 exhibited a relatively lower thermoelectric figure of merit (ZT) at ambient temperature. At higher temperatures, the ZT value experiences a substantial increase, reaching 0.89 and 0.82 for p-type and n-type materials, respectively, at 900 K. Through our analysis of the electronic band structure, we have determined that ZrSi2N4 and ZrGe2N4 exhibit indirect bandgaps (BG) of 2.74 eV and 2.66 eV, respectively, as per the Heyd-Scuseria-Ernzerhof (HSE) approximation.

The steady-state current-voltage response of ion-selective systems varies as the number of ion-selective components is varied. For the highly investigated unipolar system, including only one ion-selective component, it has been shown that above a supercritical voltage, an electroosmotic instability is triggered, leading to over-limiting currents. In contrast, the effects of this instability on the current-voltage response of the second most common system of a bipolar system, including two oppositely charged permselective regions, have yet to be reported. Using simulations, we investigate the steady-state electrical response of bipolar systems as we vary the ratio of the charge within the two oppositely charged regions. The responses are divided into those with an internal symmetry related to the surface charge and those without. In contrast to the unipolar systems, bipolar systems with the internal symmetry do not exhibit overlimiting currents, and their steady-state response is identical to the convectionless steady-state response. In contrast, the systems without the internal symmetry exhibit much more complicated behavior. For positive voltages, they have overlimiting currents, while for negative voltages, they do not have over-limiting currents. Our findings contribute to a more profound understanding of the behavior of the current-voltage response in bipolar systems.

The most common method to characterize the electrical response of a nanofluidic system is through its steady-state current-voltage response. In Part I, we demonstrated that this current-voltage response depends on the geometry, the layout of the surface charge, and the effects of advection. We demonstrated that each configuration has a unique steady-state signature. Here, we will elucidate the behavior of the time-transient response. Similar to the steady-state response, we will show that each configuration has its own unique time-transient signature when subjected to electroosmotic instability. We show that bipolar systems behave differently than unipolar systems. In unipolar systems, the instability appears only at one end of the system. In contrast, in bipolar systems the instability will either appear on both sides of the nanochannel or not at all. If it does appear on both sides, the instability will eventually vanish on one or both sides of the system. In Part I, these phenomena were explained using steady-state considerations of the behavior of the fluxes. Here, we will examine the time-transient behavior to reveal the governing principles that are, on the one hand, not so different from unipolar systems and, on the other hand, remarkably different.

Incompressible insulating phases of electronic systems at partial filling of a lattice are often associated with charge ordering that breaks lattice symmetry. The resulting phases have an enlarged unit cell with an effective integer filling. Here we explore the possibility of insulating states - which we dub "Quantum Charge Liquids" (QCL) - at partial lattice filling that preserve lattice translation symmetry. Such QCL phases must necessarily either have gapped fractionally charged excitations and associated topological order or have gapless neutral excitations. We establish some general constraints on gapped fermionic QCL phases that restrict the nature of their topological order. We prove a number of results on the minimal topological order that is consistent with the lattice filling. In particular we show that at rational fillings $\nu = p/q$ with $q$ an even integer the minimal ground state degeneracy on a torus of the fermionic QCL is $4q^2$, 4 times larger than that of the bosonic QCL at the same filling. We comment on models and physical systems which may host fermionic QCL phases and discuss the phenomenology of these phases.

We consider chiral, generally nonlinear density waves in one dimension, modelling the bosonized edge modes of a two-dimensional fermionic topological insulator. Using the coincidence between bosonization and Lie-Poisson dynamics on an affine U(1) group, we show that wave profiles which are periodic in time produce Berry phases accumulated by the underlying fermionic field. These phases can be evaluated in closed form for any Hamiltonian, and they serve as a diagnostic of nonlinearity. As an explicit example, we discuss the Korteweg-de Vries equation, viewed as a model of nonlinear quantum Hall edge modes.

Superconductivity in the recently discovered altermagnetic materials hosts large prospects for both fundamental physics and technological applications. In this work we show that a characteristic spin-sublattice locking in altermagnets puts severe constraints on possible superconducting pairing. In particular, we uncover that the most common form of superconductivity, uniform $s$-wave spin-singlet pairing is not possible to achieve in altermagnets. Considering an effective model for a $d_{x^2-y^2}$-wave altermagnet on a square lattice, we instead find that the most likely forms of spin-singlet pairing have $d_{x^2-y^2}$- or extended $s$-wave symmetry. We also find that the simplest form of equal-spin-triplet $p$-wave pairing is not allowed, but it can only exist as a mixed-spin-triplet $p$-wave state. We verify these constraints on pairing within an interaction-induced model of altermagnetism, where we also establish their validity for finite-momentum pairing. Additionally we discuss the possible pairing symmetries for odd-frequency superconducting pairing. Due to the generality of our results, they are applicable to both intrinsic superconductivity and proximity-induced superconductivity in altermagnet-superconductor hybrid junctions.

A quantum spin liquid (QSL) is an exotic insulating phase with emergent gauge fields and fractionalized excitations. However, the unambiguous demonstration of the existence of a QSL in a "non-engineered" microscopic model (or in any material) remains challenging. Here, using numerically-exact sign-problem-free quantum Monte Carlo simulations, we show that a QSL arises in a non-engineered electron-phonon model. Specifically, we investigate the ground-state phase diagram of the bond Su-Schrieffer-Heeger (SSH) model on a 2D triangular lattice at half filling (one electron per site) which we show includes a QSL phase which is fully gapped, exhibits no symmetry-breaking order, and supports deconfined fractionalized holon excitations. This suggests new routes for finding QSLs in realistic materials and high-$T_c$ superconductivity by lightly doping them.

Fusion surface models, as introduced by Inamura and Ohmori, extend the concept of anyon chains to 2+1 dimensions, taking fusion 2-categories as their input. In this work, we construct and analyze fusion surface models on the honeycomb lattice built from braided fusion 1-categories. These models preserve mutually commuting plaquette operators and anomalous 1-form symmetries. Their Hamiltonian is chosen to mimic the structure of Kitaev's honeycomb model, which is unitarily equivalent to the Ising fusion surface model. In the anisotropic limit, where one coupling constant is dominant, the fusion surface models reduce to Levin-Wen string-nets. In the isotropic limit, they are described by weakly coupled anyon chains and are likely to realize chiral topological order. We focus on three specific examples: (i) Kitaev's honeycomb model with a perturbation breaking time-reversal symmetry that realizes chiral Ising topological order, (ii) a $\mathbb{Z}_N$ generalization proposed by Barkeshli et al., which potentially realizes chiral parafermion topological order, and (iii) a novel Fibonacci honeycomb model featuring a non-invertible 1-form symmetry.

Easy-plane quantum magnets are strikingly similar to superconductors, allowing for spin supercurrent and an effective superconducting phase stemming from their $U(1)$ rotation symmetry around the $z$-axis. We uncover a generalized fractional Josephson effect with a periodicity that increases linearly with system size in one-dimensional spin-$1/2$ chains at selected anisotropies and phase-fixing boundary fields. The effect combines arbitrary integer periodicities in a single system, exceeding the $4\pi $ and $8\pi $ periodicity of superconducting Josephson effects of Majorana zero modes and other exotic quasiparticles. We reveal a universal energy-phase relation and connect the effect to the recently discovered phantom helices.

Following the recent interest in van der Waals (vdW) ferroelectrics, topologically nontrivial polar structures have been predicted to form in twisted bilayers. Due to the unconventional nature of vdW ferroelectricity, these topological polar structures have proven difficult to observe experimentally. Here we propose that these textures may be probed optically by showing that topological polarization textures result in exotic nonlinear optical responses. We derive this relationship analytically using non-Abelian Berry connections and a quantum-geometric framework, supported by tight-binding and first-principles calculations. For the case of moir\'e materials without centrosymmetry, which form networks of polar merons and antimerons, the shift photoconductivity forms a vortex-like structure in real space. For a range of frequencies where transitions occur at the Brillouin zone edge, the shift photocurrents are anti-parallel to the in-plane electronic polarization field. Our findings highlight the interplay between complex polarization textures and nonlinear optical responses in vdW materials and provide a sought-after strategy for the experimental detection of topological polarization structures.