A variational framework for spectral discretization of the density matrix in Kohn Sham density functional theory

Kohn-Sham density functional theory (KSDFT) is currently the main work-horse of quantum mechanical calculations in physics, chemistry, and materials science. From a mechanical engineering perspective, we are interested in studying the role of defects in the mechanical properties in materials. In real materials, defects are typically found at very small concentrations e.g., vacancies occur at parts per million, dislocation density in metals ranges from $10^{10} m^{-2}$ to $10^{15} m^{-2}$, and grain sizes vary from nanometers to micrometers in polycrystalline materials, etc. In order to model materials at realistic defect concentrations using DFT, we would need to work with system sizes beyond millions of atoms. Due to the cubic-scaling computational cost with respect to the number of atoms in conventional DFT implementations, such system sizes are unreachable. Since the early 1990s, there has been a huge interest in developing DFT implementations that have linear-scaling computational cost. A promising approach to achieving linear-scaling cost is to approximate the density matrix in KSDFT. The focus of this thesis is to provide a firm mathematical framework to study the convergence of these approximations. We reformulate the Kohn-Sham density functional theory as a nested variational problem in the density matrix, the electrostatic potential, and a field dual to the electron density. The corresponding functional is linear in the density matrix and thus amenable to spectral representation. Based on this reformulation, we introduce a new approximation scheme, called spectral binning, which does not require smoothing of the occupancy function and thus applies at arbitrarily low temperatures. We proof convergence of the approximate solutions with respect to spectral binning and with respect to an additional spatial discretization of the domain. For a standard one-dimensional benchmark problem, we present numerical experiments for which spectral binning exhibits excellent convergence characteristics and outperforms other linear-scaling methods.

On certain discrete inequalities and their continous analogues

In a 1955 paper, Ky Fan, Olga Taussky, and John Todd presented discrete analogues of inequalities of Wirtinger type, and by taking limits they were able to recover the continuous inequalities. We generalize their techniques to mixed and higher derivatives and inequalities with weight functions in the integrals. We have also considered analogues of inequalities of Müller and Redheffer and have used these inequalities to derive a necessary and sufficient condition on ordered pairs of numbers so that the first number is the square norm of the k^th derivative of some periodic function and the second number is the square norm of the m^th derivative of the same periodic function.

Energy inequalities and error estimates for axisymmetric torsion of thin elastic shells of revolution

The problem motivating this investigation is that of pure axisymmetric torsion of an elastic shell of revolution. The analysis is carried out within the framework of the three-dimensional linear theory of elastic equilibrium for homogeneous, isotropic solids. The objective is the rigorous estimation of errors involved in the use of approximations based on thin shell theory.

The underlying boundary value problem is one of Neumann type for a second order elliptic operator. A systematic procedure for constructing pointwise estimates for the solution and its first derivatives is given for a general class of second-order elliptic boundary-value problems which includes the torsion problem as a special case.

The method used here rests on the construction of “energy inequalities” and on the subsequent deduction of pointwise estimates from the energy inequalities. This method removes certain drawbacks characteristic of pointwise estimates derived in some investigations of related areas.

Special interest is directed towards thin shells of constant thickness. The method enables us to estimate the error involved in a stress analysis in which the exact solution is replaced by an approximate one, and thus provides us with a means of assessing the quality of approximate solutions for axisymmetric torsion of thin shells.

Finally, the results of the present study are applied to the stress analysis of a circular cylindrical shell, and the quality of stress estimates derived here and those from a previous related publication are discussed.

Part I. Approximate Hartree-Fock wavefunctions one-electron properties and electronic structure of the water molecule. Part II. Perturbation-variational calculation of the nuclear spin-spin isotope coupling constant in HD

Part I

Several approximate Hartree-Fock SCF wavefunctions for the ground electronic state of the water molecule have been obtained using an increasing number of multicenter s, p, and d Slater-type atomic orbitals as basis sets. The predicted charge distribution has been extensively tested at each stage by calculating the electric dipole moment, molecular quadrupole moment, diamagnetic shielding, Hellmann-Feynman forces, and electric field gradients at both the hydrogen and the oxygen nuclei. It was found that a carefully optimized minimal basis set suffices to describe the electronic charge distribution adequately except in the vicinity of the oxygen nucleus. Our calculations indicate, for example, that the correct prediction of the field gradient at this nucleus requires a more flexible linear combination of p-orbitals centered on this nucleus than that in the minimal basis set. Theoretical values for the molecular octopole moment components are also reported.

Part II

The perturbation-variational theory of R. M. Pitzer for nuclear spin-spin coupling constants is applied to the HD molecule. The zero-order molecular orbital is described in terms of a single 1s Slater-type basis function centered on each nucleus. The first-order molecular orbital is expressed in terms of these two functions plus one singular basis function each of the types e^-r/r and e^-r ln r centered on one of the nuclei. The new kinds of molecular integrals were evaluated to high accuracy using numerical and analytical means. The value of the HD spin-spin coupling constant calculated with this near-minimal set of basis functions is J_HD = +96.6 cps. This represents an improvement over the previous calculated value of +120 cps obtained without using the logarithmic basis function but is still considerably off in magnitude compared with the experimental measurement of J_HD = +43 0 ± 0.5 cps.

Changes in social inequalities in disability-free life expectancy in Southern Europe: the case of the Basque Country

Background: Health expectancy is a useful tool to monitor health inequalities. The evidence about the recent changes in social inequalities in healthy expectancy is relatively scarce and inconclusive, and most studies have focused on Anglo-Saxon and central or northern European countries. The objective of this study was to analyse the changes in socioeconomic inequalities in disability-free life expectancy in a Southern European population, the Basque Country, during the first decade of the 21st century. Methods: This was an ecological cross-sectional study of temporal trends on the Basque population in 1999-2003 and 2004-2008. All-cause mortality rate, life expectancy, prevalence of disability and disability free-life expectancy were calculated for each period according to the deprivation level of the area of residence. The slope index of inequality and the relative index of inequality were calculated to summarize and compare the inequalities in the two periods. Results: Disability free-life expectancy decreased as area deprivation increased both in men and in women. The difference between the most extreme groups in 2004-2008 was 6.7 years in men and 3.7 in women. Between 1999-2003 and 2004-2008, socioeconomic inequalities in life expectancy decreased, and inequalities in disability-free expectancy increased in men and decreased in women. Conclusions: This study found important socioeconomic inequalities in health expectancy in the Basque Country. These inequalities increased in men and decreased in women in the first decade of the 21st century, during which the Basque Country saw considerable economic growth.

The impact of economic crises on social inequalities in health: what do we know so far?

Since 2008, Western countries are going through a deep economic crisis whose health impacts seem to be fundamentally counter-cyclical: when economic conditions worsen, so does health, and mortality tends to rise. While a growing number of studies have presented evidence on the effect of crises on the average population health, a largely neglected aspect of research is the impact of crises and the related political responses on social inequalities in health, even if the negative consequences of the crises are primarily borne by the most disadvantaged populations. This commentary will reflect on the results of the studies that have analyzed the effect of economic crises on social inequalities in health up to 2013. With some exceptions, the studies show an increase in health inequalities during crises, especially during the Southeast Asian and Japanese crises and the Soviet Union crisis, although it is not always evident for both sexes or all health or socioeconomic variables. In the Nordic countries during the nineties, a clear worsening of health equity did not occur. Results about the impacts of the current economic recession on health equity are still inconsistent. Some of the factors that could explain this variability in results are the role of welfare state policies, the diversity of time periods used in the analyses, the heterogeneity of socioeconomic and health variables considered, the changes in the socioeconomic profile of the groups under comparison in times of crises, and the type of measures used to analyze the magnitude of social inequalities in health. Social epidemiology should further collaborate with other disciplines to help produce more accurate and useful evidence about the relationship between crises and health equity.