No Arabic abstract
The textit{Spirit} framework is designed for atomic scale spin simulations of magnetic systems of arbitrary geometry and magnetic structure, providing a graphical user interface with powerful visualizations and an easy to use scripting interface. An extended Heisenberg type spin-lattice Hamiltonian including competing exchange interactions between neighbors at arbitrary distance, higher-order exchange, Dzyaloshinskii-Moriya and dipole-dipole interactions is used to describe the energetics of a system of classical spins localised at atom positions. A variety of common simulations methods are implemented including Monte Carlo and various time evolution algorithms based on the Landau-Lifshitz-Gilbert equation of motion, which can be used to determine static ground state and metastable spin configurations, sample equilibrium and finite temperature thermodynamical properties of magnetic materials and nanostructures or calculate dynamical trajectories including spin torques induced by stochastic temperature or electric current. Methods for finding the mechanism and rate of thermally assisted transitions include the geodesic nudged elastic band method, which can be applied when both initial and final states are specified, and the minimum mode following method when only the initial state is given. The lifetime of magnetic states and rate of transitions can be evaluated within the harmonic approximation of transition-state theory. The framework offers performant CPU and GPU parallelizations. All methods are verified and applications to several systems, such as vortices, domain walls, skyrmions and bobbers are described.
We present a method for performing atomistic spin dynamic simulations. A comprehensive summary of all pertinent details for performing the simulations such as equations of motions, models for including temperature, methods of extracting data and numerical schemes for performing the simulations is given. The method can be applied in a first principles mode, where all interatomic exchange is calculated self-consistently, or it can be applied with frozen parameters estimated from experiments or calculated for a fixed spin-configuration. Areas of potential applications to different magnetic questions are also discussed. The method is finally applied to one situation where the macrospin model breaks down; magnetic switching in ultra strong fields.
We present a stochastic modeling framework for atomistic propagation of a Mode I surface crack, with atoms interacting according to the Lennard-Jones interatomic potential at zero temperature. Specifically, we invoke the Cauchy-Born rule and the maximum entropy principle to infer probability distributions for the parameters of the interatomic potential. We then study how uncertainties in the parameters propagate to the quantities of interest relevant to crack propagation, namely, the critical stress intensity factor and the lattice trapping range. For our numerical investigation, we rely on an automated version of the so-called numerical-continuation enhanced flexible boundary (NCFlex) algorithm.
A solid-state analogue of Stimulated Raman Adiabatic Passage can be implemented in a triple well solid-state system to coherently transport an electron across the wells with exponentially suppressed occupation in the central well at any point of time. Termed coherent tunneling adiabatic passage (CTAP), this method provides a robust way to transfer quantum information encoded in the electronic spin across a chain of quantum dots or donors. Using large scale atomistic tight-binding simulations involving over 3.5 million atoms, we verify the existence of a CTAP pathway in a realistic solid-state system: gated triple donors in silicon. Realistic gate profiles from commercial tools were combined with tight-binding methods to simulate gate control of the donor to donor tunnel barriers in the presence of cross-talk. As CTAP is an adiabatic protocol, it can be analyzed by solving the time independent problem at various stages of the pulse - justifying the use of time-independent tight-binding methods to this problem. Our results show that a three donor CTAP transfer, with inter-donor spacing of 15 nm can occur on timescales greater than 23 ps, well within experimentally accessible regimes. The method not only provides a tool to guide future CTAP experiments, but also illuminates the possibility of system engineering to enhance control and transfer times.
Controlling magnetism in low dimensional materials is essential for designing devices that have feature sizes comparable to several critical length scales that exploit functional spin textures, allowing the realization of low-power spintronic and magneto-electric hardware. [1] Unlike conventional covalently-bonded bulk materials, van der Waals (vdW)-bonded layered magnets [2-4] offer exceptional degrees of freedom for engineering spin textures. [5] However, their structural instability has hindered microscopic studies and manipulations. Here, we demonstrate nanoscale structural control in the layered magnet CrSBr creating novel spin textures down to the atomic scale. We show that it is possible to drive a local structural phase transformation using an electron beam that locally exchanges the bondings in different directions, effectively creating regions that have vertical vdW layers embedded within the horizontally vdW bonded exfoliated flakes. We calculate that the newly formed 2D structure is ferromagnetically ordered in-plane with an energy gap in the visible spectrum, and weak antiferromagnetism between the planes. Our study lays the groundwork for designing and studying novel spin textures and related quantum magnetic phases down to single-atom sensitivity, potentially to create on-demand spin Hamiltonians probing fundamental concepts in physics, [6-10] and for realizing high-performance spintronic, magneto-electric and topological devices with nanometer feature sizes. [11,12]
Molecular dynamics simulations are an invaluable tool in numerous scientific fields. However, the ubiquitous classical force fields cannot describe reactive systems, and quantum molecular dynamics are too computationally demanding to treat large systems or long timescales. Reactive force fields based on physics or machine learning can be used to bridge the gap in time and length scales, but these force fields require substantial effort to construct and are highly specific to a given chemical composition and application. A significant limitation of machine learning models is the use of element-specific features, leading to models that scale poorly with the number of elements. This work introduces the Gaussian multipole (GMP) featurization scheme that utilizes physically-relevant multipole expansions of the electron density around atoms to yield feature vectors that interpolate between element types and have a fixed dimension regardless of the number of elements present. We combine GMP with neural networks to directly compare it to the widely used Behler-Parinello symmetry functions for the MD17 dataset, revealing that it exhibits improved accuracy and computational efficiency. Further, we demonstrate that GMP-based models can achieve chemical accuracy for the QM9 dataset, and their accuracy remains reasonable even when extrapolating to new elements. Finally, we test GMP-based models for the Open Catalysis Project (OCP) dataset, revealing comparable performance to graph convolutional deep learning models. The results indicate that this featurization scheme fills a critical gap in the construction of efficient and transferable machine-learned force fields.