Source author record

Michele Ruggeri

Michele Ruggeri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

The mass-lumped midpoint scheme for computational micromagnetics: Newton linearization and application to magnetic skyrmion dynamics

We discuss a mass-lumped midpoint scheme for the numerical approximation of the Landau-Lifshitz-Gilbert equation, which models the dynamics of the magnetization in ferromagnetic materials. In addition to the classical micromagnetic field contributions, our setting covers the non-standard Dzyaloshinskii-Moriya interaction, which is the essential ingredient for the enucleation and stabilization of magnetic skyrmions. Our analysis also includes the inexact solution of the arising nonlinear systems, for which we discuss both a constraint preserving fixed-point solver from the literature and a novel approach based on the Newton method. We numerically compare the two linearization techniques and show that the Newton solver leads to a considerably lower number of nonlinear iterations. Moreover, in a numerical study on magnetic skyrmions, we demonstrate that, for magnetization dynamics that are very sensitive to energy perturbations, the midpoint scheme, due to its conservation properties, is superior to the dissipative tangent plane schemes from the literature.

preprint2021arXiv

Multi-scale simulation of the adsorption of lithium ion on graphite surface: from Quantum Monte Carlo to Molecular Density Functional Theory

The structure of the double-layer formed at the surface of carbon electrodes is governed by the interactions between the electrode and the electrolyte species. However, carbon is notoriously difficult to simulate accurately, even with well-established methods such as electronic Density Functional Theory and Molecular Dynamics. Here we focus on the important case of a lithium ion in contact with the surface of graphite, and we perform a series of reference Quantum Monte Carlo calculations that allow us to benchmark various electronic Density Functional Theory functionals. We then fit an accurate carbon--lithium pair potential, which is used in molecular Density Functional Theory calculations to determine the free energy of the adsorption of the ion on the surface in the presence of water. The adsorption profile in solution differs markedly from the gas phase results, which emphasize the role of the solvent on the properties of the double-layer.

preprint2019arXiv

The saturation assumption yields optimal convergence of two-level adaptive BEM

We consider the convergence of adaptive BEM for weakly-singular and hypersingular integral equations associated with the Laplacian and the Helmholtz operator in 2D and 3D. The local mesh-refinement is driven by some two-level error estimator. We show that the adaptive algorithm drives the underlying error estimates to zero. Moreover, we prove that the saturation assumption already implies linear convergence of the error with optimal algebraic rates.

preprint2018arXiv

Computational micromagnetics with Commics

We present our open-source Python module Commics for the study of the magnetization dynamics in ferromagnetic materials via micromagnetic simulations. It implements state-of-the-art unconditionally convergent finite element methods for the numerical integration of the Landau-Lifshitz-Gilbert equation. The implementation is based on the multiphysics finite element software Netgen/NGSolve. The simulation scripts are written in Python, which leads to very readable code and direct access to extensive post-processing. Together with documentation and example scripts, the code is freely available on GitLab.

preprint2018arXiv

Iterative solution and preconditioning for the tangent plane scheme in computational micromagnetics

The tangent plane scheme is a time-marching scheme for the numerical solution of the nonlinear parabolic Landau-Lifshitz-Gilbert equation (LLG), which describes the time evolution of ferromagnetic configurations. Exploiting the geometric structure of LLG, the tangent plane scheme requires only the solution of one linear variational form per time-step, which is posed in the discrete tangent space determined by the nodal values of the current magnetization. We develop an effective solution strategy for the arising constrained linear systems, which is based on appropriate Householder reflections. We derive possible preconditioners, which are (essentially) independent of the time-step, and prove that the preconditioned GMRES algorithm leads to linear convergence. Numerical experiments underpin the theoretical findings.

preprint2016arXiv

A self-consistent spin-diffusion model for micromagnetics

We propose a three-dimensional micromagnetic model that dynamically solves the Landau-Lifshitz-Gilbert equation coupled to the full spin-diffusion equation. In contrast to previous methods, we solve for the magnetization dynamics and the electric potential in a self-consistent fashion. This treatment allows for an accurate description of magnetization dependent resistance changes. Moreover, the presented algorithm describes both spin accumulation due to smooth magnetization transitions and due to material interfaces as in multilayer structures. The model and its finite-element implementation are validated by current driven motion of a magnetic vortex structure. In a second experiment, the resistivity of a magnetic multilayer structure in dependence of the tilting angle of the magnetization in the different layers is investigated. Both examples show good agreement with reference simulations and experiments respectively.

preprint2016arXiv

Low-density phases of $^3$He monolayers adsorbed on graphite

Quantum Monte Carlo simulations at zero temperature of a $^3$He monolayer adsorbed on graphite, either clean or preplated with $^4$He, unexpectedly point to a gas-liquid phase transition at a very low areal density of the order of 0.01Å$^{-2}$. This result stems from an essentially unbiased calculation of the ground state energy for an infinite, defect-free substrate which interacts with He atoms via a realistic potential, whereas the interaction between two He atoms includes two- and three-body terms. The sensitivity of the gas-liquid coexistence region on the model Hamiltonian employed is discussed.

preprint2015arXiv

A three-dimensional spin-diffusion model for micromagnetics

We implement a finite-element scheme that solves the Landau-Lifshitz-Gilbert equation coupled to a diffusion equation accounting for spin-polarized currents. The latter solves for the spin accumulation not only in magnetic materials but also in nonmagnetic conductors. The presented method incorporates the model by Slonczewski for the description of spin torque in magnetic multilayers as well as the model of Zhang and Li for the description of current driven domain-wall motion. Furthermore it is able to do both resolve the time evolution of the spin accumulation or treat it in an adiabatic fashion by the choice of sufficiently large time steps.

preprint2015arXiv

Iterative backflow renormalization procedure for many-body ground state wave functions of strongly interacting normal Fermi liquids

We show how a ground state trial wavefunction of a Fermi liquid can be systematically improved introducing a sequence of renormalized coordinates through an iterative backflow transformation. We apply this scheme to calculate the ground state energy of liquid $^3$He in two dimensions at freezing density using variational and fixed-node diffusion Monte Carlo. Comparing with exact transient estimate results for systems with small number of particles, we find that variance extrapolations provide accurate results for the true ground state together with stringent lower bounds. For larger systems these bounds can in turn be used to quantify the systematic bias of fixed-node calculations. These wave functions are size consistent and the scaling of their computational complexity with the number of particles is the same as for standard backflow wave functions.

preprint2014arXiv

Multiscale modeling in micromagnetics: existence of solutions and numerical integration

Various applications ranging from spintronic devices, giant magnetoresistance sensors, and magnetic storage devices, include magnetic parts on very different length scales. Since the consideration of the Landau-Lifshitz-Gilbert equation (LLG) constrains the maximum element size to the exchange length within the media, it is numerically not attractive to simulate macroscopic parts with this approach. On the other hand, the magnetostatic Maxwell equations do not constrain the element size, but cannot describe the short-range exchange interaction accurately. A combination of both methods allows to describe magnetic domains within the micromagnetic regime by use of LLG and also considers the macroscopic parts by a non-linear material law using the Maxwell equations. In our work, we prove that under certain assumptions on the non-linear material law, this multiscale version of LLG admits weak solutions. Our proof is constructive in the sense that we provide a linear-implicit numerical integrator for the multiscale model such that the numerically computable finite element solutions admit weak $H^1$-convergence (at least for a subsequence) towards a weak solution.

preprint2014arXiv

Quantum Monte Carlo study of a two dimensional dipolar Bose gas in a harmonic trap

We present a Quantum Monte Carlo study of the ground state properties of a two dimensional system of Bose particle with dipole moment in a harmonic trap. The direction of the dipoles is assumed to be fixed by an external field. We study how the system behaves when the direction of the dipole moments is changed. Our analysis is made using Path Integral Monte Carlo simulations in the Grand Canonical ensemble, using the Worm Algorithm. We study systems of increasing size, focusing on the spatial distribution of the atoms inside the trap and on the superfluid fraction of the sample. We find that vertical dipoles form a crystal, characterized by a shell structure, while tilting the dipoles has the effect of inducing a striped structure. We find that similar behaviour can be found in the analogue classical system. Vertical dipoles display nonclassical inertia while a strong tilting decreases superfluidity.

preprint2014arXiv

Spin-polarized transport in ferromagnetic multilayers: An unconditionally convergent FEM integrator

We propose and analyze a decoupled time-marching scheme for the coupling of the Landau-Lifshitz-Gilbert equation with a quasilinear diffusion equation for the spin accumulation. This model describes the interplay of magnetization and electron spin accumulation in magnetic and non-magnetic multilayer structures. Despite the strong nonlinearity of the overall PDE system, the proposed integrator requires only the solution of two linear systems per time-step. Unconditional convergence of the integrator towards weak solutions is proved.

preprint2013arXiv

Quasi-2D Liquid He-3

Quantum Monte Carlo simulations at zero temperature of an ensemble of He-3 atoms adsorbed on Mg and Alkali substrates yield strong evidence of a thermodynamically stable liquid He-3 monolayer on all Alkali substrates, with the possible exception of Li. The effective 2D density is close to 0.02$ inverse square Angstroms on Na, making it the lowest density liquid in Nature. Its existence is underlain by zero-point atomic motion perpendicular to the substrate, whose effect is softening the short-range repulsion of the helium interatomic potential. The monolayer films should turn superfluid at a temperature of order 1 mK. No liquid film is predicted to form on Mg, or on stronger substrates such as graphite.