Bootstrap approximation for the exchange-correlation kernel of time-dependent density functional theory

A new parameter-free approximation for the exchange-correlation kernel $f_{\rm xc}$ of time-dependent density functional theory is proposed. This kernel is expressed as an algorithm in which the exact Dyson equation for the response as well as a further approximate condition are solved together self-consistently leading to a simple parameter-free kernel. We apply this to the calculation of optical spectra for various small bandgap (Ge, Si, GaAs, AlN, TiO$_2$, SiC), large bandgap (C, LiF, Ar, Ne) and magnetic (NiO) insulators. The calculated spectra are in very good agreement with experiment for this diverse set of materials, highlighting the universal applicability of the new kernel.Comment: 4 figures 5 page

Multigrid Preconditioning for a Space-Time Spectral-Element Discontinuous-Galerkin Solver

In this work we examine a multigrid preconditioning approach in the context of a high- order tensor-product discontinuous-Galerkin spectral-element solver. We couple multigrid ideas together with memory lean and efficient tensor-product preconditioned matrix-free smoothers. Block ILU(0)-preconditioned GMRES smoothers are employed on the coarsest spaces. The performance is evaluated on nonlinear problems arising from unsteady scale- resolving solutions of the Navier-Stokes equations: separated low-Mach unsteady ow over an airfoil from laminar to turbulent ow. A reduction in the number of ne space iterations is observed, which proves the efficiency of the approach in terms of preconditioning the linear systems, however this gain was not reflected in the CPU time. Finally, the preconditioner is successfully applied to problems characterized by stiff source terms such as the set of RANS equations, where the simple tensor product preconditioner fails. Theoretical justification about the findings is reported and future work is outlined

Low-density silicon allotropes for photovoltaic applications

Silicon materials play a key role in many technologically relevant fields, ranging from the electronic to the photovoltaic industry. A systematic search for silicon allotropes was performed by employing a modified ab initio minima hopping crystal structure prediction method. The algorithm was optimized to specifically investigate the hitherto barely explored low-density regime of the silicon phase diagram by imitating the guest-host concept of clathrate compounds. In total 44 metastable phases are presented, of which 11 exhibit direct or quasi-direct band-gaps in the range of $\approx$1.0-1.8 eV, close to the optimal Shockley-Queisser limit of $\approx$1.4 eV, with a stronger overlap of the absorption spectra with the solar spectrum compared to conventional diamond silicon. Due to the structural resemblance to known clathrate compounds it is expected that the predicted phases can be synthesized

First-principles predicted low-energy structures of NaSc(BH4)4

According to previous interpretations of experimental data, sodium-scandium double-cation borohydride NaSc(BH$_4$)$_4$ crystallizes in the crystallographic space group $Cmcm$ where each sodium (scandium) atom is surrounded by six scandium (sodium) atoms. A careful investigation of this phase based on \textit{ab initio} calculations indicates that the structure is dynamically unstable and gives rise to an energetically and dynamically more favorable phase with $C222_1$ symmetry and nearly identical x-ray diffraction pattern. By additionally performing extensive structural searches with the minima-hopping method we discover a class of new low-energy structures exhibiting a novel structural motif in which each sodium (scandium) atom is surrounded by four scandium (sodium) atoms arranged at the corners of either a rectangle with nearly equal sides or a tetrahedron. These new phases are all predicted to be insulators with band gaps of $7.9-8.2$ eV. Finally, we estimate the influence of these structures on the hydrogen-storage performance of NaSc(BH$_4$)$_4$.Comment: Version publishe

Potentials of mean force in acidic proton transfer reactions in constrained geometries

Free energy barriers associated with the transfer of an excess proton in water and related to the potentials of mean force in proton transfer episodes have been computed in a wide range of thermodynamic states, from low-density amorphous ices to high-temperature liquids under the critical point for unconstrained and constrained systems. The latter were represented by set-ups placed inside hydrophobic graphene slabs at the nanometric scale allocating a few water layers, namely one or two in the narrowest case. Water–proton and carbon–proton forces were modelled with a Multi-State Empirical Valence Bond method. As a general trend, a competition between the effects of confinement and temperature is observed on the local hydrogen-bonded structures around the lone proton and, consequently, on the mean force exerted by its environment on the water molecule carrying the proton. Free energy barriers estimated from the computed potentials of mean force tend to rise with the combined effect of increasing temperatures and the packing effect due to a larger extent of hydrophobic confinement. The main reason observed for such enhancement of the free energy barriers was the breaking of the second coordination shell around the lone proton.Postprint (author's final draft

Prediction of a novel monoclinic carbon allotrope

A novel allotrope of carbon with $P2/m$ symmetry was identified during an \emph{ab-initio} minima-hopping structural search which we call $M10$-carbon. This structure is predicted to be more stable than graphite at pressures above 14.4 GPa and consists purely of $sp^3$ bonds. It has a high bulk modulus and is almost as hard as diamond. A comparison of the simulated X-ray diffraction pattern shows a good agreement with experimental results from cold compressed graphite.Comment: 3 pages, 3 figure

h-multigrid agglomeration based solution strategies for discontinuous Galerkin discretizations of incompressible flow problems

In this work we exploit agglomeration based $h$-multigrid preconditioners to speed-up the iterative solution of discontinuous Galerkin discretizations of the Stokes and Navier-Stokes equations. As a distinctive feature $h$-coarsened mesh sequences are generated by recursive agglomeration of a fine grid, admitting arbitrarily unstructured grids of complex domains, and agglomeration based discontinuous Galerkin discretizations are employed to deal with agglomerated elements of coarse levels. Both the expense of building coarse grid operators and the performance of the resulting multigrid iteration are investigated. For the sake of efficiency coarse grid operators are inherited through element-by-element $L^2$ projections, avoiding the cost of numerical integration over agglomerated elements. Specific care is devoted to the projection of viscous terms discretized by means of the BR2 dG method. We demonstrate that enforcing the correct amount of stabilization on coarse grids levels is mandatory for achieving uniform convergence with respect to the number of levels. The numerical solution of steady and unsteady, linear and non-linear problems is considered tackling challenging 2D test cases and 3D real life computations on parallel architectures. Significant execution time gains are documented.Comment: 78 pages, 7 figure