Search

Van der Waals (vdW) heterostructures have attracted intense interest worldwide that offers several routes to design materials with novel features which can be served various applications. Unfortunately, at present, this method is restricted to a small number of stackable layers. In this work, we report a bulk vdW material consisting of alternative TiS2 and TiSe2 layers vertically arranged as a potential candidate for anode applications. Moreover, with density functional theory (DFT) calculations and ab-initio molecular dynamics (AIMD) simulations, we explored the effect of high entropy on electrochemical performance of the bulk heterostructure (TSS-HS) by substituting Mo6+ and Al3+ at transition metal site (Ti4+). Solvation shell formation at the electrode-electrolyte interface (EEI) was also studied using AIMD to determine the Li-coordination. Our DFT and AIMD analysis showed considerable improvement in the performance of high entropy TSS-HS (TSS-HE) as compared to TSS-HS. Factors responsible for enhanced performance of TSS-HE were 1) Less structural deformation, 2) Strong bonding (Metal-Chalcogen), 3) Better ionic mobility, 4) Higher theoretical specific capacity, 5) Wider operational voltage window, and 6) Faster Li-ion diffusion. Our observations suggests that 'high entropy' can be an effective strategy to design new anode materials for improving electrochemical performance of Li-ion batteries.

The ability to control forces between sub-micron-scale building blocks offers considerable potential for designing new materials through self-assembly. A typical paradigm is to first identify a particular (crystal) structure that has some desired property, and then design building-block interactions so that this structure assembles spontaneously. While significant theoretical and experimental progress has been made in assembling complicated structures in a variety of systems, this two-step paradigm fundamentally fails for structurally disordered solids, which lack a well-defined structure to use as a target. Here we show that disordered solids can still be treated from an inverse self-assembly perspective by targeting material properties directly. Using the Poisson's ratio, $\nu$, as a primary example, we show how differentiable programming connects experimentally relevant interaction parameters with emergent behavior, allowing us to iteratively "train" the system until we find the set of interactions that leads to the Poisson's ratio we desire. Beyond the Poisson's ratio, we also tune the pressure and a measure of local 8-fold structural order, as well as multiple properties simultaneously, demonstrating the potential for nontrivial design in disordered solids. This approach is highly robust, transferable, and scalable, can handle a wide variety of model systems, properties of interest, and preparation dynamics, and can optimize over 100s or even 1000s of parameters. This result connects the fields of disordered solids and inverse self-assembly, indicating that many of the tools and ideas that have been developed to understand the assembly of crystals can also be used to control the properties of disordered solids.

M\"ossbauer spectroscopy is a technique employed to investigate the microscopic properties of materials using transitions between energy levels in the nuclei. Conventionally, in synchrotron-radiation-based M\"ossbauer spectroscopy, the measurement window is decided by the researcher heuristically, although this decision has a significant impact on the shape of the measurement spectra. In this paper, we propose a method for evaluating the precision of the spectral position by introducing Bayesian estimation. The proposed method makes it possible to select the best measurement window by calculating the precision of M\"ossbauer spectroscopy from the data. Based on the results, the precision of the M\"ossbauer center shifts improved by more than three times compared with the results achieved with the conventional simple fitting method using the Lorentzian function., Comment: 8 pages, 8 figures

Representation learning for the electronic structure problem is a major challenge of machine learning in computational condensed matter and materials physics. Within quantum mechanical first principles approaches, Kohn-Sham density functional theory (DFT) is the preeminent tool for understanding electronic structure, and the high-dimensional wavefunctions calculated in this approach serve as the building block for downstream calculations of correlated many-body excitations and related physical observables. Here, we use variational autoencoders (VAE) for the unsupervised learning of high-dimensional DFT wavefunctions and show that these wavefunctions lie in a low-dimensional manifold within the latent space. Our model autonomously determines the optimal representation of the electronic structure, avoiding limitations due to manual feature engineering and selection in prior work. To demonstrate the utility of the latent space representation of the DFT wavefunction, we use it for the supervised training of neural networks (NN) for downstream prediction of the quasiparticle bandstructures within the GW formalism, which includes many-electron correlations beyond DFT. The GW prediction achieves a low error of 0.11 eV for a combined test set of metals and semiconductors drawn from the Computational 2D Materials Database (C2DB), suggesting that latent space representation captures key physical information from the original data. Finally, we explore the interpretability of the VAE representation and show that the successful representation learning and downstream prediction by our model is derived from the smoothness of the VAE latent space, which also enables the generation of wavefunctions on arbitrary points in latent space. Our work provides a novel and general machine-learning framework for investigating electronic structure and many-body physics.

The study of phonon dynamics is pivotal for understanding material properties, yet it faces challenges due to the irreversible information loss inherent in powder inelastic neutron scattering spectra and the limitations of traditional analysis methods. In this study, we present a machine learning framework designed to reveal obscured phonon dynamics from powder spectra. Using a variational autoencoder, we obtain a disentangled latent representation of spectra and successfully extract force constants for reconstructing phonon dispersions. Notably, our model demonstrates effective applicability to experimental data even when trained exclusively on physics-based simulations. The fine-tuning with experimental spectra further mitigates issues arising from domain shift. Analysis of latent space underscores the model's versatility and generalizability, affirming its suitability for complex system applications. Furthermore, our framework's two-stage design is promising for developing a universal pre-trained feature extractor. This approach has the potential to revolutionize neutron measurements of phonon dynamics, offering researchers a potent tool to decipher intricate spectra and gain valuable insights into the intrinsic physics of materials., Comment: 22 pages, 9 figures

We study the phase behaviour of the Al$_x$CrFeCoNi high-entropy alloy. Our approach is based on a perturbative analysis of the internal energy of the paramagnetic solid solution as evaluated within the Korringa-Kohn-Rostoker formulation of density functional theory, using the coherent potential approximation to average over disorder. Via application of a Landau-type linear response theory, we infer preferential chemical orderings directly. In addition, we recover a pairwise form of the alloy internal energy suitable for study via atomistic simulations, which in this work are performed using the nested sampling algorithm, which is well-suited for studying complex potential energy surfaces. When the underlying lattice is fcc, at low concentrations of Al, depending on the value of $x$, we predict either an $\mathrm{L}1_2$ or $\mathrm{D}0_{22}$ ordering emerging below approximately 1000 K. On the other hand, when the underlying lattice is bcc, consistent with experimental observations, we predict $\mathrm{B}2$ ordering temperatures higher than the melting temperature of the alloy, confirming that this ordered phase forms directly from the melt. For both fcc and bcc systems, chemical orderings are dominated by Al moving to one sublattice, Ni and Co the other, while Cr and Fe remain comparatively disordered. On the bcc lattice, our atomistic modelling suggests eventual decomposition into $\mathrm{B}2$ NiAl and Cr-rich phases. These results shed light on the fundamental physical origins of atomic ordering tendencies in these intriguing materials., Comment: 17 pages, 11 figures, 1 table

Layered halide perovskites have emerged as potential alternatives to three-dimensional halide perovskites due to their improved stability and larger material phase space, allowing fine-tuning of structural, electronic, and optical properties. However, their charge carrier mobilities are significantly smaller than that of three-dimensional halide perovskites, which has a considerable impact on their application in optoelectronic devices. Here, we employ state-of-the-art ab initio approaches to unveil the electron-phonon mechanisms responsible for the diminished transport properties of layered halide perovskites. Starting from a prototypical ABX$_{3}$ halide perovskite, we model the case of $n=1$ and $n=2$ layered structures and compare their electronic and transport properties to the three-dimensional reference. The electronic and phononic properties are investigated within density functional theory (DFT) and density functional perturbation theory (DFPT), while transport properties are obtained via the ab initio Boltzmann transport equation. The vibrational modes contributing to charge carrier scattering are investigated and associated with polar-phonon scattering mechanisms arising from the long-range Fr\"ohlich coupling and deformation potential scattering processes. Our investigation reveals that the lower mobilities in layered systems primarily originates from the increased electronic density of states at the vicinity of the band edges, while the electron-phonon coupling strength remains similar. Such increase is caused by the dimensionality reduction and the break in octahedra connectivity along the stacking direction. Our findings provide a fundamental understanding of the electron-phonon coupling mechanisms in layered perovskites and highlight the intrinsic limitations of the charge carrier transport in these materials.

The crystallographic texture is a key organization feature of many technical and biological materials. In these materials, especially hierarchically structured ones, the preferential alignment of the nano constituents is heavily influencing the macroscopic behaviour of the material. In order to study local crystallographic texture with both high spatial and angular resolution, we developed Texture tomography (TexTOM). This approach allows to model the diffraction data of polycrystalline materials by using the full reciprocal space of the ensemble of crystals and describe the texture in each voxel via a orientation distribution function. This means, it provides 3D reconstructions of the local texture by measuring the probabilities of all crystal orientations. The TexTOM approach addresses limitations associated with existing models: It correlates the intensities from several Bragg reflections, thus reduces ambiguities resulting from symmetry. Further, it yields quantitative probability distributions of local real space crystal orientations without further assumptions on the sample structure. Finally, its efficient mathematical formulation enables reconstructions faster than the time-scale of the experiment. In this manuscript, we present the mathematical model, the inversion strategy and its current experimental implementation. We show characterizations of simulated data as well as experimental data obtained from a synthetic, inorganic model sample, the silica-witherite biomorph. In conclusion, Tex-TOM provides a versatile framework to reconstruct 3D quantitative texture information for polycrystalline samples. In this way, it opens the door for unprecedented insights into the nanostructural makeup of natural and technical materials.

We present an advanced scanning probe microscopy system enhanced with artificial intelligence (AI-SPM) designed for self-driving atomic-scale measurements. This system expertly identifies and manipulates atomic positions with high precision, autonomously performing tasks such as spectroscopic data acquisition and atomic adjustment. An outstanding feature of AI-SPM is its ability to detect and adapt to surface defects, targeting or avoiding them as necessary. It's also engineered to address typical challenges such as positional drift and tip apex atomic variations due to the thermal effect, ensuring accurate, site-specific surface analyses. Our tests under the demanding conditions of room temperature have demonstrated the robustness of the system, successfully navigating thermal drift and tip fluctuations. During these tests on the Si(111)-(7x7) surface, AI-SPM autonomously identified defect-free regions and performed a large number of current-voltage spectroscopy measurements at different adatom sites, while autonomously compensating for thermal drift and monitoring probe health. These experiments produce extensive data sets that are critical for reliable materials characterization and demonstrate the potential of AI-SPM to significantly improve data acquisition. The integration of AI into SPM technologies represents a step toward more effective, precise and reliable atomic-level surface analysis, revolutionizing materials characterization methods.

We derive the two-time linear response theory for out-of-equilibrium pumped systems, generic pump-probe delays and probe frequencies. Such a theory enormously simplifies the numerical calculations, for instance, of the optical conductivity with respect to the actual procedure, which requires computing the effect of the probe pulse for each time delay with respect to the pump pulse\textcolor{red}{.} The theory is given for a generic observable and pumped Hamiltonian and then specialized for a system with a quadratic Hamiltonian and its transient optical properties, exploiting the Dynamical Projective Operatorial Approach (DPOA). The theory is complemented by a set of crucial numerical guidelines that help perform actual calculations in a computationally affordable way. The optical response (differential transient reflectivity and absorption) of a prototypical three-band (core, valence, and conduction) model in the XUV regime is analyzed in detail to illustrate the theory and its application. Using some generalizations of the density of states, we provide a systematic approach to exploring the optical properties in terms of the system band structure features and the pump parameters. Such an analysis can be extremely helpful in understanding the actual results of experimental optical measurements. Moreover, we study the effects of inter-band and intra-band transitions, the local dipole coupling, and single and multi-photon processes. The latter is further investigated by varying the central frequency of the pump pulse to have different regions of the first Brillouin zone in resonance with it. We also study the effect of varying the pump pulse intensity. Finally, we study and analyze the transient optical properties in the probe pulse regime of IR and visible., Comment: 25 pages, 11 figures, 40 panels

Emerging two-dimensional (2D) materials bring unprecedented opportunities for electronic applications. The design of high-performance devices requires an accurate prediction of carrier mobility in 2D materials, which can be obtained using state-of-the-art $ab~initio$ calculations. However, various factors impact the computational accuracy, leading to contradictory estimations for the mobility. In this work, targeting accurate and efficient $ab~initio$ calculations, transport properties in III-V monolayers are reported using the Boltzmann transport equation, and the influences of pseudopotential, quadrupole correction, Berry connection, and spin-orbit coupling (SOC) on mobilities are systematically investigated. Our findings are as follows: (1) The inclusion of semi-core states in pseudopotentials is important to obtain accurate calculations. (2) The variations induced by dynamical quadrupole and Berry connection when treating long range fields can be respectively 40% and 10%. (3) The impact of SOC can reach up to 100% for materials with multi-peak bands. Importantly, although SOC notably modifies the electronic wavefunctions, it negligibly impacts the dynamical matrices and scattering potential variations. As a result, the combination of fully-relativistic electron calculation and scalar-relativistic phonon calculation can strike a good balance between accuracy and cost. This work compares computational methodologies, providing guidelines for accurate and efficient calculations of mobilities in 2D semiconductors., Comment: 11 pages, 15 figures

In this article, we present an interpolative separable density fitting (ISDF) based algorithm to calculate exact exchange in periodic mean field calculations. In the past, decomposing the two-electron integrals into tensor hypercontraction (THC) form using ISDF was the most expensive step of the entire mean field calculation. Here we show that by using a multigrid-ISDF algorithm both the memory and the CPU cost of this step can be reduced. The CPU cost is brought down from cubic scaling to quadratic scaling with a low computational prefactor which reduces the cost by almost two orders of magnitude. Thus, in the new algorithm, the cost of performing ISDF is largely negligible compared to other steps. Along with the CPU cost, the memory cost of storing the factorized two-electron integrals is also reduced by a factor of up to 35. With the current algorithm, we can perform Hartree-Fock calculations on a Diamond supercell containing more than 17,000 basis functions and more than 1,500 electrons on a single node with no disk usage. For this calculation, the cost of constructing the exchange matrix is only a factor of four slower than the cost of diagonalizing the Fock matrix. Augmenting our approach with linear scaling algorithms can further speed up the calculations.

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.

Correlated systems represent a class of materials that are difficult to describe through traditional electronic structure methods. The computational demand to simulate the structural dynamics of such systems, with correlation effects considered, is substantial. Here, we investigate the structural dynamics of $f$- and $d$-electron correlated systems by integrating quantum embedding techniques with interatomic potentials derived from graph neural networks. For Cerium, a prototypical correlated $f$-electron system, we use Density Functional Theory with the Gutzwiller approximation to generate training data due to efficiency with which correlations effects are included for large multi-orbital systems. For Nickel Oxide, a prototypical correlated $d$-electron system, advancements in computational capabilities now permit the use of full Dynamical Mean Field Theory to obtain energies and forces. We train neural networks on this data to create a model of the potential energy surface, enabling rapid and effective exploration of structural dynamics. Utilizing these potentials, we delineate transition pathways between the $\alpha$, $\alpha'$, and $\alpha''$ phases of Cerium and predict the melting curve of Nickel Oxide. Our results demonstrate the potential of machine learning potentials to accelerate the study of strongly correlated systems, offering a scalable approach to explore and understand the complex physics governing these materials., Comment: 7 pages, 4 figures

The success behind many pseudopotential methods such as Projected Augmented Waves (PAW) and Phillips-Kleinman pseudopotential methods is that these methods are all electron methods in disguise. For the Phillips Kleinman and PAW methods we show that there is an all electron formulation. In the all electron formulation there are regular low wavevector wavefunctions and several, specially chosen, high wavevector wavefunctions that are specialized to the atomic environment. Using this as inspiration here we propose a new, PAW pseudopotential inspired, basis set method for small molecules, where we use the recently introduced LO or lo (Localized orbitals) basis set is paired with a Gaussian basis for a hybrid basis for small molecules., Comment: Comments welcome, V2:Minor typos corrected

Over the past decade, experimental microscopy and spectroscopy have made significant progress in the study of the morphological, optical, electronic and transport properties of materials. These developments include higher spatial resolution, shorter acquisition times, more efficient monochromators and electron analysers, improved contrast imaging and advancements in sample preparation techniques. These advances have driven the need for more accurate theoretical descriptions and predictions of material properties. Computer simulations based on first principles and Monte Carlo methods have emerged as a rapidly growing field for modeling the interaction of charged particles, such as electron, proton and ion beams, with various systems, such as slabs, nanostructures and crystals. This report delves into the theoretical and computational approaches to modeling the physico-chemical mechanisms that occur when charged beams interact with a medium. These mechanisms encompass single and collective electronic excitation, ionization of the target atoms and the generation of a secondary electron cascade that deposits energy into the irradiated material. We show that the combined application of ab initio methods, which are able to model the dynamics of interacting many-fermion systems, and Monte Carlo methods, which capture statistical fluctuations in energy loss mechanisms by random sampling, proves to be an optimal strategy for the accurate description of charge transport in solids. This joint quantitative approach enables the theoretical interpretation of excitation, loss and secondary electron spectra, the analysis of the chemical composition and dielectric properties of solids and contributes to our understanding of irradiation-induced damage in materials, including those of biological significance., Comment: 142 pages, 45 figures, 6 tables, review paper

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.

Strong laser pulses can control superconductivity, inducing non-equilibrium transient pairing by leveraging strong-light matter interaction. Here we demonstrate theoretically that equilibrium ground-state phonon-mediated superconductive pairing can be affected through the vacuum fluctuating electromagnetic field in a cavity. Using the recently developed $\textit{ab initio}$ quantum electrodynamical density-functional theory approximation, we specifically investigate the phonon-mediated superconductive behavior of MgB$_2$ under different cavity setups and find that in the strong-light matter coupling regime its superconducting transition temperature can be enhanced by $\approx 73\%$ ($\approx 40\%$) in an in-plane (out-of-plane) polarized cavity. The results highlight that strong-light matter coupling in extended systems can profoundly alter material properties in a non-perturbative way by modifying their electronic structure and phononic dispersion at the same time. Our findings indicate a pathway to the experimental realization of light-controlled superconductivity in solid-state materials at equilibrium via cavity-material engineering., Comment: 25 pages, 3 figures

We develop a framework for on-the-fly machine learned force field molecular dynamics simulations based on the multipole featurization scheme that overcomes the bottleneck with the number of chemical elements. Considering bulk systems with up to 6 elements, we demonstrate that the number of density functional theory calls remains approximately independent of the number of chemical elements, in contrast to the increase in the smooth overlap of atomic positions scheme., Comment: 15 total pages including references and SI, 4 main figures, 1 supplemental, 6 supplemental tables. Submitted to ACS Journal of Chemical Theory and Computation. Written and reviewed by all

Predicting the charged particle transport properties of warm dense matter / hot dense plasma mixtures is a challenge for analytical models. High accuracy ab initio methods are more computationally expensive, but can provide critical insight by explicitly simulating mixtures. In this work, we investigate the transport properties and optical response of warm dense carbon-hydrogen mixtures at varying concentrations under either conserved electronic pressure or mass density at a constant temperature. We compare options for mixing the calculated pure species properties to estimate the results of the mixtures. We find that a combination of the Drude model with the Matthiessen's rule works well for DC electron transport and low frequency optical response. This breaks down at higher frequencies, where a volumetric mix of pure-species AC conductivities works better., Comment: Special Collection: Charged-Particle Transport in High Energy Density Plasmas

The insulator-to-metal transition (IMT) in vanadium dioxide (VO2) has garnered extensive attention for its potential applications in ultrafast switches, neuronal network architectures, and storage technologies. However, a significant controversy persists regarding the formation of the IMT, specifically concerning whether a complete structural phase transition from monoclinic (M1) to rutile (R) phase is necessary. Here we employ the real-time time-dependent density functional theory (rt-TDDFT) to track the dynamic evolution of atomic and electronic structures in photoexcited VO2, revealing the emergence of a long-lived monoclinic metal phase (MM) under low electronic excitation. The emergence of the metal phase in the monoclinic structure originates from the dissociation of the local V-V dimer, driven by the self-trapped and self-amplified dynamics of photoexcited holes, rather than by a pure electron-electron correction. On the other hand, the M1-to-R phase transition does appear at higher electronic excitation. Our findings validate the existence of MM phase and provide a comprehensive picture of the IMT in photoexcited VO2., Comment: 34 pages, 18 figures

Despite the widespread applications of machine learning force field (MLFF) on solids and small molecules, there is a notable gap in applying MLFF to complex liquid electrolytes. In this work, we introduce BAMBOO (ByteDance AI Molecular Simulation Booster), a novel framework for molecular dynamics (MD) simulations, with a demonstration of its capabilities in the context of liquid electrolytes for lithium batteries. We design a physics-inspired graph equivariant transformer architecture as the backbone of BAMBOO to learn from quantum mechanical simulations. Additionally, we pioneer an ensemble knowledge distillation approach and apply it on MLFFs to improve the stability of MD simulations. Finally, we propose the density alignment algorithm to align BAMBOO with experimental measurements. BAMBOO demonstrates state-of-the-art accuracy in predicting key electrolyte properties such as density, viscosity, and ionic conductivity across various solvents and salt combinations. Our current model, trained on more than 15 chemical species, achieves the average density error of 0.01 g/cm$^3$ on various compositions compared with experimental data. Moreover, our model demonstrates transferability to molecules not included in the quantum mechanical dataset. We envision this work as paving the way to a "universal MLFF" capable of simulating properties of common organic liquids.

Members of the DFT+U family of functionals are increasingly prevalent methods of addressing errors intrinsic to (semi-) local exchange-correlation functionals at minimum computational cost, but require their parameters U and J to be calculated in situ for a given system of interest, simulation scheme, and runtime parameters. The SCF linear response approach offers ab initio acquisition of the U and has recently been extended to compute the J analogously, which measures localized errors related to exchange-like effects. We introduce a renovated post-processor, the lrUJ utility, together with this detailed best-practices guide, to enable users of the popular, open-source Abinit first-principles simulation suite to engage easily with in situ Hubbard parameters and streamline their incorporation into material simulations of interest. Features of this utility, which may also interest users and developers of other DFT codes, include $n$-degree polynomial regression, error analysis, Python plotting facilities, didactic documentation, and avenues for further developments. In this technical introduction and guide, we place particular emphasis on the intricacies and potential pitfalls introduced by the projector augmented wave (PAW) method, SCF mixing schemes, and non-linear response, several of which are translatable to DFT+U(+J) implementations in other packages., Comment: 40 pages, 8 figures (1 figure in appendix), 3 tables

We introduce an electron-photon exchange-correlation functional for quantum electrodynamical density-functional theory (QEDFT). The approach, photon MBD (pMBD), is inspired by the many-body dispersion (MBD) method for weak intermolecular interactions, which is generalized to include both electronic and photonic (electromagnetic) degrees of freedom on the same footing. We demonstrate that pMBD accurately captures effects that arise in the context of strong light-matter interactions, such as anisotropic electron-photon interactions, beyond single-photon effects, and cavity modulated van der Waals interactions. Moreover, we show that pMBD is computationally efficient and allows simulations of large complex systems coupled to optical cavities., Comment: 10 pages, 3 figures

Structure-informed materials informatics is a rapidly evolving discipline of materials science relying on the featurization of atomic structures or configurations to construct vector, voxel, graph, graphlet, and other representations useful for machine learning prediction of properties, fingerprinting, and generative design. This work discusses how current featurizers typically perform redundant calculations and how their efficiency could be improved by considering (1) fundamentals of crystallographic (orbits) equivalency to optimize ordered cases and (2) representation-dependent equivalency to optimize cases of dilute, doped, and defect structures with broken symmetry. It also discusses and contrasts ways of (3) approximating random solid solutions occupying arbitrary lattices under such representations. Efficiency improvements discussed in this work were implemented within pySIPFENN or python toolset for Structure-Informed Property and Feature Engineering with Neural Networks developed by authors since 2019 and shown to increase performance from 2 to 10 times for typical inputs. Throughout this work, the authors explicitly discuss how these advances can be applied to different kinds of similar tools in the community., Comment: 4 figures; GitHub repository at https://git.pysipfenn.org; Documentation at https://pysipfenn.org

Since the implementation of the Materials Genome Project by the Obama administration in the United States, the development of various computational materials databases has fundamentally expanded the choices of industries such as materials and energy. In the field of thermoelectric materials, the thermoelectric figure of merit ZT quantifies the performance of the material. From the viewpoint of calculations for vast materials, the ZT values are not easily obtained due to their computational complexity. Here, we show how to build a database of thermoelectric materials based on first-principles calculations for the electronic and heat transport of materials. Firstly, the initial structures are classified according to the values of bandgap and other basic properties using the clustering algorithm K-means in machine learning, and high-throughput first principles calculations are carried out for narrow-bandgap semiconductors which exhibiting potential thermoelectric application. The present framework of calculations mainly includes deformation potential module, electrical transport performance module, mechanical and thermodynamic properties module. We have also set up a search webpage for the calculated database of thermoelectric materials, providing searching and viewing the related physical properties of materials. Our work may inspire the construction of more computational databases of first-principle thermoelectric materials and accelerate research progress in the field of thermoelectrics., Comment: 13 pages, 5 figures

Given the vast compositional possibilities Na$_nM_mM_{m'}$Si$_{3-p-a}$P$_p$As$_a$O$_{12}$, Na-ion superionic conductors (NASICON) are attractive but complicate for designing materials with enhanced room-temperature Na-ion conductivity $\sigma_{\rm Na,300K}$. We propose an explicit regression model for $\sigma_{\rm Na,300K}$ with easily-accessible descriptors, by exploiting density functional theory molecular dynamics (DFT-MD). Initially, we demonstrate that two primary descriptors, the bottleneck width along Na-ion diffusion paths $d_1$ and the average Na-Na distance $\langle d_{\rm Na-Na} \rangle$, modulate room-temperature Na-ion self-diffusion coefficient $D_{\rm Na,300K}$. Then, we introduce two secondary easily-accessible descriptors: Na-ion content n, which influences $d_1$, $\langle d_{\rm Na-Na} \rangle$, and Na-ion density $\rho_{\rm Na}$; and the average ionic radius $\langle r_M \rangle$ of metal ions, which impacts $d_1$ and $\langle d_{\rm Na-Na} \rangle$. These secondary descriptors enable the development of a regression model for $\sigma_{\rm Na,300K}$ with $n$ and $\langle r_M \rangle$ only. Subsequently, this model identifies a promising yet unexplored stable composition, Na$_{2.75}$Zr$_{1.75}$Nb$_{0.25}$Si$_2$PO$_{12}$, which, upon DFT-MD calculations, indeed exhibits $\sigma_{\rm Na,300K} > 10^{^3}$ S$\cdot$cm$^{-1}$. Furthermore, the adjusted version effectively fits $140$ experimental values with $R^2=0.718$.

Sulphur vacancies in MoS$_2$ on Au(111) have been shown to be negatively charged as reflected by a Kondo resonance. Here, we use scanning tunneling microscopy to show that these vacancies serve as anchoring sites for thiol-based molecules (CF$_3$-3P-SH) with two distinct reaction products, one of them showing a Kondo resonance. Based on comparisons with density-functional theory (DFT) calculations, including a random structure search and computation of energies and electronic properties at a hybrid exchange-correlation functional level, we conclude that both anchored molecules are charge neutral. One of them is an anchored intact CF$_3$-3P-SH molecule while the other one is the result of catalytically activated dehydrogenation to CF$_3$-3P-S with subsequent anchoring. Our investigations highlight a perspective of functionalizing defects with thiol-terminated molecules that can be equipped with additional functional groups, such as charge donor- or acceptor-moieties, switching units or magnetic centers.

Identifying local structural motifs and packing patterns of molecular solids is a challenging task for both simulation and experiment. We demonstrate two novel approaches to characterize local environments in different polymorphs of molecular crystals using learning models that employ either flexibly learned or handcrafted molecular representations. In the first case, we follow our earlier work on graph learning in molecular crystals, deploying an atomistic graph convolutional network, combined with molecule-wise aggregation, to enable per-molecule environmental classification. For the second model, we develop a new set of descriptors based on symmetry functions combined with a point-vector representation of the molecules, encoding information about the positions as well as relative orientations of the molecule. We demonstrate very high classification accuracy for both approaches on urea and nicotinamide crystal polymorphs, and practical applications to the analysis of dynamical trajectory data for nanocrystals and solid-solid interfaces. Both architectures are applicable to a wide range of molecules and diverse topologies, providing an essential step in the exploration of complex condensed matter phenomena.

Electronic Structure Theory (EST) describes the behavior of electrons in matter and is used to predict material properties. Conventionally, this involves forming a Hamiltonian and solving the Schr\"odinger equation through discrete computation. Here, a new perspective to EST is provided by treating a perfectly crystalline material as a Linear Translation Invariant (LTI) system. The validity of this LTI-EST formalism is demonstrated by determining band structures for a one-dimensional chain of atoms, including the phenomenon of band structure folding in super cells. The proposed formalism allows for analytical traceability of band structure folding and offers computational advantage by bypassing the O(N) eigenvalue calculations. The spike-based computing nature of the proposed LTI-EST formalism is highlighted; thereby implying potential for material simulations solely in the spiking domain.

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory (DFT) MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity auto-correlation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

Training machine learning interatomic potentials often requires optimizing a loss function composed of three variables: potential energies, forces, and stress. The contribution of each variable to the total loss is typically weighted using fixed coefficients. Identifying these coefficients usually relies on iterative or heuristic methods, which may yield sub-optimal results. To address this issue, we propose an adaptive loss weighting algorithm that automatically adjusts the loss weights of these variables during the training of potentials, dynamically adapting to the characteristics of the training dataset. The comparative analysis of models trained with fixed and adaptive loss weights demonstrates that the adaptive method not only achieves a more balanced predictions across the three variables but also improves overall prediction accuracy.

The combination of isochoric heating of solids by free electron lasers (FEL) and in situ diagnostics by X-ray Thomson scattering (XRTS) allows for measurements of material properties at warm dense matter (WDM) conditions relevant for astrophysics, inertial confinement fusion, and material science. In the case of metals, the FEL beam pumps energy directly into electrons with the lattice structure of ions being nearly unaffected. This leads to a unique transient state that gives rise to a set of interesting physical effects, which can serve as a reliable testing platform for WDM theories. In this work, we present extensive linear-response time-dependent density functional theory (TDDFT) results for the electronic dynamic structure factor of isochorically heated copper with a face-centered cubic lattice. At ambient conditions, the plasmon is heavily damped due to the presence of $d$-band excitations, and its position is independent of the wavenumber. In contrast, the plasmon feature starts to dominate the excitation spectrum and has a Bohm-Gross type plasmon dispersion for temperatures $T \geq 4~{\rm eV}$, where the quasi-free electrons in the interstitial region are in the WDM regime. In addition, we analyze the thermal changes in the $d$-band excitations and outline the possibility to use future XRTS measurements of isochorically heated copper as a controlled testbed for WDM theories.

Noble gas atoms such as Kr and Xe are byproducts of nuclear fission in nuclear plants. How to trap and confine these volatile even radioactive gases is particularly challenging. Recent studies have shown that they can be trapped in nanocages of ultrathin silica. Here, we exhibit with self-consistent phonon theory and four-phonon (4ph) scattering where the adsorption of noble gases results in an anomalous increase in lattice thermal conductivity, while the presence of Cu atoms doping leads to a reduction in lattice thermal conductivity. We trace this behavior in host-guest 2D silica to an interplay of tensile strain, rattling phonon modes, and redistribution of electrons. We also find that 4ph scatterings play indispensable roles in the lattice thermal conductivity of 2D silica. Our work illustrates the microscopic heat transfer mechanism in 2D silica nanocages with the immobilization of noble gas atoms and inspires further exploring materials with the kagome and glasslike lattice thermal conductivity., Comment: 7 pages, 4 figures

Designing polymers with high intrinsic thermal conductivity (TC) is critically important for the thermal management of organic electronics and photonics. However, this is a challenging task owing to the diversity of the chemical space and the barriers to advanced synthetic experiments/characterization techniques for polymers. In this Tutorial, the fundamentals and implementation of combining classical molecular dynamics simulation and machine learning (ML) for the development of polymers with high TC are comprehensively introduced. We begin by describing the core components of a universal ML framework, involving polymer datasets, property calculators, feature engineering and informatics algorithms. Then, the process of constructing interpretable regression algorithms for TC prediction is introduced, aiming to extract the underlying relationships between microstructures and TCs for polymers. We also explore the design of sequence-ordered polymers with high TC using lightweight and mainstream active learning algorithms. Lastly, we conclude by addressing the current limitations and suggesting potential avenues for future research on this topic.

We introduce CrystalFormer, a transformer-based autoregressive model specifically designed for space group-controlled generation of crystalline materials. The space group symmetry significantly simplifies the crystal space, which is crucial for data and compute efficient generative modeling of crystalline materials. Leveraging the prominent discrete and sequential nature of the Wyckoff positions, CrystalFormer learns to generate crystals by directly predicting the species and locations of symmetry-inequivalent atoms in the unit cell. Our results demonstrate that CrystalFormer matches state-of-the-art performance on standard benchmarks for both validity, novelty, and stability of the generated crystalline materials. Our analysis also shows that CrystalFormer ingests sensible solid-state chemistry information from data for generative modeling. The CrystalFormer unifies symmetry-based structure search and generative pre-training in the realm of crystalline materials. The simplicity, generality, and flexibility of CrystalFormer position it as a promising architecture to be the foundational model of the entire crystalline materials space, heralding a new era in materials modeling and discovery., Comment: 17 pages, 8 figures

Exploring ultrafast magnetization control in two-dimensional (2D) magnets through optically driven coherent phonons has been well-established. Yet, the microscopic interplay between spin dynamics and lattice degrees of freedom remains less explored. Employing real-time time-dependent density functional theory (rt-TDDFT) coupled with Ehrenfest dynamics, we systematically investigate laser-induced spin-nuclei dynamics with coherent phonon excitation in the 2D ferromagnet Fe3GeTe2. We found that selectively pre-exciting three typical coherent phonon modes results in up to a 53% additional spin moment loss in an out-of-plane A2 1g mode within ~50 fs. Coherent phonon control of spin dynamics is closely linked to laser pulse parameters. The underlying microscopic mechanism of this phenomenon is primarily governed by coherent phonon-induced asymmetric spin-resolved charge transfer following the disappearance of the laser pulse, thereby enabling effective control of the spin moment loss. Our findings offer a novel insight into the coupling of coherent phonons with spin systems in 2D limits on femtosecond timescales.

Phase transition problems on curved surfaces can lead to a panopticon of fascinating patterns. In this paper we consider finite element approximations of phase field models with a spatially inhomogeneous and anisotropic surface energy density. The problems are either posed in $\mathbb R^3$ or on a two-dimensional hypersurface in $\mathbb R^3$. In the latter case, a fundamental choice regarding the anisotropic energy density has to be made. Our numerical method can be employed both situations, where for the problems on hypersurfaces the algorithm uses parametric finite elements. We prove an unconditional stability result for our schemes and present several numerical experiments, including for the modelling of ice crystal growth on a sphere., Comment: 27 pages, 18 figures

Ensuring solid-state lithium batteries perform well across a wide temperature range is crucial for their practical use. Molecular dynamics (MD) simulations can provide valuable insights into the temperature dependence of the battery materials, however, the high computational cost of ab initio MD poses challenges for simulating ion migration dynamics at low temperatures. To address this issue, accurate machine-learning interatomic potentials were trained, which enable efficient and reliable simulations of the ionic diffusion processes in Li6PS5Cl over a large temperature range for long-time evolution. Our study revealed the significant impact of subtle lattice parameter variations on Li+ diffusion at low temperatures and identified the increasing influence of surface contributions as the temperature decreases. Our findings elucidate the factors influencing low temperature performance and present strategic guidance towards improving the performance of solid-state lithium batteries under these conditions.

The design and discovery of new materials is fundamental to advancing scientific and technological innovation. The recent emergence of the materials genome concept holds great promise in revolutionising materials science by enabling the systematic utilisation of data for efficient prediction and optimisation of superior materials. However, the materials genome approach can be stymied by the vast complexity of design spaces, which often demand substantial computational resources and sophisticated data processing capabilities. To address these challenges, this work introduces a novel generative design framework called the non-dominant sorting optimisation-based generative adversarial networks (NSGAN). Capitalising on the synergies of genetic algorithms (GA) and generative adversarial networks (GANs), NSGAN provides a robust and efficient approach for tackling high-dimensional multi-objective optimisation design problems. To validate the efficacy of the proposed framework, we applied the model to a comprehensive dataset of aluminium alloys. Additionally, an online tool was created as a supplementary resource, offering a brief introduction to this innovative method for the wider scientific community. This study explores the potential of a predictive and data-driven approach in material design, indicating a promising pathway for widespread applications in the field of materials science., Comment: 28 pages

Interacting defect systems are ubiquitous in materials under realistic scenarios, yet gaining an atomic-level understanding of these systems from a computational perspective is challenging - it often demands substantial resources due to the necessity of employing supercell calculations. While machine learning techniques have shown potential in accelerating materials simulations, their application to systems involving interacting defects remains relatively rare. In this work, we present a comparative study of three different methods to predict the free energy change of systems with interacting defects. We leveraging a limited dataset from Density Functional Theory(DFT) calculations to assess the performance models using materials descriptors, graph neural networks and cluster expansion. Our findings indicate that the cluster expansion model can achieve precise energetics predictions even with this limited dataset. Furthermore, with synthetic data generate from cluster expansion model at near-DFT levels, we obtained enlarged dataset to assess the demands on data for training accurate prediction models using graph neural networks for systems featuring interacting defects. A brief discussion of the computational cost for each method is provided at the end. This research provide a preliminary evaluation of applying machine learning techniques in imperfect surface systems.

Predicting physical properties of materials from their crystal structures is a fundamental problem in materials science. In peripheral areas such as the prediction of molecular properties, fully connected attention networks have been shown to be successful. However, unlike these finite atom arrangements, crystal structures are infinitely repeating, periodic arrangements of atoms, whose fully connected attention results in infinitely connected attention. In this work, we show that this infinitely connected attention can lead to a computationally tractable formulation, interpreted as neural potential summation, that performs infinite interatomic potential summations in a deeply learned feature space. We then propose a simple yet effective Transformer-based encoder architecture for crystal structures called Crystalformer. Compared to an existing Transformer-based model, the proposed model requires only 29.4% of the number of parameters, with minimal modifications to the original Transformer architecture. Despite the architectural simplicity, the proposed method outperforms state-of-the-art methods for various property regression tasks on the Materials Project and JARVIS-DFT datasets., Comment: 13 main pages, 3 figures, 4 tables, 10 appendix pages. Published as a conference paper at ICLR 2024. For more information, see https://omron-sinicx.github.io/crystalformer/

This study delves into the complex mechanisms of defect production and accumulation in nanocrystalline tungsten under irradiation, with a particular focus on the interplay between grain boundaries and free surfaces. Through molecular dynamics simulations, the research explores how grain boundaries act as sinks for irradiation-induced defects, a critical aspect in designing nanomaterials with enhanced radiation tolerance. The investigation leverages a novel Modified Wigner-Seitz Analysis to accurately quantify defect trends amidst dynamic surface reconstruction, providing a nuanced understanding of defect distribution in response to irradiation. This methodology underscores the intricate relationship between defect dynamics and the nano-scale structure of materials, specifically highlighting the role of interfaces in mediating these dynamics. The findings reveal a complex balance between defect production, surface interactions, and the influence of pre-existing defects and temperature on the primary defect production process. Surfaces are shown to amplify defect production due to biased accumulation of interstitials, alongside suppressed defect recombination, emphasizing the nuanced nature of defect dynamics in irradiated materials. This research contributes significantly to the fundamental understanding of defect formation and evolution in irradiated tungsten, offering insights that are instrumental in the development of nanomaterials poised for applications in extreme irradiation environments, such as fusion reactors, thereby advancing the field of materials science and engineering.

This study explores the impact of temperature on defect dynamics in tungsten, emphasizing its application in nuclear fusion reactors as Plasma Facing Components (PFCs). Through atomistic simulations, the research elucidates the intricate interplay of defect production, annihilation, and redistribution under irradiation at room (300K) and elevated temperatures (1000K). It demonstrates that higher temperatures significantly increase the number and mobility of defects, leading to a substantial rise in the total number of surviving Frenkel Pairs (FPs), with a notable preference for surface distribution. This redistribution is attributed to energy gradient-driven relocation processes, enhanced by the defects' increased mobility at elevated temperatures. Moreover, the study reveals that elevated temperatures promote biased accumulation of interstitial defects in grain boundaries, especially in configurations that facilitate efficient interstitial migration, indicating a strategy for minimizing lattice interstitial accumulation under irradiation. These findings underscore the critical role of temperature in modulating irradiation-induced defect dynamics in tungsten, providing valuable insights for designing and selecting materials with optimized irradiation resistance for use in extreme conditions of nuclear fusion reactors. The research suggests a material design approach that accounts for temperature effects to enhance the durability and performance of nuclear fusion materials.

This manuscript embarks on an inquiry into the influence of helium implantation on nanocrystalline tungsten, a contender for plasma-facing components (PFCs) in nuclear fusion reactors. The study underscores the inevitability of helium retention in tungsten due to the anticipated high flux of helium atoms from D-T fusion reactions in future reactors such as ITER and DEMO. This retention, potentially culminating in surface blistering and grain boundary embrittlement, propels an investigation into the helium's preferential substitutional configuration over interstitial within tungsten's lattice, leading to cavity and bubble formation. Atomistic tensile testing simulations, exploiting a helium-specific interatomic potential, dissect the ramifications of helium concentration and voids on tungsten's mechanical properties across grain sizes of 10nm and 15nm. The presence of helium at grain boundaries instigates detachment, influencing structural behaviors under strain. Notably, the study reveals a dichotomy in helium's effect: benign at concentrations up to the critical point, beyond which significant embrittlement and elastic softening occur, corroborating with theoretical predictions on helium-induced grain boundary embrittlement. This nuanced exploration delineates the intricate dance between helium implantation levels, grain size, and the ensuing mechanical property alterations in nanocrystalline tungsten, casting light on its viability and endurance as a PFC material under the rigorous conditions forecasted in fusion reactors.

This study explores the intricate interactions between grain boundaries (GBs) and irradiation-induced defects in nanocrystalline iron, highlighting the role of dopants like copper. Utilizing molecular dynamics simulations, the research delineates how GB properties, such as GB energy and defect formation energies, influence the formation and evolution of defects in low energy collision cascades. It reveals that GBs not only augment defect production but also show a marked preference for interstitials over vacancies, a behavior significantly modulated by the cascade's proximity to the GB. The presence of dopants is shown to alter GB properties, affecting both the rate and type of defect production, thereby underscoring the complex interplay between GB characteristics, dopant elements, and defect dynamics. Moreover, the investigation uncovers that the structural characteristics of GBs play a crucial role in cascade evolution and defect generation, with certain GB configurations undergoing reconfiguration in response to cascades. For instance, the reconfiguration of one pure Fe twist GB suggests that GB geometry can significantly influence defect generation mechanisms. These findings point to the potential of GB engineering in developing materials with enhanced radiation tolerance, advocating for a nuanced approach to material design. By tailoring GB properties and selectively introducing dopant elements, materials can be optimized to exhibit superior resistance to radiation-induced damage, offering insights for applications in nuclear reactors and other radiation-prone environments.

Constitutive models are essential for describing the mechanical behavior of materials under high temperatures and pressures, among which the Steinberg-Guinan (SG) model is widely adopted. Recent work has discovered a peculiar dual anomaly of compression-induced softening and heating-induced hardening in the elasticity of compressed vanadium [Phys. Rev. B 104, 134102 (2021)], which is beyond the capability of the SG model to describe. In this work, a modified SG constitutive model is proposed to embody such an anomalous behavior. Elemental vanadium is considered as an example to demonstrate the effectiveness of this improved model in describing the dual anomalies of mechanical elasticity. This new variant of the SG model can also be applied to other materials that present an irregular variation in the mechanical elasticity, and is important to faithfully model and simulate the mechanical response of materials under extreme conditions., Comment: 14 pages, 4 figures

First-principle calculations were employed to analyze the effects induced by vacancies of molybdenum (Mo) and sulfur (S) on the dielectric properties of few-layered MoS2. We explored the combined effects of vacancies and dipole interactions on the dielectric properties of few-layered MoS2. In the presence of dielectric screening, we investigated uniformly distributed Mo and S vacancies, and then considered the case of concentrated vacancies. Our results show that the dielectric screening remarkably depends on the distribution of vacancies owing to the polarization induced by the vacancies and on the interlayer distances. This conclusion was validated for a wide range of wide-gap semiconductors with different positions and distributions of vacancies, providing an effective and reliable method for calculating and predicting electrostatic screening of dimensionally reduced materials. We further provided a method for engineering the dielectric constant by changing the interlayer distance, tuning the number of vacancies and the distribution of vacancies in few-layered van der Waals materials for their application in nanodevices and supercapacitors.

Deep-learning density functional theory (DFT) shows great promise to significantly accelerate material discovery and potentially revolutionize materials research, which demands a close combination between neural networks and DFT computation. However, current research in this field primarily relies on supervised learning, making the developments of neural networks and DFT isolated from each other. In this work, we present a theoretical framework of neural-network DFT, which unifies the optimization of neural networks with the variational computation of DFT, enabling physics-informed unsupervised learning. Moreover, we develop a differential DFT code incorporated with deep-learning DFT Hamiltonian, and introduce algorithms of automatic differentiation and backpropagation to DFT, demonstrating the concept of neural-network DFT. The advanced neural-network architecture not only surpasses conventional approaches in accuracy and efficiency, but offers a new paradigm for developing deep-learning DFT methods.

Recent advances in deep learning generative models (GMs) have created unprecedented high capabilities in accessing and assessing complex high-dimensional data, allowing superior efficiency in navigating vast material configuration space in search of viable structures. Coupling such capabilities with physically significant data to construct trained models for materials discovery is crucial to moving this emerging field forward. Here, we present a universal GM for crystal structure prediction (CSP) via a conditional crystal diffusion variational autoencoder (Cond-CDVAE) approach, which is tailored to allow user-defined material and physical parameters such as composition and pressure. This model is trained on an expansive dataset containing over 670,000 local minimum structures, including a rich spectrum of high-pressure structures, along with ambient-pressure structures in Materials Project database. We demonstrate that the Cond-CDVAE model can generate physically plausible structures with high fidelity under diverse pressure conditions without necessitating local optimization, accurately predicting 59.3% of the 3,547 unseen ambient-pressure experimental structures within 800 structure samplings, with the accuracy rate climbing to 83.2% for structures comprising fewer than 20 atoms per unit cell. These results meet or exceed those achieved via conventional CSP methods based on global optimization. The present findings showcase substantial potential of GMs in the realm of CSP.