In press, SIAM J. Matrix Anal. Appl.Abstract. In this paper we develop new Newton and conjugate gradient algorithms on the Grassmann and Stiefel manifolds. These manifolds represent the constraints that arise in such areas as the symmetric eigenvalue problem, nonlinear eigenvalue problems, electronic structures computations, and signal processing. In addition to the new algorithms, we show how the geometrical framework gives penetrating new insights allowing us to create, understand, and compare algorithms. The theory proposed here provides a taxonomy for numerical linear algebra algorithms that provide a top level mathematical view of previously unrelated algorithms. It is our hope that developers of new algorithms and perturbation theories will benefit from the theory, methods, and examples in this paper.

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes. C 2016 AIP Publishing LLC. [http://dx.

Nanoelectromechanical systems (NEMS) hold promise for a number of scientific and technological applications. In particular, NEMS oscillators have been proposed for use in ultrasensitive mass detection, radio-frequency signal processing, and as a model system for exploring quantum phenomena in macroscopic systems. Perhaps the ultimate material for these applications is a carbon nanotube. They are the stiffest material known, have low density, ultrasmall cross-sections and can be defect-free. Equally important, a nanotube can act as a transistor and thus may be able to sense its own motion. In spite of this great promise, a room-temperature, self-detecting nanotube oscillator has not been realized, although some progress has been made. Here we report the electrical actuation and detection of the guitar-string-like oscillation modes of doubly clamped nanotube oscillators. We show that the resonance frequency can be widely tuned and that the devices can be used to transduce very small forces.

Electron scattering rates in metallic single-walled carbon nanotubes are studied using an atomic force microscope as an electrical probe. From the scaling of the resistance of the same nanotube with length in the low and high bias regimes, the mean free paths for both regimes are inferred. The observed scattering rates are consistent with calculations for acoustic phonon scattering at low biases and zone boundary/optical phonon scattering at high biases.

This work explores the use of joint density-functional theory, a new form of density-functional theory for the ab initio description of electronic systems in thermodynamic equilibrium with a liquid environment, to describe electrochemical systems. After reviewing the physics of the underlying fundamental electrochemical concepts, we identify the mapping between commonly measured electrochemical observables and microscopically computable quantities within an, in principle, exact theoretical framework. We then introduce a simple, computationally efficient approximate functional which we find to be quite successful in capturing a priori basic electrochemical phenomena, including the capacitive Stern and diffusive Gouy-Chapman regions in the electrochemical double layer, quantitative values for interfacial capacitance, and electrochemical potentials of zero charge for a series of metals. We explore surface charging with applied potential and are able to place our ab initio results directly on the scale associated with the Standard Hydrogen Electrode (SHE). Finally, we provide explicit details for implementation within standard density-functional theory software packages at negligible computational cost over standard calculations carried out within vacuum environments.