Reply to 'Comment on "Extending Hirshfeld-I to bulk and periodic materials" '

The issues raised in the comment by T.A. Manz are addressed through the presentation of calculated atomic charges for NaF, NaCl, MgO, SrTiO$_3$ and La$_2$Ce$_2$O$_7$, using our previously presented method for calculating Hirshfeld-I charges in Solids [J. Comput. Chem.. doi: 10.1002/jcc.23088]. It is shown that the use of pseudo-valence charges is sufficient to retrieve the full all-electron Hirshfeld-I charges to good accuracy. Furthermore, we present timing results of different systems, containing up to over $200$ atoms, underlining the relatively low cost for large systems. A number of theoretical issues is formulated, pointing out mainly that care must be taken when deriving new atoms in molecules methods based on "expectations" for atomic charges.Comment: 7 pages, 2 Tables, 2 figure

Extending Hirshfeld-I to bulk and periodic materials

In this work, a method is described to extend the iterative Hirshfeld-I method, generally used for molecules, to periodic systems. The implementation makes use of precalculated pseudo-potential based charge density distributions, and it is shown that high quality results are obtained for both molecules and solids, such as ceria, diamond, and graphite. The use of such grids makes the implementation independent of the solid state or quantum chemical code used for studying the system. The extension described here allows for easy calculation of atomic charges and charge transfer in periodic and bulk systems.Comment: 11 pages, 4 Tables, 5 Figures, pre-referee draft only, much extended post referee version only available at publishe

Three-dimensional superconductors with hybrid higher order topology

We consider three dimensional superconductors in class DIII with a four-fold rotation axis and inversion symmetry. It is shown that such systems can exhibit higher order topology with helical Majorana hinge modes. In the case of even-parity superconductors we show that higher order topological superconductors can be obtained by adding a small pairing with the appropriate $C_4$ symmetry implementation to a topological insulator. We also show that a hybrid case is possible, where Majorana surface cones resulting from non-trivial strong topology coexist with helical hinge modes. We propose a bulk invariant detecting this hybrid scenario, and numerically analyse a tight binding model exhibiting both Majorana cones and hinge modes.Comment: Published versio

A primal-dual semidefinite programming algorithm tailored to the variational determination of the two-body density matrix

The quantum many-body problem can be rephrased as a variational determination of the two-body reduced density matrix, subject to a set of N-representability constraints. The mathematical problem has the form of a semidefinite program. We adapt a standard primal-dual interior point algorithm in order to exploit the specific structure of the physical problem. In particular the matrix-vector product can be calculated very efficiently. We have applied the proposed algorithm to a pairing-type Hamiltonian and studied the computational aspects of the method. The standard N-representability conditions perform very well for this problem.Comment: 24 pages, 5 figures, submitted to the Journal of Computational Physic

Tuning of CeO2_2 buffer layers for coated superconductors through doping

The appearance of microcracks in CeO$_2$ buffer layers, as used in buffer layer architectures for coated superconductors, indicates the presence of stress between this buffer layer and the substrate. This stress can originate from the differences in thermal expansion or differences in lattice parameters between the CeO$_2$ buffer layer and the substrate. In this article, we study, by means of \textit{ab initio} density functional theory calculations, the influence of group IV doping elements on the lattice parameter and bulk modulus of CeO$_2$. Vegard's law behavior is found for the lattice parameter in systems without oxygen vacancies, and the Shannon crystal radii for the doping elements are retrieved from the lattice expansions. We show that the lattice parameter of the doped CeO$_2$ can be matched to that of the La$_2$Zr$_2$O$_7$ coated NiW substrate substrate for dopant concentrations of about $5\%$, and that bulk modulus matching is either not possible or would require extreme doping concentrations.Comment: 5 pages, 1 table, 2 figures, EMRS 2011 Fall meeting symposium on Stress, structure and stoichiometry effects on nanomaterial

A mechanism for anomalous Hall ferromagnetism in twisted bilayer graphene

We use a lowest Landau level model to study the recent observation of an anomalous Hall effect in twisted bilayer graphene. This effective model is rooted in the occurrence of Chern bands which arise due to the coupling between the graphene device and its encapsulating substrate. Our model exhibits a phase transition from a spin-valley polarized insulator to a partial or fully valley unpolarized metal as the bandwidth is increased relative to the interaction strength, consistent with experimental observations. In sharp contrast to standard quantum Hall ferromagnetism, the Chern number structure of the flat bands precludes an instability to an inter-valley coherent phase, but allows for an excitonic vortex lattice at large interaction anisotropy.Comment: 7 + 11 pages, 4 + 2 figures; (v2) 5 + 12 pages, 4 + 3 figures, added figures and references; (v3) 5 + 13 pages, minor changes, added a comment on Landau fans. Published versio