Characterization of partially accessible anisotropic spin chains in the presence of anti-symmetric exchange

We address quantum characterization of anisotropic spin chains in the presence of antisymmetric exchange, and investigate whether the Hamiltonian parameters of the chain may be estimated with precision approaching the ultimate limit imposed by quantum mechanics. At variance with previous approaches, we focus on the information that may be extracted by measuring only two neighbouring spins rather than a global observable on the entire chain. We evaluate the Fisher information (FI) of a two-spin magnetization measure, and the corresponding quantum Fisher information (QFI), for all the relevant parameters, i.e. the spin coupling, the anisotropy, and the Dzyaloshinskii Moriya (DM) parameter. Our results show that the reduced system made of two neighbouring spins may be indeed exploited as a probe to characterize global properties of the entire system. In particular, we find that the ratio between the FI and the QFI is close to unit for a large range of the coupling values. The DM coupling is beneficial for coupling estimation, since it leads to the presence of additional bumps and peaks in the FI and QFI, which are not present in a model that neglects exchange interaction and may be exploited to increase the robustness of the overall estimation procedure. Finally, we address the multiparameter estimation problem, and show that the model is compatible but sloppy, i.e. both the Uhlmann curvature and the determinant of the QFI matrix vanish. Physically, this means that the state of the system actually depends only on a reduced numbers of combinations of parameters, and not on all of them separately.Comment: 10 pages, 7 figure

Characterization of partially accessible anisotropic spin chains in the presence of anti-symmetric exchange

We address quantum characterization of anisotropic spin chains in the presence of anti-symmetric exchange, and investigate whether the Hamiltonian parameters of the chain may be estimated with precision approaching the ultimate limit imposed by quantum mechanics. At variance with previous approaches, we focus on the information that may be extracted by measuring only two neighboring spins rather than a global observable on the entire chain. We evaluate the Fisher information (FI) of a two-spin magnetization measure, and the corresponding quantum Fisher information (QFI), for all the relevant parameters, i.e. the spin coupling, the anisotropy, and the Dzyaloshinskii–Moriya (DM) parameter. Our results show that the reduced system made of two neighboring spins may be indeed exploited as a probe to characterize global properties of the entire system. In particular, we find that the ratio between the FI and the QFI is close to unit for a large range of the coupling values. The DM coupling is beneficial for coupling estimation, since it leads to the presence of additional bumps and peaks in the FI and QFI, which are not present in a model that neglects exchange interaction and may be exploited to increase the robustness of the overall estimation procedure. Finally, we address the multiparameter estimation problem, and show that the model is compatible but sloppy, i.e. both the Uhlmann curvature and the determinant of the QFI matrix vanish. Physically, this means that the state of the system actually depends only on a reduced numbers of combinations of parameters, and not on all of them separately

An automatic procedure to forecast tephra fallout

Tephra fallout constitutes a serious threat to communities around active volcanoes. Reliable short-term forecasts represent a valuable aid for scientists and civil authorities to mitigate the effects of fallout on the surrounding areas during an episode of crisis. We present a platform-independent automatic procedure with the aim to daily forecast transport and deposition of volcanic particles. The procedure builds on a series of programs and interfaces that automate the data flow and the execution and subsequent postprocess of fallout models. Firstly, the procedure downloads regional meteorological forecasts for the area and time interval of interest, filters and converts data from its native format, and runs the CALMET diagnostic model to obtain the wind field and other micro-meteorological variables on a finer local-scale 3-D grid defined by the user. Secondly, it assesses the distribution of mass along the eruptive column, commonly by means of the radial averaged buoyant plume equations depending on the prognostic wind field and on the conditions at the vent (granulometry, mass flow rate, etc). All these data serve as input for the fallout models. The initial version of the procedure includes only two Eulerian models, HAZMAP and FALL3D, the latter available as serial and parallel implementations. However, the procedure is designed to incorporate easily other models in a near future with minor modifications on the model source code. The last step is to postprocess the outcomes of models to obtain maps written in standard file formats. These maps contain plots of relevant quantities such as predicted ground load, expected deposit thickness and, for the case of or 3-D models, concentration on air or flight safety concentration thresholds

Il progetto EPLORIS: La ricostruzione virtuale dell'eruzione del Vesuvio

The main objective of the Exploris project consists in the quantitative analysis of explosive eruption risk in densely populated EU volcanic regions and the evaluation of the likely effectiveness of possible mitigation measures through the development of volcanic risk facilities (such as supercomputer models, vulnerability databases, and probabilistic risk assessment protocols) and their application to high-risk European volcanoes. Exploris’ main ambition is to make a significant step forward in the assessment of explosive eruption risk in highly populated EU cities and islands. For this project, a new simulation model, based on fundamental transport laws to describe the 4D (3D spatial co-ordinates plus time) multiphase flow dynamics of explosive eruptions has been developed and parallelized in INGV and CINECA. Moreover, CINECA developed specific tools to efficiently visualise the results of simulations. This article presents the results of the large numerical simulations, carred out with CINECA’s Supercomputers, to describe the collapse of the volcanic eruption column and the propagation of pyroclastic density currents, for selected medium scale (sub-Plinian) eruptive scenarios at Vesuvius

An application of parallel computing to the simulation of volcanic eruptions

A parallel code for the simulation of the transient 3D dispersal of volcanic particles produced by explosive eruptions is presented. The model transport equations, based on the multiphase flow theory, describe the atmospheric dynamics of the gas-particle mixture ejected through the volcanic crater. The numerics is based on a finite-volume discretization scheme and a pressure-based iterative non-linear solver suited to compressible multiphase flows. The code has been parallelized by adopting an ad hoc domain partitioning scheme that enforces the load balancing. An optimized communication layer has been built over the Message-Passing Interface. The code proved to be remarkably efficient on several high-performance platforms and makes it possible to simulate fully 3D eruptive scenarios on realistic volcano topography

Coin dimensionality as a resource in quantum metrology involving discrete-time quantum walks

We address metrological problems where the parameter of interest is encoded in the internal degree of freedom of a discrete-time quantum walker, and provide evidence that coin dimensionality is a potential resource to enhance precision. In particular, we consider estimation problems where the coin parameter governs rotations around a given axis and show that the corresponding quantum Fisher information (QFI) may increase with the dimension of the coin. We determine the optimal initial state of the walker to maximize the QFI and discuss whether, and to which extent, precision enhancement may be achieved by measuring only the position of the walker. Finally, we consider Grover-like encoding of the parameter and compare results with those obtained from rotation encoding.Comment: revised version, 14 pages, 5 figure

Continuous-time quantum walks in the presence of a quadratic perturbation

We address the properties of continuous-time quantum walks with Hamiltonians of the form H = L + \u3bbL2, with L the Laplacian matrix of the underlying graph and the perturbation \u3bbL2 motivated by its potential use to introduce next-nearest-neighbor hopping. We consider cycle, complete, and star graphs as paradigmatic models with low and high connectivity and/or symmetry. First, we investigate the dynamics of an initially localized walker. Then we devote attention to estimating the perturbation parameter \u3bb using only a snapshot of the walker dynamics. Our analysis shows that a walker on a cycle graph spreads ballistically independently of the perturbation, whereas on complete and star graphs one observes perturbation-dependent revivals and strong localization phenomena. Concerning the estimation of the perturbation, we determine the walker preparations and the simple graphs that maximize the quantum Fisher information. We also assess the performance of position measurement, which turns out to be optimal, or nearly optimal, in several situations of interest. Besides fundamental interest, our study may find applications in designing enhanced algorithms on graphs

The role of bone metastases on the mechanical competence of human vertebrae

Spine is the most common site for bone metastases. The evaluation of the mechanical competence and failure location in metastatic vertebrae is a biomechanical and clinical challenge. Little is known about the failure behaviour of vertebrae with metastatic lesions. The aim of this study was to use combined micro-Computed Tomography (microCT) and time-lapsed mechanical testing to reveal the failure location in metastatic vertebrae. Fifteen spine segments, each including a metastatic and a radiologically healthy vertebra, were tested in compression up to failure within a microCT. Volumetric strains were measured using Digital Volume Correlation. The images of undeformed and deformed specimens were overlapped to identify the failure location. Vertebrae with lytic metastases experienced the largest average compressive strains (median ± standard deviation: −8506 ± 4748microstrain), followed by the vertebrae with mixed metastases (−7035 ± 15605microstrain), the radiologically healthy vertebrae (−5743 ± 5697microstrain), and the vertebrae with blastic metastases (−3150 ± 4641microstrain). Strain peaks were localised within and nearby the lytic lesions or around the blastic tissue. Failure between the endplate and the metastasis was identified in vertebrae with lytic metastases, whereas failure localised around the metastasis in vertebrae with blastic lesions. This study showed for the first time the role of metastases on the vertebral internal deformations. While lytic lesions lead to failure of the metastatic vertebra, vertebrae with blastic metastases are more likely to induce failure in the adjacent vertebrae. Nevertheless, every metastatic lesion affects the vertebral deformation differently, making it essential to assess how the lesion affects the bone microstructure. These results suggest that the properties of the lesion (type, size, location within the vertebral body) should be considered when developing clinical tools to predict the risk of fracture in patients with metastatic lesions

Full configuration interaction approach to the few-electron problem in artificial atoms

We present a new high-performance configuration interaction code optimally designed for the calculation of the lowest energy eigenstates of a few electrons in semiconductor quantum dots (also called artificial atoms) in the strong interaction regime. The implementation relies on a single-particle representation, but it is independent of the choice of the single-particle basis and, therefore, of the details of the device and configuration of external fields. Assuming no truncation of the Fock space of Slater determinants generated from the chosen single-particle basis, the code may tackle regimes where Coulomb interaction very effectively mixes many determinants. Typical strongly correlated systems lead to very large diagonalization problems; in our implementation, the secular equation is reduced to its minimal rank by exploiting the symmetry of the effective-mass interacting Hamiltonian, including square total spin. The resulting Hamiltonian is diagonalized via parallel implementation of the Lanczos algorithm. The code gives access to both wave functions and energies of first excited states. Excellent code scalability in a parallel environment is demonstrated; accuracy is tested for the case of up to eight electrons confined in a two-dimensional harmonic trap as the density is progressively diluted and correlation becomes dominant. Comparison with previous Quantum Monte Carlo simulations in the Wigner regime demonstrates power and flexibility of the method.Comment: RevTeX 4.0, 18 pages, 6 tables, 9 postscript b/w figures. Final version with new material. Section 6 on the excitation spectrum has been added. Some material has been moved to two appendices, which appear in the EPAPS web depository in the published versio

Bone metastases do not affect the measurement uncertainties of a global digital volume correlation algorithm

Introduction: Measurement uncertainties of Digital Volume Correlation (DVC) are influenced by several factors, like input images quality, correlation algorithm, bone type, etc. However, it is still unknown if highly heterogeneous trabecular microstructures, typical of lytic and blastic metastases, affect the precision of DVC measurements. Methods: Fifteen metastatic and nine healthy vertebral bodies were scanned twice in zero-strain conditions with a micro-computed tomography (isotropic voxel size = 39 μm). The bone microstructural parameters (Bone Volume Fraction, Structure Thickness, Structure Separation, Structure Number) were calculated. Displacements and strains were evaluated through a global DVC approach (BoneDVC). The relationship between the standard deviation of the error (SDER) and the microstructural parameters was investigated in the entire vertebrae. To evaluate to what extent the measurement uncertainty is influenced by the microstructure, similar relationships were assessed within sub-regions of interest. Results: Higher variability in the SDER was found for metastatic vertebrae compared to the healthy ones (range 91-1030 με versus 222–599 με). A weak correlation was found between the SDER and the Structure Separation in metastatic vertebrae and in the sub-regions of interest, highlighting that the heterogenous trabecular microstructure only weakly affects the measurement uncertainties of BoneDVC. No correlation was found for the other microstructural parameters. The spatial distribution of the strain measurement uncertainties seemed to be associated with regions with reduced greyscale gradient variation in the microCT images. Discussion: Measurement uncertainties cannot be taken for granted but need to be assessed in each single application of the DVC to consider the minimum unavoidable measurement uncertainty when interpreting the results