Search

In this work, methods are presented to automatically generate a fully atomistic LAMMPS models of arbitrary linear multiblock polyurethane copolymers. The routine detailed here receives as parameters the number of repeat units per hard block, the number of units in a soft block, and the number of soft blocks per chain, as well as chemical formulae of three monomers which will form the hard component, soft component, and chain extender. A routine is detailed for converting the chemical structure of a free monomer to the urethane bonded repeat units in a polymer. The python package RadonPy is leveraged to assemble these units into blocks, and the blocks into copolymers. Care is taken in this work to ensure that plausible atomic charges are assigned to repeat units in different parts of the chain. The static structure factor is calculated for a variety of chemistries, and the results compared with wide angle x-ray scattering data from experiments with corresponding composition. The generated models reproduce the amorphous halo observed in the scattering data as well as some of the finer details. Structure factor calculations are decomposed into the partial structure factors to interrogate the structural properties of the two block types separately. Parametric surveys are carried out of the effects of various parameters, including temperature, soft block length, and block connectivity on the observed structure. The routine detailed here for constructing models is robust enough to be executed automatically in a high throughput workflow for material design and discovery., Comment: 6 figures, 18 pages

Atomic-scale simulations have progressed tremendously over the past decade, largely due to the availability of interatomic potentials based on machine-learning algorithms. These potentials enable the combination of the accuracy of electronic structure calculations with extensive length and time scales. In this paper, we present a new release of the i-PI code that allows the community to fully benefit from the rapid developments in the field. i-PI is a Python software that facilitates the integration of different methods and different software tools by using a socket interface for inter-process communication. The current framework is flexible enough to support rapid prototyping and the combination of various simulation techniques, while maintaining a speed that prevents it from becoming the bottleneck in most workflows. We discuss the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities. For this release, we also improved some computational bottlenecks of the implementation, reducing the overhead associated with using i-PI over a native implementation of molecular dynamics techniques. We show numerical benchmarks using widely adopted machine learning potentials, such as Behler-Parinello, DeepMD and MACE neural networks, and demonstrate that such overhead is negligible for systems containing between 100 and 12000 atoms.

Topological insulators have been studied intensively over the last decades. Among these materials, three-dimensional (3D) zirconium pentatelluride (ZrTe$_5$) stands out as one of the most intriguing for both theoretical and experimental studies because of its diverse range of distinct topological phases. In this work, we employ density functional theory to study the electronic structure and quantum oscillations exhibited by various topological phases of 3D bulk ZrTe$_5$. We have discovered that by analyzing combined patterns in band structures, Fermi surfaces, and Shubnikov-de Haas (SdH) oscillations we can determine the corresponding topological phase without relying on the conventional calculation of topological invariants or boundary state contributions. This approach facilitates the identification of topological phases in ZrTe$_5$ directly from experimental quantum oscillation measurements. Using this method, we have analyzed the entire process of topological phase transition, revealing changes in the topology of the Fermi pockets and validating the shapes deduced from the experimental data for the topological phases.

The manifold of energetically degenerate configurations arising from competing interactions in frustrated magnets gives rise to an enhanced entropy at lowest temperatures, which can be utilized for adiabatic demagnetization refrigeration (ADR). We review structural and magnetic properties of various Yb- and Gd-based oxides featuring frustration related to different triangular moment configurations and (in some cases) structural randomness. In comparison to paramagnetic hydrated salts, which have traditionally been employed for mK-ADR, these novel ADR materials enable cooling to temperatures several times lower than the magnetic interaction strength, significantly enhancing the entropy density and cooling power at a given target temperature. A further advantage is their chemical stability, allowing for a much simpler ADR pill design and ultra-high vacuum applications. For the temperature range between 0.02 and 2 K, a systematic comparison of the field-induced entropy density change is provided, that illustrates the advantages of frustrated magnets for low-temperature ADR., Comment: 17 pages, comments and suggestions are welcome

Probing the quantum geometry and topology in condensed matter systems has relied heavily on static electronic transport experiments in magnetic fields. Yet, contact-free optical measurements have rarely been explored. Magnetic dichroism (MCD), the nonreciprocal absorption of circular polarized light, was theoretically linked to the quantized anomalous Hall effect in magnetic insulators and can identify the bands and momenta responsible for the underlying Berry Curvature (BC). Detecting BC through MCD faces two challenges: First, the relevant inter-band transitions usually generate MCD in the infrared (IR) range, requiring large samples with high quality. Second, while most magnetic materials are metallic, the relation between MCD and BC in metals remains unclear. Here, we report the observation of MCD in the IR range along with the anomalous Hall effect in thin film MnBi2Te4. Both phenomena emerge with a field-driven phase transition from an antiferromagnet to a canted ferromagnet. By theoretically relating the MCD to the anomalous Hall effect via BC in a metal, we show that this transition accompanies an abrupt onset of BC, signaling a topological phase transition from a topological insulator to a doped Chern insulator. Our density functional theory calculation suggests the MCD signal mainly originates from an optical transition at the Brillouin zone edge, hinting at a potential new source of BC away from the commonly considered {\Gamma} point. Our findings demonstrate a novel experimental approach for detecting BC and identifying the responsible bands and momenta, generally applicable to magnetic materials.

Realizing topological superconductivity by integrating high-transition-temperature ($T_C$) superconductors with topological insulators can open new paths for quantum computing applications. Here, we report a new approach for increasing the superconducting transition temperature ($T_{C}^{onset}$) by interfacing the unconventional superconductor Fe(Te,Se) with the topological insulator Bi-Te system in the low-Se doping regime, near where superconductivity vanishes in the bulk. The critical finding is that the $T_{C}^{onset}$ of Fe(Te,Se) increases from nominally non-superconducting to as high as 12.5 K when $Bi_2Te_3$ is replaced with the topological phase $Bi_4Te_3$. Interfacing Fe(Te,Se) with $Bi_4Te_3$ is also found to be critical for stabilizing superconductivity in monolayer films where $T_{C}^{onset}$ can be as high as 6 K. Measurements of the electronic and crystalline structure of the $Bi_4Te_3$ layer reveal that a large electron transfer, epitaxial strain, and novel chemical reduction processes are critical factors for the enhancement of superconductivity. This novel route for enhancing $T_C$ in an important epitaxial system provides new insight on the nature of interfacial superconductivity and a platform to identify and utilize new electronic phases., Comment: 18 pages, 5 figures, accepted by Advanced Materials

Three solid solutions of [CH$_3$NH$_3$]Co$_x$Ni$_{1-x}$(HCOO)$_3$, with $x$ = 0.25 (1), 0.50 (2) and 0.75 (3), were synthesized and their nuclear structures and magnetic properties were characterized using single crystal neutron diffraction and magnetization measurements. At room temperature, all three compounds crystallize in the Pnma orthorhombic space group, akin to the cobalt and nickel end series members. Upon cooling, each compound undergoes distinct series of structural transitions to modulated structures. Compound 1 exhibits a phase transition to a modulated structure analogous to the pure nickel compound, while compound 3 maintains the behaviour observed in the pure cobalt compound, although in both cases the temperatures at which the phase transitions occur differ slightly from the pure phases. Monochromatic neutron diffraction measurements showed that the structural evolution of 2 diverges from that of either parent compound, with the competing hydrogen bond interactions which drive the modulation throughout the series producing a unique sequence of phases. It involves two modulated phases below 96(3) K and 59(3) K, with different q vectors, similar to the pure cobalt compound (with modulated phases below 128 K and 96 K), however it maintains the modulated phase below magnetic order (at 22.5(7) K), resembling the pure nickel compound (which present magnetic order below 34 K), resulting in an improper modulated magnetic structure. Despite these large scale structural changes, magnetometry data reveal that the bulk magnetic properties of these solid solutions form a linear continuum between the end members. Notably, doping of the metal site in these solid solutions allows for tuning of bulk magnetic properties, including magnetic ordering temperature, transition temperatures, and the nature of nuclear phase transitions, through adjustment of metal ratios., Comment: 21 pages (31 pages inc SI), 7 Figures, 9 SI Figures

The manipulation of magnetization via Magnetic torques is one of the most important phenomena in spintronics. In thin films, conventionally, a charge current flowing in a heavy metal is used to generate transverse spin currents and to exert torques on the magnetization of an adjacent ferromagnetic thin film layer. Here, in contrast to the typically employed heavy metals, we study spin-to-charge conversion in ferromagnetic heterostructures with large spin-orbit interaction that function as the torque-generating layers. In particular, we chose perpendicular magnetic anisotropy (PMA) multilayers [Co/Ni] and [Co/Pt] as the torque-generating layers and drive magnetization dynamics in metallic ferromagnetic thin film $\mathrm{Co_{20}Fe_{60}B_{20}}$ (CoFeB) layers with in-plane magnetic anisotropy (IMA). We investigate the spin dynamics driven by spin-orbit torque (SOT) and the concomitant charge current generation by the inverse SOT process using an inductive technique based on a vector network analyzer. In our experimental findings, we find that the SOTs generated by our multilayers are of a magnitude comparable to those produced by Pt, consistent with first-principles calculations. Furthermore, we noted a significant correlation between the SOT and the thickness of the CoFeB layer.

Lablets are autonomous microscopic particles with programmable CMOS electronics that can control electrokinetic phenomena and electrochemical reactions in solution via actuator and sensor microelectrodes. In this paper, we describe the design and fabrication of optimized singulated lablets (CMOS3) with dimensions 140x140x50 um carrying an integrated coplanar encapsulated supercapacitor as a rechargeable power supply. The lablets are designed to allow docking to one another or to a smart surface for interchange of energy, electronic information, and chemicals. The paper focusses on the digital and analog design of the lablets to allow significant programmable functionality in a microscopic footprint, including the control of autonomous actuation and sensing up to the level of being able to support a complete lablet self-reproduction life cycle, although experimentally this remains to be proven. The potential of lablets in autonomous sensing and control and for evolutionary experimentation are discussed., Comment: This article was originally submitted (2016) for review as supporting information for the final review of the EU MICREAgents Project # 318671 (2012-2016). Here it is presented in slightly revised form

The layered van der Waals material CrCl3 exhibits very strongly bound ligand field excitons that control optoelectronic applications and are connected with magnetic ordering by virtue of their d-orbital origin. Time-resolved photoluminescence of these exciton populations at room temperature shows that their relaxation is dominated by exciton-exciton annihilation and that the spontaneous decay lifetime is very long. These observations allow the rough quantification of the exciton annihilation rate constant and contextualization in light of a recent theory of universal scaling behavior of the annihilation process.

We extend the electron localization function -- a popular descriptor of molecular bonds and atomic shells -- to general non-collinearly magnetized states, open-shell solutions of the Pauli-Schr{\"o}dinger equation. The description is also improved in paradigmatic, collinear states., Comment: 4 figures

Solving the Schr\"{o}dinger equation for interacting many-body quantum systems faces computational challenges due to exponential scaling with system size. This complexity limits the study of important phenomena in materials science and physics. We develop an Artificial Neural Network (ANN)-driven algorithm to simulate fermionic systems on lattices. Our method uses Pauli matrices to represent quantum states, incorporates Markov Chain Monte Carlo sampling, and leverages an adaptive momentum optimizer. We demonstrate the algorithm's accuracy by simulating the Heisenberg Hamiltonian on a one-dimensional lattice, achieving results with an error in the order of $10^{-4}$ compared to exact diagonalization. Furthermore, we successfully model a magnetic phase transition in a two-dimensional lattice under an applied magnetic field. Importantly, our approach avoids the sign problem common to traditional Fermionic Monte Carlo methods, enabling the investigation of frustrated systems. This work demonstrates the potential of ANN-based algorithms for efficient simulation of complex quantum systems, opening avenues for discoveries in condensed matter physics and materials science.

Functionalities of solid-state materials are usually considered to be dependent on their crystal structures. The limited structural types observed in the emerged high-entropy oxides put constraints on exploration of their physical properties and potential applications. Herein, we synthesized the first high-entropy oxide in a trirutile structure, (Mn0.2Fe0.2Co0.2Ni0.2Cu0.2)Ta2O6, and investigated its magnetism. The phase purity and high-entropy nature were confirmed by powder X-ray diffraction and energy-dispersive spectroscopy, respectively. X-ray photoelectron spectroscopy indicated divalent Mn, Co and Cu along with coexistence of divalent and trivalent Fe and Ni. Magnetic properties measurements showed antiferromagnetic coupling and potential short-range magnetic ordering below ~ 4 K. The temperature-dependent heat capacity data measured under zero and high magnetic field confirmed the short-range magnetic ordering and a possible low-temperature phonon excitation. Moreover, the possible reasons for the absence of long-range magnetic ordering are concluded from a molecular orbital perspective. The discovery of the first trirutile high-entropy oxide opens a new way for studying the relationship between the highly disordered atomic arrangement and their magnetic interaction. Furthermore, it provides a new direction for exploring functionalities of high-entropy oxides., Comment: 25 pages, 11 figures, 2 tables

The simulation of quantum many-body systems poses a significant challenge in physics due to the exponential scaling of Hilbert space with the number of particles. Traditional methods often struggle with large system sizes and frustrated lattices. In this research article, we present a novel algorithm that leverages the power of deep neural networks combined with Markov Chain Monte Carlo simulation to address these limitations. Our method introduces a neural network architecture specifically designed to represent bosonic quantum states on a 1D lattice chain. We successfully achieve the ground state of the Bose-Hubbard model, demonstrating the superiority of the adaptive momentum optimizer for convergence speed and stability. Notably, our approach offers flexibility in simulating various lattice geometries and potentially larger system sizes, making it a valuable tool for exploring complex quantum phenomena. This work represents a substantial advancement in the field of quantum simulation, opening new possibilities for investigating previously challenging systems.

1T-VSe$_2$ is a narrow band transition metal chalcogenide that shows charge density wave (CDW) state below $T_{CDW}$ = 110 K. Here, we have explored the relevance of Kohler's rule and the thermal transport properties of VSe$_2$ across the CDW state. The magnetoresistance (MR) follows Kohler's rule above $T_{CDW}$, while an extended Kohler's rule is employed below $T_{CDW}$. Interestingly, we observed an anomaly in MR at T = 20 K, below which MR value decreases on lowering temperature. This anomaly is also reflected in the slope ($\kappa$) of Kohler's plots and the relative change in the thermal excitation induced carrier density ($n_T$) also. The $T_{CDW}$ remains largely unaffected in both electrical resistivity ($\rho(T)$) and longitudinal Seebeck coefficient ($\it{S_{xx}}$) even under a strong magnetic field of 14 Tesla. However, the application of magnetic field enhances the peak intensity of $\it{S_{xx}}$ at T $\sim$ 60 K. Additionally, $\it{S_{xx}(T)}$ curves measured at different fields exhibit a crossover at T = 20 K, which suggest the existence of unique feature in the CDW state of VSe$_2$ \textit{i.e.} a locally exact isosbestic point., Comment: 10 pages, 5 figures

The structure of evaporated amorphous Ge$_x$Sb$_x$Te$_{100-2x}$ ($x =$ 6, 9, 13) alloys was investigated by neutron diffraction, X-ray diffraction and extended X-ray absorption spectroscopy (EXAFS) at the Ge, Sb and Te K-edges. Large scale structural models were generated by fitting the experimental datasets (5 for each composition) simultaneously in the framework of the reverse Monte Carlo simulation technique. It was found that the alloys are chemically ordered (Ge and Sb have predominantly Te neighbors) and within the experimental uncertainty each component satisfies the $8-N$ rule. A comparison with the pair correlation functions of melt quenched Ge$_{20}$Te$_{80}$ revealed that the first minimum of $g_{\mathrm{TeTe}}(r)$ is shallower in the ternary alloys than in Ge$_{20}$Te$_{80}$. On the other hand, the separation of the first and second coordination environments of Ge atoms is stronger in the Ge-Sb-Te alloys investigated., Comment: 32 pages, submitted to Journal of Alloys and Compounds

Exploring new Hall effect is always a fascinating research topic. The ordinary Hall effect and the quantum Hall effect, initially discovered in two-dimensional (2D) non-magnetic systems, are the phenomena that a transverse current is generated when a system carrying an electron current is placed in a magnetic field perpendicular to the currents. In this work, we propose the electric counterparts of these two Hall effects, termed as electric Hall effect (EHE) and quantum electric Hall effect (QEHE). The EHE and QEHE emerge in 2D magnetic systems, where the transverse current is generated by applying an electric gate-field instead of a magnetic field. We present a symmetry requirement for intrinsic EHE and QEHE. With a weak gate-field, we establish an analytical expression of the intrinsic EHE coefficient. We show that it is determined by intrinsic band geometric quantities: Berry curvature and its polarizability which consists of both intraband and interband layer polarization. Via first-principles calculations, we investigate the EHE in the monolayer Ca(FeN)$_2$, where significant EHE coefficient is observed around band crossings. Furthermore, we demonstrate that the QEHE can appear in the semiconductor monolayer $\rm BaMn_2S_3$, of which the Hall conductivity exhibits steps that take on the quantized values $0$ and $\pm1$ in the unit of $e^2/h$ by varying the gate-field within the experimentally achievable range. Due to the great tunability of the electric gate-field, the EHE and QEHE proposed here can be easily controlled and should have more potential applications., Comment: 7 pages, 3 figures

In this letter, we derive an expression for magnetic dipole transitions that is analogous to the Lyddane-Sachs-Teller relation for dielectric polar lattice vibrations. We thereby define transverse and longitudinal optical frequencies at which paramagnetic resonance and antiresonance occurs, respectively. The relation found here thus permits non-invasive optical analysis of static magnetization properties in paramagnetic materials. Validated through terahertz electron paramagnetic resonance ellipsometry and superconducting quantum interference device measurements on Iron-doped Gallium Nitride, our findings show very good agreement between theory and experiment. We term the excitations associated with the paramagnetic transitions as paramagnetic polaritons which may find use in future photonic applications

Among the 3d transition metal ions doped polar Ni3TeO6, Mn-doped Ni3TeO6 has stimulated great interest due to its high magnetic ordering temperature and complex magnetic phases, but the mechanism of magnetoelectric (ME) coupling is far from understood. Herein we report our systematic investigation of the chemical control of magnetism, metamagnetic transition, and ME properties of Ni3-xMnxTeO6 single crystals in high magnetic field (H) up to 52 T. We present a previously unreported weak ferromagnetic behavior appeared in the ab plane below 9.5 K in addition to the incommensurate helical and commensurate collinear antiferromagnetic states. In the low-field region, a spin-flop type metamagnetic transition without any hysteresis occurs at Hc1 for H // c, while another metamagnetic transition accompanied with a change in electric polarization is observed at Hc2 in the high-field region both for H // c and H // ab above 30 K, which can be attributed to the sudden rotation of magnetic moments at Ni2 sites. The ME measurements reveal that a first-order ME effect is observed in the low-T and low-H regions, while a second-order ME coupling term appears above 30 K in the magnetic field range of Hc1 < H < Hc2 for H // c and H < Hc2 for H // ab, both becoming significant with increasing temperature. Eventually, they are dominated by the second-order ME effect near the antiferromagnetic transition temperature. The present work demonstrates that Ni3-xMnxTeO6 is an exotic magnetoelectric material compared with Ni3TeO6 and its derivatives, thereby providing insights to better understand the magnetism and ME coupling in Ni3TeO6 and its derivatives., Comment: 30 pages with 8 figures

Scientific text is a promising source of data in materials science, with ongoing research into utilising textual data for materials discovery. In this study, we developed and tested a novel approach to extract material-property relationships from scientific publications using the Question Answering (QA) method. QA performance was evaluated for information extraction of perovskite bandgaps based on a human query. We observed considerable variation in results with five different large language models fine-tuned for the QA task. Best extraction accuracy was achieved with the QA MatBERT and F1-scores improved on the current state-of-the-art. This work demonstrates the QA workflow and paves the way towards further applications. The simplicity, versatility and accuracy of the QA approach all point to its considerable potential for text-driven discoveries in materials research., Comment: The following article has been submitted to npj Computational Materials

The van der Waals material MnBi$_2$Te$_4$ and the related Sb-substituted compounds Mn(Bi$_{1-x}$Sb$_x$)$_2$Te$_4$ are prominent members of the family of magnetic topological insulators, in which rare quantum mechanical states can be realized. In this work, we study the evolution of the charge dynamics in MnBi$_2$Te$_4$ and the Sb-substituted compound Mn(Bi$_{1-x}$Sb$_x$)$_2$Te$_4$ with $x=0.26$ under hydrostatic pressure. For MnBi$_2$Te$_4$, the pressure dependence of the screened plasma frequency, the dc conductivity, and the reflectance at selected frequencies shows weak anomalies at $\sim$2 and $\sim$4~GPa, which might be related to an electronic phase transition driven by the enhanced interlayer interaction. We observe a pressure-induced decrease in the optical gap, consistent with the decrease in and closing of the energy gap reported in the literature. Both studied materials show an unusual decrease in the low-energy optical conductivity under pressure, which we attribute to a decreasing spectral weight of the Drude terms describing the free charge carrier excitations. Our results suggest a localization of conduction electrons under pressure, possibly due to hybridization effects., Comment: 9 pages, 3 figures; accepted for publication in Phys. Rev. B

While all of the polymorphs of NH$_4$I and NH$_4$Br are non-polar, a reversible electric polarization is established in the ordered $\gamma$ phases of (NH$_4$)$_{0.73}$(ND$_4$)$_{0.27}$I and (NH$_4$)$_{0.84}$(ND$_4$)$_{0.16}$Br (where D is $^2$H) via $dc$ electric fields. The presence of two groups of orbital magnetic moments appears to be responsible for the asymmetric lattice distortions. Our findings provide an alternative pathway for hydrogen-based materials to potentially add a ferroelectric functionality., Comment: 14 pages, 3 figures

While all of the polymorphs of pure NH$_4$I and KI are non-polar, we identify that (NH$_4$)$_{0.95}$K$_{0.05}$I is ferroelectric and (NH$_4$)$_{0.87}$K$_{0.13}$I and (NH$_4$)$_{0.83}$K$_{0.17}$I are pyroelectric through measurements of their pyroelectric current and complex dielectric constant. The order to disorder phase transitions occur near 245 K. Magnetic susceptibility measurements indicate that the proton orbitals of the NH$_4$$^+$ continue to become ordered in the ground state in the (NH$_4$)$_{1-x}$K$_x$I system up to x <= 0.17. The polar phases are proposed to stem from K$^+$ ions disrupting the symmetry of proton-orbital-lattice interactions between the NH$_4$$^+$ and I$^-$ ions. Our work introduces a new pathway for the ordered phases of ammonium-based compounds to potentially become ferroelectric., Comment: 13 pages, 5 figures

High-accuracy microwave sensing is widely demanded in various fields, ranging from cosmology to microwave quantum technology. Quantum receivers based on inorganic solid-state spin systems are promising candidates for such purpose because of the stability and compatibility, but their best sensitivity is currently limited to a few pT/$\sqrt{\rm{Hz}}$. Here, by utilising an enhanced readout scheme with the state-of-the-art solid-state maser technology, we develop a robust microwave quantum receiver functioned by organic molecular spins at ambient conditions. Owing to the maser amplification, the sensitivity of the receiver achieves 6.14 $\pm$ 0.17 fT/$\sqrt{\rm{Hz}}$ which exceeds three orders of magnitude than that of the inorganic solid-state quantum receivers. The heterodyne detection without additional local oscillators improves bandwidth of the receiver and allows frequency detection. The scheme can be extended to other solid-state spin systems without complicated control pulses and thus enables practical applications such as electron spin resonance spectroscopy, dark matter searches, and astronomical observations., Comment: 10 pages, 4 figures

The thickness of 2D materials not only plays a crucial role in determining the performance of nanoelectronic and optoelectronic devices but also introduces complexities in predicting volume-dependent properties such as energy storage capacity, due to the intrinsic vacuum within these materials. Although a plethora of experimental techniques, including but not limited to optical contrast, Raman spectroscopy, nonlinear optical spectroscopy, near-field optical imaging, and hyperspectral imaging, facilitate the measurement of 2D material thickness, comprehensive data for many materials remains elusive. Over the last decade, the exponential proliferation of 2D materials and their heterostructures has outstripped the capabilities of conventional experimental and computational approaches. In this evolving landscape, machine learning (ML) has emerged as an indispensable tool, offering novel avenues to augment these traditional methodologies. Addressing the critical gap, we introduce THICK2D - Thickness Hierarchy Inference and Calculation Kit for 2D Materials. This Python-based computational framework harnesses an autogenerated thickness database, developed using large language models (LLMs), and advanced ML algorithms to facilitate the rapid and scalable estimation of material thickness, relying solely on crystallographic data. To demonstrate the utility and robustness of THICK2D, we successfully employed the toolkit to predict the thickness of more than 8000 2D-based materials, sourced from two extensive 2D material databases. THICK2D is disseminated as an open-source utility, accessible on GitHub https://github.com/gmp007/THICK2D, and archived on Zenodo at https://doi.org/10.5281/zenodo.11216648}{10.5281/zenodo.11216648.

Moire superlattices formed in van der Waals heterostructures due to twisting, lattice mismatch and strain present an opportunity for creating novel metamaterials with unique properties not present in the individual layers themselves. Ferroelectricity for example, arises due to broken inversion symmetry in twisted and strained bilayers of 2D semiconductors with stacking domains of alternating out-of-plane polarization. However, understanding the individual contributions of twist and strain to the formation of topological polar nanostructures remains to be established and has proven to be experimentally challenging. Inversion symmetry breaking has been predicted to give rise to an in-plane component of polarization along the domain walls, leading to the formation of a network of topologically non-trivial merons (half-skyrmions) that are Bloch-type for twisted and Neel-type for strained systems. Here we utilise angle-resolved, high-resolution vector piezoresponse force microscopy (PFM) to spatially resolve polarization components and topological polar nanostructures in marginally twisted bilayer WSe2, and provide experimental proof for the existence of topologically non-trivial meron/antimeron structures. We observe both Bloch-type and Neel-type merons, allowing us to differentiate between moire superlattices formed due to twist or heterogeneous strain. This first demonstration of non-trivial real-space topology in a twisted van der Waals heterostructure opens pathways for exploring the connection between twist and topology in engineered nano-devices.

Medium and high entropy ceramics, with their distinctive disordered structures, exhibit ultra-low thermal conductivity and high temperature stability. These properties make them strong contenders for next generation thermal barrier coating (TBC) materials. However, predicting their thermal conductivity has been challenging, primarily due to their unique phonon scattering mechanisms. Apart from the conventional phonon-phonon scattering mechanism, the phonon-disorder scattering, comprising both mass and force constant disorder, are also expected to make significant contribution in determining the thermal conductivity of medium and high entropy ceramics. However, it remains challenging to quantify the phonon-disorder contribution, particular in the aspect of force constant disorder. Here we demonstrated a relationship between the lattice disorder, a quantity more readily calculable, with force constant disorder, rendering it possible to substitute the force constant disorder by lattice disorder. Based on this relationship and drawing inspiration from Klement's equation of static imperfection, we have developed a model that quantitatively assesses the connection between disorder and thermal conductivity. Applying our model to the medium/high entropy rocksalt and pyrochlore oxides as representatives, we found good alignment between the theoretical predictions and experimental measurements of thermal conductivities, confirming the validity of our model. The model developed offers a critical predictive tool for rapid screening of TBC materials based on medium and high entropy ceramics., Comment: 25 pages, 7 figures

Semiconducting transition metal dichalcogenides (MX$_2$) occur in 2H and rhombohedral polytypes, respectively distinguished by anti-parallel and parallel orientation of consecutive monolayer lattices. In its bulk form, rhombohedral MX$_2$ is ferroelectric, hosting an out-of-plane electric polarization, the direction of which is dictated by stacking. Here, we predict that twin boundaries, separating adjacent polarization domains with reversed built-in electric fields, are able to host two-dimensional electrons and holes with an areal density reaching $\sim 10^{13} {\rm cm}^{-2}$. Our modelling suggests that n-doped twin boundaries have a more promising binding energy than p-doped ones, whereas hole accumulation is stable at external surfaces of a twinned film. We also propose that the introduction of a twist at the twinned interface, with a `magic' angle $\theta^* \approx 3^{\circ}$, should promote a regime of strongly correlated states of electrons, such as Wigner crystals, due to commensurability between the moir\'e pattern at the interface and the accumulated carrier density., Comment: 10 pages, 7 Figures

We investigate the dynamics of two-dimensional assemblies of rod-shaped magnetic colloids under the influence of an external rotating magnetic field. Using Molecular Dynamics, we simulate the formation of patterns that emerge based on the synchronization degree between the magnetic rods and the rotating field. We then explore the structural and dynamic characteristics of the resulting steady states, examining their evolution as a function of changes in the rods' aspect ratio, the strength of the external magnetic field, and its rotation frequency. Three distinct synchronization regimes of the rods with the magnetic field are clearly observed. A detailed set of phase diagrams illustrates the complex relationship between the magnitude of the external magnetic field and its rotation frequency and how these parameters govern the formation of unique self-organized structures., Comment: 10 pages, 9 figures

CrI3 is a layered ferromagnetic insulator that has recently attracted enormous interest as it was the first example of a stand-alone monolayer ferromagnet, paving the way towards the study of two-dimensional magnetic materials and their use as building blocks of hybrid van der Waals layered heterostructures. Here we go one step down in the dimensionality ladder and report the synthesis and characterization of a tubular one-dimensional van der Waals heterostructure where CrI3 nanotubes are encapsulated within multiwall carbon nanotubes, integrating a magnetic insulator and a conductor. By means of the capillary filling of multi-wall carbon nanotubes (MWCNT), we obtained single-wall CrI3 nanotubes with diameters ranging between 2 nm and 10 nm, with an average of 5.3 nm. Using aberration corrected electron microscopy in combination with spectroscopic techniques we confirm the structure and chemical composition of the nanotubes. SQUID measurements, combined with element-specific X-ray magnetic circular dichroism (XMCD) indicate unequivocally that the Cr atoms in encapsulated CrI3 nanotubes are magnetic with a collective state compatible with a radial magnetization state predicted both by first-principles calculations and a model Hamiltonian. Our results represent a step forward in establishing 1D van der Waals heterostructures as a playground for the exploration of non-collinear magnetic states arising from the interplay between magnetic anisotropy and curvature in tubular geometries., Comment: 19 pages with 5 figures

In nuclear and/or solar applications that involve molten salts, impurities frequently enter the salt as either fission products or via corrosion. Impurities can interact and make complexes, but the impact of such complexation on the properties of the salts and corrosion rates has not been understood. Common impurities in molten salts, such as FLiBe, include Cr, Ni, and Fe. Here, we investigate the complexation of Cr and Ni in FLiBe using molecular dynamics based on a machine learning interatomic potential (MLIP) fitted using the atomic cluster expansion (ACE) method. The MLIP allows us to overcome the challenges of simultaneously needing accurate energetics and long time scale to study complexation. We demonstrate that impurity behavior is more difficult to capture than that of concentrated elements with MLIPs due to less sampling in training data, but that this can be overcome by using active learning strategies to obtain a robust fit. Our findings suggest that there is a weak but potentially significant binding free energy between CrF2 and NiF2 in eutectic FLiBe of -0.112 eV. Under certain conditions this binding creates a significant short-range order between the species and lowers the redox potential of NiF2 in the presence of CrF2 in FLiBe, making Ni dissolution more favorable in the presence of Cr as compared to its dissolution in pure FLiBe. However, we find little impact of this complexation on the diffusivity of Ni and Cr. Overall, the methodology presented here suggests an approach to modeling complexation with MLIPs and suggests that interactions between dissolved cations could be playing a significant role in some salt thermophysical properties.

We propose a method to evaluate the Gibbs free energy from constant-volume first-principles calculations. The volume integral of the pressure is performed by determining the volume and the bulk modulus in equilibrium at finite temperatures, where the pressure and its volume derivative are evaluated utilizing first-principles calculations of the Gr\"uneisen parameter without varying the volume. As an example, the validity of our method is demonstrated for fcc-Al by comparing with the conventional quasiharmonic approximation that is much more computationally-demanding.

Hydrogen has been identified as a clean, zero carbon, sustainable, and promising energy source for the future, and electrochemical water splitting for hydrogen production is an emission-free, efficient energy conversion technology. A major limitation of this approach is the unavailability of efficient, abundant, inexpensive catalysts, which prompts the need for new catalytic materials. Here, we report the synthesis and electrocatalytic properties of a novel transition metal-based ternary chalcogenide family, LaMS$_3$ (M = Mn, Fe, Co, Ni). Powder X-ray diffraction confirms the phase purity of these materials, while composition analysis using energy dispersive spectroscopy (EDS) confirms the presence of the stoichiometric ratio of elements in these compounds. X-ray photoelectron spectroscopy (XPS) and X-ray absorption spectroscopy (XAS) were used to study the chemical states on the surface and in bulk, respectively. These materials exhibit bifunctional catalytic activity towards the two half-reactions of the water-splitting process, with LaNiS$_3$ being the most active material for both hydrogen evolution reaction (HER) and oxygen evolution reaction (OER). The LaMS$_3$ compounds show long-term stability with negligible change in the overpotential at a constant current density of 10 mA cm$^{-2}$ over 18 hours of measurements. As compared to the corresponding ternary oxides, the LaMS$_3$ materials exhibit higher activity and significantly lower Tafel slopes. The ability to catalyze both half-reactions of water electrolysis makes these materials promising candidates for bifunctional catalysts and presents a new avenue to search for high-efficiency electrocatalysts for water splitting.

Plutonium's phase diagram is host to complex structures and interactions that make the description of its ground state properties elusive. Using all-electron density functional theory, we study the thermodynamic properties of $\alpha$-Pu. To do this, we build on recent work in the literature by introducing a novel noncollinear magnetic ansatz for $\alpha$-Pu's ground state. The noncollinear ansatz accurately recovers the experimental phonon density of states, heat capacity, and thermal expansion. These new results on $\alpha$-Pu along with recent results on $\delta$-Pu demonstrate the efficacy of noncollinear ansatzes for the description of plutonium.

Graphene oxide (GO) materials are widely studied, and yet their atomic-scale structures remain to be fully understood. Here we show that the chemical and configurational space of GO can be rapidly explored by advanced machine-learning methods, combining on-the-fly acceleration for first-principles molecular dynamics with message-passing neural-network potentials. The first step allows for the rapid sampling of chemical structures with very little prior knowledge required; the second step affords state-of-the-art accuracy and predictive power. We apply the method to the thermal reduction of GO, which we describe in a realistic (ten-nanometre scale) structural model. Our simulations are consistent with recent experimental findings and help to rationalise them in atomistic and mechanistic detail. More generally, our work provides a platform for routine, accurate, and predictive simulations of diverse carbonaceous materials.

Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $\mathbf{E}^{p}$ and stress $\mathbf{\sigma}$ tensors and volume fraction $c$ of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of $\mathbf{\sigma}$ and $\mathbf{E}^{p}$, there are no strict kinetic equations for them. Here, we develop combined experimental-computational approaches to determine all fields in strongly plastically predeformed Zr and kinetic equation for $\alpha$-$\omega$ PT consistent with experimental data for the entire sample. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions., Comment: 28 pages, 11 figures

Crystalline alloys and related mixed systems make up a large family of materials with high tunability which have been proposed as the solution to a large number of energy related materials design problems. Due to the presence of chemical order and disorder in these systems, neither experimental efforts nor ab-initio computational methods alone are sufficient to span the inherently large configuration space. Therefore, fast and accurate models are necessary. To this end, cluster expansions have been widely and successfully used for the past decades. Cluster expansions are generalized Ising models designed to predict the energy of any atomic configuration of a system after training on a small subset of the available configurations. Constructing and sampling a cluster expansion consists of multiple steps that have to be performed with care. In this tutorial, we provide a comprehensive guide to this process, highlighting important considerations and potential pitfalls. The tutorial consists of three parts, starting with cluster expansion construction for a relatively simple system, continuing with strategies for more challenging systems such as surfaces and closing with examples of Monte Carlo sampling of cluster expansions to study order-disorder transitions and phase diagrams., Comment: 18 pages; 12 figures

Altermagnets constitute a novel, third fundamental class of collinear magnetic ordered materials, alongside with ferro- and antiferromagnets. They share with conventional antiferromagnets the feature of a vanishing net magnetization. At the same time they show a spin-splitting of electronic bands, just as in ferromagnets, caused by the atomic exchange interaction. Interestingly, the altermagnetic spin-splitting may in principle allow for topological electronic structures far beyond collinear antiferromagnets with spin-degenerate bands. Here, using high-resolution and spin angle-resolved photoemission spectroscopy, we observe that the momentum-dependent spin-splitting in altermagnetic CrSb reaches up to 1 eV, while vanishing at high symmetry lines in the Brillouin zone. We find it causes spin-carrying altermagnetic Weyl points, and observe their protected Fermi-arc surface states. This establishes that in altermagnets the large energy scale intrinsic to the spin-splitting - orders of magnitude larger that the relativistic spin-orbit coupling - creates its own realm of robust electronic topology, with spin-polarized topological surface states that are altermagnetic., Comment: 16 pages, 4 figures, comments are welcome

We present a scalable machine learning (ML) force-field model for the adiabatic dynamics of cooperative Jahn-Teller (JT) systems. Large scale dynamical simulations of the JT model also shed light on the orbital ordering dynamics in colossal magnetoresistance manganites. The JT effect in these materials describes the distortion of local oxygen octahedra driven by a coupling to the orbital degrees of freedom of $e_g$ electrons. An effective electron-mediated interaction between the local JT modes leads to a structural transition and the emergence of long-range orbital order at low temperatures. Assuming the principle of locality, a deep-learning neural-network model is developed to accurately and efficiently predict the electron-induced forces that drive the dynamical evolution of JT phonons. A group-theoretical method is utilized to develop a descriptor that incorporates the combined orbital and lattice symmetry into the ML model. Large-scale Langevin dynamics simulations, enabled by the ML force-field models, are performed to investigate the coarsening dynamics of the composite JT distortion and orbital order after a thermal quench. The late-stage coarsening of orbital domains exhibits pronounced freezing behaviors which are likely related to the unusual morphology of the domain structures. Our work highlights a promising avenue for multi-scale dynamical modeling of correlated electron systems., Comment: 17 pages, 11 figures

A simple micromechanical model of polycrystalline materials is proposed, which enables us to swiftly produce grain-boundary-stress distributions induced by the uniform external loading (in the elastic strain regime). Such statistical knowledge of local stresses is a necessary prerequisite to assess the probability for intergranular cracking initiation. Model predictions are verified through finite element calculations for various loading configurations, material properties, and grain-boundary types specified by the properties of a bicrystal pair of grains enclosing the grain boundary.

Machine Learning (ML) potentials such as Gaussian Approximation Potential (GAP) have demonstrated impressive capabilities in mapping structure to properties across diverse systems. Here, we introduce a GAP model for low-dimensional Ni nanoclusters and demonstrate its flexibility and effectiveness in capturing the energetics, structural diversity and thermodynamic properties of Ni nanoclusters across a broad size range. Through a systematic approach encompassing model development, validation, and application, we evaluate the model's efficacy in representing energetics and configurational features in low-dimensional regimes, while also examining its extrapolative nature to vastly different spatiotemporal regimes. Our analysis and discussion shed light on the data quality required to effectively train such models. Trajectories from large scale MD simulations using the GAP model analyzed with data-driven models like Graph Neural Networks (GNN) reveal intriguing insights into the size-dependent phase behavior and thermo-mechanical stability characteristics of porous Ni nanoparticles. Overall, our work underscores the potential of ML models which coupled with data-driven approaches serve as versatile tools for studying low-dimensional systems and complex material dynamics.

The superconducting properties of layered materials can be controlled by thinning, stacking, and twisting, demanding investigation of electronic states by spectroscopic means at the nanometer scale. Here, we reveal the spatial variations of the electronic states in heterostructures of the superconducting monolayer NbSe$_2$/graphene using spectroscopic-imaging scanning tunneling microscopy. The NbSe$_2$ monolayer grown by molecular beam epitaxy is naturally twisted with respect to the graphene substrate and exhibits interference patterns of Bogoliubov quasiparticles twisted with respect to the NbSe$_2$ and graphene lattices. We find that the twisted interference patterns originate from a sextet of regions in momentum space where the Fermi surfaces of NbSe$_2$ and graphene overlap. The Fermi surface overlap is sensitive to the twist angle, providing a knob to tune superconductivity., Comment: 20 pages + 8 pages supplementary, 5 figures in main text + 6 supplementary figures

There is a growing attention given to utilizing Lagrangian and Hamiltonian mechanics with network training in order to incorporate physics into the network. Most commonly, conservative systems are modeled, in which there are no frictional losses, so the system may be run forward and backward in time without requiring regularization. This work addresses systems in which the reverse direction is ill-posed because of the dissipation that occurs in forward evolution. The novelty is the use of Morse-Feshbach Lagrangian, which models dissipative dynamics by doubling the number of dimensions of the system in order to create a mirror latent representation that would counterbalance the dissipation of the observable system, making it a conservative system, albeit embedded in a larger space. We start with their formal approach by redefining a new Dissipative Lagrangian, such that the unknown matrices in the Euler-Lagrange's equations arise as partial derivatives of the Lagrangian with respect to only the observables. We then train a network from simulated training data for dissipative systems such as Fickian diffusion that arise in materials sciences. It is shown by experiments that the systems can be evolved in both forward and reverse directions without regularization beyond that provided by the Morse-Feshbach Lagrangian. Experiments of dissipative systems, such as Fickian diffusion, demonstrate the degree to which dynamics can be reversed.

Energy dispersive X-ray (EDX) spectroscopy in the transmission electron microscope is a key tool for nanomaterials analysis, providing a direct link between spatial and chemical information. However, using it for precisely determining chemical compositions presents challenges of noisy data from low X-ray yields and mixed signals from phases that overlap along the electron beam trajectory. Here, we introduce a novel method, non-negative matrix factorisation based pan-sharpening (PSNMF), to address these limitations. Leveraging the Poisson nature of EDX spectral noise and binning operations, PSNMF retrieves high quality phase spectral and spatial signatures via consecutive factorisations. After validating PSNMF with synthetic datasets of different noise levels, we illustrate its effectiveness on two distinct experimental cases: a nano-mineralogical lamella, and supported catalytic nanoparticles. Not only does PSNMF obtain accurate phase signatures, datasets reconstructed from the outputs have demonstrably lower noise and better fidelity than from the benchmark denoising method of principle component analysis.

Recently, topolectrical circuits (TECs) boom in studying the topological states of matter. The resemblance between circuit Laplacians and tight-binding models in condensed matter physics allows for the exploration of exotic topological phases on the circuit platform. In this review, we begin by presenting the basic equations for the circuit elements and units, along with the fundamentals and experimental methods for TECs. Subsequently, we retrospect the main literature in this field, encompassing the circuit realization of (higher-order) topological insulators and semimetals. Due to the abundant electrical elements and flexible connections, many unconventional topological states like the non-Hermitian, nonlinear, non-Abelian, non-periodic, non-Euclidean, and higher-dimensional topological states that are challenging to observe in conventional condensed matter physics, have been observed in circuits and summarized in this review. Furthermore, we show the capability of electrical circuits for exploring the physical phenomena in other systems, such as photonic and magnetic ones. Importantly, we highlight TEC systems are convenient for manufacture and miniaturization because of their compatibility with the traditional integrated circuits. Finally, we prospect the future directions in this exciting field, and connect the emerging TECs with the development of topology physics, (meta)material designs, and device applications.

Empirical surface treatments, such as low-temperature baking (LTB) in a gaseous atmosphere or in vacuum, are important for the surface preparation of Nb superconducting radio frequency (SRF) cavities. These treatments inhomogeneously dope the first $\sim$50 nm of Nb's subsurface and are expected to impart depth-dependent characteristics to its Meissner response; however, direct evidence supporting this remains elusive, suggesting the effect is subtle. In this work, we revisit the Meissner profile data for several LTB treatments obtained from low-energy muon spin rotation (LE-$\mu$SR) experiments [A. Romanenko et al., Appl. Phys. Lett. 104, 072601 (2014) and R. M. L. McFadden et al., Phys. Rev. Appl. 19, 044018 (2023)], and search for signatures of inhomogeneous field screening. Using a generalized London expression with a recently proposed empirical model for a depth-dependent magnetic penetration depth $\lambda(z)$, we obtain improved fits to the Meissner data, revealing that the presence of a non-superconducting surface "dead layer" $d \geq 25$ nm is a strong indicator of a reduced supercurrent density at shallow subsurface depths. Our analysis supports the notion that vacuum annealing at 120 $^{\circ}$C for 48 h induces a depth-dependent Meissner response, which has consequences for Nb's ability to maintain a magnetic-flux-free state. Evidence of similar behavior from a "nitrogen infusion" treatment is less compelling. Suggestions for further investigation into the matter are provided., Comment: 6 pages, 3 figures, 1 table

In this work, we report on the concurrent emergence of the directional Kondo behavior and incommensurate magnetic ordering in a layered material. We employ temperature- and magnetic field-dependent resistivity measurements, susceptibility measurements, and high resolution wavelength X-ray diffraction spectroscopy to study the electronic properties of AgCrSe$_2$. Impurity Kondo behavior with a characteristic temperature of $T_\text K$ = 32 K is identified through quantitative analysis of the in-plane resistivity, substantiated by magneto-transport measurements. The agreement between our experimental data and the Schlottmann's scaling theory allows us to determine the impurity spin as $S$ = 3/2. Furthermore, we discuss the origin of the Kondo behavior and its relation to the material's antiferromagnetic transition. Our study uncovers an unusual phenomenon -- the equivalence of the N\'eel temperature and the Kondo temperature -- paving the way for further investigations into the intricate interplay between impurity physics and magnetic phenomena in quantum materials, with potential applications in advanced electronic and magnetic devices., Comment: 11 pages, 4 figures

Neutron imaging has gained significant importance as a material characterisation technique and is particularly useful to visualise hydrogenous materials in objects opaque to other radiations. Particular fields of application include investigations of hydrogen in metals as well as metal corrosion, thanks to the fact that neutrons can penetrate metals better than e.g. X-rays and are at the same time highly sensitive to hydrogen. However at interfaces for example those that are prone to corrosion, refraction effects sometimes obscure the attenuation image, which is used to for hydrogen quantification. Refraction, as a differential phase effect, diverts the neutron beam away from the interface in the image which leads to intensity gain and intensity loss regions, which are superimposed to the attenuation image, thus obscuring the interface region and hindering quantitative analyses of e.g. hydrogen content in the vicinity of the interface or in an oxide layer. For corresponding effects in X-ray imaging, a phase filter approach was developed and is generally based on transport-of-intensity considerations. Here, we compare such an approach, that has been adapted to neutrons, with another simulation-based assessment using the ray-tracing software McStas. The latter appears superior and promising for future extensions which enable fitting forward models via simulations in order to separate phase and attenuation effects and thus pave the way for overcoming quantitative limitations at refracting interfaces.

The utilization of ferrimagnets in spintronic applications purportedly offers the potential for high-speed and energy-efficient switching of magnetic states, coupled with their notable stability. This collection of characteristics is commonly observed in ferrimagnetic materials that are in close proximity to states of magnetic or angular momentum compensation. This is owing to the presence of two antiparallel ordered magnetic sublattices that exhibit differing responses to variations in temperature or composition. In the vicinity of the magnetic compensation state, an external field has the capacity to disrupt the collinearity between the sublattices, resulting in the codirectionality of the magnetic projections. The existence of the canted phase has been extensively described theoretically, but its experimental investigation remains limited. Here, a violation of antiferromagnetic ordering is detected through a change in the direction of the effective field induced by the spin-orbital torque, without altering the dominant characteristics of the ferrimagnetic structure. This effect is observed both during external heating of the sample and as a consequence of Joule heating with an increase in the transmitted current. In the examined structure of W/Co70Tb30/Ru, the canted phase region is observed at approximately room temperature and at external fields on the order of 0.1 T. Through macrospin modeling and analytical explanations, a correlation between anisotropy, interlattice exchange interaction, and the presence of the canted phase region is established.

Slow magnetization relaxation processes are an important time-dependent property of many magnetic materials. We show that magnetotransport measurements based on a well-established current-reversal method can be utilized to implement a simple and robust screening scheme for such relaxation processes. We demonstrate our approach considering the anomalous Hall effect in a Pt/Co/AlOx trilayer model system, and then explore relaxation in {\tau} -MnAl films. Compared to magnetotransport experiments based on ac lock-in techniques, we find that the dc current-reversal method is particulary sensitive to relaxation processes with relaxation time scales on the order of seconds, comparable to the current-reversal measurement time scales.