Search

In disordered, viscoelastic or viscoplastic materials a sample response exhibits a recovery phenomenon after the removal of a constant load or after creep. We study experimentally the recovery in paper, a quasi two-dimensional system with intrinsic structural disorder. The deformation is measured by using the digital image correlation (DIC) method. By the DIC we obtain accurate displacement data and the spatial fields of deformation and recovered strains. The averaged results are first compared to several heuristic models for in particular viscoelastic polymer materials. The most important experimental quantity is the permanent creep strain, and we analyze whether it is non-zero by fitting the empirical models of viscoelasticity. We then present in more detail the spatial recovery behavior results from DIC, and show that they indicate a power-law -type relaxation. We outline results on sample-to-sample variation and collective, spatial fluctuations in the recovery behaviour. An interpretation is provided of the relaxation in the general context of glassy, interacting systems with barriers.

We present statistical analysis of acoustic emission (AE) data from tensile experiments on paper sheets, loading mode I, with samples broken under strain control. The results are based on 100 experiments on unnotched samples and 70 samples with a long initial edge notch. First, AE energy release and AE event rates are considered for both cases, to test for the presence of "critical points" in fracture. For AE energy, no clear signatures are found, whereas the main finding is that the event rate diverges when a sample-dependent "critical time" of the maximum event rate is approached. This takes place after the maximum stress is reached. The results are compared with statistical fracture models of heterogenous materials. We also discuss the dependence of the AE energy and event interval distributions on average event rates., Comment: 25 pages, 15 figures

The dynamics of a "peeling front" or an elastic line is studied under creep (constant load) conditions. Our experiments show an exponential dependence of the creep velocity on the inverse force (mass) applied. In particular, the dynamical correlations of the avalanche activity are discussed here. We compare various avalanche statistics to those of a line depinning model with non-local elasticity, and study various measures of the experimental avalanche-avalanche and temporal correlations such as the autocorrelation function of the released energy and aftershock activity. From all these we conclude, that internal avalanche dynamics seems to follow "line depinning" -like behavior, in rough agreement with the depinning model. Meanwhile, the correlations reveal subtle complications not implied by depinning theory. Moreover, we also show how these results can be understood from a geophysical point of view., Comment: 22 pages, 14 figs

Decades ago, Gurzhi proposed that if momentum-conserving collisions prevail among heat-carrying phonons in insulators and charge-carrying electrons in metals, hydrodynamic features will become detectable. In this paper, we will review the experimental evidence emerging in the last few years supporting this viewpoint and raising new questions. The focus of the paper will be bulk crystals without (or with a very dilute concentration) of mobile electrons and steady-state thermal transport. We will also discuss the possible link between this field of investigation and other phenomena, such as the hybridization of phonon modes and the phonon thermal Hall effect., Comment: A contribution to the special issue of Fizika Nizkikh Temperatur commemorating Radii Gurzhi; 10 pages, 6 figures

Digital image correlation is a widely used technique in the field of experimental mechanics. In fracture mechanics, determining the precise location of the crack tip is crucial. In this paper, we introduce a universal crack tip detection algorithm based on displacement and strain fields obtained by digital image correlation. Iterative crack tip correction formulas are discovered by applying deep symbolic regression guided by physical unit constraints to a dataset of simulated cracks under mode I, II and mixed-mode conditions with variable T-stress. For the training dataset, we fit the Williams series expansion with super-singular terms to the simulated displacement fields at randomly chosen origins around the actual crack tip. We analyse the discovered formulas and apply the most promising one to digital image correlation data obtained from uniaxial and biaxial fatigue crack growth experiments of AA2024-T3 sheet material. Throughout the experiments, the crack tip positions are reliably detected leading to improved stability of the crack propagation curves.

The boson peak (BP), a low-energy excess in the vibrational density of states over the phonon Debye contribution, is usually identified as one of the distinguishing features between ordered crystals and amorphous solid materials. Despite decades of efforts, its microscopic origin still remains a mystery and a consensus on its theoretical derivation has not yet been achieved. Recently, it has been proposed, and corroborated with simulations, that the BP might stem from intrinsic localized modes which involve string-like excitations ("stringlets") having a one-dimensional (1D) nature. In this work, we build on a theoretical framework originally proposed by Lund that describes the localized modes as 1D vibrating strings, but we specify the stringlet size distribution to be exponential, as observed in independent simulation studies. We show that a generalization of this framework provides an analytically prediction for the BP frequency $\omega_{BP}$ in the temperature regime well below the glass transition temperature in both 2D and 3D amorphous systems. The final result involves no free parameters and is in quantitative agreement with prior simulation observations. Additionally, this stringlet theory of the BP naturally reproduces the softening of the BP frequency upon heating and offers an analytical explanation for the experimentally observed scaling with the shear modulus in the glass state and changes in this scaling in cooled liquids. Finally, the theoretical analysis highlights the existence of a strong damping for the stringlet modes at finite temperature which leads to a large low-frequency contribution to the 3D vibrational density of states, as observed in both experiments and simulations., Comment: We withdraw the manuscript temporally to abide by the internal review requirements of one of the authors (Jack Douglas) at the National Institute of Standards and Technology. (NIST). The withdrawal has nothing to do with the content of the paper, and we plan to resubmit our paper to the arXiv once the paper has been approved for public release by the NIST internal review board

During his time at Royal Institute of Technology (Kungliga Tekniska Hogskolan, KTH) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concept is the separation of internal and external variables of a system. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book "Phase Equilibria, Phase Diagrams and Phase Transformations" that puts equilibrium, nonequilibrium, and internal processes on equal footing. In the present paper honoring Hillert, the present author reflects his experiences with Hillert in last 35 years and expresses his gratitude for all the wisdom and support from him in terms of "Hillert nonequilibrium thermodynamics" and discusses some recent topics that the present author has been working on.

Molecular dynamics (MD) simulations play an important role in understanding and engineering heat transport properties of complex materials. An essential requirement for reliably predicting heat transport properties is the use of accurate and efficient interatomic potentials. Recently, machine-learned potentials (MLPs) have shown great promise in providing the required accuracy for a broad range of materials. In this mini review and tutorial, we delve into the fundamentals of heat transport, explore pertinent MD simulation methods, and survey the applications of MLPs in MD simulations of heat transport. Furthermore, we provide a step-by-step tutorial on developing MLPs for highly efficient and predictive heat transport simulations, utilizing the neuroevolution potentials (NEPs) as implemented in the GPUMD package. Our aim with this mini review and tutorial is to empower researchers with valuable insights into cutting-edge methodologies that can significantly enhance the accuracy and efficiency of MD simulations for heat transport studies., Comment: 25 pages, 9 figures. This paper is part of the special topic, Machine Learning for Thermal Transport

Macroscopic mechanical properties of polymers are determined by their microscopic molecular chain distribution. Due to randomness of these molecular chains, probability theory has been used to find their micro-states and energy distribution. In this paper, aided by central limit theorem and mixed Bayes rule, we showed that entropy elasticity based on Gaussian distribution is questionable. By releasing freely jointed chain assumption, we found that there is energy redistribution when each bond of a molecular chain changes its length. Therefore, we have to change Gaussian distribution used in polymer elasticity to Maxwell-Boltzmann distribution. Since Maxwell-Boltzmann distribution is only a good energy description for gas molecules, we found a mathematical path to change Maxwell-Boltzmann distribution to Fermi-Dirac distribution based on molecular chain structures. Because a molecular chain can be viewed as many monomers glued by covalent electrons, Fermi-Dirac distribution describes the probability of covalent electron occupancy in micro-states for solids such as polymers. Mathematical form of Fermi-Dirac distribution is logistic function. Mathematical simplicity and beauty of Fermi-Dirac distribution make many hard mechanics problems easy to understand. Generalized logistic function or Fermi-Dirac distribution function was able to understand many polymer mechanics problems such as viscoelasticity [1], viscoplasticity [2], shear band and necking [3], and ultrasonic bonding [4]., Comment: 3 figures

Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.

The existence of a crumpled phase for self-avoiding elastic surfaces was postulated more than three decades ago using simple Flory-like scaling arguments. Despite much effort, its stability in a microscopic environment has been the subject of much debate. In this paper we show how a crumpled phase develops reliably and consistently upon subjecting a thin spherical shell to active fluctuations. We find a master curve describing how the relative volume of a shell changes with the strength of the active forces, that applies for every shell independent of size and elastic constants. Furthermore, we extract a general expression for the onset active force beyond which a shell begins to crumple. Finally, we calculate how the size exponent varies along the crumpling curve., Comment: 9 pages and 9 figures including the appendix

In recent decades, computer experiments have led to an accurate and fundamental understanding of atomic and molecular mechanisms in fluids, such as different kinds of relaxation processes toward steady physical states. In this paper, we investigate how exactly the configuration of initial states in a molecular-dynamics simulation can affect the rates of decay toward equilibrium for the widely-known Langevin canonical ensemble. For this purpose, we derive an original expression relating the system relaxation time {\tau}_{sys} and the radial distribution function g(r) in the near-zero and high-density limit. We found that for an initial state which is slightly marginally inhomogeneous in the number density of atoms, the system relaxation time {\tau}_{sys} is much longer than that for the homogeneous case and an increasing function of the Langevin coupling constant, {\gamma}. We also found during structural equilibration, g(r) at large distances approaches 1 from above for the inhomogeneous case and from below for the macroscopically homogeneous one.

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon ($a$-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in $a$-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts -- originating from isolated and vacancy-induced clustered defects -- and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct $ab$ $initio$ molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in $a$-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of $a$-Si is found to be in the range of 488--541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487--528 K obtained from specific-heat measurements of $a$-Si at low temperatures., Comment: 14 pages, 17 figures

A simple stochastic model of solute drag by moving grain boundaries (GBs) is presented. Using a small number of parameters, the model describes solute interactions with GBs and captures nonlinear GB dynamics, solute saturation in the segregation atmosphere, and the breakaway from the atmosphere. The model is solved by kinetic Monte-Carlo (KMC) simulations with time-dependent transition barriers. The non-Markovian nature of the KMC process is discussed. In Part I of this work, the model is applied to planar GBs driven by an external force. The model reproduces all basic features of the solute drag effect, including the maximum of the drag force at a critical GB velocity. The force-velocity functions obtained depart from the scaling predicted by the classical models by Cahn and L\"ucke-St\"uwe, which are based on more restrictive assumptions. The paper sets the stage for Part II, in which the GB will be treated as a 2D solid-on-solid interface., Comment: Submitted to Physical Review Materials

In a recent paper [S. Khrapak, Molecules {\bf 25}, 3498 (2000)] the longitudinal and transverse sound velocities of conventional Lennard-Jones systems at the liquid-solid coexistence were calculated. It was shown that the sound velocities remain almost invariant along the liquid-solid coexistence boundary lines and that their magnitudes are comparable with those of repulsive soft sphere and hard sphere models at the fluid-solid phase transition. This implies that attraction does not affect the magnitude of sound velocities at the fluid-solid phase transition. This paper provides further evidence to this by examining the generalized Lennard-Jones $n$-6 fluids with $n$ ranging from $12$ to $7$ and demonstrating that the steepness of repulsive term has only a minor effect on the magnitude of the sound velocities., Comment: 4 pages, 1 figure

Granular materials such as sand, powders, and food grains are ubiquitous in civil engineering, geoscience, agriculture, and medicine. While the influence of friction between the grains on the static structure of these systems is well understood, its impact on the dynamics is an open problem. Here we use particle-based simulations of a granular pack under cyclic shear and discover that the relaxation time of the system is a non-monotonic function of friction. By introducing the concept of dynamic constraints, we reveal that this re-entrant dynamics is due to the competition between increasing frictional coupling and a concurrent change in the structure of the granular pack. Our theoretical approach, which unifies the dynamics of friction-less systems with frictional ones, is applicable to other systems that have a complex free energy landscape and a dynamics which involves time-dependent constraints, thus setting the stage for a description of the dynamic behavior of a large class of complex systems.

Topological defects are singularities in the order parameter space that are mathematically described by topological invariants and cannot be removed by continuous transformations. These defects play a significant role in various fields, ranging from cosmology to solid-state physics and biological matter. The definition of these irregularities requires an ordered reference configuration, leading to decades of debate about the existence of topological defects in disordered systems, such as glasses. Recently, it has been proposed that well-defined topological defects might emerge in the dynamical properties of glasses under deformation, potentially relating to their plastic behavior. In this study, we investigate a two-dimensional colloidal glass system composed of particles interacting by an effective magnetic repulsion. We reveal the presence of topological defects in the eigenspace of the vibrational frequencies of this experimental 2D amorphous solid. The vibrational density of states of this 2D amorphous material exhibits distinct glassy properties such as the presence of a boson peak anomaly. Through extensive numerical and theoretical quantitative analysis, we establish a robust positive correlation between the vibrational characteristics and the total number of topological defects. We furthermore show that defects of opposite charge tend to pair together and prove their local statistical correlation with the "soft spots", the regions more prone to plastic flow. This work provides the experimental confirmation for the existence of topological defects in disordered systems revealing a complex interplay between topology, disorder, and vibrational behavior.

Recoverable strain is the strain recovered once a stress is removed from a body, in the direction opposite to that in which the stress had acted. To date, the phenomenon has been understood as being elastic in origin: polymer chains stretched in the direction of an imposed stress will recoil after the stress is removed, for example. Any unrecoverable strain is instead attributed to irreversible plastic deformations. Here we study theoretically recoverable and unrecoverable strain within the soft glassy rheology model, aimed at describing the rheology of elastoplastic yield stress fluids and amorphous soft solids. We consider a material subject to the switch-on of a shear stress that is held constant before later being set back to zero, after which the strain recovery is observed. After an initially fast recoil that is indeed elastic in nature, consistent with the usual intuition, we find that significant subsequent strain recovery then arises not via further recoil of elasticity but instead by `reverse plasticity', in a sense that we discuss carefully. Also unexpectedly, although in rare parameter regimes, this plastic part of the strain post switch-off does not always in fact recover in the negative direction, counter to that of the previously imposed stress, but can sometimes continue to accumulate in the forward direction. The recovery is then non-monotonic overall, reminiscent of recent observations of non-monotonic stress relaxation after straining, which have been attributed to complex material memory effects. We elucidate the mechanisms underlying these behaviours in terms of the evolution of the SGR model's population of elastoplastic elements, and suggest connections with the notion of `reversible plasticity', as discussed in the recent physics literature., Comment: submitted for publication, 11 pages, 6 figures

Traditionally, stress fluctuations in flowing and deformed materials are overlooked, with an obvious focus on average stresses in a continuum mechanical approximation. However, these fluctuations, often dismissed as noise, hold the potential to provide direct insights into the material structure and its structure-stress coupling, uncovering detailed aspects of fluid transport and relaxation behaviors. Despite advancements in experimental techniques allowing for the visualization of these fluctuations, their significance remains largely untapped, as modeling efforts continue to target Newtonian fluids within the confines of Gaussian noise assumptions. In the present work a comparative analysis of stress fluctuations in two distinct microstructural models is carried out: the temporary network model and the dumbbell model. Despite both models conforming to the Upper Convected Maxwell Model at a macroscopic level, the temporary network model predicts non-Gaussian fluctuations. We find that stress fluctuations within the temporary network model exhibit more pronounced abruptness at local scale, with only an enlargement of the control volume leading to a gradual Gaussian-like noise, diminishing the differences between the two models. These findings underscore the heightened sensitivity of fluctuating rheology to microstructural details and the microstructure-flow coupling, beyond what is captured by macroscopically averaged stresses., Comment: 11 pages, 12 figures, The following article has been submitted to The Journal of Chemical Physics. After it is published, it will be found at Link: https://pubs.aip.org/aip/jcp?gad_source=1&gclid=Cj0KCQjwir2xBhC_ARIsAMTXk85XWFzcccAjiPwwNGQTvsW5B-DrpAbdS5nJ0D75RtXKOqHGsfyNzpgaAqVhEALw_wcB

We numerically study the relaxation dynamics of impurity-host systems, focusing on the presence of long-lived metastable states in the non-equilibrium dynamics after an initial excitation of the impurities. In generic systems, an excited impurity coupled to a large bath at zero temperature is expected to relax and approach its ground state over time. However, certain exceptional cases exhibit metastability, where the system remains in an excited state on timescales largely exceeding the typical relaxation time. We study this phenomenon for three prototypical impurity models: a tight-binding quantum model of independent spinless fermions on a lattice with two stub impurities, a classical-spin Heisenberg model with two weakly coupled classical impurity spins, and a tight-binding quantum model of independent electrons with two classical impurity spins. Through numerical integration of the fundamental equations of motion, we find that all three models exhibit similar qualitative behavior: complete relaxation for nearest-neighbor impurities and incomplete or strongly delayed relaxation for next-nearest-neighbor impurities. The underlying mechanisms leading to this behavior differ between models and include impurity-induced bound states, emergent approximately conserved local observables, and exact cancellation of local and nonlocal dissipation effects., Comment: 15 pages, 14 figures

By critically examining the traditional theory of homogeneous nucleation of precipitates in solid solutions, it is revealed that the theory's assertion regarding an increase in the nucleation free energy due to elastic strain associated with the difference in atomic volumes between the two phases is applicable to coherent precipitates, but becomes incorrect when applied to incoherent precipitates. This conclusion is obtained by accounting for thermal point defects in the matrix, which can be absorbed at the interface between an incoherent particle and the matrix during nucleation, thereby relieving elastic stresses. Accordingly, a new kinetic model based on the Reiss theory for binary nucleation is proposed for predicting the nucleation rate of incoherent precipitates by agglomeration of solute atoms and point defects, with a further extension to account for excess vacancies formed under non-equilibrium conditions of quenching experiments., Comment: Published in Russian Zhurnal Eksperimental'noi i Teoreticheskoi Fiziki 165 (2024) 214-225

The ability to understand and ultimately control the transformations and properties of various nanoscale systems, from proteins to synthetic nanomaterial assemblies, hinges on the ability to directly elucidate their dynamics on their characteristic length and time scales. Here, we use MHz X-ray photon correlation spectroscopy (XPCS) to directly elucidate the characteristic microsecond-dynamics of density fluctuations of semiconductor nanocrystals (NCs), not only in a colloidal dispersion but also in a liquid phase consisting of densely packed, yet mobile, NCs with no long-range order. By carefully disentangling X-ray induced effects, we find the wavevector-dependent fluctuation rates in the liquid phase are suppressed relative to those in the colloidal phase and to those in experiments and hydrodynamic theories of densely packed repulsive particles. We show that the suppressed rates are due to a substantial decrease in the self-diffusion of NCs in the liquid phase, which we attribute to explicit attractive interactions. Via comparison with simulations, we find that the extracted strength of the attractions explains the stability of the liquid phase, in contrast to the gelation observed via XPCS in many other charged colloidal systems. This work opens the door to elucidating fast, condensed phase dynamics in a variety of complex fluids and other nanoscale soft matter systems, such as densely packed proteins and non-equilibrium self-assembly processes., Comment: 15 pages, 4 figures

Compared to the widely investigated crystalline polymorphs of gallium oxide (Ga2O3), knowledge about its amorphous state is still limited. With the help of a machine-learning interatomic potential, we conducted large-scale atomistic simulations to investigate the glass transition and mechanical behavior of amorphous Ga2O3 (a-Ga2O3). During the quenching simulations, amorphization of gallium oxide melt is observed at ultrahigh cooling rates, including a distinct glass transition. The final densities at room temperature have up to 4% variance compared to experiments. The glass transition temperature is evaluated to range from 1234 K to 1348 K at different cooling rates. Structural analysis of the amorphous structure shows evident similarities in structural properties between a-Ga2O3 and amorphous alumina (a-Al2O3), such as radial distribution function, coordination distribution, and bond angle distribution. An amorphous gallium oxide structure that contains approximately one million atoms is prepared for the tension simulation. A highly plastic behavior is observed at room temperature in the tension simulations, comparable to amorphous alumina. With quantitative characterization methods, we show that a-Ga2O3 can possibly has a higher nucleation rate of localized plastic strain events compared to a-Al2O3, which can increase the material's resistance to shear banding formation during deformation., Comment: 15 pages, 7 figures, under review

The probabilistic cellular automaton (PCA) method is highlighted for its relatively simple numerical algorithm and low computational cost in the simulation of microstructural evolution. In this method, probabilistic state change rules are implemented to compute the evolution of cell states at each time step. The stochastic nature of this simulation method leads to non-repeatable simulation results, introducing inherent uncertainty. In this study, the uncertainty and dispersion in PCA simulations of microstructural evolution were investigated. Hence, the probabilistic transformations of cell states were meticulously considered at each time step, and discrete probability distribution functions (dPDF) were introduced to analyze the frequency distribution of simulation outcomes. To evaluate the performance of the proposed dPDFs, cellular automaton models were developed with various numbers of cells and distribution of transformation probabilities. Multiple iterations of these simulations were conducted, and the validity of the presented distribution functions was assessed through statistical analysis of the simulations outcomes. Comparisons between PCA simulation results and distribution functions demonstrate consistency, emphasizing the predictive capability of the proposed models. Also, the effects of modeling parameters on the uncertainty of simulation results in two and three-dimensional PCA modeling were studied, introducing the coefficient of variation as a measure of dispersion. Results indicate that increasing the number of boundary cells, cellular resolution, and model size reduces uncertainty, enhancing the repeatability of PCA simulation outcomes.

Metal halide perovskites have emerged as an exciting class of materials for applications in solar energy harvesting, optical devices, catalysis, and other photophysical applications. Many of the exciting properties of halide perovskites are tied to their soft, dynamic, and anharmonic lattice. In particular, the precise coupling between anharmonic lattice dynamics and electronic fluctuations is not completely understood. To build an understanding of this coupling, we use ab initio molecular dynamics simulations supplemented by the calculation of maximally localized Wannier functions to carry out a dynamic group theory analysis of local electron density fluctuations and how these fluctuations are coupled to lattice fluctuations in the model inorganic halide perovskite CsSnBr3. We detail symmetry-dependent couplings between vibrational modes, including octahedral tilting. Importantly, we suggest that the large anharmonicity of some of the vibrational modes in CsSnBr3 result from electron rotation--nuclear translation coupling, in analogy to rotation--translation coupling effects in molecular plastic crystals. We also identify electronic fluctuations in the Cs cation that couple to distortions in the surrounding Sn-Br cubic coordination environment. We anticipate that our approach and resulting insights into electronic fluctuations will aid in further understanding the role of the fluctuating lattice in determining important physical properties of halide perovskites and beyond., Comment: 15 pages, 15 figures

Modeling ion correlations in inhomogeneous polymers and soft matters with spatially varying ionic strength or dielectric permittivity remains a great challenge. Here, we develop a new theory which systematically incorporates electrostatic fluctuations into the self-consistent field theory for polymers. Applied to polyelectrolyte brushes, the theory predicts that ion correlations induce non-monotonic change of the brush height: collapse followed by reexpansion. The scaling analysis elucidates the competition between the repulsive osmotic pressure due to translational entropy and the attraction induced by ion correlations. We also clarify the absence of causal relationship between the brush collapse-reexpansion and the inversion of the surface electrostatic potential. Furthermore, strong ion correlations can trigger microphase separation, either in the lateral direction as pinned micelles or in the normal direction as oscillatory layers. Our theoretical predictions are in good agreement with the experimental results reported in the literature.

The structural, dielectric, and thermodynamic properties of the hydrogen-bonded ferroelectric crystal potassium dihydrogen phosphate ($\mathbf{KH_2PO_4}$), KDP for short, differ significantly from those of DKDP ($\mathbf{KD_2PO_4}$). It is well established that deuteration affects the interplay of hydrogen-bond switches and heavy ion displacements that underlie the emergence of macroscopic polarization, but a detailed microscopic model is missing. Here we show that all-atom path integral molecular dynamics simulations can predict the isotope effects, revealing the microscopic mechanism that differentiates KDP and DKDP. Proton tunneling in the hydrogen bonds generates phosphate configurations that do not contribute to the polarization. These dipolar defects are always present in KDP, but disappear at low temperatures in DKDP, which behaves more classically. Quantum disorder confers residual entropy to ferroelectric KDP, explaining its lower spontaneous polarization and transition entropy relative to DKDP. Tunneling should also contribute to the anomalous heat capacity observed in KDP near absolute zero. The prominent role of quantum fluctuations in KDP is related to the unusual strength of the hydrogen bonds in this system and should be equally important in the other crystals of the KDP family, which exhibit similar isotope effects.

The first-order phase transition between the liquid and gaseous phases ends at a critical point. Critical opalescence occurs at this singularity. Discovered in 1822, it is known to be driven by diverging fluctuations in the density. During the past two decades, boundaries between the gas-like and liquid-like regimes have been theoretically proposed and experimentally explored. Here, we show that fast cooling of near-critical sulfur hexafluoride (SF$_6$), in presence of Earth's gravity, favors dark opalescence, where visible photons are not merely scattered, but also absorbed. When the isochore fluid is quenched across the critical point, its optical transmittance drops by more than three orders of magnitude in the whole visible range, a feature which does not occur during slow cooling. We show that transmittance shows a dip at 2eV near the critical point, and the system can host excitons with binding energies ranging from 0.5 to 4 eV. The spinodal decomposition of the liquid-gas mixture, by inducing a periodical modulation of the fluid density, can provide a scenario to explain the emergence of this platform for coupling between light and matter. The possible formation of excitons and polaritons points to the irruption of quantum effects in a quintessentially classical context.

Lattice softness has a significant impact on charge carrier dynamics in condensed matter systems, contributing to the emergence of various properties and functions. Examples include the remarkable carrier lifetimes and defect tolerances of hybrid organic-inorganic perovskites. Recent studies suggest the contribution of quartic anharmonicity of the lattice vibrations. The quartic anharmonicity can be discussed with a double-well potential, and the transition between the two minima can be coarse-grained as a two-state jump stochastic process. Such a stochastic approach is typically employed to describe dynamic fluctuations introduced into a system by two-state transitions in the surroundings. To investigate charge transport in materials, however, it is crucial to describe not only the fluctuations but also the dynamic lattice distortion associated with charge transport. Therefore, there is a need for a theory to describe the charge carrier dynamics proceeding alongside the lattice distortion dynamics. In this study, we present a theory that describes quantum dynamics under the influence of an environment with two stable states, termed a bistable environment. The theory describes the effects of fluctuations and dissipation induced from the bistable environment in a reasonable manner, and the effects exhibit a different temperature dependence than the widely employed Gaussian environment. The physical implication of this temperature dependence is provided in terms of the environmental dynamics. The results of this study are expected to provide a step forward in describing charge carrier dynamics in materials with lattice softness and pronounced lattice anharmonicity, e.g., hybrid organic-inorganic perovskites. Moreover, these findings represent an advancement in our understanding of and capacity to predict and control the physical properties and functions of these materials.

Strain at interfaces may profoundly impact the microstructure and properties of materials; thus, it is a major consideration when designing and engineering materials. Dislocation formation is a commonly known mechanism to release mismatch strain at solid-solid interfaces. However, it is still unclear about how materials accommodate interfacial strain under drastically accelerated structural transformation kinetics, since it is extremely challenging to directly observe the atomic structure evolution of fast-propagating interfaces. Utilizing liquid phase transmission electron microscopy (TEM), we have achieved atomic-scale imaging of hydrogen-induced phase transformations of palladium nanocrystals with different transformation speeds. Our observation reveals that the fast phase transformation occurs with an expanded interface of mixed $\alpha$- and $\beta$-$\mathrm{PdH}_x$ phases, and tilting of (020) planes to accommodate mismatch strain. In contrast, slow phase transformations lead to sharp interfaces with slipping misfit dislocations. Our kinetic Monte Carlo simulations show that fast phase transformation pushes the system far-from-equilibrium, generically roughening the interface; however, a smooth boundary minimizes strain near-equilibrium. Unveiling the atomic pathways of transformations from near-equilibrium to far-from-equilibrium, which was previously possible only computationally, this work holds significant implications for engineering microstructure of materials through modulating solid-solid transformations in a wide range of kinetics., Comment: 6 pages, 4 figures, plus Methods and Supplementary Notes

The Taylor-Aris theory of shear diffusion predicts that the effective diffusivity of a tracer molecule in a sheared liquid is enhanced by a term quadratic in the shear rate. In sheared supercooled liquids, instead, the observed enhancement is linear in the shear rate. This is a fundamental observation for the physics of nonequilibrium liquids for which no theory or fundamental understanding is available. We derive a formula for the effective molecular diffusivity in supercooled liquids under shear flow based on the underlying Smoluchowski equation with shear (Smoluchowski diffusion-convection equation) with an energy barrier representing a crowded energy landscape. The obtained formula correctly recovers the effective diffusivity with a correction term linear in the shear rate, in agreement with results from numerical simulations of different liquids as well as with earlier experimental results on shear melting of colloidal glass. The theory predictions are supported by comparisons with molecular simulations of supercooled water and supercooled Lennard-Jones liquids, which confirm that the predicted enhancement of diffusivity is inversely proportional to temperature and directly proportional to the zero shear viscosity.

Solid materials that deviate from the harmonic crystal paradigm exhibit characteristic anomalies in the specific heat and vibrational density of states (VDOS) with respect to Debye's theory predictions. The boson peak (BP), a low-frequency excess in the VDOS over Debye law $g(\omega) \propto \omega^2$, is certainly the most famous among them; nevertheless, its origin is still subject of fierce debate. Recent simulation works provided strong evidence that localized one-dimensional string-like excitations (stringlets) might be the microscopic origin of the BP. In this work, we study the dynamics of acoustic phonons interacting with a bath of vibrating 1D stringlets with exponentially distributed size, as observed in simulations. We show that stringlets strongly renormalize the phonon propagator and naturally induce a BP anomaly in the vibrational density of states, corresponding to the emergence of a dispersionless BP flat mode. Additionally, phonon-stringlet interactions produce a strong enhancement of sound attenuation and a dip in the speed of sound near the BP frequency, consistent with experimental and simulation data. The qualitative trends of the BP frequency and intensity are predicted within the model and shown to be in good agreement with previous observations. In summary, our results substantiate with a simple theoretical model the recent simulation results by Hu and Tanaka claiming the origin of the BP from stringlet dynamics., Comment: v1

I define the fields that describe electrical macrostructure, and their rates of change, in terms of the microscopic charge density, electric field, electric potential, and their rates of change. To deduce these definitions, I lay some foundations of a general theory of structure homogenization, meaning a theory of how any observable macroscopic field is related to spatial averages of its microscopic counterpart. I apply homogenization theory to the fields that appear in the vacuum theory of electromagnetism and deduce that the relationships between macroscopic fields are identical in form to the relationships between their microscopic counterparts. Without invoking quantum mechanics, I derive the expressions for polarization current established by the Modern Theory of Polarization. I prove that the bulk-average electric potential, or mean inner potential, vanishes in a macroscopically-uniform charge-neutral material, and I show that when a crystal lattice lacks inversion symmetry, it does not imply the existence of macroscopic $\vec{E}$ or $\vec{P}$ fields in the crystal's bulk. I point out that symmetry is scale-dependent. Therefore, if anisotropy of the microstructure does not manifest as anisotropy of the macrostructure, it cannot be the origin of a macroscopic vector field. The macroscopic charge density is isotropic in the bulks of most materials, because it vanishes at every point. Therefore, regardless of the microstructure, a macroscopic electric field cannot emanate from the bulk. I find that all relationships between observable macroscopic fields can be expressed mathematically without introducing the polarization ($\vec{P}$) and electric displacement ($\vec{D}$) fields, neither of which is observable.

Quantum integrability emerging near a quantum critical point (QCP) is manifested by exotic excitation spectrum that is organized by the associated algebraic structure. A well known example is the emergent $E_8$ integrability near the QCP of a transverse field Ising chain (TFIC), which was long predicted theoretically and initially proposed to be realized in the quasi-one-dimensional (q1D) quantum magnet CoNb$_2$O$_6$. However, later measurements on the spin excitation spectrum of this material revealed a series of satellite peaks that cannot be described by the $E_8$ Lie algebra. Motivated by these experimental progresses, we hereby revisit the spin excitations of CoNb$_2$O$_6$ by combining numerical calculation and analytical analysis. We show that, as effects of strong interchain fluctuations, the spectrum of the system near the 1D QCP is characterized by the $D_{8}^{(1)}$ Lie algebra with robust topological soliton excitation. We further show that the $D_{8}^{(1)}$ spectrum can be realized in a broad class of interacting quantum systems. Our results advance the exploration of integrability and manipulation of topological excitations in quantum critical systems., Comment: 6 pages, 3 figures - Supplementary Material 5 pages

The structural origin of the first sharp diffraction peak (FSDP) in amorphous silica is studied by analyzing chemical and radial ordering of silicon (Si) and oxygen (O) atoms in binary amorphous networks. The study shows that the chemical order involving Si--O and O--O pairs play a major role in the formation of the FSDP in amorphous silica. This is supplemented by small contributions arising from the relatively weak Si--Si correlations in the Fourier space. A shell-by-shell analysis of the radial correlations between Si--Si, Si--O and O--O atoms in the network reveals that the position and the intensity of the FSDP are largely determined by atomic pair correlations originating from the first two/three radial shells on a length scale of about 5--8 {\AA}, whereas the fine structure of the intensity curve in the vicinity of the FSDP is perturbatively modified by atomic correlations arising from the radial shells beyond 8 {\AA}. The study leads to a simple mathematical relationship between the position of the radial peaks ($r_k$) in the partial pair-correlation functions and the diffraction peaks ($Q_k$) that can be used to obtain approximate positions of the FSDP and the principal peak. The results are complemented by numerical calculations and an accurate semi-analytical expression for the diffraction intensity obtained from the partial pair-correlation functions of amorphous silica for a given radial shell., Comment: 13 pages, 16 figures

We demonstrate an efficient and accurate, general-purpose first-principles blueprint for calculating anharmonic vibrational free energy and predicting structural phase transition temperatures of solids. Thermodynamic integration is performed without molecular dynamics using only interatomic force constants to model analogues of the true potential and generate their thermal ensembles. By replacing \textit{ab initio} molecular dynamics (AIMD) with statistical sampling of ensemble configurations and trading density-functional theory (DFT) energy calculations on each configuration for a set of matrix operations, our approach enables a faster thermodynamic integration by 4 orders of magnitude over the traditional route via AIMD. Experimental phase transition temperatures of a variety of strongly anharmonic materials with dynamical instabilities including shape-memory alloys are recovered to largely within 25% error. Such a combination of speed and accuracy enables the method to be deployed at a large-scale for predictive mapping of phase transition temperatures.

In this work, we compare existing non-polarizable force fields developed to study the solid or solution phases of hybrid organic-inorganic halide perovskites with the AMOEBA polarizable force field. The aim is to test whether more computationally expensive polarizable force fields like AMOEBA offer better transferability between solution and solid phases, with the ultimate goal being the study of crystal nucleation, growth and other interfacial phenomena involving these ionic compounds. In the context of hybrid perovskites, AMOEBA force field parameters already exist for several elements in solution and we decided to leave them unchanged and to only parameterize the missing ones (Pb\textsuperscript{2+} and CH\textsubscript{3}NH\textsubscript{3}\textsuperscript{+} ions) in order to maximise transferability and avoid over-fitting to the specific examples studied here. Overall, we find that AMOEBA yields accurate hydration free energies (within 5\%) for typical ionic species while showing the correct ordering of stability for the different crystal polymorphs of CsPbI\textsubscript{3} and CH\textsubscript{3}NH\textsubscript{3}PbI\textsubscript{3}. While the existing parameters do not accurately reproduce all transition temperatures and lattice parameters, AMOEBA offers better transferability between solution and solid states than existing non-polarizable force fields.

We consider second-order phase transitions in which the order parameter is a replicated overlap matrix. We focus on a tricritical point that occurs in a variety of mean-field models and that, more generically, describes higher order liquid-liquid or liquid-glass transitions. We show that the static replicated theory implies slowing down with a logarithmic decay in time. The dynamical equations turn out to be those predicted by schematic Mode Coupling Theory for supercooled viscous liquids at a $A_3$ singularity, where the parameter exponent is $\lambda=1$. We obtain a quantitative expression for the parameter $\mu$ of the logarithmic decay in terms of cumulants of the overlap, which are physically observable in experiments or numerical simulations., Comment: 22 pages, 2 figures

The transition from static to dynamic friction when an elastic body is slid over another is now known to result from the motion of interface rupture fronts. These fronts may be either crack-like or pulse-like, with the latter involving reattachment in the wake of the front. How and why these fronts occur remains a subject of active theoretical and experimental investigation, given its wide ranging implications for a range of problems in tribology. In this work, we investigate this question using an elastic lattice-network representation; bulk and interface bonds are simulated to deform and, in the latter case, break and reform dynamically in response to an applied remote displacement. We find that, contrary to the oft-cited rigid body scenario with Coulomb-type friction laws, the type of rupture front observed depends intimately on the location of the applied boundary condition. Depending on whether the sliding solid is pulled, pushed or sheared -- all equivalent applications in the rigid case -- distinct interface rupture modes can occur. We quantify these rupture modes, evaluate the interface stresses that lead to their formation, and and study their subsequent propagation dynamics. A strong analogy between the sliding friction problem and mode II fracture emerges from our results, with attendant wave speeds ranging from slow to Rayleigh. We discuss how these fronts mediate interface motion and implications for the general transition mechanism from static to dynamic friction., Comment: 35 pages, 15 figures

In geometrically frustrated assemblies local inter-subunit misfits propagate to intra-assembly strain gradients, giving rise to anomalous self-limiting assembly thermodynamics. Here, we use theory and coarse-grained simulation to study a recently developed class of ``curvamer'' particles, flexible shell-like particles that exhibit self-limiting assembly due to the build up of curvature deformation in cohesive stacks. To address a generic, yet poorly understood aspect of frustrated assembly, we introduce a model of curvamer assembly that incorporates both {\it intra-particle} shape deformation as well as compliance of {\it inter-particle} cohesive gaps, an effect we can attribute to a {\it finite range of attraction} between particles. We show that the ratio of intra-particle (bending elasticity) to inter-particle stiffness not only controls the regimes of self-limitation but also the nature of frustration propagation through curvamer stacks. We find a transition from uniformly-bound, curvature-focusing stacks at small size to gap-opened, uniformly curved stacks at large size is controlled by a dimensionless measure of inter- versus intra-curvamer stiffness. The finite range of inter-particle attraction determines range of cohesion in stacks are self-limiting, a prediction which is in strong agreement with numerical studies of our coarse-grained colloidal model. These predictions provide critical guidance for experimental realizations of frustrated particle systems designed to exhibit self-limitation at especially large multi-particle scales., Comment: 23 pages, 16 figures, 7 appendices

A simple and computationally efficient algorithm enables implementing negative temperature values in a spin dynamics simulation. The algorithm uses a Langevin spin dynamics thermostat with a negative damping parameter, enabling the thermalization of an arbitrary interacting spin system to the Gibbs energy distribution with a given negative temperature value. Canonical spin dynamics simulations at a negative temperature are as robust as conventional positive spin temperature simulations, providing a tool for quantitative dynamic studies of the physics of highly excited magnetic states. Two simulation case studies describing spin systems with antiferromagnetic and ferromagnetic ground states are explored. The phase transitions occurring in the negative temperature range do not necessarily exhibit similarities with their positive temperature counterparts. The transition temperatures and the character of spin alignment vary depending on the spatial range and strength of spin-spin interactions.

The spin current model of electric polarization of multiferroics is justified via the quantum hydrodynamic method and the mean-field part of the spin-orbit interaction. The spin current model is applied to derive the electric polarization proportional to the scalar product of the spins of the near by ions, which appears to be caused by the Dzylaoshinskii-Moriya interaction. Symmetric tensor spin structure of the polarization is discussed as well. We start our derivations for the ferromagnetic multiferroic materials and present the further generalization for the antiferromagnetic multiferroic materials. We rederive the operator of the electric dipole moment, which provides the macroscopic polarization obtained via the spin current model. Finally, we use the quantum average of the found electric dipole moment operator in order to derive the polarization evolution equation for the antiferromagnetic multiferroic materials. Possibility of spiral spin structures is analyzed., Comment: 12 pages

Polymers near the glass transition temperature Tg often exhibit a breakdown of time-temperature-superposition (TTS), with chain relaxation times and viscosity exhibiting a weaker temperature dependence than segmental relaxation times. The origin of this onset of thermorheological complexity has remained unsettled and a matter of debate. Here we extend the Heterogeneous Rouse Model (HRM), which generalizes the Rouse model to account for dynamic heterogeneity, to make predictions for the relaxation modulus G(t) and complex modulus G*($\omega$) of unentangled polymers near Tg. The HRM predicts that G(t) and G*($\omega$) exhibit enhanced effective scaling exponents in the Rouse regime in the presence of dynamic heterogeneity, with a more rapid decay from the glassy plateau emerging as the system becomes more dynamically heterogeneous on cooling. This behavior is predicted to emerge from a strand-length dependence of the moment of the segmental mobility distribution probed by chain dynamics. We show that the HRM predictions are in good accord with experimental complex modulus data for polystyrene, poly(methyl methacrylate), and poly(2-vinyl pyridine). The HRM also predicts the onset of distinct temperature dependences among chain scale quantities such as terminal relaxation time and viscosity in our experimental systems, apparently resolving one of the most significant standing objections to a heterogeneity-based origin of TTS-breakdown. The HRM thus provides a generalized theory of the chain-scale linear rheological response of unentangled polymers near Tg, accounting for the origin of TTS-breakdown at a molecular mechanistic level. It also points towards a new strategy of inferring the dynamic heterogeneity of glass-forming polymeric systems based on the temperature-evolution of modified scaling in the Rouse regime.

Heat flux ($J$) generally increases with temperature difference in a material. A differential coefficient of $J$ against temperature ($T$) is called differential thermal conductance ($k$), and an inverse of $k$ is differential thermal resistance ($r$). Although $k$ and $r$ are generally positive, they can be negative in a macroscopic heterojunction with positive $T$-dependent interfacial thermal resistance (ITR). The negative differential thermal resistance (NDTR) effect is an important effect that can realize thermal transistor, thermal memory, and thermal logic gate. In this paper, we examine analytical expressions of $J$, $k$, $r$, and other related quantities as a function of parameters related to thermal conductivity ($\kappa$) and ITR in a macroscopic heterojunction to precisely describe the NDTR effect., Comment: 8 pages, 9 figures, J. Stat. Mech. (in press)

We reply to the Comment by Svhaika, Ruocco, Schirmacher and collaborators. There were two accidental mistakes in our original paper (Phys. Rev. Lett. 112, 145501 (2019)), which have been now corrected. All the physical conclusions and results of the original paper, including the prediction of boson peak due to anharmonicity, remain valid in the corrected version., Comment: Reply to Ruocco et al. (arXiv:2104.13076). We thank the authors of arXiv:2104.13076 for pointing out the mathematical mistakes in (Phys. Rev. Lett. 112, 145501 (2019)) and for motivating us to revisit our analysis. Mathematica file downloadable by all interested readers who wish to check our calculations

Identification of important topics in a text can facilitate knowledge curation, discover thematic trends, and predict future directions. In this paper, we aim to quantitatively detect the most common research themes in the emerging supercapacitor research area, and summarize their trends and characteristics through the proposed unsupervised, machine learning approach. We have retrieved the complete reference entries of article abstracts from Scopus database for all original research articles from 2004 to 2021. Abstracts were processed through a natural language processing pipeline and analyzed by a latent Dirichlet allocation topic modeling algorithm for unsupervised topic discovery. Nine major topics were further examined through topic-word associations, Inter-topic distance map and topic-specific word cloud. We observed the greatest importance is being given to performance metrics (28.2%), flexible electronics (8%), and graphene-based nanocomposites (10.9%). The analysis also points out crucial future research directions towards bio-derived carbon nanomaterials (such as RGO) and flexible supercapacitors.

Controlling the connectivity and rigidity of kirigami, i.e. the process of cutting paper to deploy it into an articulated system, is critical in the manifestations of kirigami in art, science and technology, as it provides the resulting metamaterial with a range of mechanical and geometric properties. Here we combine deterministic and stochastic approaches for the control of rigidity in kirigami using the power of $k$ choices, an approach borrowed from the statistical mechanics of explosive percolation transitions. We show that several methods for rigidifying a kirigami system by incrementally changing either the connectivity or the rigidity of individual components allow us to control the nature of the explosive transition by a choice of selection rules. Our results suggest simple lessons for the design and control of mechanical metamaterials.

The aim of this paper is to present a modeling for the rheological behavior of simple liquids as a function of the amplitude of the imposed shear stress or strain. The elastic mode theory (Ref. 6) is first generalized to take into account the fact that during a flow experiment, mechanical energy is injected in a system initially at thermodynamic equilibrium. This generalized theory can be seen as a particular aspect of the general problem of perturbation by the measurement, associated with that of the coupling between fluctuation and dissipation. This generalization leads to a "finitary" character of the model. It is then combined with the inertial mode theory (Ref. 7). The formalism thus obtained allows to model the rheological behavior of liquids over a wide range of velocity gradients, including the intermediate narrow range corresponding to the Newtonian regime. As experimental tests, viscosity measurements with two kinds of moving rotor rheometers were performed. Only data obtained with liquid water at room temperature are presented and quantitatively analyzed here. It is also shown that liquid n-octane exhibits the same qualitative behaviors as those of liquid water. In the appendices, connection of this theory with quantum mechanics and turbulence phenomena are discussed, and the notion of viscous mass is introduced., Comment: 46 pages, 23 figures, 8 tables, 3 appendices

We use vortex matter in type-II superconductors as a playground to study how different types of disorder affect the long wavelength density fluctuations of the system. We find that irrespective of the vortex-vortex interaction, in the case of samples with weak and dense point defects the system presents the hidden order of hyperuniformity characterized by an algebraic suppression of density fluctuations when increasing the system size. We also reveal that, on the contrary, for samples with planar defects hyperuniformity is suppressed since density fluctuations have a tendency to unboundedness on increasing the system size. Although some of these results were known from previous works, this paper makes the fundamental discovery that the ability of planar disorder to suppress hyperuniformity grows on increasing the softness of the structure for more diluted systems., Comment: 5 pages, 2 figures, Conference LT29, Journal of the Physical Society of Japan: Conference Proceedings accepted 2022

In a recent paper [Phys. Rev. X 12, 031021], we reported experimental observations of ``ultrastable'' states in a shear-jammed granular system subjected to small-amplitude cyclic shear. In such states, all the particle positions and contact forces are reproduced after each shear cycle so that a strobed image of the stresses and particle positions appears static. In the present work, we report further analyses of data from those experiments to characterize both global and local responses of ultrastable states within a shear cycle, not just the strobed dynamics. We find that ultrastable states follow a power-law relation between shear modulus and pressure with an exponent $\beta\approx 0.5$, reminiscent of critical scaling laws near jamming. We also examine the evolution of contact forces measured using photoelasticimetry. We find that there are two types of contacts: non-persistent contacts that reversibly open and close; and persistent contacts that never open \yiqiu{and display no measurable sliding}. We show that the non-persistent contacts make a non-negligible contribution to the emergent shear modulus. We also analyze the spatial correlations of the stress tensor and compare them to the predictions of a recent theory of the emergent elasticity of granular solids, the Vector Charge Theory of Granular mechanics and dynamics (VCTG) [Phys. Rev. Lett. 125, 118002, arXiv:2204.11811]. We show that our experimental results can be fit well by VCTG, assuming uniaxial symmetry of the contact networks. The fits reveal that the response of the ultrastable states to additional applied stress is substantially more isotropic than that of the original shear-jammed states. Our results provide important insight into the mechanical properties of frictional granular solids created by shear.