All-electron magnetic response with pseudopotentials: NMR chemical shifts

A theory for the ab initio calculation of all-electron NMR chemical shifts in insulators using pseudopotentials is presented. It is formulated for both finite and infinitely periodic systems and is based on an extension to the Projector Augmented Wave approach of Bloechl [P. E. Bloechl, Phys. Rev. B 50, 17953 (1994)] and the method of Mauri et al [F. Mauri, B.G. Pfrommer, and S.G. Louie, Phys. Rev. Lett. 77, 5300 (1996)]. The theory is successfully validated for molecules by comparison with a selection of quantum chemical results, and in periodic systems by comparison with plane-wave all-electron results for diamond.Comment: 25 pages, 4 tables, submitted to Physical Review

Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations

We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyse the factors that limit the size of the integration time step in approaches based on plane-wave expansions. The maximum allowed time step is dictated by the highest frequency components of the fictitious electronic dynamics. These can result either from the large wavevector components of the kinetic energy or from the small wavevector components of the Coulomb potential giving rise to the so called {\it charge sloshing} problem. We show how to eliminate large wavevector instabilities by adopting a preconditioning scheme that is implemented here for the first-time in the context of Car-Parrinello ab-initio molecular dynamics simulations of the ionic motion. We also show how to solve the charge-sloshing problem when this is present. We substantiate our theoretical analysis with numerical tests on a number of different silicon and carbon systems having both insulating and metallic character.Comment: RevTex, 9 figures available upon request, to appear in Phys. Rev.

A Self-Consistent First-Principles Technique Having Linear Scaling

An algorithm for first-principles electronic structure calculations having a computational cost which scales linearly with the system size is presented. Our method exploits the real-space localization of the density matrix, and in this respect it is related to the technique of Li, Nunes and Vanderbilt. The density matrix is expressed in terms of localized support functions, and a matrix of variational parameters, L, having a finite spatial range. The total energy is minimized with respect to both the support functions and the elements of the L matrix. The method is variational, and becomes exact as the ranges of the support functions and the L matrix are increased. We have tested the method on crystalline silicon systems containing up to 216 atoms, and we discuss some of these results.Comment: 12 pages, REVTeX, 2 figure

Potential for Lung Recruitment and Ventilation-Perfusion Mismatch in Patients With the Acute Respiratory Distress Syndrome From Coronavirus Disease 2019

OBJECTIVES: Severe cases of coronavirus disease 2019 develop the acute respiratory distress syndrome, requiring admission to the ICU. This study aimed to describe specific pathophysiological characteristics of acute respiratory distress syndrome from coronavirus disease 2019. DESIGN: Prospective crossover physiologic study. SETTING: ICU of a university-affiliated hospital from northern Italy dedicated to care of patients with confirmed diagnosis of coronavirus disease 2019. PATIENTS: Ten intubated patients with acute respiratory distress syndrome and confirmed diagnosis of coronavirus disease 2019. INTERVENTIONS: We performed a two-step positive end-expiratory pressure trial with change of 10\u2009cm H2O in random order. MEASUREMENTS AND MAIN RESULTS: At each positive end-expiratory pressure level, we assessed arterial blood gases, respiratory mechanics, ventilation inhomogeneity, and potential for lung recruitment by electrical impedance tomography. Potential for lung recruitment was assessed by the recently described recruitment to inflation ratio. In a subgroup of seven paralyzed patients, we also measured ventilation-perfusion mismatch at lower positive end-expiratory pressure by electrical impedance tomography. At higher positive end-expiratory pressure, respiratory mechanics did not change significantly: compliance remained relatively high with low driving pressure. Oxygenation and ventilation inhomogeneity improved but arterial CO2 increased despite unchanged respiratory rate and tidal volume. The recruitment to inflation ratio presented median value higher than previously reported in acute respiratory distress syndrome patients but with large variability (median, 0.79 [0.53-1.08]; range, 0.16-1.40). The FIO2 needed to obtain viable oxygenation at lower positive end-expiratory pressure was significantly correlated with the recruitment to inflation ratio (r = 0.603; p = 0.05). The ventilation-perfusion mismatch was elevated (median, 34% [32-45%] of lung units) and, in six out of seven patients, ventilated nonperfused units represented a much larger proportion than perfused nonventilated ones. CONCLUSIONS: In patients with acute respiratory distress syndrome from coronavirus disease 2019, potential for lung recruitment presents large variability, while elevated dead space fraction may be a specific pathophysiological trait. These findings may guide selection of personalized mechanical ventilation settings

Efficient Recursion Method for Inverting Overlap Matrix

A new O(N) algorithm based on a recursion method, in which the computational effort is proportional to the number of atoms N, is presented for calculating the inverse of an overlap matrix which is needed in electronic structure calculations with the the non-orthogonal localized basis set. This efficient inverting method can be incorporated in several O(N) methods for diagonalization of a generalized secular equation. By studying convergence properties of the 1-norm of an error matrix for diamond and fcc Al, this method is compared to three other O(N) methods (the divide method, Taylor expansion method, and Hotelling's method) with regard to computational accuracy and efficiency within the density functional theory. The test calculations show that the new method is about one-hundred times faster than the divide method in computational time to achieve the same convergence for both diamond and fcc Al, while the Taylor expansion method and Hotelling's method suffer from numerical instabilities in most cases.Comment: 17 pages and 4 figure

Basis Functions for Linear-Scaling First-Principles Calculations

In the framework of a recently reported linear-scaling method for density-functional-pseudopotential calculations, we investigate the use of localized basis functions for such work. We propose a basis set in which each local orbital is represented in terms of an array of `blip functions'' on the points of a grid. We analyze the relation between blip-function basis sets and the plane-wave basis used in standard pseudopotential methods, derive criteria for the approximate equivalence of the two, and describe practical tests of these criteria. Techniques are presented for using blip-function basis sets in linear-scaling calculations, and numerical tests of these techniques are reported for Si crystal using both local and non-local pseudopotentials. We find rapid convergence of the total energy to the values given by standard plane-wave calculations as the radius of the linear-scaling localized orbitals is increased.Comment: revtex file, with two encapsulated postscript figures, uses epsf.sty, submitted to Phys. Rev.

Towards a Linear-Scaling DFT Technique: The Density Matrix Approach

A recently proposed linear-scaling scheme for density-functional pseudopotential calculations is described in detail. The method is based on a formulation of density functional theory in which the ground state energy is determined by minimization with respect to the density matrix, subject to the condition that the eigenvalues of the latter lie in the range [0,1]. Linear-scaling behavior is achieved by requiring that the density matrix should vanish when the separation of its arguments exceeds a chosen cutoff. The limitation on the eigenvalue range is imposed by the method of Li, Nunes and Vanderbilt. The scheme is implemented by calculating all terms in the energy on a uniform real-space grid, and minimization is performed using the conjugate-gradient method. Tests on a 512-atom Si system show that the total energy converges rapidly as the range of the density matrix is increased. A discussion of the relation between the present method and other linear-scaling methods is given, and some problems that still require solution are indicated.Comment: REVTeX file, 27 pages with 4 uuencoded postscript figure

Review of AdS/CFT Integrability, Chapter IV.3: N=6 Chern-Simons and Strings on AdS4xCP3

We review the duality and integrability of N=6 superconformal Chern-Simons theory in three dimensions and IIA superstring theory on the background AdS4xCP3. We introduce both of these models and describe how their degrees of freedom are mapped to excitations of a long-range integrable spin-chain. Finally, we discuss the properties of the Bethe equations, the S-matrix and the algebraic curve that are special to this correspondence and differ from the case of N=4 SYM theory and strings on AdS5xS5.Comment: 22 pages, see also overview article arXiv:1012.3982, v2: references to other chapters updated, v3: references added, v4: brief discussion of giant magnons added, further minor changes, published version, v5: union of v3 and v4 because changes made in v3 were accidentally lost in v

Measurement of Through-Going Particle Momentum By Means Of Multiple Scattering With The ICARUS T600 TPC

The ICARUS collaboration has demonstrated, following the operation of a 600 ton (T600) detector at shallow depth, that the technique based on liquid Argon TPCs is now mature. The study of rare events, not contemplated in the Standard Model, can greatly benefit from the use of this kind of detectors. In particular, a deeper understanding of atmospheric neutrino properties will be obtained thanks to the unprecedented quality of the data ICARUS provides. However if we concentrate on the T600 performance, most of the $\nu_\mu$ charged current sample will be partially contained, due to the reduced dimensions of the detector. In this article, we address the problem of how well we can determine the kinematics of events having partially contained tracks. The analysis of a large sample of atmospheric muons collected during the T600 test run demonstrate that, in case the recorded track is at least one meter long, the muon momentum can be reconstructed by an algorithm that measures the Multiple Coulomb Scattering along the particle's path. Moreover, we show that momentum resolution can be improved by a factor two using an algorithm based on the Kalman Filtering technique

Energy reconstruction of electromagnetic showers from [Pi 0] decays with the ICARUS T600 liquid argon TPC

We discuss the ICARUS T600 detector capabilities in electromagnetic shower reconstruction through the analysis of a sample of 212 events, coming from the 2001 Pavia surface test run, of hadronic interactions leading to the production of 0 mesons. Methods of shower energy and shower direction measurements were developed and the invariant mass of the photon pairs was reconstructed. The ( ) invariant mass was found to be consistent with the value of the 0 mass. The resolution of the reconstructed 0 mass was found to be equal to 27.3%. An improved analysis, carried out in order to clean the full event sample from the events measured in the crowded environment, mostly due to the trigger conditions, gave a 0 mass resolution of 16.1%, significantly better than the one evaluated for the full event sample. The trigger requirement of the coincidence of at least four photo-multiplier signals favored the selection of events with a strong pile up of cosmic ray tracks and interactions. Hence a number of candidate 0 events were heavily contaminated by other tracks and had to be rejected. Monte Carlo simulations of events with 0 production in hadronic and neutrino interactions confirmed the validity of the shower energy and shower direction reconstruction methods applied to the real data