Researcher profile

Fabiana Zama

Fabiana Zama contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2022arXiv

An automatic L1-based regularization method for the analysis of FFC dispersion profiles with quadrupolar peaks

Fast Field-Cycling Nuclear Magnetic Resonance relaxometry is a non-destructive technique to investigate molecular dynamics and structure of systems having a wide range of applications such as environment, biology, and food. Besides a considerable amount of literature about modeling and application of such technique in specific areas, an algorithmic approach to the related parameter identification problem is still lacking. We believe that a robust algorithmic approach will allow a unified treatment of different samples in several application areas. In this paper, we model the parameters identification problem as a constrained $L_1$-regularized non-linear least squares problem. Following the approach proposed in [Analytical Chemistry 2021 93 (24)], the non-linear least squares term imposes data consistency by decomposing the acquired relaxation profiles into relaxation contributions associated with 1H-1H and 1H-14N dipole-dipole interactions. The data fitting and the L1-based regularization terms are balanced by the so-called regularization parameter. For the parameters identification, we propose an algorithm that computes, at each iteration, both the regularization parameter and the model parameters. In particular, the regularization parameter value is updated according to a Balancing Principle and the model parameters values are obtained by solving the corresponding $L_1$-regularized non-linear least squares problem by means of the non-linear Gauss-Seidel method. We analyse the convergence properties of the proposed algorithm and run extensive testing on synthetic and real data. A Matlab software, implementing the presented algorithm, is available upon request to the authors.

preprint2022arXiv

Explicit Exactly Energy-conserving Methods for Hamiltonian Systems

For Hamiltonian systems, simulation algorithms that exactly conserve numerical energy or pseudo-energy have seen extensive investigation. Most available methods either require the iterative solution of nonlinear algebraic equations at each time step, or are explicit, but where the exact conservation property depends on the exact evaluation of an integral in continuous time. Under further restrictions, namely that the potential energy contribution to the Hamiltonian is non-negative, newer techniques based on invariant energy quadratisation allow for exact numerical energy conservation and yield linearly implicit updates, requiring only the solution of a linear system at each time step. In this article, it is shown that, for a general class of Hamiltonian systems, and under the non-negativity condition on potential energy, it is possible to arrive at a fully explicit method that exactly conserves numerical energy. Furthermore, such methods are unconditionally stable, and are of comparable computational cost to the very simplest integration methods (such as Stormer-Verlet). A variant of this scheme leading to a conditionally-stable method is also presented, and follows from a splitting of the potential energy. Various numerical results are presented, in the case of the classic test problem of Fermi, Pasta and Ulam, as well as for nonlinear systems of partial differential equations, including those describing high amplitude vibration of strings and plates.

preprint2022arXiv

MUPen2DTool: a Matlab Tool for 2D Nuclear Magnetic Resonance relaxation data inversion

Accurate and efficient analysis of materials properties from Nuclear Magnetic Resonance (NMR) relaxation data requires robust and efficient inversion procedures. Despite the great variety of applications requiring to process two-dimensional NMR data (2DNMR), a few software tools are freely available. The aim of this paper is to present MUPen2DTool, an open-source MATLAB based software tool for 2DNMR data inversion. The user can choose among several types of NMR experiments, and the software provides codes that can be used and extended easily. Furthermore, a MATLAB interface makes it easier to include users own data. The practical use is demonstrated in the reported examples of both synthetic and real NMR data.

preprint2020arXiv

2DNMR data inversion using locally adapted multi-penalty regularization

A crucial issue in two-dimensional Nuclear Magnetic Resonance (NMR) is the speed and accuracy of the data inversion. This paper proposes a multi-penalty method with locally adapted regularization parameters for fast and accurate inversion of 2DNMR data. The method solves an unconstrained optimization problem whose objective contains a data-fitting term, a single $L1$ penalty parameter and a multiple parameter $L2$ penalty. We propose an adaptation of the Fast Iterative Shrinkage and Thresholding (FISTA) method to solve the multi-penalty minimization problem, and an automatic procedure to compute all the penalty parameters. This procedure generalizes the Uniform Penalty principle introduced in [Bortolotti et al., \emph{Inverse Problems}, 33(1), 2016]. The proposed approach allows us to obtain accurate relaxation time distributions while keeping short the computation time. Results of numerical experiments on synthetic and real data prove that the proposed method is efficient and effective in reconstructing the peaks and the flat regions that usually characterize NMR relaxation time distributions.

preprint2020arXiv

Preliminary analysis of COVID-19 spread in Italy with an adaptive SEIRD model

In this paper we propose a Susceptible-Infected-Exposed-Recovered-Dead (SEIRD) differential model for the analysis and forecast of the COVID-19 spread in some regions of Italy, using the data from the Italian Protezione Civile from February 24th 2020. In this study investigate an adaptation of the model. Since several restricting measures have been imposed by the Italian government at different times, starting from March 8th 2020, we propose a modification of SEIRD by introducing a time dependent transmitting rate. In the numerical results we report the maximum infection spread for the three Italian regions firstly affected by the COVID-19 outbreak(Lombardia, Veneto and Emilia Romagna). This approach will be successively extended to other Italian regions, as soon as more data will be available.