Source author record

Stefano Baroni

Stefano Baroni 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

12works
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

12 published item(s)

preprint2022arXiv

SporTran: a code to estimate transport coefficients from the cepstral analysis of (multivariate) current time series

SporTran is a Python utility designed to estimate generic transport coefficients in extended systems, based on the Green-Kubo theory of linear response and the recently introduced cepstral analysis of the current time series generated by molecular dynamics simulations. SporTran can be applied to univariate as well as multivariate time series. Cepstral analysis requires minimum discretion from the user, in that it weakly depends on two parameters, one of which is automatically estimated by a statistical model-selection criterion that univocally determine the resulting accuracy. In order to facilitate the optimal refinement of these parameters, SporTran features an easy-to-use graphical user interface. A command-line interface and a Python API, easy to embed in complex data-analysis workflows, are also provided.

preprint2022arXiv

turboMagnon -- A code for the simulation of spin-wave spectra using the Liouville-Lanczos approach to time-dependent density-functional perturbation theory

We introduce turboMagnon, an implementation of the Liouville-Lanczos approach to linearized time-dependent density-functional theory, designed to simulate spin-wave spectra in solid-state materials. The code is based on the noncollinear spin-polarized framework and the self-consistent inclusion of spin-orbit coupling that allow to model complex magnetic excitations. The spin susceptibility matrix is computed using the Lanczos recursion algorithm that is implemented in two flavors - the non-Hermitian and the pseudo-Hermitian one. turboMagnon is open-source software distributed under the terms of the GPL as a component of Quantum ESPRESSO. As with other components, turboMagnon is optimized to run on massively parallel architectures using native mathematical libraries (LAPACK and FFTW) and a hierarchy of custom parallelization layers built on top of MPI. The effectiveness of the code is showcased by computing magnon dispersions for the CrI$_3$ monolayer, and the importance of the spin-orbit coupling is discussed.

preprint2022arXiv

Viscosity in water from first-principles and deep-neural-network simulations

We report on an extensive study of the viscosity of liquid water at near-ambient conditions, performed within the Green-Kubo theory of linear response and equilibrium ab initio molecular dynamics (AIMD), based on density-functional theory (DFT). In order to cope with the long simulation times necessary to achieve an acceptable statistical accuracy, our ab initio approach is enhanced with deep-neural-network potentials (NNP). This approach is first validated against AIMD results, obtained by using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional and paying careful attention to crucial, yet often overlooked, aspects of the statistical data analysis. Then, we train a second NNP to a dataset generated from the Strongly Constrained and Appropriately Normed (SCAN) functional. Once the error resulting from the imperfect prediction of the melting line is offset by referring the simulated temperature to the theoretical melting one, our SCAN predictions of the shear viscosity of water are in very good agreement with experiments.

preprint2020arXiv

Gauge optimization of time series for thermal-transport simulations

Thermal and other transport coefficients were recently shown to be largely independent of the microscopic representation of the energy (current) densities or, more generally, of the relevant conserved densities/currents. In this paper we show how this gauge invariance, which is intimately related to the intrinsic indeterminacy of the energy of isolated atoms, can be exploited to optimize the statistical properties of the current time series from which the transport coefficients can be evaluated. To this end, we make use of a variational principle that relies on the metric properties of the conserved currents, treated as elements of an abstract linear space. Different metrics would result in different variational principles. We finally show how a recently proposed data-analysis methodology based on the theory of transport in multi-component systems can be recovered by a suitable choice of this metric.

preprint2020arXiv

Heat and charge transport in H$_2$O at ice-giant conditions from ab initio molecular dynamics simulations

The impact of the inner structure and thermal history of planets on their observable features, such as luminosity or magnetic field, crucially depends on the poorly known heat and charge transport properties of their internal layers. The thermal and electric conductivities of different phases of water (liquid, solid, and super-ionic) occurring in the interior of ice giant planets, such as Uranus or Neptune, are evaluated from equilibrium ab initio molecular dynamics, leveraging recent progresses in the theory and data analysis of transport in extended systems. The implications of our findings on the evolution models of the ice giants are briefly discussed

preprint2020arXiv

Thermal Evolution of Uranus with a Frozen Interior

The intrinsic luminosity of Uranus is a factor of 10 less than that of Neptune, an observation that standard giant planetary evolution models, which assume negligible viscosity, fail to capture. Here we show that more than half of the interior of Uranus is likely to be in a solid state, and that thermal evolution models that account for this high viscosity region satisfy the observed faintness of Uranus by storing accretional heat deep in the interior. A frozen interior also explains the quality factor of Uranus required by the evolution of the orbits of its satellites.

preprint2015arXiv

Microscopic theory and quantum simulation of atomic heat transport

Quantum simulation methods based on density-functional theory are currently deemed unfit to cope with atomic heat transport within the Green-Kubo formalism, because quantum-mechanical energy densities and currents are inherently ill-defined at the atomic scale. We show that, while this difficulty would also affect classical simulations, thermal conductivity is indeed insensitive to such ill-definedness by virtue of a sort of gauge invariance resulting from energy extensivity and conservation. Based on these findings, we derive an expression for the adiabatic energy flux from density-functional theory, which allows heat transport to be simulated using ab-initio equilibrium molecular dynamics. Our methodology is demonstrated by comparing its predictions with those of classical equilibrium and ab-initio non-equilibrium (Müller-Plathe) simulations for a liquid-Argon model, and finally applied to heavy water at ambient conditions.

preprint2014arXiv

turboTDDFT 2.0 - Hybrid functionals and new algorithms within time-dependent density-functional perturbation theory

We present a new release of the turboTDDFT code featuring an implementation of hybrid functionals, a recently introduced pseudo-Hermitian variant of the Liouville-Lanczos approach to time-dependent density-functional perturbation theory, and a newly developed Davidson-like algorithm to compute selected interior eigenvalues/vectors of the Liouvillian super-operator. Our implementation is thoroughly validated against benchmark calculations performed on the cyanin (C$_{21}$O$_{11}$H$_{21}$) molecule using the Gaussian09 and turboTDDFT 1.0 codes.

preprint2013arXiv

Electron energy-loss and inelastic X-ray scattering cross sections from time-dependent density-functional perturbation theory

The Liouville-Lanczos approach to linear-response time-dependent density-functional theory is generalized so as to encompass electron energy-loss and inelastic X-ray scattering spectroscopies in periodic solids. The computation of virtual orbitals and the manipulation of large matrices are avoided by adopting a representation of response orbitals borrowed from (time-independent) density-functional perturbation theory and a suitable Lanczos recursion scheme. The latter allows the bulk of the numerical work to be performed at any given transferred momentum only once, for a whole extended frequency range. The numerical complexity of the method is thus greatly reduced, making the computation of the loss function over a wide frequency range at any given transferred momentum only slightly more expensive than a single standard ground-state calculation, and opening the way to computations for systems of unprecedented size and complexity. Our method is validated on the paradigmatic examples of bulk silicon and aluminum, for which both experimental and theoretical results already exist in the literature.

preprint2011arXiv

Thermal properties of materials from ab-initio quasi-harmonic phonons

This paper gives a short overview of the calculation of thermal properties of materials from first principles, using the Quasi-Harmonic Approximation (QHA). We first introduce some of the thermal properties of interest and describe how they can be calculated in the framework of the QHA; then we briefly recall Density-Functional Perturbation Theory as a tool for calculating phonons from first principles, and present some codes that implement it; finally we review recent applications of first-principle QHA.

preprint2010arXiv

Reptation quantum Monte Carlo for lattice Hamiltonians with a directed-update scheme

We provide an extension to lattice systems of the reptation quantum Monte Carlo algorithm, originally devised for continuous Hamiltonians. For systems affected by the sign problem, a method to systematically improve upon the so-called fixed-node approximation is also proposed. The generality of the method, which also takes advantage of a canonical worm algorithm scheme to measure off-diagonal observables, makes it applicable to a vast variety of quantum systems and eases the study of their ground-state and excited-states properties. As a case study, we investigate the quantum dynamics of the one-dimensional Heisenberg model and we provide accurate estimates of the ground-state energy of the two-dimensional fermionic Hubbard model.

preprint2010arXiv

The itinerant ferromagnetic phase of the Hubbard model

Using a newly developed quantum Monte Carlo technique, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, a discontinuous transition to a paramagnetic phase is observed, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition. The remarkable stability of itinerant ferromagnetism renews the hope to describe this phenomenon within a purely kinetic mechanism and will facilitate the validation of experimental quantum simulators with cold atoms loaded in optical lattices.