Calculating the dynamic magnetic behaviour of molecules is key for advances in spintronics and quantum information science; however, accurately modelling this process has remained challenging. A new Python package named qdmag computes the magnetization within magnetic molecules by solving a generalised Lindblad quantum master equation that incorporates spin-phonon coupling as a means of dissipation. The software accurately models how magnetism changes within molecules over extended periods; existing computational tools previously struggled with this process.

The package calculates magnetization by simulating interactions between electron spins and vibrations, known as spin-phonon coupling, using complex mathematical descriptions of magnetic behaviour. This enables investigation into dynamic processes inaccessible using older methods, benefiting both physicists and chemists studying magnetic materials. Northeastern University unveiled qdmag, new software designed to accurately simulate how magnetism evolves within molecules over time; current computational methods often struggle with this complex process.

The package functions by modelling interactions between electron spins and vibrations, known as spin-phonon coupling, using sophisticated mathematical descriptions akin to writing out all the ingredients and instructions for baking a complex cake, termed the spin Hamiltonian formalism. Importantly, qdmag accounts for tiny variations in energy levels even without an external magnetic field, similar to subtle differences in pitch when strumming different strings on a guitar, using what is called zero-field splitting. This allows scientists to investigate dynamic processes previously beyond reach. Further details regarding the underlying equations and implementation are presented below.

Millisecond magnetisation dynamics now simulated via a staircase approximation technique

A Python package called qdmag can simulate magnetization evolution up to milliseconds; existing software was limited to nanosecond timescales because of computational constraints. This breakthrough stems from a new ‘staircase’ technique approximating continuously changing magnetic fields as discrete stepped levels, dramatically reducing processing demands without sacrificing accuracy. Employing this method alongside calculations with spin Hamiltonians enables complex descriptions of molecular magnetism within qdmag, overcoming previous limitations hindering long-time dynamic studies of magnetic molecules and materials.

Lower-dimensional effective Hamiltonians for complex systems were constructed, enabling analysis even when direct calculation was impractical, further expanding the scope of accessible research. The package modelled systems containing up to three interacting spins, previously challenging due to computational limits, incorporating ‘spin-phonon coupling’, representing how vibrations influence magnetism; adjustable parameters allow exploration of sub-Ohmic, Ohmic, and super-Ohmic behaviours describing vibrational energy transfer. In particular, confirming consistent results regardless of specific values when their product remained fixed; this highlights strong performance in the method.

Molecular magnetism underpins advances in fields ranging from high-density data storage to quantum technologies; accurately predicting material behaviour demands increasingly sophisticated computational methods. While qdmag successfully simulates magnetization dynamics up to milliseconds, it currently does not account for external perturbations beyond varying magnetic fields or more intricate environmental effects which could impact real-world material behaviour. Northeastern University scientists created qdmag, a new computational tool simulating molecular magnetic behaviour with unprecedented duration and accuracy, tackling longstanding challenges in modelling complex spin dynamics within materials.

Spin-phonon coupling calculations, describing electron spins and vibrations, are now possible using qdmag allowing exploration of previously inaccessible phenomena. Representing fluctuating magnetic fields via an efficient ‘staircase’ method, approximating continuous change with discrete steps, enables calculations for systems lasting up to milliseconds, bridging theoretical models with practical numerical analysis. The team’s focus on establishing a strong computational framework handling spin-phonon coupling through a generalised Lindblad equation alongside techniques like staircase approximation offers valuable groundwork applicable across diverse materials.

Northeastern University scientists developed qdmag, a new utility that calculates the magnetization of magnetic molecules incorporating interactions between electron spins and vibrations, known as spin-phonon coupling. This tool allows researchers to simulate how magnetism evolves in complex molecular systems over timescales reaching several milliseconds, exceeding previous capabilities due to its efficient ‘staircase’ method for modelling time-varying fields.

By supporting calculations for up to three interacting spins with adjustable parameters, qdmag provides a means to connect theoretical predictions with numerical analysis of dynamic behaviour. The authors demonstrated its use through case studies involving one-, two-, and three-spin examples.

👉 More information
🗞 qdmag: A Python package for simulating nonequilibrium magnetization of molecules using quantum master equations
✍️ Shuanglong Liu, Xiao Chen, Andrew Cupo, James N. Fry and Hai-Ping Cheng
🧠 ArXiv: https://arxiv.org/abs/2609.16352

More like this